High Accuracy Buoyancy for Underwater Gliders: The Uncertainty in the Depth Control.

This paper is a section of several preliminary studies of the Underwater Drones Group of the Università degli Studi "Roma Tre" Science Department: We describe the study philosophy, the theoretical technological considerations for sizing and the development of a technological demonstrator of a high accuracy buoyancy and depth control. We develop the main requirements and the boundary conditions that design the buoyancy system and develop the mathematical conditions that define the main parameters.

1. Introduction

This paper is part of several preliminary studies by Underwater Drones Group (UDG) of the Science Department of the Università degli Studi “Roma Tre”, which is developing an advanced Autonomous Underwater Vehicle (AUV) for the exploration of the sea at high depths. The final aim of the project is to create a platform for underwater scientific research that can accommodate a wide range of different payloads.

We will examine the buoyancy system and evaluate its sizing; then we will illustrate the technological solution we have come up with in order to realize the hydraulic system to be assembled in the Underwater Glider Mk. III (see Figure 1) [1,2,3,4,5].

Figure 1

Perspective view of the Underwater Glider Mk. III. (a) Rear/port; (b) front/starboard.

1.1. The Underwater Glider

1.1.1. AUV Evolution

The exploration of the underwater world has always been one of mankind’s dreams: submarines and bathyscaphes (for extreme depths) have been developed to study the “deep blue”. Due to obvious dangers, the human exploration can take place only for very short periods and very limited areas: for these reasons, the exploration of the sea has immediately been drawn towards unmanned automatic systems [6,7,8].

An AUV is a vehicle that travels underwater without requiring input from an operator; this means that it must be equipped with a “brain” that regulates and coordinates its position, its depth and its speed: moreover, it is able to collect and store data from the payload. One of the first realizations was the Autonomous LAgrangian Circulation Explorer (ALACE) system, a buoy that was able to vary its buoyancy and therefore its depth. Although it possessed a great endurance, it only could be employed for great depths and in open sea—the consequences of these limitations are evident.

The next step was the use of Remote Operated Vehicles (ROVs). These, thanks to the constant development of electronic miniaturization, are extremely high performing vehicles for short-lasting marine operations, but the require the constant presence of a support vessel.

The need to get rid of the randomness of the currents has led to the natural development of the underwater glider concept [9,10,11,12,13,14,15].

1.1.2. The Underwater Glider

An underwater glider is a vehicle that, by changing its buoyancy, moves up and down in the ocean like a profiling float [16]. It uses hydrodynamic wings to convert vertical motion into horizontal motion, moving forward with very low power consumption [17,18,19,20,21,22]. While not as fast as conventional AUVs, the glider, using buoyancy-based propulsion, offers increased range and endurance compared to motor-driven vehicles and missions may extend to months and to several thousands of kilometres in range. An underwater glider follows an up-and-down, sawtooth-like mission profile providing data on temporal and spatial scales unavailable with previous types of AUVs [23,24,25,26,27].

1.1.3. The Mk. III Architecture

The Mk. III sub-glider has a cylindrical fuselage with a radome on the bow containing the customizable payload and, on the other end, the hydrodynamic fairing. The vehicle does not have moving surfaces: control is provided by the displacement of the battery package that varies the position of the centre of mass. The wings aerofoil is based on the Eppler E838 Hydrofoil. The aerofoil has the maximum thickness 18.4% at 37.2% chord and maximum camber 0% at 46.5% chord. The arrangement of the internal sectors is visible in Figure 2a,b. The buoyancy system is contained in the buoyancy control bay: it accommodates the buoyancy motor and the oil tank and provides longitudinal balance to the system by adjusting the level in the reservoir. The bladder is contained in the hydrodynamic fairing, in contact with the open water. The fairing is not a critical structural part—it has the task of not disturbing the hydrodynamic flow of the fuselage [28].

Figure 2

Underwater Glider Mk. III cutaway: (a) fuselage prospective section; (b) fuselage sagittal section.

1.1.4. Conventions

We introduce, for clarity, the mathematical conventions and symbols that will be used in the subsequent discussion (see Figure 3a,b):

Figure 3

The Euler angles. (a) Body frame (blue) and reference frame (red); (b) The body frame referred to the drone.

where:

α (or φ) is the angle between the x axis and the N axis.

β (or θ) is the angle between the z axis and the Z axis.

γ (or ψ) is the angle between the N axis and the X axis.

For the rotation matrix, we have: where is the parameter vector [29].

2. Materials and Methods

2.1. The Buoyancy System

2.1.1. Basic Concepts

