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

• International Journal of Highway Engineering
• /
• v.9 no.3
• /
• pp.111-121
• /
• 2007

During the whole month of December in 2005, Korea experienced both heavy snowfall and freezing temperature in southeast regions, which had caused frost related damages to many pavements laid on top of box culverts. In-situ observation revealed that the formation of ice lenses in subgrade and subsequent unbound layers led to upward heaving and transverse cracks in concrete and asphalt pavements. This has affected the long-term performance of pavements, as well as has threatened drivers' safety for a while. Recently, Korea Expressway Corporation has proposed a design guide to better protect newly constructed unbound pavement layers over culverts from frost heave. A trench drainage system has been selected to effectively draw off water and to alleviate pore-water pressure in soils during the coldest season. This paper presents experimental and analytical backgrounds behind this new design guide. Soil specimens retrieved from the sites are tested to quantify clay content and to estimate the permeability of subgrade. A 2-D ground seepage analysis has been conducted to better understand the changes in pore water pressures as a function of grain size. Finally, an optimum size of trench drainage is determined based on numerical analysis and workability in the field.

• Geophysics and Geophysical Exploration
• /
• v.24 no.4
• /
• pp.202-208
• /
• 2021

Induced polarization (IP) was first published in a Korean academic journal in 1973, and it was soon applied to coal and metal ore exploration. Then, in universities and research institutes, IP modeling studies using the finite element approach and experimental studies on IP responses for artificial samples were conducted. In the mid-1980s, the spectral IP (SIP) measurement module was introduced to Korea, and physical scale modeling and inversion approaches were developed. Due to the decline of the mineral resource industry, this method was not actively applied. However, the SIP method was not applied In the 1990s, IP exploration was applied in the investigation of hydrothermal deposits of sulfide minerals and bentonite mineralization zones, as well as to areas where the groundwater was contaminated by intruding seawater. In the 2000s, three-dimensional inversion of the IP approach was developed, and high-precision geophysical exploration was required to secure domestic and overseas mineral resources, so SIP experiments on rock samples and approaches for field exploration were developed. The SIP approach was proven useful for the exploration of metal deposits containing sulfide minerals by applying it to explore the mineralization zone of gold-silver deposits in the Haenam region. The IP method is considered to be effective in exploring critical minerals (lithium, cobalt, and nickel) in high-tech industries. It also is expected to be useful for environmental and geotechnical investigations.

• Journal of the Korean Geotechnical Society
• /
• v.26 no.1
• /
• pp.35-43
• /
• 2010

The behavior of the Top-Base foundation was investigated by carrying out 3D finite element method. Special attention is given to the settlement behavior of concrete Top-Base foundation due to the consolidation settlement of the embedding depth and the effect of footing dimensions which are not included in the practical design. To obtain the detailed informations, a series of numerical analyses were performed for different pile configurations. It is shown that as the number of piles in a group increases, the calculated settlement also increases. However, for the $7\times7$ group, there is no further increase in settlement. Based on this study, it is found that the total settlement of Top-Base foundation is highly influenced by the consolidation settlement and footing configurations. It is also found that the current design method overestimates the settlement, and thus, needs to be modified and supplemented.

• Fire Science and Engineering
• /
• v.18 no.1
• /
• pp.24-30
• /
• 2004

In this work the dynamic heat transfer occurring in a cable penetration fire stop system built in the firewall of nuclear power plants is three-dimensionally investigated to develop a test-simulator that can be used to verify effectiveness of the sealants. Here was carried out an experiment to observe the heat transfer in the cable penetration fire stop system made of DOW CORNING products. The dynamic heat transfer occurring in the fire stop system is formulated in a parabolic partial differential equation subjected to a set of initial and boundary conditions. And it was modeled, simulated, and analyzed. The simulation results were illustrated in three-dimensional graphics and were compared with experimental data. Through the simulations, it was shown clearly that the temperature distribution was influenced very much by the number, position, and temperature of the cable streams. It also was found that the dynamic heat transfer through the cable streams was one of the most dominant factors, and the feature of heat conduction could be understood as an unsteady-state process. It is certain that these numerical results are useful for making a performance-based design for the cable penetration fire stop system.

• Geophysics and Geophysical Exploration
• /
• v.21 no.3
• /
• pp.150-161
• /
• 2018

