Reduction of Hydrodynamic Noise of 3D Hydrofoil with Spanwise Microgrooved Surfaces Inspired by Sharkskin

: Loud hydrodynamic noise is not only potentially harmful to the health of organisms in the ocean, but it is also a threat to the survival of underwater vehicles. Di ﬀ erent from the general noise reduction technologies at present, a new idea for a ﬂow-induced noise reduction design with spanwise microgrooved surfaces inspired by sharkskin is introduced in this paper. Large eddy simulations (LES) combined with the Ffowcs Williams and Hawkings (FW-H) equation are adopted to simulate the hydrodynamic noise of the three-dimensional (3D) hydrofoil. The accuracy of the numerical predictions is checked against existing experimental data, achieving good agreement. With the increase of observing distance, the noise reduction e ﬀ ect at the trailing edge direction is gradually apparent, and a maximum noise reduction of up to 7.28 dB can be observed. It is seen from the noise spectra of the biomimetic hydrofoil that the main peaks are eliminated, and the noise level at high frequency is also decreased. The cause of noise reduction lies in the secondary vortex generated in the microgrooves, which hinder the process of turbulence, consume the energy of the ﬂow, and weaken the intensity of turbulent burst. The results of this study provide a new way to design low-noise underwater structures with hydrofoils.


Introduction
Hydrofoil or airfoil is an important structure applied in many underwater devices, ranging from the sail hulls of submarines and the flapping wings of underwater vehicles to ocean turbines and marine propellers. However, noise originates from various sources on these devices, such as (1) the self-noise of airfoil generated by the interaction of boundary layers and the wake and (2) the noise caused by cavitation of some ocean engineering applications, for example, marine propellers. As a kind of pollution, noise is harmful to both the health of organisms in nature and the normal operation of the marine structures. The noise problems of ocean engineering structures with airfoils is currently a great concern [1][2][3][4][5].
In 1989, Brooks et al.
[6] identified five mechanisms for airfoil self-noise due to specific boundary-layer phenomena, including the noise caused by boundary-layer separation and large-scale separation. In the following several decades, researchers have carried out a lot of studies to reduce the noise of airfoils or blades. Inspired by the wings of low-noise flying owls, Howe [7,8] theoretically analyzed the effects of serration structure in the trailing-edge of airfoil on aerodynamic noise. Since then, the noise reduction mechanism and the key parameters influencing the noise reduction effect of serrated trailing edge have been studied by a lot of scholars [9][10][11][12][13]. Some scholars have even modified the serrated trailing edge, using non-flat plate serrations [14,15], serrations with flap angles [16], iron-shaped serrations [17], and combed-sawtooth serrations [18]. Moreover, some researchers have studied the noise reduction effects of wavy leading edge inspired by biomimetic designs [19][20][21][22][23][24]. In addition, some other methods have been proposed to eliminate aerodynamic noise, including performing trailing edge perforation [25], using different blade shapes [26], introducing porous trailing edge inserts [27], and setting distributed boundary-layer suctions [28]. The above methods are all effective in reducing noise to a certain extent; however, most of these studies are still in the stage of academic research and much work has to be done to narrow the gap between the academic research and engineering application. Thus, it is necessary to introduce fresh ideas to achieve better noise reduction.
A biomimetic surface is a new technology that imitates the surface morphology of an organism to improve the performance of equipment. Research into biomimetic surface technology started in the 1970s. At first, researchers paid more attention to the drag reduction of the flow field control in the structural boundary layer [29,30]. Choi [31] pointed out that biomimetic surface technology is a smart flow control method, because the microgrooves not only control momentum, but also control flow noise due to the reduction of pressure fluctuations. The experimental results showed that the energy level at low frequency is dramatically decreased [32]. Joslin et al. [33] compared studies concerning flow and noise control. They concluded that microgrooves can lead to drag reduction, and the reduced energy in turbulent boundary layer might lead to reduced interior noise. In reality, the first study on noise reduction based on microgrooves was conducted by Gillcrist and Reidy [34]. They carried out experiments on drag reduction and noise reduction of a marine vehicle with microgrooved surface coatings. The results showed that an ideal noise reduction effect could be obtained. Then, Shi et al. [35] utilized the water tunnel at Northwestern Polytechnical University to investigate the flow-induced noise of underwater vehicles with microgrooved surfaces. When the width of the microgroove was 0.1 mm and the height of the microgrooves was 0.09 mm, there was an obvious noise reduction of up to 5 dB. Therefore, the bionic microstructural microgrooved surface can control the flow field in the boundary layer. At the same time, it can suppress the flow noise generated by the dipole sound source produced by the solid boundary reaction to the fluid to a certain extent. Figure 1 collects the present technologies to reduce the aerodynamic or hydrodynamic noise of airfoils or hydrofoils, and it also provides the most promising technology to reduce noise by biomimetic surface design, although, in this study, the greatest concern is drag reduction design.
As emphasized by Fu et al. [36], sharkskin is not a smooth surface but a kind of microgrooved surface with micro-scales called dermal denticles. Currently, research on drag reduction with microgrooved surfaces inspired by sharkskin is a hot topic [36][37][38][39]. However, with respect to airfoils or hydrofoils, most scholars have concentrated on flow field characteristics [40][41][42][43][44], and to the authors' knowledge, few studies have reported on the noise performance of airfoils with microgrooved surfaces inspired by sharkskin, which is addressed in this paper. In fact, the microgrooves of real sharkskin are not parallel to the flow direction totally. Some researchers have also investigated microgrooves or riblets across the flow direction [36,39,43,44], which is the research object in this study. The layout of this paper is as follows: Section 1 gives a brief introduction and Section 2 briefly demonstrates the theory of numerical simulation of this paper. In Section 3, the implementation of microgrooves on the surface of hydrofoil is carried out, and then, the computational model is illustrated in Section 4. The accuracy of numerical model predictions for both the flow field and sound field is validated against the experimental results in Section 5. In Section 6, results and discussions of hydrodynamic performance and sound characteristics in near-field and in far-field are demonstrated. Finally, conclusions are presented in the last section.
width of the microgroove was 0.1 mm and the height of the microgrooves was 0.09 mm, there was an obvious noise reduction of up to 5 dB. Therefore, the bionic microstructural microgrooved surface can control the flow field in the boundary layer. At the same time, it can suppress the flow noise generated by the dipole sound source produced by the solid boundary reaction to the fluid to a certain extent. Figure 1 collects the present technologies to reduce the aerodynamic or hydrodynamic noise of airfoils or hydrofoils, and it also provides the most promising technology to reduce noise by biomimetic surface design, although, in this study, the greatest concern is drag reduction design.

