• Journal of Ocean Engineering and Technology
• /
• v.32 no.4
• /
• pp.228-236
• /
• 2018

The performance of a porous swash bulkhead for the reduction of the resonant liquid motion in a swaying rectangular tank was investigated based on the assumption of linear potential theory. The Galerkin method (Porter and Evans, 1995) was used to solve the potential flow model by adding a viscous frictional damping term to the free-surface condition. By comparing the experimental results and the analytical solutions, we verified that the frictional damping coefficient was 0.4. Darcy's law was used to consider the energy dissipation at a porous bulkhead. The tool that was developed with a built-in frictional damping coefficient of 0.4 was confirmed by small-scale experiments. Using this tool, the free-surface elevation, hydrodynamic force (added mass, damping coefficient) on a wall, and the horizontal load on a bulkhead were assessed for various combinations of porosity and submergence depth. It was found that the vertical porous bulkhead can suppress sloshing motions significantly when properly designed and by selecting the appropriate porosity(${\approx}0.1$) and submergence depth.

• The Journal of Engineering Geology
• /
• v.23 no.4
• /
• pp.447-455
• /
• 2013

A micro-mechanical pipe network model with the shape of a cube was developed to simulate the behavior of fluid flow through a porous medium. The fluid-flow mechanism through the cubic pipe network channels was defined mainly by introducing a well-known percolation theory (Stauffer and Aharony, 1994). A non-uniform flow generally appeared because all of the pipe diameters were allocated individually in a stochastic manner based on a given pore-size distribution curve and porosity. Fluid was supplied to one surface of the pipe network under a certain driving pressure head and allowed to percolate through the pipe networks. A percolation condition defined by capillary pressure with respect to each pipe diameter was applied first to all of the network pipes. That is, depending on pipe diameter, the fluid may or may not penetrate a specific pipe. Once pore pressures had reached equilibrium and steady-state flow had been attained throughout the network system, Darcy's law was used to compute the resultant permeability. This study investigated the sensitivity of network size to permeability calculations in order to find out the optimum network size which would be used for all the network modelling in this study. Mean pore size and pore size distribution curve obtained from field are used to define each of pipe sizes as being representative of actual oil sites. The calculated and measured permeabilities are in good agreement.

• Tunnel and Underground Space
• /
• v.29 no.3
• /
• pp.197-213
• /
• 2019

We simulated the fault reactivation experiment conducted at 'Main Fault' intersecting the low permeability clay formations of Mont Terri Underground Research Laboratory in Switzerland using TOUGH-FLAC simulator. The fluid flow along a fault was modelled with solid elements and governed by Darcy's law with the cubic law in TOUGH2, whereas the mechanical behavior of a single fault was represented by creating interface elements between two separating rock blocks in FLAC3D. We formulate the hydro-mechanical coupling relation of hydraulic aperture to consider the elastic fracture opening and failure-induced dilation for reproducing the abrupt changes in injection flow rate and monitoring pressure at fracture opening pressure. A parametric study was conducted to examine the effects of in-situ stress condition and fault deformation and strength parameters and to find the optimal parameter set to reproduce the field observations. In the best matching simulation, the fracture opening pressure and variations of injection flow rate and monitoring pressure showed good agreement with field experiment results, which suggests the capability of the numerical model to reasonably capture the fracture opening and propagation process. The model overestimated the fault displacement in shear direction and the range of reactivated zone, which was attributed to the progressive shear failures along the fault at high injection pressure. In the field experiment results, however, fracture tensile opening seems the dominant mechanism affecting the hydraulic aperture increase.

• Coupled systems mechanics
• /
• v.7 no.6
• /
• pp.649-668
• /
• 2018

In this paper, we present a numerical model for fluid-structure interaction between structure built of porous media and acoustic fluid, which provides both pore pressure inside porous media and hydrodynamic pressures and hydrodynamic forces exerted on the upstream face of the structure in an unified manner and simplifies fluid-structure interaction problems. The first original feature of the proposed model concerns the structure built of saturated porous medium whose response is obtained with coupled discrete beam lattice model, which is based on Voronoi cell representation with cohesive links as linear elastic Timoshenko beam finite elements. The motion of the pore fluid is governed by Darcy's law, and the coupling between the solid phase and the pore fluid is introduced in the model through Biot's porous media theory. The pore pressure field is discretized with CST (Constant Strain Triangle) finite elements, which coincide with Delaunay triangles. By exploiting Hammer quadrature rule for numerical integration on CST elements, and duality property between Voronoi diagram and Delaunay triangulation, the numerical implementation of the coupling results with an additional pore pressure degree of freedom placed at each node of a Timoshenko beam finite element. The second original point of the model concerns the motion of the outside fluid which is modeled with mixed displacement/pressure based formulation. The chosen finite element representations of the structure response and the outside fluid motion ensures for the structure and fluid finite elements to be connected directly at the common nodes at the fluid-structure interface, because they share both the displacement and the pressure degrees of freedom. Numerical simulations presented in this paper show an excellent agreement between the numerically obtained results and the analytical solutions.

