• Proceedings of the Earthquake Engineering Society of Korea Conference
• /
• 2001.09a
• /
• pp.83-90
• /
• 2001

The site effects of local geological conditions on seismic ground motion are performed using 2D numerical method. For the analysis, a numerical method far ground response analysis using FE-BE coupling method is developed. The total system is divided into two parts so called far field and near field. The far field is modeled by boundary element formulation using the multi-layered dynamic fundamental solution that satisfied radiational condition of wave. And this is coupled with near field modeled by finite elements. In order to verify the seismic response analysis, the results are compared with those of commercial code. As a result, it is shown that the developed method can be an efficient numerical method to solve the seismic response analysis of the site effect in 2D problem.

• Proceedings of the Korea Concrete Institute Conference
• /
• 1996.04a
• /
• pp.247-252
• /
• 1996

In the concrete structures, cracks occur in various causes and the cracks seriously affect the functions of structures. The analysis techniques of progressive crack in the concrete have been improved with the advance of numerical techniques. The discrete crack model used in finite element program for the analysis of progressive failure is very effective, but it can not be easily implemented into numerical procedures because of difficult handing of nodal points in finite element meshes for crack growth. This paper introduces one of the techniques which skips the difficulty. In this paper, the modeling of progressive failure using finite element formulation is explained for the analysis of concrete fracture. The discontinuous element using the discontinuous shape function and the dual mapping technique in the numerical integration are implemented into finite element code for this purpose. It is shown that developed finite element program can predict the quasi-brittle behavior of concrete including ultimate load. The comparisons of the analysis results with other data are also shown.

• Proceedings of the Korea Committee for Ocean Resources and Engineering Conference
• /
• 2006.11a
• /
• pp.113-120
• /
• 2006

In this paper, generation mechanisms of ocean freak waves are briefly introduced in the context of wave-current interaction phenomena. As an accurate and efficient numerical tool, the spectral element method is presented with general features and specific treatment for the wave-current interaction problem. The present model of the fluid motion is based on the Navier-Stokes equations incorporating a velocity-pressure formulation. In order to deal with the free surface motion, an Arbitrary Lagrangian-Eulerian (ALE) description is adopted. As an intermediate stage of development, solution procedure and characteristic aspects of the present modeling and numerical method features are addressed in detail, and numerical results for wave-current interaction is left as further study.

• Journal of computational fluids engineering
• /
• v.21 no.3
• /
• pp.8-14
• /
• 2016

A numerical analysis is conducted to estimate the core expansion and the energy behaviors induced by a core disruptive accident in a sodium-cooled fast reactor. The numerical formulation based on underwater explosion theory is carried out to simulate the core explosion inside the reactor vessel. The transient pressure, temperature and expansion of the core are examined by solving the equation of state and nonlinear governing equation of momentum conservation in one-dimensional spherical coordinates. The energy balance inside the computation domain is examined during the core expansion process. Heat transfer between the core and the sodium coolant, and the bubble rise during the expansion process are briefly investigated.

• 한국전산유체공학회:학술대회논문집
• /
• 2007.04a
• /
• pp.177-181
• /
• 2007

Present study examines the numerical issues of cell structure simulation for various regimes of detonation phenomena ranging from weakly unstable to highly unstable detonations. Inviscid fluid dynamics equations with $variable-{\gamma}$ formulation and one-step Arrhenius reaction model are solved by a MUSCL-type TVD scheme and 4th order accurate Runge-Kutta time integration scheme. A series of numerical studies are carried out for the different regimes of the detonation phenomena to investigate the computational requirements for the simulation of the detonation wave cell structure by varying the reaction constants and grid resolutions. The computational results are investigated by comparing the solution of steady ZND structure to draw out the minimum grid resolutions and the size of the computational domain for the capturing cell structures of the different regimes of the detonation phenomena.

