• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2012.05a
• /
• pp.187-195
• /
• 2012

This paper deals with numerical computations of cryogenic flows around turbopump inducer. Firstly, we introduce numerical methods to compute compressible/incompressible cryogenic two-phase flow. As a validation problem, computation results of 2 dimensional/axi-symmetric cryogenic flow will be presented. In this process, various cavitation model will be compared. Finally, numerical simulation of 3 dimensional turbopump inducer will be presented.

• 한국전산유체공학회:학술대회논문집
• /
• 2005.10a
• /
• pp.169-177
• /
• 2005

A high resolution scheme for solving gas-liquid two-phase flows with cavitation is described. This scheme uses the curvilinear coordinate grid and solves the density based momentum equations for mixture of gas-liquid medium with a preconditioning method to treat both compressible and incompressible flow characteristics. The present preconditioned method is based on the Runge-Kutta explicit finite-difference scheme, and is improved by using the diagonalization, the flux difference splitting and the MUSCL-TVD schemes to save computational effort and to increase stability and resolvability, especially at gas-liquid contact surfaces. A homogeneous equilibrium cavitation model is used to treat the gas-liquid two-phase medium in cavitating flow as a locally homogeneous pseudo-single-phase medium. Therefore, it is easy to solve cavitating flow, including wave propagation, large density changes and incompressible flow characteristic at low Mach number. Some numerical results obtained by the present scheme are shown.

• International Journal of Fluid Machinery and Systems
• /
• v.9 no.1
• /
• pp.66-74
• /
• 2016

Concerning the numerical simulation of high-speed water jet with intensive cavitation this paper presents a practical compressible mixture flow method by coupling a simplified estimation of bubble cavitation and a compressible mixture flow computation. The mean flow of two-phase mixture is calculated by URANS for compressible fluid. The intensity of cavitation in a local field is evaluated by the volume fraction of gas phase varying with the mean flow, and the effect of cavitation on the flow turbulence is considered by applying a density correction to the evaluation of eddy viscosity. High-speed submerged water jets issuing from a sheathed sharp-edge orifice nozzle are treated when the cavitation number, ${\sigma}=0.1$, and the computation result is compared with experimental data The result reveals that cavitation occurs initially at the entrance of orifice and bubble cloud develops gradually while flowing downstream along the shear layer. Developed bubble cloud breaks up and then sheds downstream periodically near the sheath exit. The pattern of cavitation cloud shedding evaluated by simulation agrees experimental one, and the possibility to capture the unsteadily shedding of cavitation clouds is demonstrated. The decay of core velocity in cavitating jet is delayed greatly compared to that in no-activation jet, and the effect of the nozzle sheath is demonstrated.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.17 no.3
• /
• pp.209-219
• /
• 2013

The present paper deals with the extension of AUSMPW+ scheme into two-fluid model for multiphase flow. AUSMPW+ scheme is the improvement of a single-phase AUSM+ scheme by designing pressure-based weighting functions to prevent oscillations near a wall and shock instability after a strong shock. Recently, Kitamura and Liou assessed a family of AUSM-type schemes with two-fluid model governing equations [K. Kitamura and M.-S. Liou, Comparative study of AUSM-Family schemes in compressible multi-phase flow simulations, ICCFD7-3702 (2012)]. It was observed that the direct application of the single-phase AUSMPW+ did not provide satisfactory results for most of numerical test cases, which motivates the current study. It turns out that, by designing pressure-based weighting functions, which play a key role in controlling numerical diffusion for two-fluid model, problems reported in can be overcome. Various numerical experiments validate the proposed modification of AUSMPW+ scheme is accurate and robust to solve multiphase flow within the framework of two-fluid model.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2009.11a
• /
• pp.221-224
• /
• 2009

Since the coupling of cavitation modeling with turbulent flow is the difficulty topic, a numerical simulation for two phase flow remains as one of the challenging issues in the society. This research focuses on the development of numerical code to deal with incompressible two phase flow around conical body combined with cavitation model suggested by Kunz et al. with k-e turbulent model. The simulation results are compared to experimental data to verify the validity of the developed code. The calculation results show very good agreement with experimental observations. Also, the calculation of cavitation in cryogenic fluid is being done by implementing the temperature sensitivity in government equations and it is still in the progress. This code have been being further extended to 3D compressible two phase flow for the study on the fluid dynamics around inducers and impellers in turbo pump system.