• Geomechanics and Engineering
• /
• v.13 no.2
• /
• pp.285-300
• /
• 2017

Shallow coal resources are increasingly depleted, the mining has entered the deep stage. Due to "High stress, high gas, strong adsorption and low permeability" of coal seam, the gas drainage has become more difficult and the probability of coal and gas outburst accident increases. Based on the flow solid coupling theory of coal seam gas, the coupling model about stress and gas seepage of coal seam was set up by solid module and Darcy module in Comsol Multiphysics. The gas extraction effects were researched after applying hydraulic technology to increase permeability. The results showed that the effective influence radius increases with the expanded borehole radius and drainage time, decreases with initial gas pressure. The relationship between the effective influence radius and various factors presents in the form: $y=a+{\frac{b}{\left(1+{(\frac{x}{x_0})^p}\right)}}$. The effective influence radius with multiple boreholes is obviously larger than that of the single hole. According to the actual coal seam and gas geological conditions, appropriate layout way was selected to achieve the best effect. The field application results are consistent with the simulation results. It is found that the horizontal stress plays a very important role in coal seam drainage effect. The stress distribution change around the drilling hole will lead to the changes in porosity of coal seam, further resulting in permeability evolution and finally gas pressure distribution varies.

• International Journal of Aeronautical and Space Sciences
• /
• v.2 no.2
• /
• pp.65-72
• /
• 2001

This paper demonstrates an explicit-implicit, finite element analysis for linear as well as nonlinear hygrothermal stress problems. Additional features, such as moisture diffusion equation, crack element and virtual crack extension(VCE) method for evaluating J-integral are implemented in this program. The Linear Elastic Fracture Mechanics(LEFM) Theory is employed to estimate the crack driving force under the transient condition for an existing crack. Pores in materials are assumed to be saturated with moisture in the liquid form at the room temperature, which may vaporize as the temperature increases. The vaporization effects on the crack driving force are also studied. The ideal gas equation is employed to estimate the thermodynamic pressure due to vaporization at each time step after solving basic nodal values. A set of field equations governing the time dependent response of porous media are derived from balance laws based on the mixture theory. Darcy's law is assumed for the fluid flow through the porous media. Perzyna's viscoplastic model incorporating the Von-Mises yield criterion are implemented. The Green-Naghdi stress rate is used for the invariant of stress tensor under superposed rigid body motion. Isotropic elements are used for the spatial discretization and an iterative scheme based on the full Newton-Raphson method is used for solving the nonlinear governing equations.

• Journal of Ocean Engineering and Technology
• /
• v.15 no.2
• /
• pp.22-32
• /
• 2001

The numerical investigation of obliquely incident wave interactions with fully submerged dual buoy/porous-membrane floating breakwaters placed in parallel with spacing is studied based on linear potential theory and Darcy's law. The numerical solutions are obtained by using a discrete-membrane dynamic model and second-kind modified Bessel function distribution over the entire boundaries of fluid regions. First, numerical solutions for an idealized dual submerged system without buoys are obtained. Second, a more practical dual submerged system with membrane tension provided by buoys at its tops is investigated by the multi-domain boundary element method particularly devised for dual buoy/porous-membrane problems with gaps. The velocity potentials of wave motion are coupled with porous-membrane deformation, and solved simultaneously since the boundary condition on porous-membrane is not known in advance. The effects of varying permeability on membranes and wave characteristics are discussed for the optimum design parameters of systems previously studied. The inclusion of permeability on membrane eliminates the resonances that aggravate the breakwater performance. The system is highly efficient when waves generated by the buoys and membranes were mutually canceled and its energy at resonance frequency dissipates through fine pores on membranes.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2017.05a
• /
• pp.67-67
• /
• 2017

