Next Article in Journal
Outcome Prediction for Patients with Bipolar Disorder Using Prodromal and Onset Data
Next Article in Special Issue
Monitoring of Possible Activities of Yangsan Fault Zone Using GNSS
Previous Article in Journal
Characterization and Development of Gelatin from Cow Bones: Investigation of the Effect of Solvents Used for Soaking Beef Bones
Previous Article in Special Issue
Fusion Swarm-Intelligence-Based Decision Optimization for Energy-Efficient Train-Stopping Schemes
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analysis of a New Twin Hybrid Autonomous Underwater Vehicle

1
Department of Mechanical Engineering, Korea Maritime & Ocean University, Busan 49112, Republic of Korea
2
Interdisciplinary Major of Ocean Renewable Energy Engineering, Korea Maritime & Ocean University, Busan 49112, Republic of Korea
3
Department of Ocean Advanced Materials Convergence Engineering, Korea Maritime & Ocean University, Busan 49112, Republic of Korea
4
Maritime ICT R&D Center, Korea Institute of Ocean Science and Technology, Busan 49111, Republic of Korea
*
Author to whom correspondence should be addressed.
Appl. Sci. 2023, 13(3), 1551; https://doi.org/10.3390/app13031551
Submission received: 1 December 2022 / Revised: 10 January 2023 / Accepted: 20 January 2023 / Published: 25 January 2023
(This article belongs to the Special Issue Selected Papers from IMETI 2021)

Abstract

:
The twin hybrid autonomous underwater vehicle (THAUV) is a novel type of unmanned underwater platform that consists of a twin torpedo-shaped hull and is actuated by two buoyancy engines and two thrusters proposed in this paper. The THAUV was designed to have faster speed generated by the two buoyancy engines and two thrusters. The two buoyancy engines on each hull and the airfoil are mainly responsible for the diving and surfacing motion, and the thrusters drive the THAUV along the horizontal plane. The THAUV is capable of carrying more instrumentation and energy than a conventional hybrid autonomous underwater vehicle (HAUV) with a single buoyancy engine such that the THAUV can perform more exploration tasks and operate for a longer period in a one-time operation. Different from other unmanned underwater vehicles (UUVs) with two airfoils or wings, the THAUV has a single airfoil connecting the twin hull such that it does not require connecting bars and additional airfoils. For this reason, the structure of THAUV is more compact and simpler. In this paper, a new compact THAUV is designed and CFD simulation is used to obtain the hydrodynamic parameters of THAUV operation in water. The motion model of the THAUV is also established and the operating parameters of the THAUV are obtained by simulation.

1. Introduction

