• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.11 no.4
• /
• pp.103-112
• /
• 1991

In this study, a numerical model was established and applied to simulate the steady-state groundwater and heat flow in an isotropic, heterogeneous, three dimensional aquifer system with uniform thermal properties and no change of state. This model was developed as an aid in screening large groundwater-flow systems as prospects for underground waste storage. Driving forces on the system are external hydrologic conditions of recharge from precipitation and fixed hydraulic head boundaries. Heat flux includes geothermal heat-flow, conduction to the land surface, advection from recharge, and advection to or from fixed-head boundaries. The model uses an iterative procedure that alternately solves the groundwater-flow and heat-flow equations, updating advective flux after solution of the groundwater-flow equation, and updating hydraulic conductivity after solution of the heat-flow equation. Dierect solution is used for each equation. Travel time is determined by particle tracking through the modeled space. Velocities within blocks are linear interpolations of velocities at block faces. Applying this model to the groundwater-flow system located in Jigyung-ri. Songla-myun, Youngil-gun. Kyungsangbuk-do, the groundwater-flow system including distribution of head, temperature and travel time and flow line, is analyzed.

• Journal of the Korean Society of Groundwater Environment
• /
• v.4 no.3
• /
• pp.111-115
• /
• 1997

Discrete fracture model has become one of the alternatives for the classical continuum model to simulate the irregular aspects of the fluid flow and the solute transport in fractured rocks. It is based on the assumptions that the discharge in a single fracture is proportional to the cube of the aperture and the fractured rock can be represented by the statistical assemblage of such single fractures. This study is intended to evaluate the effect of the fracture junction on the cubic law. Numerical solution of flow in junction system was obtained by using the Boundary-Fitted Coordinate System (BFCS) method. Results with different intersection angles in crossing fractures show that the geometry of the junction affects the discharge pattern under the same simulation conditions. Therefore, strict numerical and experimental examinations on this subject are required.

• Tunnel and Underground Space
• /
• v.14 no.4
• /
• pp.287-294
• /
• 2004

To predict the temperature distribution around a underground rock storage cavern, two- and three- dimensional numerical analysis using FLAC was conducted. The effects of groundwater and latent heat on thermal properties were considered in numerical calculation. The temperature estimated by FLAC are compared with the temperature measured for 5-year operation at Gonjiam storage cavern. Estimated and measured temperatures showed great discrepancy when thermal properties from laboratory tests were used and showed good agreement when the effects from 20% of volumetric water fraction and latent heat were considered. However, the discrepancy still increased with operation time due to the heat flow from ground surface. Three-dimensional numerical models were established to closely approximate the boundary condition of the test site, and numerical results better agreement when groundwater and latent heat effects were considered.

• Journal of the Korean Society of Groundwater Environment
• /
• v.6 no.2
• /
• pp.101-106
• /
• 1999

Open dump of refuse causes groundwater and soil contamination by leachate and air pollution by LFG(Landfill Gas). In this paper, in order to perform a study about reduction of high leachate and LFG collection & control, using a 3-D, 2-phase, transient FDM model, the analysis of simultaneous flow of leachate and LFG has been carried out. In present numerical analysis it is assumed that 58 percents of LFG will evaporate to the ambient air and the recharge rate of a landfill be 12 percent of the average precipitation per year. All other data were excerpted at the point of 1995 when three refuse layers had been buried. From numerical analysis we concluded that maximum head value is approximately 26 mH2O<-에이치투오 (2.52 atm) in the center of the system and that installing venting trench plays an important role in landfill stabilization. Evan with the assumption of three layers constructed and low recharge rate applied, it is found that cumulative leachate and LFG productions will be 15.1 million 세제곱미터, 5.58 billion 세제곱미터, respectively after 40 years.

• Proceedings of the Korean Geotechical Society Conference
• /
• 2005.03a
• /
• pp.1426-1433
• /
• 2005

In this paper, reliability analysis of surface settlement by ground water drawdown is performed using a reliability-groundwater flow numerical model. The result is compared with that of the deterministic model to evaluate the influence of the uncertainty from hydraulic conductivity in the soft ground as well as to determine the range of hydraulic conductivity of grouted ground. From the analyses, it was found that probability of failure to exceed the tolerable settlement was very high, if the hydraulic conductivity of grouted ground is decided from the deterministic flow model only. Reliability analysis which evaluates variance of hydraulic conductivity should be used together with the deterministic model for grouting design of tunnels to prevent ground water drawdown.

• Journal of the Korean Geotechnical Society
• /
• v.29 no.7
• /
• pp.75-89
• /
• 2013

