• 한국전산유체공학회:학술대회논문집
• /
• 2008.10a
• /
• pp.327-333
• /
• 2008

The cavitating flow simulation is of practical importance for many engineering systems, such as marine propellers, pump impellers, nozzles, injectors, torpedoes, etc. The present work has developed a base code for simulating cavitating flows past cylinders and hydrofoils. The governing equation is the Navier-Stokes equation based on homogeneous mixture model. The momentum and energy equation is in the mixture phase while the continuity equation is solved in liquid and vapor phase, separately. The solver employs an implicit preconditioning algorithm in curvilinear coordinates. The computations have been carried out for the cylinders with spherical, 1- and 0-caliber forebody and hydrofoil of ALE and NACA cross-section and, then, compared with experiments and other numerical results. Fairly good agreements with experiments and numerical results have been achieved. The present base code has shown the feasibility to solve the cavitating flow past supercavitating torpedo after the improvement for compressibility effects and interactions with hot exhaust gas of propulsive rocket.

• 한국전산유체공학회:학술대회논문집
• /
• 2002.05a
• /
• pp.146-151
• /
• 2002

Characteristics of flow rate control has been studied for a cavitating venturi adopted in a liquid rocket propellant feed system. Numerical simulation has been peformed to give about $10\%$ discrepancy of mass flow rate to the experimental data for cavitating flow regime. Mass flow rate is confirmed to be saturated for pressure difference higher than $3\times10^5$pa when the upstream pressure is fixed to $22.8\times10^5$pa and the downstream pressure is varied. The evaporation amount depends substantially to non-condensable gas concentration. However the mass flow rate characteristic is relatively insensitive to the mass fraction of non-condensable gas. So it is reduced by only $2\%$ when the non- condensable gas concentration is increased from 1.5PPM to 150PPM. From the previous comparison the expansions of the non-condensable gas and the evaporation of liquid are verified to have same effect to pressure recovery.

• 한국전산유체공학회:학술대회논문집
• /
• 2005.04a
• /
• pp.219-222
• /
• 2005

A compressible two-fluid two-phase flow computation model using the stiffened-gas equation of state is formulated. Since the conservation equation system is of mixed type, it gives complex eigenvalues. The sonic speeds obtained from the individual single phase have been simply used in the literature for the fastest wave speeds necessary in the HLL scheme. This method has worked fine but proved to be quite diffusive according to our test. To improve the accuracy, we here propose to utilize the analytic eigenvalues evaluated from an approximate Jacobian matrix lot the fastest wave speeds. The interfacial transfer terms were dropped in constituting the Jacobian matrix for this purpose. The present scheme proved efficient, robust and accurate in comparison with other existing methods. We solved the cavitating flow problem using the present scheme. The result shows more detailed wave structure in the cavitating process caused by the strong expansion waves.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.4 no.3
• /
• pp.256-266
• /
• 2012

Cavitating flow simulation is of practical importance for many engineering systems, such as marine propellers, pump impellers, nozzles, torpedoes, etc. The present work has developed the base code to solve the cavitating flows past the axisymmetric bodies with several forebody shapes. The governing equation is the Navier-Stokes equation based on homogeneous mixture model. The momentum is in the mixture phase while the continuity equation is solved in liquid and vapor phase, separately. The solver employs an implicit preconditioning algorithm in curvilinear coordinates. The computations have been carried out for the cylinders with hemispherical, 1-caliber, and 0-caliber forebody and, then, compared with experiments and other numerical results. Fairly good agreements with experiments and numerical results have been achieved. It has been concluded that the present numerical code has successfully accounted for the cavitating flows past axisymmetric bodies. The present code has also shown the capability to simulate ventilated cavitation.

• Journal of the Society of Naval Architects of Korea
• /
• v.54 no.3
• /
• pp.215-226
• /
• 2017

This paper presents numerical results of the performance of a marin propeller in cavitating and non-cavitating flow conditions. The geometry and experimental validation data of the propeller are provided in Potsdam Propeller Test Case(PPTC) in the framework of the second International Symposium on Marine Propulsors 2011(SMP'11) workshop. The PPTC includes open water tests, velocity field measurements and cavitation tests. The present numerical analysis was carried out by using the Reynolds averaged Navier-Stokes(RANS) method on a wall-resolved grid ensuring a y+=1, where the SST k-${\omega}$ model was mainly used for turbulence closure. The influence of the turbulence model was investigated in the prediction of the wake field under a non-cavitating flow condition. The propeller tip vortex flows in both cavitating and non-cavitating conditions were captured through adaptation of additional grids. For the cavitation flows at three operation points, Schnerr-Sauer's cavitation model was used with a Volume-Of Fluid(VOF) approach to capture the two-phase flows. The present numerical results for the propeller wake and cavitation predictions including the open water performance showed a qualitatively reasonable agreement with the model test results.

• Journal of computational fluids engineering
• /
• v.16 no.1
• /
• pp.14-21
• /
• 2011

Super-cavitating flows around under-water bodies are being studied for drag reduction and dramatic speed increase. In this paper, high speed super-cavitating flow around a two-dimensional symmetric wedge-shaped body were studied using an unsteady Reynolds-averaged Navier-Stokes equations solver based on a cell-centered finite volume method. To verify the computational method, flow over a hemispherical head-form body was simulated and validated against existing experimental data. Various computational conditions, such as different wedge angles and caviation numbers, were considered for the super-cavitating flow around the wedge-shaped body. Super-cavity begins to form in the low pressure region and propagates along the wedge body. The computed cavity lengths and velocities on the cavity boundary with varying cavitation number were validated by comparing with analytic solution.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.2 no.1
• /
• pp.39-44
• /
• 2010

Recently underwater systems moving at high speed such as a super-cavitating torpedo have been studied for their practical advantage of the dramatic drag reduction. In this study we are focusing our attention on super-cavitating flows around axisymmetric cavitators. A numerical method based on inviscid flow is developed and the results for several shapes of the cavitator are presented. First using a potential based boundary element method, we find the shape of the cavtiator yielding a sufficiently large enough cavity to surround the body. Second, numerical predictions of supercavity are validated by comparing, with experimental observations carried out in a high speed cavitation tunnel at Chungnam National University (CNU CT).

• Journal of the Korean Society of Propulsion Engineers
• /
• v.18 no.3
• /
• pp.84-91
• /
• 2014

The flow rate control of a cavitating venturi has been investigated with downstream pressure variation. A set of cavitating venturies for a liquid rocket engine thrust chamber firing test facility have been designed and manufactured. The flow characteristics of the cavitating venturies have been analyzed by experimental and computational methods. Results showed that constant mass flow rate condition was established by the cavitation inside the venturi. However, upstream pressure less than the actual design pressure of the cavitating venturi could not supply a constant flow rate.

• 한국전산유체공학회:학술대회논문집
• /
• 2009.11a
• /
• pp.177-184
• /
• 2009

Cavitating flow simulation is of practical importance for many engineering systems, such as marine propellers, pump impellers, nozzles, injectors, torpedoes, etc. The present work has focused on the simulation of cavitating flow past cylinders with strong side flows. The governing equation is the Navier-Stokes equation based on the homogeneous mixture model. The momentum and energy equation is in the mixture phase while the continuity equation is solved liquid and vapor phase, separately. An implicit dual time and preconditioning method are employed for computational analysis. For the code validation, the results from the present solver have been compared with experiments and other numerical results. A fairly good agreement with the experimental data and other numerical results have been obtained. After the code validation, the strong side flow was applied to include the wake flow effects of the submarine or ocean tide.

• 한국전산유체공학회:학술대회논문집
• /
• 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.