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

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.21 no.1
• /
• pp.24-35
• /
• 1997

Melting process inside in a vertical cylinder has been investigated numerically to observe heat transfer characteristics in the latent heat storage vessel applied to the thermal storage system. The time-dependent boundary fitted coordinate system was introduced to overcome the difficulty caused by the moving boundary. The present results are in good agreement with the available previous data when the initial subcooling effect of the solid phase is not considered. It is found that the melting is promoted by the natural convection, but is delayed by the initial subcooling effect of the solid phase.

• Journal of Advanced Marine Engineering and Technology
• /
• v.20 no.3
• /
• pp.162-170
• /
• 1996

A numerical study on the Stefan problem occurred in cryosurgery is performed. Crank-Nicholson type finite difference algorithm based on the enthaly method is adapted to solve the phase change problem in this study. As it is a moving boundary problem, special emphasis is put on the estimation of the freezing front location. Two cases selected here are freezings of human tissue by disk type cryoprobe and by hemispherical one. In both cases, the heat flows are considered to be one dimensional. The calculated results using enthalpy method are compared with those using the program TRUMP and with Neumann's solution. These results agree guite well with each other. While it is pretty difficult to get accurate freezing front location by TRUMP due to the so- called "phase change knee" occured during the phase change, the algorithm based on the enthalpy method is proved to be very powerful to cope with this kind of problem.f problem.

• Proceedings of the Korean Nuclear Society Conference
• /
• 1995.05b
• /
• pp.643-650
• /
• 1995

A numerical method for treating transient diffusion Involving change of phase is presented. In other methods of dealing with this class of problems, the mass flux balance at the moving phase boundary requires explicit treatment of two distinct phases. The technique, originating from the apparent heat capacity method in transient heat conduction with the phase change, avoids the difficulty by transferring the concentration discontinuity at the boundary to smoothed physical property variations near the moving front. This technique accomodates the nonlinearities which preclude use of analytical solutions. It was tested against known analytical solutions for simple cases and turned out to be quite accurate.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.10 no.4
• /
• pp.585-592
• /
• 1986

본 연구에서는 시간의존 Boundary-Fitted 좌표계로의 좌표변환을 이용하여 자 연대류를 고려한 임의 형상에서의 2차원 상변화 문제를 해석하는 안정성 있는 수치해 석 방법을 제시하고자 한다.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.20 no.1
• /
• pp.83-106
• /
• 2016

This paper reviews and compares three different methods for modeling incompressible and immiscible ternary fluid flows: the immersed boundary, level set, and phase-field methods. The immersed boundary method represents the moving interface by tracking the Lagrangian particles. In the level set method, an interface is defined implicitly by using the signed distance function, and its evolution is governed by a transport equation. In the phase-field method, the advective Cahn-Hilliard equation is used as the evolution equation, and its order parameter also implicitly defines an interface. Each method has its merits and demerits. We perform the several simulations under different conditions to examine the merits and demerits of each method. Based on the results, we determine the most suitable method depending on the specific modeling needs of different situations.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.29 no.4 s.235
• /
• pp.425-434
• /
• 2005

In order to analyze pulsating flows in elastic blood vessels, a method based on the ALE concept and finite volume method was reformed and modulated to include wall motion of elastic vessels and impedance phase angle(phase difference between wall motion and blood flow). Our study indicated wall shear rates(WSR) were significantly influenced by the wall motion and the impedance phase angle. For larger wall motion more than $5{\%},$ the computed WSR started to deviate from the results of the perturbation theory that assumed smaller wall motion. The study showed that oscillatory shear index increased as the impedance phase angle became more negative like $-70{\circ}\;or\;-80{\circ}$ due to reduced mean WSR and increased amplitude of WSR. This result indicated that hypertensive patients are more vulnerable to atherosclerosis than normal persons because of the role of more negative impedance phase angles usually observed in these patients.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2001.05a
• /
• pp.385-392
• /
• 2001

Time-dependent frequency and energy of free vibration of the Spagetti problem, that is the axially moving continuum with time-varying length, are investigated. Exact expressions for the natural frequency and time-varying vibration energy are derived by dealing with traveling waves. When the string length is increased, the vibration period increases, but the free vibration energy varies as a function of both translating velocity and boundary velocity of the continuum. However, when the string undergoes retraction, the vibration energy increases with time, String tension together with non-zero instantaneous velocity at the moving boundary results in energy variation.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.16 no.12
• /
• pp.2385-2394
• /
• 1992

Granular solid propellants having energy and fast burning rate produce great thrusts within extremely short time intervals. Thus numerical researchs prevailed rather than experimental. Using a 2-phase fluid dynamics model among 1-dimensional 2-phase models, a numerical program was set up to describe reacting flow fields, moving boundary with oscillating pressure waves and constitutive laws research. It deserves special emphasis that correlations of convective heat transfer coefficient and viscous drag force among constitutive laws are tested and discussed because slight variations of their constants make a large influence on their results. In this calculations, some of correlations make the large difference in results. Therefore constitutive laws for convective heat transfer coefficient and viscous drag force need more considerations with experiments.

• Journal of Mechanical Science and Technology
• /
• v.20 no.6
• /
• pp.849-863
• /
• 2006

IBM (Immersed Boundary Method) with feedback momentum forcing was applied to stationary and moving bodies. The capability of IBM to treat the obstacle surfaces, especially with moving effect has been tested for two dimensional problems. Stationary and oscillating cylinders were simulated by using IBM based on finite volume method with Cartesian coordinates. For oscillating cylinder, lateral and vertical motions are considered, respectively. Present results such as time histories of drag and lift coefficients for both stationary and oscillating cases are in good agreement with previous numerical and experimental results. Also, the instantaneous wake patterns of oscillating cylinder with different oscillating frequency ratios well represented those of previous researches. More feasibility study for IBM has been carried out to two oscillating cylinders. Drag and lift coefficients are presented for two cylinders oscillating sinusoidally with phase difference of $180^{\circ}$.