Gliders are controlled through hydrostatics (vertical forces) and manipulate hydrostatic balances in order to accomplish roll and pitch of the vehicle. Stability of the vehicle is a major critical factor: a stable vehicle has the centre of gravity below the centre of buoyancy. In this configuration, the weight of the vehicle creates a restoring moment to add stability to the vehicle. Roll and pitch on the glider is accomplished by moving the battery pack. Figure 4 below displays a basic concept of a buoyancy system for the glider [30,31,32,33,34,35,36].

Figure 4

Basic scheme of the buoyancy system.

The system is extremely simple: while descending, hydraulic fluid moves from the external inflatable bladder, which produces a high pressure in the internal reservoir, which is at a low pressure through a valve: the decrease in volume of the bladder creates an increase in density, causing negative buoyancy [37,38,39,40,41,42,43,44].

While ascending, hydraulic fluid moves from internal accumulator to the external inflatable bladder through the pump. The increase in volume creates a decrease in density causing positive buoyancy. The seawater also flushes out the open hydrodynamic fairing of the vehicle, aiding it to rise to the surface. For neutral buoyancy, the vehicle must have a density equal to seawater [45,46,47,48,49,50,51,52].

2.1.2. The System Prototype

Our group has developed a technology demonstrator (see Figure 4) of the buoyancy system to validate the related technology and then to test it. To reduce the force required to actuate the oil piston, which pushes the oil in the bladder at high pressure, is necessary to reduce the piston surface (diameter) and increase the stroke: so, the buoyancy engine resembles a “shotgun”. An open-loop stepper motor was used to drive the screw inside the actuator that, in turn, pushes the piston. Two solenoid valves regulate the flow of oil into the bladder [53,54,55,56,57].

The first problem was the occurrence of actuator buckling: under the push of the engine, the probability of a part bending is high, thus deforming the thread and jamming the mechanism. The problem was solved by constructing a rigid cage with four struts that support the piston’s push load, leaving the screw only with the rolling friction load. In the early project development stages, the workgroup was oriented to use a centrifugal pump for all drives: this technology however did not allow us to create strong pressure differences; the need to use a more powerful engine was also highlighted because the prevalence was too low: this would have led to an excessive battery consumption. The second solution was to use volumetric pumps in order to obtain greater differences in pressure (even ones considerably higher than needed). Unfortunately, these would require too much power and are too heavy for our small vehicle [58,59,60,61,62].

At this stage of development, we have also thought to use the oil tank only as a passive fluid reservoir and trimmable counterweight of the payload and as an active actuator for longitudinal stability. The long travel of the piston (forward or backward) ensures the necessary variation of the bladder volume for manoeuvrability [63].

Now, the system prototype seems to work correctly and promises new developments: it is under in several fatigue cycle trials in order to investigate the parts more prone to failure as a sort of burn-in test. Our prototype (breadboard) is presented in Figure 5.

Figure 5

Buoyancy engine prototype.

2.2. Buoyancy

Archimedes’ principle is the main concept underlying the buoyancy of underwater vehicles. When a vehicle is submerged in water, a buoyant force acts on the body vertically upward due to the pressure forces below the submerged body being greater than the pressure forces above. The buoyant force results in a value equal to the weight it displaces [64].

2.2.1. Static Buoyancy

Now, consider our drone (see Figure 6): it is in steady state.

Figure 6

Drone in buoyancy Balance.

In these conditions, the total weight W of the glider is given by: where:

= Net total “weight” in the water.

= Dry Weight of the glider.

= Buoyancy of the oil bladder.

= Buoyancy of the naked glider.

The expression of dry weight is: where:

= Weight of the battery pack.

= Weight of the oil tank.

= Weight of the naked glider (without oil tank and batteries). For “dry weight”, we mean the weight of the vehicle out of the sea, without the hydrostatic force.

So Equation (2) becomes:

This series of equations will be useful later to establish the drone descent attitude.

2.2.2. Dynamic Balance on the Vertical Plane

Here, the drone dive (or emersion) is examined at constant speed: it is in steady-state gliding, the geometry of total forces is explained in Figure 7a. Figure 7b shows the Eppler E838 characteristic “Cl vs. Alpha (=angle of attack )”.

Figure 7

