Numerical Simulation of the Hydrodynamic Performance and Self-Propulsion of a UUV near the Seabed

: Unmanned underwater vehicles (UUV) face maneuverability and rapidity challenges when they are applied for detecting and repairing submarine oil and gas pipelines, and ﬁber cables near the seabed. This research establishes numerical models of the bare UUV and self-propelled UUV near the seabed using the computational ﬂuid dynamics (CFD) method. The effect of dimensionless distance Hd and Re L on the hydrodynamic performance of the vehicle and the interaction between the hull and the propeller is investigated. The range of Hd is 1.5 D –10 D , and the Re L is 9.97 × 10 5 ~7.98 × 10 6 . Findings indicate that: (1) There is an obvious strong coupling between the hydrodynamic performance of the bare UUV and Hd . With the increase of Hd , the hydrodynamic performance such as C d , the absolute value of C l and m y decreases continuously and ﬁnally tends to be stable. The absolute values of C d and C l increase with the increase of Re L . The change trend of m y is opposite to that of C l . (2) The variation trend of hydrodynamic performance of the self-propelled UUV with Hd is consistent with those of the bare UUV. Additionally, it increases to some extent, respectively, compared with the bare UUV. (3) The self-propelled characteristics such as t , η H , w and η i are weakly related to Hd . The t and η i increase with the increasing of Re L , while η H and w decrease with the increasing of Re L .


Introduction
UUV is a kind of underwater vehicle, which integrates technologies of underwater communication, automatic control, mode recognition, and artificial intelligence.It is capable of performing underwater tasks autonomously without any operator and is important equipment for exploring the ocean.The UUV has been widely used in marine environment exploration, undersea oil and gas pipeline and cable inspection, hydrological research, etc.The operability and rapidity of UUV are receiving more and more attention, which reflects the safe operations and navigation efficiency of the vehicles.They play an important role in traveling distance and duration of travel.The operability is mainly represented by the hydrodynamic characteristics of the vehicle, and the rapidity generally needs to consider interferences between hull and thruster.In recent years, a great number of studies have studied the hydrodynamic characteristics of UUV and the interaction between hull and thruster through experiments and CFD methods.
The drag and self-propulsion analysis of SUBOFF with E1619 propeller was carried out using CFD method by Sezen et al. [1,2].The propellers of self-propelled submarines were modeled using the actuated disc and real propellers, respectively.The numerical results are compared with the available experimental data.The suitability of CFD methods for the prediction of self-propulsion performance of underwater vehicles is discussed.Posa et al. [3] reported the results of self-propelled SUBOFF and compared with the towing test results.They find that the bimodal distribution of turbulent stresses in the towing case is different with the axial peak of turbulent kinetic energy in self-propulsion due to the backpressure gradient induced by the stern section.Carria et al. [4] studied the effect of mesh refinement on propeller performance and wake flow for E1619 propeller and SUBOFF by CFD method.The results showed the applicability of different turbulence models for propeller wake flow and blade tip vortex analysis, and analyzed the self-propulsion performance of SUBOFF.
The near-surface and subsea conditions are two extreme conditions for UUV, which put the maneuverability of UUV to severe challenges.The UUV need to navigate near the surface in order to share, transmit information, obtain energy and update their position.In the near-subsea condition, the UUV detects submarine oil and gas pipelines and fiber optic cables.These two conditions need to consider the nonlinear effect of waves on the vehicles near the surface and the blocking effect of the seabed.
Studies on the interaction between the hull and the propeller near the surface and fully submerged were carried out in a towed tank [5].The experimental results show that the drag coefficient of UUV in the surface condition is larger than that in the submerged condition.Besides, the reliability of CFD method is verified by the experimental results.It provides a new method to study the interference between hull and propeller.Tian et al. [6] used the two-phase flow calculation CFD method to study the performance of UUV under the wave-current coupling.It's found that the lift of the UUV is greatly influenced by waves, and the average drag increases with the increase of wave height and the decrease of submerged depth.Carria [7] studied the hydrodynamic performance of the SUBOFF with an E1658 propeller near the water surface by experiments and CFD overlapping mesh method.The interaction between the vessel and free surface strongly influences the inflow and wake of the propeller, generating high local propulsion coefficients and blade loads near the propeller surface.Strong surface fluctuations at small axial depths lead to instability and collapse of propeller tip vortices in the near field.The performance of the UUV with a propeller near water surface was simulated by Zhang et al. [8] using the slip grid technique and the VOF method.Results show that the effect of the free surface on the drag force is greater than that of the self-propulsion.
Hydrodynamic characteristics and self-propulsion performance of self-propelled of UUV near the water surface have been carried out by a large number of researchers using CFD and experimental methods.However, there are few papers studying the nearsubsea conditions to the authors' knowledge.Song et al. [9][10][11] proposed two subsea stationing strategies, hydraulically supported and mooring chain types, for UUV that require permanent fixed-point exploration.Two simulation models of the landing process of UUV were established.The results show that the lateral force of UUV with vertical thrust is 5~15% larger than that without vertical thrust.Zhang et al. [12] developed a mathematical model of critical instability considering the complex coupling between the glider and the seabed as well as the flow field in which the underwater glider resides stably.The results show that the center of buoyancy is 0.3 m from the seabed and the angle of attack of the glider has a strong influence on the stability of the system.
However, the UUV analyzed in these papers are stationary without motion.Zhu et al. [13] conducted a model resistance test for UUV at different distances near the seabed and near the sea surface, and established a model test and analysis method to obtain the influence rules of UUV drag near the seabed and near the sea surface.The results demonstrate that considering only the effect of the free surface, it can be assumed that there is no wave resistance when the vessel travels at a depth of over 1/3 of its length.Taking into account only the influence of the seabed wall, the effect of the seabed is no longer considered when the submerged body travels at a distance of over 1/4 of its length from the seabed.When the influences of both free water surface and sea bed wall exist, to completely avoid the interference of restricted water, it is required that submerged hull dive depth exceeds 1/3 of the hull length, and the distance between the seabed and the axis of the submerged hull should not be less than 1/2 of the hull length.
Du et al. [14] developed a numerical model of the near-subsea UUV maneuverability and investigated the effects of distance from subsea, Re and angle of attack on the hydrody-namic performance of the UUV.The results show that the drag coefficient of UUV increases with decreasing distance from the seabed and with increasing angle of attack.The absolute value of the lift coefficient increases with decreasing distance and with the increment of the angle of attack.The absolute value of the pitching moment coefficient increases as the distance decreases and the angle of attack increases.Yan et al. [15] established a numerical model of a near-subsea UUV in order to obtain the effects of subsea boundary effects and complex terrain on the operational stability and maneuverability.The effects of subsea distance, Re, angle of attack, and inclined terrain on the hydrodynamic performance of the UUV were analyzed.
In summary, there is little experimental and CFD literature on UUV near the seabed.The studies focus on the hydrodynamic characteristics of UUV, while those of self-propelled UUV are still relatively few.Therefore, it is meaningful to investigate the hydrodynamic performance and self-propulsion performance of the UUV near the seabed through CFD method in this paper.Two models, bare UUV and self-propelled UUV near the seabed, were established and analyzed the effects of Hd and Re L on the performance of vessel.

