• 한국추진공학회:학술대회논문집
• /
• 한국추진공학회 2008년 영문 학술대회
• /
• pp.659-665
• /
• 2008

The 1-D numerical study of the interior ballistics has been conducted. The unsteady compressible 1-D CFD code using SIMPLER algorithm and QUICK scheme has been developed. The mathematical model of the two-phase flow has been established for the behavior of the interior ballistics. The moving boundary due to the projectile motion as the physical phenomena of the interior ballistics results in the varied control volume. In order to analyze the moving boundary, the numerical codes, which apply the ghost-cell extrapolation method and the Lagrangian method respectively, have been developed. The ghost-cell extrapolation method has been used in the Eulerian coordinate system. The Lagrangian method has been used in Non-Eulerian coordinate system. These codes have been verified through the analysis of the free piston motion problem in the tube. Through this study, the basic techniques of the numerical code for the multi-dimensional two-phase flow of the interior ballistics have been obtained.

• Nuclear Engineering and Technology
• /
• 제53권10호
• /
• pp.3264-3274
• /
• 2021

In a severe accident of light water reactor (LWR), molten core material (corium) can be released into the wet cavity, and a fuel-coolant interaction (FCI) can occur. The molten jet with high speed is broken and fragmented into small debris, which may cause a steam explosion or a molten core concrete interaction (MCCI). Since the premixing stage where the jet breakup occurs has a large impact on the severe accident progression, the understanding and evaluation of the jet breakup phenomenon are highly important. Therefore, in this study, the jet breakup simulations were performed using the Smoothed Particle Hydrodynamics (SPH) method which is a particle-based Lagrangian numerical method. For the multi-fluid system, the normalized density approach and improved surface tension model (CSF) were applied to the in-house SPH code (single GPU-based SOPHIA code) to improve the calculation accuracy at the interface of fluids. The jet breakup simulations were conducted in two cases: (1) jet breakup without structures, and (2) jet breakup with structures (control rod guide tubes). The penetration depth of the jet and jet breakup length were compared with those of the reference experiments, and these SPH simulation results are qualitatively and quantitatively consistent with the experiments.

• 한국해양환경ㆍ에너지학회지
• /
• 제10권4호
• /
• pp.218-226
• /
• 2007

본 연구는 생존성 측면에서 잠수함의 내충격 설계시 고려해야 하는 수중폭발 충격응답해석 기법에 관하여 고찰하였다. 생존성을 고려하는 것은 현대의 모든 함정설계에 필수요건이며, 특히 잠수함의 수중폭발에 대한 내충격 설계는 반드시 고려되어야 할 사항이다. 기존의 수중함에 대한 내충격 설계는 가정된 이론에 의한 충격가속도를 정적해석의 결과를 통하여 단순하게 적용하였으나 본 연구에서는 직접해석법에 수중폭발 충격응답해석을 수행함으로써 기존의 방법보다 매우 합리적이고 신뢰성 높은 결과를 제시하였다. 특히 본 연구에서는 직접해석법 중에서 기존 수상함에서 널리 사용하고 있는 LS-DYNA/USA code일 적용보다 LS-DYNA code만으로 사용자가 보다 손쉽게 접근할 수 있는 ALE기법을 이용하여 잠수함 액화산소탱크를 대상으로 수중폭발해석 방식을 제안하고자 하였다.

• Nuclear Engineering and Technology
• /
• 제56권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.

• 천문학회보
• /
• 제39권2호
• /
• pp.65.2-65.2
• /
• 2014

We present a new hydrodynamic simulation code based on the Voronoi tessellation for estimating the density precisely. The code employs both of Lagrangian and Eulerian description by adopting the movable mesh scheme, which is superior to the conventional SPH (smoothed particle hydrodynamics) and AMR (adaptive mesh refinement) schemes. The code first generates unstructured meshes by the Voronoi tessellation at every time step, and then solves the Riemann problem for all surfaces of each Voronoi cell so as to update the hydrodynamic states as well as to move current meshes. Besides, the IEM (incremental expanding method) is devised to compute the Voronoi tessellation to desired degree of speed, thereby the CPU time is turned out to be just proportional to the number of particles, i.e., O(N). We discuss the applications of our code in the context of cosmological simulations as well as numerical experiments for galaxy formation.

• 소성∙가공
• /
• 제14권1호
• /
• pp.43-48
• /
• 2005

