# Layout Optimization of Two Autonomous Underwater Vehicles for Drag Reduction with a Combined CFD and Neural Network Method

Authors: Wenlong Tian; Zhaoyong Mao; Fuliang Zhao; Zhicao Zhao
Journal: Complexity
 (2017)

---

## Abstract

This paper presents an optimization method for the design of the layout of an autonomous underwater vehicles (AUV) fleet to minimize the drag force. The layout of the AUV fleet is defined by two nondimensional parameters. Firstly, three-dimensional computational fluid dynamics (CFD) simulations are performed on the fleets with different layout parameters and detailed information on the hydrodynamic forces and flow structures around the AUVs is obtained. Then, based on the CFD data, a back-propagation neural network (BPNN) method is used to describe the relationship between the layout parameters and the drag of the fleet. Finally, a genetic algorithm (GA) is chosen to obtain the optimal layout parameters which correspond to the minimum drag. The optimization results show that(1) the total drag of the AUV fleet can be reduced by 12% when the follower AUV is located directly behind the leader AUV and (2) the drag of the follower AUV can be reduced by 66% when it is by the side of the leader AUV.


---

## Body

## 1. Introduction

Autonomous underwater vehicles (AUVs) are a kind of self-sailing, self-executing underwater robots, which play an important role in expanding people’s knowledge of the ocean. AUVs are capable of a wide range of applications, such as pipeline inspection [1, 2], underwater search and rescue, mine-sweeping [3], and oceanographic exploration [4]. The tasks of current AUVs are relatively single, and in the future, AUVs will be able to perform a variety of tasks. Therefore, AUVs will be composed of more complex structures and control systems, as well as more sensors and functional modules. Another solution for complex tasks is to use a fleet of simple AUVs, among which each one has a specific function, so that the fleet can complete complex tasks with simple AUVs and at a lower cost.When the AUVs travel in a fleet, the layout of the fleet becomes a key problem. Observations of animal motions such as birds in flocks suggest that some energy benefit may be obtained by certain fleet shapes [5]. It is known to all that migrating geese fly in a V or I formation for long distance migrations (Figure 1). The leader goose produces high speed updrafts in wake, which help the follower geese to save energy. Similarly, the leader AUV in the fleet produces low speed wake and low pressure side flows. Therefore, the follower AUVs could possibly use the velocity difference or pressure difference to minimize the drag and save energy.Figure 1
(a) Migrating geese and (b) possible layout of an AUV fleet.
(a)
(b)The drag of an individual AUV can be reduced by optimizing the shape of the hull, propeller, and surface control [6, 7]. Alvarez et al. introduced an optimization method for the shape optimization of an AUV moving near the water surface [8]. They used a panel model to predict the wave resistance of a revolution body and a simulated annealing algorithm to optimize the geometric parameters that minimize the wave resistance. Kim et al. proposed a Computational Fluid Dynamics (CFD) method for the optimization of monohull ships with the minimum drag [9]. Joung et al. proposed a procedure using the CFD method and a Design of Experiments (DOE) optimization for the hull optimization [10]. Sun et al. established an energy consumption model for an underwater glider and optimized the shape of the glider with a combined CFD and Efficient Global Optimization (EGO) method [11].CFD simulations can predict the hydrodynamic performance of the AUV with high accuracy, however, at the cost of time. Combined optimization method based on the CFD, Back-Propagation Neural Network (BPNN), and Genetic Algorithm (GA) can efficiently save the simulation time to obtain the optimal target. BPNN and GA are used at some of the nonlinear data to find the optimal value [12, 13]. Safikhani et al. studied the multiobjective optimization of nanofluid flow in flat tubs, which combined with a CFD, artificial neural networks, and genetic algorithms [14]. They got important design information about nanofluids and flat tubes. Avcı et al. studied the optimization of the deign parameters of a home refrigerator using CFD and artificial neural network [15]; they showed that CFD simulation and the ANN can determine the best value for the refrigerator design.Inspired by the above research, this paper introduces an optimization procedure of an AUV fleet, which contains two AUVs, to minimize its drag. The layout of the fleet is defined by two nondimensional parameters. Three-dimensional CFD simulations are performed to find the drag of each case. Then BPNN and GA optimization method is used to find the optimal layout which has the minimum drag.


## 2. Geometry Configuration

### 2.1. The AUV Hull Design

The AUV considered in this paper has a torpedo-like axisymmetric shape, which is characterized with a length ofL=1.3 m and a maximum diameter of D=0.2 m, resulting in a slenderness ratio (the ratio of the length to the diameter) of 6.5. The two-dimensional sketch of the AUV is shown in Figure 2. The designed AUV is composed of a nose section, a middle section, and a tail section. The lengths of the three sections are 0.2 m, 0.7 m, and 0.4 m, respectively. In order to obtain a low-drag shape, the Myring Equation [17] is used to design the shape of the AUV nose. The Myring Equation is widely used in the design of AUV hulls [10]. The expression for the nose shape is determined by(1)rx=12D1-x-Lhx21/1.8,where x is the position along the rotation axis, r is the radius at a specific x, and Lh is the length of the nose. The tail section is smoothly transited from the mid-section to a circular cone, which has a cone angle of 30°. The diameter of the rear surface is 0.04 m. Table 1 shows the values of the dimensional parameters.Table 1
Dimensional parameters of the AUV.
D (m)
D
t (m)
L (m)
L
h (m)
L
c (m)
L
t (m)
θ (°)
0.20
0.04
1.30
0.20
0.70
0.40
30Figure 2
Shape of the AUV.


### 2.2. Parameters of the AUV Fleet

As shown in Figure3, a coordinate frame is assigned to the nose of the leader AUV, with its horizontal axis pointing right and its vertical axis pointing to the side of the follower AUV. The layout of the AUV fleet can be characterized by two parameters,a andb, which represent the relative location of two adjacent AUVs in the fleet. To eliminate the dimensional effects,a andb are normalized by dividing the length of the AUV,(2)a=yfL,b=xfL,wherexf andyf are the position of the follower AUV in the proposed coordinate frame.Figure 3
Parameter definition for the AUV fleet.A series of fleet layouts can be described by changing the two design parameters,a andb. Before the optimization procedure, 53 groups of CFD simulations are performed as sample points witha varying from 0.00 to 0.60 andb varying from 0.00 to 2.00. It should be noted that simulations for a≤0.2 and b≤1.00 are not performed because of the geometric interference between the two AUVs in these situations.This study aims optimize the drag of the AUV fleet. The drag is normalized in the following way:(3)Cdl=Fxl1/2ρSU2,Cdf=Fxf1/2ρSU2,Cdsum=Cdl+Cdf,where Cdl, Cdf, and Cdsum are the drag coefficient of the leader AUV, the follower AUV, and the fleet, respectively, ρ is the water density and is equal to 998 kg/m3, S is the area of the cross-section of the AUV, and U is the velocity of the AUV.


## 2.1. The AUV Hull Design