UUV Geometry Parameters
The self-propelled UUV used in this research includes three parts: the main body of the vessel (hull), the fin plates and the propellers, as shown in Figure 1.The hull profile is generated using the theoretical equation proposed in reference [6].Its total length is 2 m and the maximum diameter is 0.2 m.It can be divided into three parts: nose, parallel section and tail end, with the lengths of 0.3 m, 1.2 m and 0.5 m, respectively.The main parameters of the UUV hull are shown in Table 1.To match the hull, a scaled-down 0.6× DTMB4119 model propeller is used for the propulsion section, with a diameter of 0.1829 m.Detailed parameters are shown in Table 2.The hub of the propeller is suitably altered to match the stern of the hull.
Du et al. [14] developed a numerical model of the near-subsea UUV maneuverability and investigated the effects of distance from subsea, Re and angle of attack on the hydrodynamic performance of the UUV.The results show that the drag coefficient of UUV increases with decreasing distance from the seabed and with increasing angle of attack.The absolute value of the lift coefficient increases with decreasing distance and with the increment of the angle of attack.The absolute value of the pitching moment coefficient increases as the distance decreases and the angle of attack increases.Yan et al. [15] established a numerical model of a near-subsea UUV in order to obtain the effects of subsea boundary effects and complex terrain on the operational stability and maneuverability.The effects of subsea distance, Re, angle of attack, and inclined terrain on the hydrodynamic performance of the UUV were analyzed.
In summary, there is little experimental and CFD literature on UUV near the seabed.The studies focus on the hydrodynamic characteristics of UUV, while those of self-propelled UUV are still relatively few.Therefore, it is meaningful to investigate the hydrodynamic performance and self-propulsion performance of the UUV near the seabed through CFD method in this paper.Two models, bare UUV and self-propelled UUV near the seabed, were established and analyzed the effects of Hd and on the performance of vessel.

UUV Geometry Parameters
The self-propelled UUV used in this research includes three parts: the main body of the vessel (hull), the fin plates and the propellers, as shown in Figure 1.The hull profile is generated using the theoretical equation proposed in reference [6].Its total length is 2 m and the maximum diameter is 0.2 m.It can be divided into three parts: nose, parallel section and tail end, with the lengths of 0.3 m, 1.2 m and 0.5 m, respectively.The main parameters of the UUV hull are shown in Table 1.To match the hull, a scaled-down 0.6× DTMB4119 model propeller is used for the propulsion section, with a diameter of 0.1829 m.Detailed parameters are shown in Table 2.The hub of the propeller is suitably altered to match the stern of the hull.

Definition of Dimensionless Parameters
In order to classify the computational cases and analyze the effects of Hd and Re L on the rapidity and maneuverability of the UUV more conveniently, the dimensionless parameter Hd, which is the ratio of the distance between the seabed and the axis of UUV to the maximum diameter of the UUV, is introduced in this paper.Additionally, the values of Hd are 1.5, 2, 2.5, 3, 5, 7, 9 and 10, which cover the extreme conditions of subsea exploration in practices.The oncoming velocity takes values in the range of 0.5 m/s, 1 m/s, 2 m/s, and 4 m/s, and the corresponding Re L ranges from 9.97 × 10 5 to 7.98 × 10 6 depending on the length dimension of the vessel.The 3D dimensionless coefficients of drag, lift and pitch moment coefficients are used to describe the hydrodynamic characteristics of the UUV, respectively, which are calculated using following formula. (2) C d , C l and m y describe the hydrodynamic characteristics of the UUV.The formulas F x , F z and M y the drag, lift and pitch moment of the UUV.S A is the maximum cross-sectional area of the UUV.U is the speed of the vehicle.L A is the total length of the UUV.ρ is the fluid density taken as 998 kg/m 3 .
J 0 , K T and K Q are the propeller open water performance parameters.J 0 is an important parameter for propeller's performance, which is similar to the angle of attack in wing theory.V A is the velocity at the propeller disc; n is the rotational speed of propeller in rps; D is the propeller diameter.K T and K Q are the propeller thrust and torque coefficients, with T and Q being the propeller thrust and torque.
The dimensionless parameters t, w, η i and η H depict the interaction between the hull and propeller.t describes the effects of the propeller on the drag characteristics of the hull.R is the drag of the self-propelled UUV.R 0 is the drag of the bare UUV at the corresponding speed.Additionally, T is the thrust provided by the propeller on the aft of the UUV.In the absence of scaling effects, R and T are equal.w and η i define the effects of the wake field of the self-propelled UUV on the propeller performance.J and K Q are the feed ratios and Appl.Sci.2022, 12, 6975 5 of 21 torque coefficients of the propeller of the self-propelled UUV.J Q and K Q 0 are the feed ratios and torque coefficients of the propeller in open water condition derived by using the equal thrust method.η H describes the hull efficiency and represents the combined effect of thrust reduction and wake fraction.
The calculation cases in this paper are classified into two main categories, as shown in Table 3.

