Noise Characteristics Analysis of the Horizontal Axis Hydrokinetic Turbine Designed for Unmanned Underwater Mooring Platforms

: Operating horizontal axis hydrokinetic turbine (HAHT) generates noise a ﬀ ecting the ocean environment adversely. Therefore, it is essential to determine the noise characteristics of such types of HAHT, as large-scale turbine sets would release more noise pollution to the ocean. Like other rotating machinery, the hydrodynamic noise generated by the rotating turbine has been known to be the most important noise source. In the present work, the transient turbulent ﬂow ﬁeld of the HAHT is obtained by incompressible large eddy simulation, thereafter, the Ffowcs Williams and Hawkings acoustic analogy formulation is carried out to predict the noise generated from the pressure ﬂuctuations of the blade surface. The coe ﬃ cient of power is compared with the experimental results, with a good agreement being achieved. It is seen from the pressure contours that the 80% span of the blade has the most severe pressure ﬂuctuations, which concentrate on the region of leading the edge of the airfoil and the suction surface of the airfoil. Then, the noise characteristics around a single turbine are systematically studied, in accordance with the results of the ﬂow ﬁeld. The noise characteristics around the whole turbine are also investigated to determine the directionality of the noise emission of HAHT.


Introduction
Unmanned underwater mooring platforms (UUMPs), which are anchored to the seabed using mooring cables, could perform various functions in the ocean with a long period of time, typically ranging from several months to several years. In the present stage, most of the battery-powered UUMPs face the problem of an insufficiency if energy supply during the service period, because of the limited battery capacity and uninterrupted energy consumption by the electronic devices. Thus, the energy problem of the UUMPs is a significant issue to be considered in the design process.
Marine renewable energy is a huge potential energy resource, including ocean thermal energy [1], ocean surface solar energy [2], ocean wave energy [3], and ocean current energy [4], which could potentially recharge the batteries of the UUMPs. As most of the UUMPs are expected to carry out fix-point missions for long durations, ocean thermal energy and ocean surface solar energy are not suitable for continuous energy supply of UUMPs. For the utilization of ocean wave energy, conversion devices are also expected to be located near the ocean surface, which is infeasible for deep-water UUMPs. As a matter of fact, UUMPs are often placed where ocean currents are consistently available [5], which could provide an effective energy supply to recharge the UUMPs.
Hydrokinetic turbines are the primary energy conversion devices for the utilization of ocean current energy, which can be classified into two categories, namely vertical axis hydrokinetic turbine (VAHT) and horizontal axis hydrokinetic turbine (HAHT; they cloud be also be called cross-flow turbines and axial-flow turbines) [6]. In the early stages, Tian et al. [4] designed a miniature VAHT to capture ocean current energy and recharge a moored autonomous underwater vehicle (AUV). However, the maximum power coefficient of the turbine is very low (only 0.1) because of the disturbance of the hull. Subsequently, they investigated several types of HAHT [5,7,8], and the results showed that the maximum power coefficient has significantly improved. The investigation of HAHT is still a hot topic all over the world [9,10].
Although extensive research and investigations over the last decade on ocean current energy have been carried out, there are still many challenges to be encountered. One of the most important aspects is that the installation of the ocean current turbines has to fully consider the interactions with the ambient marine environment. For example, several researchers have conducted research about the impact of hydrokinetic turbines on the spatial ecology of wild fishes [11], fish migration [12], phytoplankton dynamics [13], etc. Another annoying issue that has to be faced is the hydrodynamic noise when the turbine is in operational mode. Loud flow-induced noise not only affects the health and existence of the organisms in the ocean, but also threatens the safety of the high-precision mechanical devices of the UUMPs. Recently, some scholars have carried out the investigation of noise produced by hydrokinetic turbines, as follows: Hafla et al. [14] modelled underwater noise propagation from marine hydrokinetic power devices through a time-domain, velocity-pressure solution' Lossent et al. [15] measured the underwater operational noise level emitted by an ocean current turbine and assessed its potential impact on marine fauna, which indicated that the behavioral disturbance might occur up to 1 km around the device for harbor porpoises; and Pine et al. [16] found that the turbine-generated noise can reduce the listening space of marine fauna. The results showed that listening space reductions (LSRs) for harbor seals were in excess of 80% within 60 m, whilst for harbor porpoises, they were in excess of 55% within 10 m of the devices. Like wind turbines, the vertical axis turbines are much quieter compared with the horizontal types, and the horizontal axis turbines often generate low-frequency noise at blade passing frequency and higher order harmonics, which are harmful to the environment [17]. This phenomenon has also been confirmed by the research group from Northwestern Polytechnical University, with a series of continuous research [4,5,7,8,18]. In fact, the HAHT investigated in this paper and the horizontal axis wind turbines are similar to the rotating machinery of renewable energy. For wind turbines, several researchers have investigated the noise characteristics. A technical procedure [19,20] was proposed to simultaneously estimate the emission and the residual noise components measured nearby a wind farm when the residual noise is mainly generated by wind. The experimental tests of the noise emission characteristics of large wind turbines, with emphasis on infrasound and low-frequency noise [21], were also carried out. Additionally, the annoyance of people because of wind turbine noise [22] and the wind turbine noise reduction [23] have been investigated, individually. To the author's knowledge, the specific noise characteristics of the ocean current turbines have not been systematically studied. Because of the loud deleterious hydrodynamic noise observed around HAHT, it is necessary to consider the noise emission problem of ocean current energy generating units, which not only help to build a harmonious marine ecological environment, but also enhances the competitiveness of current generating units in the international market.
In this paper, the noise characteristics of a HAHT (rated power = 30.0 W, diameter D = 1.2 m) designed for UUMPs are investigated based on large eddy simulations (LES) combined with the Ffowcs Williams and Hawkings (FW-H) equation, which is adopted to simulate the rotating noise of turbines [24]. The paper is organized as follows. In Section 2, the model of the HAHT investigated in this paper is introduced. Then, the numerical model is illustrated in Section 3, and the accuracy of the numerical model predictions is verified and validated in Section 4. In Section 5, the results and discussions of the flow field characteristics and the noise field characteristics are demonstrated. Finally, the conclusions are presented in Section 6.

