• Journal of the Korean Geotechnical Society
• /
• v.29 no.7
• /
• pp.75-89
• /
• 2013

This paper describes a GIS-based geohydrologic methodology, called YSGWF (YonSei GroundWater Flow) for predicting the rainfall infiltration-groundwater flow of slopes. This physical-based model was developed by the combination of modified Green-Ampt model that considers the unsaturated soil parameters and GIS-based raster model using Darcy's law that reflects the groundwater flow. In the model, raster data are used to simulate the three dimensional inclination of bedrock surface as actual topographic data, and the groundwater flow is governed by the slope. Also, soil profile is ideally subdivided into three zones, i.e., the wetting band zone, partially saturated zone, and fully saturated zone. In the wetting band and partially saturated zones the vertical infiltration of water (rainfall) from surface into ground is modeled. When the infiltrated water recharges into the fully saturated zone, the horizontal flow of groundwater is introduced. A comparison between the numerical calculation and real landslide data shows a reasonable agreement, which indicate that the model can be used to simulate real rainfall infiltration-groundwater flow.

• 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.

• Journal of the Korean Institute of Gas
• /
• v.16 no.3
• /
• pp.35-41
• /
• 2012

Manifold researches for carbon capture and storage (CCS) have been developed and large scale-carbon capture system can be performed recently. Hence, the technologies for $CO_2$ sequestration or storage become necessary to handle the captured $CO_2$. Among them, enhanced oil recovery using $CO_2$ can be a solution since it guarantees both oil recovery and $CO_2$ sequestration. In this study, the miscible flow of oil and $CO_2$ in porous media is modeled to analyze the effect of enhanced oil recovery and $CO_2$ sequestration. Based on Darcy-Muskat law, the equation is modified to consider miscibility of oil and $CO_2$ and the change of viscosity. Finite volume method is used for numerical modeling. As results, the pressure and oil saturation changes with time can be predicted when oil, water, and $CO_2$ are injected, respectively, and $CO_2$ injection is more efficient than water injection for oil recovery.

• Journal of the Korean Society of Safety
• /
• v.11 no.4
• /
• pp.143-150
• /
• 1996

This is 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 and 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 Ship and Ocean Technology
• /
• v.1 no.1
• /
• pp.35-47
• /
• 1997

