• Proceedings of the Korean Society for Rock Mechanics Conference
• /
• 2000.09a
• /
• pp.105-115
• /
• 2000

As large underground projects such as radioactive waste disposal, hot water and heat storage, and geothermal energy become influential, the study, which consider all aspects of thermics, hydraulics and mechanics would be needed. Thermo Hydro-Mechanical coupling analysis is one of the most complex numerical technique because it should be implemented with the combined three governing equations to analyze the behavior of rock mass. In this study, finite element code, which is based on Biot's consolidation theory, was developed to analyze the thermo-hydro-mechanical coupling in continuum rock mass. To verify the implemented program, one-dimensional consolidation model under the isothermal and non-isothermal conditions was analyzed and was compared with the analytic solution. The parametric study on two-dimensional consolidation was also performed and the effects of several factors such as poisson's ratio and hydraulic anisotropy on rock mass behavior were investigated. In the future, this program would be revised to be used for analysis of general discontinuous media with incorporating discrete joint model.

• Tunnel and Underground Space
• /
• v.10 no.3
• /
• pp.355-365
• /
• 2000

As large underground projects such as radioactive waste disposal, hot water and heat storage, and geothermal energy become influential, the study, which consider all aspects of thermics, hydraulics and mechanics would be needed. Thermo-Hydro-Mechanical coupling analysis is one of the most complex numerical technique because it should be implemented with the combined three governing equations to analyze the behavior of rock mass. In this study, finite element code, which is based on Biot's consolidation theory, was developed to analyze the thermo-hydro-mechanical coupling in continuum rock mass. To verify the implemented program, one-dimensional consolidation model under the isothermal and non-isothermal conditions was analyzed and was compared with the analytic solution. The parametric study on two-dimensional consolidation was also performed and the effects of several factors such as poisson's ratio and hydraulic anisotropy on rock mass behavior were investigated. In the future, this program would be revised to be used for analysis of general discontinuous media with incorporating discrete joint model.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.9 no.3
• /
• pp.97-106
• /
• 1989

This paper is on the modeling of two-dimensional groundwater flow, which is the first step of the development of Dynamic System Model for groundwater flow and pollutant transport in subsurface porous media. The particular features of the model are its versatility and flexibility to deal with as many real-world problems as possible. Points as well as distributed sources/sinks are included to represent recharges/pumping and rainfall infiltrations. All sources/sinks can be transient or steady state. Prescribed hydraulic head on the Dirichlet boundaries and fluxes on Neumann or Cauchy boundaries can be time-dependent or constant. Sources/sinks strength over each element and node, hydraulic head at each Dirichlet boundary node and flux at each boundary segment can vary independently of each other. Either completely confined or completely unconfined aquifers, or partially confined and partially unconfined aquifers can be dealt with effectively. Discretization of a compound region with very irregular curved boundaries is made easy by including both quadrilateral and triangular elements in the formulation. Large-field problems can be solved efficiently by including a pointwise iterative solution strategy as an optional alternative to the direct elimination solution methed for the matrix equation approximating the partial differential equation of groundwater flow. The model also includes transient flow through confining leaky aquifers lying above and/or below the aquifer of interest. The model is verified against three simple cases to which analytical solutions are available. The groundwater flow model shall be combined with the model of pollutant transport in subsurface porous media. Then the combined model, with the applications of the Eigenvalue technique and the Dynamic system theory, shall be improved to the Dynamic System Model which can simulate the real groundwater flow and the pollutant transport accurately and effectively for the analyses and predictions.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.18 no.4
• /
• pp.383-389
• /
• 2006

In three dimensional wave field, a direct numerical simulation model, which is able to handle free-fall and slope-fall of a waterdrop due to overtopping, is proposed to evaluate the overtopping rate on the coastal structures like an absorbing revetment. A comparison between the numerical model and existing experimental results for overtopping rate was made to validate the proposed numerical model's accuracy, and showed fairly good agreement between them. It is confirmed in numerical and hydraulic(existing) model test that the overtopping quantity on a absorbing revetment becomes larger with an increase in Ursell number. Also, the overtopping rate estimated by 3-D numerical model is compared with it obtained by 2-D numerical model.

• Journal of Korea Water Resources Association
• /
• v.44 no.3
• /
• pp.249-262
• /
• 2011

In this study, in order to reduce the numerical oscillation due to the unbalance between source and flux terms as the HLLC scheme is applied to the flow analysis on the irregular bed topography, a unstructured finite volume model based on the well-balanced HLLC scheme and the shallow water equations is developed and applied to problems of dam-break waves. The well-balanced HLLC scheme considers directly the gradient of bed topography as the flux terms is calculated. This scheme provides the good numerical balance between the source and flux terms in the case of the application to the steady-state transcritical flow. To verify the numerical model developed in this study, it is applied to three cases of hydraulic model experiments and a field case study of Mapasset dam failure (France). As a result of the verification, the predicted numerical results agree relatively well with available laboratory and field measurements. The model provides slightly more accurate results compared with the existing models.

• Journal of Korea Water Resources Association
• /
• v.31 no.3
• /
• pp.347-360
• /
• 1998

Numerical simulations for fluid flow and solute transport in a fracture rock masses are performed by using a transient flow model, which is based on the three-dimensional stochastic and discrete fracture network model (DFN model) and is coupled hydraulic model with mechanical model. In the numerical simulations of the solute transport, we used to the particle following algorithm which is similar to an advective biased random walk. The purpose of this study is to predict the response of the tracer test between two deep bore holes (GPK1 and GPK2) implanted at Soultz sous Foret in France, in the context of the geothermal researches.l The data sets used are obtained from in situcirculating experiments during 1995. As the result of the transport simulation, the mean transit time for the non reactive particles is about 5 days between two bore holes.

