• Journal of computational fluids engineering
• /
• v.11 no.1 s.32
• /
• pp.29-35
• /
• 2006

A code is developed to analyze a spherically symmetric underwater explosion. The arbitrary Lagrangian-Eulerian(ALE) Godunov scheme for two-phase flow is used to calculate numerical fluxes through moving control surfaces. For detonation gas of TNT and liquid water, the Jones-Wilkins-Lee(JWL) equation of states and the isentropic Tait relation are used respectively. It is suggested to use the Godunov variable to estimate the velocity of a material interface. The code is validated through comparisons with other results on the gas-water shock tube problem. It is shown that the code can handle generation of discontinuity and recovering of continuity in the normal velocity near the material interface during shock waves interact with the material interface. The developed code is applied to analyze a spherically symmetric underwater explosion. Repeated transmissions of shock waves are clearly captured. The calculated period and maximum radius of detonation gas bubble show good agreements with experimental and other numerical results.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.37 no.1
• /
• pp.1-16
• /
• 2009

In this paper, we introduce two-phase versions of RoeM and AUSMPW+ schemes. Both schemes are originally developed for the gas dynamic problems, and have shown superior accuracy, efficiency and robustness. A new shock discontinuity sensing term is derived from the mixture equation of state, which is commonly used in the RoeM and AUSMPW+ schemes for the stable numerical flux calculation. The developed two-phase versions of the schemes are applied to several liquid-gas, large property discrepancy two-phase test problems, including several shock stability test problems. The results show that both schemes maintain the merits exhibited in gas dynamic problems even in two-phase flows.

• 한국전산유체공학회:학술대회논문집
• /
• 2010.05a
• /
• pp.365-367
• /
• 2010

In this paper, we presented the exact Riemann solver for the air-water two-phase shock tube problems where the strength of the propagated sock wave is moderately weak. The shock tube has a diaphragm in the middle which separates water medium in the left and air medium in the right. By rupturing the diaphragm, various waves such as rarefaction wave, shock wave and contact discontinuity are propagated into water and air. Both fluids are treated as compressible, with the linearized equations of state. We used the isentropic relations for the air and water assuming a weak shock wave. We solved the shock tube problem considering a high pressure in the water and a low pressure in the air. The numerical results cleary showed a left-traveling rarefaction wave in the water, a right-traveling shock wave in the air, and the right-traveling material interface.

• Journal of Aerospace System Engineering
• /
• v.12 no.2
• /
• pp.15-22
• /
• 2018

A compressible multiphase flow was investigated using a DIM consisting of seven equations, including the fifth-order MLP and a modified HLLC Riemann solver to achieve a precise interface structure of liquid and gas. The numerical methods were verified by comparing the flow structures of the high-pressure water and low-pressure air in the shock tube. A 2D air-helium shock-bubble interaction at the incident shock wave condition (Mach number 1.22) was numerically solved and verified using the experimental results.

• Journal of Mechanical Science and Technology
• /
• v.20 no.2
• /
• pp.286-294
• /
• 2006

Flows in a ramjet inlet is simulated for the study of the rocket-ramjet transition. The flow is unsteady, two-dimensional axisymmetric, compressible and turbulent. Double time marching method is used for the unsteady calculation and HLLC method is used as a higher order MUSCL method. As for turbulent calculation, $\kappa-\omega$ SST model is used for more accurate viscous calculations. Sinusoidal pressure perturbation is given at the exit and the flow fields at the inlet is studied. The cruise condition as well as the ground test condition are considered. The pressure level for the ground test condition is relatively low and the effect of the pressure perturbation at the combustion chamber is small. The normal shock at the cruise condition is very sensitive to the pressure perturbation and can be easily detached from the cowl when the exit pressure is relatively high. The sudden decrease in the mass flux is observed when the inlet flow becomes subcritical, which can make the inlet incapable. The amplitude of travelling pressure waves becomes larger as the downstream pressure increases, and the wavelength becomes shorter as Mach number increases. The phase difference of the travelling perturbed pressure wave in space is 180 degree.

• 한국전산유체공학회:학술대회논문집
• /
• 2001.05a
• /
• pp.121-126
• /
• 2001