Computational Methods
Large eddy simulations (LES) combined with the Ffowcs Williams and Hawkings (FW-H) equation were applied to simulate the hydrodynamic noise of 3D hydrofoil. LES were performed first to calculate the hydrodynamic characteristics of the flow field with a transient Computational Fluid Dynamics (CFD) method [45]. Then, the sound pressure levels were obtained from the fluctuating surface pressure on the hydrofoils based on the FW-H equation [46]. The flowchart of the numerical simulation is given in Figure 2. As emphasized by Fu et al. [36], sharkskin is not a smooth surface but a kind of microgrooved surface with micro-scales called dermal denticles. Currently, research on drag reduction with microgrooved surfaces inspired by sharkskin is a hot topic [36][37][38][39]. However, with respect to airfoils or hydrofoils, most scholars have concentrated on flow field characteristics [40][41][42][43][44], and to the authors' knowledge, few studies have reported on the noise performance of airfoils with microgrooved surfaces inspired by sharkskin, which is addressed in this paper. In fact, the microgrooves of real sharkskin are not parallel to the flow direction totally. Some researchers have also investigated microgrooves or riblets across the flow direction [36,39,43,44], which is the research object in this study. The layout of this paper is as follows: Section 1 gives a brief introduction and Section 2 briefly demonstrates the theory of numerical simulation of this paper. In Section 3, the implementation of microgrooves on the surface of hydrofoil is carried out, and then, the computational model is illustrated in Section 4. The accuracy of numerical model predictions for both the flow field and sound field is validated against the experimental results in Section 5. In Section 6, results and discussions of hydrodynamic performance and sound characteristics in near-field and in far-field are demonstrated. Finally, conclusions are presented in the last section.

Computational Methods
Large eddy simulations (LES) combined with the Ffowcs Williams and Hawkings (FW-H) equation were applied to simulate the hydrodynamic noise of 3D hydrofoil. LES were performed first to calculate the hydrodynamic characteristics of the flow field with a transient Computational Fluid Dynamics (CFD) method [45]. Then, the sound pressure levels were obtained from the fluctuating surface pressure on the hydrofoils based on the FW-H equation [46]. The flowchart of the numerical simulation is given in Figure 2.

Hydrodynamic Calculation
Sound Source Analogy

Acoustic Calculation
Acoustic Results

Large Eddy Simulation Model
Turbulence is the three-dimensional unsteady random motion observed in fluids at moderate to high Reynolds numbers [47]. There are three main numerical methods, including direct numerical simulation (DNS), Reynolds-averaged Navier-Stokes (RANS) equation, and LES. The most accurate results could be obtained with DNS. However, in most situations, it is not feasible to resolve the wide range of scales in time and space by using DNS, as the central processing unit (CPU) requirements

Large Eddy Simulation Model
Turbulence is the three-dimensional unsteady random motion observed in fluids at moderate to high Reynolds numbers [47]. There are three main numerical methods, including direct numerical simulation (DNS), Reynolds-averaged Navier-Stokes (RANS) equation, and LES. The most accurate results could be obtained with DNS. However, in most situations, it is not feasible to resolve the wide range of scales in time and space by using DNS, as the central processing unit (CPU) requirements would by far exceed the available computing power. For this reason, averaging procedures are applied to the Navier-Stokes equations to obtain the flow field characteristics equivalently and rapidly. The most well-known one is the RANS equation, through which a smooth variation of the averaged velocity and pressure fields can be obtained. With low spatial resolution and a small workload, it cannot reflect the fluctuation of the flow field precisely. As a result, LES is proposed to overcome the shortcomings of the other two above methods. It not only greatly reduces the computational complexity, but also obtains the low-frequency evolution information in the instantaneous turbulent flow field with high numerical accuracy. Therefore, this method can simulate the detailed characteristics of the pulsation in the flow field better, so as to calculate the acoustic field more accurately. An incompressible simulation was used for the LES model in this paper.
In the theory of LES, there are two parts for turbulent vortices: the low-pass filtered quantity ϕ, representing large scale vortices before the filter cut-off, and the sub-grid scale quantity ϕ , representing small-scale vortices after the filter cut-off. The variable ϕ after filtering is the average in the spatial domain, which is obtained directly by solving the governing equation of LES, while the filtered variable ϕ(x) is solved by establishing the turbulence model.
where D is the flow area, x is the spatial coordinate of the actual flow, x is the large-scale spatial coordinate after filtering, G(x, x ) is the filtering function, which can be written as: where V is the volume of a computational cell. The substitution of the filtering function G(x, x ) of Equation (2) into Equation (1) results in: Then, the governing equations of LES can be obtained by substituting the mass conservation equation and the transient Navier-Stokes equation into Equation (3): where τ ij is the term of sub-grid stress, also called the filtered stress tensor, represented as: The sub-grid stress τ ij denotes the momentum transport between small-scale vortex ϕ and large-scale vortex ϕ. Additionally, the effect of small-scale fluctuations on the overall flow field is quantified. To make Equation (5) enclosed, the corresponding sub-grid model has to be constructed.
At present, the sub-grid model is based largely on the hypothesis proposed by Boussinesq [48], and the sub-grid stress is determined as: where δ ij is the sub-grid scale of Reynolds stress, µ t is the coefficient of the sub-grid eddy viscosity, and S ij is the deformation rate tensor corresponding to the scale of solution, which can be written as: The eddy viscosity hypothesis above will transform solving the sub-grid stress into solving for the sub-grid scale eddy viscosity coefficient. The Smagorinsky-Lilly model is developed to solve this problem. The coefficient of the sub-grid eddy viscosity, µ t , can be obtained as: where L s is the mixing length of sub-grid, and it can be expressed as: where k is Karman constant, d is the distance from the node to the nearest wall, C s is Smagorinsky constant, and ∆ is the filter scale.