• Journal of Environmental Science International
• /
• v.14 no.1
• /
• pp.41-51
• /
• 2005

Using three-dimensional non-hydrostatical numerical model with one way double nesting technique, atmo­spheric circulation in the mountainous coastal region in summer was investigated from August 13 through 15, 1995. During the day, synoptic westerly wind blows over Mt. Mishrung in the west of a coastal city, Sokcho toward the East Sea, while simultaneously, easterly upslope wind combined with both valley wind from plain (coast) toward mountain and sea-breeze from sea toward inland coast blows toward the top of the mountain. Two different directional wind systems confront each other in the mid of eastern slope of the mountain and the upslope wind goes up to the height over 2 km, becoming an easterly return flow in the upper level over the sea and making sea-breeze front with two kinds of sea-breeze circulations of a small one in the coast and a large one in the open sea. Convective boundary layer is developed with a thickness of about 1km over the ground in the upwind side of the mountain in the west and a thickness of thermal internal boundary layer from the coast along the eastern slope of the mountain is only confined to less than 200 m. On the other hand, after sunset, no prohibition of upslope wind generated during the day and downward wind combined with mountain wind from mountain towardplain and land-breeze from land toward under nocturnal radiative cooling of the ground surfaces should intensify westerly downslope wind, resulting in the formation of wind storm. As the wind storm moving down along the eastern slop causes the development of internal gravity waves with hydraulic jump motion in the coast, bounding up toward the upper level of the coastal sea, atmospheric circulation with both onshore and offshore winds like sea-breeze circulation forms in the coastal sea within 70 km until midnight and after that, westerly wind prevails in the coast and open seas.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2016.05a
• /
• pp.65-65
• /
• 2016

Various types of flow conditions are developed in the region just downstream of hydraulic structures such as weir and drop structures. One of distinct flow conditions occurred downstream of drop structures is the wave type flow with undular hydraulic jump formation. We present three-dimensional numerical simulations of a wave type flow formed downstream of a stepped weir which were experimentally investigated by Kang et al. (2010). The turbulent flow over the weir structure is modeling using the unsteady Reynolds-averaged Navier-Stokes (URANS) simulation employing the Spalart-Allmaras one equation model and the detached eddy simulation. Numerical modeling and the performance of turbulence modeling approaches are evaluated by comparing with the experimental measurements in terms of the free surface variation, the shapes and sizes of undular wave, roller near at free surface, recirculation zone near the channel bottom downstream of the structures, and streamwise velocity profiles at selected longitudinal locations.

• Proceedings of the Korean Environmental Sciences Society Conference
• /
• 2003.11a
• /
• pp.19-26
• /
• 2003

The dispersion of recycled particulates in the complex coastal terrain containing Kangnung city, Korea was investigated using a three-dimensional non-hydrostatic numerical model and lagrangian particle model (or random walk model). The results show that particulates at the surface of the city that float to the top of thermal internal boundary layer (TIBL) are then transported along the eastern slope of the mountains with the passage of sea breeze and nearly reach the top of the mountains. Those particulates then disperse eastward at this upper level over the coastal sea and finally spread out over the open sea. Total suspended particulate (TSP) concentration near the surface of Kangnung city is very low. At night, synoptic scale westerly winds intensify due to the combined effect of the synoptic scale wind and land breeze descending the eastern slope of the mountains toward the coast and further seaward. This increase in speed causes development of internal gravity waves and a hydraulic jump up to a height of about 1km above the surface over the city. Particulate matter near the top of the mountains also descends the eastern slope of the mountains during the day, reaching the central city area and merges near the surface inside the nocturnal surface inversion layer (NSIL) with a maximum ground level concentration of TSP occurring at 0300 LST. Some particulates were dispersed following the propagation area of internal gravity waves and others in the NSIL are transported eastward to the coastal sea surface, aided by the land breeze. The following morning, particulates dispersed over the coastal sea from the previous night, tend to return to the coastal city of Kangnung with the sea breeze, developing a recycling process and combine with emitted surface particulates during the morning. These processes result in much higher TSP concentration. In the late morning, those particulates float to the top of the TIBL by the intrusion of the sea breeze and the ground level TSP concentration in the city subsequently decreases.

• Geomechanics and Engineering
• /
• v.10 no.5
• /
• pp.599-615
• /
• 2016

The present research presents experimental and finite element studies to investigate the behavior of piled raft-tunnel system in a sandy soil. In the experimental work, a small scale model was tested in a sand box with load applied vertically to the raft through a hydraulic jack. Five configurations of piles were tested in the laboratory. The effects of pile length (L), number of piles in the group and the clearance distance between pile tip and top of tunnel surface (H) on the load carrying capacity of the piled raft-tunnel system are investigated. The load sharing percent between piles and rafts are included in the load-settlement presentation. The experimental work on piled raft-tunnel system yielded that all piles in the group carry the same fraction of load. The load carrying capacity of the piled raft-tunnel model was increased with increasing (L) for variable (H) distances and decreased with increasing (H) for constant pile lengths. The total load carrying capacity of the piled raft-tunnel model decreases with decreasing number of piles in the group. The total load carrying capacity of the piles relative to the total applied load (piles share) increases with increasing (L) and the number of piles in the group. The increase in (L/H) ratio for variable (H) distance and number of piles leads to an increase in piles share. ANSYS finite element program is used to model and analyze the piled raft-tunnel system. A three dimensional analysis with elastoplastic soil model is carried out. The obtained results revealed that the finite element method and the experimental modeling are rationally agreed.