• Explosives and Blasting
• /
• v.34 no.2
• /
• pp.31-38
• /
• 2016

The hydrodynamics code is a numerical tool developed for modeling high velocity impacts where the materials are assumed to behave like fluids. The hydrodynamics code is widely used for solving impact problems, such as rock blasting using explosives. For a realistic simulation of rock blasting, it is necessary to model explosives numerically so that the interaction problem between rock and explosives can be solved in a fully coupled manner. The equation of state of explosives, which describes the state of the material under given physical conditions, should be established. In this paper, we introduced the hydrodynamics code used for explosion process modeling, the equation of state of explosives, and the determination of associated parameters.

• Journal of The Korean Astronomical Society
• /
• v.27 no.1
• /
• pp.13-29
• /
• 1994

We have constructed a 3-dim hydrodynamics code called BTSPH. The fluid dynamics part of the code is based on the smoothed particle hydrodynamics (SPH), and for its Poisson solver the binary tree (BT) scheme is employed. We let the smoothing length in the SPH algorithm vary with space and time, so that resolution of the calculation is considerably enhanced over the version of SPH with fixed smoothing length. The binary tree scheme calculates the gravitational force at a point by collecting the monopole forces from neighboring particles and the multipole forces from aggregates of distant particles. The BTSPH is free from geometric constraints, does not rely on grids, and needs arrays of moderate size. With the code we have run the following set of test calculations: one-dim shock tube, adiabatic collapse of an isothermal cloud, small oscillation of an equilibrium polytrope of index 3/2, and tidal encounter of the polytrope and a point mass perturber. Results of the tests confirmed the code performance.

• The Bulletin of The Korean Astronomical Society
• /
• v.45 no.1
• /
• pp.55.1-55.1
• /
• 2020

In an attempt to investigate the nonlinear dynamics such as shock, shear, and turbulence associated with ultra-relativistic jets, we develop a new relativistic hydrodynamics (RHD) code based on the weighted essentially non-oscillatory (WENO) scheme. It is a 5th-order accurate, finite-difference scheme, which has been widely used for solving hyperbolic systems of conservation equations. The code is parallelized with MPI and OpenMP. Through an extensive set of tests, the accuracy and efficiency of different WENO reconstructions, and different time discretizations are assessed. Different implementations of the equation of state (EOS) for relativistic fluid are incorporated, As the fiducial setup for simulations of ultra-relativistic jets, we adopt the EOS in Ryu et al. (2006) to treat arbitrary adiabatic index of relativistic fluid, the WENO-Z reconstructions to minimize numerical dissipation without loss of stability, and the strong stability preserving Runge-Kutta (SSPRK) method to achieve stable time stepping with large CFL numbers. In addition, the code includes a high-order flux averaging along the transverse directions for multi-dimensional problems, and the modified eigenvalues for the acoustic modes to effectively control the carbuncle instability. We find that the new code performs satisfactorily simulations of ultra-relativistic jets.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.43 no.1
• /
• pp.8-22
• /
• 2015

Most jet spray and atomization simulations are done with the Eulerian method which has inherent disadvantage in representing jet breakups and droplets. Full Lagrangian particles method called Smoothed Particle Hydrodynamics(SPH) is used in this work. We develop the SPH code and perform validations that confirm the suitability of our SPH method for simulating liquid jet atomization problem. Then, we conduct the simulation of liquid-liquid swirl coaxial injector for comparison against the experimental data.

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

• 한국연소학회:학술대회논문집
• /
• 2012.11a
• /
• pp.15-16
• /
• 2012

The present study aims to investigate quantitatively the hydrodynamic effect on combustion process of pulverized coal particles in large scale combustion chamber using computational analysis, with a general purpose computational fluid dynamics code. Burner hydrodynamics include swirl and turbulence intensity from the burner. To understand the phenomena which are difficult to observe how flow has influence on the combustion process, comparative effect of combustion related coal properties and hydrodynamics is evaluated on flame formation and development in burner flames.

• Journal of The Korean Astronomical Society
• /
• v.34 no.4
• /
• pp.251-254
• /
• 2001

To examine the structure and dynamics of thick accretion disks, we use a two-dimensional viscous hydrodynamic code coupled with radiation transport. The $\alpha$-model and the full viscous stress-tensor description for the kinematic viscosity are used. The radiation transport is treated in the gray, flux-limited diffusion approximation. The finite difference methods used are based on an explicit-implicit method. We apply the numerical code to the Super-Eddington black-hole model for SS 433.@The result for a very small viscosity parameter a reproduces well the characteristic features of SS 433, such as the relativistic jets with $\~$0.26c, the small collimation degree of the jets, the mass-outflow rate of ${\ge}5{\times}10^{-7}M{\bigodot}yr^{-1}$, and the formation of the X-ray iron emission lines.

• Journal of The Korean Astronomical Society
• /
• v.29 no.2
• /
• pp.195-205
• /
• 1996

The evolution of the stellar debris after tidal disruption due to the super massive black hole's tidal force is difficult to solve numerically because of the large dynamical range of the problem. We developed an SPH (Smoothed Particle Hydrodynamics) - TVD (Total Variation Diminishing) hybrid code in which the SPH is used to cover a widely spread debris and the TVD is used to compute the stream collision more accurately. While the code in the present form is not sufficient to obtain desired resoultion, it could provide a useful tool in studying the aftermath of the stellar disruption by a massive black hole.

• Journal of the Korea Institute of Military Science and Technology
• /
• v.14 no.3
• /
• pp.532-537
• /
• 2011

This paper describes the development of SPH(Smooth Particle Hydrodynamics) scheme and integration into the multi-material shock physics code(ExLO) for the purpose of the application to the extreme large deformation problems. SPH numerical scheme has been extended into the fluid dynamics and the high-speed impact events, such as space structure protection against space debris and meteorite catering. Like other hydrocodes, SPH scheme also solves the conservation equations with the constitutive equation including equation of state. The benchmark problem, Taylor-Impact test, was simulated and the predictions show good agreements with both the published numerical data and experimental data. Currently, the contact treatment between materials is under development.

• Journal of The Korean Astronomical Society
• /
• v.34 no.4
• /
• pp.255-259
• /
• 2001

PMDSPH is a combined 3D particle-mesh and SPH code aimed to simulate the self-consistent dynamical evolution of spiral galaxies including live stellar and collisionless dark matter components, as well as an isothermal gas component. This paper describes some aspects of this code and shows how its application to the Milky Way helps to recover the gas flow within the Galactic bar region from the observed HI and CO longitude-velocity distributions.