Numerical Methods
To evaluate the hydrodynamic performance of the bare UUV and the self-propelled UUV, as well as the self-propulsion performance near subsea, the ANSYS FLUENT 17.0 was used to solve the governing equations.The Moving Reference Frame (MRF) model and SST − kw turbulence model are used.The governing equations are discretized by the finitevolume method, where convective and diffusive terms in the equations are discretized with second-order-upwind differential format, central differential format, respectively.To speed up the convergence, the Couple method is used to solve the discrete equations.

Governing Equations
The Reynolds average equation includes: the continuum equation, and the momentum equation, as shown in Equations ( 12) and ( 13).
∂ρ ∂t u i is the fluid velocity component.P is the total pressure.F i is the volume force term.Additionally, µ is the dynamic viscosity coefficient.The SST − kw turbulence model is used to simulate the turbulence terms in the RANS equations.This model is able to simulate turbulent shear stress transport and accurately predict the flow separation under the inverse pressure gradient.Turbulence kinetic energy and turbulence dissipation rate are calculated by the following equations.
To investigate the effect of Hd and Re L on the hydrodynamic performance of the selfpropelled UUV and the coupling effect of the hull and the propeller, firstly, it is necessary to obtain the rotational speed of the propeller under different operating conditions, and then solve the hydrodynamic performance and self-propulsion factor of the UUV.The rotational speed of propeller is corrected by the DEFINE_ZONE_MOTION macro once every 1000 iterations.The rotational speed for self-propelled point is considered to be obtained when the rotational speed value of the previous step does not change more than 0.01% with respect to the previous rotational speed correction value and the relative error between the drag of the UUV and the thrust of the propeller does not exceed 0.01%.The propeller rotational speed is calculated according to the following equation.
N is the currently rotational speed.F uuv is the drag of the UUV without propeller.F blade is the thrust of the propeller.δ n is the convergence coefficient.

Computational Domain and Boundary Conditions
The computational domain contains two parts: the stationary domain and the rotational domain.The cross section of the stationary domain is a square with a length and width of 20 D A , which is the diameter of the maximum cross section of the UUV.The length of the flow field domain is 10 L A , which is the length of UUV to ensure the full development of both the upstream incoming flow and wake flow, and the nose is 4 L A from the inlet boundary.The rotating domain is a cylinder with a diameter of 200 mm and a length of 60 mm.The detailed sizes of the domain and the boundary condition settings are shown in Figure 2. Information about the flow field is transferred between the stationary and rotating domains through an interface.By reference to the findings of other researchers, the size of this flow field domain satisfies the requirements of the calculation.Additionally, the results of the validation cases that adopt the same strategy show that the size of the computational domain is reasonable.
The propeller rotational speed is calculated according to the following equation.
, is the currently rotational speed. is the drag of the UUV without propeller.is the thrust of the propeller.
is the convergence coefficient.

Computational Domain and Boundary Conditions
The computational domain contains two parts: the stationary domain and the rotational domain.The cross section of the stationary domain is a square with a length and width of 20 DA, which is the diameter of the maximum cross section of the UUV.The length of the flow field domain is 10 LA, which is the length of UUV to ensure the full development of both the upstream incoming flow and wake flow, and the nose is 4 LA from the inlet boundary.The rotating domain is a cylinder with a diameter of 200 mm and a length of 60 mm.The detailed sizes of the domain and the boundary condition settings are shown in Figure 2. Information about the flow field is transferred between the stationary and rotating domains through an interface.By reference to the findings of other researchers, the size of this flow field domain satisfies the requirements of the calculation.Additionally, the results of the validation cases that adopt the same strategy show that the size of the computational domain is reasonable.

Meshing Settings
Ansys fluent mesh is applied to discrete the computational domain.The stationary domain is a mixed mesh containing polyhedral mesh and hexahedral mesh, and the rotating domain is a polyhedral mesh, as shown in Figure 3.In order to effectively and accurately predict the hydrodynamic properties of the UUV and propeller, two layers of encrypted regions are set up around the UUV and propeller, as shown in Figure 3a. Figure 3b shows the nose grid and Figure 3c shows the stern and rotation domain grid.To capture the shape of the propeller, the grid size of the leading and trailing edges of the propeller is 0.05 mm. Figure 3d shows the grid of the propeller blade surface.Prismatic boundary layer meshes are set for all the non-slip boundaries in the flow field domain, and the height Y+ of the mesh is controlled to be around 1 to meet the needs of the turbulence model.A total of 15 boundary layers are set up.In order to capture the effects of the seabed on the hydrodynamics of the vehicle, the bottom wall of the stationary domain is arranged with 100 nodes in the length of 1L A .
tating domain is a polyhedral mesh, as shown in Figure 3.In order to effectively and accurately predict the hydrodynamic properties of the UUV and propeller, two layers of encrypted regions are set up around the UUV and propeller, as shown in Figure 3a. Figure 3b shows the nose grid and Figure 3c shows the stern and rotation domain grid.To capture the shape of the propeller, the grid size of the leading and trailing edges of the propeller is 0.05 mm. Figure 3d shows the grid of the propeller blade surface.Prismatic boundary layer meshes are set for all the non-slip boundaries in the flow field domain, and the height Y+ of the mesh is controlled to be around 1 to meet the needs of the turbulence model.A total of 15 boundary layers are set up.In order to capture the effects of the seabed on the hydrodynamics of the vehicle, the bottom wall of the stationary domain is arranged with 100 nodes in the length of 1 .A grid independent validation is performed in this paper to confirm that the numerical results do not depend on the quantity of grids.Three grid resolutions, 3.2 M (coarse grid), 6 M (medium grid) and 8.4 M (fine grid), are set for the self-propelled UUV case with Hd of 2 and oncoming velocity of 1 m/s.Table 4 summarizes the numerical results for different grid resolutions.It can be seen that the numerical results of the fine and medium meshes are close to each other, which predicts that further increment of the number of mesh will not significantly affects the numerical results of CFD.Therefore, the medium resolution grid is used to complete the subsequent numerical simulations.In order to eliminate the influence of grid resolution on the numerical results of different simulation cases, the medium grid resolution is used for the bare UUV simulation cases and the selfpropelled UUV cases.A grid independent validation is performed in this paper to confirm that the numerical results do not depend on the quantity of grids.Three grid resolutions, 3.2 M (coarse grid), 6 M (medium grid) and 8.4 M (fine grid), are set for the self-propelled UUV case with Hd of 2 and oncoming velocity of 1 m/s.Table 4 summarizes the numerical results for different grid resolutions.It can be seen that the numerical results of the fine and medium meshes are close to each other, which predicts that further increment of the number of mesh will not significantly affects the numerical results of CFD.Therefore, the medium resolution grid is used to complete the subsequent numerical simulations.In order to eliminate the influence of grid resolution on the numerical results of different simulation cases, the medium grid resolution is used for the bare UUV simulation cases and the self-propelled UUV cases.