The AUV considered in this paper has a torpedo-like axisymmetric shape, which is characterized with a length ofL=1.3 m and a maximum diameter of D=0.2 m, resulting in a slenderness ratio (the ratio of the length to the diameter) of 6.5. The two-dimensional sketch of the AUV is shown in Figure 2. The designed AUV is composed of a nose section, a middle section, and a tail section. The lengths of the three sections are 0.2 m, 0.7 m, and 0.4 m, respectively. In order to obtain a low-drag shape, the Myring Equation [17] is used to design the shape of the AUV nose. The Myring Equation is widely used in the design of AUV hulls [10]. The expression for the nose shape is determined by(1)rx=12D1-x-Lhx21/1.8,where x is the position along the rotation axis, r is the radius at a specific x, and Lh is the length of the nose. The tail section is smoothly transited from the mid-section to a circular cone, which has a cone angle of 30°. The diameter of the rear surface is 0.04 m. Table 1 shows the values of the dimensional parameters.Table 1
Dimensional parameters of the AUV.
D (m)
D
t (m)
L (m)
L
h (m)
L
c (m)
L
t (m)
θ (°)
0.20
0.04
1.30
0.20
0.70
0.40
30Figure 2
Shape of the AUV.


## 2.2. Parameters of the AUV Fleet

As shown in Figure3, a coordinate frame is assigned to the nose of the leader AUV, with its horizontal axis pointing right and its vertical axis pointing to the side of the follower AUV. The layout of the AUV fleet can be characterized by two parameters,a andb, which represent the relative location of two adjacent AUVs in the fleet. To eliminate the dimensional effects,a andb are normalized by dividing the length of the AUV,(2)a=yfL,b=xfL,wherexf andyf are the position of the follower AUV in the proposed coordinate frame.Figure 3
Parameter definition for the AUV fleet.A series of fleet layouts can be described by changing the two design parameters,a andb. Before the optimization procedure, 53 groups of CFD simulations are performed as sample points witha varying from 0.00 to 0.60 andb varying from 0.00 to 2.00. It should be noted that simulations for a≤0.2 and b≤1.00 are not performed because of the geometric interference between the two AUVs in these situations.This study aims optimize the drag of the AUV fleet. The drag is normalized in the following way:(3)Cdl=Fxl1/2ρSU2,Cdf=Fxf1/2ρSU2,Cdsum=Cdl+Cdf,where Cdl, Cdf, and Cdsum are the drag coefficient of the leader AUV, the follower AUV, and the fleet, respectively, ρ is the water density and is equal to 998 kg/m3, S is the area of the cross-section of the AUV, and U is the velocity of the AUV.


## 3. CFD Method

The drag coefficients and the flow structures for each simulation case are obtained by the CFD approaches separately. A group of transient representations of the three-dimensional flow are carried out to improve the efficiency of the optimization process. Further, the results of this numerical model will be compared with existing experimental data for the model validation, which will be discussed in Section3.4.


### 3.1. Governing Equations

The prediction of the flow around the AUVs is based on the incompressible Navier-Stokes equations [18],(4)∂ρ∂t+∇ρv→=0,∂∂tρv→+∇ρv→v→=-∇p+∇τ--,where t represents time, v→ is the vector of velocity, p is pressure, and τ-- refers to the stress tensor.The shear stress transport (SST)k-ω turbulence model was selected to model the turbulence terms of the RANS equations [18]. The SSTk-ω turbulence model is able to model the transport of turbulent shear stress and gives accurate predictions on the onset and amount of flow separation under adverse pressure gradients and has been successfully used in the CFD simulation of complex flows [11, 19–22].


### 3.2. Computational Domain and Boundary Conditions

The computational domain is a finite space used to simulate the flow around the AUV. In order to minimize the effect of block caused by the AUV, a rectangular domain with a dimension of8L×20D×15D is selected, resulting in a maximum block ratio of 0.67%, as shown in Figure 4. The length of the domain is 8L and the leader AUV is placed in the mid-plane of the domain and 3L to the left boundary. The follower AUV is in the same horizontal plane with the leader AUV and is offset in the coordinate frame.Figure 4
Computational domain and boundary conditions.Steady and uniform velocity inlet boundary is set at the left surface of the domain. The magnitude of the inflow velocity is 6 m/s and a moderate turbulence intensity of 5% is chosen. Pressure outlet boundary is chosen for the right surface of the domain, which has the same turbulence intensity with the inlet boundary. Smooth wall conditions are imposed at the four-side surfaces, where the shear effects are neglected to minimize the influences of the walls. Standard wall conditions are applied to the surfaces of both AUVs.


### 3.3. Mesh Generation

In order to obtain more accurate results with smaller number of grids, the computational domain is meshed with structured hexahedral grids, as shown in Figure5. Because the flow around the AUVs is more complex and has higher velocity gradients, a higher density of grid is set in this place to capture the development of flow with accuracy. Particularly, prism layer grid elements (Figure 5(c)) are generated from the wall surfaces to refine the mesh quality near the boundary layer flows. The initial height of the prism layer above the wall surfaces is an important indicator for the quality of the mesh. In this study, we use the nondimensional y+ value to evaluate the quality of the prism layer and a maximum value of 1 is chosen for all meshes. For a smoother transition of the adjacent elements, a general mesh growth rate of 1.2 is chosen for all grids.Figure 5
Details of the computation mesh: (a) overview; (b) grids around the vehicle; and (c) prism layer grid elements.
(a)
(b)
(c)


### 3.4. Verification and Validation

A mesh verification study is carried out to determine the proper density of mesh. The verification simulations are performed on a single AUV with different meshes. Three meshes, with approximately 3.76 million (fine mesh), 2.80 million (mediate mesh), and 1.67 million (coarse mesh) elements, respectively, are generated for the grid resolution verification. The drag coefficients obtained by the three meshes are listed in Table2. It is found that the coarse mesh predicts slightly higher Cd and the other two give almost the same results. Therefore, the mediate mesh is chosen for the following simulations. It should be noted that, for the simulations of the AUV fleet, the grid density is kept the same with that of the single AUV and the number of elements is about 5.5 million due to the increased number of AUVs.Table 2
Results of the mesh verification study.
Mesh
Grid number
C
d
Coarse
1670000
0.0786
Mediate
2800000
0.0774
Fine
3760000
0.0771A validation study is performed for the validation of the proposed numerical model. Simulations are performed according to an experiment of twin parallel bare prolate spheroids with transverse separations [16] (a varies and b=0 in current study), as shown in Figure 6. The length of the model was 1.2 m and the maximum diameter was 0.2 m, which is very close to the size of the AUV in this study. The testing inlet velocity was 40 m/s, corresponding to a Reynolds number of 3.2×106 according to the length of the model. This Reynolds number was also close to the current simulation. The calculated coefficients of drag and side force of the top spheroid are then compared with those obtained by experiment [16] and a previous CFD study [23] (Figure 7). It is noticed that the current CFD model provides more accurate results than that presented by Molland and Utama [23]. The maximum relative errors of the predicted drag and side force are 3.3% and 9.1%, respectively. This proves that the CFD model in this study can predict the forces of two AUVs with acceptable accuracy.Figure 6
Wind tunnel test of the twin spheroids [16].Figure 7
Results of the validation study: (a) coefficients of drag and (b) coefficients of side force.
(a)
(b)


## 3.1. Governing Equations

The prediction of the flow around the AUVs is based on the incompressible Navier-Stokes equations [18],(4)∂ρ∂t+∇ρv→=0,∂∂tρv→+∇ρv→v→=-∇p+∇τ--,where t represents time, v→ is the vector of velocity, p is pressure, and τ-- refers to the stress tensor.The shear stress transport (SST)k-ω turbulence model was selected to model the turbulence terms of the RANS equations [18]. The SSTk-ω turbulence model is able to model the transport of turbulent shear stress and gives accurate predictions on the onset and amount of flow separation under adverse pressure gradients and has been successfully used in the CFD simulation of complex flows [11, 19–22].