Underwater gliders (UG), autonomous underwater vehicles (AUV), and remotely operated vehicles (ROV) [1,2,3] have been developed for decades as major underwater monitoring platforms. Compared with ROV, AUV are untethered, pre-programmed operating platforms, which means they have a wider working envelope and longer mission cycles [4]. UG is a monitoring system driven by buoyancy, which makes UG possess the characteristics of low energy consumption and longer working time. In the literature, [5] studied a vertically aligned two-compartment autonomous underwater vehicle and performed motion simulations. The hybrid autonomous underwater vehicle (HAUV) combines the advantages of low energy consumption of UG and strong mobility of AUV, alternating the working modes of UG and AUV to complete a wide range of underwater data collection and subsea monitoring. In reference [6], a nonlinear kinematic and dynamic model of a hybrid autonomous underwater vehicle (HAUV) is derived for the two modes of locomotion, namely, propelled (AUV) and gliding (UG) modes. The HAUV adopts the UG operation mode in the process of underwater data acquisition, i.e., by adjusting the magnitude and direction of buoyancy to make it alternately operate underwater diving and surfacing. It can also adjust the buoyancy to make it hover horizontally and use thrusters to drive it to complete the acquisition of data such as seafloor topography. After finishing the sampling work in a certain water area, the AUV mode is used to change the course to the next collection region. The THAUV is a kind of HAUV with two buoyant bodies. The fixed wing plays a vital role in the work of the UG and the HAUV; the lift generated by the wing drives the UG and HAUV forward during the dive and surfacing states and CFD numerical calculations are used to study the hydrodynamic characteristics of torpedo-type underwater gliders. [7]. The acquisition of hydrodynamic parameters is crucial to the simulation prediction of the UG and the HAUV. Common methods include the experimental method, parameter identification method, and computational fluid dynamic (CFD) simulation. The hydrodynamic performance and parameters of the underwater glider were calculated using experimental methods as described in reference [8]. In reference [9], the authors described the identification of model parameters for SLOCUM using experimental flight test data. In particular, the focus is on data from a steady linear glide. CFD is widely used for the calculation of hydrodynamic parameters of underwater equipment because of its low cost and fast result acquisition. The application of CFD methods for modeling and hydrodynamic analysis of underwater gliders is described in detail in Reference [10].
The operating modes of this THAUV can be divided into two types according to the operating mechanism, i.e., glide mode and AUV mode. In the glide mode, the net buoyancy and the position of the center of gravity are regulated by the displacement change of the piston, which drives the glider to move downward (gravity over buoyancy) or move upward (buoyancy over gravity) in the vertical plane. In AUV mode, the piston is adjusted to keep the glider in a neutral state, and the thruster is used to drive the glider to move in the horizontal plane for tasks such as target search and seafloor scanning, while the elevator is responsible for adjusting the attitude of the glider to keep it horizontal (the angle of attack is 0). In this paper, the gliding mode is taken as the research objective, and suitable airfoil selection is obtained through simulation of the gliding motion.

2. Structure of THAUV

The twin hybrid autonomous underwater vehicle (THAUV) is a new type of underwater monitoring platform with two buoyant bodies, which allows it to carry more energy and a variety of measurement equipment. The two buoyant bodies are connected by an aluminum plate, which can also be considered as a fixed wing of the THAUV. In order to reduce the resistance to underwater movement, the main body of the THAUV is a torpedo-shaped rotating body. The rotating body is modeled using the MYRING hull profile equations [11], and the overall structural diagram of the THAUV is shown in Figure 1.
The mass of the entire THAUV system consists of two parts: the dynamic mass, including the mass of the piston ( m p = 1.022 ), the mass of the ballast water ( m b ), and the static mass ( m s = 65.33 ), which is the sum of the mass of the roller ( m r = 1.04 ), and the mass of the remaining static accessories ( m s b = 64.29 ). Due to the limited rotation range of the elevator, the effect of the elevator on the floating center and center of gravity is negligible. The mass distribution is illustrated in Figure 2.
The total mass of the THAUV is expressed in Equation (1):
m = m p + m s = m p r + m p l + m r + m s b
The derived equation for the center of gravity is shown in Equation (2):
r C G = m p r r p r + m p l r p l + m s r s m
r s = m r r r + m s b r s b m s
Here r p r and r p l are the positions of the right and left pistons in the body coordinates, respectively ;   r s is the position of static mass (including roller); r r is the position of the roller in the body coordinates; and r s b is the equivalent position of the remaining static accessories in the body coordinates (it is a fixed value).
Since the roller only rotates around the axis of the right body shown in Figure 2 and there is no displacement change in the x-axis direction, the position in the x-axis direction is fixed. Figure 3 shows the schematic diagram after the roller rotation angle γ .
At this time, the center of gravity is shifted in the y-axis direction, resulting in the tilting of the body in the y-axis, and the THAUV will produce a rotary hovering motion around a certain point due to the action of hydrodynamic lift. When the rotation angle of the roller is γ , the position coordinates of the roller can be expressed as:
r r = r r x i 1 0 0 + r r y i + R r s i n γ 0 1 0 + r r z i R r R r c o s γ 0 0 1
Here, R r is the radius of rotation of the roller center of mass and r r x i ,   r r y i , and r r z i are the initial coordinates of the roller in the body coordinates along the x,y,z axis, respectively, when the roller does not rotate.
In this design, the buoyancy engine adopts the structure of piston, which changes the magnitude and direction of the net buoyancy of the THAUV through the movement of the piston to absorb or drain water. When the piston absorbs water inward, the gravity is greater than the buoyancy, and the THAUV is diving; when the piston drains water outward, the buoyancy is greater than the gravity, and the THAUV is surfacing. In order to keep the glider balanced and stable, the ballast water in the two buoyancy engines is always equal in the absence of special requirements. The ballast water is expressed with Equation (5):
m b = ρ π r 2 d p ,
Here, ρ is the density of the water, r is the radius of the piston ( r = 0.05   m ), and   d p is the displacement of the piston along the x-axis. Therefore, the expression of the buoyancy force of the THAUV is shown below:
m B = M B 2 m b
Here, m B is the maximum displacement mass of THAUV.
The buoyancy center is derived with Equation (7):
r C B = m b r p r + r p l m B

