• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.23 no.8
• /
• pp.1063-1071
• /
• 1999

Fine grid calculations are reported for the developing turbulent flow in a curved duct of square cross-section with a radius of curvature to hydraulic diameter ratio ${\delta}=Rc/D_H=3.357 $ and a bend angle of 720 deg. A sequence of modeling refinements is introduced; the replacement of wall function by a fine mesh across the sublayer and a low Reynolds number algebraic second moment closure up to the near wall sublayer in which the non-linear return to isotropy model and the cubic-quasi-isotropy model for the pressure strain are adopted; and the introduction of a multiple source model for the exact dissipation rate equation. Each refinement is shown to lead to an appreciable improvement in the agreement between measurement and computation.

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

This paper presents groundwater flow characteristics in discontinuous rock mass using fracture network program(NAPSAC) by statistical approach. Equivalent permeability coefficients are estimated from borehole data around Mabuk test tunnel site and fracture map on the arch of the tunnel. The reliability of fracture network model is obtained from determination of input data for statistical fracture network analysis from the real data(data of fracture network, data of hydraulic tests). The variation of permeability and mean anisotropic permeability coefficients are calculated from the realized model by increasing the size. As a result of analysis, a strong anisotropy of permeability is observed according to the direction of the fracture sets around the test tunnel.

• Tunnel and Underground Space
• /
• v.10 no.3
• /
• pp.378-386
• /
• 2000

This paper presents groundwater flow characteristics in discontinuous rock mass using fracture network program(NAPSAC) by statistical approach. Equivalent permeability coefficients are estimated from borehole data around Mabuk test tunnel site and fracture map on the arch of the tunnel. The reliability of fracture network model is obtained from determination of input data for statistical fracture network analysis from the real data(data of fracture network, data of hydraulic tests). The variation of permeability and mean anisotropic permeability coefficients are calculated from the realized model by increasing the size. As a result of analysis, a strong anisotropy of permeability is observed according to the direction of the fracture sets around the test tunnel.

• Economic and Environmental Geology
• /
• v.50 no.1
• /
• pp.61-71
• /
• 2017

The effects of joint density and size distribution on the hydrogeologic characteristics of jointed rock masses are addressed through numerical experiments based on the 2-D DFN (discrete fracture network) fluid flow analysis. Using two joint sets, a total of 51 2-D joint network system were generated with various joint density and size distribution. Twelve fluid flow directions were chosen every $30^{\circ}$ starting at $0^{\circ}$, and total of 612 $20m{\times}20m$ DFN blocks were prepared to calculate the directional block conductivity. Also, the theoretical block conductivity, principal conductivity tensor and average block conductivity for each generated joint network system were determined. The directional block conductivity and chance for the equivalent continuum behavior of the 2-D DFN system were found to increase with the increase of joint density or size distribution. However, the anisotropy of block hydraulic conductivity increases with the increase of density discrepancy between the joint sets, and the chance for the equivalent continuum behavior were found to decrease. The smaller the intersection angle of the two joint sets, the more the equivalent continuum behavior were affected by the change of joint density and size distribution. Even though the intersection angle is small enough that it is difficult to have equivalent continuum behavior, the chance for anisotropic equivalent continuum behavior increases as joint density or size distribution increases.

• The Journal of Engineering Geology
• /
• v.10 no.1
• /
• pp.13-25
• /
• 2000

A three-dimensional discrete fracture network model based on probabilistic characteristics of fracture geometry and transmissivity was designed to calculate the conductivity tensor and to estimate theanisotropy of conductivity. The conductivities, $K_p$, obtained from the numerical simulation of single-holepacker test corresponded well to those from the field tests. From this, it can be concluded that thefracture network model designed in this study can represent hydraulic characteristics of in-situ fractured rock mass. Block-scale conductivities, $K_b$, estimated from the modelling of steady-state flow through the REV-scale block were ranged between the arithmetic mean and harmonic mean of theconductivity estimates from packer tests. The conductivity along north-south direction was 1.4 timesgreater than that along the east-west direction. It was concluded that the anisotropy of conductivitywas insignificant. It was also found that there was a little correlation between $K_b$ and $K_p$. This would be to that the conductivities from the packer test simulation was strongly dependent on thetransmissivity and the number of fractures within the packer test intervals.

• 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.

• Tunnel and Underground Space
• /
• v.30 no.3
• /
• pp.242-255
• /
• 2020

The numerical simulations of Heater Experiment-D (HE-D) at the Mont Terri rock laboratory in Switzerland were performed to investigate an applicability of FLAC3D to reproduce the coupled thermo-hydro-mechanical (THM) behaviour in Opalinus Clay, as part of the DECOVLEX-2015 project Task B. To investigate the reliability of numerical simulations of the coupled behaviour using FLAC3D code, the simulation results were compared with the observations from the in-situ experiment, such as temperature at 16 sensors, pore pressure at 6 sensors, and strain at 22 measurement points. An anisotropic heat conduction model, fluid flow model, and transversely isotropic elastic model in FLAC3D successfully represented the coupled thermo-hydraulic behaviour in terms of evolution for temperature and pore pressure, however, performance of the models for mechanical behavior is not satisfactory compared with the measured strain.

• Tunnel and Underground Space
• /
• v.26 no.6
• /
• pp.516-523
• /
• 2016

Roughness of a joint surface is one of the most important parameters that affects the mechanical and hydraulic behavior of rock mass. Therefore, various studies on making constitutive model and/or roughness quantification have been conducted in experimental and empirical manners. Advances in recent 3D printing technology can be utilized to generate a joint surface with a specific roughness. In this study, a reliable technique to generate a rough joint surface was introduced and its quantitative assessment was made. Random midpoint displacement method was applied to generate a joint surface and the distribution of $Z_2$ was investigated to assess its roughness. As a result, a certain roughness can be embodied by controlling input parameters and furthermore it was able to generate a joint surface with specific roughness anisotropy.

• The Journal of the Petrological Society of Korea
• /
• v.23 no.2
• /
• pp.75-103
• /
• 2014

The Quaternary faults are extensively observed along major inherited fault zones (i.e. Yangsan Fault System, Ulsan Fault, Yeonil Tectonic Line, Ocheon Fault System) in SE Korea. Their geometry and kinematics provide a very useful piece of information about the Quaternary crustal deformation and stress field in and around Korean Peninsula. Using magnetic fabrics (AMS), we attempted to determine the slip senses of Jinti, Mohwa, Suseongji2, and Wangsan faults and then interpreted the fabric development process of fault gouge and the characteristics of stress field during the Quaternary. All the magnetic fabrics of the faults, except the Wangsan Fault, consistently indicate a dominant reverse-slip sense with weak strike-slip component. Most of the oblate fabrics are nearly parallel to the fault surface and the anisotropy degrees generally increase in proportion to the oblatenesses. These results suggest that the fabrics of the fault gouges resulted from a progressive deformation due to continuous simple shear during the last reactivation stage as reverse faulting. It is also interpreted that the pre-existing fabrics were overwhelmed and obliterated by the re-activated faulting. Paleostress field calculated from the fault slip data indicates an ENE-WNW compressive stress, which is in accord with those determined from previous fault tectonic analysis, focal mechanism solution, and hydraulic fracturing test in and around Korean Peninsula.