• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.31 no.5B
• /
• pp.483-489
• /
• 2011

In this study, the coupled-numerical model has been newly developed to investigate numerically scouring and deposition around a coastal structure like a submerged breakwater using a numerical wave model and a lagrangian particle model for sand transport. As a numerical wave model, LES-WASS-2D (Hur and Choi, 2008) is adopted. The model is able to consider the flow through a porous midium with inertial, laminar and turbulent resistance term and determine the eddy viscosity with LES turbulence model. Distinct element method (Cundall and Strack, 1979), which is able to apply to many dynamical analysis of particulate media, as a lagrangian particle model for sand transport is newly coupled to the numerical wave model. The numerical simulation has been carried out to examine the scour problem in front of an impermeable submerged breakwater using the newly coupled-numerical model. The numerical results has been compared qualitatively with an existing experimental data and then its applicability has been discussed.

• Interaction and multiscale mechanics
• /
• v.2 no.1
• /
• pp.45-68
• /
• 2009

An aim of the study is to develop an efficient numerical simulation technique that can handle the two-scale analysis of fluid permeation filters fabricated by the partial sintering technique of small spherical ceramics. A solid-fluid mixture homogenization method is introduced to predict the mechanical characters such as rigidity and permeability of the porous ceramic filters from the micro-scale geometry and configuration of partially-sintered particles. An extended finite element (X-FE) discretization technique based on the enriched interpolations of respective characteristic functions at fluid-solid interfaces is proposed for the non-interface-fitted mesh solution of the micro-scale analysis that needs non-slip condition at the interface between solid and fluid phases of the unit cell. The homogenization and localization performances of the proposed method are shown in a typical two-dimensional benchmark problem whose model has a hole in center. Three-dimensional applications to the body-centered cubic (BCC) and face-centered cubic (FCC) unit cell models are also shown in the paper. The 3D application is prepared toward the computer-aided optimal design of ceramic filters. The accuracy and stability of the X-FEM based method are comparable to those of the standard interface-fitted FEM, and are superior to those of the voxel type FEM that is often used in such complex micro geometry cases.

• Water for future
• /
• v.27 no.3
• /
• pp.107-113
• /
• 1994

The analysis of Richards eq. requires data of the soil water retention function and the unsaturated hydraulic conductivity. The soil water retention function was measured through the use of the developed apparatus and the saturated hydraulic conductivity was measured by the constant head method for each soil sample corresponding to the A, B, C types of SCS. In order to obtain one water retention function and one unsaturated hydraulic conductivity which represent each soil group, van Genuchten's eq. and Mualem's pore-structure model was chosen respectively. Parameters of van Genuchten's eq. are estimated for each soil group using data obtained in the experiments, and estimated values give a basis to analyze the unsaturated flow in the non-measured region efficiently.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.7 no.4
• /
• pp.41-50
• /
• 1987

The present study was conducted to measure and predict the leachate generation and to establish the methods of leachate control and handling by both field and lab-lysimeter studies. The change of biodegradable matter, field capacity, dry density, quantity and quality of leachate as a function of time of landilling were measured. The model based on the theory of unsaturated flow and contaminant transport in porous media was developed and simulation model was used for the prediction of movement of landfill moisture and for the strength in leachate in terms of organic and inorganic contaminant. The model constructed from the present study in terms of unsaturated decomposition and transport combined with the evapotranspiration and surface runoff showed good agreement with the results obtained by lysimeter studies.

• 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 Wetlands Research
• /
• v.14 no.1
• /
• pp.75-87
• /
• 2012

Hyporheic zone is a region beneath and alongside a stream, river, or lake bed, where there is mixing of shallow groundwater and surfacewater. Hyporheic exchange controls a variety of physical, biogeochemical and thermal processes, and provides unique ecotones in a aquatic ecosystem. Field and experimental observations, and modeling studies indicate that hyporheic exchange is mainly in response to pressure gradients driven by the geomorphological features of stream beds. In the reach scale of a stream, pool-riffle structures dominate the exchange patterns. Flow over a pool-riffle sequence develops recirculation zones and stagnation points, and this flow structures make irregular pressure gradient which is driving force of the hyporheic exchange. In this study, 3 D hydro-dynamic model solves the Reynolds-averaged Navier-Stokes equations for the surface water and Darcy's Law and the continuity equation for ground water. The two sets of equations are coupled via the pressure distribution along the interface. Simulation results show that recirculation zones and stagnation points in the pool-riffle structures dominantly control the upwelling and downwelling patterns. With decrease of recirculation zones, length of donwelling zone formed in front of riffles is reduced and position of maximum downwelling point moves downward. The numerical simulation could successfully predict the behavior of hyporheic exchange and contribute the field study, river management and restoration.

• Journal of the Korean Geotechnical Society
• /
• v.36 no.11
• /
• pp.149-156
• /
• 2020