Recently, unmanned aircraft EM (electromagnetic) survey based on ICT (Information and Communication Technology) has been widely utilized because of the efficiency in regional survey. We performed the theoretical study on the unmanned airship EM system developed by KIGAM (Korea Institute of Geoscience and Mineral resources) as part of the practical application of unmanned aircraft EM survey. Since this system has different configurations of transmitting and receiving loops compared to the conventional aircraft EM systems, a new technique is required for the appropriate interpretation of measured responses. Therefore, we proposed a method to calculate the EM field for the arbitrary shaped transmitter and verified its validity through the comparison with analytic solution for circular loop. In addition, to simulate the magnetic responses by three-dimensionally (3D) distributed anomalies, we have adapted our algorithm to 3D frequency-domain EM modeling algorithm based on the edge-FEM (finite element method). Though the analysis on magnetic field responses from a subsurface anomaly, it was found that the response decreases as the depth of the anomaly increases or the flight altitude increases. Also, it was confirmed that the response became smaller as the resistivity of the anomaly increases. However, a nonlinear trend of the out-of-phase component is shown depending on the depth of the anomaly and the used frequency, that makes it difficult to apply simple analysis based on the mapping of the magnitude of the responses and can cause the non-uniqueness problem in calculating the apparent resistivity. Thus, it is a prerequisite to analyze the appropriate frequency band and flight altitude considering the purpose of the survey and the site conditions when conducting a survey using the unmanned aircraft EM system.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.26 no.3
• /
• pp.174-183
• /
• 2014

Seabed beneath and near coastal structures may undergo large excess pore water pressure composed of oscillatory and residual components in the case of long durations of high wave loading. This excess pore water pressure may reduce effective stress and, consequently, the seabed may liquefy. If liquefaction occurs in the seabed, the structure may sink, overturn, and eventually increase the failure potential. In this study, to evaluate the liquefaction potential on the seabed, numerical analysis was conducted using the expanded 2-dimensional numerical wave tank to account for an irregular wave field. In the condition of an irregular wave field, the dynamic wave pressure and water flow velocity acting on the seabed and the surface boundary of the composite breakwater structure were estimated. Simulation results were used as input data in a finite element computer program for elastoplastic seabed response. Simulations evaluated the time and spatial variations in excess pore water pressure, effective stress, and liquefaction potential in the seabed. Additionally, the deformation of the seabed and the displacement of the structure as a function of time were quantitatively evaluated. From the results of the analysis, the liquefaction potential at the seabed in front and rear of the composite breakwater was identified. Since the liquefied seabed particles have no resistance to force, scour potential could increase on the seabed. In addition, the strength decrease of the seabed due to the liquefaction can increase the structural motion and significantly influence the stability of the composite breakwater. Due to limitations of allowable paper length, the studied results were divided into two portions; (I) focusing on the dynamic response of structure, acceleration, deformation of seabed, and (II) focusing on the time variation in excess pore water pressure, liquefaction, effective stress path in the seabed. This paper corresponds to (II).

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.26 no.3
• /
• pp.160-173
• /
• 2014

Seabed beneath and near coastal structures may undergo large excess pore water pressure composed of oscillatory and residual components in the case of long durations of high wave loading. This excess pore water pressure may reduce effective stress and, consequently, the seabed may liquefy. If liquefaction occurs in the seabed, the structure may sink, overturn, and eventually increase the failure potential. In this study, to evaluate the liquefaction potential on the seabed, numerical analysis was conducted using the expanded 2-dimensional numerical wave tank to account for an irregular wave field. In the condition of an irregular wave field, the dynamic wave pressure and water flow velocity acting on the seabed and the surface boundary of the composite breakwater structure were estimated. Simulation results were used as input data in a finite element computer program for elastoplastic seabed response. Simulations evaluated the time and spatial variations in excess pore water pressure, effective stress, and liquefaction potential in the seabed. Additionally, the deformation of the seabed and the displacement of the structure as a function of time were quantitatively evaluated. From the results of the analysis, the liquefaction potential at the seabed in front and rear of the composite breakwater was identified. Since the liquefied seabed particles have no resistance to force, scour potential could increase on the seabed. In addition, the strength decrease of the seabed due to the liquefaction can increase the structural motion and significantly influence the stability of the composite breakwater. Due to limitations of allowable paper length, the studied results were divided into two portions; (I) focusing on the dynamic response of structure, acceleration, deformation of seabed, and (II) focusing on the time variation in excess pore water pressure, liquefaction, effective stress path in the seabed. This paper corresponds to (I).