Computational investigation of cavitating flow around two dimensional NACA 4424 and MHKF240 hydrofoil
Srijna Singh^{1} , Mohammad Danish^{2} , Kaushik Saha^{3}
^{1, 2}Bennett University, Greater Noida, India
^{3}Indian Institute of Technology, New Delhi, India
^{1}Corresponding author
Vibroengineering PROCEDIA, Vol. 29, 2019, p. 159164.
https://doi.org/10.21595/vp.2019.21121
Received 24 October 2019; accepted 31 October 2019; published 28 November 2019
JVE Conferences
This study focuses on the comparison of the performance of two unsymmetrical hydrofoils, NACA 4424 and MHKF240 at 60 angle of attack under cavitation. The Schnerr and Sauer cavitation model along with Realizable $k$$\epsilon $ turbulence model is used for numerical computation in commercial software ANSYS Fluent. The lift, drag and pressure coefficients for different cavitation numbers were studied. Among both the hydrofoils MHKF240 gives a higher lift coefficient which is the parameter of better performance.
Keywords: NACA 4424, MHKF240, realizable kε, Schnerr and Sauer model.
1. Introduction
Cavitation is a phenomenon of formation and subsequent collapse of vapor bubbles that form in the liquid flow when the local pressure falls below the saturated vapor pressure. The inception of cavitation in flow causes severe structural damage due to noise, vibration, and erosion making it complex and interesting area for research [1]. Cavitation instabilities may be of different types such as cavitation surge, rotating cavitation, asymmetric cavitation, and high order instabilities. The cavitation phenomenon can be commonly observed around the hydrofoils and propellers in a marine application.
The hydrofoil is the streamlined body designed in such a way that it can provide sufficient lift required by the ship to move on water. This lift force produced on the hydrofoil is due to the pressure difference between the lower and upper surface of the hydrofoil as the upper side of hydrofoil is having lower pressure. This lower pressure can cause the problem of cavitation. By nature, a cavitating flow is always an unsteady threedimensional flow and to study the cavitation on the hydrofoil various experimental and numerical studies have been reported in the literature. Fabula [2] has applied ThinAerofoil theory and closuresingularity cavity termination method on thin hydrofoil for steady, plane irrotational and incompressible flow to know the phenomenon of cavity termination along with the wake thickness for a different angle of attack. His theoretical results for lift and drag force on the cavitating aerofoil are found to be in good agreement with the experimental results. Similarly, Woods [3] studied the unsteady cavitating flow past a finite length twodimensional flat plate hydrofoil and symmetrical wedge to understand the feature of cavity growth and required sink for the accommodation of displaced liquid.
The occurrence of cavitation can be quantitatively measured by a nondimensional parameter, termed as cavitation number ($\sigma $), which is essentially a measure of how close the pressure in liquid flow is to vapor pressure. Mathematically, the cavitation number is defined as $\sigma =\left({p}_{\infty}{p}_{v}\right)/\left(0.5\rho {U}_{\infty}^{2}\right)$, where ${p}_{\infty}$ is the reference pressure, ${p}_{v}$ is the vapor pressure of the liquid, $\rho $ is the density and ${U}_{\infty}$ is the free stream velocity of liquid flow [4]. Cavitation number is also an important parameter to determine the cavity length on the hydrofoil, as lower the cavitation number more will be the cavity length and below certain value of $\sigma $ the complete hydrofoil is engulfed in the cavity bubble. This phenomenon is known as supercavitation which greatly reduces the skin friction drag on hydrofoil. Green [5] extended thin aerofoil theory on two supercavitating flatplate hydrofoils at an angle ${\gamma}_{1}$ for upper wetted surface and ${\gamma}_{2}$ on the lower wetted surface in the presence of a free surface using twodimensional, incompressible, irrotational, linearized flow model. In addition to cavitation number the cavitation behaviour is also influenced by the Reynold’s number as in many instances the local pressure reduces as the result of accelerating high velocity.
The NACA aerofoil developed by National Advisory Committee for Aeronautics (NACA) are the profiles used for aircraft wings, ship propellers and hydrofoils. The computational investigation of the performance of different NACA profiles under cavitation has been the topic of intense research in the past. To improve the performance and structural stability of hydrofoils, some researchers have modified these NACA profiles. Recently, Shiu et al. [6] has modified different NACA44 hydrofoils by adding finite thickness to the trailing edge and named the hydrofoil as Marine and Hydrokinetic (MHKF1) hydrofoils and then compared the performance of MHKF1 for different angle of attack with NACA 44 series using both experimental and numerical method. They concluded that these modified hydrofoils give higher lift coefficients. They have focused on cavitating phenomenon with respect to change in angle of attack but not mentioned the cavitation behaviour and hydrofoil performance related to different cavitation number.
In this paper, the comparison of the performance of NACA 4424 and MHKF240 hydrofoil for different cavitation numbers ($\sigma )$ has been presented. The angle of attack used is 6° as internal parametric studies gives the stall condition for NACA 44 series at 6°. The MHKF240 hydrofoil is designed by adding finite thickness to the sharp training edge of NACA 4424 hydrofoil. The twophase single fluid homogeneous mixture approach with Schnerr and Sauer cavitation model [7] is used and to study the turbulent flow dynamics the Realizable $k$$\epsilon $ turbulence model is applied [8]. The chord length of both the hydrofoils is 100 mm and span is considered as a unity. The results of NACA4424 and MHKF240 are compared and it was found that the lift coefficient of MHKF240 has higher value as expected.
2. Governing equations
The mass conservation equation and momentum equation are the governing equations for two phase flow [9]:
where mixture density ${(\rho}_{m})$ is related to vapor volume fraction (${\alpha}_{v}$) as:
where $p$ is the pressure, ${\rho}_{v}$ and ${\rho}_{l}$ are the vapor and liquid phase densities. The effective viscosity in Eq. (2) is expressed as ${\mu}_{eff}=\mu +{\mu}_{t}$ where $\mu $ is the molecular viscosity and ${\mu}_{t}$ is the turbulent viscosity which is modelled using different turbulence model. For $k$$\epsilon $ turbulence models:
In this equation $k$ is the turbulent kinetic energy and $\mathrm{\epsilon}$ is dissipation of kinetic energy and ${C}_{\mu}$ is the model constant.
3. Computational method
The commercial CFD simulation software ANSYS Fluent is used to numerically solve twophase cavitating flow on the NACA 4424 and MHKF240 hydrofoil at 6° angle of attack for different cavitation numbers ($\sigma )$. NACA 4424 hydrofoil is having chord length of 100 mm, maximum camber of 4 % at 0.4 chords from the leading edge with maximum thickness of 24 % of the chord. The MHKF240 is having the same chord length but with the finite thickness (2.5455 mm) at the trailing edge [6]. Fig. 1 shows the profile of both the hydrofoils.
To know the hydrodynamic behaviour of the hydrofoils under the cavitation phenomenon the C type domain is made around the hydrofoil having a size of 1800 mm×1500 mm, extending 13.5 times the chord length ahead of leading edge and 12 times the chord length behind the leading edge [10]. The computational domain with the required dimensions and boundary conditions is shown in Fig. 2. The Dirichlet boundary condition is given to the domain with C type boundary in the upstream as velocity inlet with freestream velocity in $x$ direction as 7.03 m/s. The outlet is given as pressure outlet and all other boundaries including hydrofoil is given noslip condition. The Reynold’s number ($R{e}_{c}$) based on the chord length is 0.69×10^{6} [11]. The domain is meshed using unstructured quadrilateral elements and to create the denser mesh around the hydrofoil the sphere of influence is provided as shown is the Fig. 3. To obtain higher computational accuracy the inflation layers are generated around the hydrofoil as shown in Fig. 4.
The pressure based transient mixture model is used for solving the governing equations. The realizable $k$$\epsilon $ turbulence model [12] and Schnerr and Sauer cavitation model is used to calculate the ${\mu}_{t}$ and ${\alpha}_{v}$ respectively. For pressure and velocity coupling Coupled method is used. Here the momentum, volume fraction, turbulent kinetic energy and turbulent dissipation rate are solved using QUICK scheme whereas pressure is computed using PRESTO!
Fig. 1. NACA 4424 and MHKF240 hydrofoil profile
Fig. 2. Computational domain with boundary conditions
Fig. 3. Meshing in computational domain
Fig. 4. Inflation layers around the hydrofoil
4. Results and discussions
The independency of mesh is an important parameter for optimum computational cost without compromising the accuracy of the solution. The grid independence test was performed for both the hydrofoils at $\sigma =\text{1.5}$ and $\alpha =$6°. The timeaveraged lift coefficient and percentage error for all grids concerning the finest grid E are presented in Table 1. and it can be seen that error between grid D and grid E for both the hydrofoils is minimum. Hence, grid D with 96000 elements has been selected as the final mesh for studying the hydrodynamic behaviour of both hydrofoils at different cavitation numbers. The plots for lift coefficient for different grid size are shown in Fig. 5.
Table 1. Time averaged lift coefficient for different grids
S. No