## 3.2. Computational Domain and Boundary Conditions

The computational domain is a finite space used to simulate the flow around the AUV. In order to minimize the effect of block caused by the AUV, a rectangular domain with a dimension of8L×20D×15D is selected, resulting in a maximum block ratio of 0.67%, as shown in Figure 4. The length of the domain is 8L and the leader AUV is placed in the mid-plane of the domain and 3L to the left boundary. The follower AUV is in the same horizontal plane with the leader AUV and is offset in the coordinate frame.Figure 4
Computational domain and boundary conditions.Steady and uniform velocity inlet boundary is set at the left surface of the domain. The magnitude of the inflow velocity is 6 m/s and a moderate turbulence intensity of 5% is chosen. Pressure outlet boundary is chosen for the right surface of the domain, which has the same turbulence intensity with the inlet boundary. Smooth wall conditions are imposed at the four-side surfaces, where the shear effects are neglected to minimize the influences of the walls. Standard wall conditions are applied to the surfaces of both AUVs.


## 3.3. Mesh Generation

In order to obtain more accurate results with smaller number of grids, the computational domain is meshed with structured hexahedral grids, as shown in Figure5. Because the flow around the AUVs is more complex and has higher velocity gradients, a higher density of grid is set in this place to capture the development of flow with accuracy. Particularly, prism layer grid elements (Figure 5(c)) are generated from the wall surfaces to refine the mesh quality near the boundary layer flows. The initial height of the prism layer above the wall surfaces is an important indicator for the quality of the mesh. In this study, we use the nondimensional y+ value to evaluate the quality of the prism layer and a maximum value of 1 is chosen for all meshes. For a smoother transition of the adjacent elements, a general mesh growth rate of 1.2 is chosen for all grids.Figure 5
Details of the computation mesh: (a) overview; (b) grids around the vehicle; and (c) prism layer grid elements.
(a)
(b)
(c)


## 3.4. Verification and Validation

A mesh verification study is carried out to determine the proper density of mesh. The verification simulations are performed on a single AUV with different meshes. Three meshes, with approximately 3.76 million (fine mesh), 2.80 million (mediate mesh), and 1.67 million (coarse mesh) elements, respectively, are generated for the grid resolution verification. The drag coefficients obtained by the three meshes are listed in Table2. It is found that the coarse mesh predicts slightly higher Cd and the other two give almost the same results. Therefore, the mediate mesh is chosen for the following simulations. It should be noted that, for the simulations of the AUV fleet, the grid density is kept the same with that of the single AUV and the number of elements is about 5.5 million due to the increased number of AUVs.Table 2
Results of the mesh verification study.
Mesh
Grid number
C
d
Coarse
1670000
0.0786
Mediate
2800000
0.0774
Fine
3760000
0.0771A validation study is performed for the validation of the proposed numerical model. Simulations are performed according to an experiment of twin parallel bare prolate spheroids with transverse separations [16] (a varies and b=0 in current study), as shown in Figure 6. The length of the model was 1.2 m and the maximum diameter was 0.2 m, which is very close to the size of the AUV in this study. The testing inlet velocity was 40 m/s, corresponding to a Reynolds number of 3.2×106 according to the length of the model. This Reynolds number was also close to the current simulation. The calculated coefficients of drag and side force of the top spheroid are then compared with those obtained by experiment [16] and a previous CFD study [23] (Figure 7). It is noticed that the current CFD model provides more accurate results than that presented by Molland and Utama [23]. The maximum relative errors of the predicted drag and side force are 3.3% and 9.1%, respectively. This proves that the CFD model in this study can predict the forces of two AUVs with acceptable accuracy.Figure 6
Wind tunnel test of the twin spheroids [16].Figure 7
Results of the validation study: (a) coefficients of drag and (b) coefficients of side force.
(a)
(b)


## 4. Optimal Design Method

This paper establishes an optimization method based on the combination of CFD method, BPNN and GA. The drag coefficients and the flow structures for 53 different layouts are obtained by the CFD approaches separately. The results are then used to build the agent model betweenCd and the design parameters (a and b) using BPNN. Finally, GA is applied to search for the optimal Cd together with the values for the corresponding design parameters based on the agent model.


### 4.1. BP Neural Network

The artificial neural network is widely accepted as an alternative to providing solutions to complex and ambiguous problems [24, 25]. In this study, a back-propagation algorithm that optimizes weighted connections by allowing errors to propagate from the output layer to the input layer is used to train the network. The topology of the BPNN in present study is shown in Figure 8. The structure has two input layers, five hidden layers and one output layer. Each node of the structure represents a neuron, and the nodes between the layers and the layers are connected by weights.Figure 8
The topology of the BPNN.The BPNN is working based on the following principle:(
1
) According to CFD simulation results, determine the network input data and output data (X,Y). H is the data matrix of the hidden layer, as can be seen in(5)Xp=xp1,xp2,…,xpn,Hp=hp1,hp2,…,hpl,Yp=yp1,yp2,…,ypm,where n,m,l, and p are the numbers of inputs, outputs, hidden layers, and the training samples, respectively.(
2
) The outputs of the hidden layer are obtained by(6)hpi=f∑i=1nωijxpi-aj,j=1,2,…,l,where ωij is the connection between the input layer and the hidden layer.(
3
) The predicted outputs of the BPNN can be obtained by(7)y^pk=∑j=1lhpiω^jk-bk,k=1,2,…,m, where ω^jk is the connection between the hidden layer and the output layer.(
4
) According to the network forecast output y^k and the expected output yk, calculate the network prediction error ek. (8)ek=ypk-y^pk.(
5
) The weights are updated according to the network prediction error (ek).(9)ωij=ωij+ηhpi1-hpixi∑k=1mω^jkek,ω^jk=ω^jk+ηhpjek,where η is the learning efficiency of BPNN.(
6
) The thresholds are updated according to the network prediction error (ek).(10)aj=aj+ηhpj1-hpj∑k=1mω^jkek,bk=bk+ηek.(7) Return to step(2) if the optimization is still running.It is noted that all the iterative processes of the BP neural network have a mean square error (MSE). In addition, the average accuracy of the prediction(R) is defined by Kreith (2000). In other words, R represents the goodness of fit, which is used to measure the correlation between prediction output data and training samples data. The closer it is to the value of 1, the better the training network is.(11)MSE=1mp∑p=1p∑j=1my^pj-ypj2,R=1P∑i=1PRi=1P∑i=1Pypjy^pj.


### 4.2. Genetic Algorithm

In this paper, the GA tool in MATLAB is used for the optimization, which runs in the following principle:(1)
Initialize the population, calculate the fitness value, and find the best chromosome from the population.(2)
Iterative optimization:
(a)
Select: first, the solution to the problem is encoded by using the floating-point encoding. This function selects the chromosomes in each generation population for subsequent crossover and mutation. The method used is the roulette selection method.(b)
Crossover: this function is a random selection of two chromosomes, according to determine the crossover probability to determine whether the cross, and the cross position is also random.(c)
Mutation: this function performs the mutation operation. The mutation chromosomes and mutated positions are randomly selected. Finally, it will check the feasibility of chromosomes; otherwise, it will be recompiled.(d)
Result analysis.The optimal solution can be found after several generations. The establishment of BPNN approximate model and the numerical optimization of the genetic algorithm are shown in Figure9.Figure 9
Neural network and genetic algorithm flowchart.


