• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.23 no.6
• /
• pp.675-682
• /
• 2010

The flood impact pressure acting on a vertical wall resulting from a dam-breaking problem is simulated using a navier-Stokes(N-S) solver. The N-S solver uses Eulerian Finite Volume Method(FVM) along with Volume Of Fluid(VOF) method for 2-D incompressible free surface flows. A Split Lagrangian Advection(SLA) scheme for VOF method is implemented in this paper. The SLA scheme is developed based on an algorithm of Piecewise Linear Interface Calculation(PLIC). The coupling between the continuity and momentum equations is affected by using a well-known Semi-Implicit Method for Pressure-Linked Equations (SIMPLE) algorithm. Several two-dimensional numerical simulations of the dam-breaking problem are presented to validate the accuracy and demonstrate the capability of the present algorithm. The significance of the time step and grid resolution are also discussed. The computational results are compared with experimental data and with computations by other numerical methods. The results showed a favorable agreement of water impact pressure as well as the global fluid motion.

• Ocean Systems Engineering
• /
• v.7 no.3
• /
• pp.299-318
• /
• 2017

The paper aims at illustrating several key issues and ongoing efforts for development of a reliable fully-Lagrangian particle-based solver for simulation of hydroelastic slamming. Fluid model is founded on the solution of Navier-Stokes along with continuity equations via an enhanced version of a projection-based particle method, namely, Moving Particle Semi-implicit (MPS) method. The fluid model is carefully coupled with a structure model on the basis of conservation of linear and angular momenta for an elastic solid. The developed coupled FSI (Fluid-Structure Interaction) solver is applied to simulations of high velocity impact of an elastic aluminum wedge and hydroelastic slammings of marine panels. Validations are made both qualitatively and quantitatively in terms of reproduced pressure as well as structure deformation. Several remaining challenges as well as important key issues are highlighted. At last, a recently developed multi-scale MPS method is incorporated in the developed FSI solver towards enhancement of its adaptivity.

• Journal of Mechanical Science and Technology
• /
• v.18 no.8
• /
• pp.1470-1478
• /
• 2004

In this paper, a new 'presumed' Probability Density Function (PDF) approach coupled with a Lagrangian tracking method is proposed for turbulent combustion modeling. The test and the investigation of the model are conducted by comparing the model results with DNS data for a premixed flame subjected in a decaying turbulent field. The newly constructed PDF, which incorporates the instantaneous chemical reaction term, demonstrates consistent improvement over conventional assumed PDF models. It has been found that the time evolution of the mean scalar, the variance and the mean reaction rate are strongly influenced by a parameter deduced by a Lagrangian equation which takes into account explicitly the local reaction rate. Tests have been performed for a moderate Damkohler number, and it is expected the model may cover a broader range of Damkohler number. The comparison with the DNS data demonstrates that the proposed model may be promising and affordable for implementation in a moment-equation solver.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2009.04a
• /
• pp.264-270
• /
• 2009

선진국 국립연구소에서만 개발하여 적용중인 hydrocode 또는 wavecode를 최근에 국내에서 C++언어로 ExLO를 개발하였다. 3차원 FEM을 적용한 본 코드의 가장 두드러진 특징은 Lagrangian/ALE/Eulerian solver 들을 모두 한 프레임에 내포하고 있어 적용범위에 따라 즉 변형양상에 따라 선택적인 해석 solver의 적용이 가능하다는 것이다. 즉 문제에 따른 (변형의 양상에 따른) 최적의 시뮬레이션 해석 solver 적용할 수 있다. 본 논문에서는 ExLO를 이용하여 고속 물체의 충돌, 파편의 구조물 관통해석 및 대기 중 Air-blast 충격파 전파해석, 물속에서의 수중폭발 충격파 생성 및 전파해석 등의 예제를 소개하고 그 신뢰성을 확인해 보고자 한다. 대체로 군사적인 적용이 많으나 차츰 민간분야의 다양한 방재현상 시뮬레이션에 적용이 가능할 것이다.

• Nuclear Engineering and Technology
• /
• v.56 no.1
• /
• pp.9-18
• /
• 2024