Permeation grouting effectively enhances soil strength and decreases permeability of soil; however, the flow of grout is heavily affected by anisotropy of hydraulic conductivity in layers. Therefore, this study investigates the effect of permeability anisotropy on the effective radius of horizontal permeation grout using computational fluid dynamics (CFD). We modeled the horizontal permeation grout flow as a two-phase viscous fluid flow in porous media, and the model incorporated the chemical diffusion and the viscosity variation due to hardening. The numerical simulation reveals that the permeability anisotropy shapes the grout bulb to be elliptic and the dissolution-driven diffusion causes a gradual change in grout pore saturation at the edge of the grout bulb. For the grout pore saturations of 10%, 50% and 90%, the horizontal and vertical radii of grout bulb are estimated when the horizontal-to-vertical permeability ratio varies from 0.01 to 100, and the predictive model equations are suggested. This result contributes to more efficient design of injection strategy in formation layers with permeability anisotropy.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.7 no.3
• /
• pp.111-119
• /
• 1987

Scale effects of stability, run-up, run-down and reflection of layered coastal structures are investigated through the experiments with 7 kinds of hydraulic scale models. The occurrence mechanism and the control method of scale effects are also discussed. As a result, it is found that the similarity of permeability of inner layers plays an important role in the occurrence of scale effects, which has been neglected in the most of conventional model tests. To assure the best scale effects for permeable coastal structures, control of Reynolds numbers of the porous media flow in each layer is recommended. It is also found that Reynolds numbers in revetment, filter, and core layer must be greater than $2{\times}10^4$, $3{\times}10^3$, and $1{\times}10^3$, respectively.

• Geophysics and Geophysical Exploration
• /
• v.18 no.4
• /
• pp.223-231
• /
• 2015

Self-potential (SP) is sensitive to groundwater flow and there are many causes to generate SP. Among many mechanisms of SP, pore-fluid flow in porous media can generate potential without any external current source, which is referred to as electrokinetic potential or streaming potential. When calculating SP responses on the surface due to geothermal fluid within an engineered geothermal system (EGS) reservoir, SP anomaly is usually considered to be generated by fluid injection or production within the reservoir. However, SP anomaly can also result from geothermal water fluid within EGS reservoirs experiencing temperature changes between injection and production wells. For more precise simulation of SP responses, we developed an algorithm being able to take account of SP anomalies produced by not only water injection and production but also the fluid of geothermal water, based on three-dimensional finite-element-method employing tetrahedron elements; the developed algorithm can simulate electrical potential responses by both point source and volume source. After verifying the developed algorithm, we assumed a simple geothermal reservoir model and analyzed SP responses caused by geothermal water injection and production. We are going to further analyze SP responses for geothermal water in the presence of water production and injection, considering temperature distribution and geothermal water flow in the following research.

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

최근 전 세계적으로 심각해지는 물 부족 현상과 수질오염으로 대량의 원수를 빠른 시간 내에 여과하기 위한 여과장치의 개발 및 효율성 향상을 위한 연구의 필요성이 증가되고 있다. 특히 여과필터의 내부구조에 의해 유동이 편중되는 현상이 발생하면 여과효율 및 여과필터 유지관리에 문제가 발생되기 때문에 최적의 여과필터를 설계하는 것이 중요하다. 이에 본 연구에서는 수리모형실험으로 검토하기에 어려움이 있는 여과기 내부구조에 대한 유동특성을 수치해석을 이용하여 검토하였다. 수치해석은 유한요소법 기반의 수치모형으로 여과기 내부를 상세하게 모의할 수 없기 때문에 유한체적법 기반인 ANSYS CFX 모형을 이용하였다. 여과기 내 여과필터는 두께 2.0 mm, 공극율 25%로 가정하고 다공성 기법(porous media)을 적용하였다. 검토를 위한 경계조건은 유입부에 목표 취수량, 유출부에 대기압 조건을 적용하였으며, 여과기에 비해 매우 작게 구성된 여과필터 내부의 유동특성을 검토하기 위해 여과기는 최소 3.0 mm, 여과필터는 1.0 mm의 격자를 적용하였다. 현재 실제 여과시설에 적용되고 있는 여과기 제품 형상을 기준으로 여과기 내부 흐름공간의 크기 및 각도 조정에 따른 유동특성을 검토하여 여과효율을 비교하였으며 통과유량, 유속, 유속벡터 등을 검토하여 균등한 유량과 유속이 발생되는 최적의 여과장치 구조를 도출하였다. 본 연구에서 여과기 내부 구조 변경에 따른 유동특성 검토를 통해 도출된 최적의 여과기 내부크기 및 각도에 대한 설계인자는 여과기 내 여과필터의 효율을 증가시킬 뿐만 아니라 내구성 증진에 도움이 될 것으로 예상된다.