Description of the HAHT
The HAHT designed in this paper is a three-blade turbine with a diameter of D = 1.2 m. As shown in Figure 1, the UUMPs are anchored to the sea floor with mooring cables, and the HAHT is mounted on the UUMPs with a hub, gearing system, bearing system, and sealing system. As the ocean current flows, the torque of the turbine is transmitted to the generator via the magnetic coupling unit, where the energy conversion process could be accomplished from kinetic energy to electric energy. Figure 1 simply illustrates the diagrammatic sketch of the working principle of the HAHT designed for UUMPs. As this paper is focused on the noise characteristics of the HAHT, therefore, the details of the magnetic coupling unit, generator, and gearbox are not introduced here, which could be found in the previous study [7,8].

Description of the HAHT
The HAHT designed in this paper is a three-blade turbine with a diameter of D = 1.2 m. As shown in Figure 1, the UUMPs are anchored to the sea floor with mooring cables, and the HAHT is mounted on the UUMPs with a hub, gearing system, bearing system, and sealing system. As the ocean current flows, the torque of the turbine is transmitted to the generator via the magnetic coupling unit, where the energy conversion process could be accomplished from kinetic energy to electric energy. Figure 1 simply illustrates the diagrammatic sketch of the working principle of the HAHT designed for UUMPs. As this paper is focused on the noise characteristics of the HAHT, therefore, the details of the magnetic coupling unit, generator, and gearbox are not introduced here, which could be found in the previous study [7,8]. The blade of the turbine is generated from the FF-77-W airfoil shape. More concretely, it is created from the fx77w121, fx77w153, and fx77w258 airfoils, which are available in the literature [25]. The details of the geometrical parameters of the airfoil are summarized in Table 1, where R is the radius of the turbine, r is the spanwise location of the blade, c is the chord length, t is the thickness of the airfoil, and φ is the twist angle of the airfoil. To give a closer view of the geometric appearance of the blade, the three-dimensional view of the blade, and the view of the blade from the tip of the blade to blade root are shown in Figure 2.  The blade of the turbine is generated from the FF-77-W airfoil shape. More concretely, it is created from the fx77w121, fx77w153, and fx77w258 airfoils, which are available in the literature [25]. The details of the geometrical parameters of the airfoil are summarized in Table 1, where R is the radius of the turbine, r is the spanwise location of the blade, c is the chord length, t is the thickness of the airfoil, and ϕ is the twist angle of the airfoil. To give a closer view of the geometric appearance of the blade, the three-dimensional view of the blade, and the view of the blade from the tip of the blade to blade root are shown in Figure 2.

Numerical Model
The procedure of noise calculation for the HAHT could be represented in four steps, as in Figure  3. The geometry modeling was carried out to build the flow field of the turbine initially, and then the mesh generation was conducted to meet the grid distribution requirements for the subsequent flow field analysis and noise field analysis. The flow field characteristics of the rotational turbine were simulated using LES with a transient computational fluid dynamics (CFD) method. Thereafter, the FW-H acoustic analogy was employed with the pressure fluctuations near the surface of the turbine as the input to calculate the hydrodynamic noise result at the designated sound observing point.

LES Formulation
A three-dimensional unsteady incompressible simulation was used for the LES model. In the procedure of LES calculation, the low-pass filtered quantity of turbulent vortices was resolved by establishing the turbulence model, while the sub-grid scale quantity of the turbulent vortices, which were smaller than the size of the computational mesh, was modeled. The filtered variable () x is expressed as follows: where x is the large-scale spatial coordinate after filtering, V is the volume of a computational cell, and x ' is the spatial coordinate of the actual flow. The unsteady filtered mass conservation equation and the filtered incompressible transient three-dimensional Navier-Stokes equations are as follows:

Numerical Model
The procedure of noise calculation for the HAHT could be represented in four steps, as in Figure 3. The geometry modeling was carried out to build the flow field of the turbine initially, and then the mesh generation was conducted to meet the grid distribution requirements for the subsequent flow field analysis and noise field analysis. The flow field characteristics of the rotational turbine were simulated using LES with a transient computational fluid dynamics (CFD) method. Thereafter, the FW-H acoustic analogy was employed with the pressure fluctuations near the surface of the turbine as the input to calculate the hydrodynamic noise result at the designated sound observing point.

Numerical Model
The procedure of noise calculation for the HAHT could be represented in four steps, as in Figure  3. The geometry modeling was carried out to build the flow field of the turbine initially, and then the mesh generation was conducted to meet the grid distribution requirements for the subsequent flow field analysis and noise field analysis. The flow field characteristics of the rotational turbine were simulated using LES with a transient computational fluid dynamics (CFD) method. Thereafter, the FW-H acoustic analogy was employed with the pressure fluctuations near the surface of the turbine as the input to calculate the hydrodynamic noise result at the designated sound observing point.

LES Formulation
A three-dimensional unsteady incompressible simulation was used for the LES model. In the procedure of LES calculation, the low-pass filtered quantity of turbulent vortices was resolved by establishing the turbulence model, while the sub-grid scale quantity of the turbulent vortices, which were smaller than the size of the computational mesh, was modeled. The filtered variable () x is expressed as follows: where x is the large-scale spatial coordinate after filtering, V is the volume of a computational cell, and x ' is the spatial coordinate of the actual flow. The unsteady filtered mass conservation equation and the filtered incompressible transient three-dimensional Navier-Stokes equations are as follows:

