• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.14 no.5
• /
• pp.1137-1143
• /
• 1994

Through the numerical analysis of the capillary pressure-based Richards equation with and without the effect of the capillary hysteresis under the boundary condition having an antecedent rainfall. the moving tendency of the wetting front, the redistribution of the moisture content, infiltration rate, cumulative infiltration etc, were computed. The effect of the capillary hysteresis cannot be neglected in analyzing an unsaturated flow, and the more accurate results may be obtained by the consideration of the hysteresis effect. If the effect of the hysteresis cannot be considered, the analysis by the use of the main wetting curve may give more reliable result than that of the main drying curve.

• Journal of Korean Society of Disaster and Security
• /
• v.10 no.1
• /
• pp.19-27
• /
• 2017

In recent years, ground subsidence has been frequently occurred by underground cavities due to the excessive groundwater inflow, caused by poor construction and management, during tunnel excavation and underground structure construction. In this study, a numerical model (SEEFLOW3D) was developed to estimate groundwater fluctuations for saturated-unsaturated poros media, evaluates the impact on ground excavation with open cut and non-open cut scenarios. In addition, the visual MODFLOW was applied to demonstrate the verification of the model compared with both results. Our results indicated that the RMSE and NRMSE was obtained to range over -3.95~5.7% and 0.56~4.62%, respectively. The developed model was expected to estimate groundwater discharges and apply analysis tool for optimum design of waterproof wall in future.

• Water Engineering Research
• /
• v.2 no.3
• /
• pp.187-196
• /
• 2001

An approximate analytical solution of a nonlinear hydro-thermo coupled diffusion equation is derived using the dimensionless form of the equation and transformation method. To derive an analytical solution, it is drastically assumed that the product of first order derivatives in the non-dimensionalized governing equation has little influence on the solution of heat and moisture behavior problem. The validity of this drastic assumption is demonstrated. Some numerical simulation is performed to investigate the applicability of a derived approximate analytical solution. The results show a good agreement between analytical and numerical solutions. The proposed solution may provide a useful tool in the verification process of the numerical models. Also, the solution can be used for the analysis of one-dimensional coupled heat and moisture movements in unsaturated porous media.

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

토양내에서 오염유기물질이 불포화토양내에 유입될 때의 dispersion coefficient를 adsorption과 desorption과정에 대해 알아보았다. apparent dispersion coefficient를 측정하기 위해 일상적인 상대습도(46%)조건에서 parametric analysis를 행하였다. 실험에 사용된 토양은 fine sand와 silt-clay혼합시료였고, 흐름방향은 상향과 하향으로 하였다. 그리고, Freon gas를adsorbing solute로 사용하였다. 오염물질로는 DCM, TCE, DCB를 사용하였다. 분석을 위해서 linear와 probability scale의 breakthrough curve를 사용하였다. 공기에서의 diffusion coefficient의 예측을 위하여 Graham's law를 계산에 사용하였고, DCM diffusion coefficient는 0.098$\textrm{cm}^2$/s로 계산되었다. 연구결과, adsorption과 desorption의 속도는 차이가 있는 것으로 나타났으며, diffusion이 flow regime을 좌우하는 것으로 나타났다. 그리고, desorption에서의 D$^{a}$ D$^{o}$ 는 1보다 클수도 있다. 또한, dispersion은 silt-clay혼합시료에서의 속도와 함께 증가한다. dispersion은 Freon의 sorption방향에 크게 의존한다.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.12 no.6
• /
• pp.599-608
• /
• 2000

One dimensional numerical modeling was carried out for the melting behavior of dry snow and the unsaturated flow when heat was supplied from the bottom surface. Discrepancy between the previous experimental data and the present numerical results is substantially reduced by considering the density change of water permeation layer due to the infiltration of meltwater. In the parametric study for effective thermal conductivity, it was found that the effect of this parameter to the behavior of snow melting is minor. Sensitivity analysis showed that the melting time is most sensitive to changes in supplied heat flux, snow temperature, and bulk density, whereas snow bulk density and residual saturation have a significant effect on the height of water permeation layer in snow.

