1. Introduction
These days, water-powered aerial systems are becoming popular. Especially in hydroflight sports, it is very likely that readers have seen spectacular performances performed by these systems either online or in person. Structurally, each of these systems features three main components: a jet ski that pumps water, a water-conveying hose, and a head assembly that carries a rider. In general, their operation basics is that the water jetted from the nozzles underneath the head assembly generates thrusts and lifts off. The rider controls the system either by leaning their body (as is the case for flyboards [
1,
2] and hoverboards [
3]) or by rotating the direction of the water nozzle (as is the case for jet packs [
4] and jet bikes [
3,
5]) to ride. Thereby, the rider can perform complex movements, such as turning, flipping, or spinning. With this maneuverability, water-powered systems have the potential for many aerial operations. However, their inevitable drawbacks make them currently only applicable to hydroflight sports or show performances. Firstly, the systems require a substantial amount of water for propulsion. Secondly, their flight distance is limited by the length of the water-conveying hose. Especially, most of the water-powered systems are manual or semiautomated such that only trained and skilled people are able to operate them. Newly developed products can operate autonomously, thanks to the use of sensors, actuators, and automatic controllers. For example, the Zapata Aquadrone [
3] is equipped with an electric pump and four electric motors rotating four thrust nozzles for self-manipulation. Nevertheless, it is designed primarily for demonstration purposes.
Water-powered aerial systems have shown potential in firefighting missions, particularly in scenarios where fires occur in water-based environments, such as crowded harbors, bridges with vehicle traffic, or boats in the middle of the sea. These scenarios present a significant challenge due to limited accessibility and potential safety hazards. For example, the current method used to put out a boat fire on the sea is demonstrated in
Figure 1. It involves the use of one or more vessels to approach the scene and spray water onto the burning boat, which can be inaccurate and potentially harmful to property and people still on the scene. Moreover, terrestrial robots are not deployable, and the firefighting methods of using commercial aerial vehicles, such as helicopters and fixed-wing aircraft, also face similar challenges in these scenarios. Nevertheless, small-scale aerial systems offer several advantages, including aerial perspective, remote control capabilities, and access to difficult-to-reach areas. Hence, drones can carry fire retardants [
6,
7,
8,
9] or attach with water-conveying hoses and fire nozzles [
10,
11,
12,
13] to battle the fire. Cooperative use of multiple drones, as in [
14,
15], can cover a large object. However, their adoption is still limited due to multiple constraints, such as their limited load-carrying capacity, the potential for generated airflow to exacerbate the fire, and the presence of dangerous high-speed rotating blades.
Meanwhile, in these scenarios, the water supply concern for water-powered systems is no longer a limiting factor. The systems can transport to get close to the accident areas, and their head assembly only needs to take off from there, thus overcoming their limitation on flight distance. By flying above the flames and taking advantage of high ground, water-powered systems can effectively extinguish fires. One demonstration is that firefighters in Dubai in 2007 used a flyboard to bring them to the height needed to extinguish fires on a bridge over the sea [
16,
17]. However, firefighters need skills to operate these because of the manual systems. Additionally, firefighting at a close distance puts a threat to firefighters’ lives. In [
18,
19,
20,
21], Dragon Firefighter, an aerial robot with a water-conveying hose and water-jet propulsion, was developed for firefighting. The robot features four rotatable nozzles that are controlled to maneuver the robot itself and suppress the fire. Experimental studies have demonstrated the robot’s successful and safe suppression of flames in some simple scenarios but have also revealed its limitations in reachability and maneuverability. Therefore, the development of a different method for actuation and control is needed to enable the autonomous operation of water-powered systems.
Three different actuating methods have been proposed in the literature for controlling the motion of water-powered aerial systems. The first method involves altering the weight distribution of the head assembly in order to change its attitude angles, such as the weight-shifting mechanism proposed by Tri et al. [
22]. This mimics the way a person rides his/her flyboard or hoverboard. Unfortunately, this method has some drawbacks, including an increase in weight, nonminimum phase phenomena, and uncontrollable yaw rotation. The second, and the most common actuating method, is to rotate the thrust nozzles only instead of the whole head assembly to change the direction of the thrust force generated by each nozzle. Liu et al. [
23] used two actuators to drive a nozzle rotating 360 degrees, two nozzles in total, for their unmanned water-powered aerial vehicle. The flying robots with water jets in [
24,
25,
26,
27,
28] used a module consisting of two fixed nozzles and two biaxial active nozzles. Additionally, the Zapata Aquadrone [
3] is equipped with four rotatable thrust nozzles. These systems are fully controllable in three-dimensional space but come with several drawbacks, such as nonminimum phase characteristics [
23], very limited operational space in comparison to firefighting scenarios [
24,
25,
26,
27,
28], or complex structure and high requirement of waterproofing [
3], etc. Finally, the last method is to regulate the water flow rate out of the thrust nozzles, for example, the flow-regulated, unmanned, water-powered, aerial vehicle in [
23]. Thereby, the amplitude of the thrust forces is controllable. This mechanism is simpler and easier to control in comparison to the two other methods. The available design in [
23] is also not possible to control the yaw rotation though.
In terms of control strategy, several linear ones have been proposed to achieve the desired motion of the water-powered aerial systems, for example, proportional control with speed feedback in [
24,
25,
26,
27,
28], proportional–derivative control in [
23], and linear quadratic integrator in [
22]. These control systems are simple and only require feedback from the motion of the head assembly. However, simulations and experiments show performance is poor and greatly influenced by the water-conveying hose movement. Controllers similar to those for multi-propeller drones have also been suggested, for example, cascade proportional–integral–derivative control and cascade sliding mode control [
29]. Furthermore, numerous advanced control schemes have been developed in the literature and considered different aspects of aerial systems, such as decoupled control [
30] and cooperative control of multiple systems [
31,
32]. However, it is worth noting that the hose dynamics are not fully represented in the model of the water-powered systems in the mentioned studies, which does not help validate the control performance and preserve control system stability. For instance, [
22,
29] only considered that the gravitation of the hose varies with the altitude of the system. The force from the hose was formulated via a heuristic assumption in [
23]. In [
26,
27,
28], the hose was approximated by multiple links connected by elastic joints and dampers.
The hose in the abovementioned systems can be considered as a cantilevered pipe that conveys fluid under deformation. Hence, a dynamical model of it, and subsequently of the whole system, can be obtained. Theoretical models of cantilevered pipes were introduced in [
33,
34,
35] in the form of nonlinear partial differential equations. Although the resulting models are typically complex, they can be used to determine the pipe’s deformation and vibration characteristics. Thereby, a control system that takes the partial states of the hose into consideration is able to achieve effective and robust motion performances. He et al. [
36,
37,
38,
39,
40] applied this principle to stabilize flexible marine risers, which share some characteristics with a water-conveying hose. In detail, their proposed controllers need to know up to third-order partial derivatives of the riser at the control position. However, these values are challenging to obtain in practice. Alternatively, the conditions for the hose stability can be deduced from its model, allowing a control system to be designed with only the feedback from the head assembly while still preserving the stability of the whole system. Ambe et al. [
24] derived the model of their complete water-jet system in the form of partial differential equations, but only for small deformation of the hose in the lateral direction. The authors proposed a proportional–derivative controller for the motion of the head and ensured energy dissipation of the system. In [
41], it is shown that if the rate of water flow inside the hose is below a critical value, the water-conveying hose can be kept stable, and the total effect of the hose on the head assembly is bounded. Hence, a state-feedback controller that considers the saturation of the control inputs can preserve the system’s stability. Unfortunately, the designed control law is still insufficient in dealing with disturbances effectively.
Therefore, in this study, we first introduce a water-powered aerial system with a novel flow-regulating mechanism that makes it fully controllable in three-dimensional space. The system is aimed at firefighting missions, especially at or in bodies of water. It also can be extended for land missions by using the water source from the available fire hydrant system. The dynamics of the designed system are derived in the form of coupled partial ordinary differential equations. We then proposed a composite controller that ensures the desired altitude and attitude of the head while still preserving the overall stabilization of the whole system. Distinct characteristics of the system, such as input coupling, saturation, and motion of the hose, were considered in the control system design process. The result is that uniform ultimate boundedness is achieved. Comparative studies with two previous control strategies validated the efficiency and robustness of the proposed system. Hence, the contributions of the paper can be summarized as follows:
A novel water-powered aerial system with a flow-regulating mechanism is introduced and modeled.
A composite controller is proposed to fully drive the system in three-dimensional space and preserve the system’s stability. The uniform ultimate boundedness of the control system is proven.
Computational studies in comparison with two previous control strategies validated the efficiency and robustness of the proposed control system.
Accordingly, the remainder of the paper is organized as follows.
Section 2 describes the design and operation of the system. Mathematical models of the water-conveying hose, the head assembly, and the whole system are derived in
Section 3. The proposed controller is presented in
Section 4. In
Section 5, the comparative results are shown and discussed. Finally, conclusions are drawn in
Section 6.
2. System Description
As indicated by its name, the proposed system uses water as a source of propulsion and firefighting. At or in bodies of water, the enormous amount of water is used. In other fire sites, water can be taken from fire hydrants. The proposed system utilizes a water pump that pulls water through a flexible hose and delivers it to a head assembly, where it flows through a manifold before being distributed to nozzles for generating thrust and performing maneuvering and firefighting actions. The head assembly, which contains major components of the system, such as nozzles, landing gear, instruments, actuators, etc., connects to the hose via a swivel ball joint to isolate the torques from the hose. The head is designed such that the water inlet is located at its mass center, while the water-conveying hose has a uniform cross-section to ensure a consistent flow rate at any point. The design of the water-powered firefighting system is illustrated in
Figure 2a.
On the head assembly, there are four nozzles for generating thrust. Each nozzle has an outlet with a much smaller cross-sectional area than the head inlet. This structure accelerates the water flow, hence increasing its momentum and generating thrust. The four nozzles are numbered 1 to 4, with nozzle 1 being located at the front left and the remaining nozzles following in a counterclockwise direction. Each nozzle injects water downward and slightly out to the lateral sides, as shown in the figure. This arrangement gives several advantages over the other configuration in [
1,
3,
4,
5,
23,
42], including increased flight stability, ability to achieve yaw rotation, and decoupling of control input in the longitudinal and pitch directions. The system also features another nozzle located at the front of the head that sprays the water forward in fine droplets to cover and suppress a flame. The posture of this fire nozzle is slightly downward, and its axial line goes through the head’s center to prevent unwanted torque generation.
Three motorized valves, labeled as
a,
b, and
c, regulate the water flow distribution among the four thrust nozzles. An additional on/off valve, labeled as valve
d, is used to fully open or close the fire hose nozzle. Apparently, the fire hose nozzle only needs to open after approaching the fire. The corresponding flow diagram is shown in
Figure 2b. At the neutral position of the motorized valves, the water flow is equally distributed among the four thrust nozzles. Manipulating one or more of these valves results in a difference in the flow rates and thus creates torques. For instance, opening valve
a toward valve
b increases the water flow rate and thrust of nozzles 1 and 3 and decreases those of nozzles 2 and 4. This results in a torque around the vertical axis passing through the center of the head, thereby producing a yaw rotation.
3. System Modeling
In order to realize the maneuverability of the proposed system, its mathematical model is derived in this section. Initially, a system of coordinates is defined, and the kinematic relationships among them are derived. Next, the nonlinear dynamics of the water-conveying hose and the head part are considered. Finally, a governor equation with initial and boundary conditions represents the completed model of the system.
3.1. Coordinates and System Kinematics
Two Cartesian coordinates and a curvilinear one are defined regarding the motion of the water hose and the system head. One Cartesian coordinate is the Earth-fixed reference frame
XYZ, and the other is the body-fixed frame of the head assembly
XbYbZb. The latter is placed at the mass center of the head with the water inside it. The
Xb-,
Yb-, and
Zb-axes follow the longitudinal, lateral, and vertical of the head, with the
Xb-axis pointing forward and the
Zb-axis pointing upward. The curvilinear coordinate starts at the pump outlet and is along the center line of the water-conveying hose. This is depicted in
Figure 3. Forces acting on the head part are also shown in this figure.
In the reference frame, the hose element at a coordinate
s at a time
t can be described by the partial vector
. Let the dots be denoted for time derivatives, and the primes stand for derivatives with respect to coordinate
s. Then, the velocity of a hose element is
, and the velocity of the water flowing through that is computed by:
where
is the axial velocity of the water flow in the hose.
It is noted that
defines the position of the hose inlet, and
, with
L being the length of the hose, is for the hose outlet. In other words,
is also the position of the head assembly. Additionally, the completed motion of the head can be represented by:
In which:
: vector of translational velocity with respect to the body-fixed frame.
: vector of three Euler angles in XYZ order.
: vector of angular velocities.
, with and , are the transformation matrix from the head-fixed frame to the reference frame. is that from the angular velocities to the Euler angles’ rate.
: m-by-n zero matrix.
For the water flow at an inlet/outlet of the head assembly, its velocity can be obtained as follows:
with
indicating the head inlet, the outlet of the fire hose nozzle, and the outlet of the
ith thrust nozzle, respectively. The second subscript
b means that the value is in the body-fixed frame of the head assembly.
is the position vector of the
j inlet/outlet, and
is the unit-length normal vector of the corresponding cross-sectional area.
Remark 1: is the vector tangent with the hose center line at coordinate s, so it is also the normal vector of the cross-section at this location. Subsequently, .
3.2. Water-Conveying Hose Model
In order to realize the dynamics of the water-conveying hose, let the hose be isolated from the head assembly and replaced by corresponding reaction forces. Then, it can be considered that the hose has a free end. A modified version of Hamilton’s principle for the assembly of the hose and the water inside is written as follows:
where
is the Hamilton’s variational operator.
is the virtual work performed by external forces on the hose. They are listed as follows:
In which , with being a matrix of positive definite coefficients, is the damping on the hose. and are the reaction forces due to isolating the head and the total force holding the hose’s first end. is the total interactive force between the hose and its surroundings. For the sake of simplicity, complex aerodynamic forces, such as drag, wind, ground effect, etc., were neglected. Hence, contains the contact forces (normal and frictional forces) with the ground in onshore applications and includes buoyance force, hydrodynamic force, and seabed interaction force in offshore applications. is the magnitude of the force due to the exit pressure at the hose outlet. The Dirac delta functions indicate the locations that the forces act on, i.e., .
The second term in Equation (4) is the Lagrangian including kinetic energy
T and potential energy thanks to the hose’s bending strain
B and gravitation
G. That is:
In detail, they are obtained as follows:
where
and
are the masses per unit length of the hose and the water inside, respectively.
is the flexural rigidity of the hose, and
is the curvature of the deformed centerline.
g is the gravitational acceleration.
The right-hand term in Equation (4) is the virtual momentum transport across the open surface at the end of the hose. Additionally, we have the following relations with the mass flow rate of the water
inside the hose:
with
being the water density and
being the hose’s cross-sectional area.
Substitute Equations (5)–(8) into Equation (4) and neglect the acceleration/deceleration of the water flow inside the hose. After a complicated but straightforward calculation, the governor equation of the water-conveying hose is obtained:
with the boundary and initial conditions given by:
Remark 2: The three-dimensional dynamics (9) of the water-conveying hose can be reduced to two partial motions thanks to the inextensible condition of the hose as follows [33,34,35]: Remark 3: Consider the system operating in onshore missions, and the interactive force consists of the normal contact force with the ground and frictional forces generated by the motion of the hose on the ground surface. They can be modeled by: with
being the frictional coefficient and
being the signum functional.
3.3. Head Assembly Model
For the sake of modeling, the head assembly with its contained water was treated as a rigid body. This was performed thanks to the isolation of the hose and the nondeforming control volume from the water inlet to the nozzle outlet. The Newton–Euler formulation gives the dynamics of the head with respect to the body-fixed frame as follows:
is the mass of the head assembly, including its contained water.
is its matrix of inertia. The cross-product terms appear thanks to time derivatives with respect to the rotating frame. The total force acting on the head is due to the momentum change of the water flowing through it
; the gravitation
; the water pressure on the inlet cross-section
; and the reaction due to the hose isolation
. Meanwhile, the torque acting on the head is only due to the thrust at the nozzles’ outlets
, because the water inlet is at the mass center of the head. The force
and the torque
were computed by:
Moreover, let us define a mapping variable
,
, between the opening of the three-way valve
j and the resulting flow rates after it.
means that valve
j is at a neutral position and the water is distributed equally to its two outlets.
indicates that the valve is fully opened to the left side in the diagram in
Figure 1b. For example,
means all the water entering valve
a will pass to valve
b, and if
, the water entering valve
b will flow to nozzle 1 completely. In contrast,
represents that the valve is fully opened to the right side in the diagram. Hence, the complete mapping can be written as follows:
In order to obtain the representation of the head dynamics in the reference frame, substitute Equations (14) and (15), the boundary conditions (10), and the kinematics (1) into Equation (13). By neglecting higher-order terms and applying the approximation of small pitch and roll angles, the mathematical model of the head assembly is derived as in Equations (16) and (17).
In these equations, l, w, and h are the distances in the longitudinal, lateral, and vertical directions, from the coordinate origin to the outlets of the thrust nozzles, respectively. Similarly, ll, wl, and hl are those regarding the fire extinguishing nozzle. is the magnitude of the acute angle between the thrust nozzles and the vertical axis. is the acute angle between the fire extinguishing nozzle and the longitudinal axis Xb. Finally, and are the remaining disturbances and higher-order terms acting on the head’s translational and rotational movements, respectively.
Remark 4: The system is underactuated because four control inputs control the 6-DoF motion of the head, as well as maintain the stability of the flexible hose. Moreover, there is a coupling of input in the lateral motion and roll rotation, with respect to the body-fixed coordinate, of the head.
4. Control System Design
In a firefighting mission, the head assembly follows a given trajectory to approach its working position at a proper distance to a flame and then sprays the water to suppress the flame. During the process, the motion of the head should be kept stabilized, and the fluctuation of the water-conveying hose should be minimized. Furthermore, because the system is underactuated, it is not possible to control every DoF independently. In fact, the head’s longitudinal and lateral translations are coupled with its pitch and roll rotations, respectively. Hence, a composite control law is designed in this section and takes the above characteristics into account. Because it is difficult to obtain the partial states of the hose in practical operation, the proposed control law uses the feedback from the head assembly and the water supply pump. Additionally, a command filter is introduced to ensure the proposed control law is applicable.
4.1. Control Law Design
The control law to be designed consists of the pump’s required flow rate
, the opening of the valves
,
, and two virtual functions
and
of the head’s desired roll and pitch rotations. For the system head, the references of the trajectory and attitude are given by
and
. Firstly, filtered references were introduced to guarantee they are twice differentiable.
with
. In Equation (18), the hats indicate the filtered values, and
is the vector of difference between them and the unfiltered values.
is the diagonal matrix of damping ratios, and
is that of the natural frequencies of the filters. Then, a vector of tracking errors of the head was denoted with respect to the newly filtered references as follows:
To compensate for this error, a virtual control law for the head’s velocity was designed as in (20):
with
,
and
are the positive definite and diagonal matrices. In addition, a composite error was considered as follows:
with
In and
On being the identity matrix and zero matrix, respectively, of size
n. In addition to the difference between the actual and desired velocities of the head assembly (i.e.,
and
), this error representation also takes into account coupling and the saturation of the control inputs. In detail,
in Equation (21) is aimed to decouple the control input
that appears in the lateral motion.
is the first-order filtered error due to saturated control inputs. They are given by:
and in the dynamics of are positive definite diagonal matrices. is chosen such that . The bar above a variable indicates the computed values from the controller, i.e., without saturations, to distinguish from the true value (without the bars) performed on the system.
Then, the dynamics of the error
can be obtained by substituting the head dynamics in Equation (16) and the terms in Equation (22) into its time derivative:
with
Then, a saturation-free control law was designed as follows:
with
as positive definite and diagonal matrices and
as a positive constant value. From the above control vector, we derived the control inputs defined in (17), as well as the required roll and pitch rotations of the head, as follows:
As a result, the required flow rate from the pump and the opening of the flow-regulating valves were obtained by mapping (27). The overall structure of the proposed control system is depicted in
Figure 4.
4.2. Stability Analysis
Theorem. For the water-powered aerial system whose dynamics are represented by Equations (9), (10), (16), and (17), the control law (25)–(27), with a proper choice of gains and , , guarantees the uniform ultimate boundedness of the entire control system.
Proof of the Theorem. To begin with, we considered the following Lyapunov function candidate:
For proof of the Theorem, the following Lemma 1 and Lemma 2 were considered.
Lemma 1: The Lyapunov function candidate in (28) has the upper and lower boundaries established as follows.
In Equation (30), the notations and stand for the eigenvalue and Euclidean norm of , respectively.
Proof of Lemma 1: From the Cauchy–Schwarz inequality and Holder’s inequality for integral, one can derive:
Substituting inequalities (31) into Equation (28) gives the result as follows:
Moreover:
Hence, Lemma 1 is obtained. □
Lemma 2: With the proposed control law, the derivative of in (28) along the trajectory of the system satisfies:with as a positive definite matrix, as a positive definite constant, and: Proof of Lemma 2: By taking the integration of both sides of the governor Equation (9) and adopting the initial conditions of the hose, the third-order differential term
can be obtained as follows:
The total weight of the hose that the system must lift increases with its height and becomes a noticeable influence. From the properties of a line integral in a vector field and by adopting the mean value theorem for definite integrals [
43], this weight can be represented in the following way:
with
. Let us write
, where
is the nonzero constant of the mean value of
, and
is the remainder to make the equality (39) true.
Taking the time derivative of the Lyapunov candidate function, substituting the system’s dynamics and above control law, and taking integrals by parts result in (38).
with
and
Here, one can see that except , all the elements in Equation (38) are negative definite. In Equations (39) and (40), is the remaining interactive force between the hose and its surroundings after being compensated by the adaptive element in the control law, and is the difference between the head’s speed and the average speed of the hose, due to its flexibility. Then, by adopting Young’s inequality for in (39), the differential was obtained as follows:
By choosing the control matrices such that is positive definite, we thereby completed the proof of Lemma 2. □
Lemma 1 and 2 imply the uniform ultimate boundedness of the control system. From the right inequality of (31), we have:
From the left inequality of (31):
Then, the ultimate bound
b with respect to
can be taken as (44) regardless of either the hose’s curvature or the head’s fly height.
This completes the proof for the Theorem. □