LES Formulation
A three-dimensional unsteady incompressible simulation was used for the LES model. In the procedure of LES calculation, the low-pass filtered quantity of turbulent vortices was resolved by establishing the turbulence model, while the sub-grid scale quantity of the turbulent vortices, which were smaller than the size of the computational mesh, was modeled. The filtered variable ϕ(x) is expressed as follows: where x is the large-scale spatial coordinate after filtering, V is the volume of a computational cell, and x is the spatial coordinate of the actual flow. The unsteady filtered mass conservation equation and the filtered incompressible transient three-dimensional Navier-Stokes equations are as follows: where ρ is the density, u i and u j are the velocity components, p is the pressure, and τ ij is the term of sub-grid stress, which includes the residual stresses, represented as follows: The sub-grid stress τ ij indicates the momentum transport between large-scale vortex and small-scale vortex. What is more, the effect of small-scale fluctuations on the overall flow field is quantified. To make Equation (3) enclosed, the corresponding sub-grid model has to be reconstructed. At present, the sub-grid model is based largely on the hypothesis proposed by Boussinesq [26], expressed as follows: where µ t is the coefficient of the sub-grid eddy viscosity, δ ij is the sub-grid scale of Reynolds stress, and S ij is the deformation rate tensor corresponding to the scale of solution, written as follows: The eddy viscosity hypothesis of Equation (5) would transform solving the sub-grid stress into solving the coefficient of sub-grid scale eddy viscosity. The Smagorinsky-Lilly model is developed to solve this problem, where the coefficient of the sub-grid eddy viscosity could be obtained as follows: where L s is the mixing length of sub-grid, which can be expressed as follows: 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, which is computed according to the volume of the computational cell, as follows:

Acoustic Simulation FW-H Equation
Lighthill [27] was the earliest researcher concentrated on aerodynamic noise, whose idea was rearranging the Navier-Stokes equations into an inhomogeneous wave equation, as follows: where c 0 is the sound speed, Q ij is Lighthill's stress tensor. It was used to predict the noise generated from jet noise initially. Unfortunately, Lighthill only considered the quadrupole sources. Thereafter, Ffowcs-Williams and Hawkings [28] extended Equation (11) to predict the noise from a moving surface, where the monopole and dipole sources are also considered, as follows: where H( f ) is the Heaviside function, n j is the unit normal vector pointing toward the exterior region, u i is the fluid velocity in the i direction, u n is the fluid velocity component normal to the surface f = 0 (f = 0 represents a mathematical surface introduced to "embed" the exterior flow problem in an unbounded space [29]), v n is the surface velocity normal to the surface f = 0, δ( f ) is the Dirac delta function, and P ij is the compressive stress tensor, expressed as follows: Lighthill's stress tensor Q ij could be written as follows Equation (12) could be integrated analytically on the assumption that the flow space is free and there are no obstacles between the noise sources and the receivers. With the generalized function theory [30] and free space Green's function, the overall sound pressure at the observing points could be obtained, as follows: where p T and p L denote the sound pressure of the thickness noise (monopole source) and loading noise (dipole source), respectively, expressed as follows [31]: where where r is the distance to receiver, t is the observer time, and M r is the Mach number of a point on the moving surface the Mach number vector. M i relates to the motion of the integration surface: The dot above a variable represents the source-time differentiation of that variable, the subscripts represent the inner products of a vector and a unit vector implied by the subscript, and the subscript ret indicates that the integrals are calculated at the retarded time, defined as follows: After a Fast Fourier Transform (FFT) algorithm, the sound pressure data at a noise observing point could be obtained in the considered frequency range. The sound pressure level (SPL) is obtained as follows: SPL(dB) = 20 log 10 p rms p re f (21) where p rms is the effective acoustic pressure amplitude of the sound wave and p re f is the reference acoustic pressure. For noise transmission in water, the reference pressure is normally defined as 1µPa(1 × 10 −6 Pa). The overall sound pressure level (OASPL) is defined as follows: 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.

Computation Domains and Grid Generation
To consider the full development of the incoming flow and to eliminate the blockage effect, the overall computational domain is a circular cylinder with a diameter of 5D and a length of 15D. The HAHT was placed at the centerline of the cylinder at a distance of 5D from the inlet boundary, as shown in Figure 4. The overall computational domain was divided into three subdomains. The first domain was called the rotating flow domain, which contained the grid elements surrounding the HAHT. The second domain was called wake flow domain, which was a little larger than the first domain, to simulate the flow characteristics of the wake flow after the HAHT. The third domain which contains the cells in the outer region was called outer flow domain.  (21) where rms p is the effective acoustic pressure amplitude of the sound wave and ref p is the reference acoustic pressure. For noise transmission in water, the reference pressure is normally defined as 1 μPa (  -6 1 10 Pa ). The overall sound pressure level (OASPL) is defined as follows: OASPL(dB) 10 log 10 (22) where i SPL( ) 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.

Computation Domains and Grid Generation
To consider the full development of the incoming flow and to eliminate the blockage effect, the overall computational domain is a circular cylinder with a diameter of 5D and a length of 15D. The HAHT was placed at the centerline of the cylinder at a distance of 5D from the inlet boundary, as shown in Figure 4. The overall computational domain was divided into three subdomains. The first domain was called the rotating flow domain, which contained the grid elements surrounding the HAHT. The second domain was called wake flow domain, which was a little larger than the first domain, to simulate the flow characteristics of the wake flow after the HAHT. The third domain which contains the cells in the outer region was called outer flow domain.  Figure 4. A flow with a steady and uniform velocity of 0.5 m/s, which is considered here according to the marine current energy resources distribution in China [32], was imposed at the inlet boundary of the computational domain. A pressure outlet boundary was applied at the outlet boundary of the computational domain. No-slip boundary conditions were applied at the surface of the blades. To enhance the stability of the numerical simulations, symmetry boundary conditions were imposed at the side wall of the outer flow domain, which could help to avoid the effects of walls. Two pairs of interfaces, interface groups 1 and 2, were applied at the overlap faces between adjacent subdomains to consider the data transmission of the flow properties.
To simulate the flow characteristics of the HAHT accurately, fine meshes were generated for the region of the rotating flow domain and wake flow domain. Figure   The boundary conditions are shown in Figure 4. A flow with a steady and uniform velocity of 0.5 m/s, which is considered here according to the marine current energy resources distribution in China [32], was imposed at the inlet boundary of the computational domain. A pressure outlet boundary was applied at the outlet boundary of the computational domain. No-slip boundary conditions were applied at the surface of the blades. To enhance the stability of the numerical simulations, symmetry boundary conditions were imposed at the side wall of the outer flow domain, which could help to avoid the effects of walls. Two pairs of interfaces, interface groups 1 and 2, were applied at the overlap faces between adjacent subdomains to consider the data transmission of the flow properties.
To simulate the flow characteristics of the HAHT accurately, fine meshes were generated for the region of the rotating flow domain and wake flow domain. Figure 5a shows the mesh on the surface of the computational model from A-A plane. A hex dominant method was applied in the wake flow domain and outer flow domain, while the tetrahedrons mesh was employed in the region of the rotating flow domain. Fine mesh in the wake flow domain could capture the flow information in the wake flow. Extremely refined mesh in the rotating flow domain could obtain the flow fluctuations near the blades. To capture the pressure fluctuations of the surface of the blade, which lead to the generation of the flow-induced noise, 20 layers of prisms were imposed on the surface of the blades with a mesh growth rate of 1.12. The mesh of the blade subdomain in three different cross sections with the view from blade tip to blade root is shown in Figure 5b. To observe the mesh near the blade clearly, Figure 5c gives a view of the mesh near the blade from the cross section along the spanwise direction. generation of the flow-induced noise, 20 layers of prisms were imposed on the surface of the blades with a mesh growth rate of 1.12. The mesh of the blade subdomain in three different cross sections with the view from blade tip to blade root is shown in Figure 5b. To observe the mesh near the blade clearly, Figure 5c gives a view of the mesh near the blade from the cross section along the spanwise direction.