Ffowcs-Williams and Hawkings Acoustics Model
In the direct method of predicting hydrodynamically generated noise, both generation and propagation of sound waves are directly computed by solving the appropriate fluid dynamics equations. It is thus computationally expensive and difficult inasmuch as it requires very fine computational meshes, highly accurate numbers all the way to receivers, and acoustically nonreflecting boundary conditions [49], whereas for the methods based on acoustic analogy, such as the FW-H equation, the propagation of sound from its generation is essentially decoupled, allowing one to separate the flow solution process form the acoustic analysis. By solving the FW-H equation, the sound pressure level spectrum and total sound pressure level at the observer points are then calculated from the fluctuating surface pressure on the hydrofoil.
The Lighthill equation is the foundation of hydrodynamic noise, which is derived by the RANS equation and the continuity equation, expressed as: where Q is the radiation of unit volume source per unit time and F is the generalized force per unit volume. Subtracting the result of differentiating Equation (12) with respect to time t from the result of differentiating Equation (13) with respect to x i , results in: where T ij is the Lighthill stress tensor, represented as: It can be seen from Equation (14) that, if radiation Q changes with time, generalized force F changes with time, or Lighthill stress tensor T ij changes, hydrodynamic noise will be produced.
For solid wall turbulence, the general solution for Equation (14) can be written as: where l j is cosine in the normal direction of source plane area S and the direction is from inside to outside for the boundary-layer flow of a flat plate, l 1 = l 3 = 0, l 2 = 1. In this way, P 1 = −σ 12 , P 2 = p, and P 3 = −σ 32 , where σ 12 and σ 32 represent the shear stress p 12 and p 32 , separately. Substituting the value of P i into Equation (16) yields the following sound pressure equation: Then, a further transformation is conducted into a form suitable for engineering calculation of hydrodynamic sound pressure: From Equation (18), it can be seen that the hydrodynamic noise is composed of monopole, dipole, and quadrupole sound sources. Then, the sound pressure in the time domain is obtained. After Fourier transformation, the noise characteristics in the frequency domain are calculated. The sound pressure level (SPL) is defined as SPL(dB) = 20 log 10 p p re f (20) where p is the acoustic pressure and p re f is the reference acoustic pressure. For sound transmission in water, the reference pressure is defined as 1 µPa(1 × 10 −6 Pa). Correspondingly, the overall sound pressure level (OASPL) is obtained as OASPL(dB) = 10 log 10 where SPL(i) is the octave band sound pressure level in the broadband frequency domain and n is the number of octave bands from the lowest frequency to the highest frequency concerned.

Design of Spanwise Microgrooves on the Surface of Hydrofoils
The airfoil shapes used in this paper are created from the FF-77-W airfoil shape, which is available at [50]. As is shown in Figure 3, the relative thickness of the airfoil is 14.8% of the length of the chord. When the fluid flows around the hydrofoil, it will produce drag in the flow direction and lift in the vertical direction. Boundary layer separation occurs when the flow angle is larger than a certain value. The location of the separation point and the flow pattern in the boundary layer are closely related to the surface structure of the airfoil, and they also determine the resistance of the airfoil. The turbulent boundary layer has stronger fluid mixing, more uniform fluid kinetic energy distribution, and smaller shape resistance compared with that of laminar boundary layer. Therefore, increasing the surface roughness of the structure is often an effective measure to improve the flow field.
According to the results of Walsh [51][52][53], the optimal drag reduction effect could be obtained with microgrooves when its non-dimensional height h + and non-dimensional spacing s + are both 15. The expression of h + and s + , according to Walsh, are defined as where h is the microgroove height, s is the microgroove spacing, v is the kinematic viscosity, u ∞ is the free-stream velocity, and c f is the local skin friction coefficient. In this way, the dimension parameters of the microgrooves are determined. As a preliminary design inspired by [43] and [44], partial coverage is applied to the upper surface from 40% of the airfoil chord to 80% of the airfoil chord. The overall 3D hydrofoil partially covered by microgrooves and details of the microgrooves are shown in Figure 4. increasing the surface roughness of the structure is often an effective measure to improve the flow field.
According to the results of Walsh [51][52][53], the optimal drag reduction effect could be obtained with microgrooves when its non-dimensional height + h and non-dimensional spacing + s are both 15. The expression of + h and + s , according to Walsh, are defined as where h is the microgroove height, s is the microgroove spacing, v is the kinematic viscosity,  u is the free-stream velocity, and f c is the local skin friction coefficient. In this way, the dimension parameters of the microgrooves are determined. As a preliminary design inspired by [43] and [44], partial coverage is applied to the upper surface from 40% of the airfoil chord to 80% of the airfoil chord. The overall 3D hydrofoil partially covered by microgrooves and details of the microgrooves are shown in Figure 4.

Computational Model
The transient simulations were performed utilizing the FLUENT code in the ANSYS v15.0 environment. In the simulation, the computational domain was sized to decrease the blockage effect and to allow for full development of the upstream flow. The hydrofoil, with a fixed chord c , stood in the middle position of the vertical direction of the computational domain. The upstream inlet was c 5 from the leading edge of the hydrofoil, and the outlet boundary was placed c 10 downstream of the trailing edge of the hydrofoil. Spanwise length in the z -direction was set to 50 mm (0.99c). The computational domain and the corresponding boundary conditions are shown in Figure 5. increasing the surface roughness of the structure is often an effective measure to improve the flow field.
According to the results of Walsh [51][52][53], the optimal drag reduction effect could be obtained with microgrooves when its non-dimensional height + h and non-dimensional spacing + s are both 15. The expression of + h and + s , according to Walsh, are defined as where h is the microgroove height, s is the microgroove spacing, v is the kinematic viscosity,  u is the free-stream velocity, and f c is the local skin friction coefficient. In this way, the dimension parameters of the microgrooves are determined. As a preliminary design inspired by [43] and [44], partial coverage is applied to the upper surface from 40% of the airfoil chord to 80% of the airfoil chord. The overall 3D hydrofoil partially covered by microgrooves and details of the microgrooves are shown in Figure 4.

Computational Model
The transient simulations were performed utilizing the FLUENT code in the ANSYS v15.0 environment. In the simulation, the computational domain was sized to decrease the blockage effect and to allow for full development of the upstream flow. The hydrofoil, with a fixed chord c , stood in the middle position of the vertical direction of the computational domain. The upstream inlet was c 5 from the leading edge of the hydrofoil, and the outlet boundary was placed c 10 downstream of the trailing edge of the hydrofoil. Spanwise length in the z -direction was set to 50 mm (0.99c). The computational domain and the corresponding boundary conditions are shown in Figure 5.

Computational Model
The transient simulations were performed utilizing the FLUENT code in the ANSYS v15.0 environment. In the simulation, the computational domain was sized to decrease the blockage effect and to allow for full development of the upstream flow. The hydrofoil, with a fixed chord c, stood in the middle position of the vertical direction of the computational domain. The upstream inlet was 5c from the leading edge of the hydrofoil, and the outlet boundary was placed 10c downstream of the trailing edge of the hydrofoil. Spanwise length in the z-direction was set to 50 mm (0.99c). The computational domain and the corresponding boundary conditions are shown in Figure 5. Due to the small dimensions of microgrooves (less than 1 mm), unstructured grids were applied in the computational domain. In order to make the calculation results more accurate, (1) mesh refinement was carried out at the microgrooves, and (2) boundary layer mesh was applied at the surface of the hydrofoil with more layers in the boundary layers and fine element size of the boundary layers. The details of the grids near the leading edge region, the trailing edge region, and the microgrooves are shown in Figure 6.
A flow with uniform and steady velocity of 1 m/s was applied at the inlet of the computational domain, and a pressure outlet boundary was imposed at the outlet of the computational domain. A no-slip boundary condition was applied at the surface of the hydrofoils. To improve the reliability of the numerical model, symmetry boundary conditions were used at the four side walls of the domain, which allowed us to consider the computational domain in a larger domain, avoiding the effects of walls.

