• Journal of the Korean Geotechnical Society
• /
• v.30 no.2
• /
• pp.77-85
• /
• 2014

The paper presents comparison between numerical and experimental results to verify finite element algorithm of air diffusion in three-phase porous media. The theory of two-phase flow in partially saturated soil is a fundamental study to predict the safety for $CO_2$ sequestration. Geotechnical engineering problems exposed to atmospheric pressure generally assume that air pressure is equal to zero in three-phase porous media, but the geotechnical engineering problem at a depth of about over 1 km needs to be considered in connection with deformable solid skeleton due to air pressure. Therefore, the objective of this study is to validate the numerical algorithm by comparing with results obtained from measurement of air diffusion and dissipation through drainage test.

• Journal of Korean Society of Environmental Engineers
• /
• v.30 no.10
• /
• pp.984-991
• /
• 2008

This paper reports experimental results, demonstrating the feasibility of horseradish peroxidase(HRP) and H$_2$O$_2$ to reduce phenol transport in saturated porous media. A laboratory-scale packed column reactor(ID: 4.1 cm, sand-bed height 12 cm) column was utilized to simulate injection of HRP and H$_2$O$_2$ into an aquifer contaminated with phenol. Effluent concentrations of phenol and polymerization products were monitored before and after enzyme addition under various experimental conditions(enzyme dose: 0$\sim$2 AU/mL, [ionic strength]: 5$\sim$100 mM, pH: 5$\sim$9). The concentration of phenol in the column effluent was found to decrease by nearly 90% in the presence of HRP(2 AU/mL) and H$_2$O$_2$ in the continuous flow system at pH 7 and ionic strength 20 mM. The influent phenol was converted in the system to insoluble precipitate, which deposited in pore spaces. The remains were discharged as soluble oligomers. About 8% of total pore volume in column system was decreased by deposition of polymer produced.

• Computers and Concrete
• /
• v.4 no.2
• /
• pp.157-166
• /
• 2007

The paper considers a theoretical model to study sulfate ion diffusion in saturated porous media - cement based mineral composites, accounting for simultaneous effects, such as filling micro-capillaries (pores) with ions and chemical products and liquid push out of them. Pore volume change and its effect on the distribution of ion concentration within the specimen are investigated. Relations for the distribution of the capillary relative radius and volume within the composite under consideration are found. The numerical algorithm used is further completed to consider capillary size change and the effects accompanying sulfate ion diffusion. Ion distribution within the cross section and volume of specimens fabricated from mineral composites is numerically studied, accounting for the change of material capillary size and volume. Characteristic cases of 2D and 3D diffusion are analyzed. The results found can be used to both assess the sulfate corrosion in saturated systems and predict changes occurring in the pore structure of the composite as a result of sulfate ion diffusion.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.13 no.5
• /
• pp.109-122
• /
• 1993

ln this paper, a finite element analysis procedure is proposed for the incremental multi-step excavations in a fluid-saturated porous medium such as saturated soil ground. As the basis of derivation, Biot's equation was used. The proposed procedure was applied to some one- and two-dimensional problems under incremental excavations. Unsaturated cases as well as saturated cases were considered for comparison. Through numerical tests, the effects of permeability and excavation speed on the deformation history was investigated. Results showed that pore pressure built up during incremental excavation has a significant effect on the deformation and stresses of solid skeleton and validated the use of the present procedure for the analysis of multi-step excavations in fluid-saturated media such as in saturated shallow ground.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.19 no.2
• /
• pp.197-215
• /
• 2015

Modeling water flow in variably saturated, porous media is important in many branches of science and engineering. Highly nonlinear relationships between water content and hydraulic conductivity and soil-water pressure result in very steep wetting fronts causing numerical problems. These include poor efficiency when modeling water infiltration into very dry porous media, and numerical oscillation near a steep wetting front. A one-dimensional finite element formulation is developed for the numerical simulation of variably saturated flow systems. First order backward Euler implicit and second order Crank-Nicolson time discretization schemes are adopted as a solution strategy in this formulation based on Picard and Newton iterative techniques. Five examples are used to investigate the numerical performance of two approaches and the different factors are highlighted that can affect their convergence and efficiency. The first test case deals with sharp moisture front that infiltrates into the soil column. It shows the capability of providing a mass-conservative behavior. Saturated conditions are not developed in the second test case. Involving of dry initial condition and steep wetting front are the main numerical complexity of the third test example. Fourth test case is a rapid infiltration of water from the surface, followed by a period of redistribution of the water due to the dynamic boundary condition. The last one-dimensional test case involves flow into a layered soil with variable initial conditions. The numerical results indicate that the Crank-Nicolson scheme is inefficient compared to fully implicit backward Euler scheme for the layered soil problem but offers same accuracy for the other homogeneous soil cases.

• Advances in environmental research
• /
• v.3 no.1
• /
• pp.45-69
• /
• 2014