Verification of Numerical Method
Due to the lack of experimental data for the self-propelled model, the articles adopt the same validation strategy as the others [1][2][3][4].Numerical simulations of the UUV and propeller from the available experimental study data are carried out separately to verify the feasibility of the numerical method used in the paper.
The SUBOFF is an axisymmetric vehicle and the hydrodynamic test is carried out by David Taylor Research Center.It is aimed at providing a data base for CFD validation and other flow field analysis related to submarines.A large number of academics adopt SUBOFF as a research object for UUV.The AFF-3 model in SUBOFF is chosen in this paper to verify the feasibility of the numerical method.Huang et al. [16] experimentally studied the hydrodynamic performance of the AFF-8 model at a Re L of 12 × 10 6 .The wake characteristics of the AFF-2 model at Re L 1.1 × 10 6 were investigated by Jiménez et al. [17] through wind tunnel tests.Posa et al. [18] investigated the turbulent boundary layer on the surface of the AFF-8 model with a Re of 1.2 × 10 6 by CFD method.
The distribution of pressure coefficient (C p ) with Re L of 1.2 × 10 6 is obtained by CFD method.C p is calculated by Equation (17).The comparison of the numerical results is given in Figure 4.It can be seen that the distribution of C p on the hull surface agrees well with the results of other researchers, except that field near the fin plate of the model.The main reason for the phenomenon is the absence of a bridge in the AFF-3 model.The differences in the geometric models lead to the divergence in the flow field distribution at the leading edge of the stern appendages.
the surface of the AFF-8 model with a Re of 1.2 × 10 6 by CFD method.
The distribution of pressure coefficient ( ) with of 1.2 × 10 6 is ob method.
is calculated by Equation (17).The comparison of the nume given in Figure 4.It can be seen that the distribution of on the hull surfa with the results of other researchers, except that field near the fin plate of t main reason for the phenomenon is the absence of a bridge in the AFF-3 m ferences in the geometric models lead to the divergence in the flow field the leading edge of the stern appendages.
There are two reasons for the discrepancy.
1.The geometric profile of the model used in the article is different from The geometric model in the article contains the stern appendages, bu water, while [16,18] contain both the stern appendages and the fair-w has neither the stern appendages nor the fair-water.2. The most significant reason is that the cross section of [16][17][18] from w obtained did not pass through the stern appendages, while the cross paper passed through it.There are two reasons for the discrepancy.

1.
The geometric profile of the model used in the article is different from that of [16][17][18].The geometric model in the article contains the stern appendages, but not the fairwater, while [16,18] contain both the stern appendages and the fair-water, and [17] has neither the stern appendages nor the fair-water.

2.
The most significant reason is that the cross section of [16][17][18] from which the C p is obtained did not pass through the stern appendages, while the cross section of this paper passed through it.
The DTMB 4119 propeller has been used to verify the reliability of numerical methods and has been widely applied to research the cavitation performance of the propeller.The open water characteristics of this propeller were obtained by CFD method and compared with the experimental results.As shown in Figure 5, the numerical results agree well with the experimental results.The largest discrepancy appears in the efficiency variable with J of 0.833, with a relative error of 2.57%, which can be accepted.The experimental results for this propeller were obtained from reference [19].
The open water characteristics of this propeller were obtained by CFD met pared with the experimental results.As shown in Figure 5, the numerical well with the experimental results.The largest discrepancy appears in the e able with J of 0.833, with a relative error of 2.57%, which can be accepted.The results for this propeller were obtained from reference [19].

The Hydrodynamic Characteristics of Near-Seabed bare UUV
The hydrodynamic characteristics of the bare UUV near the seabed incl and pitching moment characteristics.The variation of the hydrodynamic c of the bare UUV near the seabed with Hd and are given in Figures shows the distribution of UUV cross-section (a) for each Hd case when th m/s.Cross-section (a) is defined by Figure 1, which is the longitudinal sectio Figure 10 shows the distribution of UUV at cross-section (a) for Hd of 1 1 × 10 6 ~8 × 10 6 .The is defined by Equation (17), where is the statics p at the inlet and is the fluid velocity at the inlet.

Drag Characteristics
The drag coefficient of the bare UUV decreases continuously with the i and finally tends to be stable, and the variation rate of the coefficient is th the seabed.As shown in Figure 6, the drag coefficient of UUV is greater ne When Hd is greater than 5, the drag coefficient tends to be stable and rem This blocking effect between the UUV and the seabed causes an increase in t pressure in the X-direction on the UUV surface.Figure 9 shows that th greater than 0 until 3.5% , and decreases with increasing Hd.When Hd i 5, the is almost constant.The value near the nose of the UUV decre ously with the increase of Hd, that is, the pressure near the nose along th decreases with the increase of Hd.The pressure component in the X-directi allel section of the UUV is 0, so it is ignored.The value near 83.45% and the absolute value of decreases continuously with the increase of H

