ORIGINAL RESEARCH article
Hydrodynamical Fingerprint of a Neighbour in a Fish Lateral Line
- 1Center for Mathematical Science and Advanced Technology, Japan Agency for Marine-Earth Science and Technology (JAMSTEC), Yokohama, Japan
- 2Center for Design, Manufacturing and Materials (CDMM), Skolkovo Institute of Science and Technology, Moscow, Russia
- 3Graduated School of Engineering, Chiba University, Chiba, Japan
- 4Laboratoire de Physique et Mécanique des Milieux Hétérogènes (PMMH), CNRS UMR 7636, ESPCI Paris—PSL University, Sorbonne Université, Université de Paris, Paris, France
For fish, swimming in group may be favorable to individuals. Several works reported that in a fish school, individuals sense and adjust their relative position to prevent collisions and maintain the group formation. Also, from a hydrodynamic perspective, relative-position and kinematic synchronisation between adjacent fish may considerably influence their swimming performance. Fish may sense the relative-position and tail-beat phase difference with their neighbors using both vision and the lateral-line system, however, when swimming in dark or turbid environments, visual information may become unavailable. To understand how lateral-line sensing can enable fish to judge the relative-position and phase-difference with their neighbors, in this study, based on a verified three-dimensional computational fluid dynamics approach, we simulated two fish swimming adjacently with various configurations. The lateral-line signal was obtained by sampling the surface hydrodynamic stress. The sensed signal was processed by Fast Fourier Transform (FFT), which is robust to turbulence and environmental flow. By examining the lateral-line pressure and shear-stress signals in the frequency domain, various states of the neighboring fish were parametrically identified. Our results reveal that the FFT-processed lateral-line signals in one fish may potentially reflect the relative-position, phase-differences, and the tail-beat frequency of its neighbor. Our results shed light on the fluid dynamical aspects of the lateral-line sensing mechanism used by fish. Furthermore, the presented approach based on FFT is especially suitable for applications in bioinspired swimming robotics. We provide suggestions for the design of artificial systems consisting of multiple stress sensors for robotic fish to improve their performance in collective operation.
Collective behavior has been proposed to benefit animals in many aspects, such as predator avoidance, feeding, reproduction, migration and social learning (Partridge, 1982; Pavlov and Kasumyan, 2000; Vicsek and Zafeiris, 2012). Particularly in fish, collective behavior has been recognized as a means to improve energetic performance in swimming (Weihs, 1973; Abrahams and Colgan, 1985; Chen et al., 2016; Ashraf et al., 2017; Filella et al., 2018).
In a fish school, relative position and kinematic synchronisation between neighboring fish are major factors that influence the energetic expenditure of individual fish. Weihs proposed that a diamond-shaped formation would enable one fish to take advantage of the low relative speed passage between the vortex streets formed by two preceding fish (Weihs, 1973). Observations on sea bass by Herskin and Steffensen suggest that swimming at the front of a school was significantly costlier than swimming at the rear (Herskin and Steffensen, 1998). Marras et al. reported that in a fish school, individuals in any position had reduced energetic expenditure: fish swimming behind their neighbors save the most energy, and those fish swimming ahead of their nearest neighbor gain a minor net energetic benefit over swimming solely (Marras et al., 2015). Ashraf et al., by contrast, suggested that a simple “phalanx” formation (swimming in a line) may provide a significant energetical benefit (Ashraf et al., 2017). The choice of relative phase is associated with the relative position between fish: preference of synchronization between nearby fish (in-phase or anti-phase) is observed in fish swimming experiments (Ashraf et al., 2016, 2017), while the beneficial outcome of synchronization are confirmed numerically (Li et al., 2019a) and experimentally (Godoy-Diana et al., 2019). Besides synchronization, recent studies also suggest that matching the vortex phase of a neighboring fish is an effective means to improve the energetic efficiency (Daghooghi and Borazjani, 2015; Khalid et al., 2016; Li et al., 2020). In that case, the optimal phase between neighbors is dynamic and should be adjusted according to their relative position and wake morphology.
Effective sensing is important to maintain the group configuration (Puckett et al., 2018). The sensory basis of fish schooling includes multiple systems such as the lateral line and vision (Partridge and Pitcher, 1980; Partridge, 1981). Attenuated vision may alter schooling behavior: fish are observed to swim slower in low light environment (Berdahl et al., 2013), and keep larger relative distance (Hunter, 1968). Nevertheless, fish can school even under a blindfolding condition, while blinding had little effect on the position that experimental fish took up with respect to their neighbors within the school (Pitcher et al., 1976; Partridge and Pitcher, 1980). This indicates the significant role of non-visual sensing in fish schooling behavior. On the other hand, ablating the lateral-line also alters schooling behavior (Partridge and Pitcher, 1980; Faucher et al., 2010; Mekdara et al., 2018). Fish may still school when the posterior lateral-line system is disabled, but they reduce the distance to neighboring fish and may occasionally collide (Partridge and Pitcher, 1980). A more recent study reports that, without functioning superficial neuromasts, schooling behavior was disrupted under both photopic and scotopic conditions and the ability to detect stationary objects decreased (Middlemiss et al., 2017). Experiment on rummy-nose tetra fish shows that, when the whole lateral-line system is inactivated, fish cannot maintain the schooling behavior and swim with greater distances between neighbors (Faucher et al., 2010). Mekdara et al. found that tail beat synchronization during schooling requires a functional posterior lateral-line system in Giant Danios, and they hypothesize that the anterior branch may be more important for regulating position within the school, while the posterior branch may be more important for synchronizing tail movements (Mekdara et al., 2021).
Lateral-line sensing is realized by a mechanosensory system, comprised of arrays of sensors called neuromasts, which respond to the motion of the surrounding water relative to the skin of the fish, and to pressure and shear stress changes (Coombs et al., 1988; Montgomery et al., 1995; Coombs, 2001; Ren and Mohseni, 2012). There are two main types of lateral-line organs in lower vertebrates: superficial neuromasts (SN), with a cupula that protrudes in the surrounding water, and canal neuromasts (CN), located in the lateral-line canal. The scales of the trunk lateral-line canal of fish contain SNs as well as CNs (Kroese and Schellart, 1992; Coombs, 2001). While the neural response of the hair cell is proportional to the displacement of the cupula and the underlying ciliary bundles (van Netten and Kroese, 1987), cupula displacement is largely proportional to the velocity of water flowing past it, coupling the motion of the surrounding water to the underlying cilia through viscous forces. On the other hand, inertial forces are required for the fluid to break through the boundary layer and move into the small-diameter canal. Hence, flow velocity inside the canal is more or less proportional to the net acceleration between the fish and the surrounding water, CNs respond to changes in external flow acceleration and to net pressure differences between the two surrounding canal pores (Denton and Gray, 1983; Kroese and Schellart, 1992). CNs and SNs may have different response characteristics: SNs best respond to the direct current and low-frequency components of the incoming flow, whereas CNs respond best to high-frequency components of the flow (Coombs et al., 2001).
Artificial lateral lines with biomimetic neuromasts have been developed and applied to bioinspired swimming robots (Fan et al., 2002; Yang et al., 2010; Xu and Mohseni, 2017; Liu et al., 2020). Optimal sensor locations for artificial swimmers have been also investigated to improve the efficiency (Verma et al., 2019; Weber et al., 2020), especially when the number of biomimetic neuromasts is limited. The sensing of external hydrodynamic stress may improve the control robustness in single-fish self-organized undulatory swimming (Thandiackal et al., 2021), as well as the propulsive efficiency (Tytell et al., 2010). In collective swimming, the interaction of the flow field and the surface stress between neighboring fish have been established by several recent studies (Daghooghi and Borazjani, 2015; Chen et al., 2016; Peng et al., 2018, Li et al., 2019a; 2019b; Verma et al., 2019). This interaction is the basis for perceiving the collective swimming configuration by stress-sensing. It is reasonable to expect that fish tune the collective swimming spatio-temporal patterns and improve the collective swimming efficiency based on stress information (Novati et al., 2017; Verma et al., 2018).
The aim of this study is to investigate, using a numerical model, the stress signal features induced by a nearby swimming fish, and examine how to identify the state of a neighbor fish based on lateral-line sensing. The general objective is to assess to what extent the sensing of flow fluctuations through lateral line sensors can be used to determine decision making in schooling fish. Since fish deformation and deformation-induced stress are periodic, we apply a Fourier Transform to identify the signal changes with respect to different frequency components. The zero-frequency component (constant in Fourier Transform), the tail-beat frequency component, and the higher-frequency ones correspond, respectively, to large-scale environmental flow, to the fish own tail-beat behavior, and to turbulence. Previous biological studies provide basis for such decomposition: in the fish peripheral nervous system, superficial neuromasts respond best to slow unidirectional flows, whereas canal neuromasts respond best to more rapidly fluctuating flows (e.g., 30-100 Hz); in the fish central nervous system, there is also a filter mechanism to increase signal-to-noise ratios by filtering-out the fish’s own movements (see Coombs, 2001, for a review). This analysis based on Fourier Transform provides a new angle to quantitatively understand the hydrodynamic interaction in collective swimming, and is especially suitable to for implementation in swimming robotic control.
Fish Morphology and Kinematics
We simulated a virtual swimmer with body shape based on the profile of an adult trout (Supplemenatary Figure S12). The body length is set as 4 cm to limit the simulation within laminar flow, and the body wave was prescribed by sinusoidal functions:
Three Dimensional Navier-Stokes Solver
The hydrodynamic solution is based on a validated three-dimensional Navier-Stokes (NS) solver (Liu, 2009; Li et al., 2012; 2014; Bomphrey et al., 2017). The governing equations for fluid solution are three-dimensional, incompressible and unsteady NS equations written in strong conservation form for mass and momentum (Liu, 2009). To accelerate the computation and improve the robustness during iteration, the artificial compressibility method is adopted by adding a pseudo time derivative of pressure to the continuity equation. The solving process is implemented in multiple grids using the finite volume method (FVM). Navier-Stokes equations were solved in each grid and results were interpolated at the grid interfaces.
The approach comprises surface models of the changing fish shape (dimension: 133 × 97), and local fine-scale body-fitted grids (dimension: 133 × 97 × 20) plus a large stationary global grid (dimension: various) to calculate the flow patterns around the fish with sufficient resolution. The boundary conditions are set as follows: 1) in the fish-body-fitted grid, the non-slip condition is applied to the cells on the surface of the fish body; 2) in the global grid, an incoming flow
We implemented pre-simulations on a single fish to determine an equilibrium speed U of 2.375L s−1 (9.5 cm s−1), at which the drag balances the propulsion by the fish. In all the simulations in this study, the centers of mass of fish were fixed, while the incoming flow was set as U, so that the hydrodynamics of the tethered fish were similar to those in realistic steady swimming. The Reynolds number of the simulations is defined as
The computational accuracy in terms of stress and flow around a fluctuating rigid cylinder is validated in §B-8 in the Supplemenatary Materials of (Li et al., 2021)). The validation of flow simulation accuracy around a deforming fish body is provided in Figure 9 and Fig. 18 in (Li et al., 2012), which can also be an indirect validation of stress accuracy on a deforming fish body, since surface stress is coupled with near-surface flow motion. Further details of the NS solver are provided in Part C of Supplemenatary Materials.
Virtual Lateral Line
We assume that the model fish possesses a continuous virtual lateral line at mid-body height, from snout to tail tip, on both sides of the body. The virtual lateral-line may accurately sense both local pressure- and shear-stresses. The lateral-line stress signal X (s, b, l) during an arbitrary tail-beat cycle is sampled 40 times (sampling rate = 320 Hz), where s represents the stress type (s = pressure or shear stress), b indicates the body side (b = left or right); and l represents the sensor position along the body (0 < l < L).
Stress Signal Processing by Fast Fourier Transform (FFT)
Due the specific nature of the fish-like swimming kinematics, the pressure- and shear stress-over the fish surface shows strong periodicity: a frequency-domain analysis is thus ideal to examine the lateral-line stress signal. The lateral-line stress signal X during an arbitrary tail-beat cycle is sampled and processed by a fast Fourier transform (FFT) algorithm (function: fft, MATLAB R2020b, The Mathworks).
Signals are processed and decomposed with respect to the dimensionless frequency k used in this study, computed as k = f/fref, where the k = 1 component corresponds to a sinusoidal signal in tail-beat frequency, and k = 0 corresponds to the direct current component. We will use in the paper the direct current (DC), alternating current (AC) nomenclature. The single-sided spectrum is computed for all components. Note that for an AC component, the single-sided spectrum only shows its magnitude, while for the DC component, we preserve the sign of signal since the negativity of the DC component represents hydrodynamic information of the stress direction. The phase of each component is computed by the MATLAB function angle.
Further details of the Fast Fourier Transform are provided in §C-5 of Electronical Supplemenatary Materials.
Basal Signal in a Swimming Fish Alone
The lateral-line signal of a swimming fish alone provides a basal value to estimate the relative stress change associated to collective swimming. Since the deformation is governed by a sinusoidal function, the curvature time-series periodically repeats a same pattern. We choose an arbitrary tail-beat cycle (from t0 to t0+T) and the curvature map is demonstrated in Figure 1A.
FIGURE 1. Deformation and lateral-line stress in solely swimming fish. (A) curvature map of fish axis in an arbitrary tail-beat cycle (from t0 to t0+T); (B) Flow field surrounding a solely swimming fish; (C) pressure-stress map of the left-side lateral-line; (D) Frequency domain analysis of pressure stress signal at five locations along the lateral-line; (E) Distribution of single-sided spectrum of left-side pressure, at k = 0, 1, and 2; (F) shear-stress map of the left-side lateral-line; (G) Frequency domain analysis of shear-stress signal at five locations along the lateral-line; (H) Distribution of single-sided spectrum of left-side shear-stress signal, at k = 0, 1, and 2. One stress unit = 10 Pa.
The left-side lateral-line pressure signal map during one tail-beat cycle is displayed in Figure 1C. The pressure-signal map sensed by the right-side lateral-line is omitted since it is of an identical pattern to the left-side but with a half-tailbeat-cycle phase difference (
Figure 1F shows the shear-stress signal map sensed by the left lateral-line from t0 to t0+T. Unlike the pressure stress that shifts between positive and negative values, the shear-stress signal map is overall negative, indicating that shear stress consistently hinders the forward motion of the fish. The magnitude of the DC signal is usually stronger than the
FIGURE 2. Phases of
Presence of an Alongside, In-phase Neighbor
The presence of an alongside neighbor may influence the lateral-line signal in several aspects. Figure 3 shows the signal features sensed by the protagonist fish when an alongside, in-phase neighbor is 0.35 L away on its right side. In terms of pressure stress, the DC term of both left and right sides changed: on the right side (neighbor-ward) the DC pressure signal becomes greater at the head and tail zones, and weakens in the mid-region of the body. On the left side (free side) the DC pressure signal of the anterior body slightly decreases, suggesting that the free side may still be influenced by the presence of the neighbor fish. A similar phenomenon (stress influence across the body) is also reported in Verma et al. (2019). For k = 1 AC signal, the pressure on the free side is consistent with the reference signal, while it is overall weak on the neighbor-ward side. The latter is considered as a key feature in this collective swimming configuration. In terms of shear stress, on the neighbor-ward side the shear stress DC signal is reduced in the anterior part and increased in the posterior part. The k = 1 AC shear signal shows an increment in the middle of body (0.3–0.5L), and since the basal signal is relatively low at this location, there is an approximately 30% increment in proportion.
FIGURE 3. Signal features sensed by the protagonist fish when an alongside, in-phase neighbor is 0.35 L away on its right side. Full results are presented in Supplementary Figure S1, ESM. One stress unit = 10 Pa.
Dependence on the Distance of the Alongside In-phase Neighbor
As the alongside, in-phase (
FIGURE 4. Signal features sensed by the protagonist fish when an alongside, in-phase neighbor is 0.5 L,0.75 and 1 L away on its right side. Full results are presented in Supplementary Figures S2-S4, ESM. One stress unit = 10 Pa.
The alongside, in-phase neighbor may also cause changes in higher order (k > 2) AC components. These results are shown in Supplemenatary Figure S2-S4, ESM.
Phase Difference With the Alongside Neighbor
The phase difference between two adjacent fish also induces changes in the stress signal. As shown in Figure 5, when a neighbor fish swims alongside the protagonist with various phase differences
FIGURE 5. Signal features sensed by the protagonist fish when an alongside neighbor is 0.35 L away on its right side, with
The phase of the lateral line signal is also sensitive to the phase difference with a neighbor fish. As shown in Figure 5, at Location 2, the k = 1 AC pressure signal phase on the neighbor-ward side shifts from the reference (i.e., the gap between the red and black solid lines). The phase-shift is most significant at
In contrast to the neighbor-ward side, the pressure stress signal on the free-side is very similar to the reference value (Figure 5: the difference between the green and blue solid lines is small), suggesting fish may sense a phase shift between its neighbor-ward and free sides, and may further judge the phase of the neighbor fish based on the relative change of magnitude and phase of the k = 1 AC pressure signal between its two sides.
In terms of the DC component (for full results, see Supplemenatary Figure S1, Supplemenatary Figure S5-S7, ESM), for both pressure and shear stress at
When the neighbor fish is located in the front diagonal of the protagonist fish (see Figure 6, LHS), the posterior body of the protagonist fish is exposed to the wake of the neighbor fish and the protagonist fish may easily sense the associated flow fluctuations. As shown in the LHS of Figure 6, on the neighbor-ward side, both DC and AC components of pressure- and shear-stress signals in the posterior body undergo noticeable changes. For the pressure, the DC and the k = 1 and 2 AC components show a magnitude fluctuation mainly in posterior body, while the anterior part signal is less influenced. The k = 3 AC components of pressure signal undergo complex magnitude changes along the entire body, with strengthened and weakened zones in a staggered pattern. Concerning the shear-stress signal, the DC and k = 1 AC components show a magnitude reduction in the posterior part of the body. The k = 2 and 3 AC components of the pressure signal fluctuate along the entire body, being strengthened in some locations, and weakened in others. In terms of the phase, a shift is observed in the k = 1 AC components in pressure- and shear-signals: the k = 1 AC pressure-signal of the posterior part generally shifts to a delayed phase, while the k = 1 AC shear-signal of the posterior part shows a fluctuation of phase modification. These signal changes concentrate on the neighbor-ward side of the body. On the free-side, the signal modification is insignificant.
FIGURE 6. Signal features sensed by the protagonist fish when an in-phase neighbor is in diagonal front position (LHS), and diagonal rear position (RHS). Full results are presented in Supplementary Figures S8,S9, ESM. One stress unit = 10 Pa.
On the contrary, when the protagonist fish is located in the front diagonal of the neighbor fish (see Figure 6, RHS), the protagonist fish senses very minor changes in both pressure- and shear-components, The slight change concentrates in the posterior body, still suggesting a possibility to sense a fish behind. The phase of k = 1 AC pressure- and shear-signals hardly have any notable change.
Right Behind a Frontal Neighbor
When the neighbor fish is right in front of the protagonist fish, the DC pressure signal in the head zone of the protagonist fish decreases dramatically (e.g., DC pressure at around 0.05 L almost vanished), and the k = 1 AC component also shows moderate decrement in the head zone (see Figure 7, pressure). The higher order pressure signals, k = 2 and k = 3 AC, undergo irregular fluctuations. The shear stress on the entire body is reduced. Meanwhile, irregular extra fluctuations are observed in all AC shear signals (see Figure 7, shear). When the protagonist fish is right behind a neighbor fish, the stress modification on the left- and right-side body is generally same, which is regarded as a feature to identify such symmetric configurations.
FIGURE 7. Signal features sensed by the protagonist fish when an in-phase neighbor is in front. Full results are presented in Supplementary Figure S10, ESM. One stress unit = 10 Pa.
Since the frontally swimming fish forms a rear zone with slower flow speed (see Figure 8B in Li et al. (2019a)), the protagonist fish in this zone may confront relatively slower oncoming flow, thus its frontal pressure and overall shear stress all attenuate, which reduces both form drag and friction. Meanwhile, the relatively slow oncoming flow will reduce the local slip ratio (the ratio of local flow speed to the fixed body-wave speed) at the posterior body, and the posterior body may generate stronger propulsion, expressed as the increment of the DC and k = 1 AC pressure-stress components at the posterior body.
FIGURE 8. Signal features sensed by the protagonist fish when an alongside neighbor uses double frequency as the protagonist fish. Full results are presented in Supplementary Figure S11, ESM. One stress unit = 10 Pa.
When the Neighbor Fish Uses a Different Tail-Beat Frequency
In all aforementioned simulations, both neighbor and protagonist fish use the same tail-beat frequency. In order to investigate if the changes in the tail-beat frequency of the neighbor fish can also be detected by the lateral-line stress signal, we doubled the tail-beat frequency of the neighbor fish with respect to that of the protagonist fish and simultaneously reduced the tail-beat amplitude to maintain a gross net-force balance. As shown in Figure 8, the doubled frequency used by the neighbor fish causes a dramatical change in the k = 2 AC pressure component, while the influences on k = 1 AC pressure component is much smaller than that observed in the case with the two fish using the same frequency (Figure 3). Interestingly, the DC pressure signal hardly reflects the change of the tail-beat frequency of the neighbor fish, retaining a similar pattern to that in Figure 3. As to the shear-stress signal, all components are insignificantly influenced. The shift of the AC pressure signal magnitude from k = 1 to k = 2 component suggests that the frequency of the neighbor fish can be sensed according to the frequency domain analysis. In contrast, DC signals of pressure and shear stress seem robust against a tail-beat frequency change of the neighbor fish.
This study provides a unique approach to understand the stress signals sensed by a fish lateral line. Based on Fourier Transform, the magnitude and phase in the frequency domain show characteristics that may be used to identify the relative position, phase and other status of a neighbor fish, even in the absence of a visual signal. It is worth to note that stress signals in the solo vs pair conditions to a large extent coincide with each other, and whether their difference can be distinguished by a fish or robot should still depend on a sufficiently accurate sensory system. Previous observations by Kanter and Coombs (Kanter and Coombs, 2003) suggest that the fish sensory system is rather accurate and robust to background flow where fish typically swim. They found that sculpin (Cottus bairdi) were able to detect relatively weak, prey-like signals in the presence of a strong ambient background flow, and further estimated that fish may sense flow changes several orders of magnitude below the background flow levels. The sensory system is also reported to concentrate at locations where changes in pressure are greatest during motion (Ristroph et al., 2015) to further strengthen sensing accuracy. In this study, the characteristic stress difference caused by a neighbour fish is usually one order of magnitude below the basal stress. Therefore, the sensitivity of a fish lateral line system, according to the literature, seems enough to detect the difference. Furthermore, since the water movements caused by a specific animal species usually show unique features (e.g., see Figure 1 in Mogdans (2019)), schooling fish could be instinctively sensitive to the characteristic signal formed by a neighbour fish of the same species, and less interfered by environmental noise.
Overall, the DC component of shear and pressure stresses can be used to judge the lateral distance of an alongside neighbor, while the magnitude and phase of the k = 1 AC pressure component can reflect the phase of an alongside neighboring fish. Similar to the pressure signal, the shear stress signal on the neighbor-ward side will change and indicates the lateral distance of the neighbor, but unlike the pressure signal, the shear signal is robust to the phase difference between the two fish (see Figure 5), and can provide a clear signal when the fish enters the wake of another fish (see Figures 6, 7). Since multiple vortices are generated during one tail-stroke, and further break down into small vortices, the wake flow causes lateral-line stress fluctuations with a frequency higher than the tail-beat frequency. Exposure to the wake of the neighbor cause changes in k
Regarding the influence caused by the phase difference in a fish pair, we only discussed one specific relative position (
There are limitations due to the numerical approach used in this study. Firstly, we simulated laminar flow, thus turbulence formed on the body surface and during the breakdown of the wake are beyond our simulation. Turbulence may disturb the lateral-line signal (Coombs et al., 1999), especially the higher terms in our analysis (k
It remains to be rigorously confirmed whether the sensing and nervous system of fish can process stress signals effectively filtering frequencies in a similar fashion to the Fourier Transform used in our study. The experimental-numerical study by Van Trump and McHenry, (2008) reports a 30-fold range in the amplitude of sensitivity and more than a 200-fold range of variation in cut-off frequency in lateral line sensing in zebrafish larvae (Danio rerio), which suggests that natural variation in cupular height within a species is capable of generating large differences in their mechanical filtering and dynamic range. On the other hand, even if a biological sensory system does not really perform FFT decomposition, we consider that this approach is especially suitable to engineering systems such as swimming robotics. Our study reminds that the stress signal sensed by an artificial lateral line should not be simply processed by averaging, and that the signal oscillating components carry useful information and should not be discarded. FFT is suitable for analyzing unsteadiness: high-frequency turbulence, low frequency ambient currents and self-induced flow can all be processed using FFT, leaving important DC and low-order AC components to indicate the state of a neighbor fish. Also, FFT can be implemented in artificial electronical systems with high efficiency. The stress signal regulation detected by FFT in this study must be applied flexibly in robotic applications, since robotic swimmers with variable morphological and kinematical properties may not possess universal signal features. We recommend researchers to implement preliminary tests to construct a database to connect the FFT stress-signal features with group status for their own robot swimmers, so that by comparing to the database, the regulation detected by FFT processing can be used to judge the schooling status in robot fish and assist to maintain group formation.
Data Availability Statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.
Conceptualization: GL Project administration: GL Methodology and Validation: GL, HL Simulation: GL Visualization: GL Results interpretation: GL, DK, RG-D, and BT Writing and revision: All authors.
GL is funded by the Japan Society for the Promotion of Science (20K14978).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
We would like to thank Dr(s) Liang Li, Chen Wang and Sridhar Ravi for hosting the article collection “Robotics to Understand Animal Behaviour”.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/frobt.2022.825889/full#supplementary-material
Ashraf, I., Bradshaw, H., Ha, T.-T., Halloy, J., Godoy-Diana, R., and Thiria, B. (2017). Simple Phalanx Pattern Leads to Energy Saving in Cohesive Fish Schooling. Proc. Natl. Acad. Sci. USA. 114, 9599–9604. doi:10.1073/pnas.1706503114
Ashraf, I., Godoy-Diana, R., Halloy, J., Collignon, B., and Thiria, B. (2016). Synchronization and Collective Swimming Patterns in Fish ( Hemigrammus Bleheri ). J. R. Soc. Interf. 13, 20160734. doi:10.1098/rsif.2016.0734
Berdahl, A., Torney, C. J., Ioannou, C. C., Faria, J. J., and Couzin, I. D. (2013). Emergent Sensing of Complex Environments by Mobile Animal Groups. Science. 339, 574–576. doi:10.1126/science.1225883
Chen, S.-Y., Fei, Y.-H. J., Chen, Y.-C., Chi, K.-J., and Yang, J.-T. (2016). The Swimming Patterns and Energy-Saving Mechanism Revealed from Three Fish in a School. Ocean Eng. 122, 22–31. doi:10.1016/j.oceaneng.2016.06.018
Coombs, S., Janssen, J., and Webb, J. F. (1988). “Diversity of Lateral Line Systems: Evolutionary and Functional Considerations,” in Sensory Biology of Aquatic Animals (New York, NY: Springer), 553–593.
Fan, Z., Chen, J., Zou, J., Bullen, D., Liu, C., and Delcomyn, F. (2002). Design and Fabrication of Artificial Lateral Line Flow Sensors. J. Micromech. Microeng. 12, 655–661. doi:10.1088/0960-1317/12/5/322
Faucher, K., Parmentier, E., Becco, C., Vandewalle, N., and Vandewalle, P. (2010). Fish Lateral System Is Required for Accurate Control of Shoaling Behaviour. Anim. Behav. 79, 679–687. doi:10.1016/j.anbehav.2009.12.020
Herskin, J., and Steffensen, J. F. (1998). Energy Savings in Sea Bass Swimming in a School: Measurements of Tail Beat Frequency and Oxygen Consumption at Different Swimming Speeds. J. Fish Biol. 53, 366–376. doi:10.1111/j.1095-8649.1998.tb00986.x
Khalid, M. S. U., Akhtar, I., and Dong, H. (2016). Hydrodynamics of a Tandem Fish School with Asynchronous Undulation of Individuals. J. Fluids Structures. 66, 19–35. doi:10.1016/j.jfluidstructs.2016.07.008
Li, G., Liu, H., Müller, U. K., Voesenek, C. J., and van Leeuwen, J. L. (2021). Fishes Regulate Tail-Beat Kinematics to Minimize Speed-specific Cost of Transport. Proc. R. Soc. B. 288. doi:10.1098/rspb.2021.1601
Li, G., Müller, U. K., Van Leeuwen, J. L., and Liu, H. (2014). Escape Trajectories Are Deflected when Fish Larvae Intercept Their Own C-Start Wake. J. R. Soc. Interf. 11, 20140848. doi:10.1098/rsif.2014.0848
Li, L., Nagy, M., Graving, J. M., Bak-Coleman, J., Xie, G., and Couzin, I. D. (2020). Vortex Phase Matching as a Strategy for Schooling in Robots and in Fish. Nat. Commun. 11, 1–10. doi:10.1038/s41467-020-19086-0
Liu, G., Hao, H., Yang, T., Liu, S., Wang, M., Incecik, A., et al. (2020). Flow Field Perception of a Moving Carrier Based on an Artificial Lateral Line System. Sensors. 20, 1512–1522. doi:10.3390/s20051512
Marras, S., Killen, S. S., Lindström, J., McKenzie, D. J., Steffensen, J. F., and Domenici, P. (2015). Fish Swimming in Schools Save Energy Regardless of Their Spatial Position. Behav. Ecol. Sociobiol. 69, 219–226. doi:10.1007/s00265-014-1834-4
Mekdara, P. J., Nasimi, F., Schwalbe, M. A. B., and Tytell, E. D. (2021). Tail Beat Synchronization During Schooling Requires a Functional Posterior Lateral Line System in Giant Danios, Devario Aequipinnatus. Integr. Comp. Biol. 61, 427–441. doi:10.1093/icb/icab071
Mekdara, P. J., Schwalbe, M. A. B., Coughlin, L. L., and Tytell, E. D. (2018). The Effects of Lateral Line Ablation and Regeneration in Schooling Giant Danios. J. Exp. Biol. 221. doi:10.1242/jeb.175166
Middlemiss, K. L., Cook, D. G., Jerrett, A. R., and Davison, W. (2017). Morphology and Hydro-Sensory Role of Superficial Neuromasts in Schooling Behaviour of Yellow-Eyed Mullet (Aldrichetta Forsteri). J. Comp. Physiol. A. 203, 807–817. doi:10.1007/s00359-017-1192-6
Mogdans, J. (2019). Sensory Ecology of the Fish Lateral‐line System: Morphological and Physiological Adaptations for the Perception of Hydrodynamic Stimuli. J. Fish. Biol. 95, 53–72. doi:10.1111/jfb.13966
Novati, G., Verma, S., Alexeev, D., Rossinelli, D., Van Rees, W. M., and Koumoutsakos, P. (2017). Synchronisation Through Learning for Two Self-Propelled Swimmers. Bioinspir. Biomim. 12, 036001. doi:10.1088/1748-3190/aa6311
Ristroph, L., Liao, J. C., and Zhang, J. (2015). Lateral Line Layout Correlates with the Differential Hydrodynamic Pressure on Swimming Fish. Phys. Rev. Lett. 114, 1–5. doi:10.1103/PhysRevLett.114.018102
Thandiackal, R., Melo, K., Paez, L., Herault, J., Kano, T., Akiyama, K., et al. (2021). Emergence of Robust Self-Organized Undulatory Swimming Based on Local Hydrodynamic Force Sensing. Sci. Robot. 6, eabf6354. doi:10.1126/scirobotics.abf6354
Tytell, E. D., Hsu, C.-Y., Williams, T. L., Cohen, A. H., and Fauci, L. J. (2010). Interactions Between Internal Forces, Body Stiffness, and Fluid Environment in a Neuromechanical Model of Lamprey Swimming. Proc. Natl. Acad. Sci. 107, 19832–19837. doi:10.1073/pnas.1011564107
van Netten, S. M., and Kroese, A. B. A. (1987). Laser Interferometric Mea- Surement on the Dynamic Behavior of the Cupula in the Fish Lateral Line. Hear. Res. 29, 55–61. doi:10.1016/0378-5955(87)90205-x
Van Trump, W. J., and McHenry, M. J. (2008). The Morphology and Mechanical Sensitivity of Lateral Line Receptors in Zebrafish Larvae (Danio rerio). J. Exp. Biol. 211, 2105–2115. doi:10.1242/jeb.016204
Verma, S., Novati, G., and Koumoutsakos, P. (2018). Efficient Collective Swimming by Harnessing Vortices Through Deep Reinforcement Learning. Proc. Natl. Acad. Sci. USA. 115, 5849–5854. doi:10.1073/pnas.1800923115
Xu, Y., and Mohseni, K. (2017). A Pressure Sensory System Inspired by the Fish Lateral Line: Hydrodynamic Force Estimation and Wall Detection. IEEE J. Oceanic Eng. 42, 532–543. doi:10.1109/JOE.2016.2613440
Yang, Y., Nguyen, N., Chen, N., Lockwood, M., Tucker, C., Hu, H., et al. (2010). Artificial Lateral Line with Biomimetic Neuromasts to Emulate Fish Sensing. Bioinspir. Biomim. 5, 016001. doi:10.1088/1748-3182/5/1/016001
Zheng, X., Wang, M., Zheng, J., Tian, R., Xiong, M., and Xie, G. (2019). “Artificial Lateral Line Based Longitudinal Separation Sensing for Two Swimming Robotic Fish with Leader-Follower Formation,” in IEEE/RSJ International Conference on Intelligent Robots and Systems (IROS) (IEEE), 2539–2544.
Keywords: sensing, surface stress, swimmer interactions, fish school, computational fluid dynamics (CFD), fast fourier transform (FFT), lateral line
Citation: Li G, Kolomenskiy D, Liu H, Thiria B and Godoy-Diana R (2022) Hydrodynamical Fingerprint of a Neighbour in a Fish Lateral Line. Front. Robot. AI 9:825889. doi: 10.3389/frobt.2022.825889
Received: 30 November 2021; Accepted: 20 January 2022;
Published: 11 February 2022.
Edited by:Liang Li, Max Planck Institute of Animal Behaviour, Germany
Reviewed by:Xingwen Zheng, University of Groningen, Netherlands
Peng Zhang, New York University, United States
Copyright © 2022 Li, Kolomenskiy, Liu, Thiria and Godoy-Diana. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Gen Li, email@example.com