The production of metal strip with uniform thickness and flatness are two important problems associated thin strip rolling. The thickness and flatness of strip are affected by the flattening of contact surface between strip and roll, the elastic recovery and the bending of roll. Especially, the flatness of the strip is greatly affected by bending deflection of roll. The roll must be designed considered the elastic deformation of roll. This study describes the measurement of thickness and flatness of strip and shows the crown roll for producing flat strip. But it is difficult to produce the crown roller. The cross rolling that is a simple method which can produce the flat strip is introduced and it is found the optimal cross angle for improvement of flatness of plate. These problems are solved by the MARC code on the basis of elastic-plastic material and the updated Lagrangian formulation.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2003년도 The Fifth Asian Computational Fluid Dynamics Conference
• /
• pp.233-234
• /
• 2003

Fly ash enters axial compressor when a turbomachinery is operated in an adverse environment. We have numerically investigated erosion of the blade and shroud in the turbulent compressor passage flow under the influence of gas-particle two-phase interaction. There have appeared quasi-three dimensional calculations on this subject but not the complete three-dimensional gas-particle interaction as done in the present work. Lagrangian particle tracing technique is used on the base of parallel processing for efficient calculation. Accuracy of the present code is tested using the benchmark lPL nozzle. In the DFVLR compressor blades, we have shown that a large number of particles passing through the tip clearance make impact on the blade tip and on the shroud. Higher degree of erosion is resulted by the heavier particles due to the centrifugal force.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2003년도 The Fifth Asian Computational Fluid Dynamics Conference
• /
• pp.203-204
• /
• 2003

Fly ash enters axial compressor when a turbomachinery is operated in an adverse environment. We have numerically investigated erosion of the blade and shroud in the turbulent compressor passage flow under the influence of gas-particle two-phase interaction. There have appeared quasi-three dimensional calculations on this subject but not the complete three-dimensional gas-particle interaction as done in the present work. Lagrangian particle tracing technique is used on the base of parallel processing for efficient calculation. Accuracy of the present code is tested using the benchmark JPL nozzle. In the DFVLR compressor blades, we have shown that a large number of particles passing through the tip clearance make impact on the blade tip and on the shroud. Higher degree of erosion is resulted by the heavier particles due to the centrifugal force.

• 대한기계학회논문집A
• /
• 제26권3호
• /
• pp.576-585
• /
• 2002

A three-dimensional Lagrangian explicit time-integration finite element code for analyzing the dynamic impact phenomena was developed. It uses four node tetrahedral elements. In order to consider the effects of strain rate hardening, strain hardening and thermal softening, which are frequently observed in high-velocity deformation phenomena, Johnson-Cook model is used as constitutive model. For more accurate and robust contact force computation, the defense node contact algorithm was adopted and implemented. In order to evaluate the performance of the newly developed three-dimensional hydrocode NET3D, numerical simulations of the oblique impact of mild steel plate by mild steel sphere were carried out. Ballistic limit about various oblique angle between 0 degree and 80 degree was estimated through a series of simulations with different initial velocities of sphere. Element eroding by equivalent plastic strain was applied to mild steel spheres and targets. Ballistic limits and fracture characteristics obtained from simulation were compared with experimental results conducted by Finnegan et al. From numerical studies, the following conclusions were reached. (1) Simulations could successfully reproduce the key features observed in experiment such as tensile failure termed "disking"at normal impacts and outwards bending of partially formed plus segments termed "hinge-mode"at oblique impacts. (2) Simulation results fur 60 degrees oblique impact at 0.70 km/s and 0.91 km/s were compared with experimental results and Eulerian hydrocode CTH simulation results. The Lagrangian code NET3D is superior to Eulerian code CTH in the computational accuracy. Agreement with the experimentally obtained final deformed cross-sections of the projectile is excellent. (3) Agreement with the experimental ballistic limit data, particularly at the high-obliquity impacts, is reasonably good. (4) The simulation result is not very sensitive to eroding condition but slightly influenced by friction coefficient.

• 대한기계학회논문집A
• /
• 제30권10호
• /
• pp.1320-1327
• /
• 2006

Parallel procedure and performance of two-step Eulerian code have not been reported sufficiently yet even though it was developed and utilized widely in the impact simulation. In this study, parallel strategy of two-step Eulerian code was proposed and described in detail. The performance was evaluated in the self-made linux cluster computer. Compared with commercial code, a relatively good performance is achieved. Through the performance evaluation of each computation stage, remap is turned out to be the most time consuming part among the other part such as FE processing, communication, time marching etc.