• Explosives and Blasting
• /
• v.22 no.4
• /
• pp.7-15
• /
• 2004

An explosion modeling technique was developed by using the spherical discrete element code, PFC3D, which can be used to model the dynamic stress wave propagation phenomenon. The modeling technique is simply based on an idea that the explosion pressure should be applied to a PFC3D particle assembly not in the form of an external force (body force), but in the form of a contact force (surface force). According to this concept, the explosion pressure is applied to the wall particles by the scheme of radius expansion/contraction of inner-hole particles. The output wall force is compared to the input hole pressure in every time step, and a correction routine is activated to control the radius multiplier of the inner-hole particles. A comparative blast simulation far a cement mortar block of $80\times90\times80mm$ was conducted by using the conventional explosion modeling method and the new one. The results of the simulation are presented in a qualitative fashion.

• Journal of the Korean Society of Groundwater Environment
• /
• v.6 no.4
• /
• pp.159-170
• /
• 1999

The hydrogeochemical study on the $CO_2$-rich water in the Chojeong area was carried out. The $CO_2$-rich water of Ca-$HCO_3$type is characterized by low pH (5.0~5.8). high $CO_2$concentration ($Pco_2$＜$10^{0.31}$atm) and high TDS. The water chemistry indicates that the $CO_2$-rich water was probably evolved by the local suppy of deep seated $CO_2$gas resulting in the enhanced water/rock (granite) interaction under low pH conditions. High $NO_3$concentration indicates that the $CO_2$water was mixed and diluted with low $CO_2$groundwater in the vicinity of the area, in which the extensive groundwater abstraction occurred during the past years. The evoiution of the $CO_2$-rich water in the Chojeong area for the process of $CO_2$injection water/rock interaction and mixing processes was thermodynamically simulated by PHREEQC. Although the simulation was limited to water/plagioclase interaction, the results show the feasible explanation about the observed trend of pH and Ca and Na concentrations of the $CO_2$-rich water.

• Economic and Environmental Geology
• /
• v.49 no.5
• /
• pp.371-380
• /
• 2016

This study aims to identify the mineraloical and petrographical characteristics of caprock from drilling cores of Pohang basin as a potential $CO_2$ storage site. Experiments and modeling were conducted in order to investigate the geochemical and mineralogical caprock effects of carbon dioxide. A series of autoclave experiments were conducted to simulate the interaction in the $scCO_2$-caprock-brine using a high pressure and temperature cell at $50^{\circ}C$ and 100 bar. Geochemical and mineralogical alterations after 15 days of $scCO_2$-caprock-brine sample reactions were quantitatively examined by XRD, XRF, ICP-OES investigation. Results of mineralogical studies, together with petrographic data of caprock and data on the physicochemical parameters of brine were used for geochemical modeling. Modelling was carried out using the The Geochemist's Workbench 11.0.4 geochemical simulator. Results from XRD analysis for caprock sample showed that major compositional minerals are quartz, plagioclase, and K-feldspar, and muscovite, pyrite, siderite, calcite, kaolinite and montnorillonite were included on a small scale. Results from ICP-OES analysis for brine showed that concentration of $Ca^{2+}$, $Na^+$, $K^+$ and $Mg^{2+}$ increased due to dissolution of plagioclase, K-feldspar and muscovite. Results of modeling for the period of 100 years showed that the recrystallization of kaolinite, dawsonite and beidellite, at the expense of plagioclase and K-feldspar is characteristic. Volumes of newly precipitation minerals and minerals passing into brine were balanced, so the porosity remained nearly unchanged. Experimental and modeling results indicate the interaction between caprock and $scCO_2$ during geologic carbon sequestration can exert significant impacts in brine pH and solubility/stability of minerals.

• Proceedings of the Korean Geotechical Society Conference
• /
• 1998.05a
• /
• pp.35-81
• /
• 1998

Distinct Element Method(DEM) has a great advantage to model the discontinuous behaviour of jointed rock masses such as rotation, sliding, and separation of rock blocks. Geometrical data of joints by a field monitoring is not enough to model the jointed rock mass though the results of DE analysis for the jointed rock mass is most sensitive to the distributional properties of joints. Also, it is important to use a properly joint law in evaluating the stability of a jointed rock mass because the joint is considered as the contact between blocks in DEM. In this study, a stochastic modelling technique is developed and the dilatant rock joint is numerically modelled in order to consider th geometrical and mechanical properties of joints in DE analysis. The stochastic modelling technique provides a assemblage of rock blocks by reproducing the joint distribution from insufficient joint data. Numerical Modelling of joint dilatancy in a edge-edge contact of DEM enable to consider not only mechanical properties but also various boundary conditions of joint. Preprocess Procedure for a stochastic DE model is composed of a statistical process of raw data of joints, a joint generation, and a block boundary generation. This stochastic DE model is used to analyze the effect of deviations of geometrical joint parameters on .the behaviour of jointed rock masses. This modelling method may be one tool for the consistency of DE analysis because it keeps the objectivity of the numerical model. In the joint constitutive law with a dilatancy, the normal and shear behaviour of a joint are fully coupled due to dilatation. It is easy to quantify the input Parameters used in the joint law from laboratory tests. The boundary effect on the behaviour of a joint is verified from shear tests under CNL and CNS using the numerical model of a single joint. The numerical model developed is applied to jointed rock masses to evaluate the effect of joint dilation on tunnel stability.

• Tunnel and Underground Space
• /
• v.26 no.4
• /
• pp.266-273
• /
• 2016