급속한 산업화와 도시화로 인하여, 투수지역은 감소함으로써, 개발전과 다른 지표, 지표하 유출이 나타난다. 이에 대한 대안으로 최근 저영향개발 (LID)이 수문학적 및 환경, 생태적 개선으로 대안으로 대두 되고 있다. 이에 많은 연구자들이 EPA SWMM 모델의 이용하여 LID 설치 전, LID를 모의하였으나, 불포화토양 및 토양 내의 matric head에 대한 고려가 없어 정확한 LID 모의가 힘든 실정이다. 이에 본 연구에서 상세한 토양 모의가 가능 HYDRUS를 이용하여, SWMM-HYDRUS 모델을 개발하였다. EPA SWMM 모델의 경우, 가장 상단의 layer에서 green ampt equation을 이용하여 침투량을 계산 후, 다음 layer에서 Darcy eqation을 이용하여 토양 물이동을 계산되어진다. 하지만 기존의 SWMM모델의 경우, 불포화토양내의 물 흐름에 대한 고려와 Matric head와 Pond depth에 대한 고려가 없어, LID 모의 시 한계점이 나타났다. 이에 본 연구에서는 이러한 한계점을 개선하기 위하여, 기존의 EPA SWMM의 LID 모듈을 Van Genuchten's equaton과 Richard Equation을 이용하여 정확한 토양 물 흐름을 계산하는 HYDRUS을 SWMM 모델에 결합하여, 더욱 정확한 LID 모의를 실시하였다. 개선된 SWMM-HYDRUS 모델의 모의 결과, 기존의 SWMM에서 한계점을 보여주는 Metric head를 고려하여 불포화 침투가 이루어지며, 또한 포화 후 LID 위에 존재하는 Pond depth를 고려해주는 결과가 나타났다. 향후 개발된 SWMM-HYDRUS모델를 이용하여 LID를 검증 시 기존의 모델보다 정확한 모의가 가능하다.

• Journal of the Korean GEO-environmental Society
• /
• v.8 no.5
• /
• pp.5-14
• /
• 2007

The quantity of noxious wastes generated by the growth in industrialization and population in all over the world and its potential hazards in subsurface environments are becoming increasingly significant. The extraction of the contaminant from the soil and movement of the water are restricted due to the low permeability and adsorption characteristics of the reclaimed soils. Incorporated technique with PVDs have been used for dewatering from fine-grained soils for the purpose of ground improvement by means of soil flushing and soil vapor extraction systems. This paper is to evaluate several key parameters that affected to the performance of the PVDs specifically with regard to: well resistance of PVD, zone of influence, and smear effects. In the feasibility of contaminant remediation was evaluated in pilot-scale laboratory experiments. Well resistance is affected on the vertical discharge capacity of the PVDs under the various vacuum pressures. The discharge capacity increases consistently in areal extents with higher applied vacuum up to a limiting vacuum pressure. The head values for each piezometer at different vacuum pressures show that the largest head loss occurs within 14 cm of the PVD. Air flow rates and head losses were measured for the PVD placed in the model test box and the gas permeability of the silty soils was calculated. Increasing the equivalent diameter results in a decrease in the calculated gas permeability. It is concluded that the gas permeability determined over the 1,500 to 2,000 $cm^3/s$ flow rates are the most accurate values which yields gas permeability of about 3.152 Darcy.

• Journal of Korea Water Resources Association
• /
• v.45 no.1
• /
• pp.39-52
• /
• 2012

Surface-subsurface interactions are an intrinsic component of the hydrologic response within a watershed. In general, these interactions are considered to be one of the most difficult areas of the discipline, particularly for the modeler who intends simulate the dynamic relations between these two major domains of the hydrological cycle. In essence, one major complexity is the spatial and temporal variations in the dynamically interacting system behavior. The proper simulation of these variations requires the need for providing an appropriate coupling mechanism between the surface and subsurface components of the system. In this study, an approach for modelling surface-subsurface flow and transport in a fully intergrated way is presented. The model uses the 2-dimensional diffusion wave equation for sheet surface water flow, and the Boussinesq equation with the Darcy's law and Dupuit-Forchheimer's assumption for variably saturated subsurface water flow. The coupled system of equations governing surface and subsurface flows is discretized using the finite volume method with central differencing in space and the Crank-Nicolson method in time. The interactions between surface and subsurface flows are considered mass balance based on the continuity conditions of pressure head and exchange flux. The major module consists of four sub-module (SUBFA, SFA, IA and NS module) is developed.