An accelerating flow field through a compressor cascade is studied numerically by unsteady computational simulation. The two-dimensional Navier-Stokes equations for compressible flow is used for the study of unsteady high incidence angle flow, with preconditioning scheme to cover the wide range of Mach number and $\kappa-\omega$ model for the turbulent viscous flow analysis. A DCA(double circular arc) compressor blade is accelerated artificially in this study to understand the unsteady effect by comparing the present results with the existing steady-state experimental and computational results. Also, the accelerating flow field during the starting phase of gas turbine is studied with actual experimental data for the understanding of flow field and performance characteristics at off-design condition.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2008.11a
• /
• pp.78-81
• /
• 2008

A numerical study of the unsteady flow in an over-expanded thrust optimized contour and compressed truncated perfect rocket nozzle is carried out in present paper. These rocket nozzles are subject to flow separation in transient phase at engine start-up and/or engine shut-down. The separation flow structures at different pressure ratios are observed. The start-up process exhibits two different shock structures such as FSS (Free Shock Separation) and RSS (Restricted Shock Separation). For a range of pressure ratios, hysteresis phenomenon occurs between these two separation patterns. A three-dimension compressible Navier-Stokes solver is used for the present study. One equation Spalart-Allmaras turbulence model is selected. The computed nozzle wall pressures show a good agreement with the experimental measurements. Present results have shown that present code can be used for the analysis of the transient flows in nozzle.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.24 no.9 s.180
• /
• pp.2159-2166
• /
• 2000

Recently, high performance machining center requires special type of sealing mechanism that prevent a leakage of oil jet or oil mist lubrication system. Sealing of oil-air mixture plays important r oles to have an enhanced lubrication for performance machining center. Current work emphasizes on investigations of the air jet effect on the protective collar type labyrinth seal. To improve sealing capabilities of conventional labyrinth seals, air jet is injected against the leakage flow. In this study, an adapted model is introduced to improve sealing capability of conventional non-contact type seals. It has a combined geometry of a protective collar type and an air jet type. Both of a numerical analysis by CFD (Computational Fluid Dynamics) and experimental measurements are carried out to verify sealing improvement. The sealing effects of the leakage clearance and the air jet magnitude aic studied in various parameters. Gas or liquid has been used as a working fluid for most of nori-contact types seals including the labyrinth seal. However, it is more reasonable to regard two-phase flows because oil mist or oil jet are used for high performance spindle's lubrication. In this study, working fluid is regarded as two phases that are mixed flow of oil and air phase. Both of turbulence and compressible flow model are also introduced in a CFD analysis to represent an isentropic process. Estimation of non-leaking property is determined by amount of pressure drop in the leakage path. Results of pressure drop in the experiment match reasonably to those of the simulation by introducing a flow coefficient. Effect of the sealing improvement is explained as decreasing of leakage clearance by air jetting. Thus, sealing effect is improved by amount of air jetting even though clearance becomes larger

• Proceedings of the KSME Conference
• /
• 2007.05b
• /
• pp.2929-2934
• /
• 2007

Prediction methods for cavitation noise are presented. At first, direct numerical simulation of cavitating flow noise has been performed, and acoustic analogy equation based on the cavitation noise modeling is derived. For the direct numerical simulation, a density based homogenous equilibrium model is employed to simulate cavitating two-phase flow and the governing equations are solved with high-order numerical schemes to resolve cavitation noise. The compressible Navier-Stokes equations for mixture fluids are discretized with a sixth-order central compact scheme, and the steep gradient of flow variables and supersonic regions are treated with the selective spatial filtering technique. The direct simulation of cavitating flow noise is performed for a 2D circular cylinder at cavitation number 0.7 and 1. The far-field noise is also predicted with the derived analogy equation. Noise spectrum predicted with the equation is well compared with the result of direct numerical simulation and also agree well with the theory.

• Proceedings of the KSME Conference
• /
• 2003.04a
• /
• pp.1821-1826
• /
• 2003

In this study, we discuss CIP method, which can treat compressible/incompressible problem and multi-phase problem. We can apply this method to the general equations by using CCUP. In this paper, non-staggered grid arrangement and predictor-corrector method are used to enhance later development and the solution accuracy and convergence performance. To validate the numerical algorithm proposed in this paper, the two-dimensional unsteady flow in the driven cavity is simulated. The numerical results of this subject using the present algorithm are compared with other numerical results. As a result, it is prived that the present scheme gives stable and improved solutions, and the results even coarse grid are in good agreement with other result using a fine grid arrangement.