3. Mathematic Model

In order to analyze the motion of the THAUV underwater, the coordinate system of the THAUV's motion was established as shown in Figure 4. The coordinate system includes a reference coordinate system E-XYZ and body coordinate system O-xyz, with the buoyancy center of the THAUV as the coordinate origin, in addition to the velocity coordinate system π 0 of the THAUV's motion.
The position and attitude of the THAUV in the inertial coordinate system can be expressed as b = x y z T and Ω = p q r T . The linear and angular velocities of the THAUV in the body coordinate system can be expressed as V = u v w T and ω = p q r T . The rotation matrix from the body frame to the inertial frame is R .
R = c o s ψ cos θ s i n ψ c o s ϕ + c o s ψ s i n θ s i n ϕ s i n ψ s i n ϕ + c o s ψ c o s ϕ s i n θ s i n ψ c o s θ c o s ψ c o s ϕ + s i n ϕ s i n θ s i n ψ c o s ψ s i n ϕ + s i n θ s i n ψ c o s ϕ s i n θ c o s θ s i n ϕ c o s θ c o s ϕ
Then, the kinematic equation of the THAUV can be expressed as:
b ˙ = R V
R ˙ = R ω ^
A detailed derivation of the kinetic equations of the underwater glider is given in the literature [12,13]. Referring to the above literature, the equations of motion governing a THAUV propagating in a three-dimensional space are expressed as:
R ˙ b ˙ v ˙ ω ˙ r ˙ p P p P ˙ p m ˙ 0 = R ω ^ R v M 1 F J 1 τ d ˙ p m p v + ω × r p + r ˙ p u p u 0 d ˙ p
The total moment and force are expressed as follows:
τ = J ω + r ^ p P p × ω + M v × v + ω × r p × P p + m p r ^ p g R T k + τ e x t r ^ p u p
F = M v + P p × ω + m 0 g R T k + F e x t u p
Here, τ e x t includes the hydrodynamic moment and thruster torque, and F e x t includes the hydrodynamic forces and thrust of thrusters.
The THAUV in this paper can not only perform pitching motion in the vertical plane, but it can also change the heading in the horizontal direction by using the thrusters and achieve gyratory motion by controlling the deflection angle of the roller. In this paper, we focus on the pitch motion attitude of the THAUV, where the thruster, roller, and elevator remain stationary during the pitch motion. Then, the equations of the motion of the THAUV in the vertical plane can be simplified as:
x ˙ = v 1 c o s θ + v 3 s i n θ
z ˙ = v 1 s i n θ + v 3 c o s θ
θ ˙ = ω 2
v ˙ 1 = 1 m 1 m 3 v 3 ω 2 P p 3 ω 2 m 0 g s i n θ + L s i n α Dcos α u p 1
v ˙ 3 = 1 m 3 m 1 v 1 ω 2 + P p 1 ω 2 + m 0 g c o s θ L c o s α + Dsin α u p 3
ω 2 ˙ = 1 J 2 m 3 m 1 v 1 v 3 r p 1 P p 1 + r p 3 P p 3 ω 2 m g r r x cos θ + r r z sin θ + M D L r p 3 u p 1 + r P 1 u P 3
r ˙ p 1 = d ˙ p 1
P p 1 = m p v 1 + ω 2 r p 3 + d ˙ p 1
P p 3 = m p v 3 ω 2 r p 1
P ˙ p 1 = u p 1
P ˙ p 3 = u p 3
m ˙ 0 = u 0 d ˙ p 1
α =   acrtan v 3 v 1
ζ =   θ α
Here, x and z are the horizontal and vertical displacements of the THAUV in the inertial frame, respectively; α is the angle of attack; θ is the pitching angle; ζ is the gliding angle; v 1 and   v 3 are the velocity components of the THAUV along the x- and z-axes in the body frame, respectively; ω 2 represents the pitch angular velocity of the THAUV rotating around the y-axis of the body frame; L is the hydrodynamic lift; D is the hydrodynamic drag; M D L represents the pitching moment yielded by the hydrodynamic analysis; m 1 and   m 3 represent the added mass along the x- and z-axes, respectively; J 2 represents the moment of inertia rotating around the y-axis; and P P 1 ,   P P 3 represent the momenta of the piston along the x- and z-axes, respectively.
Table 1 shows a detailed description of the parameters used in the previous equations and Table 2 presents the values of each case.

