• Journal of Advanced Marine Engineering and Technology
• /
• 제20권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.

• 설비공학논문집
• /
• 제3권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.

• 한국전산구조공학회논문집
• /
• 제20권2호
• /
• pp.95-103
• /
• 2007

본 논문에서는 포화된 다공성매체에서 파동의 전파속도와 감쇠를 구할 수 있는 해석적 이론해를 유도하여 제시하였다. 이론해의 유도를 위하여 압축성의 고체입자와 간극수를 고려하는 완전 연계 Field모델을 사용하였다. 완전 포화된 다공성 매체의 해석을 위한 공학적인 접근방법이 채택되었으며, 균질 영역에서 1차원 파동의 전파를 위한 이론해가 유도되었다. 본 논문에서 유도한 이론해는 고체입자의 압축성, 간극수의 압축성, 다공성입자의 변형, 공간의 감쇠(Spatial damping) 등을 고려할 수 있어 매우 다양하게 사용될 수 있다. 또한 다양한 지반체에서 두 가지 종류의 파속(Wavespeed)과 감쇠계수를 계산하는데 이용 가능하다. 본 논문에서 제시한 이론해를 전산코드화하여 파동의 전파속도와 감쇠에 대한 파라미터연구를 수행한 결과는 본 연구의 II부에 제시할 예정이다.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• 제20권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
• /
• 제2권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
• /
• 제4권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
• /
• 제19권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
• /
• 제3권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.

• 지질공학
• /
• 제14권2호
• /
• pp.179-187
• /
• 2004

Frequency Domain Reflectometry(FDR) system과 길이 7cm의 측정센서를 사용하여 에탄을 믹싱리퀴드(EML)를 주입함으로씨 포화된 흙의 유전율상수 변화를 파악하기 위한 실내실험을 실시하였다. 측정센서에 대한 유전율 측정 범위의 확인 및 물과 동일한 비중을 갖는 EML이 포화된 흙의 공극을 통하여 이동 경로를 파악하기 위함이다. EML 확산실험에서는 포화 흙 칼럼 상단부의 배출구가 EML 주입구로부터 이격될수록 포화 흙의 공극을 통한 EML의 확산범위가 확대됨을 확인하였다. 그리고 EML 유동실험의 결과에서는 포화된 흙 칼럼에 대한 주입과 동시에 일정한 간격으로 설치된 모든 측정센서에서 유전율의 변화를 파악하였다. 따라서 EML 주입에 따른 포화 흑의 공극 내에 존재하는 물의 치환이 충분히 가능하며 공극을 통한 이동성도 함께 확인하였다.

• 한국지반공학회:학술대회논문집
• /
• 한국지반공학회 2005년도 춘계 학술발표회 논문집
• /
• 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.