## 4.1. BP Neural Network

The artificial neural network is widely accepted as an alternative to providing solutions to complex and ambiguous problems [24, 25]. In this study, a back-propagation algorithm that optimizes weighted connections by allowing errors to propagate from the output layer to the input layer is used to train the network. The topology of the BPNN in present study is shown in Figure 8. The structure has two input layers, five hidden layers and one output layer. Each node of the structure represents a neuron, and the nodes between the layers and the layers are connected by weights.Figure 8
The topology of the BPNN.The BPNN is working based on the following principle:(
1
) According to CFD simulation results, determine the network input data and output data (X,Y). H is the data matrix of the hidden layer, as can be seen in(5)Xp=xp1,xp2,…,xpn,Hp=hp1,hp2,…,hpl,Yp=yp1,yp2,…,ypm,where n,m,l, and p are the numbers of inputs, outputs, hidden layers, and the training samples, respectively.(
2
) The outputs of the hidden layer are obtained by(6)hpi=f∑i=1nωijxpi-aj,j=1,2,…,l,where ωij is the connection between the input layer and the hidden layer.(
3
) The predicted outputs of the BPNN can be obtained by(7)y^pk=∑j=1lhpiω^jk-bk,k=1,2,…,m, where ω^jk is the connection between the hidden layer and the output layer.(
4
) According to the network forecast output y^k and the expected output yk, calculate the network prediction error ek. (8)ek=ypk-y^pk.(
5
) The weights are updated according to the network prediction error (ek).(9)ωij=ωij+ηhpi1-hpixi∑k=1mω^jkek,ω^jk=ω^jk+ηhpjek,where η is the learning efficiency of BPNN.(
6
) The thresholds are updated according to the network prediction error (ek).(10)aj=aj+ηhpj1-hpj∑k=1mω^jkek,bk=bk+ηek.(7) Return to step(2) if the optimization is still running.It is noted that all the iterative processes of the BP neural network have a mean square error (MSE). In addition, the average accuracy of the prediction(R) is defined by Kreith (2000). In other words, R represents the goodness of fit, which is used to measure the correlation between prediction output data and training samples data. The closer it is to the value of 1, the better the training network is.(11)MSE=1mp∑p=1p∑j=1my^pj-ypj2,R=1P∑i=1PRi=1P∑i=1Pypjy^pj.


## 4.2. Genetic Algorithm

In this paper, the GA tool in MATLAB is used for the optimization, which runs in the following principle:(1)
Initialize the population, calculate the fitness value, and find the best chromosome from the population.(2)
Iterative optimization:
(a)
Select: first, the solution to the problem is encoded by using the floating-point encoding. This function selects the chromosomes in each generation population for subsequent crossover and mutation. The method used is the roulette selection method.(b)
Crossover: this function is a random selection of two chromosomes, according to determine the crossover probability to determine whether the cross, and the cross position is also random.(c)
Mutation: this function performs the mutation operation. The mutation chromosomes and mutated positions are randomly selected. Finally, it will check the feasibility of chromosomes; otherwise, it will be recompiled.(d)
Result analysis.The optimal solution can be found after several generations. The establishment of BPNN approximate model and the numerical optimization of the genetic algorithm are shown in Figure9.Figure 9
Neural network and genetic algorithm flowchart.


## 5. CFD Results

### 5.1. Force Performance

In order to more intuitively show the influences of the fleet layout on the drags of the two AUVs, the coefficients of drags are expressed in a nondimensional drag ratio by dividingCd0 of a single AUV (0.0774 in Table 2):(12)rl=CdlCd0,rf=CdfCd0,rsum=Cdl+CdfCd0.The drag ratios for the 53 different fleet layouts are shown in Figure10. The configuration of the AUV fleet can be divided into four regions depending on the drag ratios: the Parallel Region, the Pull Region, the Push Region, and the Tandem Region.Figure 10
Drag ratios for different AUV fleet layouts.


### 5.2. Parallel Region

The Parallel Region locates atb=0. AUVs in this region generate a higher drag than that in an infinite domain. Besides, the drag ratios for both AUVs are almost the same as each other due to the symmetrical geometric positions. For example, the drag ratios of AUVs at a=0.20 and b=0.00 are 5.52% higher than that of a single AUV. It can also be observed that as the latitudinal offset increases, the drag ratios gradually reduce. At a=0.60 and b=0.00, the drags of the AUVs are only approximately 0.6% higher than that of a single AUV. These results are consistent with previous studies [16, 26] which suggested a minimum latitudinal offset of 0.5L to preclude the interactions between two AUVs.The contours of pressure and velocity around the two parallel AUVs are shown in Figure11. Three different latitudinal offsets, a=0.2, 0.4, and 0.6, are presented and compared. It can be seen from the figure that the velocity of the fluid is obviously increased in the range between the two AUVs. According to Bernoulli’s Equation, larger fluid velocity corresponds to lower pressure, which can be verified by the pressure distributions in Figure 11(b). For example, in the pressure diagram of a=0.2, the pressure on the nose of the AUV is significantly reduced, which eventually leads to an increase in the drag. As a increases, the interactions between the AUVs get weaker and the drag of the AUV is gradually restored to normal levels.Figure 11
Contours of (a) velocity and (b) pressure around two parallel AUVs at different latitudinal offsets.
(a)
(b)


### 5.3. Tandem Region

The Tandem Region is where the follower AUV locates just behind the leader AUV (a=0). The drag ratios of the AUVs in this region are greatly reduced, especially for the leader AUV. For example, at a=0.00 and b=1.50, the drag reduction is 6.18% for the leader AUV and is 2.52% for the follower AUV. It should be noted that the follower AUV experiences an 9.72% increase in drag at a=0.00 and b=1.25, but the total drag of the AUV fleet reduces.The contours of pressure and velocity around the two tandem AUVs are shown in Figure12. Three different longitudinal offsets, b=1.25, 1.50, and 1.75, are presented and compared. When the follower AUV is located in tandem in the wake of the leader AUV, a certain amount of pressure is recovered on the tail of the leader AUV, which contributes to the reduction of the pressure drag of the leader AUV, while for the follower AUV, the velocity of the upstream flow is reduced by the leader AUV, resulting in a reduction of the nose pressure and a decrease in final resistance.Figure 12
Contours of (a) velocity and (b) pressure around two parallel AUVs at different latitudinal offsets.
(a)
(b)


### 5.4. Pull Region

The Pull Region is where the drag of the leader AUV is increased and that of the follower AUV is reduced. The Pull Region mainly locates between0.20≤a≤0.60 and 0.25≤b≤1.00. In this region, it is likely that the leader AUV pulls the follower AUV to move, so that the drag of the leader AUV is increased and the drag of the follower AUV changes oppositely.The contours of pressure and velocity around the two AUVs in the Pull Region are shown in Figure13. Three different locations with a=0.2, 0.3, and 0.4, b=0.25, 0.50, and 0.75, respectively, are presented and compared. The longitudinal offset is within one AUV length and the nose of the follower AUV locates at the lower pressure region by the side of the leader AUV. Under the effects of the follower AUV, the pressure on the side surface of the leader AUV is reduced, resulting in an augment in the drag. Besides, as the a increases, the interactions between the two AUVs are weaker and the variation in the drags of both AUVs is smaller.Figure 13
Contours of (a) velocity and (b) pressure around two AUVs in the Pull Region.
(a)
(b)


