• Geomechanics and Engineering
• /
• v.11 no.1
• /
• pp.77-94
• /
• 2016

A forked tunnel, as a special complicated underground structure, is composed of big-arch tunnel, multi-arch tunnel, neighborhood tunnels and separate tunnels according to the different distances between two separate tunnels. Due to the complicated process of design and construction, surrounding jointed rock mass stability of the big-arch tunnel which belongs to the forked tunnel during excavation is a hot issue that needs special attentions. In this paper, an elasto-plastic damage constitutive model for jointed rock mass is proposed based on the coupling method considering elasto-plastic and damage theories, and the irreversible thermodynamics theory. Based on this elasto-plastic damage constitutive model, a three dimensional elasto-plastic damage finite element code (D-FEM) is implemented using Visual Fortran language, which can numerically simulate the whole excavation process of underground project and perform the structural stability of the surrounding rock mass. Comparing with a popular commercial computer code, three dimensional fast Lagrangian analysis of continua (FLAC3D), this D-FEM has advantages in terms of rapid computing process, element grouping function and providing more material models. After that, FLAC3D and D-FEM are simultaneously used to perform the structural stability analysis of the surrounding rock mass in the forked tunnel considering three different computing schemes. The final numerical results behave almost consistent using both FLAC3D and D-FEM. But from the point of numerically obtained damage softening areas, the numerical results obtained by D-FEM more closely approach the practical behaviors of in-situ surrounding rock mass.

• Nuclear Engineering and Technology
• /
• v.42 no.6
• /
• pp.620-635
• /
• 2010

In this paper we describe current activities on the project Multi-Scale Modeling and Analysis of convective boiling (MSMA), conducted jointly by the Paul Scherrer Institute (PSI) and the Swiss Nuclear Utilities (Swissnuclear). The long-term aim of the MSMA project is to formulate improved closure laws for Computational Fluid Dynamics (CFD) simulations for prediction of convective boiling and eventually of the Critical Heat Flux (CHF). As boiling is controlled by the competition of numerous phenomena at various length and time scales, a multi-scale approach is employed to tackle the problem at different scales. In the MSMA project, the scales on which we focus range from the CFD scale (macro-scale), bubble size scale (meso-scale), liquid micro-layer and triple interline scale (micro-scale), and molecular scale (nano-scale). The current focus of the project is on micro- and meso-scales modeling. The numerical framework comprises a highly efficient, parallel DNS solver, the PSI-BOIL code. The code has incorporated an Immersed Boundary Method (IBM) to tackle complex geometries. For simulation of meso-scales (bubbles), we use the Constrained Interpolation Profile method: Conservative Semi-Lagrangian $2^{nd}$ order (CIP-CSL2). The phase change is described either by applying conventional jump conditions at the interface, or by using the Phase Field (PF) approach. In this work, we present selected results for flows in complex geometry using the IBM, selected bubbly flow simulations using the CIP-CSL2 method and results for phase change using the PF approach. In the subsequent stage of the project, the importance of effects of nano-scale processes on the global boiling heat transfer will be evaluated. To validate the models, more experimental information will be needed in the future, so it is expected that the MSMA project will become the seed for a long-term, combined theoretical and experimental program.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.14 no.2
• /
• pp.245-255
• /
• 1994

A two-dimensional Lagrangian finite-difference computer program is developed for the wave propagation analysis of impact phenomena. The numerical scheme is the standard method originally proposed by Von Neuman and Richtmyer, using artificial viscosity to smooth shock fronts. The material model used in the study is the standard hydrodynamic-elastic-plastic relations with Von-Mises yield criterion. A test configuration consisted of a target and a projectile were calculated to understand the response of a colliding event. However, the computer code is in plane strain, the calculations were intended for generating the qualitative features of the model behaviors. Nevertheless, the computational results were consistent with the experimental observations and provided a rational basis to interpret the modes of failures.

• Korea-Australia Rheology Journal
• /
• v.18 no.2
• /
• pp.109-118
• /
• 2006

A finite-element method was employed to analyze axisymmetric unsteady motion of a deformable bubble near the wall. In the present study a deformable bubble in a Newtonian medium near the wall was considered. In solving the governing equations a structured mesh generator was used to describe the collapse of highly deformed bubbles with the Arbitrary Lagrangian Eulerian (ALE) method being employed in order to capture the transient bubble boundary effectively. In order to check the accuracy of the present FE analysis we compared the results of our FE solutions with the result of the collapse of spherical bubbles in a large body of fluid in which solutions can be obtained using a 1D FE analysis. It has been found that 1D and 2D bubble deformations are in good agreement for spherically symmetric problems confirming the validity of the numerical code. Non-spherically symmetric problems were also solved for the collapse of bubble located near a plane solid wall. We have shown that a microjet develops at the bubble boundary away from the wall as already observed experimentally. We have discussed the effect of Reynolds number and distance of the bubble center from the wall on the transient collapse pattern of bubble.

• Journal of computational fluids engineering
• /
• v.17 no.2
• /
• pp.35-41
• /
• 2012