The Hydrodynamic Characteristics of Near-Seabed Bare UUV
The hydrodynamic characteristics of the bare UUV near the seabed includes drag, lift and pitching moment characteristics.The variation of the hydrodynamic characteristics of the bare UUV near the seabed with Hd and Re L are given in Figures 6-8. Figure 9 shows the C p distribution of UUV cross-section (a) for each Hd case when the speed is 0.5 m/s.Cross-section (a) is defined by Figure 1, which is the longitudinal section of the UUV. Figure 10 shows the C p distribution of UUV at cross-section (a) for Hd of 1.5 and Re L of 1 × 10 6 ~8 × 10 6 .The C p is defined by Equation (17), where P 0 is the statics pressure value at the inlet and V 0 is the fluid velocity at the inlet.The drags of the bare UUV are mainly contributed from the sum of the ues of the nose and tail.The variation of the tail pressure value is one order smaller than that of the nose pressure.When Hd is greater than 5, both ar order of magnitude.When Hd is less than 5, the drag of the UUV is mainly d the component of the pressure at the nose of the UUV along the X-directi creases, the drag of UUV decreases.When Hd is greater than 5, the drag of U constant.

•
Lifting Characteristics Absolute value of lift coefficient of bare UUV decreases continuously crease of Hd and finally tends to be zero, and the variation rate of the coe           coefficient of the UUV decreases with increasing .The values of UUV are less tha 0 behind 3.5% , and the values decrease with increasing .This explains that th absolute values of UUV lift and pitching moment increase with increasing .The variation rule of pitch moment coefficient with is not consistent with the lif coefficient because the increment of pitch moment with velocity is smaller than the sec ond-squares increment of velocity.The drag coefficient of the self-propelled UUV decreases continuously with the in crease of Hd and finally tends to be stable, and the variation rate of the coefficient is th largest near the seabed.Figure 11 shows that the drag coefficient of self-propelled UUV near the seabed is greater than that away from seabed.When Hd is greater than 5, the dra coefficient tends to be stable and remains constant.The drag coefficient of the self-pro pelled UUV increased by 16.74~17.88%compared to that of the bare UUV.The drag coefficient of the bare UUV decreases continuously with the increase of Hd and finally tends to be stable, and the variation rate of the coefficient is the largest near the seabed.As shown in Figure 6, the drag coefficient of UUV is greater near the seabed.When Hd is greater than 5, the drag coefficient tends to be stable and remains constant.This blocking effect between the UUV and the seabed causes an increase in the differential pressure in the X-direction on the UUV surface.Figure 9 shows that the C p value is greater than 0 until 3.5% L A , and decreases with increasing Hd.When Hd is greater than 5, the C p is almost constant.The value C p near the nose of the UUV decreases continuously with the increase of Hd, that is, the pressure near the nose along the X-direction decreases with the increase of Hd.The pressure component in the X-direction of the parallel section of the UUV is 0, so it is ignored.The C p value near 83.45%L A is less than 0, and the absolute value of C p decreases continuously with the increase of Hd.That is, the pressure value there increases continuously with the increase of Hd increases.
The drags of the bare UUV are mainly contributed from the sum of the pressure values of the nose and tail.The variation of the tail pressure value is one order of magnitude smaller than that of the nose pressure.When Hd is greater than 5, both are of the same order of magnitude.When Hd is less than 5, the drag of the UUV is mainly determined by the component of the pressure at the nose of the UUV along the X-direction.As Hd increases, the drag of UUV decreases.When Hd is greater than 5, the drag of UUV is almost constant.

• Lifting Characteristics
Absolute value of lift coefficient of bare UUV decreases continuously with the increase of Hd and finally tends to be zero, and the variation rate of the coefficient is the greatest near the seabed.As shown in Figure 7, the absolute value of lift coefficient of UUV near seabed is larger than that away from seabed, which is shown as the ground diving phenomenon near seabed without angle of attack.This phenomenon benefits from the asymmetric distribution of pressure on the upper and lower surfaces of the UUV, as shown in Figure 15.A further reason is due to the blocking effect between the UUV and the seabed.This leads to an increase in the fluid velocity between them and consequently a decrease in the static pressure value at that position, which eventually leads to a lift of the UUV less than 0. Figure 9 shows that the pressure at the nose and parallel segments varies with Hd.When X/L A is greater than 3.5%, the surface pressure coefficient of the UUV is less than 0. Additionally, the absolute value of the pressure coefficient decreases continuously with the increase of Hd.The absolute value of the parallel segment C p decreases continuously with the increase of Hd.When Hd is greater than 7, the absolute value of parallel segment C p is almost unaffected by the seabed.Therefore, the absolute value of UUV lift keeps decreasing with the increase of Hd, and when Hd is greater than 7, the UUV lift is almost constant.

• Pitch Moment Characteristics
The center of buoyancy of the bare UUV is treated as the origin of the pitching moment.The pitching moment coefficient of the bare UUV decreases with the increase of Hd and finally tends to be zeros.Additionally, the variation rate of the coefficient is the greatest near the seabed.As shown in Figure 8, the pitching moment coefficient of UUV near seabed is greater than that away from seabed.There is a tendency of the UUV to dive in the near-seabed case.The lift coefficients on the UUV surface are unevenly distributed along the X-direction.It can be seen from Figure 9 that the lift coefficient before the buoyancy center is larger than that after the buoyancy center.The rule of variation of the bare UUV pitching moment with Hd is consistent with that of the lift coefficient.

