• Water Engineering Research
• /
• v.1 no.3
• /
• pp.243-256
• /
• 2000

This physically-based hydrologic model is developed to calculate the soil-moisture balance on paddy fields. This model consists of three modules; the first is the unsaturated module, the second is the rice evapotranspiration module with SPAC(soil-plant-atmospheric-continuum), and the third is the groundwater and open channel flows based upon the interrehtionship module. The model simulates the hydrlogical processes of infiltration, soil water storage, deep perocolation or echarge to the shallow water table, transpiration and evaporation from the soil surface and also the interrelationship of the groundwater and river flow exchange. To verify the applicability of the developed model, it was applied to the Kimjae Plains, located in the center of the Dongjin river basin in Korea, during the most serious drought season of 1994. The result shows that the estimated water net requirement was 757mm and the water deficit was about 5.9% in this area in 1994. This model can easily evaluate the irrigated water quantity and visualize the common crop demands and soil moisture conditions.

• Proceedings of the Korea Water Resources Association Conference
• /
• 1988.07a
• /
• pp.85-92
• /
• 1988

A model for the transient one-dimensional moisture movement in the saturated-unsaturated zone using a finite difference method is developed. Hysteresis in the soil water retention is incorporated. The model considers layered geologic formations. Monte Carlo simulation, together with the nearest neighbor model is used. Outputs of the model include pressure head, water content, and the water table elevation. Two Monte Carlo simulations of 100 realizations each are made for a 12-day simulation period with different input values. The simulation results show that the S.D. of the outputs increases with an increase in the input, the S.D. of the log K$$. The model is applied to predict a long term water table fluctuation, and the predicted water table agress well with the observed one.

• Journal of Soil and Groundwater Environment
• /
• v.10 no.3
• /
• pp.9-15
• /
• 2005

Soil aquifer treatment is a water reuse technology that secondary or tertiary treated wastewater is infiltrated into the aquifer in which physical and biochemical reactions occur. Major consideration in SAT is the removal and transport of DOC and nitrogen species. In this study, reaction mechanism in SAT was examined considering nitrification, denitrification and organic oxidation. In addition, SAT modeling system was developed as the reaction mechanism was applied to groundwater flow and transport model. In verification of the reaction module by 1-dimensional unsaturated soil column test, the experimental data of all of the species, ammonium, nitrate, DOC and DO, were well matched with the simulation results. In sensitivity analysis, ammonium partition coefficient, dissolved oxygen inhibition constant and biomass decay rate affect ammonium, DOC and DO concentration of effluent, respectively.

• Journal of Soil and Groundwater Environment
• /
• v.27 no.spc
• /
• pp.75-91
• /
• 2022

Differences in subsurface migration of LNAPL/DNAPL contaminants caused by a selection of 3-phase (aqueous, NAPL, and gas) relative permeability function (RPF) models in numerical modeling were investigated. Several types of RPF models developed from both experimental and theoretical backgrounds were introduced prior to conducting numerical modeling. Among the RPF models, two representative models (Stone I and Parker model) were employed to simulate subsurface LNAPLs/DNAPLs migration through numerical calculation. For each model, the spatiotemporal distribution of individual phases and the mole fractions of 6 NAPL components (4 LNAPL and 2 DNAPL components) were calculated through a multi-phase and multi-component numerical simulator. The simulation results indicated that both spilled LNAPLs and DNAPLs in the unsaturated zone migrated faster and reached the groundwater table sooner for Stone I model than Parker model while LNAPLs migrated faster on the groundwater table under Parker model. This results signified the crucial effect of 3-phase relative permeability on the prediction of NAPL contamination and suggested that RPF models should be carefully selected based on adequate verification processes for proper implementation of numerical models.

• Geomechanics and Engineering
• /
• v.23 no.2
• /
• pp.139-149
• /
• 2020

Piles are widely used in structural foundations of engineering projects. However, the deformation of the soil around the pile caused by driving process has an adverse effect on adjacent existing underground buildings. Many previous studies have addressed related problems in sand and saturated clay. Nevertheless, the failure mechanism of pile driving in unsaturated soil remains scarcely reported, and this issue needs to be studied. In this study, a modeling test system based on particle image velocimetry (PIV) was developed for studying deformation characteristics of pile driving in unsaturated silt with different water contents. Meanwhile, a series of direct shear tests and soil-water characteristic curve (SWCC) tests also were conducted. The test results show that the displacement field shows an apparent squeezing effect under the pile end. The installation pressure and displacement field characteristics are sensitive to the water content. The installation pressure is the largest and the total displacement field is the smallest, for specimens compacted at water content of 11.5%. These observations can be reasonably interpreted according to the relevant unsaturated silt theory derived from SWCC tests and direct shear tests. The variation characteristics of the soil displacement field reflect the macroscopic mechanical properties of the soil around the pile.

• Journal of the Korean Geotechnical Society
• /
• v.29 no.1
• /
• pp.93-107
• /
• 2013