• Journal of Mechanical Science and Technology
• /
• v.18 no.3
• /
• pp.460-470
• /
• 2004

A computational analysis of engineering problems with moving domain or/and boundary according to either Lagrangian or Eulerian approach may encounter inherent numerical difficulties, the extreme mesh distortion in the former and the material boundary indistinctness in the latter. In order to overcome such defects in classical numerical approaches, the ALE(arbitrary Lagrangian Eulerian) method is widely being adopted in which the finite element mesh moves with arbitrary velocity. This paper is concerned with the ALE finite element formulation, aiming at the dynamic response analysis of baffled fuel-storage container in yawing motion, for which the coupled time integration scheme, the remeshing and smoothing algorithm and the mesh velocity determination are addressed. Numerical simulation illustrating theoretical works is also presented.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.27 no.9
• /
• pp.1327-1334
• /
• 2003

Flow characteristics of the fluid laden with many particles in the two-dimensional channel are investigated using the Navier-Stokes equations coupled with the equation of motion of particles by direct numerical simulation. A four-step fractional step method with Crank-Nicolson scheme and ALE technique is used for P2P1 mixed finite element method. The motion and distribution of particles in the fluid is virtually described as a result of direct numerical simulation and the increase of viscosity is compared with theoretical equations. The effect of channel height on the relative viscosity and the tubular pinch effect are discussed.

• 한국연소학회:학술대회논문집
• /
• 2006.04a
• /
• pp.205-212
• /
• 2006

Three-dimensional numerical study was carried out for the investigation of the detonation wave structures propagating in tubes. Fluid dynamics equations and conservation equation of reaction progress variable were analyzed by a MUSCL-type TVD scheme and four stage Runge-Kutta time integration. Chemical reaction was modeled by using a simplified one-step irreversible kinetics model. The variable gas properties between unburned and burned states were considered by using variable specific heat ratio formulation. The unsteady computational results in three-dimension show the detailed mechanisms of rectangular and diagonal mode of detonation wave instabilities resulting same cell length but different cell width in smoked-foil record. The results for the small reaction constant shows the spinning mode of three-dimensional detonation wave dynamics, which was rarely observed in the previous numerical simulation of the detonation waves.

• Transactions on Electrical and Electronic Materials
• /
• v.16 no.3
• /
• pp.107-111
• /
• 2015

Our work is based on the development of a numerical model to develop a methodology for predicting the aging and breakdown in insulation due to the dynamics of space charge packets. The model of bipolar charge transports is proposed to simulate space charge dynamic for high DC voltage in law-density polyethylene (LDPE), taking into account the trapping and detrapping of recombination phenomena, this model has been developed and experimentally validation. Theoretical formulation of the physical problem is based on the Poisson, the continuity and the transport equations as well as on the appropriate models for injection. Numerical results provide temporal and local distributions of the electric field, the space charge density for the different kinds of charges, conduction and displacement current densities, and the external current.

• 한국전산유체공학회:학술대회논문집
• /
• 1999.05a
• /
• pp.11-22
• /
• 1999

A numerical prediction method has been proposed to predict non-linear free surface oscillation in a three-dimensional container. The fluid motions are numerically predicted with Navier-Stokes equations discretized in a Lagrangian scheme with sufficient numerical accuracy. The profile of a free surface is precisely represented with three-dimensional body-fitted coordinates (BFC), which are regenerated in each computational step on the basis of the arbitrary Lagrangian-Eulerian (ALE) formulation. In order to confirm the reliability of the computational method, it was firstly applied to three-dimensional flows within complicated-shaped rigid boundaries, such as curved pipes and ducts. Than it was applied to benchmark computations related to free surface oscillations. Following these basic verifications, non-linear sloshings in a cylindrical tank and transitions from sloshing to swirling motions were numerically predicted. Throughout these computations, the applicability of the present computational method has been confirmed and some of the predicted free surface motions were visualized as sequential images and animations to understand their dynamic futures.