Numerical Method Verification and Validation
Because there is no available experimental data about the hydrofoils studied in this paper, airfoil NACA0012 was selected to test the reliability and accuracy of the present numerical model. The verification and validation of the present numerical method were as follows. The Reynolds numbers of the NACA0012 airfoil were 1.0 × 10 5 for grid independence verification and acoustical validation and 2.88 × 10 6 for the flow field validation. The transient pressure-based solver was adopted with the drop of all scaled residuals below 1e-5, which was employed as the convergence criterion. The time step size was 2e-5, and the maximum number of iterations was set as 5e4, which enabled the sound frequency range calculated to be 1 Hz~20 kHz with an interval of about 1 Hz. The non-iterative timeadvancement scheme, which significantly speeds up transient simulations by performing only a single outer iteration per time step, was applied in the simulation. The fractional-step method was adopted in ANSYS FLUENT as a velocity-coupling scheme in a non-iterative time-advancement algorithm for the solution method in this paper.

Verification
A grid resolution study was performed to evaluate the influence of grid resolution on the drag coefficient of five different meshing schemes. The angle of attack of the airfoil was zero. The details of the five mesh schemes, including the total number of elements, the parameters of boundary layers, Due to the small dimensions of microgrooves (less than 1 mm), unstructured grids were applied in the computational domain. In order to make the calculation results more accurate, (1) mesh refinement was carried out at the microgrooves, and (2) boundary layer mesh was applied at the surface of the hydrofoil with more layers in the boundary layers and fine element size of the boundary layers. The details of the grids near the leading edge region, the trailing edge region, and the microgrooves are shown in Figure 6.
A flow with uniform and steady velocity of 1 m/s was applied at the inlet of the computational domain, and a pressure outlet boundary was imposed at the outlet of the computational domain. A no-slip boundary condition was applied at the surface of the hydrofoils. To improve the reliability of the numerical model, symmetry boundary conditions were used at the four side walls of the domain, which allowed us to consider the computational domain in a larger domain, avoiding the effects of walls. Due to the small dimensions of microgrooves (less than 1 mm), unstructured grids were applied in the computational domain. In order to make the calculation results more accurate, (1) mesh refinement was carried out at the microgrooves, and (2) boundary layer mesh was applied at the surface of the hydrofoil with more layers in the boundary layers and fine element size of the boundary layers. The details of the grids near the leading edge region, the trailing edge region, and the microgrooves are shown in Figure 6.
A flow with uniform and steady velocity of 1 m/s was applied at the inlet of the computational domain, and a pressure outlet boundary was imposed at the outlet of the computational domain. A no-slip boundary condition was applied at the surface of the hydrofoils. To improve the reliability of the numerical model, symmetry boundary conditions were used at the four side walls of the domain, which allowed us to consider the computational domain in a larger domain, avoiding the effects of walls.

Numerical Method Verification and Validation
Because there is no available experimental data about the hydrofoils studied in this paper, airfoil NACA0012 was selected to test the reliability and accuracy of the present numerical model. The verification and validation of the present numerical method were as follows. The Reynolds numbers of the NACA0012 airfoil were 1.0 × 10 5 for grid independence verification and acoustical validation and 2.88 × 10 6 for the flow field validation. The transient pressure-based solver was adopted with the drop of all scaled residuals below 1e-5, which was employed as the convergence criterion. The time step size was 2e-5, and the maximum number of iterations was set as 5e4, which enabled the sound frequency range calculated to be 1 Hz~20 kHz with an interval of about 1 Hz. The non-iterative timeadvancement scheme, which significantly speeds up transient simulations by performing only a single outer iteration per time step, was applied in the simulation. The fractional-step method was adopted in ANSYS FLUENT as a velocity-coupling scheme in a non-iterative time-advancement algorithm for the solution method in this paper.

Verification
A grid resolution study was performed to evaluate the influence of grid resolution on the drag coefficient of five different meshing schemes. The angle of attack of the airfoil was zero. The details of the five mesh schemes, including the total number of elements, the parameters of boundary layers,

Numerical Method Verification and Validation
Because there is no available experimental data about the hydrofoils studied in this paper, airfoil NACA0012 was selected to test the reliability and accuracy of the present numerical model. The verification and validation of the present numerical method were as follows. The Reynolds numbers of the NACA0012 airfoil were 1.0 × 10 5 for grid independence verification and acoustical validation and 2.88 × 10 6 for the flow field validation. The transient pressure-based solver was adopted with the drop of all scaled residuals below 1.0 × 10 −5 , which was employed as the convergence criterion. The time step size was 2.0 × 10 −5 , and the maximum number of iterations was set as 5.0 × 10 4 , which enabled the sound frequency range calculated to be 1 Hz~20 kHz with an interval of about 1 Hz. The non-iterative time-advancement scheme, which significantly speeds up transient simulations by performing only a single outer iteration per time step, was applied in the simulation. The fractional-step method was adopted in ANSYS FLUENT as a velocity-coupling scheme in a non-iterative time-advancement algorithm for the solution method in this paper.

Verification
A grid resolution study was performed to evaluate the influence of grid resolution on the drag coefficient of five different meshing schemes. The angle of attack of the airfoil was zero. The details of the five mesh schemes, including the total number of elements, the parameters of boundary layers, and the element size along the spanwise direction are listed in Table 1. The drag coefficient at the angle of attack of zero is calculated by where F d is the drag force and l is the spanwise length of the hydrofoil.
where d F is the drag force and l is the spanwise length of the hydrofoil.

Validation
In order to validate the accuracy of the numerical method used in this paper, the flow field characteristics and sound field performance were evaluated with the experimental results simultaneously.
First, the surface pressure coefficients' distribution on the NACA0012 airfoil with an angle of attack at 0° was compared with the experiment results [54], as shown in Figure 8. In the simulation, the Reynolds number was 2.88 × 10 6 , which was in accordance with the experimental test. It was shown that the numerical results have fairly good agreement with the experimental results tested in the wind tunnel.