Dynamic balance of the forces (a) Drone geometrical balance of the forces; (b) Cl vs. Alpha (=angle of attack, diagram for Eppler E838 aerofoil.

At equilibrium, for the dynamics on the vertical plane at constant speed we have:

The expression for the lift is:

According to the Eppler E838 characteristic “Cl vs. Alpha” (Figure 7b) when the angle of attack = 0° (is null) the is zero so that the lift force is null. This shows that the drone cannot progress horizontally at constant speed (straight and level): the only mission profile allowed is a sawtooth curve [65].

2.2.3. Glider Trajectory

Because its motion is due to the difference between the forces of weight and buoyancy, the glider is unable to proceed straight and level, thus being forced to follow a dive/climb trajectory made smooth by the wings. Moreover, unlike gliders in air, AUVs can have ascending glide slopes if the net buoyancy is positive, producing a negative sink rate.

The buoyancy engine of the glider allows changing its net buoyancy into alternating positive and negative states, thereby imparting it with the ability to string together a succession of descending and ascending glide slopes referred to as a sawtooth glide.

The behaviour of the vehicle is considered in case of a simple glide slope (refer to Figure 7a):

where:

: is the dynamic pressure.

: is the characteristic area.

: is the mean aerodynamic chord.

: is the angle of attack.

For the other coefficients, we have:

The coefficient of drag is composed of two members: the first is insensitive to the angle of attack and is constant; the second one () is instead a function of the square of the angle. Note that the zero lift coefficient because the wing profile that has been chosen for our project is symmetrical (type Eppler 883). Now is necessary to separate the contributions of the fuselage (body) and of the wings, for the three factors of lift, friction and pitching moment; so the expression is:

According to the Navier-Stokes (approximated) equations and the simplifications above cited, the previous system of equations becomes: where:

is the wing area.

is the distance between the mean aerodynamic and the center of buoyancy.

is a non-dimensional coefficient.

This parameter is necessary to know the exact attitude and therefore the to obtain a constant descent profile [66].

2.2.4. Gliding Forces

In order to allow the vehicle to glide, it is necessary to create a differential buoyancy force and therefore the balance is not null: from the buoyant force expression, the change in volume needed for the buoyancy engine for a full dive is calculated as:

In which is the buoyancy force due to the difference of volume of the bladder: where:

= Seawater density (average 1.025 kg/l).

Gravity approx. to 9.81 m/s2.

Volume difference of the bladder.

The relationship between difference volume of the bladder and buoyancy force is:

2.2.5. Restoring Moment on the Vertical Plan

The position of the center of mass is given by: where:

: is the density.

: is the considered volume.

r: is the distance considered from the reference frame.

The gliders masses define above, is:

The vectorial balance of the forces becomes: where

is the volume of the drone (submerged).

is the rotation function around the reference frame.

is the upward direction.

When the geometrical centre of the body is offset from the CG frame, the resulting torque is given by:

The expression of is:

The expression of is:

Therefore, we have:

According to which, when the body frame is coincident with the centre of the figure:

Moreover, the resulting balancing torque is

This parameter is necessary to know the exact torque and therefore the position to which the servomechanism to move the battery pack will have to obtain a constant descent profile [67].

2.2.6. Uncertainty in the Depth Control

In this part, the problem of uncertainty in controlling the depth of the drone is examined. In the case where the drone has a very low vertical speed, due to the high viscosity of the water, the vehicle stops (and therefore stabilizes itself in the vertical plan) in a few centimeters of water, so this is not a problem.

Consider the case in which the drone is gliding, at a constant speed (): at a depth of +5 m compared to the “target depth”, the bladder swells up to make the drone assume a neutral buoyancy. From a dynamic point of view, our simulations show us that the drone will behave like a mass-spring-damper system, whose temporal behavior is described in Equation (23) and visible in Figure 8 (blue line).

Figure 8

Drone depth dynamic behavior: in the blue line the normal (uncompensated) damping; the azure line the compensated behavior; dotted lines are the upper and lower limit. The red line (out of scale) is the bladder’s “hydrostatic chute” action.

The maximum precision in the measurement of depth is not given by the precision of the instrument, in this case a piezoelectric depth meter can be accurate up to 2.5 cm, but in the architecture of the drone. Because of all the possible attitudes of the vehicle, it is not known if the transducer is placed higher or lower than the centre of gravity: so, unless more than one transducer is installed (which is impractical) the only consideration is that the maximum distance at which the transducer should be placed from the centre of gravity is precisely 49.7 cm, so, the maximum range of precision obtainable from the instrument is:

As shown in Figure 8 (blue line) the natural dynamic stability of the vehicle, i.e., the oscillation within the error band, is obtained after a time which is not acceptable (see Equation (25)). It is not acceptable for two reasons: first of all because of the long damping time of the oscillation; secondly for the overshoot of the “target depth”: the depth limit may have been placed not only by the type of mission but also by the nature (orography) of the seabed. Therefore, such behavior could lead to a collision of the vehicle with the seabed itself or a known obstacle, thus leading to damage:

The solution is to use the bladder as a “hydrostatic parachute” that changes the damping conditions of the system. In this case the simulation have evidenced that, if you place the bladder at maximum buoyancy (Figure 8—red line) shortly after reaching “+5” depth, the resulting behavior of the vehicle is visible in Figure 8 (azure line).

In this case the tolerance band is reached in s (see Equation (26)), half the previous time and there is no danger of overshooting the depth, thus keeping us always in safe conditions:

From a mathematical point of view, the system changes the damping factor, increasing it considerably in fact, . The model parameters are shown in Table 1 below.

Table 1

Model parameters.

f_n	ζ_normal	ζ_compensated	Z
1	0.0794	0.15779	4.972

2.2.7. Dynamic Simulation

In the section, we simulate the behavior of the Eppler 838 aerofoil and the fuselage to evaluate, before a detailed and expensive but more sophisticated 3D simulation, if the proportions and dimensions of the drone fall within the range of measures evaluated as a requirement. For this purpose, the program “JavaFoil—Analysis of Airfoils” [68] is used. JavaFoil is a program which uses several traditional methods for aerofoil analysis. The backbone of the program consists of two methods:

The evaluation of the potential flow. The analysis is done with a higher order panel method (linear varying vorticity distribution) and it calculates the local, inviscid flow velocity along the surface of the aerofoil for any desired angle of attack.

The evaluation of the boundary layer. The analysis is steps along the upper and the lower surfaces of the aerofoil, starting at the stagnation point, solves a set of differential equations to find the various boundary layer parameters, according to the integral method.

The equations and criteria for transition and separation are based on the procedures described by Eppler. A standard compressibility correction according to Karman and Tsien has been implemented to take moderate Mach number effects into account. Usually this means Mach numbers between zero and 0.5 [69,70,71,72].

The simulation shows the critical parameters for the hydrodynamic behaviour of the model: the results are given in Table 2 below [73,74,75,76,77,78].

Table 2

Body and wing critical parameters [Speed = 0.52 kts, Re = 104 and ].

CLb_α	CLw_α	CDb_0	CDw_0	CDb_α	CDw_α	Cmb_α	c_m
0.3997	0.4085	0.01746	0.01676	0.00018	0.001515	0.005	0.4752

The trend of the flow lines with ° around the body and the wing are visible in Figure 9a,b, below.

Figure 9

Flow field and pressure gradient of conditions: Speed = 0.52 kts, Re = 104 and , (a) Body profile of the drone; (b) Eppler 838 Aerofoil.

For completeness, in order to be able to express the induced resistance and the stall delay, the behavior of the aerofoil at different angles of attack is provided in Figure 10.

Figure 10

Flow field and pressure gradient of conditions (Speed = 0.52 kts, Re = 104) are shown at the following : (a) 0.5° (b) 4.0° (c) 15.0°.

3. Conclusions

This paper reports part of several preliminary studies of the Underwater Drones Group of the Università degli Studi “Roma Tre” Science Department and follows the route traced on several conference papers presented at the IEEE International Workshop on Metrology for the Sea (MetroSea); this part is dedicated to the design and engineering study of the part relating to the UAV buoyancy.

This paper highlights the large series of considerations and structural dimensions, going down in great detail, of the Underwater Glider Mk. III that is currently in an advanced development phase. The real novelties of this work are highlighted in the development due to two strong constraints that our group inserted during the design of this AUV: the first is to evaluate the project always under the most conservative (pessimistic) operating conditions; the second is to evaluate how any changes made to the subsystem in development (in our case the buoyancy s/s) is reflected (and forced) on all the other parts (or subsystems). All the results both from the partial simulations and from the construction and testing of subsystems will then be used in the operating vehicle.

The first section is dedicated to the design and engineering study of the part relating to the UAV buoyancy. In the first section, the architecture, the internal arrangement of the sub glider, the type of mission profile and the maximum requirements for the performances are broadly described. In the second part, the buoyancy system is described from an engineering-construction point of view: the solutions developed and implemented in a working prototype were illustrated.

The last part describes the mathematical requirements for sizing the vehicle. Firstly, the static requirements that are used to determine the mechanical and dimensional sizing of the buoyancy engine are examined. Then the dynamic stability (on the vertical plane) of the vehicle is analysed: this quantifies the forces involved during the “glide”. The trajectory is analysed to decide the attitude and the angle of attack: the latter is necessary, in stationary conditions, to determine the work point of the profile or the position of the profile in the diagram “Cl vs. Alpha”. In order to have a constant attitude, it is necessary to balance the moments in the vertical plane so that, once the wing profile is “started” and in progressive acceleration (i.e., while is nearly to reach the terminal velocity) the pitch up effects must be compensated by the movement of the centre of gravity. An evaluation of the uncertainty in the depth control is also provided.

Lastly, a simplified simulation is introduced in order to observe the hydrodynamic behaviour of the fuselage (limited to the profile) and of the aerofoil at different angles of attack, to highlight the stall characteristics.