• Journal of the Society of Naval Architects of Korea
• /
• v.42 no.4 s.142
• /
• pp.330-340
• /
• 2005

Compressible two-phase flow is analyzed based on the arbitrary Lagrangian-Eulerian (ALE) formulation. For water, Tamman type stiffened equation of state is used. Numerical fluxes are calculated using the ALE two-phase Godunov scheme which assumes only that the speed of sound and pressure can be provided whenever density and internal energy are given. Effects of the approximations of a material interface speed are Investigated h method Is suggested to assign a rigid body boundary condition effectively To validate the developed code, several well-known problems are calculated and the results are compared with analytic or other numerical solutions including a single material Sod shock tube problem and a gas/water shock tube problem The code is applied to analyze the refraction and transmission of shock waves which are impacting on a water-gas interface from gas or water medium.

• 한국전산유체공학회:학술대회논문집
• /
• 2007.10a
• /
• pp.249-257
• /
• 2007

A high resolution numerical method aimed at solving gas-liquid two-phase flow is proposed and applied to gas-liquid two-phase shock tube problem. The present method employs a finite-difference 4th-order Runge-Kutta method and Roe's flux difference splitting approximation with the MUSCL TVD scheme. By applying the homogeneous equilibrium cavitation model, the present density-based numerical method permits simple treatment of the whole gas-liquid two-phase flow field, including wave propagation and large density changes. The speed of sound for gas-liquid two-phase media is derived on the basis of thermodynamic relations and compared with that by eigenvalues. By this method, a Riemann problem for Euler equations of one dimensional shock tube was computed. Numerical results such as detailed observations of shock and expansion wave propagations through the gas-liquid two-phase media and some data related to computational efficiency are made. Comparisons of predicted results and exact solutions are provided and discussed.

• Journal of computational fluids engineering
• /
• v.15 no.4
• /
• pp.25-31
• /
• 2010

Interface tracking of two phase is significant to analyze multi-phase phenomena. The VOF(Volume of Fluid) and level set are well known interface tracking method. However, they have limitations to solve compressible flow and incompressible flow at the same time. CIP(Cubic Interpolate Propagation) method is appropriate for considering compressible and incompressible flow at once by solving the governing equation which is divided up into advection and non-advection term. In this article, we analyze the droplet impingement according to various We number using improved CIP method which treats nonlinear term once more comparison with original CIP method. Furthermore, we compare spread radius after droplet impingement on the wall with the experimental data and original CIP method. The result using improved CIP method shows the better result of the experiments, comparison with result of original CIP method, and it reduces the mass conservation error which is generated in the numerical analysis comparison with original CIP method.

• Journal of the Korean Society for Precision Engineering
• /
• v.14 no.4
• /
• pp.56-63
• /
• 1997

Sealing of lubricant-air mixture in the high performance machining center is one of most the important characteristics to carry out enhanced lubrication. High speed spindle requires non-contact type of sealing mechanism. Evaluating an optimum seal design to minimize leakage is concerned in the aspect of flow control. Effect of geometry and leakage path are evaluated according to variation of sealing geometry. Velocity, pressure, turbulence intensity of profile is calculated to find more efficient geometry and variables. This offers a methodological way of enhancement seal design for high speed spindle. The working fluid is regarded as two phases that are mixed flow of oil phase and air phase. It is more reasonable to simulate an oil jet or oil mist type high speed spindle lubrication. Turbulence and compressible flow model are used to evaluate a flow characteristic. This paper considers a design effect of sealing capability of non- contact type seals for high speed spindle and analyzes leakage characteristics to minimize a leakage 7 on the same sealing area.

• Proceedings of the Korean Society of Tribologists and Lubrication Engineers Conference
• /
• 1997.04a
• /
• pp.176-180
• /
• 1997

Evaluating an optimum seal design to minimize leakage is concerned in the aspect of flow control. Flow is characterized into five categories according to its leakage path. Effect of geometry and leakage path are evaluated according to variation of sealing geometry. To simulate an oil jet or oil mist type high speed spindle lubrication, the working fluid is regarded as two phases that are mixed flow of oil phase and air phase. Both of the turbulence and the compressible flow model were introduced in CFD(Computational Fluid Dynamics) analysis. This offers a methodological way of enhancement seal design for high speed spindle.