Correlation between the Hydrodynamic Characteristics of Bare UUV and Re L
From Figures 6-8, some remarks can be concluded.The UUV drag coefficient decreases with the increase of Re L .The absolute of lift coefficient increases with the increase of Re L .When Hd is greater than 7, the relationship between lift coefficient and pitching moment coefficient and Re L becomes weaker.
From Figure 10, it can be found that the C p near the nose of the UUV decreases with increasing Re L , while the absolute value of C p near the tail increases with increasing Re L .As a result, the differential pressure drag of the UUV decreases with increasing Re L , while the frictional drag coefficient remains constant in this Re L range.Finally, the drag coefficient of the UUV decreases with increasing Re L .The C p values of UUV are less than 0 behind 3.5% L A , and the C p values decrease with increasing Re L .This explains that the absolute values of UUV lift and pitching moment increase with increasing Re L .
The variation rule of pitch moment coefficient with Re L is not consistent with the lift coefficient because the increment of pitch moment with velocity is smaller than the second-squares increment of velocity.The drag coefficient of the self-propelled UUV decreases continuously with the increase of Hd and finally tends to be stable, and the variation rate of the coefficient is the largest near the seabed.Figure 11 shows that the drag coefficient of self-propelled UUV near the seabed is greater than that away from seabed.When Hd is greater than 5, the drag coefficient tends to be stable and remains constant.The drag coefficient of the self-propelled UUV increased by 16.74~17.88%compared to that of the bare UUV.
The absolute value of the lift coefficient of the self-propelled UUV decreases continuously with the increase of Hd and finally tends to be around 0.008, and the variation rate of the coefficient is the largest near the seabed.As can be seen from Figure 12, the absolute value of the lift coefficient of self-propelled UUV near the seabed is greater than that away from seabed.When Hd is greater than 7, the lift coefficient of UUV tends to be constant.It is consistent with the trend of the lift coefficient of bare UUV with Hd.However, the lift coefficient of the self-propelled UUV increases by 69~130% compared to that of the bare UUV when the Hd less than 3.The absolute value of the lift coefficient of the self-propelled uously with the increase of Hd and finally tends to be around 0.00 of the coefficient is the largest near the seabed.As can be seen from value of the lift coefficient of self-propelled UUV near the seabed i from seabed.When Hd is greater than 7, the lift coefficient of UUV is consistent with the trend of the lift coefficient of bare UUV wi coefficient of the self-propelled UUV increases by 69~130% comp UUV when the Hd less than 3.The absolute value of the lift coefficient of the self-propelled uously with the increase of Hd and finally tends to be around 0.008 of the coefficient is the largest near the seabed.As can be seen from value of the lift coefficient of self-propelled UUV near the seabed is from seabed.When Hd is greater than 7, the lift coefficient of UUV is consistent with the trend of the lift coefficient of bare UUV wit coefficient of the self-propelled UUV increases by 69~130% comp UUV when the Hd less than 3. Similar to the bare UUV, the self-propelled UUV takes the UUV with a propeller as the origin of the pitching moment.The p of the self-propelled UUV decreases with the increase of Hd and fin 0.005, and the variation rate of the coefficient is the largest near the from Figure 13, the pitching moment coefficient near the seabed i from the seabed, which causes the dive phenomenon near the s attack.The pitching moment curve of the self-propelled UUV v according to the same trend as that of the bare UUV.However, th efficient of the self-propelled UUV increased by 100~300% relative Similar to the bare UUV, the self-propelled UUV takes the buoyancy center of the UUV with a propeller as the origin of the pitching moment.The pitch moment coefficient of the self-propelled UUV decreases with the increase of Hd and finally tends to be around 0.005, and the variation rate of the coefficient is the largest near the seabed.As can be seen from Figure 13, the pitching moment coefficient near the seabed is larger than that away from the seabed, which causes the dive phenomenon near the seabed without angle of attack.The pitching moment curve of the self-propelled UUV varies with Hd and Re L according to the same trend as that of the bare UUV.However, the pitching moment coefficient of the self-propelled UUV increased by 100~300% relative to that of the bare UUV when the Hd less than 3.The correlation of hydrodynamics coefficients of the self-pro consistent with that of the bare UUV.The drag coefficient decrea .The absolute of lift coefficient increases with the increase of of pitch moment coefficient with is not consistent with the lif increment of pitch moment with velocity is smaller than the secon velocity.When Hd is greater than 7, the relationship between lift moment coefficient and become weaker.

Effects of Hd and on UUV Self-Propelled Performance
The relevant self-propulsion coefficients are derived from Eq alyze the effects of Hd and on the interactions between the parameters such as thrust reduction, hull efficiency, effective wake efficiency are selected for analysis.The first two of these paramet thruster on the UUV, and the last two reflect the effect of the UUV

•
The Variation of Thrust Reduction describes the effect of propeller at the wake of the hull o rotation of the thruster at the tail of the UUV, velocity of flow nea the static pressure here decreases relative to the static pressure of t in an increase in the drag of the UUV.For the increased UUV dra erally attributed to the reduced thrust of the propellor.
As can be seen from Figure 14, decreases firstly with the increases slowly and finally tends to stable.Near the seabed, the the seabed increases the pressure difference between the lead and The correlation of hydrodynamics coefficients of the self-propelled UUV and Re L is consistent with that of the bare UUV.The drag coefficient decreases with the increase of Re L .The absolute of lift coefficient increases with the increase of Re L .The variation rule of pitch moment coefficient with Re L is not consistent with the lift coefficient because the increment of pitch moment with velocity is smaller than the second-squares increment of velocity.When Hd is greater than 7, the relationship between lift coefficient and pitching moment coefficient and Re L become weaker.

Effects of Hd and Re L on UUV Self-Propelled Performance
The relevant self-propulsion coefficients are derived from Equations ( 8)- (11).To analyze the effects of Hd and Re L on the interactions between the UUV and the thruster, parameters such as thrust reduction, hull efficiency, effective wake fraction, and rotational efficiency are selected for analysis.The first two of these parameters are the effects of the thruster on the UUV, and the last two reflect the effect of the UUV on the thruster.

•
The Variation of Thrust Reduction t describes the effect of propeller at the wake of the hull on the UUV.Due to the rotation of the thruster at the tail of the UUV, velocity of flow near the tail increases and the static pressure here decreases relative to the static pressure of the bare UUV, resulting in an increase in the drag of the UUV.For the increased UUV drag, its convention is generally attributed to the reduced thrust of the propellor.
As can be seen from Figure 14, t decreases firstly with the increasing of Hd, then increases slowly and finally tends to stable.Near the seabed, the wall blocking effect of the seabed increases the pressure difference between the lead and tail of the UUV.Additionally, there is an abrupt change in the t curve for Hd equal to 2.5.This is because that although the trends of R 0 , R and T with Hd are consistent, the inconsistency of the rapidity of their changes leads to an abrupt change in t when Hd is 2.5.In addition, comparing the t values for different Hd, it is found that the relative change of t with respect to that at Hd of 10 is less than 1%.It suggests that the effect of the seabed wall on t is relatively weak.It is difficult to capture the small changes by numerical methods.t increases with the increase of Re L .   Figure 16 shows the pressure distribution of UUV at a cross section with inco      Figure 16 shows the pressure distribution of UUV at a cross section with incoming flow velocity of 0.5 m/s and Hd of 1.5, 2, 2.5, 3, 5, 7, 9, 10, where (a) column is the pressure distribution of bare UUV and (b) column is the pressure distribution of UUV together with propeller.Comparing the two columns (a) and (b), it can be seen that the pressure at the trailing end of the UUV is reduced due to the suction of the propeller, which in turn causes an increase in the drag of the UUV.• Hull Efficiency depicts the hull efficiency and represents the combined effect of thrust reduction and wake flow.As can be seen from Figure 17, this efficiency decreases with increasing .As Hd increases, this efficiency gradually increases and remains stable when Hd is greater than 3.