We have numerically modeled thermal evolution of Ulleung Island after an emplacement of magma chamber. The disk-shape magma chamber is assumed to locate at 2.9 km beneath the island and has a diameter and a thickness of 10 km and 300 (or 600) m, respectively. The geothermal gradients evaluated from the numerical modeling coincide well with the range of the geotherms (${\sim}95^{\circ}C/km$) observed from the well logging. Although there are limitations in the application of the numerical results directly to the interpretation of the observed geotherms, we believe that an existence of a hot magma chamber in molten or in solidified state is the most plausible explanation for the observed geotherms.

• Explosives and Blasting
• /
• v.23 no.1
• /
• pp.47-54
• /
• 2005

An explosion modeling technique was developed by using the spherical discrete element code, $PFC^{3D}$, which can be used to model the dynamic stress wave propagation phenomenon. The modeling technique is simply based on an idea that the explosion pressure should be applied to a $PFC^{3D}$ particle assembly not in the form of an external force (body force), but in the form of a contact force (surface force). A test blast was conducted for a RC column, whose dimension was $600\times300\times1800$ in millimeters. The initial velocities of the surface movements were measured to be in the range of $14\~18\;m/s$ with the initiation times of $1.5\~2.0m$. Then the blasting procedure was simulated by using the modeling technique. The particle assembly representing the concrete was made of cement mortar and coarse aggregates, whose mirco-properties were obtained from the calibration processes. As a result, the modeling technique developed in this study made it possible for the burden to move with the velocity of $17\~24\;m/s$, which are slightly higher values compared to those of the test blast.

• Journal of the Korean earth science society
• /
• v.28 no.3
• /
• pp.367-373
• /
• 2007

A geoacoustic modeling has been developed to predict sound transmission through the submarine layers of sediment and rock. It demands a geoacoustic model with the measured, extrapolated, and predicted values of geoacoustic parameters controlling acoustic propagation. In the coastal areas of Okgye and Bukpyeong, the East Sea, the marine succession consists of Quaternary/Tertiary deposits and acoustic basement. The basement of Okgye coastal area is indicative of siliciclastics of the Pyeongan Group in Paleozoic, and the average velocities of P-wave and S-wave are 4276 m/s and 2400 m/s, respectively. The basement of Bukpyeong coastal area is indicative of limestone of the Joseon Supergroup in early Paleozoic, and the average velocities of P-wave and S-wave are 5542 m/s and 2742 m/s, respectively.

• Tunnel and Underground Space
• /
• v.19 no.5
• /
• pp.440-449
• /
• 2009

LCM test is one of the most powerful and reliable methods of experiment for the cutter head design and the performance prediction of TBM. In many cases, however, the predicted design model can be directly applied to the field design, because this test may have an uppermost limit in preparation and/or transportation of the large size rock samples and the test for the jointed rock mass is not easy. When the proper and reasonable numerical modeling is considered to overcome this limit, the most adequate cutter head design for TBM could be presented without any complicate preconsideration in the field. In this study, the crack propagation patterns dependent on the contact point of disc cutter and the angle of rock joint are analyzed for the rock specimen with a single joint using the UDEC. The authors could derive the appropriate contact points of disc cutters and their space with respect to the joint angle in rock mass thru the numerical analysis.

• Journal of the Korean earth science society
• /
• v.41 no.2
• /
• pp.137-146
• /
• 2020

A great number of geological research data have been produced by individually conducted researchers and then personally stored in domestic universities and research institutes. However, it is difficult to share data with other researchers owing to low and limited accessibility. The purpose of this study is to systematically establish metadata for inaccessible data, to manage them collectively and to provide opportunities for utilizing the data to those who require efficient research methods. Approximately 1,000 geological specimens (900 rocks and fossils, 100 thin sections) were gathered, along with their metadata such as high-resolution photographs, classification, name, owner, address, and geographical coordinates of the sample site, to establish their features. Additionally, 3D modeling data for 100 rocks and fossils were generated. On the basis of this study, it is possible for researchers to access and share crucial geological data that have a high potential to be lost and have been neglected in restricted spaces; by avoiding the wasted time, energy, and costs caused by repetitive collection of data, researchers may perform effective research and achieve qualified and competitive research results. Moreover, vulnerable and important geological data in the field can be protected from damage caused by indiscriminate, repetitive collection of specimens that have previously been secured. Through the establishment of additional metadata concerning the diversity of rocks, fossils, and thin sections kept at other universities and research institutes, much more data can be recognized, leading to advanced research results. Furthermore, specialized comparison and analysis of basic mineralogy and petrology knowledge are anticipated, based on the use of the metadata.

• Tunnel and Underground Space
• /
• v.32 no.6
• /
• pp.568-585
• /
• 2022

In the present study, we proposed a numerical method for simulating thermally induced fracture slip using a grain-based distinct element model (GBDEM). As a part of DECOVALEX-2023, the thermo-mechanical loading test on a saw-cut rock fracture conducted at the Korea Institute of Civil Engineering and Building Technology was simulated. In the numerical model, the rock sample including a saw-cut fracture was represented as a group of random Voronoi polyhedra. Then, the coupled thermo-mechanical behavior of grains and their interfaces was calculated using 3DEC. The key concerns focused on the temperature evolution, thermally induced principal stress increment, and fracture normal and shear displacements under thermo-mechanical loading. The comparisons between laboratory experimental results and the numerical results revealed that the numerical model reasonably captured the heat transfer and heat loss characteristics of the rock specimen, the horizontal stress increment due to constrained displacement, and the progressive shear failure of the fracture. However, the onset of the fracture slip and the magnitudes of stress increment and fracture displacement showed discrepancies between the numerical and experimental results. We expect the numerical model to be enhanced by continuing collaboration and interaction with other research teams of DECOVALEX-2023 Task G and validated in further study.