4. Hydrodynamic Analysis

When underwater equipment works in the water, it is mainly affected by buoyancy and the hydrodynamic force generated by the relative movement of the water around the body in the process of movement. Compared with the buoyancy, the calculation of hydrodynamic force is not easy; however, it is very critical. Currently, three methods are usually used to calculate the hydrodynamic force of water-related equipment, i.e., for the target with a regular shape, the empirical formula is mostly used to predict the hydrodynamic force; for the equipment with a simple structure and low manufacturing cost, the method of the proportional model experiment is often used to obtain the hydrodynamic parameters directly; however, for the equipment with a complex structure and high manufacturing cost, the above method cannot accurately and effectively calculate the hydrodynamic force of the target. For such cases, ANSYS Fluent [14] software is used to calculate the hydrodynamic parameters of the target.
The hydrodynamic drag force of the THAUV was simulated using Fluent in this paper. The simulation domain is a rectangular body and the dimensions of the simulation domain are: 14 m in length, 6m in both width and height, and the inlet of the simulation domain is 5 m away from the floating center of the THAUV. The k-omega SST model [15] is selected as the turbulence model for this simulation and the SST model is widely used for underwater equipment simulation analysis. Due to the complexity of the THAUV tail structure, a non-structural grid is used for hydrodynamic simulation in this paper, with a total grid number of about 3,790,000 and a Yplus less than 1 [16]. The created grid is shown in Figure 5.
In order to ensure that the number of meshes can meet the accuracy of the simulation results, the comparison of the simulation results of different meshes in the state of the THAUV with an angle of attack of +2 degrees and a gliding speed of 0.5 m/s can be obtained from the verification results of the mesh independence [17] shown in Figure 6. The grid numbers used in the figure are 2,250,000, 3,330,000, 3,790,000, and 75,200,000. In view of the simulation results, simulation time, and computer performance, the number of meshes used in this paper is about 3,790,000.
In this paper, the effect of the hydrodynamics of three different sizes of fixed wing on the movement of the THAUV underwater is studied, and the three cases are shown in Figure 7 below.
The hydrodynamic lift, drag, and moment are shown in Figure 8, Figure 9 and Figure 10 below. The trend of the values shown in the graphs shows that when the angle of attack and the velocity increase, the hydrodynamic lift, drag, and moment will also increase; the lift and moment increase almost linearly, while the growth trend of the hydrodynamic drag is in line with the quadratic growth. Since the three cases have similar hydrodynamic trends, only the graphs of hydrodynamic changes for Case 01 are shown here.
From the trend of hydrodynamic forces and moment shown above, it can be seen that the hydrodynamic lift ( L ) and moment ( M D L ) can be expressed as a linear function of the angle of attack ( α ), and the hydrodynamic resistance ( D ) can be expressed as a quadratic function equation of the angle of attack. These expressions are shown as follows:
L = K L 0 + K L α V 2
D = K D 0 + K D α 2 V 2
M D L = K M 0 + K M α V 2
Here, K s are the corresponding hydrodynamic parameters and their values are shown in the following Table 3.
A comparison of hydrodynamic lift and drag relative to cruising velocity of the THAUV when it is in the AUV mode of horizontal cruise is shown in the figure below. Figure 11 and Figure 12 show the relationship between the hydrodynamic lift and drag of the THAUV in the horizontal plane relative to the cruising speed, respectively. It can be concluded from the plot of drag that the wing area has a relatively small effect on drag as it decreases, since the buoyant body is the main factor in drag generation.