The groundwater has been a major source of water supply throughout the ages. Around 50% of the rural as well as urban population in the developing countries like India depends on groundwater for drinking. The groundwater is also an important source in the agriculture and industrial sector. In many parts of the world, groundwater resources are under increasing threat from growing demands, wasteful use and contamination. A good planning and management practices are needed to face this challenge. A key to the management of groundwater is the ability to model the movement of fluids and contaminants in the subsurface environment. It is obvious that the contaminant source activities cannot be completely eliminated and perhaps our water bodies will continue to serve as receptors of vast quantities of waste. In such a scenario, the goal of water quality protection efforts must necessarily be the control and management of these sources to ensure that released pollutants will be sufficiently attenuated within the region of interest and the quality of water at points of withdrawal is not impaired. In order to understand the behaviour of contaminant transport through different types of media, several researchers are carrying out experimental investigations through laboratory and field studies. Many of them are working on the analytical and numerical studies to simulate the movement of contaminants in soil and groundwater of the contaminant transport. With the advent of high power computers especially, a numerical modelling has gained popularity and is indeed of particular relevance in this regard. This paper provides the state of the art of contaminant transport and reviews the allied research works carried out through experimental investigation or using the analytical solution and numerical method. The review involves the investigation in respect of both, saturated and unsaturated, porous media.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.23 no.4
• /
• pp.361-368
• /
• 2010

In this paper, the parallel algorithm using MPI(Message-Passing Interface) library is introduced in order to improve numerical efficiency for the staggered method for nearly incompressible and impermeable porous media which was introduced by Park and Tak(2010). The porous media theory and the staggered method are also briefly introduced in this paper. Moreover, we account for MPI library for blocking, non-blocking, and collective communication, and propose combined the staggered method with the blocking and nonblocking MPI library. And then, we present how to allocate CPUs on the staggered method and the MPI library, which is related with the numerical efficiency in order to solve unknown variables on nearly incompressible and impermeable porous media. Finally, the results comparing serial solution with parallel solution are verified by 2 dimensional saturated porous model according to the number of FEM meshes.

• Computers and Concrete
• /
• v.4 no.4
• /
• pp.273-284
• /
• 2007

The paper considers a theoretical model for the study of the process of transfer of sulfate ions in saturated porous media - mineral composites. In its turn, the model treats diffusion of sulfate ions into cement based composites, accounting for simultaneous effects such as filling of micro-capillaries with ions and chemical products and liquid push out of them. The proposed numerical algorithm enables one to account for those simultaneous effects, as well as to model the diffusive behavior of separate sections of the considered volume, such as inert fillers. The cases studied illustrate the capabilities of the proposed model and those of the algorithm developed to study diffusion, considering the specimen complex configuration. Computations show that the theoretical assumptions enable one to qualitatively estimate the experimental evidence and the capabilities of the studied composite. The results found can be used to both assess the sulfate corrosion in saturated systems and predict and estimate damage of structures built of cement-based mineral composites.

• Proceedings of the Korean Geotechical Society Conference
• /
• 2005.03a
• /
• pp.354-361
• /
• 2005

Physical parameters such as porosity and effective porosity are important physical parameters that determine the transfer and movement of water and solutes in porous media. Various methods of determining these parameters have been developed, with varying degrees of accuracy and applicability. Most of the existing methods produce static results. They do not produce instantaneous and real time of porosity and effective porosity in a porous media. In this study, a new permittivity method called Frequency Domain Reflectometry with Vector analyzer (FDR-V) is proposed to determine the porosity and effective porosity of some sand samples in the laboratory. The advantage of the FDR-V method is that it instantaneously determines the temporal variation of dielectric constants of porous media. Then, the porosity and the effective porosity of porous media are computed using well established empirical equations. Results obtained from the FDR-V method compared favorably with results from other permittivity methods such as gravimetric, injection and replacement tests. The ratio of effective porosity to porosity was $85{\sim}92%$, when FDR-V was used. This value compared favourably with 90%, which has been usually quoted in previous studies. Considering the convenience and its applicability, the measurement system of FDR-V permittivity holds a great potential in porous media and contaminant transport studies.

• Geomechanics and Engineering
• /
• v.2 no.2
• /
• pp.141-156
• /
• 2010

In the framework of meshfree methods, a new methodology is developed based on radial point interpolation method (RPIM). This methodology is applied to a one-dimensional contaminant transport modelling in the saturated porous media. The one-dimensional form of advection-dispersion equation involving reactive contaminant is considered in the analysis. The Galerkin weak form of the governing equation is formulated using 1D meshfree shape functions constructed using thin plate spline radial basis functions. MATLAB code is developed to obtain the numerical solution. Numerical examples representing various phenomena, which occur during migration of contaminants, are presented to illustrate the applicability of the proposed method and the results are compared with those obtained from the analytical and finite element solutions. The proposed RPIM has generated results with no oscillations and they are insensitive to Peclet constraints. In order to test the practical applicability and performance of the RPIM, three case studies of contaminant transport through the landfill liners are presented. A good agreement is obtained between the results of the RPIM and the field investigation data.