4.1. State-of-the-Art
One of the first experimental demonstrations of a physical reservoir computer was a system that exploited the ripples on the surface of water contained in a tank placed on an electric motor-driven platform [
84] (for a review see, e.g., Refs. [
23,
38,
47]). The ripples caused by judiciously controlled vibrations were recorded using a digital camera. Then, the images were processed on a computer and then the processed data were used as the input of the reservoir computing algorithm.
Although that pioneering work clearly demonstrated the potential of liquids to perform RC calculations, for about two decades it was considered to be mostly of fundamental interest to researchers working on physics-inspired AI systems [
47,
48,
85]. A similar idea was exploited in the recent theoretical [
25,
82] and experimental works [
83], where it was suggested that the ripples created by motors can be replaced by the dynamical properties of solitary waves—nonlinear, self-reinforcing and localised wave packets that can move for long distances without changing their shape [
86].
Importantly, those novel works demonstrated the ability of water-based RC system to do complex calculations relying only on low computational power microcontrollers (e.g., Arduino models [
87]) that can, in principle, continuously operate for several months without the need to recharge the battery. Therefore, such systems are especially attractive for application in drones, robots and other autonomous platforms.
4.2. Towards Reservoir Computing with Water Waves Created by an ROV
Figure 3a shows the photographs of an experimental ROV designed to test different kinds of autonomous AI systems. The hull of ROV consists of two wing-like fins that hold two electric motors that rotate the propellers: one in the counterclockwise direction and another one in the clockwise direction [
90]. The fins are connected to a cylindrical frame that is also used as a compartment where sensing equipment, including lasers and photodetectors, can be placed. As with many commercial ROV designs [
91], the movement of the experimental ROV is controlled via an electric cable.
When both propellers rotate at the same speed, the ROV moves forward (
Figure 3b). However, the ROV will turn left (right) when the speed of the left (right) motor is reduced. Thus, a moving ROV creates a wake wave whose pattern depends on the relative speed of the two motors. The temporal dynamics of the wave pattern can be sensed using two laser-photodiode pairs that measures the light reflected from the waves created by the left and right propellers (a similar detection mechanism was used in Ref. [
87]). The resulting optical signals can then be converted in electric signals that are transmitted to the control unit via the cable. The waveform of the so-generated electric signals will be a function of time and it will correlate with the commands (left, right, forward and so on) given to the ROV by the operator. Indeed, as schematically shown in
Figure 3b, when the ROV moves forward the waves created by both propellers are approximately the same. However, the wave patterns change when the ROV turns left or right.
Similarly to the seminar work Ref. [
84], the differences in the wave pattern of the moving ROV can be employed to perform reservoir computations. Such computations can help predict an optimal trajectory of the ROV based on the previous inputs of the human operator and taking into account the environmental factors, including the speed of wind, temperature and salinity of water as well as the presence of obstacles such as rocks and debris [
92,
93]. Since the ROV creates the waves that are used to predict the trajectory, the on-board reservoir computer does not consume significant energy apart from the need to power the sensors, thereby satisfying the requirements for autonomous AI systems [
34].
4.3. Physical Reservoir Computing Using Fluid Flow Disturbances
The research work involving the ROV shown in
Figure 3 is still in early development. However, at the time of writing, to the best of our knowledge it is the only experimental attempt to build a prototype of a water drone that employs physical processes in its environment as a means of computation (a UGV designed using similar principles is discussed in
Section 5). Nevertheless, similar ideas have been explored theoretically and also tested in laboratory settings [
36,
83,
87,
88,
89,
94]. Although the results presented in the cited papers do not involve any autonomous vehicle, in the following we will analyse them in the context of the ROV model discussed above.
Before discussing those results, we also highlight the proposals of physical RC systems designed to predict the trajectory of a flying drone [
93,
95]. Although the idea of using physical processes taking place in the surrounding environment of the drone were not expressed in Ref. [
93,
95], the approach presented there can be extended to implement the concepts illustrated in
Figure 1a, b. In this context, we remind that the discipline of fluid dynamics is concerned with the flow of both liquids and gases [
96]. This means that AI systems designed for flying and underwater drones can, in principle, exploit the same physical processes. In fact, as mentioned above, both moving airborne and underwater vehicles create vortices as they move in the atmosphere [
97,
98,
99] and water [
21].
Vortices are ubiquitous in nature (e.g., they can be observed in whirlpools, smoke rings and winds surrounding tropical cyclones and tornados) and they exhibit interesting nonlinear dynamical that has been the subject of fundamental and applied research [
100,
101,
102,
103]. For instance, it is well-established that when fluid flows around a cylindrical object, the physical effects known as von Kármán vortex street can be observed in wide a range of flow velocities [
100,
102]. The shedding of such vortices imparts a periodic force on the object. In many situations, this force is not significant enough to accelerate the object. However, in some practical cases the object can vibrate about a fixed position, undergoing harmonic motion. Yet, when the frequency of the periodic driving force matches the natural frequency of the oscillation of the object, the resonance processes come into play and the amplitude of the oscillations can increase dramatically [
101,
103].
Therefore, it has been theoretically demonstrated that the physical properties of vortices can be used in reservoir computing [
88,
94]. In Ref. [
88], the authors conducted a rigorous numerical analysis of a vortex-based RC system based on a von Kármán vortex street (
Figure 4a). A periodic pattern of numerical sensors located across the computational domain was used to monitor the nonlinear dynamics and collect data for further processing following the traditional RC algorithm. The flow of fluid was used as the input. For example, to create an input that corresponds to a signal that varies in time, the velocity of the flow was modulated such that it follows the time-varying shape of the signal (which can readily be done using an electric pump [
83,
87]).
The dynamics of the resulting computational reservoir depends on the value of the Reynolds number
. Subsequently, different operating regimes were tested using inputs corresponding to the values of
below and above the threshold of the formation of a von Kármán vortex street [
88]. Those tests confirmed a high memory capacity of the reservoir and its ability to learn from input data and generalise them. Further tests also revealed the ability of the reservoir to make accurate predictions of time series datasets.
Figure 4.
(
a) Illustration of the vortex shedding taking place when a fluid such as air or water flows past a cylinder. As theoretically shown in Ref. [
88], modulating the flow velocity and monitoring the vortex dynamics using a set of virtual sensors one can create an efficient physical RC system. An experimental implementation of this computational approach was discussed in Ref. [
89]. (
b) Photograph of vortices and other water flow effects created by the ROV in a lab setting.
Figure 4.
(
a) Illustration of the vortex shedding taking place when a fluid such as air or water flows past a cylinder. As theoretically shown in Ref. [
88], modulating the flow velocity and monitoring the vortex dynamics using a set of virtual sensors one can create an efficient physical RC system. An experimental implementation of this computational approach was discussed in Ref. [
89]. (
b) Photograph of vortices and other water flow effects created by the ROV in a lab setting.
We note that in an experimental attempt to test the theoretical vortex-based RC system [
89] the virtual sensors shown in
Figure 4 can be replaced by real ones. Alternatively, one can use a digital camera to film the vortex and then process different pixels of the individual frames extracted from the video file [
82,
83,
87,
89]. It is also noteworthy that propellers of the ROV can also create vortices [
21,
104] that can be used for reservoir computing purposes (
Figure 4b).
4.4. Acoustic-Based Reservoir Computing
Another promising approach to physical reservoir computing employs acoustic waves, vibrations and adjacent physical processes [
105,
106,
107]. Although high-frequency (MHz-range) acoustic waves have been used in the cited papers, the ideas presented in those works can be implemented using the acoustic phenomena observed in a wide range of frequencies. For example, the temporal dynamics of vortices and other disturbances created by ROVs (
Figure 4b) has a spectral signature in the frequency range that spans from several tens of Hz to several hundreds of kHz. Such acoustic signals can be detected using hydrophones, sonar technologies and other well-established acoustic location techniques [
108,
109,
110,
111] and then processed using an RC algorithm [
112].
As with the sound radiated by a moving ship [
113], ROVs and similar autonomous vehicles can produce a tonal (related to the blade pass frequency) acoustic disturbance and broadband noise associated with the presence of unsteadiness in the flow. The tonal disturbances can be further categorised into the contributions related to such technical parameters of the propellers as the blade thickness parameter and blade loading [
113]. These acoustic processes also exhibit significant nonlinear effects [
114] that can be exploited in a computational reservoir [
38,
47,
48].
Moreover, as the propeller rotates it pushes the ROV through the water, causing a positive pressure on the face of the blade and a negative pressure on its back. The negative pressure causes any gas in solution in the water to evolve into bubbles [
115,
116,
117,
118,
119]. These bubbles collapse via the process called cavitation, causing hammer-like impact loads on the blades and damaging their surface [
117,
118,
119]. The cavitation also causes significant acoustic noise that originates both from oscillations and collapse of bubbles [
117] and the formation of vortices [
20,
92]. This physical picture is sketched in
Figure 5.
In particular, the underwater acoustic noise is associated with the highly-nonlinear oscillation of the bubble volumes that typically occur in the range of frequencies from several hundred Hz to approximately 40 kHz [
116,
117,
120,
121] (microscopic bubbles oscillate at higher frequencies [
122,
123,
124]; however, the physics of their acoustic response remains essentially the same). Indeed, considering an idealised scenario of a sinusoidal acoustic pressure wave, when a wave moves through water, its initial waveform changes so that its initial monochromatic spectrum acquires higher harmonic frequencies [
121]. The more nonlinear the medium in which the wave propagates, the stronger the enrichment of the spectrum with the peaks corresponding to harmonics.
The degree of acoustic nonlinearity can be characterised by the acoustic parameter
, which is the ratio of coefficients
B and
A of quadratic and linear terms in the Taylor series expansion of the equation of state of the medium (see [
121] and references therein). The larger the value of
, the more nonlinear the medium, the stronger a distortion of the acoustic spectrum from the initial monochromatic state. For instance, water with
is more acoustically nonlinear than air with
(we note that the degree of nonlinearity is considered to be moderate in both media).
However, when air bubbles are present in water, the value of
increases to around 5000 [
121]. This is because liquids are dense and have little free space between molecules, which explains their low compressibility. In contrast, gases are easily compressible. When an acoustic pressure wave propagating in water reaches a bubble, due to the high compressibility of the gas trapped in it, its volume changes dramatically, causing substantial local acoustic wavefront deformations that result in strong variation of the initial acoustic spectrum.
The evolution of the acoustic spectrum of an oscillating bubble trapped in the bulk of water is illustrated in
Figure 6 (for the computational details and model parameters see, e.g., [
118,
119,
123,
125]). The units of the
x-axis of this figure correspond to the normalised frequency
, where
is the frequency of the incident sinusoidal acoustic pressure waves. The
y-axis corresponds to the peak pressure of the incident wave (in kPa units) but the false colour encodes the amplitude (in dB units) of the acoustic pressure scattered by the bubbles. The bright traces with the amplitude of approximately 0 dB correspond to the frequency peaks in the spectra of the bubble forced at one particular value of the peak pressure of the incident wave. We can see that at a relatively low pressure of the incident wave the spectrum contains the frequency peaks at the normalised frequencies
and so forth. However, an increase in the pressure results in the generation of the subharmonic frequency
and its ultraharmonic components
,
and so on. Further increase in the peak pressure of the incident waves leads to a cascaded generation of subharmonics frequency peaks, resulting in a comb-like spectrum [
119,
125].
The nonlinear response of a cluster of oscillating bubbles trapped in water was used to create a computational reservoir in theoretical work Ref. [
36]. Each bubble in the cluster oscillates at a certain frequency when the entire cluster is irradiated with an acoustic wave. The oscillation frequency of each bubble in the cluster depends on the equilibrium radius of the bubble and the strength of its interaction with the other bubbles in the cluster. Importantly, the cluster maintains its structural stability (i.e. the bubble do not merge) when the pressure of the incident acoustic wave remains relatively low [
125].
Thus, when the temporal profile of the incident pressure wave is modulated to encode the input data (e.g., a time series that needs to be learnt and then forecast by the RC system), the cluster of bubbles acts as a network of artificial neurons. Since each oscillating bubble emits sound waves, these waves can be detected using either a hydrophone or a laser beam [
118,
119] and then process following to the traditional RC algorithm. As demonstrated in Ref. [
36], this computational procedure enables the RC system to predict the future evolution of highly nonlinear and complex time series with the efficiency of the traditional RC algorithm while using low energy consuming computational resources.
Subsequently, it is plausible that bubbles created by a moving ROV (
Figure 5) can be used to construct an onboard RC system. The input of such a reservoir will be the local pressure variations caused by the propellers of the ROV. It can be shown that these variations correlate with the control signals received by the ROV from the operator and/or its onboard control unit [
127,
128]. While the implementation of such a computational scheme requires resolving several technological problems, it enables the researchers to access a rich spectrum of fascinating nonlinear effects associated with oscillating bubbles [
117,
121,
123,
124,
129,
130]. Interestingly enough, as demonstrated in the following section, a similar, from the point of view of nonlinear dynamics, idea was proposed in the domain of autonomous ground vehicles.