5. Results

Based on the modeling of the aforementioned Equations (14)–(27) of motion in Section 3, this chapter mainly shows the simulation results of the THAUV within the vertical water column. The prerequisites for the simulations shown in this paper are that the roller, thruster, and elevator are at rest, mainly to show the effect of the buoyancy engine on the glide angles, and the simulation results can also be used as a reference for the actual underwater experiments. In this paper, the motion of the piston is used as the only control input.
The simulation results of the THAUV for the three cases in two working cycles (dive and then surfacing, then dive and surfacing again) are shown in Figure 13, Figure 14 and Figure 15. During the first stage of diving and surfacing, the displacement of the piston is plus or minus 0.03 m, i.e., the change in the mass of the ballast water is plus or minus 483 g. The positive sign indicates that the buoyancy engine absorbs water inward and the buoyancy of the THAUV decreases. During the second stage of diving and surfacing at the time of 1060 s and 1430s, the displacement of the piston is plus or minus 0.05m. Figure 16 presents the comparison of the angle of attack (AoA) for the three cases.
From the above comparison results, it is clear that the gliding performance of the THAUV is enhanced when the fixed wing size decreases. From the trend of the gliding performance, it can be inferred that the gliding performance will not change any more when the fixed wing size decreases to a specific value. From the curve of glide velocity, we can find that the THAUV reached the stable dive state at around 170s, and the dive velocities of Cases 01, 02, and 03 were 0.25m/s, 0.32m/s, and 0.34m/s, respectively. In the second dive stage of the THAUV, when the THAUV reached the stable dive state, the dive velocity was 0.43m/s, 0.54m/s, and 0.57m/s, respectively. By simple calculation, the gliding speed of Case 02 increased by 25.58% and 28% at the first and second dive states, respectively, compared with Case 01, and the increase became bigger when the control input was larger. After more simulations, the optimal size of the fixed wing can be obtained. It is important to note that the fixed wing in the middle of the THAUV in this paper not only provides lift but also connects two buoyant bodies. While obtaining the maximum gliding performance, it is also necessary to pay attention to the rigid strength of the whole system to avoid the unstable failure of the system in the unstable underwater environment. In the above three cases, the gliding performance of Case 01 is not as good as that of Case 02, although Case 01 meets the system strength requirement. While the gliding performance of Case 03 is better than that of Case 02, the performance improvement is not obvious and the smaller airfoil size does not satisfy the rigid body strength requirement. Therefore, through the comparison of these three cases, the airfoil in Case 02 can not only achieve the best gliding performance but also meet the rigidity requirements of the system.

6. Conclusions

In this paper, a novel twin hybrid autonomous underwater vehicle (THAUV) that consists of a twin torpedo-shaped hull and is actuated by two buoyancy engines and two thrusters is proposed. The THAUV was designed to have faster speed generated by the two buoyancy engines and two thrusters. The two buoyancy engines on each hull and the airfoil are mainly responsible for the diving and surfacing motion, and the thrusters drive the THAUV along the horizontal plane.
In this paper, the modeling and motion simulation of a newly designed THAUV were performed and the gliding parameters of the THAUV in the vertical plane were obtained by simulation under the condition that only the buoyancy engine was operating. The effect of fixed wing size on the gliding performance of the THAUV was analyzed through the simulation comparison of three case designs.
From the cruising drag in the horizontal plane, it can be seen that the size change of the wing has little effect on the hydrodynamic drag compared with the effect of the buoyant body. The gliding parameters of the THAUV are obtained by simulation, and in the comparison of Case 01 with Case 02 it is found that the downward gliding speed of Case 02 increases by about 26% compared with that of Case 01 for the same input. It can be seen that among the three cases, Case 02 showed the best performance in gliding speed performance and in the structural strength requirements of the system.