•
Hull Efficiency η H depicts the hull efficiency and represents the combined effect of thrust reduction and wake flow.As can be seen from Figure 17, this efficiency decreases with increasing Re L .As Hd increases, this efficiency gradually increases and remains stable when Hd is greater than 3.

•
Effective Wake Fraction of Self-Propelled UUV As shown in Figure 18, increases slowly with the increasing of Hd and bilizes, and decreases with the increase of .Figure 19 shows the distribu at the UUV propeller disc with Hd of 1.5 and oncoming velocities from top of 0.5 m/s, 1 m/s, 2 m/s, and 4 m/s.(a) is the nominal wake at the bare UUV pro and (b) is the effective wake at the propeller disc.The in the cloud plot is from Equation (18).
is the speed at the propeller disc.Additionally, rep UUV's speed.

•
Effective Wake Fraction of Self-Propelled UUV As shown in Figure 18, w increases slowly with the increasing of Hd and finally stabilizes, and w decreases with the increase of Re L . Figure 19 shows the distribution of the w at the UUV propeller disc with Hd of 1.5 and oncoming velocities from top to bottom of 0.5 m/s, 1 m/s, 2 m/s, and 4 m/s.(a) is the nominal wake at the bare UUV propeller disc, and (b) is the effective wake at the propeller disc.The w in the cloud plot is calculated from Equation (18).V A is the speed at the propeller disc.Additionally, V represents the UUV's speed.• Effective Wake Fraction of Self-Propelled UUV As shown in Figure 18, increases slowly with the increasing of Hd an bilizes, and decreases with the increase of .Figure 19 shows the distrib at the UUV propeller disc with Hd of 1.5 and oncoming velocities from to of 0.5 m/s, 1 m/s, 2 m/s, and 4 m/s.(a) is the nominal wake at the bare UUV pr and (b) is the effective wake at the propeller disc.The in the cloud plot from Equation (18).
is the speed at the propeller disc.Additionally, re UUV's speed.20 shows that the effective wake flow is smaller than t wake flow, and the effective wake flow is weak with the variation of Hd, which be seen in the cloud diagram.

Efficiency of Propeller Rotation
The efficiency increases with the increase of Hd at low and it is a pendent of Hd when Hd is greater than 7.It increases with the increase of , the efficiency is weakly correlated with Hd, as shown in Figure 21. Figure 19 shows that the effective wake is smaller than the nominal wake, and the effective wake decreases with increasing Re L , which is consistent with Figure 18.The shear layer formed at the trailing edge of the UUV fin plate and the turbulent boundary layer grows along the tail of the UUV together form the cross-shaped wake field at the propeller disc in Figure 19a.
Figure 20 shows the distribution of the wake fraction w at the propeller disc of UUV with 8 kinds of different Hd at the oncoming velocity of 0.5 m/s, (a) the nominal wake at the propeller disc of bare UUV and (b) the effective wake at the propeller disc of self-propelled UUV. Figure 20 shows that the effective wake flow is smaller than the nominal wake flow, and the effective wake flow is weak with the variation of Hd, which can hardly be seen in the cloud diagram.

• Efficiency of Propeller Rotation
The efficiency η i increases with the increase of Hd at low Re L and it is almost independent of Hd when Hd is greater than 7.It increases with the increase of Re L .At large Re L , the efficiency is weakly correlated with Hd, as shown in Figure 21.

Conclusions
The variation trend of hydrodynamic performance of bare UUV and self-propelled UUV with Hd and Re L is similar.The coefficient of drag decreases with the increase of Hd and finally tends to be stable.The absolute value of lift coefficient decreases with the increase of Hd and finally tends to be zeros.The coefficient of pitching moment decreases with the increase of Hd and finally tends to be zeros.The hydrodynamic coefficients of UUV have a significant rate of variation near the seabed.The drag coefficient decreases with increasing Re L .The absolute value of the lift coefficient increases with increasing Re L .The pitching moment coefficient changes in the opposite trend of the lift coefficient.The drag coefficient of the self-propelled UUV increases by 16.74~17.88%.Additionally, lift coefficient, pitching moment coefficient of the self-propelled UUV increase by 69~130% and 100~300%, respectively, compared with that of bare UUV when UUV near the seabed.Thrust reduction, hull efficiency, wake fraction and rotational efficiency are weakly correlated with Hd.It is found that the thrust reduction and rotational efficiency increase with increasing Re L , and the hull efficiency and associated flow fraction decrease with increasing Re L .
The case of UUV for straight sailing were discussed for the paper.In the future, it is necessary to study the effect of seabed on the hydrodynamic performance of the self-propelled UUV with variations of the angle of attack.
Inlet: Uniform incoming flow; = ; Outlet: Pressure outlet; = 0; Non-slip wall surfaces: For the hull, fin and propellers and subsea walls.= 0; Slip wall surfaces: For the top of the flow field.= 0; Symmetric boundary: For the side of the flow fields.= • = 0.

Figure 2 .
Figure 2. Schematic of flow field.Inlet: Uniform incoming flow; → u = → u c ; Outlet: Pressure outlet; p = 0; Non-slip wall surfaces: For the hull, fin and propellers and subsea walls.

→ u = 0 ;
Slip wall surfaces: For the top of the flow field.∂ → u ∂x i = 0; Symmetric boundary: For the side of the flow fields.u n = → u c • → n = 0.

Figure 3 .
Figure 3. Grid settings.(a) The longitudinal profile grid of the domain with HD of 1.5; (b) UUV nose grid; (c) UUV tail grid; (d) propeller grid.