Solution Sets
Fluent 14.5 was used as the CFD software in the simulations. According to the result of Tian et al. [8], the HAHT could obtain an optimal value of the power coefficient when the tip speed ratio (TSR) was nearly equal to 5, which was adopted in the present simulations. TSR denotes the ratio of the blade tip speed to the incoming flow speed, which is defined as follows: where  is the angular velocity of the HAHT and U is the flow speed at the inlet. For the LES model, the Smagorinsky-Lilly model was chosen as its subgrid-scale model, where the Cs value was 0.1, to yield the best results for a wide range of flows [29]. To obtain a fast convergence, the Pressure-Implicit with Splitting of Operators (PISO) algorithm was applied in pressure-velocity coupling in the solver. Additionally, the second order implicit was chosen as the transient formulation in the transient process. The drops of all of the scaled residuals (continuity,

Solution Sets
Fluent 14.5 was used as the CFD software in the simulations. According to the result of Tian et al. [8], the HAHT could obtain an optimal value of the power coefficient when the tip speed ratio (TSR) was nearly equal to 5, which was adopted in the present simulations. TSR denotes the ratio of the blade tip speed to the incoming flow speed, which is defined as follows: where ω is the angular velocity of the HAHT and U is the flow speed at the inlet.
For the LES model, the Smagorinsky-Lilly model was chosen as its subgrid-scale model, where the Cs value was 0.1, to yield the best results for a wide range of flows [29]. To obtain a fast convergence, the Pressure-Implicit with Splitting of Operators (PISO) algorithm was applied in pressure-velocity coupling in the solver. Additionally, the second order implicit was chosen as the transient formulation in the transient process. The drops of all of the scaled residuals (continuity, x-velocity, y-velocity, and z-velocity) below 1.0 × 10 −4 were used as the convergence criterion. The time step for each simulation in previous research was 1 • of the HAHT rotation [8], while the time step for the present simulation was set as 0.25 • of the HAHT rotation in order to obtain a wider frequency range, and to ensure a small Courant number and accurate temporal discretization at the same time. After the dynamic stabilization of the LES calculation, another 1800 time steps were set to enable the acoustic frequency range calculated to be 0~477 Hz with an interval of about 0.5 Hz. The surface pressure fluctuations of the HAHT were used as the noise source field data, and then the noise spectrum in the frequency domain could be obtained using FFT.

Numerical Method Verification and Validation
A mesh resolution study was carried out to evaluate the effect of mesh resolution on the torque of the HAHT and its coefficient of the power. The coefficient of the power is calculated as follows where C p is coefficient of the power, M is the torque of the HAHT, and ω is the angular velocity of the HAHT. The results with different mesh resolutions are summarized in Table 2. It is seen that three different kinds of mesh configurations were investigated. The main different settings for the mesh generation were as follows:  Unfortunately, there are no published results for the noise validation of the horizontal axis hydrokinetic turbine considered in this study, but the hydrodynamic results for this case are available. To confirm the accuracy of the present simulation, the coefficient of power was compared with the experimental data of Bahaj et al. [33].
In the experimental test, a three-blade marine current turbine with 800 mm in diameter was tested in a 2.4 m × 1.2 m cavitation tunnel and a 60 m towing tank. The blades of the turbines were developed from the profile shape of NACA 63-8XX. The inflow velocity of 1.73 m/s was chosen to be compared with. The experimental results provided a series of results with a TSR from 4.1 to 11.3, and the TSR of integers from 4 to 10 were simulated by the LES method to investigate the torque output. Figure 6 shows the comparison of the coefficient of power between the present LES model and the experimental tests [33]. It can be seen that the results of the two methods followed the same trend, and the relative errors were within 4% for the cases of TSR = 5, TSR = 6, TSR = 7, and TSR = 8. The errors seemed to be larger when TSR = 9 and TSR = 10. Overall, the comparison between the present LES method and the experimental result shows an acceptable accuracy for this LES computation.

Flow Field Characteristics
The HAHT was designed for optimal efficiency at a constant revolution per minute (RPM) of 40, relating to the TSR of five at an incoming flow of 0.5 m/s. Therefore, the noise characteristics of the HAHT working at TSR = 5 were investigated preliminarily. To complete the numerical simulations, the material properties of the water utilized in the analysis are given in Table 3. Table 3. Material properties of water adopted in this paper.

Description
Parameter After the LES computation, the vortices generated around the HAHT are shown in Figure 7. It can be seen that plenty of vortices were generated after the turbine when the turbine rotated. As the hydrodynamic noise was closely related to the pressure fluctuations in the flow field, the pressure distribution around the turbine had to be emphasized. To observe the pressure characteristics clearly, the pressure contours at several cross sections of the blade are given in Figure 8. In Figure 8, the pressure range of the five cross sections is identical. It can be seen that the pressure fluctuations were gradually becoming more serious with the increase of the spanwise length from a 20% span to 80% span. At the cross sections of the 100% span, the pressure fluctuations weakened for both the barotropic region (red areas) and the negative pressure region (blue areas). From the pressure contours of the five cross sections of the blade, it is noticeable that the maximum absolute values of the pressure concentrate on the region of the leading edge of the airfoil and the suction surface of the airfoil. Unlike the other four cross sections of the airfoil, the negative pressure region of the 100% span spread over the suction surface from the leading edge to the trailing edge. To better investigate the hydrodynamic noise of the blade of the HAHT designed for UUMPs, its noise characteristics were subsequently investigated.

Flow Field Characteristics
The HAHT was designed for optimal efficiency at a constant revolution per minute (RPM) of 40, relating to the TSR of five at an incoming flow of 0.5 m/s. Therefore, the noise characteristics of the HAHT working at TSR = 5 were investigated preliminarily. To complete the numerical simulations, the material properties of the water utilized in the analysis are given in Table 3. Table 3. Material properties of water adopted in this paper.

Description Parameter and Value
Temperature T = 20 • C Density ρ = 998.2 kg/m 3 Viscosity µ = 0.001003 kg/(m·s) Acoustic speed v = 1483 m/s Reference sound pressure p re f = 1 × 10 −6 Pa After the LES computation, the vortices generated around the HAHT are shown in Figure 7. It can be seen that plenty of vortices were generated after the turbine when the turbine rotated. As the hydrodynamic noise was closely related to the pressure fluctuations in the flow field, the pressure distribution around the turbine had to be emphasized. To observe the pressure characteristics clearly, the pressure contours at several cross sections of the blade are given in Figure 8. In Figure 8, the pressure range of the five cross sections is identical. It can be seen that the pressure fluctuations were gradually becoming more serious with the increase of the spanwise length from a 20% span to 80% span. At the cross sections of the 100% span, the pressure fluctuations weakened for both the barotropic region (red areas) and the negative pressure region (blue areas). From the pressure contours of the five cross sections of the blade, it is noticeable that the maximum absolute values of the pressure concentrate on the region of the leading edge of the airfoil and the suction surface of the airfoil. Unlike the other four cross sections of the airfoil, the negative pressure region of the 100% span spread over the suction surface from the leading edge to the trailing edge. To better investigate the hydrodynamic noise of the blade of the HAHT designed for UUMPs, its noise characteristics were subsequently investigated.

Arrangement of Noise Observation Location
First of all, the noise characteristics of a single blade were investigated. Figure 9 illustrates the noise observing points arrangement around a single blade. The acoustic monitoring points for the SPL computation were located at five different spanwise positions with h = 0.2 R, h = 0.4 R, h = 0.6 R, h = 0.8 R, and h = 1.0 R. For each spanwise position, there were two groups of observing points forming two circles with a radius r1 = 40.0 mm and r2 = 100.0 mm. The centers of the circles were at the aerodynamic center of the airfoil. For each circle, 36 acoustic monitoring points were distributed evenly with an interval of 10°. The direction of trailing edge of the airfoil was defined as 0°, the direction of the leading edge of the airfoil was defined as 180°, the direction of suction surface was 90°, and the direction of the pressure surface was 180°. To describe the noise spectrum of the monitoring points in a convenient way afterwards, the noise monitoring points at a 90° direction with radius r1 were named Point 1, Point 2, Point 3, Point 4, and Point 5 at the five different spanwise positions.

Arrangement of Noise Observation Location
First of all, the noise characteristics of a single blade were investigated. Figure 9 illustrates the noise observing points arrangement around a single blade. The acoustic monitoring points for the SPL computation were located at five different spanwise positions with h = 0.2 R, h = 0.4 R, h = 0.6 R, h = 0.8 R, and h = 1.0 R. For each spanwise position, there were two groups of observing points forming two circles with a radius r1 = 40.0 mm and r2 = 100.0 mm. The centers of the circles were at the aerodynamic center of the airfoil. For each circle, 36 acoustic monitoring points were distributed evenly with an interval of 10°. The direction of trailing edge of the airfoil was defined as 0°, the direction of the leading edge of the airfoil was defined as 180°, the direction of suction surface was 90°, and the direction of the pressure surface was 180°. To describe the noise spectrum of the monitoring points in a convenient way afterwards, the noise monitoring points at a 90° direction with radius r1 were named Point 1, Point 2, Point 3, Point 4, and Point 5 at the five different spanwise positions.

Arrangement of Noise Observation Location
First of all, the noise characteristics of a single blade were investigated. Figure 9 illustrates the noise observing points arrangement around a single blade. The acoustic monitoring points for the SPL computation were located at five different spanwise positions with h = 0.

Results and Discussions of Noise Characteristics
The comparisons of the predicted OASPL of the single blade for noise observing points with different radii are shown in Figure 10. It can be observed that the OASPLs were gradually increasing from the spanwise location of 0.2 R to 0.8 R, except for the observing positions near the directions of 0° and 180°, while the overwhelming majority of noise observing points of the spanwise location of 1.0R had smaller OASPLs than the spanwise location of 0.8 R. It is seen from Figure 10 that the noise observing points with r1 and r2 had similar results. Combined with the results of Figure 8, it is understandable that the hydrodynamic noise was closely related to the pressure fluctuations near the surface of the turbine.  Figure 11 gives the comparisons of the predicted OASPL of the single blade for noise observing points at different cross sections. To observe the results clearly, the directions from 0° to 180° were

Results and Discussions of Noise Characteristics
The comparisons of the predicted OASPL of the single blade for noise observing points with different radii are shown in Figure 10. It can be observed that the OASPLs were gradually increasing from the spanwise location of 0.2 R to 0.8 R, except for the observing positions near the directions of 0 • and 180 • , while the overwhelming majority of noise observing points of the spanwise location of 1.0R had smaller OASPLs than the spanwise location of 0.8 R. It is seen from Figure 10 that the noise observing points with r 1 and r 2 had similar results. Combined with the results of Figure 8, it is understandable that the hydrodynamic noise was closely related to the pressure fluctuations near the surface of the turbine.

Results and Discussions of Noise Characteristics
The comparisons of the predicted OASPL of the single blade for noise observing points with different radii are shown in Figure 10. It can be observed that the OASPLs were gradually increasing from the spanwise location of 0.2 R to 0.8 R, except for the observing positions near the directions of 0° and 180°, while the overwhelming majority of noise observing points of the spanwise location of 1.0R had smaller OASPLs than the spanwise location of 0.8 R. It is seen from Figure 10 that the noise observing points with r1 and r2 had similar results. Combined with the results of Figure 8, it is understandable that the hydrodynamic noise was closely related to the pressure fluctuations near the surface of the turbine.  Figure 11 gives the comparisons of the predicted OASPL of the single blade for noise observing points at different cross sections. To observe the results clearly, the directions from 0° to 180° were  Figure 11 gives the comparisons of the predicted OASPL of the single blade for noise observing points at different cross sections. To observe the results clearly, the directions from 0 • to 180 • were chosen because of the almost symmetrical display. It can be easily seen that the OASPL of the observing points with a larger radius r 2 had smaller values than that of the observing points with a smaller radius r 1 , except for the direction of 0 • and 180 • . Overall, the directionality of noise with different observing distances was identical. chosen because of the almost symmetrical display. It can be easily seen that the OASPL of the observing points with a larger radius r2 had smaller values than that of the observing points with a smaller radius r1, except for the direction of 0° and 180°. Overall, the directionality of noise with different observing distances was identical. Correspondingly, the predicted SPL spectra were analyzed. The predicted SPL spectra of Point A is shown in Figure 12. It can be easily seen that the high-value peaks were concentrated on the low frequency range, and there was broadband noise that existed after approximately 100 Hz. To exhibit the low-frequency peaks clearly, the frequency-axis was displayed in logarithmic form. It can be seen that the frequencies of the peaks were 2.0 Hz, 4.0 Hz, 6.0 Hz, 8.0 Hz, and so on, which can be calculated through the following expression:  60 knZ f (25) where n is the speed of rotation (r/min); Z is the number of blades; k is positive integer (f represents the blade passing frequency (BPF) for k = 1, and f represents the higher order harmonics for k =2, 3, 4…). Thereafter, the following frequency-axis of the noise spectra as a function of frequency were all expressed in logarithmic form. Correspondingly, the predicted SPL spectra were analyzed. The predicted SPL spectra of Point A is shown in Figure 12. It can be easily seen that the high-value peaks were concentrated on the low frequency range, and there was broadband noise that existed after approximately 100 Hz. To exhibit the low-frequency peaks clearly, the frequency-axis was displayed in logarithmic form. It can be seen that the frequencies of the peaks were 2.0 Hz, 4.0 Hz, 6.0 Hz, 8.0 Hz, and so on, which can be calculated through the following expression: where n is the speed of rotation (r/min); Z is the number of blades; k is positive integer (f represents the blade passing frequency (BPF) for k = 1, and f represents the higher order harmonics for k =2, 3, 4 . . . ). Thereafter, the following frequency-axis of the noise spectra as a function of frequency were all expressed in logarithmic form. The comparisons of the predicted SPL of the single blade for noise observing points in different directions at the spanwise position of 0.6R are shown in Figure 13. It can be seen that the sound spectra were slightly different among the different directions. For the directions of 0° and 180°, there were more peaks compared with that of 90° and 270°. However, these different directions followed the same BPF and the same higher order harmonics. Among these four observing points, the OASPL of Point A was the largest, while the OASPL of Point D was the smallest. What is more, the comparisons of the predicted SPL of the single blade for noise observing points at different cross sections are plotted in Figure 14. It is noticeable that the observing position along the spanwise direction affected the noise characteristics a lot. At a 20% span, there were only two peaks at the frequency of the BPF and two times that of the BPF. With the increase in the spanwise length, the peaks became more and more abundant. As for the OASPL, it gradually increased with the increase in the spanwise length from Point 1 to Point 4; however, the OASPL of Point 5 was smaller than that of Point 4. The comparisons of the predicted SPL of the single blade for noise observing points in different directions at the spanwise position of 0.6R are shown in Figure 13. It can be seen that the sound spectra were slightly different among the different directions. For the directions of 0 • and 180 • , there were more peaks compared with that of 90 • and 270 • . However, these different directions followed the same BPF and the same higher order harmonics. Among these four observing points, the OASPL of Point A was the largest, while the OASPL of Point D was the smallest. What is more, the comparisons of the predicted SPL of the single blade for noise observing points at different cross sections are plotted in Figure 14. It is noticeable that the observing position along the spanwise direction affected the noise characteristics a lot. At a 20% span, there were only two peaks at the frequency of the BPF and two times that of the BPF. With the increase in the spanwise length, the peaks became more and more abundant. As for the OASPL, it gradually increased with the increase in the spanwise length from Point 1 to Point 4; however, the OASPL of Point 5 was smaller than that of Point 4.
To explore the generation of the hydrodynamic noise essentially, the relationship between the flow field and the hydrodynamic noise field was further discussed. The hydrodynamic noise emission characteristics of a single blade were closely related to the pressure distributions near the surface of the turbine. It is obvious that the region near the blade had a higher OASPL value, and the value gradually increased from the blade root to the blade tip region. The maximal OASPL value was observed at the location between the blade root and the blade tip, which was closer to the blade tip region. From the point with a maximal OASPL value to the blade tip, the angle of the attack of the airfoil was extremely small, which resulted in the scatter and attenuation of the barotropic region and the negative pressure region. Thus, the OASPL at the 100% span of the blade become smaller. The noise spectra at the designated noise observing points showed that the maximum noise intensity can be obtained at 2.0 Hz, 4.0 Hz, and 6.0 Hz. The noise peak at the BPF was high, and peaks at the secondary and third frequencies were also high, which indicated that the hydrodynamic noise of the HAHT was mainly rotational noise.  To explore the generation of the hydrodynamic noise essentially, the relationship between the flow field and the hydrodynamic noise field was further discussed. The hydrodynamic noise emission characteristics of a single blade were closely related to the pressure distributions near the surface of  To explore the generation of the hydrodynamic noise essentially, the relationship between the flow field and the hydrodynamic noise field was further discussed. The hydrodynamic noise emission characteristics of a single blade were closely related to the pressure distributions near the surface of the turbine. It is obvious that the region near the blade had a higher OASPL value, and the value

Arrangement of Noise Observation Location
The investigation of the noise characteristics around a single blade could provide a closer look into the noise generation near the surface of the blade. To investigate the noise emission characteristics of the HAHT designed for underwater mooring platforms, a group of noise observing points forming a spherical surface with radius of 6500.0 mm was chosen. Figure 15 shows the observing points' arrangement around the whole turbine. The rotating center of the turbine was taken as the center of the sphere. The z-axis passed through the hydrodynamic centerline of one blade, and the y-axis, the flow direction, was the centerline of the UUMPs. There were 24 observing points evenly distributed along each circle, with an interval of 15 • . the negative pressure region. Thus, the OASPL at the 100% span of the blade become smaller. The noise spectra at the designated noise observing points showed that the maximum noise intensity can be obtained at 2.0 Hz, 4.0 Hz, and 6.0 Hz. The noise peak at the BPF was high, and peaks at the secondary and third frequencies were also high, which indicated that the hydrodynamic noise of the HAHT was mainly rotational noise.

Arrangement of Noise Observation Location
The investigation of the noise characteristics around a single blade could provide a closer look into the noise generation near the surface of the blade. To investigate the noise emission characteristics of the HAHT designed for underwater mooring platforms, a group of noise observing points forming a spherical surface with radius of 6500.0 mm was chosen. Figure 15 shows the observing points' arrangement around the whole turbine. The rotating center of the turbine was taken as the center of the sphere. The z-axis passed through the hydrodynamic centerline of one blade, and the y-axis, the flow direction, was the centerline of the UUMPs. There were 24 observing points evenly distributed along each circle, with an interval of 15°.

Results and Discussions of Noise Characteristics
As shown in Figure 16, p is an arbitrary observing point in the spherical surface with rotational center o of the three blades as the center of sphere, and line pp1 is perpendicular to plane xoy at point p1.
The angle between line op1 and the x-axis is defined as  1 , and the angle between radius op and z-axis is defined as  2 . For each observing point in the spherical surface, its position could be determined by the coordinate (  1 , 2 ), where  1 = 0°, 15°, 30°, …, 345°, and  2 = 15°, 30°, 45°,…, 165°. After the FFT analysis, the OASPL of point p could be obtained. For the noise observing points in the spherical

Results and Discussions of Noise Characteristics
As shown in Figure 16, p is an arbitrary observing point in the spherical surface with rotational center o of the three blades as the center of sphere, and line pp 1 is perpendicular to plane xoy at point p 1 . The angle between line op 1 and the x-axis is defined as α 1 , and the angle between radius op and z-axis is defined as α 2 . For each observing point in the spherical surface, its position could be determined by the coordinate (α 1 ,α 2 ), where α 1 = 0 • , 15 • , 30 • , . . . , 345 • , and α 2 = 15 • , 30 • , 45 • , . . . , 165 • . After the FFT analysis, the OASPL of point p could be obtained. For the noise observing points in the spherical surface, the length of op was equal to the radius of the sphere. To investigate the directionality of the hydrodynamic noise produced by the HAHT, the length of op had to be proportionally transformed according to the value of OASPL with the direction (α 1 ,α 2 ) unchanged. The comparison of the predicted OASPL of the whole turbine for noise observing points in a spherical surface is shown in Figure 17.
surface, the length of op was equal to the radius of the sphere. To investigate the directionality of the hydrodynamic noise produced by the HAHT, the length of op had to be proportionally transformed according to the value of OASPL with the direction (  1 , 2 ) unchanged. The comparison of the predicted OASPL of the whole turbine for noise observing points in a spherical surface is shown in Figure 17.  It is seen from Figure 17 that there were obvious directionalities for noise transmission. It is noticeable that the dipole noise source dominated the sound sources in the current model. The noise surface, the length of op was equal to the radius of the sphere. To investigate the directionality of the hydrodynamic noise produced by the HAHT, the length of op had to be proportionally transformed according to the value of OASPL with the direction (  1 , 2 ) unchanged. The comparison of the predicted OASPL of the whole turbine for noise observing points in a spherical surface is shown in Figure 17.  It is seen from Figure 17 that there were obvious directionalities for noise transmission. It is noticeable that the dipole noise source dominated the sound sources in the current model. The noise It is seen from Figure 17 that there were obvious directionalities for noise transmission. It is noticeable that the dipole noise source dominated the sound sources in the current model. The noise was more likely to transmit along the flow direction (y-axis) than the x-axis direction and z-axis direction. To investigate the noise transmitted along the six directions of the three coordinate axes, including the x-axis (+), x-axis (−), y-axis (+), y-axis (−), z-axis (+), and z-axis (−), the comparison of the predicted OASPL is shown in Figure 18. It can be seen that the noise reduction rate along the y-axis was smaller than the noise reduction rate of the x-axis and z-axis. The comparisons of the cumulative sum of the sound power spectrum E(f ) of the noise observing points along the x-axis and y-axis are exhibited in Figure 19, where f represents the frequency in the considered frequency range here. According to Wang [34], the cumulative sum of the sound power spectrum E(f ) helps to identify major noise energy distributions in a way that is better than the more well-known SPL curve versus frequency, by observing the slope dE(f )/df. The cumulative sum of the sound power spectrum is defined as the integration from f = f 1 to f = f 2 , expressed as follows: According to Wang [34], the cumulative sum of the sound power spectrum E(f) helps to identify major noise energy distributions in a way that is better than the more well-known SPL curve versus frequency, by observing the slope dE(f)/df. The cumulative sum of the sound power spectrum is defined as the integration from f= f1 to f = f2, expressed as follows: It can be seen from Figure 19 that the major noise energies were all concentrated in a low frequency range (within 10.0 Hz). With the observing distance from 1000.0 mm to 6000.0 mm, the range of E(f) along the y-axis was from 1e-4 to 1; however, the smallest E(f) along the x-axis was lower than 1e-7, and its change range was wider. Therefore, the noise reduction degree along the x-axis was larger in the observing distance from 1000.0 mm to 6000.0 mm, which shows that the directionality of the noise emission of the HAHT was obvious. For researchers dealing with the noise study of the HAHT, the flow direction (y-axis) had to be emphasized because of the larger noise pressure levels compared with the other directions.  It can be seen from Figure 19 that the major noise energies were all concentrated in a low frequency range (within 10.0 Hz). With the observing distance from 1000.0 mm to 6000.0 mm, the range of E(f ) along the y-axis was from 1e-4 to 1; however, the smallest E(f ) along the x-axis was lower than 1e-7, and its change range was wider. Therefore, the noise reduction degree along the x-axis was larger in the observing distance from 1000.0 mm to 6000.0 mm, which shows that the directionality of the noise emission of the HAHT was obvious. For researchers dealing with the noise study of the HAHT, the flow direction (y-axis) had to be emphasized because of the larger noise pressure levels compared with the other directions.
For the noise emissions around the whole turbine, the dipole noise source dominated the noise sources for the noise observing points, forming a circle parallel to the xoy surface. For each circle, the maximum OASPL could be obtained in the y-axis direction. The noise emission pattern of the HAHT investigated in this paper had similar characteristics to the results of Wasala et al. [35], concentrating on horizontal axis wind turbines, which is one typical kind of rotating machinery in the region of renewable energy. It can be seen that the OASPL did not change linearly with the observing distance, which implied that doubling the noise observing distance did not result in halving the OASPL. The cumulative sum of the sound power spectra showed that the noise energy was concentrated at a low frequency range, which is in accordance with the noise SPL spectra with main noise peaks at BPF and the first several higher order harmonics. For the noise emissions around the whole turbine, the dipole noise source dominated the noise sources for the noise observing points, forming a circle parallel to the xoy surface. For each circle, the maximum OASPL could be obtained in the y-axis direction. The noise emission pattern of the HAHT investigated in this paper had similar characteristics to the results of Wasala et al. [35], concentrating on horizontal axis wind turbines, which is one typical kind of rotating machinery in the region of renewable energy. It can be seen that the OASPL did not change linearly with the observing distance, which implied that doubling the noise observing distance did not result in halving the OASPL. The cumulative sum of the sound power spectra showed that the noise energy was concentrated at a low frequency range, which is in accordance with the noise SPL spectra with main noise peaks at BPF and the first several higher order harmonics.

Conclusions
As a renewable energy source, marine current energy has been paid increasing attention all over the world. However, the noise emission problem of the energy conversion devices cannot be ignored, and the hydrodynamic noise generated by the rotating turbine has been known to be the most important noise source. In this study, the noise characteristics of the HAHT designed for UUMPs are evaluated based on LES, combined with the FW-H equation. In general, the coefficient of power of the hydrokinetic turbine after the LES calculation shows good agreement with the experimental data at different TSRs, denoting that the present flow field model is capable of accurately simulating the flow field characteristics. After verifying and validating the numerical method, the flow field characteristics and the noise characteristics are investigated. The main conclusions from the perspective of flow field characteristics, single blade noise characteristics, and whole turbine noise characteristics are as follows: Figure 19. Comparison of the cumulative sum of the sound power spectrum E(f ) of the noise observing points along x-axis and y-axis.

Conclusions
As a renewable energy source, marine current energy has been paid increasing attention all over the world. However, the noise emission problem of the energy conversion devices cannot be ignored, and the hydrodynamic noise generated by the rotating turbine has been known to be the most important noise source. In this study, the noise characteristics of the HAHT designed for UUMPs are evaluated based on LES, combined with the FW-H equation. In general, the coefficient of power of the hydrokinetic turbine after the LES calculation shows good agreement with the experimental data at different TSRs, denoting that the present flow field model is capable of accurately simulating the flow field characteristics. After verifying and validating the numerical method, the flow field characteristics and the noise characteristics are investigated. The main conclusions from the perspective of flow field characteristics, single blade noise characteristics, and whole turbine noise characteristics are as follows: (1) Along the spanwise direction of the blade, the pressure fluctuations are gradually becoming serious, with the increase of spanwise length from 20% span to 80% span. What is more, it is noticeable that the pressure fluctuations are concentrated in the region of the leading edge of the airfoil and the suction surface of the airfoil. With the increase in spanwise length, the negative pressure region gradually moves from the leading edge to the trailing edge on the suction pressure surface. (2) From the OASPL directionality around a single blade, it can be seen that the overall noise level of the spanwise location of 0.8 R has largest OASPL than that of the other spanwise locations in general, which is in accordance with the results of the flow field characteristics near the blade region. The sound spectra are slightly different among the different observing directions at the same spanwise location, and the observing location along the spanwise direction affects the noise spectra characteristics to a very great extent. (3) The investigation of the noise around the whole turbine indicates that the dipole noise source dominates the sound sources in the current analysis. The noise is more likely to transmit along the flow direction (y-axis in this paper) than the other directions. Furthermore, the cumulative sum of the sound power spectrum of the noise observing points along the x-axis and y-axis show the change relationship clearly.
In the present research, the HAHT is modeled as rotating about the fixed axis, and the interaction movement between the UUMPs is ignored. Therefore, further studies still need to be carried out on the following issues: First of all, the effect of the interaction movement of the UUMPs on noise performance should be carried out. In actual working conditions, the UUMPs will move back and forth with the current; thus, more complex working environment and boundary conditions of the UUMPs and HAHT have to be conducted.
Secondly, the noise reduction design will be carried out according to the noise results around a single blade. Among the present noise reduction technologies of the turbines, the most suitable technology has to be determined with the modification of the blade.
Additionally, the experimental test of the noise performance of the HAHT designed for the UUMPs has to be carried out.