• Journal of Advanced Marine Engineering and Technology
• /
• v.20 no.5
• /
• pp.68-75
• /
• 1996

The vacuum drying characteristics of water-saturated porous media were studied experimentally. The water-saturated porous media, water-saturated sand layer, was heated by the isothermal bottom wall of the rectangular vessel. The vacuum drying rate and temperature distribution of the sand layer were measured and calculated under a variety of conditions of heated wall temperature, vacuum rate, and thickness of the test material. It was found that the drying rate due to the heat and mass teansfer is greatly influenced by the heated wall temperature, vacuum rate, and thickness of the test material.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.3 no.5
• /
• pp.386-394
• /
• 1991

Freezing of saturated porous media contained in a rectangular cavity has been studied experimentally. Water and different diameter glass beads consitituted the liquid and porous media. Solidification front shape, the effects of bead diameter and initial liquid temperature was investigated. When the hot wall temperature was below $4^{\circ}C$, the freezing rate was higher at the top than at the bottom due to the density inversion, but with increasing the hot wall temperature the freezing rate at the top was effected by the liquid temperature and was lower than at the bottom. With increasing the bead diameter, the difference of freezing rate between top and bottom was increased and depends on thermal conductivity. When the liquid temperature was low in the beginning, the freezing rate was high, but with increasing the time almost the same with those of high temperature liquid.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.20 no.2
• /
• pp.95-103
• /
• 2007

An analytical closed-form solution for wave propagation velocity and damping in saturated porous media is presented in this paper The fully coupled field model with compressible solid Brains and pore water were used to derive this solution. An engineering approach for the analysis of fully saturated porous media was adopted and closed-form solutions for one dimensional wave propagation in a homogeneous domain were derived. The solution is highly versatile in that it considers compression of the solid grains, compression of the pore water, deformation of the porous skeleton, and spatial damping and can be used to compute wavespeeds of first and second kind and damping coefficients in various geologic materials. This solution provides a means of analyzing the influence of material property variations on wavespeed and attenuation. In Part 2 of this work the theoretical solution is incorporated into the numerical code and the code is used in a parametric study on wave propagation velocity and damping.

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

• Coupled systems mechanics
• /
• v.2 no.4
• /
• pp.375-387
• /
• 2013

Soils consist of an assemblage of particles with different sizes and shapes which form a skeleton whose voids are filled with water and air. Hence, soil behaviour must be analyzed by incorporating the effects of the transient flow of the pore-fluid through the voids, and therefore requires a two-phase continuum formulation for saturated porous media. The present paper presents briefly the Biot's basic theory of dynamics of saturated porous media with u-P formulation to determine the responses of pore fluid and soil skeleton during cyclic loading. Kelvin elements are attached to transmitting boundary. The Pastor-Zienkiewicz-Chan model has been used to describe the inelastic behavior of soils under isotropic cyclic loadings. Newmark-Beta method is employed to discretize the time domain. The response of fluid-saturated porous media which are subjected to time dependent loads has been simulated numerically to predict the liquefaction potential of a semi-infinite saturated sandy layer using finite-infinite elements. A settlement of 17.1 cm is observed at top surface. It is also noticed that liquefaction occurs at shallow depth. The mathematical advantage of the coupled finite element analysis is that the excess pore pressure and displacement can be evaluated simultaneously without using any empirical relationship.

• Coupled systems mechanics
• /
• v.4 no.1
• /
• pp.67-98
• /
• 2015

Numerical prediction of dynamic behavior of fully coupled saturated porous media is of great importance in many engineering problems. Specifically, static and dynamic response of soils - porous media with pores filled with fluid, such as air, water, etc. - can only be modeled properly using fully coupled approaches. Modeling and simulation of static and dynamic behavior of soils require significant Verification and Validation (V&V) procedures in order to build credibility and increase confidence in numerical results. By definition, Verification is essentially a mathematics issue and it provides evidence that the model is solved correctly, while Validation, being a physics issue, provides evidence that the right model is solved. This paper focuses on Verification procedure for fully coupled modeling and simulation of porous media. Therefore, a complete Solution Verification suite has been developed consisting of analytical solutions for both static and dynamic problems of porous media, in time domain. Verification for fully coupled modeling and simulation of porous media has been performed through comparison of the numerical solutions with the analytical ones. Modeling and simulation is based on the so called, u-p-U formulation. Of particular interest are numerical dispersion effects which determine the level of numerical accuracy. These effects are investigated in detail, in an effort to suggest a compromise between numerical error and computational cost.

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

• The Journal of Engineering Geology
• /
• v.14 no.2
• /
• pp.179-187
• /
• 2004

This study conducted a laboratory work to determine the change of dielectric constant of saturated soil porous media, which is injecting EML to use Frequency Domain Reflectometry(FDR) system and length 7cm-length measurement sensor. It is purpose of study that estimate a movement path through pore of soil for an ethanol mixing liquid(EML) which have the same specific gravity of water at $20{\;}^{\circ}C$, and determine to a dielectric measuring range for the measurement sensor. As an outflow point of saturated soil column upper part recedes from an EML inject point in EML diffusion test, the diffusion extent of EML through pore of saturated soil was expanded. And results of ail EML flow experiment were measured the change of dielectric constant at all measurement sensors which are placed to constant interval at the same travel time for saturated soil column. Therefore, the displace process of water that exist in pore of saturated soil by EML injection is enough available and confirm together mobility through pores.

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