This paper presents incompressible Navier-Stokes solution algorithm for 2D Free-surface flow problems on the Cartesian mesh, which was implemented to run on Graphics Processing Units(GPU). The INS solver utilizes the variable arrangement on the Cartesian mesh, Finite Volume discretization along Constrained Interpolation Profile-Conservative Semi-Lagrangian(CIP-CSL). Solution procedure of incompressible Navier-Stokes equations for free-surface flow takes considerable amount of computation time and memory space even in modern multi-core computing architecture based on Central Processing Units(CPUs). By the recent development of computer architecture technology, Graphics Processing Unit(GPU)'s scientific computing performance outperforms that of CPU's. This paper focus on the utilization of GPU's high performance computing capability, and presents an efficient solution algorithm for free surface flow simulation. The performance of the GPU implementations with double precision accuracy is compared to that of the CPU code using an representative free-surface flow problem, namely. dam-break problem.

• Journal of ILASS-Korea
• /
• v.17 no.2
• /
• pp.64-70
• /
• 2012

Water mist fire suppression systems which use relatively small droplets of water with high injection pressure are increasingly being used in wider applications because of its greater efficiency, low flooding damage and low toxicity. However, the performance of the system significantly relies on the water mist characteristics and it requires better understanding of fire suppression mechanism of water mist. In the present study, computational fluid dynamics simulations were carried out to investigate cooling performance of water mist in heated chamber. The gas phase was prepared with natural convection heat transfer model for incompressible ideal case and then the effects of water mist injection characteristics on cooling capabilities were investigated upon the basis of the pre-determined temperature field. For the simulation of water mist behavior, Lagrangian discrete phase model was employed by using a commercial code, FLUENT. Smaller droplet sizes, greater injection angles and higher flow rates provided relatively higher cooling performance.

• Journal of computational fluids engineering
• /
• v.12 no.2
• /
• pp.37-45
• /
• 2007

The hybrid Cartesian/immersed boundary method is applied to simulate fluid-structure interaction of a two-dimensional orbiting flexible foil. The elastic deformation of the flexible foil is modelled based on the dynamic equation of a thin-plate. At each time step, the locations and velocities of the Lagrangian control points on the flexible foil are used to reconstruct the boundary conditions for the flow solver based on the hybrid staggered/non-staggered grid. To test the developed code, the flow fields around a flapping elliptical wing are calculated. The time history of the vertical force component and the evolution of the vorticity fields are compared with recent other computations and good agreement is achieved. For the orbiting flexible foil, the vorticity fields are compared with those of the case without the deformation. The combined effects of the angle of attack and the orbit on the deformation are investigated. The grid independency study is carried out for the computed time history of the deformation at the tip.

• Coupled systems mechanics
• /
• v.2 no.1
• /
• pp.85-110
• /
• 2013

Generally, mass concrete structural behavior is governed by the strain components. However, relevant guidelines in dam engineering evaluate the structural behavior of concrete dams using stress-based criteria. In the present study, strain-based criteria are proposed for the first time in a professional manner and their applicability in seismic failure evaluation of an arch dam are investigated. Numerical model of the dam is provided using NSAD-DRI finite element code and the foundation is modeled to be massed using infinite elements at its far-end boundaries. The coupled dam-reservoir-foundation system is solved in Lagrangian-Eulerian domain using Newmark-${\beta}$ time integration method. Seismic performance of the dam is investigated using parameters such as the demand-capacity ratio, the cumulative inelastic duration and the extension of the overstressed/overstrained areas. Real crack profile of the dam based on the damage mechanics approach is compared with those obtained from stress-based and strain-based approaches. It is found that using stress-based criteria leads to conservative results for arch action while seismic safety evaluation using the proposed strain-based criteria leads to conservative cantilever action.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2014.10a
• /
• pp.411-414
• /
• 2014

Numerical modeling was conducted to estimate the amount of dislocation that may occur across a frictionless fracture during an earthquake using commercial code FLAC3D (Fast Lagrangian Analysis of Continua in 3 Dimensions). The applied motion was calculated to represent a Richter 6.0 magnitude earthquake at distances of 2 km from the fracture. The velocity-time history was generated from Svensk $K{\ddot{a}}arnbr{\ddot{a}}anslehantering$ AB report. In the report, The velocity field resulting from an earthquake on a fault located in the near-field (2 km distance) was modelled using a finite difference program, WAVE. The stress-time history was substituted for velocity-time history to perform dynamic analysis using FLAC3D. During the earthquake, the maximum dislocation and change of shear stress were about 1 cm and 2MPa, respectively. Because the fracture is frictionless in this study, all dislocations relax to zero after the earthquake motions have ceased.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2017.05a
• /
• pp.451-455
• /
• 2017

The reactive Eulerian-Lagrangian code is utilized to simulate combustion of ZPP/THPP in a closed vessel. In the paper, ignition delay of THPP is mainly studied since ZPP and THPP are isolated by a boron nitride wall. Only a numerical case study is conducted as experimental observation is inaccessible. Results showed THPP ignition delay affects initial shock strength thus not only the first peak become weak, but also the frequency of a pressure oscillation is slowed.