### 5.5. Push Region

The Push Region is where the drag of the leader AUV is reduced and that of the follower AUV is augmented. The Push Region mainly locates between0.10≤a≤0.60 and b≥1.25. In this region, it is likely that the follower AUV pushes the leader AUV to move, so that the drag of the leader AUV is reduced and the drag of the follower AUV changes oppositely.The contours of pressure and velocity around the two AUVs in the Push Region are shown in Figure14. Three different locations with a=0.2, 0.3, and 0.4, b=1.25, 1.50, and 1.75, respectively, are presented and compared. The reason for the variation of the drags on the two AUVs can be explained. The obstacle flow around the nose of the follower AUV hinders the separation on the tail of the leader AUV; therefore, a certain amount of pressure is recovered on the tail of the leader AUV, which contributes to a reduction of the pressure drag of the leader AUV. A reaction force is then applied to the follower AUV, leading to an increase in the drag on the follower AUV.Figure 14
Contours of (a) velocity and (b) pressure around two AUVs in the Pull Region.
(a)
(b)


### 5.6. Data for the BPNN Training Samples

The predicted drag ratios for the AUV fleet obtained by CFD are given in Table3. It should be noted that 10 samples for a≤0.20 and b≤1.00 are added to complete an orthographical sample domain. Since it is nonphysical to arrange the two AUVs in this area, a high drag ratio of 2.0000 is assigned to the AUV which locates in this region. Therefore, there are 63 samples in total for the training of the BPNN.Table 3
CFD data for the training samples.
Sample
a
b
r
l
r
f
r
s
u
m
(1)
0.00
0.00
2.0000
2.0000
4.0000
(2)
0.00
0.25
2.0000
2.0000
4.0000
(3)
0.00
0.50
2.0000
2.0000
4.0000
(4)
0.00
0.75
2.0000
2.0000
4.0000
(5)
0.00
1.00
2.0000
2.0000
4.0000
(6)
0.00
1.25
0.7910
1.0972
1.8882
(7)
0.00
1.50
0.9392
0.9748
1.9141
(8)
0.00
1.75
0.9757
0.9602
1.9359
(9)
0.00
2.00
0.9885
0.9629
1.9514
(10)
0.10
0.00
2.0000
2.0000
4.0000
(11)
0.10
0.25
2.0000
2.0000
4.0000
(12)
0.10
0.50
2.0000
2.0000
4.0000
(13)
0.10
0.75
2.0000
2.0000
4.0000
(14)
0.10
1.00
2.0000
2.0000
4.0000
(15)
0.10
1.25
0.8245
1.1725
1.9970
(16)
0.10
1.50
0.9432
1.0604
2.0036
(17)
0.10
1.75
0.9765
1.0296
2.0061
(18)
0.10
2.00
0.9887
1.0188
2.0075
(19)
0.20
0.00
1.0552
1.0552
2.1104
(20)
0.20
0.25
1.7136
0.3598
2.0733
(21)
0.20
0.50
1.5153
0.5364
2.0517
(22)
0.20
0.75
1.4027
0.5841
1.9868
(23)
0.20
1.00
0.7591
1.2507
2.0098
(24)
0.20
1.25
0.8671
1.1107
1.9778
(25)
0.20
1.50
0.9450
1.0398
1.9848
(26)
0.20
1.75
0.9728
1.0150
1.9879
(27)
0.20
2.00
0.9841
1.0050
1.9891
(28)
0.30
0.00
1.0247
1.0247
2.0494
(29)
0.30
0.25
1.4041
0.6326
2.0367
(30)
0.30
0.50
1.3595
0.6702
2.0297
(31)
0.30
0.75
1.1989
0.7884
1.9873
(32)
0.30
1.00
0.9257
1.0920
2.0177
(33)
0.30
1.25
0.9135
1.0671
1.9806
(34)
0.30
1.50
0.9557
1.0297
1.9853
(35)
0.30
1.75
0.9761
1.0118
1.9879
(36)
0.30
2.00
0.9852
1.0037
1.9889
(37)
0.40
0.00
1.0146
1.0146
2.0292
(38)
0.40
0.25
1.2386
0.7859
2.0245
(39)
0.40
0.50
1.2489
0.7770
2.0259
(40)
0.40
0.75
1.1302
0.8617
1.9920
(41)
0.40
1.00
0.9947
1.0283
2.0229
(42)
0.40
1.25
0.9487
1.0354
1.9841
(43)
0.40
1.50
0.9669
1.0200
1.9869
(44)
0.40
1.75
0.9800
1.0086
1.9886
(45)
0.40
2.00
0.9871
1.0025
1.9896
(46)
0.50
0.00
1.0094
1.0095
2.0189
(47)
0.50
0.25
1.1488
0.8699
2.0187
(48)
0.50
0.50
1.1753
0.8506
2.0259
(49)
0.50
0.75
1.0950
0.9000
1.9951
(50)
0.50
1.00
1.0243
1.0032
2.0275
(51)
0.50
1.25
0.9723
1.0146
1.9868
(52)
0.50
1.50
0.9769
1.0118
1.9887
(53)
0.50
1.75
0.9843
1.0057
1.9900
(54)
0.50
2.00
0.9891
1.0016
1.9907
(55)
0.60
0.00
1.0062
1.0065
2.0127
(56)
0.60
0.25
1.0963
0.9170
2.0132
(57)
0.60
0.50
1.1263
0.8991
2.0253
(58)
0.60
0.75
1.0724
0.9245
1.9969
(59)
0.60
1.00
1.0327
0.9907
2.0235
(60)
0.60
1.25
0.9873
1.0018
1.9891
(61)
0.60
1.50
0.9850
1.0052
1.9903
(62)
0.60
1.75
0.9883
1.0029
1.9912
(63)
0.60
2.00
0.9911
1.0000
1.9911


## 5.1. Force Performance

In order to more intuitively show the influences of the fleet layout on the drags of the two AUVs, the coefficients of drags are expressed in a nondimensional drag ratio by dividingCd0 of a single AUV (0.0774 in Table 2):(12)rl=CdlCd0,rf=CdfCd0,rsum=Cdl+CdfCd0.The drag ratios for the 53 different fleet layouts are shown in Figure10. The configuration of the AUV fleet can be divided into four regions depending on the drag ratios: the Parallel Region, the Pull Region, the Push Region, and the Tandem Region.Figure 10
Drag ratios for different AUV fleet layouts.


## 5.2. Parallel Region

