• Tunnel and Underground Space
• /
• v.31 no.4
• /
• pp.270-288
• /
• 2021

We proposed a numerical method to simulate the hydro-mechanical behavior of rock fracture using a grain-based distinct element model (GBDEM) in the paper. As a part of DECOVALEX-2023 Task G, we verified the method via benchmarks with analytical solutions. DECOVALEX-2023 Task G aims to develop a numerical method to estimate the coupled thermo-hydro-mechanical processes within the crystalline rock fracture network. We represented the rock sample as a group of tetrahedral grains and calculated the interaction of the grains and their interfaces using 3DEC. The micro-parameters of the grains and interfaces were determined by a new methodology based on an equivalent continuum approach. In benchmark modeling, a single fracture embedded in the rock was examined for the effects of fracture inclination and roughness, the boundary stress condition and the applied pressure. The simulation results showed that the developed numerical model reasonably reproduced the fracture slip induced by boundary stress condition, the fracture opening induced by fluid injection, the stress distribution variation with fracture inclination, and the fracture roughness effect. In addition, the fracture displacements associated with the opening and slip showed good agreement with the analytical solutions. 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 experiments.

• Journal of the Korean Magnetics Society
• /
• v.23 no.1
• /
• pp.31-35
• /
• 2013

The objective of this study is numerically analyzed the behavior characteristics of the magnetic fluid in a closed rectangular container using finite element method (FEM). The governing equations are solved with magnetization and Maxwell equations for consideration of rotating effect of the magnetite particle. Then the discretized equations are solved with boundary conditions of the velocity and temperature. The developed model is validated with the results of Davis (1983) and Fusegi et al. (1991) has a good agreement within 5.5 % and 2.7 %, respectively.

• Tunnel and Underground Space
• /
• v.30 no.6
• /
• pp.573-590
• /
• 2020

This study presents the current status of DECOVALEX-2023 project Task G and our research results so far. Task G, named 'Safety ImplicAtions of Fluid Flow, Shear, Thermal and Reaction Processes within Crystalline Rock Fracture NETworks (SAFENET)' aims at developing a numerical method to simulate the fracture creation and propagation, and the coupled thermohydro-mechanical processes in fracture in crystalline rocks. The first research step of Task G is a benchmark simulation, which is designed for research teams to make their modelling codes more robust and verify whether the models can represent an analytical solution for displacements of a single rock fracture. We reproduced the mechanical behavior of rock and embedded single fracture using a three-dimensional grain-based distinct element model for the simulations. In this method, the structure of the rock was represented by an assembly of rigid tetrahedral grains moving independently of each other, and the mechanical interactions at the grains and their contacts were calculated using 3DEC. The simulation results revealed that the stresses induced along the embedded fracture in the model were relatively low compared to those calculated by stress analysis due to stress redistribution and constrained fracture displacements. The fracture normal and shear displacements of the numerical model showed good agreement with the analytical solutions. The numerical model will be enhanced by continuing collaboration and interaction with other research teams of DECOVALEX-2023 Task G and validated using various experiments in a further study.

• Journal of the Korean GEO-environmental Society
• /
• v.13 no.10
• /
• pp.55-62
• /
• 2012

In this research, the algorithm for simulating specific grain size distribution(GSD) with large diameter granular material was developed using the distinct element analysis program $PFC^{3D}$(Particle Flow Code). This modeling approach can generate the initial distinct elements without clump logic or cluster logic and prevent distinct element from escaping through the confining walls during the process. Finally the proposed distinct element model is used to simulate large triaxial compression test of the rockfill material and we compared the simulation output with lab test results. Simulation results of Assembly showed very well agreement with the GSD of the test sample and numerical modeling of granular material would be possible for various stress conditions using this application through the calibration.

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

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.36 no.11
• /
• pp.1119-1126
• /
• 2012

The collective behavior of bacteria plays an important role in biofilm development. In this study, the fluidic properties of biofilms formed in Escherichia coli (E. coli) and Bacillus subtilis (B. subtilis) colonies were compared by visualizing 200-nm fluorescent beads that were initially embedded in an agar plate and distributed spontaneously on the upper surface of the growing colonies. We conducted experiments to measure the three-dimensional profile of the E. coli colony using fluorescent microbeads that did not flow in the colony. Vortical flow patterns near the edge of the B. subtilis colony were observed clearly by tracking the movement of the beads in the biofilm of the colony. The present study should be the first step toward determining the effect of fluidic biofilms on the growth and swarming dynamics of bacteria.

• Journal of the Korean Society of Safety
• /
• v.13 no.4
• /
• pp.3-8
• /
• 1998

The vortex particle method, which includes viscous effects, consists of diffusion of boundary vorticity and creation of the vortex particles, convection, particle strength exchange, and particle redistribution. Accuracy of the boundary element method is very important since it creates the particles around the body at every time step. A boundary element method based on source panel was investigated as part of computation of unsteady separated flows by rising the vortex particle method. The potential flows were computed around a circular cylinder and a square cylinder. The results around the circular cylinder were compared with the exact solution, and the distribution of vorticity, in particular near the sharp comers of the square cylinder, is scrutinized for different number of panels.