Validation
In order to validate the accuracy of the numerical method used in this paper, the flow field characteristics and sound field performance were evaluated with the experimental results simultaneously.
First, the surface pressure coefficients' distribution on the NACA0012 airfoil with an angle of attack at 0 • was compared with the experiment results [54], as shown in Figure 8. In the simulation, the Reynolds number was 2.88 × 10 6 , which was in accordance with the experimental test. It was shown that the numerical results have fairly good agreement with the experimental results tested in the wind tunnel. Then, the predictive SPL by experimental test for the NACA0012 airfoil at the observing point was used to verify the accuracy of the acoustic signals. The observing point was perpendicular to and 1.22 m from the trailing edge and the model midspan. In the validation, the model of smooth surface airfoil where the boundary layers were untripped was chosen, the inlet flow velocity was 31.7m/s, and the angle of attack was 0°. Figure 9 shows the comparison of SPL in 1/3 octave frequency bands of 30.48-cm-chord NACA0012 airfoil in the experimental test [6] and the numerical simulation. It can be observed that the numerical results are in close agreement with the experimental results at low and medium frequencies, while for the high-frequency range, the trend was consistent. Considering that the most concern was concentrated on the main peaks of the noise spectrum at lower frequency in the subsequent analysis, which dominated the main noise energy, a conclusion could be drawn that the numerical method applied in this paper is able to study the noise reduction of biomimetic hydrofoil.

Results and Discussions
The numerical calculations of the hydrofoil were conducted at an angle of attack of 0° and a Reynolds number of 5.0 × 10 4 , based on using the hydrofoil chord as the characteristic length. The hydrofoil was created from the FF-77-W airfoil shape, whose relative thickness was 14.8% of the length of the chord, as mentioned in Section 3. The dimensions of the computational domain are shown in Section 4. A similar density grid (Mesh 4), as shown in Section 5.1, was adopted to calculate both the flow field characteristics and sound field performance of the original hydrofoil model and the biomimetic hydrofoil model. To complete the numerical calculation, the relevant material properties of water used in this paper are listed in Table 2. Then, the predictive SPL by experimental test for the NACA0012 airfoil at the observing point was used to verify the accuracy of the acoustic signals. The observing point was perpendicular to and 1.22 m from the trailing edge and the model midspan. In the validation, the model of smooth surface airfoil where the boundary layers were untripped was chosen, the inlet flow velocity was 31.7m/s, and the angle of attack was 0 • . Figure 9 shows the comparison of SPL in 1/3 octave frequency bands of 30.48-cm-chord NACA0012 airfoil in the experimental test [6] and the numerical simulation. It can be observed that the numerical results are in close agreement with the experimental results at low and medium frequencies, while for the high-frequency range, the trend was consistent. Considering that the most concern was concentrated on the main peaks of the noise spectrum at lower frequency in the subsequent analysis, which dominated the main noise energy, a conclusion could be drawn that the numerical method applied in this paper is able to study the noise reduction of biomimetic hydrofoil. Then, the predictive SPL by experimental test for the NACA0012 airfoil at the observing point was used to verify the accuracy of the acoustic signals. The observing point was perpendicular to and 1.22 m from the trailing edge and the model midspan. In the validation, the model of smooth surface airfoil where the boundary layers were untripped was chosen, the inlet flow velocity was 31.7m/s, and the angle of attack was 0°. Figure 9 shows the comparison of SPL in 1/3 octave frequency bands of 30.48-cm-chord NACA0012 airfoil in the experimental test [6] and the numerical simulation. It can be observed that the numerical results are in close agreement with the experimental results at low and medium frequencies, while for the high-frequency range, the trend was consistent. Considering that the most concern was concentrated on the main peaks of the noise spectrum at lower frequency in the subsequent analysis, which dominated the main noise energy, a conclusion could be drawn that the numerical method applied in this paper is able to study the noise reduction of biomimetic hydrofoil.

Results and Discussions
The numerical calculations of the hydrofoil were conducted at an angle of attack of 0° and a Reynolds number of 5.0 × 10 4 , based on using the hydrofoil chord as the characteristic length. The hydrofoil was created from the FF-77-W airfoil shape, whose relative thickness was 14.8% of the length of the chord, as mentioned in Section 3. The dimensions of the computational domain are shown in Section 4. A similar density grid (Mesh 4), as shown in Section 5.1, was adopted to calculate both the flow field characteristics and sound field performance of the original hydrofoil model and the biomimetic hydrofoil model. To complete the numerical calculation, the relevant material properties of water used in this paper are listed in Table 2.

Results and Discussions
The numerical calculations of the hydrofoil were conducted at an angle of attack of 0 • and a Reynolds number of 5.0 × 10 4 , based on using the hydrofoil chord as the characteristic length. The hydrofoil was created from the FF-77-W airfoil shape, whose relative thickness was 14.8% of the length of the chord, as mentioned in Section 3. The dimensions of the computational domain are shown in Section 4. A similar density grid (Mesh 4), as shown in Section 5.1, was adopted to calculate both the flow field characteristics and sound field performance of the original hydrofoil model and the biomimetic hydrofoil model. To complete the numerical calculation, the relevant material properties of water used in this paper are listed in Table 2. Table 2. Material properties of water used in the simulation.