Thermal conduction of the particulate composites or granular materials can be widely used in porous materials and geotechnical engineering. And it has continued to develop "effective thermal conductivity" of medium by modeling energy relationship among particles in medium. This study focuses on the development of the effective thermal conductivity at the unsaturated conditions of soils using the modified network model approach assisted by synthetic 3D random packed systems (DEM method, Discrete Element Method) at the particle scale. To verify the network model, three kinds of glass beads and the Jumunjin sand are used to obtain experimental values at various unsaturated conditions. The PPE (Pressure Plate Extractor) test is then performed to obtain SWCC (Soil-Water Characteristic Curve) of soil samples. In the modified network model, SWCC is used to adjust the equivalent radius of thermal cylinder at contact area between particles. And cutoff range parameter to define the effective zone is also adjusted according to the SWCC at given conditions. From a series of laboratory tests and the proposed network model, the modified network model which adopts a SWCC shows a good agreement in modeling thermal conductivity of granular soils at given conditions. And an empirical correlation between the fraction of the mean radius (${\chi}$) and thermal conductivity at given saturated condition is provided, which can be used to expect thermal conductivity of the granular soils, to estimate thermal conductivity of granular soils.

• Computers and Concrete
• /
• v.26 no.2
• /
• pp.127-136
• /
• 2020

There is an inherent randomness for concrete microstructure even with the same manufacturing process. Meanwhile, the concrete material under the aqueous environment is usually not fully saturated by water. This study aimed to develop a stochastic micromechanical framework to investigate the probabilistic behavior of the unsaturated concrete from microscale level. The material is represented as a multiphase composite composed of the water, the pores and the intrinsic concrete (made up by the mortar, the coarse aggregates and their interfaces). The differential scheme based two-level micromechanical homogenization scheme is presented to quantitatively predict the concrete's effective properties. By modeling the volume fractions and properties of the constituents as stochastic, we extend the deterministic framework to stochastic to incorporate the material's inherent randomness. Monte Carlo simulations are adopted to reach the different order moments of the effective properties. A distribution-free method is employed to get the unbiased probability density function based on the maximum entropy principle. Numerical examples including limited experimental validations, comparisons with existing micromechanical models, commonly used probability density functions and the direct Monte Carlo simulations indicate that the proposed models provide an accurate and computationally efficient framework in characterizing the material's effective properties. Finally, the effects of the saturation degrees and the pore shapes on the concrete macroscopic probabilistic behaviors are investigated based on our proposed stochastic micromechanical framework.

• Environmental Engineering Research
• /
• v.11 no.5
• /
• pp.241-249
• /
• 2006

A mathematical model is described for the prediction of convective upward transport of an organic solvent driven by evaporation at the surface, which is known as the major transport mechanism in the in-situ photolysis of a soil contaminated with 2,3,7,8-tetrachlorodibenzo-p-dioxin(TCDD). A finite-element model was proposed to incorporate the effects of multiphase flow on the distribution of each fluid, gravity as a driving force, and the use of hysteretic models for more accurate description of k-S-p relations. Extensive numerical calculations were performed to study fluid flow through three types of soils under different water table conditions. Predictions of relative permeability-saturation-pressure (k-S-p) relations and fluids distribution for an illustrative soil indicate that hysteresis effects may be quite substantial. This result emphasizes the need to use hysteretic models in performing flow simulations including reversals of flow paths. Results of additional calculations accounting for hysteresis on the one-dimensional unsaturated soil columns show that gravity affects significantly on the flow of each fluid during gravity drainage, solvent injection, and evaporation, especially for highly permeable soils. The rate and duration of solvent injection also have a profound influence on the fluid saturation profile and the amount of evaporated solvent. Key factors influencing water drainage and solvent evaporation in soils also include hydraulic conductivity and water table configuration.

• Journal of the Korean Geotechnical Society
• /
• v.19 no.5
• /
• pp.99-107
• /
• 2003

Moisture flow and heat flow within pavement systems have been recognized as coupled processes with complex interactions between them. The distribution of moisture and temperature within pavement due to the moisture flow and heat flow varies not only seasonally but also vertically and horizontally. This paper presents an analysis model by the finite element method for the two-dimensional coupled moisture and heat flow in unsaturated soils. To test the model the analysis result by the model is compared with the analysis result by the software, GEO-SLOPE developed by GEO-SLOPE International Ltd. in Alberta, Canada. And a sensitivity analysis using ASTM method is performed to identify how model inputs affect the modeling analysis.

• Proceedings of the Korean Society of Soil and Groundwater Environment Conference
• /
• 2003.09a
• /
• pp.140-143
• /
• 2003

Air permeabilities in unsaturated fractured tuff at the Apache Leap Research Site (ALRS) near Superior, Arizona, exhibit a self-affine behavior, thus renders a field random fractal. Based up fractal scaling, the observed scale effect has been interpreted [Hyun et al., 2002]. Recently, Frantziskonis and Hansen [2000] presented that fractal scaling can be represented based on wavelets. This study deals with the way of using wavelets for fractal scaling. A numerical study is presented to examine the applicability of wavelet-based approach to determining upscaled air permeability values on various data supports at the site. To characterize the scaling property of self-affine fields generated based upon wavelets, Hurst coefficient, H. was inferred by applying the average wavelet coefficient (AWC) method. The result yielded H = 0.24, which is very close to the result of geostatistical analysis using a power variogram (H = 0.22). The study concludes that wavelet-based scaling is a useful way of determining parameter values on different data supports, which is an essential task for modeling of subsurface flow and mass transport in a numeric grid with different resolutions (grid size).