• International Journal of Fluid Machinery and Systems
• /
• v.8 no.4
• /
• pp.240-253
• /
• 2015

The objective of this paper is to evaluate the applicability of mass transfer cavitation models and determine appropriate numerical parameters for cavitating flow simulations. CFD simulations were performed for a NACA66 hydrofoil at cavitation numbers of 1.49 and 1.00, corresponding to steady sheet and unsteady sheet/cloud cavitating regimes using the Kubota and Merkle cavitation models. The Merkle model was implemented into CFX by User Fortran code. The Merkle cavitation model is found to give some improvements for cavitating flow simulation results for these cases. Turbulence modeling is also found to have an important contribution to the prediction quality of the simulations. The relationship between the turbulence viscosity modification, in order to take into account the local compressibility at the vapor/liquid interfaces, and the predicted numerical results is clarified. The limitations of current cavitating flow simulation techniques are discussed throughout the paper.

• Journal of computational fluids engineering
• /
• v.16 no.3
• /
• pp.37-46
• /
• 2011

Unsteady sheet cavitation on a three-dimensional twisted hydrofoil was studied using an unsteady Reynolds-averaged Navier-Stokes equations solver based on a cell-centered finite volume method. As a verification test of the computational method, non-cavitating and cavitating flows over a modified NACA66 foil section were simulated and validated against existing experimental data. The numerical uncertainties of forces and pressure were evaluated for three levels of mesh resolution. The computed pressure on the foil and the cavity shedding behavior were validated by comparing with existing experimental data. The cavity shedding dynamics by re-entrant jets from the end and sides of the cavity were investigated.

• Journal of the Society of Naval Architects of Korea
• /
• v.34 no.1
• /
• pp.11-23
• /
• 1997

This paper contains a new approach to blade section design method for marine propellers. The hydrodynamic characteristics of 2-D section are highly influenced by its geometrical parameters i.e., thickness and camber distributions and leading edge radius etc. To consider fully turbulent flow field near 2-D section. the finite volume method with k-${\varepsilon}$ turbulent model which solve Reynolds time averaged Navier-Stokes(RANS) equation is applied. In this study, O-type grid system that can provide many calculation points on blade surface is used. The results were compared with those of the experiment of NACA0012 to confirm the accuracy of the developed codes. The goal of this study is the development of a blade section with high efficiency and low drag. To achieve this, we carried out the tests of lift, drag and cavitation characteristics in cavitation tunnel. The results of experiment were compared with numerical results in order to validate the proposed blades design method. By comparing the numerical results with the experiments, we found that the new blade section, KH28 allows superior performance in efficiency and cavitation avoidance characteristics. We further investigated the blade section design method and an application study of this section, KH28 to apply to the marine propeller. In order to improve the accuracy of numerical results on prediction of lift and drag, we conclude here that the 2-layer boundary model must be used.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.6 no.2
• /
• pp.406-417
• /
• 2014

This paper presents a numerical optimization method to improve the performance of the propeller with Turbo-Ring using real-coded genetic algorithm. In the presented method, Unimodal Normal Distribution Crossover (UNDX) and Minimal Generation Gap (MGG) model are used as crossover operator and generation-alternation model, respectively. Propeller characteristics are evaluated by a simple surface panel method "SQCM" in the optimization process. Blade sections of the original Turbo-Ring and propeller are replaced by the NACA66 a = 0.8 section. However, original chord, skew, rake and maximum blade thickness distributions in the radial direction are unchanged. Pitch and maximum camber distributions in the radial direction are selected as the design variables. Optimization is conducted to maximize the efficiency of the propeller with Turbo-Ring. The experimental result shows that the efficiency of the optimized propeller with Turbo-Ring is higher than that of the original propeller with Turbo-Ring.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.26 no.2
• /
• pp.141-147
• /
• 2016

