• Journal of The Korean Society of Agricultural Engineers
• /
• v.48 no.5
• /
• pp.83-93
• /
• 2006

In unsaturated porous media, the air-water interface (AWI) plays an important role in removing of biocolloids such as bacteria, viruses, and protozoan (oo)cysts. In this study, four models related to calculation of specific AWI area are analyzed to determine the appropriate model, and the selected models are verified using the previously reported experimental data. The results indicate that the modified model from Niemet et al. (2002) is the most appropriate tool for calculating the specific AWI area using the van Genuchten (1980) parameters obtained from the water retention curve. Hence, it is expected that this model could be used to quantitatively determine the attachment of biocolloids to AWI in the transport modeling of biocolloids in unsaturated porous media.

• Ocean Systems Engineering
• /
• v.12 no.1
• /
• pp.63-86
• /
• 2022

A refined projection-based purely Lagrangian meshfree method is presented towards reliable numerical analysis of fluid flow interactions with saturated/unsaturated porous media of uniform/spatially-varying porosities. The governing equations are reformulated on the basis of two-phase mixture theory with incorporation of volume fraction. These principal equations of mixture are discretized in the context of Incompressible SPH (Smoothed Particle Hydrodynamics) method. Associated with the consideration of governing equations of mixture, a new term arises in the source term of PPE (Poisson Pressure Equation), resulting in modified source term. The linear and nonlinear force terms are included in momentum equation to represent the resistance from porous media. Volume increase of fluid particles are taken into consideration on account of the presence of porous media, and hence multi-resolution ISPH framework is also incorporated. The stability and accuracy of the proposed method are thoroughly examined by reproducing several numerical examples including the interactions between fluid flow and saturated/unsaturated porous media of uniform/spatially-varying porosities. The method shows continuous pressure field, smooth variations of particle volumes and regular distributions of particles at the interface between fluid and porous media.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2006.04a
• /
• pp.381-388
• /
• 2006

A porous medium is composed of solids, fluids, and gas which have different physical and chemical properties. In addition, these constituents have a relative velocity between each other. So far, in order to analyze porous media using finite element method, Lagrangian or Eulerian method has been used. However, the numerical analyses for porous media have a defect that the methods do not describe the movements of constituents. In this paper, numerical analysis for unsaturated porous media was performed in frame of ALE method which has advantages of Lagrangian and Eulerian. Namely, the Lagrangian description was used in solid phase, and the Eulerian description was used in fluid or gas phase in a porous medium Then the relationship between each other was controlled by the convective term in ALE method. Finally, the numerical results of ALE were compared with tile results of Lagrangian analysis.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.20 no.3
• /
• pp.243-259
• /
• 2016

The Richards equation for water movement in unsaturated soil is highly nonlinear partial differential equations which are not solvable analytically unless unrealistic and oversimplifying assumptions are made regarding the attributes, dynamics, and properties of the physical systems. Therefore, conventionally, numerical solutions are the only feasible procedures to model flow in partially saturated porous media. The standard Finite element numerical technique is usually coupled with an Euler time discretizations scheme. Except for the fully explicit forward method, any other Euler time-marching algorithm generates nonlinear algebraic equations which should be solved using iterative procedures such as Newton and Picard iterations. In this study, lumped mass and distributed mass in the frame of Picard and Newton iterative techniques were evaluated to determine the most efficient method to solve the Richards equation with finite element model. The accuracy and computational efficiency of the scheme and of the Picard and Newton models are assessed for three test problems simulating one-dimensional flow processes in unsaturated porous media. Results demonstrated that, the conventional mass distributed finite element method suffers from numerical oscillations at the wetting front, especially for very dry initial conditions. Even though small mesh sizes are applied for all the test problems, it is shown that the traditional mass-distributed scheme can still generate an incorrect response due to the highly nonlinear properties of water flow in unsaturated soil and cause numerical oscillation. On the other hand, non oscillatory solutions are obtained and non-physics solutions for these problems are evaded by using the mass-lumped finite element method.

• Proceedings of the Korean Geotechical Society Conference
• /
• 2010.09a
• /
• pp.808-812
• /
• 2010

A great deal of attention is focused on coupled Thermo-Hydro-Mechanical (THM) behavior of multiphase porous media in diverse geo-mechanical and geo-environmental areas. This paper presents general governing equations for coupled THM processes in unsaturated porous media. Coupled partial differential equations are derived from 3 mass balances equations (solid, water, and air), energy balance equation, and force equilibrium equation. Finite element code is developed from the Galerkin formulation and time integration of these governing equations for 4 main variables (displacement $\underline{u}$, gas pressure $P_g$, liquid pressure $P_l$), and temperature T). The code is validated with theoretical solutions for linear material with simple boundary conditions.

• The Journal of Engineering Geology
• /
• v.14 no.2
• /
• pp.169-177
• /
• 2004

The permeation movements of groundwater recharge and contaminate materials receive a eat effect due to porosity and effective porosity of porous media which is composing underground consisted of saturation and unsaturated states. This study developed Frequency Domain Reflectometry(FDR) system and measurement sensor, and then carried out the laboratory experiments to measure effective porosity for unsaturated porous media. Also, I suggested dielectric mixing models(DMMs) which can calculate the effective porosity from relation of measured dielectric constants. In the experimental results the extent range of effective porosity of standard sand and river sand which are unsaturated soil sample were measured in about 65∼85 % for porosity. In relation of effective porosity and porosity, especially, effective porosity confirmed that displays decreasing a little tendency as porosity increases. This is because unsaturated soil did not reach in saturation enough by air of very small amount that exist in pore between soil particles.

• Proceedings of the Korean Society of Soil and Groundwater Environment Conference
• /
• 1996.11a
• /
• pp.58-62
• /
• 1996

To check the variability of the effective diffusion coefficient in the unsaturated porous soil media. a Monte Carlo simulation was done for the equation suggested by Millington and Quirk(1961). The results shows that the probability density function of D$_{c}$/D$_{o}$ is positively skewed. It means the chance of having less effective diffusion coefficient values in the soil media than mean value is high. Also, the distribution types of D$_{c}$/D$_{o}$ are about same regardless of assumed distribution types of input parameters.ers.ers.

• Environmental Engineering Research
• /
• v.10 no.2
• /
• pp.71-78
• /
• 2005

A series of column experiments was conducted to develop a monitoring system for in-situ and realtime measurement of ozone transport in unsaturated porous media using a fiber optic sensor. The calibration of the fiber optic transflection dip probe (FOTDP) system was successfully carried out at various ozone concentrations using a column with length of 30 cm and diameter of 5 cm packed with glass beads, which don't react with gaseous ozone. The breakthrough curves (BTCs) of ozone were obtained by converting the normalized intensity into ozone concentration. The FOTDP system worked well for in-situ monitoring of gas phase ozone at various water saturations and in presence of soil organic matter (SOM). However, the FOTDP system did not measure the ozone concentration at more than 70% water saturation.

• Proceedings of the Korean Society of Soil and Groundwater Environment Conference
• /
• 1997.05a
• /
• pp.112-115
• /
• 1997

• Proceedings of the Korean Geotechical Society Conference
• /
• 1991.10a
• /
• pp.235-254
• /
• 1991

In the analysis of seepage behavior through the porous media, unsaturated domain should be included for obtaining more accurate results. The seepage programs considering its domain have been developed and practically used widely. The programs are reviewed and two case studies are presented in this paper. The case studies are performed using PC-SEEP and obtained good results for the cases of saturated-unsaturated transient flow.