• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2008.03a
• /
• 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
• /
• v.53 no.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.

• Journal of the Korean Society for Marine Environment & Energy
• /
• v.10 no.4
• /
• pp.218-226
• /
• 2007

In modern naval ships, the design of submarines has continually evolved to improve survivability and it is also important to design ship against shock response. Exiting underwater ship design has been peformed due to results of static analysis considering shock acceleration by simple method. However, it can not be anticipated good assesment. The present study applied the Arbitrary Lagrangian-Eulerian (ALE) technique, a fluid-structure interaction approach, to simulate an underwater explosion and investigate the survival capability of a damaged submarine liquefied oxygen tank. The Lagrangian-Eulerian coupling algorithm and the equations of state for explosives and seawater were also reviewed. It is shown that underwater explosion analysis using the ALE technique can accurately evaluate structural damage after attack. This procedure could be applied quantitatively to real structural design.

• 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.

• The Bulletin of The Korean Astronomical Society
• /
• v.39 no.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.

• Transactions of Materials Processing
• /
• v.14 no.1 s.73
• /
• 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.10a
• /
• 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.10a
• /
• 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.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.26 no.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.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.30 no.10 s.253
• /
• 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.