• Journal of Soil and Groundwater Environment
• /
• v.12 no.6
• /
• pp.60-69
• /
• 2007

Understanding the behavior of gas phase VOCs (volatile organic compounds) in unsaturated soils is of a great environmental importance for public health concerns. Moment analysis for the breakthrough curves (BTCs) during transport of chemicals in porous media was known to be a useful tool to evaluate the velocity, spreadness, and the skewness of the plume of the chemicals. In this study, the temporal moments of the BTCs of a group of VOCs were analyzed for the gaseous transport in an unsaturated soil. BTCs were measured using lab-scale column experiments for four different VOCs at the water saturation range of 0.04-0.46, and for eleven VOCs at a water saturation of 0.21. The central second and third moments of the VOCs were compared with the water saturation and the first moment. It was found that both central second and third moments increased with the first moment. The central third moment was, however, found to be more sensitive to the first moment.

• The Journal of Engineering Geology
• /
• v.19 no.2
• /
• pp.205-215
• /
• 2009

The conventional numerical models to analyze flow in subsurface porous media under the transient state usually generate numerical oscillation and unstability due to local flux domain for critical cases such as infiltration into initially dry soil during rainfall period. In this case, it is required refined mesh and small time step, but it decrease efficiency of computation. In this study, numerical unstability in discontinuity domain is removed by applying particle tracking algorithm to simulate unsteady subsurface flow with inflow boundary condition. Finally the hybrid LE FEM improving numerical stability is proposed. The hypothetical domains with unsteady uniform and nonuniform flow field were used to demonstrated algorithm verification. In comparison with analytic solution, we obtained reasonable results and conducted simulation of hypothetical 3-D recharge/pumping area. The proposed algorithm can simulate saturated/unsaturated porous media with more practical problems and will greatly contribute to accuracy and stability of numerical computation.

• 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 Nuclear Fuel Cycle and Waste Technology(JNFCWT)
• /
• v.17 no.3
• /
• pp.329-338
• /
• 2019

Suction of unsaturated bentonite buffers is a very important input parameter for hydro-mechanical performance assessment and design of an engineered barrier system. This study analyzed suction measurement techniques and constitutive models of unsaturated porous media reported in the literature, and suggested suction measurement techniques and constitutive models suitable for bentonite buffer in an HLW repository. The literature review showed the suction of bentonite buffer to be much higher than that of soil, as measured by total suction including matric suction and osmotic suction. The measurement methods (RH-Cell, RH-Cell/Sensor) using a relative humidity sensor were suitable for suction measurement of the bentonite buffer; the RH-Cell /Sensor method was more preferred in consideration of the temperature change due to radioactive decay heat and measurement time. Various water retention models of bentonite buffers have been proposed through experiments, but the van Genuchten model is mainly used as a constitutive model of hydro-mechanical performance assessment of unsaturated buffers. The water characteristic curve of bentonite buffers showed different tendencies according to bentonite type, dry density, temperature, salinity, sample state and hysteresis. Selection of water retention models and determination of model input parameters should consider the effects of these controlling factors so as to improve overall reliability.

• Water for future
• /
• v.28 no.4
• /
• pp.215-222
• /
• 1995

The water retention function which has the hysteresis phenomena is required to analyze the Richards equation which is a governing equation of the unsaturated flow, and its hysteresis phenomena has influence upon the characteristics of the unsaturated flow. The accuracy of the published hysteresis models is compared by using experimental data of the water retention function. The apparatus to experiment the hysteresis phenomena on the soil is developed, and experimental data for the main wetting process and the main drying process of the water retention function are obtained. The parameters of the van Genuchten equation are calibrated by using experimentally obtained data. As a result of the comparison of the selected hysteresis models which simulate the main drying curve from the main wetting curve, the Model I-1(Mualem) overestimates and the Model II-1(Mualem) underestimates but the Model III-2(Park and Sonu) similarly estimates the experimental data of the main drying curve.