• Tunnel and Underground Space
• /
• v.33 no.1
• /
• pp.57-69
• /
• 2023

In this study, we simulated induced seismicity in the Groningen natural gas reservoir using 2D hydro-mechanical coupled discrete element modelling (DEM). The code used is PFC2D (Particle Flow Code 2D), a commercial software developed by Itasca, and in order to apply to this study we further developed 1)initialization of inhomogeneous reservoir pressure distribution, 2)a non-linear pressure-time history boundary condition, 3)local stress field monitoring logic. We generated a 2D reservoir model with a size of 40 × 50 km² and a complex fault system, and simulated years of pressure depletion with a time range between 1960 and 2020. We simulated fault system failure induced by pressure depletion and reproduced the spatiotemporal distribution of induced seismicity and assessed its failure mechanism. Also, we estimated the ground subsidence distribution and confirmed its similarity to the field measurements in the Groningen region. Through this study, we confirm the feasibility of the presented 2D hydro-mechanical coupled DEM in simulating the deformation of a complex fault system by hydro-mechanical coupled processes.

• 한국방재학회:학술대회논문집
• /
• 2011.02a
• /
• pp.75-75
• /
• 2011

본 논문에서는 진동중인 교량 거더에 작용하는 풍하중을 산정하고 그에 따른 플러터 발생풍속을 예측하기 위하여 분산형 전산환경을 활용한 수치해석 연구를 수행하였다. 분산형 전산환경은 웹 포탈을 기반으로 수치해석 환경을 제공하는 수치풍동 시스템으로서, 전산유체역학(CFD : Computational Fluid Dynamics)에 대한 전문지식이 부족한 사용자들도 격자생성, 수치해석자를 이용한 계산, 가시화 등의 전 과정을 편리하게 수행할 수 있는 차세대 토목분야 연구 환경이다. 본 시스템은 그리드스피어(GfidSphere)를 기반으로 구성되었으며, 기본적으로 사용자 관리, 세션 관리, 그룹 관리, 레이아웃 관리 등을 제공하여 사용자가 포탈을 통해서 다양한 서비스를 쉽게 사용할 수 있는 환경을 구축하도록 도와준다. 수치해석을 위한 유체 지배방정식은 2차원 비정상 비압축성 RANS(Reynolds-Averaged Navier-Stokes) 방정식이며, pseudo compressibility 방법을 적용하였다. 비정상 유동장을 해석하기 위하여 이중시간 전진법(dual time stepping)을 사용하였으며, 수렴가속화를 위해 Multi-grid 기법을 적용하였다. 또한 난류 유동장 해석을 위해서 $k-{\omega}$ SST 난류 모델을 사용하였으며, 난류 천이 과정에서의 유동을 모사하기 위하여 Total stress limitation 방법을 적용하였다. 교량 거더의 연직과 회전방향의 2자유도 움직임을 모사하기 위하여 동적격자 기법을 도입하였다. 교량 거더 주변의 비정상 유동해석 결과를 통해, 거더 표면에서 떨어져나가는 크고 작은 와류의 영향으로 양력 및 모멘트 계수 그래프가 중첩된 진폭과 주기를 갖고 주기적으로 나타나는 것을 확인할 수 있었다. 또한 계산된 비정상 공기력을 적용한 2자유도 플러터 방정식을 통하여 플러터 발생풍속을 산정하였다. 최종적으로 본 연구에서 계산된 결과의 타당성을 검증하기 위하여 수치적으로 구한 플러터 발생풍속과 기존의 실험 및 수치해석 결과를 비교하였으며, 결과는 잘 일치하였다.

• Journal of computational fluids engineering
• /
• v.1 no.1
• /
• pp.71-80
• /
• 1996

An iterative time marching procedure for solving incompressible turbulent flow has been applied to the flows around a high speed train including cross-wind effects. This procedure solves three-dimensional unsteady incompressible Reynolds-averaged Navier-Stokes equations on a non-orthogonal curvilinear coordinate system using first-order accurate schemes for the time derivatives and third/second-order accurate schemes for the spatial derivatives. Turbulent flows have been modeled by Baldwin-Lomax turbulent model. To validate present procedure, the flow around a high speed train at zero yaw angle was simulated and compared with experimental data. Generally good agreement with experiments was achieved. The flow fields around the high speed train at 9.2°, 16.7°, and 45° of yaw angle were also simulated.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.41 no.7
• /
• pp.481-487
• /
• 2017

The interaction between a shock wave and a boundary layer causes boundary layer separation, shock train, and in some cases, strong unsteadiness in the flow field. Such a situation is also observed in a shock tube, where the reflected shock wave interacts with the unsteady boundary layer. However, only a few studies have been conducted to investigate the shock train phenomenon in a shock tube. In the present study, numerical studies were conducted using the two-dimensional axisymmetric domain of a shock tube, and compressible Navier-Stokes equations were solved to clarify the flow characteristics of shock train phenomenon inside a shock tube. A detailed wave diagram was developed based on the present computational results, which were validated with existing experimental data.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.3 no.3
• /
• pp.31-39
• /
• 1999