The Parallel Region locates atb=0. AUVs in this region generate a higher drag than that in an infinite domain. Besides, the drag ratios for both AUVs are almost the same as each other due to the symmetrical geometric positions. For example, the drag ratios of AUVs at a=0.20 and b=0.00 are 5.52% higher than that of a single AUV. It can also be observed that as the latitudinal offset increases, the drag ratios gradually reduce. At a=0.60 and b=0.00, the drags of the AUVs are only approximately 0.6% higher than that of a single AUV. These results are consistent with previous studies [16, 26] which suggested a minimum latitudinal offset of 0.5L to preclude the interactions between two AUVs.The contours of pressure and velocity around the two parallel AUVs are shown in Figure11. Three different latitudinal offsets, a=0.2, 0.4, and 0.6, are presented and compared. It can be seen from the figure that the velocity of the fluid is obviously increased in the range between the two AUVs. According to Bernoulli’s Equation, larger fluid velocity corresponds to lower pressure, which can be verified by the pressure distributions in Figure 11(b). For example, in the pressure diagram of a=0.2, the pressure on the nose of the AUV is significantly reduced, which eventually leads to an increase in the drag. As a increases, the interactions between the AUVs get weaker and the drag of the AUV is gradually restored to normal levels.Figure 11
Contours of (a) velocity and (b) pressure around two parallel AUVs at different latitudinal offsets.
(a)
(b)


## 5.3. Tandem Region

The Tandem Region is where the follower AUV locates just behind the leader AUV (a=0). The drag ratios of the AUVs in this region are greatly reduced, especially for the leader AUV. For example, at a=0.00 and b=1.50, the drag reduction is 6.18% for the leader AUV and is 2.52% for the follower AUV. It should be noted that the follower AUV experiences an 9.72% increase in drag at a=0.00 and b=1.25, but the total drag of the AUV fleet reduces.The contours of pressure and velocity around the two tandem AUVs are shown in Figure12. Three different longitudinal offsets, b=1.25, 1.50, and 1.75, are presented and compared. When the follower AUV is located in tandem in the wake of the leader AUV, a certain amount of pressure is recovered on the tail of the leader AUV, which contributes to the reduction of the pressure drag of the leader AUV, while for the follower AUV, the velocity of the upstream flow is reduced by the leader AUV, resulting in a reduction of the nose pressure and a decrease in final resistance.Figure 12
Contours of (a) velocity and (b) pressure around two parallel AUVs at different latitudinal offsets.
(a)
(b)


## 5.4. Pull Region

The Pull Region is where the drag of the leader AUV is increased and that of the follower AUV is reduced. The Pull Region mainly locates between0.20≤a≤0.60 and 0.25≤b≤1.00. In this region, it is likely that the leader AUV pulls the follower AUV to move, so that the drag of the leader AUV is increased and the drag of the follower AUV changes oppositely.The contours of pressure and velocity around the two AUVs in the Pull Region are shown in Figure13. Three different locations with a=0.2, 0.3, and 0.4, b=0.25, 0.50, and 0.75, respectively, are presented and compared. The longitudinal offset is within one AUV length and the nose of the follower AUV locates at the lower pressure region by the side of the leader AUV. Under the effects of the follower AUV, the pressure on the side surface of the leader AUV is reduced, resulting in an augment in the drag. Besides, as the a increases, the interactions between the two AUVs are weaker and the variation in the drags of both AUVs is smaller.Figure 13
Contours of (a) velocity and (b) pressure around two AUVs in the Pull Region.
(a)
(b)


## 5.5. Push Region

The Push Region is where the drag of the leader AUV is reduced and that of the follower AUV is augmented. The Push Region mainly locates between0.10≤a≤0.60 and b≥1.25. In this region, it is likely that the follower AUV pushes the leader AUV to move, so that the drag of the leader AUV is reduced and the drag of the follower AUV changes oppositely.The contours of pressure and velocity around the two AUVs in the Push Region are shown in Figure14. Three different locations with a=0.2, 0.3, and 0.4, b=1.25, 1.50, and 1.75, respectively, are presented and compared. The reason for the variation of the drags on the two AUVs can be explained. The obstacle flow around the nose of the follower AUV hinders the separation on the tail of the leader AUV; therefore, a certain amount of pressure is recovered on the tail of the leader AUV, which contributes to a reduction of the pressure drag of the leader AUV. A reaction force is then applied to the follower AUV, leading to an increase in the drag on the follower AUV.Figure 14
Contours of (a) velocity and (b) pressure around two AUVs in the Pull Region.
(a)
(b)


## 5.6. Data for the BPNN Training Samples

The predicted drag ratios for the AUV fleet obtained by CFD are given in Table3. It should be noted that 10 samples for a≤0.20 and b≤1.00 are added to complete an orthographical sample domain. Since it is nonphysical to arrange the two AUVs in this area, a high drag ratio of 2.0000 is assigned to the AUV which locates in this region. Therefore, there are 63 samples in total for the training of the BPNN.Table 3
CFD data for the training samples.
Sample
a
b
r
l
r
f
r
s
u
m
(1)
0.00
0.00
2.0000
2.0000
4.0000
(2)
0.00
0.25
2.0000
2.0000
4.0000
(3)
0.00
0.50
2.0000
2.0000
4.0000
(4)
0.00
0.75
2.0000
2.0000
4.0000
(5)
0.00
1.00
2.0000
2.0000
4.0000
(6)
0.00
1.25
0.7910
1.0972
1.8882
(7)
0.00
1.50
0.9392
0.9748
1.9141
(8)
0.00
1.75
0.9757
0.9602
1.9359
(9)
0.00
2.00
0.9885
0.9629
1.9514
(10)
0.10
0.00
2.0000
2.0000
4.0000
(11)
0.10
0.25
2.0000
2.0000
4.0000
(12)
0.10
0.50
2.0000
2.0000
4.0000
(13)
0.10
0.75
2.0000
2.0000
4.0000
(14)
0.10
1.00
2.0000
2.0000
4.0000
(15)
0.10
1.25
0.8245
1.1725
1.9970
(16)
0.10
1.50
0.9432
1.0604
2.0036
(17)
0.10
1.75
0.9765
1.0296
2.0061
(18)
0.10
2.00
0.9887
1.0188
2.0075
(19)
0.20
0.00
1.0552
1.0552
2.1104
(20)
0.20
0.25
1.7136
0.3598
2.0733
(21)
0.20
0.50
1.5153
0.5364
2.0517
(22)
0.20
0.75
1.4027
0.5841
1.9868
(23)
0.20
1.00
0.7591
1.2507
2.0098
(24)
0.20
1.25
0.8671
1.1107
1.9778
(25)
0.20
1.50
0.9450
1.0398
1.9848
(26)
0.20
1.75
0.9728
1.0150
1.9879
(27)
0.20
2.00
0.9841
1.0050
1.9891
(28)
0.30
0.00
1.0247
1.0247
2.0494
(29)
0.30
0.25
1.4041
0.6326
2.0367
(30)
0.30
0.50
1.3595
0.6702
2.0297
(31)
0.30
0.75
1.1989
0.7884
1.9873
(32)
0.30
1.00
0.9257
1.0920
2.0177
(33)
0.30
1.25
0.9135
1.0671
1.9806
(34)
0.30
1.50
0.9557
1.0297
1.9853
(35)
0.30
1.75
0.9761
1.0118
1.9879
(36)
0.30
2.00
0.9852
1.0037
1.9889
(37)
0.40
0.00
1.0146
1.0146
2.0292
(38)
0.40
0.25
1.2386
0.7859
2.0245
(39)
0.40
0.50
1.2489
0.7770
2.0259
(40)
0.40
0.75
1.1302
0.8617
1.9920
(41)
0.40
1.00
0.9947
1.0283
2.0229
(42)
0.40
1.25
0.9487
1.0354
1.9841
(43)
0.40
1.50
0.9669
1.0200
1.9869
(44)
0.40
1.75
0.9800
1.0086
1.9886
(45)
0.40
2.00
0.9871
1.0025
1.9896
(46)
0.50
0.00
1.0094
1.0095
2.0189
(47)
0.50
0.25
1.1488
0.8699
2.0187
(48)
0.50
0.50
1.1753
0.8506
2.0259
(49)
0.50
0.75
1.0950
0.9000
1.9951
(50)
0.50
1.00
1.0243
1.0032
2.0275
(51)
0.50
1.25
0.9723
1.0146
1.9868
(52)
0.50
1.50
0.9769
1.0118
1.9887
(53)
0.50
1.75
0.9843
1.0057
1.9900
(54)
0.50
2.00
0.9891
1.0016
1.9907
(55)
0.60
0.00
1.0062
1.0065
2.0127
(56)
0.60
0.25
1.0963
0.9170
2.0132
(57)
0.60
0.50
1.1263
0.8991
2.0253
(58)
0.60
0.75
1.0724
0.9245
1.9969
(59)
0.60
1.00
1.0327
0.9907
2.0235
(60)
0.60
1.25
0.9873
1.0018
1.9891
(61)
0.60
1.50
0.9850
1.0052
1.9903
(62)
0.60
1.75
0.9883
1.0029
1.9912
(63)
0.60
2.00
0.9911
1.0000
1.9911


