• The Journal of Engineering Geology
• /
• v.23 no.3
• /
• pp.227-234
• /
• 2013

The vertical shaft of a selective intake structure, which is constructed in a large reservoir, is required to be impermeable and to employ a grouting technology to prevent water inflow from the reservoir or surrounding ground. In this study, groundwater inflow is estimated using a numerical model for two cases (i.e., grouting or non-grouting cases at the exterior of a vertical shaft) and compared with data measured during an excavation at the construction site of a selective intake structure in the Soyang reservoir, Korea. Groundwater inflow is estimated to range from 444 to 754 $m^3/d$ in the case of non-grouting and from 58 to 95 $m^3/d$ in the case of grouting. The groundwater inflow measured in a vertical shaft, which ranges from 30 to 100 $m^3/d$, is similar to the simulated amount. It is recommended that before the excavation of a shaft, water inflow is estimated using a numerical model and a grouting test to ensure excavation stability and improve excavation efficiency.

• Clean Technology
• /
• v.25 no.3
• /
• pp.206-222
• /
• 2019

Recently, various multiphase flows have been developed, and among them some models have been commercialized. However, most of them have been developed based on a pressure-based approach; therefore, various numerical difficulties were involved inherently. Accordingly, in order to overcome these numerical difficulties, a multiphase flow model, MultiPhaSe flow (MPS), following a fractional-flow based approach was developed. In this study, by combining a contaminant transport module describing an enhanced dissolution effect of a surfactant with MPS, a MultiPhaSe flow and TranSport (MPSTS) model was developed. The developed model was verified using the analytical solution of Clement. The MPSTS model can simulate the process of surfactant enhanced aquifer remediation including interphase mass transfer and contaminant transport in multiphase flow by using the coupled particle tracking method and Lagrangian-Eulerian method. In this study, a surfactant was used in a non aqueous phase liquid (NAPL) contaminated area, and the effect of hydro-geological heterogeneity in the layered media on remediation efficiency was studied using the developed model. According to the numerical simulation, when hydraulic conductivity in a lower layer is 10 times, 20 times, and 50 times larger than that in an upper layer, the concentration of dissolved diesel in the lower layer is much higher than that in the upper layer because the surfactant moves faster along the lower layer owing to preferential flow; thus, the surfactant enhances dissolution of residual non aqueous phase liquid in the lower layer.

• Tunnel and Underground Space
• /
• v.33 no.4
• /
• pp.265-280
• /
• 2023

Water leakage in tunnels is a defect that can affect tunnel stability and the ground movement by changing the stress and pore water pressure of the surrounding ground. Long-term or large-scale water leaks may lead to damage of tunnel structure and the surrounding environment, such as tunnel lining instability and ground surface settlement. The present study numerically investigated the effects of water leakage on the structural stability of a tunnel and the ground behavior. The tunnel was assumed to be under undrained conditions for preventing the inflow of the surrounding water and leaks occurred in the concrete lining after completion of the tunnel construction. A coupled hydro-mechanical analysis using a TOUGH-FLAC simulator developed in Python was conducted for assessing the leakage induced-behavior of the tunnel structure and ground under different conditions of the amount and location of water leak. Additionally, the effect of hydro-mechanical coupling terms on the results of coupled response was investigated and discussed.

• Nuclear Engineering and Technology
• /
• v.24 no.3
• /
• pp.228-235
• /
• 1992

Some soluble species may not be solubility-limited or congruent-released with the matrix species. For example, during the operation of the nuclear reactor, the fission products can be accumulated in the fuel-cladding gap, voids, and grain boundaries of the fuel rods. In the waste package for spent-fuel placed in a geologic repository, the high solubility species of these fission products accumulated in the“gap”, e.g. cesium or iodine are expected to dissolve rapidly when ground water penetrates fuel rods. The time and space dependent mass transport for high solubility nuclides in the gap is analyzed, and its numerical illustrations are demonstrated. The approximate solution that is valid for all times is developed, and validated by comparison with an asymptotic solution and the solution obtained by the numerical inversion of Laplace transform covering the entire time span.

• Journal of Korea Water Resources Association
• /
• v.34 no.5
• /
• pp.449-458
• /
• 2001

In order to verify the numerical model, which was developed to simulate the behavior of the two-phase fluid flow in a single fracture, the characteristic equation of relative permeability was incorporated into the developed numerical model, and the computed results were compared with the experimental results of the model test. As results of the sensitivity analysis on the roughness and the aperture size of fracture, the gas velocity was inversely proportional to the fracture roughness, and not proportional to the square of aperture size which is usually observed in single phase flow in a single fracture. The numerical model was applied to the underground LPG storage terminal in order to check the field applicability. The simultaneous flow of water and gas in accordance with the operation pressures in a single fracture near cavern was simulated by the model. It was shown that the leaked gas was able to be controlled in a single fracture neither by the pressure of operation nor by that of groundwater in case the fracture became smoother in roughness and smaller in aperture size.

• Tunnel and Underground Space
• /
• v.20 no.2
• /
• pp.82-91
• /
• 2010