Figure 3 .
Figure 3. Grid settings.(a) The longitudinal profile grid of the domain with HD of 1.5; (b) UUV nose grid; (c) UUV tail grid; (d) propeller grid.

Figure 5 .
Figure 5. Verification of open water performance of DTMB4119.

Figure 5 .
Figure 5. Verification of open water performance of DTMB4119.

Figure 9 .
Figure 9.distribution of bare UUV at cross section a for different Hd.

Figure 9 .
Figure 9.distribution of bare UUV at cross section a for different Hd.

Figure 9 .
Figure 9. distribution of bare UUV at cross section a for different Hd.Figure 9. C p distribution of bare UUV at cross section a for different Hd.

Figure 9 .
Figure 9. distribution of bare UUV at cross section a for different Hd.Figure 9. C p distribution of bare UUV at cross section a for different Hd.

Figure 10 .
Figure 10.distribution of bare UUV at cross section a for different .

4. 2 .
Hydrodynamic Characteristics of Self-Propelled UUV near the Seabed and the Self-Propelled Performences 4.2.1.Hydrodynamic Characteristics of the Self-Propelled UUV in Relation to the Hd and

Figure 10 .
Figure 10.C p distribution of bare UUV at cross section a for different Re L .4.1.1.Correlations between Hydrodynamic Characteristics of Bare UUV and Hd • Drag Characteristics

4. 2 .
Hydrodynamic Characteristics of Self-Propelled UUV near the Seabed and the Self-Propelled Performences 4.2.1.Hydrodynamic Characteristics of the Self-Propelled UUV in Relation to the Hd and Re L

Figure 15
Figure 15 shows the pressure distribution of self-propelled UUV at cross-secti with Hd = 1.5 and oncoming flow velocities of 0.5, 1, 2, 4 m/s, where (a) is the pre distribution of bare UUV and (b) is the pressure distribution of UUV with propeller.paring the two columns (a) and (b), it can be seen that the pressure at the tail end UUV decreases due to the suction effect of the propeller, which in turn causes an in in the drag of the vehicle, corresponding to a reduction in the thrust of the UUV.

Figure 15 .
Figure 15.The pressure clouds at the bare UUV as shown columns (a) and self−propelled UU shown columns (b) at cross−section a.

Figure 15 21 Figure 14 .
Figure 15 shows the pressure distribution of self-propelled UUV at cross-section (a) with Hd = 1.5 and oncoming flow velocities of 0.5, 1, 2, 4 m/s, where (a) is the pressure distribution of bare UUV and (b) is the pressure distribution of UUV with propeller.Comparing the two columns (a) and (b), it can be seen that the pressure at the tail end of the UUV decreases due to the suction effect of the propeller, which in turn causes an increase in the drag of the vehicle, corresponding to a reduction in the thrust of the UUV.

Figure 15
Figure 15 shows the pressure distribution of self-propelled UUV at cross-section (a) with Hd = 1.5 and oncoming flow velocities of 0.5, 1, 2, 4 m/s, where (a) is the pressure distribution of bare UUV and (b) is the pressure distribution of UUV with propeller.Comparing the two columns (a) and (b), it can be seen that the pressure at the tail end of the UUV decreases due to the suction effect of the propeller, which in turn causes an increase in the drag of the vehicle, corresponding to a reduction in the thrust of the UUV.

Figure 15 .
Figure 15.The pressure clouds at the bare UUV as shown columns (a) and self−propelled UUV as shown columns (b) at cross−section a.

Figure 16
Figure16shows the pressure distribution of UUV at a cross section with incoming flow velocity of 0.5 m/s and Hd of 1.5, 2, 2.5, 3, 5, 7, 9, 10, where (a) column is the pressure distribution of bare UUV and (b) column is the pressure distribution of UUV together with propeller.Comparing the two columns (a) and (b), it can be seen that the pressure at

Figure 15 .
Figure 15.The pressure clouds at the bare UUV as shown columns (a) and self−propelled UUV as shown columns (b) at cross−section a.

21 Figure 16 .
Figure 16.The pressure clouds at the bare UUV as shown columns (a) and self-propelled UUV as shown columns (b) at a cross-section.

Figure 16 .
Figure 16.The pressure clouds at the bare UUV as shown columns (a) and self-propelled UUV as shown columns (b) at a cross-section.

Figure 19
Figure19shows that the effective wake is smaller than the nominal w effective wake decreases with increasing , which is consistent with Fig shear layer formed at the trailing edge of the UUV fin plate and the turbule

Figure 19 .
Figure 19.Cloud diagram of the wake fraction at the propeller disc for various wh 1.5.Columns (a) are for the bare UUV, and columns (b) are for the self-propelled UUV Figure20shows the distribution of the wake fraction at the propeller d with 8 kinds of different Hd at the oncoming velocity of 0.5 m/s, (a) the nomin the propeller disc of bare UUV and (b) the effective wake at the propeller d propelled UUV.Figure20shows that the effective wake flow is smaller than t wake flow, and the effective wake flow is weak with the variation of Hd, which be seen in the cloud diagram.

Figure 19 .
Figure 19.Cloud diagram of the wake fraction at the propeller disc for various Re L when Hd is 1.5.Columns (a) are for the bare UUV, and columns (b) are for the self-propelled UUV.

Figure 20 .
Figure 20.Cloud diagram of the wake fraction at the propeller disc for various Hd for oncoming velocity of 0.5 m/s.Columns (a) are for the bare UUV, and columns (b) are for the self-propelled UUV.

Figure 20 .
Figure 20.Cloud diagram of the wake fraction at the propeller disc for various Hd for oncoming velocity of 0.5 m/s.Columns (a) are for the bare UUV, and columns (b) are for the self-propelled UUV.

Figure 20 .
Figure 20.Cloud diagram of the wake fraction at the propeller disc for various Hd for oncomin velocity of 0.5 m/s.Columns (a) are for the bare UUV, and columns (b) are for the self-propelle UUV.

Table 1 .
Main parameters of UUV hull.

Table 2 .
Detailed parameters of the propeller DTMB4119.

Table 1 .
Main parameters of UUV hull.

Table 2 .
Detailed parameters of the propeller DTMB4119.