Description Parameter and Value
Temperature T = 20 • C Density ρ = 998.2 kg/m 3 Viscosity µ = 0.001003 kg/(m·s) Sound speed v = 1483 m/s Reference acoustic pressure p re f = 1.0 × 10 −6 Pa 6.1. Analysis of Hydrodynamic Performance Figure 10 shows time histories of the hydrodynamic force coefficients of the original hydrofoil and biomimetic hydrofoil at 0.8-2.0 s, with a running average over 5 times of calculation. It can be seen that the curves of the lift coefficient and drag coefficient for the original hydrofoil fluctuated periodically. Further, it can be easily seen that there was a small pulse near the peak of each period of drag coefficient, which was largely due to the instability of the flow velocity gradient, while the for biomimetic model, there was not an evident periodicity of properties for the lift coefficient curve, but the values of the lift coefficient were very close. The drag coefficient curve for the biomimetic model was even much more arbitrary than the lift coefficient curve. Moreover, it is shown in Figure 10 that the lift coefficient of the biomimetic model was almost unchanged, and the drag coefficient of the biomimetic model was slightly decreased. In total, the amplitudes of fluctuations with variation of time were significantly lower than that of the original hydrofoil.
In addition, the sound pressure with variation of time of the original hydrofoil at three different observing points was depicted in Figure 11. The three points, named Receiver 1, Receiver 2, and Receiver 3, were perpendicular to the trailing edge and 0.1 c, 0.3 c, and 5.5 c from the model midspan, respectively. From Figure 11, it can be seen that the amplitudes of the sound waves gradually decreased with the sound transmission with distance. Further, one can observe that the period of the three curves was the same.   Figure 10 shows time histories of the hydrodynamic force coefficients of the original hydrofoil and biomimetic hydrofoil at 0.8-2.0 s, with a running average over 5 times of calculation. It can be seen that the curves of the lift coefficient and drag coefficient for the original hydrofoil fluctuated periodically. Further, it can be easily seen that there was a small pulse near the peak of each period of drag coefficient, which was largely due to the instability of the flow velocity gradient, while the for biomimetic model, there was not an evident periodicity of properties for the lift coefficient curve, but the values of the lift coefficient were very close. The drag coefficient curve for the biomimetic model was even much more arbitrary than the lift coefficient curve. Moreover, it is shown in Figure  10 that the lift coefficient of the biomimetic model was almost unchanged, and the drag coefficient of the biomimetic model was slightly decreased. In total, the amplitudes of fluctuations with variation of time were significantly lower than that of the original hydrofoil.
In addition, the sound pressure with variation of time of the original hydrofoil at three different observing points was depicted in Figure 11. The three points, named Receiver 1, Receiver 2, and Receiver 3, were perpendicular to the trailing edge and 0.1 c, 0.3 c, and 5.5 c from the model midspan, respectively. From Figure 11, it can be seen that the amplitudes of the sound waves gradually decreased with the sound transmission with distance. Further, one can observe that the period of the three curves was the same.

Noise Characteristics along the Circumferential Direction
A total of 16 observing points were defined to evaluate the acoustic performance of near-field hydrofoil noise along the circumferential direction, which are presented in Figure 12. It can be seen that the observing points were evenly distributed in a circle with a radius of 30.00 mm. As is shown in Figure 12, the origin of the coordinate was at the midspan of the leading edge, and the center of the circle was obtained by translating the origin of the coordinate along the x-axis direction. The predicted SPL for each observing point of original hydrofoil as a function of frequency is shown in Figure 13. It can be seen that all of the curves exhibited similar peaks and dips and followed the same trend, which first increased in the lower frequency range and then decreased in the higher frequency range, except for the points near the trailing edge p1 and leading edge p9. To show the noise characteristics of the original hydrofoil and biomimetic models, four representative observing points were chosen, that is, p1, p5, p9, and p13. The comparisons of the predicted SPL spectra generated from the original and biomimetic model for the four observing points are illustrated in Figure 14. Additionally, the OASPL of the 16 observing points for the two models are presented in Figure 15.

Noise Characteristics along the Circumferential Direction
A total of 16 observing points were defined to evaluate the acoustic performance of near-field hydrofoil noise along the circumferential direction, which are presented in Figure 12. It can be seen that the observing points were evenly distributed in a circle with a radius of 30.00 mm. As is shown in Figure 12, the origin of the coordinate was at the midspan of the leading edge, and the center of the circle was obtained by translating the origin of the coordinate along the x-axis direction. The predicted SPL for each observing point of original hydrofoil as a function of frequency is shown in Figure 13. It can be seen that all of the curves exhibited similar peaks and dips and followed the same trend, which first increased in the lower frequency range and then decreased in the higher frequency range, except for the points near the trailing edge p 1 and leading edge p 9 . To show the noise characteristics of the original hydrofoil and biomimetic models, four representative observing points were chosen, that is, p 1 , p 5 , p 9 , and p 13 . The comparisons of the predicted SPL spectra generated from the original and biomimetic model for the four observing points are illustrated in Figure 14. Additionally, the OASPL of the 16 observing points for the two models are presented in Figure 15.

Noise Characteristics along the Circumferential Direction
A total of 16 observing points were defined to evaluate the acoustic performance of near-field hydrofoil noise along the circumferential direction, which are presented in Figure 12. It can be seen that the observing points were evenly distributed in a circle with a radius of 30.00 mm. As is shown in Figure 12, the origin of the coordinate was at the midspan of the leading edge, and the center of the circle was obtained by translating the origin of the coordinate along the x-axis direction. The predicted SPL for each observing point of original hydrofoil as a function of frequency is shown in Figure 13. It can be seen that all of the curves exhibited similar peaks and dips and followed the same trend, which first increased in the lower frequency range and then decreased in the higher frequency range, except for the points near the trailing edge p1 and leading edge p9. To show the noise characteristics of the original hydrofoil and biomimetic models, four representative observing points were chosen, that is, p1, p5, p9, and p13. The comparisons of the predicted SPL spectra generated from the original and biomimetic model for the four observing points are illustrated in Figure 14. Additionally, the OASPL of the 16 observing points for the two models are presented in Figure 15.   It can be observed from Figure 14 that the noise spectra generated from the original hydrofoil and biomimetic hydrofoil were totally different in these four observing points. The SPL generated from biomimetic model increased at low frequencies, and the frequency of the main peaks of the SPL curve all increased a little. As for the noise spectrum at higher frequency, it can obviously be seen that the biomimetic model would generate more noise than the original hydrofoil at observing point p1. Combined with the results of Figure 15, it was shown that the OASPLs of the biomimetic hydrofoil at observing points p1 and p9 were larger than those of the original model. Except for observing points p1 and p9, all the other observing points had positive effects on noise attenuation. The maximal noise reduction effect was found at observing point p2, with an overall noise reduction of 3.37 dB.  It can be observed from Figure 14 that the noise spectra generated from the original hydrofoil and biomimetic hydrofoil were totally different in these four observing points. The SPL generated from biomimetic model increased at low frequencies, and the frequency of the main peaks of the SPL curve all increased a little. As for the noise spectrum at higher frequency, it can obviously be seen that the biomimetic model would generate more noise than the original hydrofoil at observing point p 1 . Combined with the results of Figure 15, it was shown that the OASPLs of the biomimetic hydrofoil at observing points p 1 and p 9 were larger than those of the original model. Except for observing points p 1 and p 9 , all the other observing points had positive effects on noise attenuation. The maximal noise reduction effect was found at observing point p 2 , with an overall noise reduction of 3.37 dB. It can be observed from Figure 14 that the noise spectra generated from the original hydrofoil and biomimetic hydrofoil were totally different in these four observing points. The SPL generated from biomimetic model increased at low frequencies, and the frequency of the main peaks of the SPL curve all increased a little. As for the noise spectrum at higher frequency, it can obviously be seen that the biomimetic model would generate more noise than the original hydrofoil at observing point p1. Combined with the results of Figure 15, it was shown that the OASPLs of the biomimetic hydrofoil at observing points p1 and p9 were larger than those of the original model. Except for observing points p1 and p9, all the other observing points had positive effects on noise attenuation. The maximal noise reduction effect was found at observing point p2, with an overall noise reduction of 3.37 dB.