## 6. Neural Network and Genetic Algorithm Optimization

### 6.1. Optimization of the AUV Fleet

When the AUV fleet executes tasks underwater, a minimum drag is expected so that the fleet will work for a longer time and a wider range. Therefore, the layout of the AUV fleet should be optimized to obtain a minimum drag.


#### 6.1.1. Neural Network Training

Using the data in Table3, training of the BPNN is carried out until the network meets the intended target. A comparison between the predicted data and the expected data is shown in Figure 15. The training network predicts drag ratios with good agreement with the CFD data. Therefore, it can be said that the BPNN successfully model and predict the drag ratio of the AUV fleet.Figure 15
Comparisons of the drag ratio of the AUV fleet by BPNN and CFD.


#### 6.1.2. Comparisons of the Optimal Drag between CFD and BPNN

A CFD simulation of the AUV fleet for the optimal parameters is carried out for the comparisons between the two methods. The optimal results are shown in Table4. The results of the BPNN and the CFD method are very close, with a relative error of 0.66%, which further validates the accuracy of the BPNN. The optimal drag ratio is obtained at a=0.00 and b=1.21, where the follower AUV is directly behind the Leading AUV, and the optimal drag ratio of the AUV fleet is 1.8825, which means that the total drag of the AUV fleet is reduced by approximately 12%.Table 4
Comparisons of the optimal drag ratio of the AUV fleet between CFD and BPNN.
a
b
BPNN
CFD
Relative error (%)
0.00
1.21
1.8949
1.8825
0.66%To further investigate the flow structures of the optimal fleet layout, Figure16 compares the pressure contours around the optimal AUV fleet and that of a single AUV. The most obvious difference between these two cases is the pressure distribution on the tail of the leader AUV. Due to the block effect of the follower AUV, the pressure on the tail of the leader AUV is recovered. This variation can be clearly observed in the pressure distributions in Figure 17.Figure 16
Contours of pressure around the AUV: (a) single AUV and (b) optimal AUV fleet.
(a)
(b)Figure 17
Pressure distributions along the longitudinal direction of the AUV: (a) single AUV and (b) optimal AUV fleet.


### 6.2. Optimization of the Follower AUV

Although the optimal fleet layout analyzed in Section6.2 characterizes with the minimum drag, the drag of the follower AUV is not the optimal value. Among the AUVs in a fleet, it is possible that leader AUV is designed with high propulsive performance and the other follower AUVs obtain lower propulsion costs at the expense of the leader AUV. Therefore, the optimal layout of the fleet for the best performance of the follower AUV should be studied.


#### 6.2.1. Neural Network Training

Using the data in Table3, training of the BPNN is carried out until the network meets the intended target. A comparison between the predicted data and the expected data for the follower AUV is shown in Figure 18. As can be seen from the figure, the training network predicts drag ratios with good agreement with the CFD data.Figure 18
Comparisons of the drag ratio of the follower AUV by BPNN and CFD.


#### 6.2.2. Comparisons of the Optimal Drag between CFD and BPNN

Based on the network after the training process, the optimal output drag ratio of the follower AUV is predicted using the BPNN model. Typically, a CFD simulation for the optimal parameters is carried out for the comparisons between the two methods. The optimal results are shown in Table5. Table 5 suggests that the results of the BPNN and the CFD method are very close, with a relative error of 1.87%, which further validates the accuracy of the BPNN. The optimal drag ratio is obtained at a=0.00 and b=1.21, where the AUV fleet is in the Pull Region. The optimal drag ratio of the follower AUV is 0.3417, which means that the drag of the follower AUV is reduced by approximately 66%.Table 5
Comparisons of the optimal drag ratio of the follower AUV between CFD and BPNN.
a
b
BPNN
CFD
Relative error (%)
0.22
0.29
0.3481
0.3417
1.87%Figure19 compares the pressure contours around the optimal AUV fleet and that of a single AUV. The most obvious difference between these two cases is the pressure distribution around the nose of the follower AUV. At this position, the nose of the follower AUV is influenced by the negative pressure on the inside surface of the leader AUV, so that a smaller pressure is obtained. This can also be observed in the pressure distributions in Figure 20.Figure 19
Contours of pressure around the AUV: (a) single AUV and (b) optimal AUV fleet.
(a)
(b)Figure 20
Pressure distributions along the longitudinal direction of the AUV: (a) single AUV and (b) optimal AUV fleet.


## 6.1. Optimization of the AUV Fleet

When the AUV fleet executes tasks underwater, a minimum drag is expected so that the fleet will work for a longer time and a wider range. Therefore, the layout of the AUV fleet should be optimized to obtain a minimum drag.


### 6.1.1. Neural Network Training

Using the data in Table3, training of the BPNN is carried out until the network meets the intended target. A comparison between the predicted data and the expected data is shown in Figure 15. The training network predicts drag ratios with good agreement with the CFD data. Therefore, it can be said that the BPNN successfully model and predict the drag ratio of the AUV fleet.Figure 15
Comparisons of the drag ratio of the AUV fleet by BPNN and CFD.


### 6.1.2. Comparisons of the Optimal Drag between CFD and BPNN

A CFD simulation of the AUV fleet for the optimal parameters is carried out for the comparisons between the two methods. The optimal results are shown in Table4. The results of the BPNN and the CFD method are very close, with a relative error of 0.66%, which further validates the accuracy of the BPNN. The optimal drag ratio is obtained at a=0.00 and b=1.21, where the follower AUV is directly behind the Leading AUV, and the optimal drag ratio of the AUV fleet is 1.8825, which means that the total drag of the AUV fleet is reduced by approximately 12%.Table 4
Comparisons of the optimal drag ratio of the AUV fleet between CFD and BPNN.
a
b
BPNN
CFD
Relative error (%)
0.00
1.21
1.8949
1.8825
0.66%To further investigate the flow structures of the optimal fleet layout, Figure16 compares the pressure contours around the optimal AUV fleet and that of a single AUV. The most obvious difference between these two cases is the pressure distribution on the tail of the leader AUV. Due to the block effect of the follower AUV, the pressure on the tail of the leader AUV is recovered. This variation can be clearly observed in the pressure distributions in Figure 17.Figure 16
Contours of pressure around the AUV: (a) single AUV and (b) optimal AUV fleet.
(a)
(b)Figure 17
Pressure distributions along the longitudinal direction of the AUV: (a) single AUV and (b) optimal AUV fleet.