In this paper, numerical experiment is attempted to analyze and compare the combustion efficiency of the burning sprays due to OFO, FOF triplet / FOOF split doublet injectors. Preconditioned Wavier-Stokes equation system with low Reynolds number $\kappa$-$\varepsilon$ model for turbulence closure, is LU-SGS time-integrated. Spray processes are modeled by DSF analysis with experimentally determined injection characteristics. n-heptane/air global reaction model approximates the combustion for simplicity, and the influence of turbulence on the chemical reaction is included using eddy dissipation model. The results showed the FOF triplet injector of highest combustion efficiency, whereas the OFO type of poet performance. It was also observed that the droplet mean diameter and the average gas temperature due to the mixing efficiency, are the representative parameters for the performance design of combustion.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.16 no.1
• /
• pp.33-42
• /
• 2003

This paper deals with the FEM analysis of nonlinear sloshing of incompressible, invicid and irrotational flow in two dimensional rectangular tank. We use laplace equation based on potential theory as governing equation. For large amplitude sloshing motion, kinematic and dynamic free surface conditions derived from Bernoulli equation are applied. This problem is solved by FEM using 9-node elements. For the time integration and accurate velocity calculation, we introduce predictor-corrector time marching scheme and least square method. Also, numerical stability in tracking of free surface is obtained by direct calculation of free surface location to time variation. Numerical results of sloshing induced by harmonic excitations, while comparing with those of linear theory and references, prove the accuracy and stability. After verification of our program, we analyze sloshing response characteristics to the fluid height and the excitation amplitude.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.19 no.2
• /
• pp.495-508
• /
• 1995

The self-excited vibrations of airfoil is related to the classical flutter problems, and it has been studied as a system with linear stiffness and small damping. However, since the actual aircraft wing and the many mechanical elements of airfoil type have various design variables and parameters, some of these could have strong nonlinearities, and the nonlinearities could be unexpectedly strong as the parameters vary. This abrupt chaotic behavior undergoes ordered routes, and the behaviors after these routes are uncontrollable and unexpectable since it is extremely sensitive to initial conditions. In order to study the chaotic behavior of the system, three parameters are considered, i.e., free-stream velocity, elastic distance and zero-lift angle. If the chaotic parameter region can be identified from the mathematically modeled nonlinear differential equation system, the designs which avoid chaotic regions could be suggested. In this study, by using recently developed dynamically system methods, and chaotic regions on the parameter plane will be found and the safe design variables will be suggested.

• Journal of the Society of Naval Architects of Korea
• /
• v.35 no.4
• /
• pp.1-10
• /
• 1998

This paper presents a vorticity-based numerical method for analyzing an incompressible Newtonian viscous flow around an impulsively started cylinder. The Navier-Stockes equations have a natural Helmholtz decomposition. The vorticity transport equation and the pressure equation are derived from this decoupled form. The associated boundary conditions are dynamic for the vorticity and pressure variables representing the coupling relation between them and the force balance on the wall. The various numerical treatments for solving the governing equations are introduced. According to Wu et al.(1994), the boundary conditions are decoupled, keeping the dynamic relation between vorticity and pressure. The vorticity transport equation is formulated by FVM and TVD(Total Variation Diminishing) scheme is used for the convection term. An integral approach similar to the panel method is used to obtain the velocity field for a given vorticity field and the pressure field, instead of the conventional differential approaches. In the numerical process, the structured grid is generated. The results are compared to existing numerical and analytic results for the validity of the present method.

• Journal of Energy Engineering
• /
• v.2 no.3
• /
• pp.251-257
• /
• 1993

The purpose of this study is to obtain an improved interpretation of heat transfer phenomena between blocks and fluids in the parallel conducting plates. Flow is two-dimensional, incompressible steady laminar flow over rectangular blocks, representing finite heat source on parallel plate. Heat transfer phenomena, temperature of blocks and heat transfer into the flow field are investigated for different spacings between blocks and Reynolds numbers. Results indicate that Nusselt number on the far upstream corner of the block was higher than that of any part of the block. As Reynolds number and spacings of blocks increased, Nusselt number increased. The distribution of local Nusselt number on the top surface of the conducting plate is similar to the case with insulated plate. Temperature of the block which has heat source in half cubage was approximately twice as high as temperature of the block which has heat source in whole cubage. As Reynolds number and spacings of blocks increased, overall temperature decreased. The peak value of block temperature occurred at position shifted to the right or upper right from center. The maximum temperature of block can be expressed as a function of Reynolds number, spacings between blocks, position of maximum temperature of each block and then it is possible to predict the maximum temperature of blocks.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.14 no.5
• /
• pp.1033-1042
• /
• 1994

Dynamic behavior of flexible rectangular liquid containers is analyzed by a two-dimensional coupled boundary element-finite element method. The irrotational motion of inviscid and incompressible ideal fluid is modeled by boundary elements and the motion of structure by finite elements. A singularity free integral formulation is employed for the implementation of boundary element method. Coupling is performed by using compatibility and equilibrium conditions along the interface between the fluid and structure. The fluid-structure interaction effects are reflected into the coupled equation of motion as added fluid mass matrix and sloshing stiffness matrix. By solving the eigen-problem for the coupled equation of motion, natural frequencies and mode shapes of coupled system are obtained. The free surface sloshing motion and hydrodynamic pressure developed in a flexible rectangular container due to horizontal and vertical ground motions are computed in time domain.