Noise Characteristics along the Radial Direction
Due to the high change rate of sound pressure in the near sound field, which can be found in Figure 11, the noise characteristics along the radial direction were investigated. Another six observing points were chosen, named p17, p18, p19, p20, p21, and p22, as shown in Figure 16. The interval between the neighboring observing points was 0.2 c. Correspondingly, the predicted SPL spectra and predicted OASPLs of the original hydrofoil and biomimetic hydrofoil at the six observing points are shown in Figure 17; Figure 18, respectively. It can be seen from Figure 17 that the main peaks in the SPL curve gradually faded away in the biomimetic design with the increase of observing distance. Further, the noise reduction effect in the high-frequency range gradually became more intense. As is shown in Figure 18, the overall noise reduction for the six observing points was 2.27 dB, 4.65 dB, 5.67 dB, 6.23 dB, 6.56 dB, and 6.78 dB, respectively. The six observing points had increasingly positive noise reduction effects, which was different from that of observing point p1.

Noise Characteristics along the Radial Direction
Due to the high change rate of sound pressure in the near sound field, which can be found in Figure 11, the noise characteristics along the radial direction were investigated. Another six observing points were chosen, named p 17 , p 18 , p 19 , p 20 , p 21 , and p 22 , as shown in Figure 16. The interval between the neighboring observing points was 0.2 c. Correspondingly, the predicted SPL spectra and predicted OASPLs of the original hydrofoil and biomimetic hydrofoil at the six observing points are shown in Figure 17; Figure 18, respectively.

Noise Characteristics along the Radial Direction
Due to the high change rate of sound pressure in the near sound field, which can be found in Figure 11, the noise characteristics along the radial direction were investigated. Another six observing points were chosen, named p17, p18, p19, p20, p21, and p22, as shown in Figure 16. The interval between the neighboring observing points was 0.2 c. Correspondingly, the predicted SPL spectra and predicted OASPLs of the original hydrofoil and biomimetic hydrofoil at the six observing points are shown in Figure 17; Figure 18, respectively. It can be seen from Figure 17 that the main peaks in the SPL curve gradually faded away in the biomimetic design with the increase of observing distance. Further, the noise reduction effect in the high-frequency range gradually became more intense. As is shown in Figure 18, the overall noise reduction for the six observing points was 2.27 dB, 4.65 dB, 5.67 dB, 6.23 dB, 6.56 dB, and 6.78 dB, respectively. The six observing points had increasingly positive noise reduction effects, which was different from that of observing point p1. It can be seen from Figure 17 that the main peaks in the SPL curve gradually faded away in the biomimetic design with the increase of observing distance. Further, the noise reduction effect in the high-frequency range gradually became more intense. As is shown in Figure 18, the overall noise reduction for the six observing points was 2.27 dB, 4.65 dB, 5.67 dB, 6.23 dB, 6.56 dB, and 6.78 dB, respectively. The six observing points had increasingly positive noise reduction effects, which was different from that of observing point p 1 .

Analysis of Far-Field Hydrodynamic Noise
To quantify the noise reduction of the biomimetic hydrofoil in the far sound field, another 16 observing points were chosen. As is shown in Figure 19, the observing points were evenly distributed in a circle with a radius of 300.00 mm. In order to observe the comparisons of the two models clearly, the observing points at 0 • , 90 • , 180 • , and 270 • were chosen, which was similar to the conduction of the near sound field in Section 6.2.1. The comparisons of the predicted SPL spectra generated from the original hydrofoil and biomimetic models are shown in Figure 20. And the OASPLs of all the far sound field observing points are demonstrated in Figure 21. To quantify the noise reduction of the biomimetic hydrofoil in the far sound field, another 16 observing points were chosen. As is shown in Figure 19, the observing points were evenly distributed in a circle with a radius of 300.00 mm. In order to observe the comparisons of the two models clearly, the observing points at 0°, 90°, 180°, and 270° were chosen, which was similar to the conduction of the near sound field in Section 6.2.1. The comparisons of the predicted SPL spectra generated from the original hydrofoil and biomimetic models are shown in Figure 20. And the OASPLs of all the far sound field observing points are demonstrated in Figure 21. It is easily observed from Figure 20 that the main peaks of the original model were eliminated after carrying out biomimetic design. Combined with the result of Figure 21, it can be seen that all of the observing points have achieved a certain degree of noise reduction effect. It is worth noting that the observing points at 0° and 180°, where negative effects could be found in the near sound field, had significant noise reduction effects in the far sound field. Among all of the 16 observing points, the maximal noise reduction was 7.28 dB, obtained at observing point p23.  It is easily observed from Figure 20 that the main peaks of the original model were eliminated after carrying out biomimetic design. Combined with the result of Figure 21, it can be seen that all of the observing points have achieved a certain degree of noise reduction effect. It is worth noting that the observing points at 0 • and 180 • , where negative effects could be found in the near sound field, had significant noise reduction effects in the far sound field. Among all of the 16 observing points, the maximal noise reduction was 7.28 dB, obtained at observing point p 23 . To quantify the noise reduction of the biomimetic hydrofoil in the far sound field, another 16 observing points were chosen. As is shown in Figure 19, the observing points were evenly distributed in a circle with a radius of 300.00 mm. In order to observe the comparisons of the two models clearly, the observing points at 0°, 90°, 180°, and 270° were chosen, which was similar to the conduction of the near sound field in Section 6.2.1. The comparisons of the predicted SPL spectra generated from the original hydrofoil and biomimetic models are shown in Figure 20. And the OASPLs of all the far sound field observing points are demonstrated in Figure 21. It is easily observed from Figure 20 that the main peaks of the original model were eliminated after carrying out biomimetic design. Combined with the result of Figure 21, it can be seen that all of the observing points have achieved a certain degree of noise reduction effect. It is worth noting that the observing points at 0° and 180°, where negative effects could be found in the near sound field, had significant noise reduction effects in the far sound field. Among all of the 16 observing points, the maximal noise reduction was 7.28 dB, obtained at observing point p23.