This paper describes a GIS-based geohydrologic methodology, called YSGWF (YonSei GroundWater Flow) for predicting the rainfall infiltration-groundwater flow of slopes. This physical-based model was developed by the combination of modified Green-Ampt model that considers the unsaturated soil parameters and GIS-based raster model using Darcy's law that reflects the groundwater flow. In the model, raster data are used to simulate the three dimensional inclination of bedrock surface as actual topographic data, and the groundwater flow is governed by the slope. Also, soil profile is ideally subdivided into three zones, i.e., the wetting band zone, partially saturated zone, and fully saturated zone. In the wetting band and partially saturated zones the vertical infiltration of water (rainfall) from surface into ground is modeled. When the infiltrated water recharges into the fully saturated zone, the horizontal flow of groundwater is introduced. A comparison between the numerical calculation and real landslide data shows a reasonable agreement, which indicate that the model can be used to simulate real rainfall infiltration-groundwater flow.

• Tunnel and Underground Space
• /
• v.16 no.4 s.63
• /
• pp.281-287
• /
• 2006

Intake of groundwater by tunneling in a mountainous area mostly results from groundwater flow through fractured parts of total rock mass. For reasonable analysis of this phenomenon the representative joint groups 1, 2, and 3 have been selected by previous investigations, geological/geophysical field tests and boring works. Three dimensional fractures were generated by the FracMan and MAFIC which is a three dimensional finite element model has been used to analyse a groundwater flow through fractured media. Monte Carlo simulation was applied to reduce the uncertainty of this study. The numerical results showed that the average and deviation of amounts of groundwater intaked into tunnel per unit length were $5.40{\times}10^{-1}$ and $3.04{\times}10^{-1}m^3/min/km$. It is concluded that tunnel would be stable on impact of groundwater environment by tunneling because of the lower value than $2.00{\sim}3.00m^3/min/km$ as previous and present standard on the application of tunnel construction.

• Journal of Soil and Groundwater Environment
• /
• v.7 no.4
• /
• pp.10-16
• /
• 2002

The multi-aquifer system separated by semipervious leaky beds was analyzed. The finite difference scheme of the Crank-Nicolson method is applied to obtain the solution for this system. The solution of this scheme was compared with the analytical solution for two-layer aquifer systems with one-dimensional steady state. The results showed a good agreement between analytical and numerical solution for two-layer aquifer systems. So, the numerical scheme can be extended to multi-aquifer system. When the pumping is tried for single or multi aquifer, the computation of the groundwater heads was possible for each aquifer in the multi-aquifer with two-dimensional system. So, it might be helpful for the effective groundwater management.

• Geomechanics and Engineering
• /
• v.15 no.3
• /
• pp.851-859
• /
• 2018

There are two primary causes of the ground movement due to tunnelling in urban areas; firstly the lost ground and secondly the groundwater depression during construction. The groundwater depression was usually not considered as a cause of settlement in previous research works. The main purpose of this study is to analyze the combined effect of these two phenomena on the transverse settlement trough. Centrifuge model tests and numerical analysis were primarily selected as the methodology. The characteristics of settlement trough were analyzed by performing centrifuge model tests where acceleration reached up to 80g condition. Two different types of tunnel models of 180 mm diameter were prepared in order to match the prototype of a large tunnel of 14.4 m diameter. A volume loss model was made to simulate the excavation procedure at different volume loss and a drainage tunnel model was made to simulate the reduction in pore pressure distribution. Numerical analysis was performed using FLAC 2D program in order to analyze the effects of various groundwater depression values on the settlement trough. Unconfined fluid flow condition was selected to develop the phreatic surface and groundwater level on the surface. The settlement troughs obtained in the results were investigated according to the combined effect of excavation and groundwater depression. Subsequently, a new curve is suggested to consider elastic settlement in the modified Gaussian curve. The results show that the effects of groundwater depression are considerable as the settlement trough gets deeper and wider compared to the trough obtained only due to excavation. The relationships of maximum settlement and infection point with the reduced pore pressure at tunnel centerline are also suggested.

• The Journal of Engineering Geology
• /
• v.29 no.2
• /
• pp.113-122
• /
• 2019

In this study, the flow analysis to evaluate the extent of groundwater decline and the effect of the small valleys caused by the decrease of groundwater level in the construction of road tunnel, and the pollutant movement analysis to evaluate pollution of nearby water source by pollutant discharge during tunnel construction, respectively. The decrease of the groundwater during the 30 month tunnel excavation period was maximum 27 m and it was found to be the largest within 50 m from the tunnel center. The flow of groundwater is shown in the form of flowing into the tunnels and the effects of groundwater level decline were observed up to a tunnel radius of 200 m. As a result of the numerical modeling of the contaminant transport to examine the influence of the polluted water discharge from the tunnel, the range of the turbid water generated at the end of the tunnel is up to 120 m and it is estimated that the risk of contamination of the small river is not large.