• 화약ㆍ발파
• /
• 제34권2호
• /
• pp.31-38
• /
• 2016

유체 동역학 코드는 고속 충돌을 모델링하는 수치해석 툴로서 재료가 유체처럼 거동한다고 가정하며, 화약을 이용한 암반발파와 같은 충돌 문제를 푸는 데 광범위하게 사용된다. 암반발파를 현실적으로 모사하기 위해서는 화약을 수치해석적으로 모델링할 필요가 있으며, 이를 통해 암반과 화약의 상호작용 문제를 완전 연계된 방식으로 풀 수 있다. 화약을 수치 모델링하기 위해서는 특정 물리적 조건에서 재료의 상태를 나타내는 상태 방정식이 수립되어야 한다. 본 고에서는 발파 과정을 수치 모델링하기 위한 유체 동역학 코드, 화약의 상태 방정식과 관련 매개변수의 결정방법에 대해 소개하였다.

• 천문학회지
• /
• 제27권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.

• 천문학회보
• /
• 제45권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.

• 한국항공우주학회지
• /
• 제43권1호
• /
• pp.8-22
• /
• 2015

지금까지 인젝터의 수치적 시뮬레이션은 대부분 Eulerian 기법의 바탕위에서 이루어져 왔다. 그러나 액체제트의 미립화현상과 복잡한 공기와의 경계면 변화를 나타내는데 있어 기존의 기법들이 갖는 선천적인 단점이 존재하며 따라서 본 연구에서는 새로운 Smoothed Particle Hydrodynamics(SPH)라는 입자 기법을 도입하였다. 수치적 시뮬레이션을 위해 먼저 해석을 위한 SPH 코드를 개발하였으며 본 논문에서는 인젝터 문제를 정확하게 나타내는데 있어 필수적인 알고리즘중 하나인 다상유동모사에 대한 검증문제가 제시 되어 있다. 마지막으로 다양한 인젝터 종류 중 하나인 액체-액체 동축형 스월 인젝터에 대한 시뮬레이션이 수행되었으며 실제실험과의 비교를 진행하였다.

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

• 한국연소학회:학술대회논문집
• /
• 한국연소학회 2012년도 제45회 KOSCO SYMPOSIUM 초록집
• /
• 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.

• 천문학회지
• /
• 제34권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.

• 천문학회지
• /
• 제29권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.

• 한국군사과학기술학회지
• /
• 제14권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.

• 천문학회지
• /
• 제34권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.