Noise Reduction Mechanism
The hydrofoil self-noise is produced by the interaction between the structure and the turbulence produced in its own boundary layer and the near wake [6]. To give a clear view of the flow field, the two-dimensional contours of the flow field at the mid plane of the hydrofoil at 2.0 s are given in this section.
The comparison of the pressure contours between the original hydrofoil and biomimetic hydrofoil is shown in Figure 22. It can be seen that the pressure of the biomimetic model was much more stable than the original one, especially near the region of suction surface. Moreover, the contours of velocities in the y direction and in the flow direction at the trailing edge of the two models investigated in this paper are illustrated in Figure 23; Figure 24, respectively. It is noticeable that the velocity fluctuations in both the y direction and in the flow direction at the trailing edge improved a lot after adding spanwise microgrooved surfaces. The velocity gradient of the flow in the wake of the trailing edge was significantly decreased. To give more explanations for the mechanism of the noise reduction demonstrated in this paper, the streamlines in the groove are exhibited in Figure 25.

Noise Reduction Mechanism
The hydrofoil self-noise is produced by the interaction between the structure and the turbulence produced in its own boundary layer and the near wake [6]. To give a clear view of the flow field, the two-dimensional contours of the flow field at the mid plane of the hydrofoil at 2.0 s are given in this section.
The comparison of the pressure contours between the original hydrofoil and biomimetic hydrofoil is shown in Figure 22. It can be seen that the pressure of the biomimetic model was much more stable than the original one, especially near the region of suction surface. Moreover, the contours of velocities in the y direction and in the flow direction at the trailing edge of the two models investigated in this paper are illustrated in Figure 23; Figure 24, respectively. It is noticeable that the velocity fluctuations in both the y direction and in the flow direction at the trailing edge improved a lot after adding spanwise microgrooved surfaces. The velocity gradient of the flow in the wake of the trailing edge was significantly decreased. To give more explanations for the mechanism of the noise reduction demonstrated in this paper, the streamlines in the groove are exhibited in Figure 25.

Noise Reduction Mechanism
The hydrofoil self-noise is produced by the interaction between the structure and the turbulence produced in its own boundary layer and the near wake [6]. To give a clear view of the flow field, the two-dimensional contours of the flow field at the mid plane of the hydrofoil at 2.0 s are given in this section.
The comparison of the pressure contours between the original hydrofoil and biomimetic hydrofoil is shown in Figure 22. It can be seen that the pressure of the biomimetic model was much more stable than the original one, especially near the region of suction surface. Moreover, the contours of velocities in the y direction and in the flow direction at the trailing edge of the two models investigated in this paper are illustrated in Figure 23; Figure 24, respectively. It is noticeable that the velocity fluctuations in both the y direction and in the flow direction at the trailing edge improved a lot after adding spanwise microgrooved surfaces. The velocity gradient of the flow in the wake of the trailing edge was significantly decreased. To give more explanations for the mechanism of the noise reduction demonstrated in this paper, the streamlines in the groove are exhibited in Figure 25. As is shown in Figure 25, two pairs of vortices rotating reversely existed in the interior of the microgrooves. It is noticeable that the rotation direction of the larger vortex near the surface of the hydrofoil was favorable to the flow of the external flow field. The "micro-air bearing" introduced by Pan [55] could provide explanations for this phenomenon in this paper. The advantages of surfaces with spanwise microgrooves for retarding the development of the turbulent boundary layer include the following: (1) weakening of the streamwise vortices that spawn the secondary vortices and (2) retention of low-speed fluid within the grooves. For the present cases, the spanwise microgrooves, which were introduced onto the surface of the hydrofoils, consumed the energy of the flow. The existence of these secondary flow vortices plays a similar role to "roller bearings". These secondary flow vortices are stable in the groove and will not appear on the surface, which would lead to the phenomenon of vortices scattering. In short, the basic idea of underwater noise reduction of hydrofoils with spanwise microgrooved surface is to hinder the process of turbulence, reduce the flow energy, and weaken the intensity of turbulent burst.

Conclusions
A biomimetic hydrofoil is proposed in this study with spanwise microgrooved surfaces inspired by sharkskin on the FF-77-W airfoil profile. Firstly, LES were carried out to simulate hydrodynamic characteristics of the hydrofoil; then, the sound pressure at the observing points was calculated based on the FW-H equation. After verifying and validating the numerical method, the hydrodynamic As is shown in Figure 25, two pairs of vortices rotating reversely existed in the interior of the microgrooves. It is noticeable that the rotation direction of the larger vortex near the surface of the hydrofoil was favorable to the flow of the external flow field. The "micro-air bearing" introduced by Pan [55] could provide explanations for this phenomenon in this paper. The advantages of surfaces with spanwise microgrooves for retarding the development of the turbulent boundary layer include the following: (1) weakening of the streamwise vortices that spawn the secondary vortices and (2) retention of low-speed fluid within the grooves. For the present cases, the spanwise microgrooves, which were introduced onto the surface of the hydrofoils, consumed the energy of the flow. The existence of these secondary flow vortices plays a similar role to "roller bearings". These secondary flow vortices are stable in the groove and will not appear on the surface, which would lead to the phenomenon of vortices scattering. In short, the basic idea of underwater noise reduction of hydrofoils with spanwise microgrooved surface is to hinder the process of turbulence, reduce the flow energy, and weaken the intensity of turbulent burst. As is shown in Figure 25, two pairs of vortices rotating reversely existed in the interior of the microgrooves. It is noticeable that the rotation direction of the larger vortex near the surface of the hydrofoil was favorable to the flow of the external flow field. The "micro-air bearing" introduced by Pan [55] could provide explanations for this phenomenon in this paper. The advantages of surfaces with spanwise microgrooves for retarding the development of the turbulent boundary layer include the following: (1) weakening of the streamwise vortices that spawn the secondary vortices and (2) retention of low-speed fluid within the grooves. For the present cases, the spanwise microgrooves, which were introduced onto the surface of the hydrofoils, consumed the energy of the flow. The existence of these secondary flow vortices plays a similar role to "roller bearings". These secondary flow vortices are stable in the groove and will not appear on the surface, which would lead to the phenomenon of vortices scattering. In short, the basic idea of underwater noise reduction of hydrofoils with spanwise microgrooved surface is to hinder the process of turbulence, reduce the flow energy, and weaken the intensity of turbulent burst.

Conclusions
A biomimetic hydrofoil is proposed in this study with spanwise microgrooved surfaces inspired by sharkskin on the FF-77-W airfoil profile. Firstly, LES were carried out to simulate hydrodynamic characteristics of the hydrofoil; then, the sound pressure at the observing points was calculated based on the FW-H equation. After verifying and validating the numerical method, the hydrodynamic

Conclusions
A biomimetic hydrofoil is proposed in this study with spanwise microgrooved surfaces inspired by sharkskin on the FF-77-W airfoil profile. Firstly, LES were carried out to simulate hydrodynamic characteristics of the hydrofoil; then, the sound pressure at the observing points was calculated based