No of grids

Time averaged lift coefficient (${C}_{l}$)
NACA 4424

Time averaged lift coefficient (${C}_{l}$)
MHKF240

Percentage error related to 115000 grids (${C}_{l}$)
NACA 4424

Percentage error related to 115000 grids (${C}_{l}$)
MHKF240

A

29000

0.0791

0.0978

1.4

1.91

B

49000

0.0797

0.0989

2.18

0.80

C

69000

0.0784

0.1006

0.512

0.90

D

96000

0.0781

0.0999

0.128

0.20

E

115000

0.0780

0.0997

–

–

Fig. 5. Time averaged lift coefficient against the grid numbers
The pressure coefficient ${C}_{p}$ as a result of pressure distribution on the hydrofoil is represented in Fig. 6 on the suction side of both NACA 4424 and MHKF240 hydrofoil for $\sigma =$ 1.28 and $\alpha =$ 6°$.$ The straight line on both the plots represents the cavity length on the hydrofoil. The cavity length of MHKF240 foil is more as compared to NACA 4424 hydrofoil as more bubble formation takes place towards the leading edge due to more pressure drop on the suction side.
Fig. 6. Pressure coefficient at the suction side of hydrofoil along the chord length
Fig. 7. Time averaged (3 sec) lift coefficient for different cavitation numbers
The comparison of lift coefficient on NACA 4424 and MHKF240 hydrofoil is plotted in Fig. 7 for the range of cavitation number from 0.17 to 4.12. For calculating the lift coefficient, we normalized our results with the chord length of NACA 4412 in the literature [11]. Here, it can be observed that lift coefficient of MHKF240 is higher as compared to NACA 4424 resulting in higher hydrodynamic performance. The other parameter to check the performance of hydrofoil is the ${\left(l/d\right)}_{max}$ value which is computed based on the wetted area of hydrofoil [13]. MHKF240 have higher ${\left(l/d\right)}_{max}$ giving better performance as compared to NACA 4424 profile. In Fig. 8 and Fig. 9 vapor volume fraction on NACA 4424 and MHKF240 hydrofoil is shown for $\alpha =$6° and $\sigma =$ 1.28. We observed that cavity on MHKF240 starts nearer to the leading edge and had larger cavity length and more vapor fraction as compared to NACA 4424 which results in the lower skin friction drag coefficient on MHKF240 which is described next section.
Fig. 8. Vapor volume fraction on NACA 4424 hydrofoil at $\sigma =$ 1.28
Fig. 9. Vapor volume fraction on MHKF240 hydrofoil
The skin friction drag is the resisting force due to the viscosity of fluid on the object. Fig. 10 represents the skin friction coefficient on the upper side of both the hydrofoils. As the flow separation on the MHKF240 foil near the leading edge is more there is less viscous resistance due to fluid. In Fig. 11 the time evolution of lift coefficients between 1.9 sec to 3$$for both the hydrofoils at$\sigma =\text{1.28}$ is shown which is later used to compute the Strouhal number using Fast Fourier Transformation (FFT) of the lift coefficient.
Fig. 10. Skin friction coefficient at the suction side of hydrofoil
Fig. 11. Time evolution of lift coefficient
To understand the nature of the shedding of the cloud cavity the Strouhal number (St) is used as an important parameter that defines the impact of unsteady oscillations on the flow [14]. To calculate the Strouhal number FFT algorithm is used which converts the signal from a time domain to frequency domain. Fig. 11 shows the power spectral density of the lift coefficient for$\sigma =$ 1.28 giving the first peak value for NACA 4424 at 13.3 Hz and for MHKF240 at 12.3 Hz. Based on the chord length of hydrofoil and peak value of frequency we get St for NACA 4424 as 0.189 and MHKF240 as 0.175.
Fig. 12. Power spectral density of lift coefficient versus frequency
5. Conclusions
In this paper, the cavitation occurrence on NACA 4424 and MHKF240 hydrofoil is been examined for different cavitation numbers at 6° angle of attack. The optimum mesh size on both the hydrofoils was determined and numerical simulation was performed using commercial software ANSYS Fluent to obtain lift coefficient, skin friction drag coefficient, pressure coefficients and Strouhal number. From these results, it is concluded that:
1) Lift coefficient of MHKF240 hydrofoil is more as compared to NACA 4424 giving higher performance.
2) Vapor volume fraction of MHKF240 is higher resulting in lower skin friction drag.
3) Strouhal number for NACA 4424 is 0.189 and for MHKF240 is 0.175 due to lower vapor shedding frequency on MHKF240.
References
 Ganz S. Cavitation: Causes, Effects, Mitigation and Application. Rensselaer Polytechnic Institute, 2012. [CrossRef]
 Fabula A. G. Thinairfoil theory applied to hydrofoils with a single finite cavity and arbitrary freestreamline detachment. Journal of Fluid Mechanics, Vol. 12, Issue 2, 1962, p. 227240. [Publisher]
 Woods L. C. On the theory of growing cavities behind hydrofoils. Journal of Fluid Mechanics, Vol. 19, Issue 1, 1964, p. 124136. [Publisher]
 Franc J. P., Michel J. M. Unsteady attached cavitation on an oscillating hydrofoil. Journal of Fluid Mechanics, Vol. 193, 1988, p. 171189. [Publisher]
 Green T., Street R. L. Two supercavitating hydrofoils near a free surface. Journal of Fluid Mechanics, Vol. 27, Issue 1, 1967, p. 128. [Publisher]
 Shiu H., van Dam C. P., Johnson E., Barone M., Phillips R., Straka W., Fontaine A., Jonson M. A design of a hydorofoil family for currentdriven marinehydorkinetic turbines. 20th International Conference on Nuclear Engineering and the ASME Power Conference, 2012. [CrossRef]
 Niedzwiedzka A. Homogeneous cavitation modelinganalysis of basics of homogeneous cavitation modeling – analysis of basics of mathematical. Modelowanie Inżynierskie, Vol. 65, 2017, p. 7582. [CrossRef]
 Huang S., He M., Wang C., Chang X. Simulation of cavitating flow around a 2D hydrofoil. Journal of Marine Science and Application, Vol. 9, Issue 1, 2010, p. 6368. [Publisher]
 Saha K. Modelling of cavitation in nozzles for diesel injection applications. UWSpace, 2014, http://hdl.handle.net/10012/8628. [CrossRef]
 Kumar T. R. S., Vidyapeetham A. V., Ramakrishnananda B., Vidyapeetham A. V., Arulkrishnan A. Numerical investigation of two element camber morphing airfoil in low Reynolds number flows. Journal of Engineering Science and Technology, Vol. 12, Issue 7, 2017, p. 19391955. [CrossRef]
 Kermeen R. W. Water tunnel tests of NACA 4412 and Walchner profile 7 hydrofoils in noncavitating and cavitating flows. Rep. No 475, 1956. [CrossRef]
 Shih T. H., Liou W. W., Shabbir A., Yang Z., Zhu J. A new KΕ{lunate} eddy viscosity model for high Reynolds number turbulent flows. Computers and Fluids, Vol. 24, Issue 3, 1995, p. 227238. [Publisher]
 Loftin Laurence K. Quest for Performance: The Evolution of Modern Aircraft. National Aeronautics and Space Administration, 2012. [CrossRef]
 Dular M., Bachert R. The issue of Strouhal number definition in cavitating flow. Strojniski Vestnik/Journal of Mechanical Engineering, Vol. 55, Issue 11, 2009, p. 666674. [CrossRef]