Grouting technology is one of the ground improvement methods used in water controlling and reinforcement of rock mass in underground structure construction. It is necessarily required to find out the characteristics of grout flow through discontinuities in a rock mass for an adequate grout design and performance assessment. Laminar flow is not always applicable in simulating a grout flow in a rock mass, since the rock joints usually have apertures at a micro-scale and the flow through these joints is affected by the joint roughness and the velocity profile of the flow changes partially near the roughness. Thus, the influence of joint roughness and aperture on the grout flow in rough rock joint was numerically investigated in this study. The commercial computational fluid dynamics code, FLUENT, was applied for this purpose. The computed results by embedded Herschel-Bulkley model and VOF (volume of fluid) model, which are applicable to simulate grout flow in a narrow rock joint that is filled with air and water, were well compared with that of analytical results and previously published laboratory test for the verification. The injection pressure required to keep constant injection rate of grout was calculated in a variety of Joint Roughness Coefficient (JRC) and aperture conditions, and the effect of joint roughness and aperture on grout flow were quantified.

• Journal of the Korean Society of Groundwater Environment
• /
• v.6 no.1
• /
• pp.23-32
• /
• 1999

A new numerical method using solution-adaptive grids (SAG) is developed to solve the Richards' equation (RE) for unsaturated flow in porous media. Using a grid generation technique, the SAG method automatically redistributes a fixed number of grid points during the flow process, so that more grid points are clustered in regions of large solution gradients. The method uses the coordinate transformation technique to employ a new transformed RE, which is solved with the standard finite difference method. The movement of grid points is incorporated into the transformed RE, and therefore all computation is performed on fixed grid points of the transformed domain without using any interpolation techniques. Thus, numerical difficulties arising from the movement of the wetting front during the infiltration process have been substantially overcome by the new method. Numerical experiments for an one-dimensional infiltration problem are presented to compare the SAG method to the modified Picard method using a fixed grid. Results show that accuracy of a SAG solution using 41 nodes is comparable with the solution of the fixed grid method using 201 nodes, while it requires only 50% of the CPU time. The global mass balance and the convergence of SAG solutions are strongly affected by the time step size (Δt) and the weighting parameter (${\gamma}$) used for generating solution-adaptive grids. Thus, the method requires automated readjustment of Δt and ${\gamma}$ to yield mass-conservative and convergent solutions, although it may increase computational costs. The method can be effective especially for simulating unsaturated flow and other transport problems involving the propagation of a sharp-front.

• Journal of the Korean GEO-environmental Society
• /
• v.15 no.6
• /
• pp.5-15
• /
• 2014

This study is to quantitatively evaluate the impact on surrounding groundwater of waterway tunnel excavation and cofferdam construction in which A-dam and B-dam, so prediction of groundwater fluctuation and tunnel lining installation was studied. As a result, drawdown of groundwater level during tunnel excavation and cofferdam construction occurred about 3.58 m in the tunnel shaft. The initial condition of groundwater level recovered by up to 90 % was simulated after the completed the construction of the tunnel and lining installation. Groundwater inflow in the tunnel evaluated was analyzed to have exceeding water design criteria of the tunnel. The groundwater inflow is reduced to maximum $0.006m^3/min/km$ after lining installation done in the tunnel, so effect of lining installation was evaluated as 93 % or more. Drawdown of about 0.04~0.31 m occurs in the houses and temples analysis of groundwater system of the surrounding area from construction. Drawdown has occurred nearly by considering annual groundwater level fluctuation of National Groundwater Observation Network.

• Journal of the Korean GEO-environmental Society
• /
• v.16 no.5
• /
• pp.35-42
• /
• 2015

This is conducted to observe underground temperature and to analyze its change affected by climate condition and soil infiltration in the mountainous area, Yesan region, Chungcheong-namdo province. Additionally, underground temperature change is also simulated using air temperature and soil thermal properties with a numerical model. Soil temperature monitoring data acquired from each depth, 20 cm, 50 cm, and 100 cm, indicates that the data within 50 cm in depth shows peak-shaped big fluctuation directly affected by air temperature and it at 100 cm has open-shaped small fluctuation. Underground temperature variation, a difference between high and low values, during monitoring period is weakly proportional to hydraulic conductivity of the sediment and it is assumed that water plays a part in delivering air temperature in soil. The underground temperature estimated by a numerical model is very similar to the observed data with an average value of 0.99 cross-correlation coefficient. From the result of this study, the aquifer unsaturated hydraulic conductivity of the soil and the groundwater recharge is likely to be able to estimate with underground temperature profile calculated using a numerical model.

• Tunnel and Underground Space
• /
• v.13 no.1
• /
• pp.20-32
• /
• 2003

To understand the behavior of migration of contaminants in a fractured porous medium is a key to assure the overall safety of a potential radwaste repository. The feasible retention mechanism of contaminant transport in a tinctured medium are sorption of contaminants on solid surface and matrix diffusion of contaminants from a fracture into an adjacent porous medium. The acceleration mechanisms are the migration of contaminants in the form of pseudo-colloids and the limit of a volume f3r matrix diffusion. In this paper, the effects of these two acceleration mechanisms are studied mathematically, then semi-analytically computed by the application of the Talbot theorem and verified. Results indicate that the acceleration processes cannot be neglected in the modeling of contaminant transport in a fractured porous medium.