In-vessel retention through external reactor vessel cooling (IVR-ERVC) is a severe accident management (SAM) strategy that has been adopted and used in many nuclear reactors such as AP1000, APR1400, and light water reactor etc. Some reactor accidents have raised concerns about nuclear reactors among residents, leading to a decrease in residents' acceptability and many studies on SAM are being conducted. Experiments on IVR-ERVC are almost impossible due to its specificity, so fluid characteristics are analyzed through BALI experiments with similar condition. In this study, computational fluid dynamics (CFD) via Reynolds-averaged Navier-Stokes (RANS) and large eddy simulation (LES) for BALI experiments were performed. Steady-state CFD analysis was performed on three turbulence models, and SST k-ω model was in good agreement with the experimental measurement temperature within the maximum error range of 1.9%. LES CFD analysis was performed based on the RANS analysis results and it was confirmed that the temperature and wall heat flux for depth was consistent within an error range of 1.0% with BALI experiment. The LES CFD analysis results were compared with those of the Lagrangian-based solver. LES matched the temperature distribution better than SOPHIA, but SOPHIA calculated the position of boundary between stratified layer and convective layer more accurately. On the other hand, Lagrangian-based solver predicted several small eddy behaviors of the convective layer and LES predicted large vortex behavior. The vibration characteristics near the cooling part of the BALI experimental device were confirmed through Fast Fourier Transform (FFT) investigation. It was found that the power spectral density for pressure at least 10 times higher near the side cooling than near the top cooling.

• 한국전산유체공학회:학술대회논문집
• /
• 2010.05a
• /
• pp.71-77
• /
• 2010

As a hybrid model of continuum motion description which combines the advantages of classical kinematical descriptions i.e. Lagrangian and Eulerian description, the ALE (Arbitrary Lagrangian Eulerian) description is adopted for the simulation of a fluid-structure interaction of solid propellant rocket interior. The fluid-structure interaction phenomenon with the deformation of solid domain during the simulation. The developed solver is applied flow and propellant structure. The computed results show complex flow physics in the combustion chamber and the behavior of a solid propellant deformation.

• Journal of the Korea Institute of Military Science and Technology
• /
• v.14 no.5
• /
• pp.932-940
• /
• 2011

This study describes computational simulation results in 2-dimensional and 3-dimensional space concerning large scale gap test(LSGT) by using commercial hydrocode such as AUTODYN and LS-DYNA to analyze the detonation phenomenons of high explosives. To consider the possibilities of LSGT simulation, we used Lee - Tarver reaction rate model of PBX-9404 and Comp-B which were implemented AUTODYN's material library. Also we have tried the diverse numerical schemes such as Lagrangian, Eulerian and ALE(Arbitary Lagrangian Eulerian), SPH(Smoothed Particle Hydrodynamics) in LSGT simulations. After LSGT simulations, we compared the simulation results with published results to verify the LSGT simulations. According to the LSGT simulations, we have concluded as follows. In 2-dimensional and 3-dimensional space, Lagrangian solver provided the most reliable results based on analysis time and accuracy. When using two hydrocodes in 2-dimensional space, the simulation results are almost same except one explosive model. We have verified the modeling method and simulation results of the LSGT by using the commenrcial hydrocode in this study.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.9 no.2
• /
• pp.286-290
• /
• 2008

This paper summarizes the components of an explicit aeroelastic solver developed especially for the simulation of dynamic fracture events occurring during the flight of solid propellant rockets. The numerical method combines an explicit Arbitrary Lagrangian Eulerian (ALE) version of the Cohesive Volumetric Finite Element (CVFE) scheme, used to simulate the spontaneous motion of one or more cracks propagating dynamically through a domain with regressing boundaries, and an explicit unstructured finite volume Euler code to follow the flow field during the failure event. A key feature of the algorithm is the ability to adaptively repair and expand the fluid mesh to handle the large geometrical changes associated with grain deformation and crack motion.

• Transactions of the Korean Society of Automotive Engineers
• /
• v.18 no.5
• /
• pp.108-114
• /
• 2010

The injection nozzle of an electro-hydraulic injector is being opened and closed by movement of a injector's needle which is balanced by pressure at the nozzle seat and at the needle control chamber, at the opposite end of the needle. In this study, the effects of needle movement in a piezo-driven injector on unsteady cavitating flows behavior inside nozzle were investigated by cavitation numerical model based on the Eulerian-Lagrangian approach. Aimed at simulating the 3-D two-phase flow behavior, the three dimensional geometry model along the central cross-section regarding of one injection hole with real design data of a piezo-driven diesel injector has been used to simulate the cavitating flows for injection time by at fully transient simulation with cavitation model. The cavitation model incorporates many of the fundamental physical processes assumed to take place in cavitating flows. The simulations performed were both fully transient and 'pseudo' steady state, even if under steady state boundary conditions. As this research results, we found that it could analyze the effect the pressure drop to the sudden acceleration of fuel, which is due to the fastest response of needle, on the degree of cavitation existed in piezo-driven injector nozzle.

• Proceedings of the KWS Conference
• /
• 2002.05a
• /
• pp.161-163
• /
• 2002

An implicit finite element implementation for Leblond's transformation plasticity constitutive equations, which are widely used in welded steel structure is proposed in the framework of parallel computing. The implementation is based upon the updated Lagrangian formulation. We examine the efficiency of parallel compuatation for the finite element analysis of a welded structure using multi-frontal solver.