The consequences associated with ships running aground depend very much on the soil characteristics of the sea bed and the geometrical shape of the ship bow. The penetration into the sea bed depends on these factors and the penetration is an important factor for the ship motion because it influences the ship heave and pitch motions as well as the friction between the ship and the soil. In this paper a rational calculation model is presented for the sea bed soil reaction forces on the ship bottom. The model is based on the assumption that the penetration of the ship bow generates a flow of pore water through the grain skeleton of the soil. The flow is governed by Darcy\`s law and it is driven by the pressure of the pore water at the bow. In addition to this pore water pressure, the bow is subjected to the effective stresses in the grain skeleton at the bow surface. These stresses are determined by the theory of frictional soils in rupture. Frictional stresses on the bow surface are assumed to be related to the normal pressure by a simple Coulomb relation. The total soil reaction as a function of velocity and penetration is found by integration of normal pressure and frictional stresses over the surface of the bow. The analysis procedure is implemented in a computer program for time domain rigid body analysis of ships running aground and it is verified in the paper through a comparison of calculated stopping lengths, effective coefficients of friction, and sea bed penetrations with corresponding experimental results obtained by model tests as well as large, scale tests.

• Transactions of Materials Processing
• /
• v.6 no.3
• /
• pp.239-249
• /
• 1997

The behaviour of alloys in the semi-solid state strongly depends on the imposed stress state and on the morphology of the phase which can vary from dendritic to globular. To optimal net shape forging of semi-solid materials, it is important to investigate for filling phenomena in forging process of arbitrarily shaped dies. To produce a automotive part which has good mechanical property, the filling pattern according to die velocity and solid fraction distribution has to be estimated for arbitrarily shaped dies. Therefore, the estimation of filling characteristic in the forging simulation with arbitrarily shaped dies of semi-solid materials are calculated by finite element method with proposed algorithm. The proposed theoretical model and a various boundary conditions for arbitrarily shaped dies is investigated with the coupling calculation between the liquid phase flow and the solid phase deformation. The simulation process with arbitrarily shaped dies is performed to the isothermal conditions of two dimensional problems. To analysis of forging process by using semi-solid materials, a new stress-strain relationship is described, and forging analysis is performed by viscoelastic model for the solid phase and the Darcy's law for the liquid flow. The calculated results for forging force and filling limitations will be compared to experimental data. The filling simulation of simple products performed with the uniform billet temperature(584$^{\circ}C$) from the induction heating by the commercial package MAGMAsoft. The initial step of computation is the touching of semi-solid material with the end of die gate and the initial concept of proposed system just fit with the capability of MAGMAsoft.

• Journal of Soil and Groundwater Environment
• /
• v.16 no.5
• /
• pp.53-66
• /
• 2011

Groundwater which infiltrated in recharge areas discharges in the forms of evapotranspiration, baseflow to streams, groundwater abstraction and eventually flows into the sea. This study characterized radon-222 concentration and electrical conductivity (EC) in coastal groundwater discharge, well groundwater, Ilkwang Stream water, and seawater in the coastal area of Busan Metropolitan City and subsequently estimated groundwater discharge rate to the sea. The median value of Rn-222 concentration is highest in well groundwater (18.36 Bq/L), and then decreases in the order of coastal groundwater discharge (15.92 Bq/L), Ilkwang Stream water (1.408 Bq/L), and seawater (0.030 Bq/L). The relationship between Rn-222 concentration and EC values is relatively strong in well groundwater and then in seawater. However, the relationship is not visible between coastal groundwater discharge and Ilkwang Stream water. The groundwater discharge rate to the sea is estimated as $3,130m^3$/day by using radon mass budget model and $16,788m^3$/day by using Darcy's law.

• Journal of The Korean Society of Agricultural Engineers
• /
• v.63 no.5
• /
• pp.63-71
• /
• 2021

Windbreak forests, which have a windproof effect against strong winds, are known to be effective in reducing the spread of odors and dust emitted from livestock farms. The effect of reducing the spread of odors and dust can be estimated through numerical models such as computational fluid dynamics, which require aerodynamic coefficients of the windbreaks for accurate prediction of their performance. In this study, we aimed to evaluate the aerodynamic coefficients, Co, C₁, C₂, and α, of two windbreaks, Thuja occidentalis and a mesh net, through wind tunnel experiments. The aerodynamic coefficients were derived by the relation between the incoming wind speed and the pressure loss due to the windbreaks which was measured by differential pressure sensors. In order to estimate the change in the aerodynamic coefficient concerning various leaf density, the experiments were conducted repeatedly by removing the leaves gradually in various stages. The results showed that the power law regression model more suitable for coefficient evaluation compared to the Darcy-Forchheimer model.

• 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.

• Tunnel and Underground Space
• /
• v.28 no.5
• /
• pp.400-425
• /
• 2018

This study presents the research results and current status of the DECOVALEX-2019 project Task B. Task B named 'Fault slip modelling' is aiming at developing a numerical method to simulate the coupled hydro-mechanical behavior of fault, including slip or reactivation, induced by water injection. The first research step of Task B is a benchmark simulation which is designed for the modelling teams to familiarize themselves with the problem and to set up their own codes to reproduce the hydro-mechanical coupling between the fault hydraulic transmissivity and the mechanically-induced displacement. We reproduced the coupled hydro-mechanical process of fault slip 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. A methodology to formulate the hydro-mechanical coupling relations of two different hydraulic aperture models and link the solid element of TOUGH2 and the interface element of FLAC3D was suggested. In addition, we developed a coupling module to update the changes in geometric features (mesh) and hydrological properties of fault caused by water injection at every calculation step for TOUGH-FLAC simulator. Then, the transient responses of the fault, including elastic deformation, reactivation, progressive evolutions of pathway, pressure distribution and water injection rate, to stepwise pressurization were examined during the simulations. The results of the simulations suggest that the developed model can provide a reasonable prediction of the hydro-mechanical behavior related to fault reactivation. The numerical model will be enhanced by continuing collaboration and interaction with other research teams of DECOLVAEX-2019 Task B and validated using the field data from fault activation experiments in a further study.