Author Contributions

Conceptualization, H.-S.C., and J.H.; methodology, J.H.; software, J.H. and P.H.N.A.; validation, J.H., H.C., J.-H.P., K.-B.C., D.-W.J., P.H.N.A., R.Z., J.-Y.K., and D.J.; formal analysis, J.H.; investigation, J.H., J.-Y.K., and D.J.; resources, J.H. and H.-S.C.; data curation, J.H.; writing—original draft preparation, J.H.; writing—review and editing, D.-W.J., H.C., J.-H.P., K.-B.C., R.Z., H.-S.C., J.-Y.K., and D.J.; supervision, H.-S.C.; project administration, H.-S.C.; funding acquisition, H.-S.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the Unmanned Vehicles Core Technology Research and Development Program through the National Research Foundation of Korea (NRF) and Unmanned Vehicle Advanced Research Center (UVARC), funded by the Ministry of Science and ICT, the Republic of Korea, (NRF-2020M3C1C1A02086321) and (NRF-2020M3C1C1A0208632612).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors acknowledge all members of the Korea Maritime University Intelligent Robot & Automation Lab.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Webb, D.C.; Simonetti, P.J.; Jones, C.P. SLOCUM: An underwater glider propelled by environmental energy. IEEE J. Ocean. Eng. 2001, 26, 447–452. [Google Scholar] [CrossRef]
  2. Joung, T.-H.; Sammut, K.; He, F.; Lee, S.-K. Shape optimization of an autonomous underwater vehicle with a ducted propeller using computational fluid dynamics analysis. Int. J. Nav. Archit. Ocean. Eng. 2012, 4, 44–56. [Google Scholar] [CrossRef] [Green Version]
  3. Choi, H.; Chung, S.; Park, H.; Seo, J. Design and control of a convertible ROV. In 2012 Oceans-Yeosu; IEEE: Piscataway, NJ, USA, 2012; pp. 1–7. [Google Scholar] [CrossRef]
  4. Cadena, A.; Teran, P.; Reyes, G.; Lino, J.; Yaselga, V.; Vera, S. Development of a hybrid autonomous underwater vehicle for benthic monitoring. In Proceedings of the 4th International Conference on Control, Automation and Robotics (ICCAR), Auckland, New Zealand, 20–23 April 2018; IEEE: Piscataway, NJ, USA, 2018; pp. 437–440. [Google Scholar] [CrossRef]
  5. He, M.-Q.; Williams, C.D.; Crocker, P.; Shea, D.; Riggs, N.; Bachmayer, R. A simulator developed for a twin-pod AUV, the Marport SQX-500. J. Hydrodyn. 2010, 22, 184–189. [Google Scholar] [CrossRef]
  6. Makdah, A.A.R.A.; Daher, N.; Asmar, D.; Shammas, E. Three-dimensional trajectory tracking of a hybrid autonomous underwater vehicle in the presence of underwater current. Ocean Eng. 2019, 185, 115–132. [Google Scholar] [CrossRef]
  7. Le, T.-L.; Hong, D.-T. Computational fluid dynamics study of the hydrodynamic characteristics of a torpedo-shaped underwater glider. Fluids 2021, 6, 252. [Google Scholar] [CrossRef]
  8. Liu, Y.; Ma, J.; Ma, N.; Huang, Z. Experimental and Numerical Study on Hydrodynamic Performance of an Underwater Glider. Math. Probl. Eng. 2018, 2018, 8474389. [Google Scholar] [CrossRef]
  9. Graver, J.G.; Bachmayer, R.; Leonard, N.E.; Fratantoni, D.M. Underwater glider model parameter identification. In Proceedings of the 13th International Symposium on Unmanned Untethered Submersible Technology, Durham, NH, USA, 24–27 August 2003. [Google Scholar]
  10. Singh, Y.; Bhattacharyya, S.; Idichandy, V. CFD approach to modelling, hydrodynamic analysis and motion characteristics of a laboratory underwater glider with experimental results. J. Ocean. Eng. Sci. 2017, 2, 90–119. [Google Scholar] [CrossRef]
  11. Ji, D.-H.; Choi, H.-S.; Kang, J.-I.; Cho, H.-J.; Joo, M.-G.; Lee, J.-H. Design and control of hybrid underwater glider. Adv. Mech. Eng. 2019, 11, 1687814019848556. [Google Scholar] [CrossRef]
  12. Huang, J.; Choi, H.-S.; Jung, D.-W.; Lee, J.-H.; Kim, M.-J.; Choo, K.-B.; Cho, H.-J.; Jin, H.-S. Design and Motion Simulation of an Underwater Glider in the Vertical Plane. Appl. Sci. 2021, 11, 8212. [Google Scholar] [CrossRef]
  13. Graver, J.G. Underwater Gliders: Dynamics, Control and Design. Ph.D. Thesis, Princeton University, Candidacy, NJ, USA, 2005. [Google Scholar]
  14. Meyers, L.M.; Msomi, V. Hydrodynamic analysis of an underwater glider wing using ANSYS fluent as an investigation tool. Mater. Today Proc. 2021, 45, 5456–5461. [Google Scholar] [CrossRef]
  15. Qin, D.; Huang, Q.; Pan, G.; Shi, Y.; Li, F.; Han, P. Numerical simulation of hydrodynamic and noise characteristics for a blended-wing-body underwater glider. Ocean. Eng. 2022, 252, 111056. [Google Scholar] [CrossRef]
  16. Du, X.; Ali, N.; Zhang, L. Numerical simulations for predicting wave force effects on dynamic and motion characteristics of blended winged-body underwater glider. Ocean. Eng. 2021, 235, 109312. [Google Scholar] [CrossRef]
  17. Javaid, M.Y.; Ovinis, M.; Hashim, F.B.; Maimun, A.; Ahmed, Y.M.; Ullah, B. Effect of wing form on the hydrodynamic characteristics and dynamic stability of an underwater glider. Int. J. Nav. Archit. Ocean. Eng. 2017, 9, 382–389. [Google Scholar] [CrossRef]