Underwater cavitation is one of the most important issues because it causes not only vibration and erosion of submerged bodies but also significant flow noise problems. In this paper, flow noise due to cavitation flows around the NACA66 MOD hydrofoil is numerically investigated. The cavitation flow simulation is conducted using the Reynolds-Averaged Navier-Stokes equations based on finite difference methods. To capture the cavitation phenomena accurately and effectively, the homogeneous mixture model with the Merkle's cavitation model is applied. The predicted results are compared with available experimental data in terms of pressure coefficients and volume fraction, which confirms the validity of numerical results. Based on flow field analysis results, hydro-acoustic noise field due to the cavitation flow is predicted using the Ffowcs-Williams and Hawkings equation derived from the Lighthill's acoustic analogy. The typical lift dipole propagation patterns are identified.

• Journal of computational fluids engineering
• /
• v.17 no.1
• /
• pp.94-100
• /
• 2012

In this paper, the cavitating flows around a hydrofoil have been numerically investigated by using a 2-d multi-phase RANS flow solver based on pseudo-compressibility and a homogeneous mixture model on unstructured meshes. For this purpose, a vertex-centered finite-volume method was utilized in conjunction with 2nd-order Roe's FDS to discretize the inviscid fluxes. The viscous fluxes were computed based on central differencing. The Spalart-Allmaras one equation model was employed for the closure of turbulence. A dual-time stepping method and the Gauss-Seidel iteration were used for unsteady time integration. The phase change rate between the liquid and vapor phases was determined by Merkle's cavitation model based on the difference between local and vapor pressure. Steady state calculations were made for the modified NACA66 hydrofoil at several flow conditions. Good agreements were obtained between the present results and the experiment for the pressure coefficient on a hydrofoil surface. Additional calculation was made for cloud cavitation around the hydrofoil. The observation of the vapor structure, such as cavity size and shape, was made, and the flow characteristics around the cavity were analyzed. Good agreements were obtained between the present results and the experiment for the frequency and the Strouhal number of cavity oscillation.

• Journal of computational fluids engineering
• /
• v.22 no.1
• /
• pp.59-66
• /
• 2017

The guideline of selecting the number of snapshot dataset, $N_s$ in proper orthogonal decomposition(POD) was presented via the analysis of Eigen values based on the singular value decomposition(SVD). In POD, snapshot datasets from the solutions of Euler or Navier-Stokes equations are utilized to SVD and a reduced order model(ROM) is constructed as the combination of Eigen vectors. The ROM is subsequently applied to reconstruct the flowfield data with new set of flow conditions, thereby enhancing the computational efficiency. The overall computational efficiency and accuracy of POD is dependent on the number of snapshot dataset; however, there is no reliable guideline of determining $N_s$. In order to resolve this problem, the order of maximum to minimum Eigen value ratio, O(R) from SVD was analyzed and presented for the decision of $N_s$; in case of steady flow, $N_s$ should be determined to make O(R) be $10^9$. For unsteady flow, $N_s$ should be increased to make O(R) be $10^{11\sim12}$. This strategy of selecting the snapshot dataset was applied to two dimensional NACA0012 airfoil and vortex flow problems including steady and unsteady cases and the numerical accuracies according to $N_s$ and O(R) were discussed.

• Journal of the Society of Naval Architects of Korea
• /
• v.59 no.3
• /
• pp.183-191
• /
• 2022

A design method for a propeller type rim-driven axial-flow turbine for a micro-hydropower system is presented. The turbine consists of pre-stator, impeller and post-stator, where the pre-stator plays a role as a guide vane to provide circumferential velocity to the on-coming flow, and the impeller as a rotational power generator by absorbing angular momentum of the flow. BEM(Blade Element Method), which is based on the turbine Euler equation, is employed to design the pre-stator and impeller blades. NACA 66 thickness form and a=0.8 mean camber line, which is widely accepted as a marine propeller blade section, is used for the pre-stator and turbine blade section. A CFD method, derived from the discretization of the RANS equations, is applied for the analysis of the designed turbine system. The design conditions of the turbine is confirmed by the CFD calculation. Turbine characteristic curve is calculated by the CFD method, in order to provide the performance characteristics at off-design operation conditions. The proposed procedures for the design of a propeller type rim-driven axial-flow turbine are established and confirmed by the CFD analysis.