• Journal of the Korean Society for Marine Environment & Energy
• /
• v.9 no.2
• /
• pp.85-91
• /
• 2006

This study deals with the flow over a flat plate with repeated roughness elements of 2-dimensional rectangular shape, which can be applied into the study on the natural geographical roughness and the turbulent flow on roughened solid surface. Part 1 of the study showed that the ratio between the spacing and the height of roughness elements plays a crucial role in developing the flow pattern near wall surface. The present study complements the turbulence characteristics, the utilization of friction velocity as well as integral parameters. Results confirmed that k-type roughness(s/H=7 or 14) is certainly a more effective means than d-type roughness (s/H=3.5) in thickening the viscous region.

• Journal of the Korean GEO-environmental Society
• /
• v.25 no.4
• /
• pp.21-28
• /
• 2024

In this study, a polyhedral domain decomposition method for Smoothed Particle Hydrodynamics (SPH) analysis is introduced. SPH which is one of meshless methods is a numerical analysis method for fluid flow simulation. It can be useful for analyzing fluidic soil or fluid-structure interaction problems. SPH is a particle-based method, where increased particle count generally improves accuracy but diminishes numerical efficiency. To enhance numerical efficiency, parallel processing algorithms are commonly employed with the Cartesian coordinate-based domain decomposition method. However, for parallel analysis of complex geometric shapes or fluidic problems under dynamic boundary conditions, the Cartesian coordinate-based domain decomposition method may not be suitable. The introduced polyhedral domain decomposition technique offers advantages in enhancing parallel efficiency in such problems. It allows partitioning into various forms of 3D polyhedral elements to better fit the problem. Physical properties of SPH particles are calculated using information from neighboring particles within the smoothing length. Methods for sharing particle information physically separable at partitioning and sharing information at cross-points where parallel efficiency might diminish are presented. Through numerical analysis examples, the proposed method's parallel efficiency approached 95% for up to 12 cores. However, as the number of cores is increased, parallel efficiency is decreased due to increased information sharing among cores.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2019.05a
• /
• pp.92-92
• /
• 2019

최근 기존의 격자방식의 해석 방법을 벗어나 해석 영역에 대한 분할을 별도로 고려치 않는 수치기법의 실무적 적용사례가 증가 하고 있으며, 이러한 방식중 SPH(Smoothed Particle Hydrodynamics) 방식이 근자에 수자원 분야에서도 도입되어 관수로 및 개수로 해석 또는 복합해석 등에 활용된 바 있다. 최초 도입된 무격자방식의 모형들은 주로 복잡한 형상을 지니는 유체기계 등에 활용성이 높았던 바, 큰 규모의 해석 도메인을 가지는 수자원 분야에서의 SPH의 실무적용 평가와 효율성의 확보를 위해서, 본 연구는 국내 댐을 시범 대상으로 하여 SPH 수치해석 툴을 적용하고자 하였다. 분석 대상댐은 국내 P댐으로서 관리수위의 변동은 크지 않으나, 댐 직상류의 만곡이 심하고 다수의 대규모 취수구를 가진 곳으로 상시 발전방류 및 수시 댐 수문방류에 의해 유체의 흐름이 2,3차원의 복잡성을 띄고있기 때문에, 3차원 전산유체역학 Tool의 적용이 적절한 것으로 판단하였다. 해석을 위해 하류경계조건을 댐축과 문비로 설정하였고 상류 1km까지를 해석의 도메인으로 설정하였다. 소요시간을 줄이기 위해 여러 번의 모의를 거쳐 입자의 평균 입경은 0.6m로 제안하고 시격은 1초 미만(평균 0.5초)로 설정하였다. 수문 및 발전방류는 해당댐의 1~2년 빈도 수준에 해당하는 $5,000m^3/s$ 이하의 유량을 기준으로 하여 모의를 수행하였다. 모의의 안정성을 확보한 이후에는 해당 댐지역의 하류영향을 고려한 문비제어를 반영한 다양한 방식의 수문운영 및 취수지점의 순간 수위 영향을 검토하였다. 그 결과로 본 모의에서는 특정한 수문의 운영 조건에서는 댐수위 계측지점과 인근 취수지점 간에도 0.2m 수준의 순간 수위차가 발생할 수 있음을 보였으며, 이는 경우에 따라 취수시설의 일시적 장애요소로 발생할 수 있음을 의미한다. 따라서, 현재의 취수구조물과 문비운영 특성에 따라 발생가능한 취수장애를 줄일 수 있는 운영조건의 탐색을 위해서 수치모의를 추가로 하였으며, 이 때 댐축 상류의 유속분포에 대한 추가 검토도 수행하였다. 다만, 댐에서 방류시 하류조건에 대한 검토는 추후 보강되어야 할 것으로 판단된다.