Figure 1. The structure of the THAUV.
Figure 1. The structure of the THAUV.
Applsci 13 01551 g001
Figure 2. Mass distribution and descriptions.
Figure 2. Mass distribution and descriptions.
Applsci 13 01551 g002
Figure 3. Schematic diagram of roller rotation.
Figure 3. Schematic diagram of roller rotation.
Applsci 13 01551 g003
Figure 4. Coordinate system of the THAUV.
Figure 4. Coordinate system of the THAUV.
Applsci 13 01551 g004
Figure 5. Non-structural grid of the CFD simulation.
Figure 5. Non-structural grid of the CFD simulation.
Applsci 13 01551 g005
Figure 6. Grid independence verification.
Figure 6. Grid independence verification.
Applsci 13 01551 g006
Figure 7. Three cases with different wing sizes.
Figure 7. Three cases with different wing sizes.
Applsci 13 01551 g007
Figure 8. Hydrodynamic lift.
Figure 8. Hydrodynamic lift.
Applsci 13 01551 g008
Figure 9. Hydrodynamic drag.
Figure 9. Hydrodynamic drag.
Applsci 13 01551 g009
Figure 10. Hydrodynamic moment.
Figure 10. Hydrodynamic moment.
Applsci 13 01551 g010
Figure 11. Variation of hydrodynamic lift with velocity.
Figure 11. Variation of hydrodynamic lift with velocity.
Applsci 13 01551 g011
Figure 12. Variation of hydrodynamic drag with velocity.
Figure 12. Variation of hydrodynamic drag with velocity.
Applsci 13 01551 g012
Figure 13. Comparison of gliding velocity.
Figure 13. Comparison of gliding velocity.
Applsci 13 01551 g013
Figure 14. Comparison of horizontal displacement.
Figure 14. Comparison of horizontal displacement.
Applsci 13 01551 g014
Figure 15. Comparison of diving depth.
Figure 15. Comparison of diving depth.
Applsci 13 01551 g015
Figure 16. Comparison of AOA.
Figure 16. Comparison of AOA.
Applsci 13 01551 g016
Table 1. Descriptions of parameters of the above equations.
Table 1. Descriptions of parameters of the above equations.
TermDescription
x The horizontal displacement in the vertical plane
z The vertical displacement in the vertical plane
v 1 The velocity of the glider along the x-axis in the body-fixed frame
  v 3 The velocity of the glider along the z-axis in the body-fixed frame