## 6.1.1. Neural Network Training

Using the data in Table3, training of the BPNN is carried out until the network meets the intended target. A comparison between the predicted data and the expected data is shown in Figure 15. The training network predicts drag ratios with good agreement with the CFD data. Therefore, it can be said that the BPNN successfully model and predict the drag ratio of the AUV fleet.Figure 15
Comparisons of the drag ratio of the AUV fleet by BPNN and CFD.


## 6.1.2. Comparisons of the Optimal Drag between CFD and BPNN

A CFD simulation of the AUV fleet for the optimal parameters is carried out for the comparisons between the two methods. The optimal results are shown in Table4. The results of the BPNN and the CFD method are very close, with a relative error of 0.66%, which further validates the accuracy of the BPNN. The optimal drag ratio is obtained at a=0.00 and b=1.21, where the follower AUV is directly behind the Leading AUV, and the optimal drag ratio of the AUV fleet is 1.8825, which means that the total drag of the AUV fleet is reduced by approximately 12%.Table 4
Comparisons of the optimal drag ratio of the AUV fleet between CFD and BPNN.
a
b
BPNN
CFD
Relative error (%)
0.00
1.21
1.8949
1.8825
0.66%To further investigate the flow structures of the optimal fleet layout, Figure16 compares the pressure contours around the optimal AUV fleet and that of a single AUV. The most obvious difference between these two cases is the pressure distribution on the tail of the leader AUV. Due to the block effect of the follower AUV, the pressure on the tail of the leader AUV is recovered. This variation can be clearly observed in the pressure distributions in Figure 17.Figure 16
Contours of pressure around the AUV: (a) single AUV and (b) optimal AUV fleet.
(a)
(b)Figure 17
Pressure distributions along the longitudinal direction of the AUV: (a) single AUV and (b) optimal AUV fleet.


## 6.2. Optimization of the Follower AUV

Although the optimal fleet layout analyzed in Section6.2 characterizes with the minimum drag, the drag of the follower AUV is not the optimal value. Among the AUVs in a fleet, it is possible that leader AUV is designed with high propulsive performance and the other follower AUVs obtain lower propulsion costs at the expense of the leader AUV. Therefore, the optimal layout of the fleet for the best performance of the follower AUV should be studied.


### 6.2.1. Neural Network Training

Using the data in Table3, training of the BPNN is carried out until the network meets the intended target. A comparison between the predicted data and the expected data for the follower AUV is shown in Figure 18. As can be seen from the figure, the training network predicts drag ratios with good agreement with the CFD data.Figure 18
Comparisons of the drag ratio of the follower AUV by BPNN and CFD.


### 6.2.2. Comparisons of the Optimal Drag between CFD and BPNN

Based on the network after the training process, the optimal output drag ratio of the follower AUV is predicted using the BPNN model. Typically, a CFD simulation for the optimal parameters is carried out for the comparisons between the two methods. The optimal results are shown in Table5. Table 5 suggests that the results of the BPNN and the CFD method are very close, with a relative error of 1.87%, which further validates the accuracy of the BPNN. The optimal drag ratio is obtained at a=0.00 and b=1.21, where the AUV fleet is in the Pull Region. The optimal drag ratio of the follower AUV is 0.3417, which means that the drag of the follower AUV is reduced by approximately 66%.Table 5
Comparisons of the optimal drag ratio of the follower AUV between CFD and BPNN.
a
b
BPNN
CFD
Relative error (%)
0.22
0.29
0.3481
0.3417
1.87%Figure19 compares the pressure contours around the optimal AUV fleet and that of a single AUV. The most obvious difference between these two cases is the pressure distribution around the nose of the follower AUV. At this position, the nose of the follower AUV is influenced by the negative pressure on the inside surface of the leader AUV, so that a smaller pressure is obtained. This can also be observed in the pressure distributions in Figure 20.Figure 19
Contours of pressure around the AUV: (a) single AUV and (b) optimal AUV fleet.
(a)
(b)Figure 20
Pressure distributions along the longitudinal direction of the AUV: (a) single AUV and (b) optimal AUV fleet.


## 6.2.1. Neural Network Training

Using the data in Table3, training of the BPNN is carried out until the network meets the intended target. A comparison between the predicted data and the expected data for the follower AUV is shown in Figure 18. As can be seen from the figure, the training network predicts drag ratios with good agreement with the CFD data.Figure 18
Comparisons of the drag ratio of the follower AUV by BPNN and CFD.


## 6.2.2. Comparisons of the Optimal Drag between CFD and BPNN

Based on the network after the training process, the optimal output drag ratio of the follower AUV is predicted using the BPNN model. Typically, a CFD simulation for the optimal parameters is carried out for the comparisons between the two methods. The optimal results are shown in Table5. Table 5 suggests that the results of the BPNN and the CFD method are very close, with a relative error of 1.87%, which further validates the accuracy of the BPNN. The optimal drag ratio is obtained at a=0.00 and b=1.21, where the AUV fleet is in the Pull Region. The optimal drag ratio of the follower AUV is 0.3417, which means that the drag of the follower AUV is reduced by approximately 66%.Table 5
Comparisons of the optimal drag ratio of the follower AUV between CFD and BPNN.
a
b
BPNN
CFD
Relative error (%)
0.22
0.29
0.3481
0.3417
1.87%Figure19 compares the pressure contours around the optimal AUV fleet and that of a single AUV. The most obvious difference between these two cases is the pressure distribution around the nose of the follower AUV. At this position, the nose of the follower AUV is influenced by the negative pressure on the inside surface of the leader AUV, so that a smaller pressure is obtained. This can also be observed in the pressure distributions in Figure 20.Figure 19
Contours of pressure around the AUV: (a) single AUV and (b) optimal AUV fleet.
(a)
(b)Figure 20
Pressure distributions along the longitudinal direction of the AUV: (a) single AUV and (b) optimal AUV fleet.


## 7. Conclusions

The drag of the AUV fleet significantly influences the operating time and range of underwater vehicles. In this study, an optimization method is proposed to find the optimal layout of the AUV fleet which has the minimum drag. A combined method of CFD simulation, BPNN, and GA method is utilized for the optimization of the AUV fleet. Important conclusions of this study include the following:(
1
) The CFD results show that the drags of the AUVs in the fleet are determined by the relative position of the two AUVs. The layout of the fleet can be categorized into four different regions based on the drag behavior of the two AUVs: the Parallel Region, the Tandem Region, the Pull Region, and the Push Region. The variation of the drags of the AUVs is the result of the interacting flow and pressure change on the two AUVs.(
2
) The optimal layout parameters for the minimum fleet drag are a=0.00 and b=1.21, where the fleet is in the Tandem Region and predicts a drag approximately 12% lower than that of a single AUV.(
3
) The optimal layout parameters for the minimum drag of the follower AUV are a=0.22 and b=0.29, where the fleet is in the Pull Region and predicts a drag approximately 66% lower than that of a single AUV.



---

*Source: 5769794-2017-12-07.xml*