θ The angle of pitch
ω 2 The pitch angular velocity rotating around the y-axis of the body-fixed frame
m 1 The mass of the glider (with added mass) in the x-axis
m 3 The mass of the glider (with added mass) in the z-axis
p p 1 The momenta of the piston along the x-axis
p p 3 The momenta of the piston along the z-axis
m 0 Net buoyancy mass
L Hydrodynamic lift
D Hydrodynamic drag
α Angle of attack
u p 1 The force acting on the piston in the x -direction
u p 3 The force acting on the piston in the z -direction
J 2 The moment of inertia rotating around the y-axis (without moveable mass and piston)
r p 1 Position of the piston in the x direction
r p 3 Position of the piston in the z direction
m Mass of the glider
r r x Position of the gravity center in the x direction
r r z Position of the gravity center in the z direction
M D L Hydrodynamic moment
d p 1 The displacement of the piston along the x direction
u 0 The ballast water per meter
ζ Gliding angle
m p Mass of the piston
J The total moment of inertia of the glider rotating around the y-axis
m s Mass of static equilibrium
ρ Density of the water
r Radius of the piston
M B Maximum buoyancy of the glider
Table 2. Values of the THAUV.
Table 2. Values of the THAUV.
TermCase 01Case 02Case 03
m   k g 67.3767.0266.3
m p   k g 1.0221.0221.022
m s   k g 65.3364.9864.256
r   (m)0.050.050.05
Table 3. Values of hydrodynamic parameters.
Table 3. Values of hydrodynamic parameters.
K L 0 K L K D 0 K D K M 0 K M
Case 01−1.91699.528111.11590.2151.1025−1.4179
Case 02−1.84619.835111.55940.24450.3720−2.4335
Case 03−1.2658.094811.28050.205250.0985−2.8938
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Huang, J.; Choi, H.-S.; Jung, D.-W.; Choo, K.-B.; Cho, H.; Anh, P.H.N.; Zhang, R.; Kim, J.-Y.; Ji, D.; Park, J.-H. Analysis of a New Twin Hybrid Autonomous Underwater Vehicle. Appl. Sci. 2023, 13, 1551. https://doi.org/10.3390/app13031551

AMA Style

Huang J, Choi H-S, Jung D-W, Choo K-B, Cho H, Anh PHN, Zhang R, Kim J-Y, Ji D, Park J-H. Analysis of a New Twin Hybrid Autonomous Underwater Vehicle. Applied Sciences. 2023; 13(3):1551. https://doi.org/10.3390/app13031551

Chicago/Turabian Style

Huang, Jiafeng, Hyeung-Sik Choi, Dong-Wook Jung, Ki-Beom Choo, Hyunjoon Cho, Phan Huy Nam Anh, Ruochen Zhang, Joon-Young Kim, Daehyeong Ji, and Jung-Hyeun Park. 2023. "Analysis of a New Twin Hybrid Autonomous Underwater Vehicle" Applied Sciences 13, no. 3: 1551. https://doi.org/10.3390/app13031551

APA Style

Huang, J., Choi, H. -S., Jung, D. -W., Choo, K. -B., Cho, H., Anh, P. H. N., Zhang, R., Kim, J. -Y., Ji, D., & Park, J. -H. (2023). Analysis of a New Twin Hybrid Autonomous Underwater Vehicle. Applied Sciences, 13(3), 1551. https://doi.org/10.3390/app13031551

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop