Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
Biomimetic Nanocarrier Targeting Drug(s) to Upstream-Receptor Mechanisms in Dementia: Focusing on Linking Pathogenic Cascades
Next Article in Special Issue
Flow Interactions Between Low Aspect Ratio Hydrofoils in In-line and Staggered Arrangements
Previous Article in Journal / Special Issue
The Ground Effect in Anguilliform Swimming
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Optimal Flow Sensing for Schooling Swimmers

by
Pascal Weber
1,
Georgios Arampatzis
1,2,
Guido Novati
1,
Siddhartha Verma
3,4,
Costas Papadimitriou
5 and
Petros Koumoutsakos
1,2,*
1
Computational Science and Engineering Laboratory, ETH Zürich, Clausiusstrasse 33, 8092 Zürich, Switzerland
2
Collegium Helveticum, 8092 Zurich, Switzerland
3
Department of Ocean and Mechanical Engineering, Florida Atlantic University, Boca Raton, FL 33431, USA
4
Harbor Branch Oceanographic Institute, Florida Atlantic University, Fort Pierce, FL 34946, USA
5
Department of Mechanical Engineering, University of Thessaly, Pedion Areos, GR-38334 Volos, Greece
*
Author to whom correspondence should be addressed.
Biomimetics 2020, 5(1), 10; https://doi.org/10.3390/biomimetics5010010
Submission received: 22 October 2019 / Revised: 17 February 2020 / Accepted: 26 February 2020 / Published: 9 March 2020
(This article belongs to the Special Issue Fluid Dynamic Interactions in Biological and Bioinspired Propulsion)

Abstract

:
Fish schooling implies an awareness of the swimmers for their companions. In flow mediated environments, in addition to visual cues, pressure and shear sensors on the fish body are critical for providing quantitative information that assists the quantification of proximity to other fish. Here we examine the distribution of sensors on the surface of an artificial swimmer so that it can optimally identify a leading group of swimmers. We employ Bayesian experimental design coupled with numerical simulations of the two-dimensional Navier Stokes equations for multiple self-propelled swimmers. The follower tracks the school using information from its own surface pressure and shear stress. We demonstrate that the optimal sensor distribution of the follower is qualitatively similar to the distribution of neuromasts on fish. Our results show that it is possible to identify accurately the center of mass and the number of the leading swimmers using surface only information.

1. Introduction

Fish navigate in their habitats by processing visual and hydrodynamic cues from their aqueous environment. Such cues may serve to provide awareness of their neighbors as fish adapt their swimming gaits in groups. Early studies have shown that vision is a critical factor for fish schooling [1]. However, more recent studies have shown that even blinded fish can keep station in a school [2]. Such capabilities are of particular importance in flow environments where vision capabilities may be limited [3]. The flow environment is replete with mechanical disturbances (pressure, shear) that can convey information about the sources that generated them. Fish swimming in groups have been found to process such hydrodynamic cues and balance them with social interactions [4,5]. In order to detect mechanical disturbances in terms of surface pressure and shear stresses fish have developed a specialized organ, the lateral line system. The mechanoreceptors in the lateral line—allowing the sensing of the disturbances in water—are called neuromasts. A number of studies and experiments have shown that the functioning of the lateral line is crucial for several tasks [6,7]. Experiments with trout in the vicinity of objects have shown its importance for Kármán gaiting and bow wake swimming as well as energy efficient station keeping [8,9]. Using the information contained in the flow, the cylinder diameter, the flow velocity, and the position relative to the generated Kármán vortex street were quantified [10,11]. Using blind cave fish, several studies have shown the importance of the lateral line to detect the location and the shape of surrounding objects and avoid obstacles [12,13,14,15]. In another study, the feeding behavior of blinded mottled sculpin was tested and it was found that they use their lateral line system to detect prey [16]. It was also found that blind fish manage to keep their position in schools and lose this ability with a disabled lateral line organ [17]. The importance of the lateral line was also shown for enhanced communication [18], the selection of habitats [19] and rheotaxis [20].
In this work, we mimic the mechanosensory receptors, more specifically the sub-surface ‘canal’ neuromasts and superficial neuromasts [21,22]. The neuromast on the fish skin are used to detect shear stresses, where the ones residing in the lateral line canals are used to detect pressure gradients [23,24,25,26,27]. Due to the filtering nature of the canals, the detection of small hydrodynamic stimuli against background noise is improved for the subsurface neuromasts [28].
The effectiveness and versatility of the lateral line organ has yielded several bio-inspired artificial flow sensors [29,30,31,32,33]. Arranging these sensors in arrays on artificial swimmers has attracted attention to transform underwater sensing [3,34,35,36,37,38,39]. Here, leveraging the intelligent distributed sensing inspired by the lateral line showed to be effective in robots moving in aquatic environments [40,41,42,43,44,45].
In order to better use and understand the capabilities of the artificial sensors several studies regarding the information content in the flow and optimal harvesting of this information were performed: The prevalence of information on the position of a vibrating source was shown to be linearly coded in the pressure gradients measured by the subsurface neuromasts [46]. Furthermore, it was shown that the variance of the pressure gradient is correlated with the presence of lateral line canals [47]. In [48], fish robots equipped with distributed pressure sensors for flow sensing were combined with Bayesian filtering in order to estimate the flow speed, the angle of attack, and the foil camber. Other studies have focused on dipole sources in order to develop methods that extract information and optimize the parameters of the sensing devices [49,50]. In a recent study artificial neural networks were employed to classify the environment using flow-only information [51,52,53,54]. In order to find effective sensor positions weight analysis algorithms were employed [55].
Following an earlier work for detection of flow disturbances generated from single obstacles [56], we examine the optimality of the spatial distribution of sensors in a self-propelled swimmer that infers the size and the relative position of the leading school. We combine numerical simulations of the two-dimensional Navier–Stokes equation and Bayesian optimal sensor placement to examine the extraction of flow information by pressure gradients and shear stresses and the optimal positioning of associated sensors. The present work demonstrates the capability of sensing a rather complex system using information of shear and pressure. Such information is available both, to biological organisms and artificial swimmers. We remark that the present work does not aim to reproduce biological systems but rather reveal algorithms that may be applicable to robotic systems. At the same time, we find that the identified optimal sensor locations for the two-dimensional artificial swimmers have similarities to biological systems indicating common governing physical mechanisms for the hydrodynamics of natural and artificial swimmers.
The paper is organised as follows: In Section 2.1 we describe the numerical simulations and in Section 2.2 the process of Bayesian optimal experimental design. We present our results in Section 3 and conclude in Section 4.

2. Materials and Methods

2.1. Flow Simulations

The swimmers are modeled by slender deforming bodies of length L which are characterized by their half-width w ( s ) along the midline [57,58]
w ( s ) = 2 w h s s 2 , 0 s < s b , w h ( w h w t ) s s b s t s b , s b s < s t , w t L s L s t , s t s L .
A sketch of the parametrization is presented in Figure 1. Following [59], we use w h = s b = 0.04 L , s t = 0.95 L and w t = 0.01 L . The swimmers propel themselves by performing sinusoidal undulations of their midline. This motion is described by a time dependent parameterization of the curvature,
k ( s , t ) = A ( s ) sin 2 π t T p 2 π s L .
Here T p = 1 is the tail-beat period and A is the undulation amplitude which linearly increases from A ( 0 ) = 0.82 / L to A ( L ) = 5.7 / L to replicate the anguilliform swimming motion described by [60]. Given the curvature along s and a center of mass, the coordinates r ( s , t ) of the swimmer’s midline can be computed by integrating the Frenet–Serret formulas [59]. In turn, the half-width w ( s ) and the coordinates r ( s , t ) characterize the swimmer’s surface.
The flow environment is described by numerical simulations of the two-dimensional incompressible Navier–Stokes equations (NSE) in velocity-pressure ( u -p) formulation. The NSE are discretized with second order finite differences and integrated in time with explicit Euler time stepping. The fluid-structure interaction is approximated with Brinkman penalization [58,61,62] by extending the fluid velocity u inside the swimmers’ bodies and by including in the NSE a penalization term to enforce no-slip and no-through boundary conditions,
u k + 1 u k δ t = p k ( u k · ) u k + ν Δ u k + i = 1 N s λ χ i ( u s , i k u k ) .
Here, ν is the kinematic viscosity, λ = 1 / δ t is the penalization coefficient, N s is the number of swimmers, u s , i k is the velocity field imposed by swimmer i (composed of translational, rotational and undulatory motions), and χ i is its characteristic function which takes value 1 inside the body of swimmer i and value 0 outside. The characteristic function χ i is computed, given the distance of each grid-point from the surface of swimmer i, by a second-order accurate finite difference approximation of a Heaviside function [63]. The pressure field is computed by pressure-projection [58,64],
Δ p k = 1 δ t · u ˜ k 1 δ t i = 1 N s χ i · u ˜ k ,
where u ˜ k = u k ( u k · ) u k + ν Δ u k . The terms inside the summation in Equation (4) are due to the non-divergence free deformation of the swimmers.

2.1.1. Schooling Formation

The tail-beating motion that propels forward a single swimmer generates in its wake a sequence of vortices. The momentum contained in the flow field induces forces which swimmers in schooling formation must overcome to maintain their positions in the group [65]. In this study, we maintain the schooling formation for multiple swimmers by employing closed-loop parametric controllers. The tail beating frequency T p , i of each swimmer i is increased or decreased if it lags behind or surpasses respectively a desired position Δ x i in the direction of the school’s motion,
T p , i = T p ( 1 Δ x i ) .
The mean school trajectory is adjusted by imposing an additional uniform curvature k C , i along each swimmer’s midline in order to minimize its lateral deviation Δ y i and its angular deflection Δ θ i ,
k C , i = [ Δ y i , Δ θ i ] + [ Δ y i , Δ θ i ] + [ Δ y i , Δ θ i ] .
Here, · defines an exponential moving average with weight δ t / T p , which approximates the integral term found in PI controllers and
[ a , b ] = | a | b , if a b < 0 , 0 , otherwise .
The formulation in Equation (6) indicates that if both the lateral displacement and the angular deviation are positive (or both negative) the swimmer will gradually revert to its position in the formation. Conversely, if Δ y i and Δ θ i have different signs the displacement has to be corrected by adding (or subtracting) curvature to the swimmer’s midline.

2.1.2. Flow Sensors

We distinguish two types of sensors on the swimmer body. The superficial neuromasts detect flow stresses and the subcanal neuromasts pressure gradients [31,66,67]. From the numerical solution of the 2D Navier–Stokes equation we obtain the flow velocity u = ( u , v ) and the pressure p at every point of the computational grid. The surface values of these quantities are obtained through a bi-linear interpolation from the nearest grid points. We perform offline analysis by recording the interpolated pressure p and flow velocity u in the vicinity of the body. We remark that we have neglected points near the end of the body to reduce the influence of large flow gradients that are generated by the motion and sharp geometry of the tail. The shear stresses are computed on the body surface using the local tangential velocity in the two nearest grid points. Moreover, we compute pressure gradients along the surface by first smoothing these pressure along the surface using splines implemented in SCIPY [68,69].

2.2. Optimal Sensor Placement Based on Information Gain

In the present work, a swimmer is equipped with sensors that are used to identify the size and location of a nearby school. The optimal sensor locations are identified using Bayesian experimental design [70] so that the information obtained from the collected measurements is maximized. We define the information gain as the distance between the prior belief on the quantities of interest and the posterior belief after obtaining the measurements. Here, we choose as measure of the distance the Kullback–Leibler divergence between the prior and the posterior distribution.

2.2.1. Bayesian Estimation of Swimmers

In the present experiment setup, we consider a group of swimmers followed by a single swimmer. The follower needs to identify (i) the relative location r of the center of mass and (ii) the population n f of the leading group. We denote with ϑ = r or ϑ = n f these unknown quantities and allow the follower to update its prior belief p ( ϑ ) about the leading group of swimmers by collecting measurements on its sensors. These sensors are distributed symmetrically on both sides of the swimmer and are represented by a single point on its mid-line. We denote the k-th measurement location at the upper and the lower part with x 1 ( s k ) and x 2 ( s k ) , respectively. The corresponding measurements are denoted by y k 1 and y k 2 , respectively (see Figure 2 for a sketch of the setup).
We denote by F ( ϑ ; s ) R 2 n the output of the flow simulation and include an error term ε to account for inaccuracies such as as numerical errors and imperfections in the sensors. The measurements on the swimmer body can be expressed as,
y = F ( ϑ ; s ) + ε .
We model the error term by a multivariate Gaussian distribution ε N ( 0 , Σ ( s ) ) with zero mean and covariance matrix Σ ( s ) R 2 n × 2 n . In this case the likelihood of a measurement is given by,
p y | ϑ , s = 1 ( 2 π ) 2 n det ( Σ ( s ) ) exp 1 2 y F ( ϑ ; s ) Σ 1 ( s ) y F ( ϑ ; s ) .
The covariance matrix depends on the sensor positions s and we assume that the prediction errors are correlated for measurements on the same side of the swimmer and uncorrelated if they originate from opposite sides. Finally, we assume that the correlation is decaying exponentially with the distance of the measurement locations. The functional form of the resulting covariance matrix is given by,
Σ i j ( s ) = σ 2 exp x 1 ( s i ) x 1 ( s j ) , if 1 i , j n , σ 2 exp x 2 ( s i n ) x 2 ( s j n ) , if n < i , j 2 n , 0 , otherwise ,
where > 0 is the correlation length and σ is the correlation strength. For all the cases described in this work, the correlation length is set to one tenth of the swimmer length = 0.1 L . The correlation strength is set to be two times the average of the signals coming from the simulations,
σ = 1 n N ϑ j = 1 2 n i = 1 N ϑ | F ( ϑ ( i ) ; s j ) | ,
where ϑ ( i ) are samples from the distribution p ( ϑ ) . We remark that the covariance matrix must be symmetric and positive definite. To ensure positive definiteness we have to take special care to the case where we pick a sensor location twice. Notice that when s i = s j for i j , a non-diagonal entry equals the diagonal entry and positive definiteness is violated. We handle this case by setting the argument of the exponential in Equation (10) to 10 7 when s i = s j . This form of the correlation error reduces the utility when sensors are placed too close together and prevents excessive clustering of the sensors [71,72].
We wish to identify the locations s yielding the largest information gain about the unknown parameter ϑ of the disturbance. A measure for information gain is defined through the Kullback–Leibler (KL) divergence between the prior belief of the parameter values and the posterior belief, i.e., after measuring the environment. The prior and posterior beliefs are represented through the density functions p ( ϑ ) and p ( ϑ | y , s ) , respectively. We denote by T the support of p ( ϑ ) . The two densities are connected through Bayes’ theorem,
p ( ϑ | y , s ) = p ( y | ϑ , s ) p ( ϑ | s ) p ( y | s ) ,
where p ( y | ϑ , s ) is the likelihood function defined in Equation (9) and p ( y | s ) is the normalization constant. We assume that the prior belief on the parameters ϑ does not depend on the sensor locations, p ( ϑ | s ) p ( ϑ ) .
The utility function is defined as [73],
u ( s , y ) : = D KL ( p ( ϑ | y , s ) | | p ( ϑ ) ) = T ln p ( ϑ | y , s ) p ( ϑ ) p ( ϑ | y , s ) ϑ .
The expected utility is defined as the average value over all possible measurements,
U ( s ) : = E y | s u ( s , y ) = Y u ( s , y ) p ( y | s ) d y = Y T ln p ( ϑ | y , s ) p ( ϑ ) p ( ϑ | y , s ) d ϑ p ( y | s ) d y ,
where Y is the domain of all possible measurements. Using Equation (12) the expected utility can be expressed as,
U ( s ) = Y T ln p ( y | ϑ , s ) p ( y | s ) p ( y | ϑ , s ) p ( ϑ ) d ϑ d y .

2.2.2. Estimated Expected Utility for Continuous Random Variables: School Location

When ϑ = r is a continuous random variable and ϑ Ω R 2 . The estimator for the expected utility in this case can be obtained by approximating the two integrals by Monte Carlo integration using N ϑ samples from p ( ϑ ) and N y samples from p ( y | ϑ , s ) [70]. The resulting estimator is given by,
U ( s ) U ^ ( s ) = 1 N ϑ N y j = 1 N y i = 1 N ϑ ln p ( y ( i , j ) | ϑ ( i ) , s ) ln 1 N ϑ k = 1 N ϑ p ( y ( i , j ) | ϑ ( k ) , s ) ,
where ϑ ( i ) p ϑ ( · ) for i = 1 , , N ϑ and y ( i , j ) p y ( · | ϑ ( i ) , s ) for j = 1 , , N y . We remark that the computational complexity of this procedure is mainly determined by the number of Navier–Stokes simulations N ϑ . There is no additional computational burden to compute the N y samples following the measurement error model in Equation (8).

2.2.3. Estimated Expected Utility for Discrete Random Variables: School Size

When ϑ is a discrete random variable with finite support taking values in the set { ϑ 1 , , ϑ N ϑ } the expected utility in Equation (15) is given by,
U ( s ) = i = 1 N ϑ p ( ϑ i ) Y ln p ( y | ϑ i , s ) p ( y | s ) p ( y | ϑ i , s ) d y .
Here, ϑ = n f represents the number of swimmers in the leading group. An estimator of the given utility can be obtained by Monte Carlo integration using N y samples from the likelihood distribution p ( y | ϑ i , s ) . The estimator is given by
U ( s ) U ^ ( s ) = 1 N y j = 1 N y i = 1 N ϑ p ( ϑ i ) ln p ( y ( i , j ) | ϑ i , s ) ln k = 1 N ϑ p ( ϑ k ) p ( y ( i , j ) | ϑ k , s ) .
where y ( i , j ) p y ( · | ϑ ( i ) , s ) for j = 1 , , N y . Let φ be the random variable representing one of the group configurations. Each group configuration is associated with a unique number φ i , for = 1 , , n i , where n i is the total number of configurations containing i swimmers. With this notation, φ takes values in the set { φ i , | i = 1 , , 8 , = 1 , , n i } . For examples of different configurations see Appendix A.
Using the fact that for i = 1 , , N ϑ ,
p ( y , φ = φ k , | ϑ = ϑ i , s ) = 0 , for k i ,
and
p ( y | ϑ = ϑ i , φ = φ i , , s ) = p ( y | φ = φ i , , s ) , for = 1 , , n i ,
and the assumption
p ( φ = φ i , | ϑ = ϑ i , s ) = 1 n i , for = 1 , , n i ,
the likelihood function can be written as,
p ( y | ϑ = ϑ i , s ) = k = 1 N ϑ = 1 n i p ( y , φ = φ k , | ϑ = ϑ i , s ) = = 1 n i p ( y , φ = φ i , | ϑ = ϑ i , s ) = = 1 n i p ( y | ϑ = ϑ i , φ = φ i , , s ) p ( φ = φ i , | ϑ = ϑ i ) = 1 n i = 1 n i p ( y | φ = φ i , , s ) .
Notice that the likelihood function for fixed ϑ i , is a mixture of Gaussian distributions with equal weights and that p ( y | φ = φ i , , s ) = N ( y | F ( φ i , ; s ) , Σ ( s ) ) . In order to draw a sample from the likelihood, first we draw an integer * with equal probability from 1 to n i and then draw y p y ( · | φ i , * , s ) .
The final form of the estimator is given by
U ^ ( s ) = 1 N y j = 1 N y i = 1 N ϑ p ( ϑ i ) ln 1 n i = 1 n i p ( y ( i , j ) | φ i , , s ) ln 1 n i k = 1 N ϑ p ( ϑ ( k ) ) = 1 n i p ( y ( i , j ) | φ i , , s ) .

2.2.4. Optimization of the Expected Utility Function

In order to determine the optimal sensor arrangement we maximize the utility estimator U ^ ( s ) described in Equation (16). It has been observed that the expected utility for many sensors often exhibit many local optima [71,74]. Heuristic approaches, such as the sequential sensor placement algorithm described by [75], have been demonstrated to be effective alternatives. Here, following [75], we perform the optimization iteratively, placing one sensor after the other. We start by placing one sensor s 1 by a grid search in the interval [ 0 , L ] , where L is the length of the swimmer. In the next step we compute the location of the second sensor by setting s = ( s 1 , s ) and repeating the grid search for the new optimal location s 2 . This procedure is then continued by defining
s i = argmax s U ^ ( s ) where s = ( s 1 , , s i 1 , s ) .
We note that the scalar variable s denotes the mid-line coordinate of a single sensor-pair, whereas the vector s holds the mid-line coordinates of all sensor-pairs. Besides the mentioned advantages, sequential placement allows to quantify the importance of each sensor placed and provides further insight into the resulting distribution of sensors.

3. Results

We examine the optimal arrangement of pressure gradient and shear stress sensors on the surface of a swimmer trailing a school of self-propelled swimmers. We consider two sensing objectives: (a) the size of the leading school and (b) the relative position of the school. The simulations correspond to a Reynolds number Re = L 2 ν = 2000 . In all experiments, we use 4096 points to discretize the horizontal direction x [ 0 , 1 ] and all artificial swimmers have a length of L = 0 . 1 .
For the “size of the leading school” experiment, where the aim is to determine the size of the group, we chose the school-sizes to be ϑ i = 1 , , 8 . First we consider one configuration per group-size. In this case inferring the configuration is equivalent to inferring the number of swimmer in the group. To increase the difficulty we consider n i different initial configurations. In each configuration we assign a number φ i , for i = 1 , , 8 and = 1 , , n i . In total, we consider N t o t = i n i = 61 distinct configurations each having the same prior probability 1 / N t o t . In Appendix A we present the initial condition for all configurations. The center of mass of the school is located at x = 0.3 and in the y-axis in the middle of the vertical extent of the domain. We use a controller to fix the distance between x and y coordinates of two swimmers to Δ x = Δ y = 0.15 , see Section 2.1.1.
For the “relative position” experiment, where the aim is to determine the relative location of the follower to the center of mass of the leading group, we consider three independent experiments with one, four and seven leading swimmers. Snapshots of the pressure field for these simulations are presented in Figure 3. The prior probability for the position of the group is uniform in the domain [ 0.6 , 0.8 ] × [ 0.1 , 0.4 ] . The support of the prior probability is discretized with 21 × 31 gridpoints. Since the experiments are independent, the total expected utility function for the three cases is the sum of the expected utility of each experiment [56].
For both experiments we record the pressure gradient and shear stress on the surface of the swimmer using the methods discussed in Section 2.1.2. The motion of the swimmer introduces disturbances on its own surface. In order to distinguish the self-induced from the environment disturbances we freeze the movement of the following swimmer and set its curvature to zero. The freezing time is selected by evolving the simulation until the wakes of the leading group are sufficiently mixed and passed the following swimmer. We found that this is the case for T = 22 . The transition from swimming to coasting motion takes place during the time interval [ T , T + 1 ] . Finally, we record the pressure gradient and the shear stress at time T + 2 . The resulting sensor-signal associated to the midline coordinates s for a given configuration ϑ is denoted F ( ϑ ; s ) , see Equation (8).

3.1. Utility Function for the First Sensor

In this section we discuss the optimal location of a single pressure gradient sensor using the estimators in Equations (16) and (20). Recall that we estimate the expected KL divergence between the prior and the posterior distribution for different sensor locations s . The KL divergence can be understood as a measure of distance between two probability distributions. Thus, higher values of divergence correspond to preferable locations for the sensor, leading to higher information gain. The resulting utilities are plotted in Figure 4. For all experiments we find that the tip of the head ( s = 0 ) exhibits the largest utility independent of the number of swimmer in the leading group.
At the tip of the head, the two symmetrically placed sensors have the smallest distance. In Equation (10) we have assumed that the two swimmer halves are symmetric and uncorrelated. Due to the small distance of the sensors at the head, spatial correlation between the sensors across the swimmer halves would decrease the utility of this location. In order to test whether the utility for sensors at the head is influenced by this symmetry assumption, we perform experiments where we place a single sensor on one side of the swimmer. Again, in this case the location at the head is found to have the highest expected utility.
There is evidence that the head experiences the largest variance of pressure gradients F ( ϑ ; s ) . The same observations can be made for the density of the sub-canal neuromasts, which is also highest in the front of the fish [47]. To check the presence of this correlation in our study, we examine the variance of the values obtained from our numerical solution of the Navier–Stokes equation. We confirm that our simulations are consistent with this experimental observation. We find that independent of the number of swimmers, the variance in the sensor signal var ϑ ( F ( ϑ ; s ) ) is largest at s = 0 .

3.2. Sequential Sensor Placement

In this section we discuss the results of the sequential sensor placement described in Section 2.2.4. For the “size of the leading school” experiment we present the results in Figure 5. In Figure 5a the utility curve for the first five sensors is shown. We observe that the utility curve becomes flatter as the number of sensors increase. Furthermore, we observe that the location where the previous sensor was placed is a minimum for the utility for the next sensor. Figure 5b shows the utility estimator at the optimal sensor for up to 20 sensors and it is evident that the value of the expected utility reaches a plateau. In Figure 5c the found optimal location of the sensors on the skin of the swimmer is presented. The numbers correspond to the iteration in the sequential procedure that the sensor was placed. Note that the sensors are being placed symmetrically.
The optimal sensor placement results for the “relative position” experiment can be found in Figure 6. Similar to the other experiment the utility curves become flatter after every placed sensor and the location for the previous sensor is a minimum for the utility for the next sensor (see Figure 6a). We plot the maximum of the utility for up to 20 sensors (see Figure 6b) and observe a convergence to a constant value. In Figure 6c the found optimal location of the first 20 sensors is presented.
For both experiments, it is evident that the utility of the optimal sensor location approaches a constant value. This fact can be explained by recalling that the expected utility in Equation (15) is a measure of the averaged distance between the prior and the posterior distribution. Increasing the number of sensors leads to an increase in the number of measurements. By the Bayesian central limit theorem, increasing the number of measurements leads to convergence of the posterior to a Dirac distribution. As soon as the posterior has converged, the expected distance from the prior, and thus the expected utility, remains constant.
The found sensor distributions for the two objectives are similar, having clusters at the head and uniform distribution along the body. In order to underpin the biological relevance of the observed sensor distribution we compare our results to [47]. Given that the canals display significant 3D branching in the head a direct comparison is difficult. However, the found cluster of sensors at the head agrees qualitatively with the high canal density reported in [47].

3.3. Inference of the Environment

In this section we demonstrate the importance of the optimal sensor locations and examine the convergence of the posterior distribution. We compute the posterior distribution via Bayes’ theorem given in Equation (12). We set y = F ( ϑ , s ) and compute the posterior for different values of ϑ in the prior region. We consider measurements collected at: (a) the optimal and (b) the worst sensors location.
The posterior probability for the “size of the leading school” experiment is shown in Figure 7. We observe that the worst sensor location implies an almost uniform posterior distribution, reflecting that measurements at this sensor carry no information. On the other hand, the posterior distribution for the optimal sensor is more informative. We observe that for groups with small size the follower is able to identify the size with more confidence, as opposed to larger groups. We compare the posterior for an experiment with only one configuration per group-size to an experiment with multiple configurations. For multiple configurations the posterior is less informative. This indicates that the second case occurs to be a more difficult problem. Finally, notice that the posterior for one configuration is symmetric, where when adding multiple configurations this symmetry is broken. This fact is discussed in Appendix B.
The posterior density for the “relative position” with one leading swimmer is presented in Figure 8. The posterior for the configuration with three and seven swimmers is similar. We compute the posterior for measurements at the best and the worst location for one and three sensors. For the three sensors the worst location has been selected in all three phases of the sequential placement. The results for the normalized densities are shown in Figure 8. We observe that one sensor at the optimal location gives a very peaked posterior. Three optimal sensors can infer the location with low uncertainty. This is not the case for the worst sensors, where adding more sensors does not immediately lead to uncertainty reduction.

3.4. Shear Stress Sensors

In this section, we discuss the results for the optimal positioning of shear stress sensors. We follow the same procedure as in Section 3.1 and Section 3.2. Here, we omit the presentation of all the results and focus on the similarities and differences to the pressure gradient sensors.
The optimal location for a single sensor for the “size of the leading school” experiment is at s * = 3.01 × 10 4 . For the “relative position” experiment we find the optimal location s * = 3.84 × 10 4 . In contrast to the optimal location for one pressure gradient sensor, the found sensor is not at the tip of the head and is at different positions for the two experiments. Examining the variance in the shear signal shows quantitatively the same behaviour as the utility. Comparing the location of the maxima in variance shows that they do not coincide with the found maxima for the expected utility for shear sensors.
We perform sequential placement of 15 sensors. The resulting distribution of sensors is shown in Figure 9. In Section 3.2 we argue that the expected utility must reach a plateau when placing many sensors using the Bayesian central limit theorem. For shear stress sensors we observe that the convergence is slower compared to the pressure gradient sensors. We conclude that the information gain per shear stress sensor placed is lower as for the pressure gradient sensors.
The posterior density obtained for both experiments is less informative when using the same number of sensors. Also this indicates that shear is a less informative quantity yielding a slower convergence of the posterior. This is in agreement with the observation that the subcanal neuromasts associated with pressure gradient sensing are more robust to noise [28]. For multiple fish in schools the resulting flow field is disturbed, thus suggesting the use of pressure gradient sensors.

4. Discussion

We present a study of the optimal sensor locations on a self-propelled swimmer for detecting the size and location of a leading group of swimmers. This optimization combines Bayesian experimental design with large scale simulations of the two dimensional Navier–Stokes equations. Mimicking the function of sensory organs in real fish, we used the shear stress and pressure gradient on the surface of the swimmers to determine the sensor feedback generated by a disturbance in the flow field.
The optimization was performed for different configurations of swimmers, ranging from a simple leader-follower configuration with two swimmers, to a group of up to eight swimmers leading a single follower. We considered two types of information: the number of swimmers in the leading group and the relative location of the leading group. We find that, although the general shape of the utility function varies between the two objectives, the preferred location of the first sensor on the head of the swimmer is consistent. Furthermore, we find that the objective is only weakly influenced when varying the number of members in the leading group.
We perform a sequential sensor placement and find that the utility converges to a constant value and thus we can conclude that few sensors suffice to infer the quantities of the surrounding flow. Indeed, we find that the optimal sensor locations correspond to a posterior distribution that is strongly peaked around the true value of the quantity of interest. In summary, we find, that for the group sizes under examination, changing the number of swimmers in the leading group does not influence the follower’s ability to infer the mean school location. Furthermore, we were able to show that choosing the locations for the measurements in a systematic way we are able to infer the number of swimmer in the leading group and the location of our agent to high accuracy.
We envision that the presented methodology can provide guidance in developing autonomous systems of schooling artificial swimmers. While biological organisms have distinct flow fields from those examined in the present two-dimensional simulations, we believe that the algorithms presented here in can be extended to 3D flows. Moreover, while we draw a distinction between fish and the studied artificial swimmers, we note the capability of identifying neighboring swimmers using shear and pressure information on the body of the swimmers, indicating sufficiency of such type of information for flow sensing.

Author Contributions

Conceptualization, G.A. and P.K.; Data curation, P.W.; Formal analysis, P.W. and G.A.; Funding acquisition, P.K.; Investigation, P.W.; Methodology, G.A., S.V., C.P. and P.K.; Project administration, P.K.; Resources, P.K.; Software, P.W. and G.N.; Supervision, G.A. and P.K.; Validation, P.W. and G.N.; Visualization, P.W. and G.A.; Writing—original draft, P.W.; Writing—review & editing, P.W., G.A., G.N., S.V., C.P. and P.K. All authors have read and agreed to the published version of the manuscript.

Funding

We would like to acknowledge the computational time at Swiss National Supercomputing Center (CSCS) under the project s929. We gratefully acknowledge support from the European Research Council (ERC) Advanced Investigator Award (No. 341117).

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Configurations

The configuration used for the “size of the leading school” experiment. For the configurations with three rows the vertical extent y [ 0 , 0.5 ] was discretized using 2048 gridpoints, for the ones with four rows it was extended to y [ 0 , 0.75 ] and discretized using 3072 gridpoints.
Figure A1. Configurations for two leading swimmers.
Figure A1. Configurations for two leading swimmers.
Biomimetics 05 00010 g0a1
Figure A2. Configurations for three leading swimmers.
Figure A2. Configurations for three leading swimmers.
Biomimetics 05 00010 g0a2
Figure A3. Configurations for four leading swimmers.
Figure A3. Configurations for four leading swimmers.
Biomimetics 05 00010 g0a3
Figure A4. Configurations for five leading swimmers.
Figure A4. Configurations for five leading swimmers.
Biomimetics 05 00010 g0a4
Figure A5. Configurations for six leading swimmers.
Figure A5. Configurations for six leading swimmers.
Biomimetics 05 00010 g0a5
Figure A6. Configurations for seven leading swimmers.
Figure A6. Configurations for seven leading swimmers.
Biomimetics 05 00010 g0a6
Figure A7. Configurations for eight leading swimmers.
Figure A7. Configurations for eight leading swimmers.
Biomimetics 05 00010 g0a7

Appendix B. The Posterior Is Not Symmetric

The estimated posterior in Figure 7b is not symmetric with respect to the ϑ true = ϑ diagonal. This observation indicates that the posterior is not symmetric with respect to an exchange of ϑ and ϑ true , the parameter we try to infer and the one used in the simulation. Here, we want to show that this observation is true in general. In order to lighten the notation, we neglect the dependence of the distributions on the sensor location s .
In Section 2.2.3 we showed that the distribution of ϑ i conditioned on measurements y , under the assumption of uniform prior, is proportional to
p ( ϑ i | y ) p ( y | ϑ i ) p ( ϑ i ) = 1 8 1 n i = 1 n i 1 2 π σ 2 exp ( y F ( φ i , ) ) 2 2 σ 2 .
We want to show that for any i j ,
p ( ϑ i | ϑ j ) = p ( ϑ i | y = F ( φ j , k ) ) p ( ϑ j | y = F ( φ i , ) ) = p ( ϑ j | ϑ i ) ,
for any configurations φ j , k and φ i , corresponding to a school of size ϑ j and ϑ i , respectively. From Equation (A1) it is easy to see that (A2) is true due to the fact that
1 n i = 1 n i exp ( F ( φ j , k ) F ( φ i , ) ) 2 2 σ 2 1 n j k = 1 n j exp ( F ( φ i , ) F ( φ j , k ) ) 2 2 σ 2 .
Finally, we note that in the case where we have only one configuration per group size, i.e., n i = 1 for all i, the statement in (A3) is not true and the posterior is symmetric.

References

  1. Morrow, J.E. Schooling Behavior in Fishes. Q. Rev. Biol. 1948, 23, 27–38. [Google Scholar] [CrossRef] [PubMed]
  2. Partridge, B.L.; Pitcher, T. The sensory basis of fish schools: Relative roles of lateral line and vision. J. Comp. Physiol. 1980, 135, 315–325. [Google Scholar] [CrossRef]
  3. Triantafyllou, M.S.; Weymouth, G.D.; Miao, J. Biomimetic Survival Hydrodynamics and Flow Sensing. Annu. Rev. Fluid Mech. 2016, 48, 1–24. [Google Scholar] [CrossRef] [Green Version]
  4. Ward, A.J.W.; Sumpter, D.J.T.; Couzin, I.D.; Hart, P.J.B.; Krause, J. Quorum decision-making facilitates information transfer in fish shoals. Proc. Natl. Acad. Sci. USA 2008, 105, 6948–6953. [Google Scholar] [CrossRef] [Green Version]
  5. Puckett, J.G.; Pokhrel, A.R.; Giannini, J.A. Collective gradient sensing in fish schools. Sci. Rep. 2018, 8, 7587. [Google Scholar] [CrossRef]
  6. Dykgraaf, S. Untersuchungen über die Funktion der Seitenorgane an Fischen. Zeitschrift für Vergleichende Physiologie 1933, 20, 162–214. [Google Scholar] [CrossRef]
  7. Dykgraaf, S. The functioning and significance of the lateral-line organs. Biol. Rev. Camb. Philos. Soc. 1963, 38, 51–105. [Google Scholar] [CrossRef]
  8. Bleckmann, H.; Przybilla, A.; Klein, A.; Schmitz, A.; Kunze, S.; Brücker, C. Station Holding of Trout: Behavior, Physiology and Hydrodynamics. In Nature-Inspired Fluid Mechanics: Results of the DFG Priority Programme 1207 ”Nature-Inspired Fluid Mechanics” 2006–2012; Springer: Berlin/Heidelberg, Germany, 2012; pp. 161–177. [Google Scholar] [CrossRef]
  9. Sutterlin, A.; Waddy, S. Possible Role of the Posterior Lateral Line in Obstacle Entrainment by Brook Trout (Salvelinus fontinalis). J. Fish. Res. Board Can. 2011, 32, 2441–2446. [Google Scholar] [CrossRef]
  10. Akanyeti, O.; Venturelli, R.; Visentin, F.; Chambers, L.; Megill, W.M.; Fiorini, P. What information do Kármán streets offer to flow sensing? Bioinspir. Biomim. 2011, 6, 036001. [Google Scholar] [CrossRef]
  11. Chambers, L.D.; Akanyeti, O.; Venturelli, R.; Ježov, J.; Brown, J.; Kruusmaa, M.; Fiorini, P.; Megill, W.M. A fish perspective: Detecting flow features while moving using an artificial lateral line in steady and unsteady flow. J. R. Soc. Interface 2014, 11. [Google Scholar] [CrossRef] [Green Version]
  12. von Campenhausen, C.; Riess, I.; Weissert, R. Detection of stationary objects by the blind Cave FishAnoptichthys jordani (Characidae). J. Comp. Physiol. 1981, 143, 369–374. [Google Scholar] [CrossRef]
  13. Hassan, E.S. Hydrodynamic Imaging of the Surroundings by the Lateral Line of the Blind Cave Fish Anoptichthys jordani. In The Mechanosensory Lateral Line; Coombs, S., Görner, P., Münz, H., Eds.; Springer: New York, NY, USA, 1989; pp. 217–227. [Google Scholar]
  14. Windsor, S.P.; Norris, S.E.; Cameron, S.M.; Mallinson, G.D.; Montgomery, J.C. The flow fields involved in hydrodynamic imaging by blind Mexican cave fish (Astyanax fasciatus). Part I: Open water and heading towards a wall. J. Exp. Biol. 2010, 213, 3819–3831. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Windsor, S.P.; Norris, S.E.; Cameron, S.M.; Mallinson, G.D.; Montgomery, J.C. The flow fields involved in hydrodynamic imaging by blind Mexican cave fish (Astyanax fasciatus). Part II: Gliding parallel to a wall. J. Exp. Biol. 2010, 213, 3832–3842. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Hoekstra, D.; Janssen, J. Non-visual feeding behavior of the mottled sculpin, Cottus bairdi, in Lake Michigan. Environ. Biol. Fishes 1985, 12, 111–117. [Google Scholar] [CrossRef]
  17. Pitcher, T.; Partridge, B.; Wardle, C. A blind fish can school. Science 1976, 194, 963–965. [Google Scholar] [CrossRef]
  18. Satou, M.; Takeuchi, H.A.; Nishii, J.; Tanabe, M.; Kitamura, S.; Okumoto, N.; Iwata, M. Behavioral and electrophysiological evidences that the lateral line is involved in the inter-sexual vibrational communication of the himé salmon (landlocked red salmon, Oncorhynchus nerka). J. Comp. Physiol. A 1994, 174, 539–549. [Google Scholar] [CrossRef]
  19. Huijbers, C.M.; Nagelkerken, I.; Lössbroek, P.A.C.; Schulten, I.E.; Siegenthaler, A.; Holderied, M.W.; Simpson, S.D. A test of the senses: Fish select novel habitats by responding to multiple cues. Ecology 2012, 93, 46–55. [Google Scholar] [CrossRef] [Green Version]
  20. Montgomery, J.C.; Baker, C.F.; Carton, A.G. The lateral line can mediate rheotaxis in fish. Nature 1997, 389, 960–963. [Google Scholar] [CrossRef]
  21. Coombs, S.; Janssen, J.; Webb, J.F. Diversity of lateral line systems: Evolutionary and functional considerations. In Sensory Biology of Aquatic Animals; Springer: Berlin/Heidelberg, Germany, 1988; pp. 553–593. [Google Scholar]
  22. Coombs, S.; Görner, P.; Münz, H. A Brief Overview of the Mechanosensory Lateral Line System and the Contributions to This Volume. In The Mechanosensory Lateral Line; Coombs, S., Görner, P., Münz, H., Eds.; Springer: New York, NY, USA, 1989; pp. 3–5. [Google Scholar]
  23. Denton, E.J.; Gray, J.A.B. Some Observations on the Forces Acting on Neuromasts in Fish Lateral Line Canals. In The Mechanosensory Lateral Line; Coombs, S., Görner, P., Münz, H., Eds.; Springer: New York, NY, USA, 1989; pp. 229–246. [Google Scholar]
  24. Coombs, S.; Braun, C.B. Information Processing by the Lateral Line System. In Sensory Processing in Aquatic Environments; Springer: New York, NY, USA, 2003; pp. 122–138. [Google Scholar] [CrossRef]
  25. Coombs, S.; Netten, S.V. The Hydrodynamics and Structural Mechanics of the Lateral Line System. In Fish Physiology; Elsevier: Amsterdam, The Netherlands, 2005; Volume 23, pp. 103–139. [Google Scholar] [CrossRef]
  26. Bleckmann, H. Peripheral and central processing of lateral line information. J. Comp. Physiol. A 2008, 194, 145–158. [Google Scholar] [CrossRef]
  27. Jiang, Y.; Ma, Z.; Zhang, D. Flow field perception based on the fish lateral line system. Bioinspir. Biomim. 2019, 14, 041001. [Google Scholar] [CrossRef]
  28. Engelmann, J.; Hanke, W.; Mogdans, J.; Bleckmann, H. Hydrodynamic stimuli and the fish lateral line. Nature 2000, 408, 1476–4687. [Google Scholar] [CrossRef] [PubMed]
  29. Kottapalli, A.G.P.; Asadnia, M.; Miao, J.M.; Barbastathis, G.; Triantafyllou, M.S. A flexible liquid crystal polymer MEMS pressure sensor array for fish-like underwater sensing. Smart Mater. Struct. 2012, 21, 115030. [Google Scholar] [CrossRef] [Green Version]
  30. Tao, J.; Yu, X. Hair flow sensors: From bio-inspiration to bio-mimicking—A review. Smart Mater. Struct. 2012, 21, 113001. [Google Scholar] [CrossRef]
  31. Asadnia, M.; Kottapalli, A.G.P.; Miao, J.; Warkiani, M.E.; Triantafyllou, M.S. Artificial fish skin of self-powered micro-electromechanical systems hair cells for sensing hydrodynamic flow phenomena. J. R. Soc. Interface 2015, 12. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Kottapalli, A.G.P.; Bora, M.; Sengupta, D.; Miao, J.; Triantafyllou, M.S. Hydrogel-CNT Biomimetic Cilia for Flow Sensing. In Proceedings of the 2018 IEEE SENSORS, New Delhi, India, 28–31 October 2018; pp. 1–4. [Google Scholar] [CrossRef]
  33. Wolf, B.J.; Morton, J.A.S.; MacPherson, W.N.; van Netten, S.M. Bio-inspired all-optical artificial neuromast for 2D flow sensing. Bioinspir. Biomim. 2018, 13, 026013. [Google Scholar] [CrossRef] [PubMed]
  34. Yang, Y.; Chen, J.; Engel, J.; Pandya, S.; Chen, N.; Tucker, C.; Coombs, S.; Jones, D.L.; Liu, C. Distant touch hydrodynamic imaging with an artificial lateral line. Proc. Natl. Acad. Sci. USA 2006, 103, 18891–18895. [Google Scholar] [CrossRef] [Green Version]
  35. Yang, Y.; Nguyen, N.; Chen, N.; Lockwood, M.; Tucker, C.; Hu, H.; Bleckmann, H.; Liu, C.; Jones, D.L. Artificial lateral line with biomimetic neuromasts to emulate fish sensing. Bioinspir. Biomim. 2010, 5, 016001. [Google Scholar] [CrossRef] [Green Version]
  36. Strokina, N.; Kämäräinen, J.; Tuhtan, J.A.; Fuentes-Pérez, J.F.; Kruusmaa, M. Joint Estimation of Bulk Flow Velocity and Angle Using a Lateral Line Probe. IEEE Trans. Instrum. Meas. 2016, 65, 601–613. [Google Scholar] [CrossRef]
  37. Xu, Y.; Mohseni, K. A Pressure Sensory System Inspired by the Fish Lateral Line: Hydrodynamic Force Estimation and Wall Detection. IEEE J. Ocean. Eng. 2017, 42, 532–543. [Google Scholar] [CrossRef]
  38. Sengupta, D.; Chen, S.H.; Kottapalli, A.G.P. Nature-Inspired Self-Powered Sensors and Energy Harvesters. In Self-Powered and Soft Polymer MEMS/NEMS Devices; Springer: Cham, Swizerland, 2019; pp. 61–81. [Google Scholar] [CrossRef]
  39. Zhang, X.; Shan, X.; Shen, Z.; Xie, T.; Miao, J. A New Self-Powered Sensor Using the Radial Field Piezoelectric Diaphragm in d33 Mode for Detecting Underwater Disturbances. Sensors 2019, 19, 962. [Google Scholar] [CrossRef] [Green Version]
  40. Kruusmaa, M.; Fiorini, P.; Megill, W.; de Vittorio, M.; Akanyeti, O.; Visentin, F.; Chambers, L.; El Daou, H.; Fiazza, M.; Ježov, J.; et al. FILOSE for Svenning: A Flow Sensing Bioinspired Robot. IEEE Robot. Autom. Mag. 2014, 21, 51–62. [Google Scholar] [CrossRef]
  41. DeVries, L.; Lagor, F.D.; Lei, H.; Tan, X.; Paley, D.A. Distributed flow estimation and closed-loop control of an underwater vehicle with a multi-modal artificial lateral line. Bioinspir. Biomim. 2015, 10, 025002. [Google Scholar] [CrossRef] [PubMed]
  42. Ježov, J.; Akanyeti, O.; Chambers, L.D.; Kruusmaa, M. Sensing oscillations in unsteady flow for better robotic swimming efficiency. In Proceedings of the 2012 IEEE International Conference on Systems, Man, and Cybernetics (SMC), Seoul, Korea, 14–17 October 2012; pp. 91–96. [Google Scholar] [CrossRef]
  43. Yen, W.; Sierra, D.M.; Guo, J. Controlling a Robotic Fish to Swim Along a Wall Using Hydrodynamic Pressure Feedback. IEEE J. Ocean. Eng. 2018, 43, 369–380. [Google Scholar] [CrossRef]
  44. Krieg, M.; Nelson, K.; Mohseni, K. Distributed sensing for fluid disturbance compensation and motion control of intelligent robots. Nat. Mach. Intell. 2019, 1, 216–224. [Google Scholar] [CrossRef]
  45. Zheng, X.; Wang, M.; Zheng, J.; Tian, R.; Xiong, M.; Xie, G. Artificial lateral line based longitudinal separation sensing for two swimming robotic fish with leader-follower formation. In Proceedings of the 2019 IEEE/RSJ International Conference on Intelligent Robots and Systems (IROS), Macau, China, 4–8 November 2019; pp. 2539–2544. [Google Scholar] [CrossRef]
  46. Ćurčić-Blake, B.; van Netten, S.M. Source location encoding in the fish lateral line canal. J. Exp. Biol. 2006, 209, 1548–1559. [Google Scholar] [CrossRef] [Green Version]
  47. Ristroph, L.; Liao, J.C.; Zhang, J. Lateral Line Layout Correlates with the Differential Hydrodynamic Pressure on Swimming Fish. Phys. Rev. Lett. 2015, 114, 018102. [Google Scholar] [CrossRef] [Green Version]
  48. Zhang, F.; Lagor, F.; Yeo, D.; Washington, P.; Paley, D. Distributed Flow Sensing Using Bayesian Estimation for a Flexible Fish Robot. In Proceedings of the ASME 2015 Dynamic Systems and Control Conference, Columbus, OH, USA, 28–30 October 2015. [Google Scholar] [CrossRef]
  49. Ahrari, A.; Lei, H.; Sharif, M.A.; Deb, K.; Tan, X. Design optimization of artificial lateral line system under uncertain conditions. In Proceedings of the 2015 IEEE Congress on Evolutionary Computation (CEC), Sendai, Japan, 25–28 May 2015; pp. 1807–1814. [Google Scholar] [CrossRef]
  50. Ahrari, A.; Lei, H.; Sharif, M.A.; Deb, K.; Tan, X. Reliable underwater dipole source characterization in 3D space by an optimally designed artificial lateral line system. Bioinspir. Biomim. 2017, 12, 036010. [Google Scholar] [CrossRef]
  51. Boulogne, L.H.; Wolf, B.J.; Wiering, M.A.; van Netten, S.M. Performance of neural networks for localizing moving objects with an artificial lateral line. Bioinspir. Biomim. 2017, 12, 056009. [Google Scholar] [CrossRef] [Green Version]
  52. Colvert, B.; Alsalman, M.; Kanso, E. Classifying vortex wakes using neural networks. Bioinspir. Biomim. 2018, 13, 025003. [Google Scholar] [CrossRef] [Green Version]
  53. Wolf, B.J.; Pirih, P.; Kruusmaa, M.; Van Netten, S.M. Shape Classification Using Hydrodynamic Detection via a Sparse Large-Scale 2D-Sensitive Artificial Lateral Line. IEEE Access 2020, 8, 11393–11404. [Google Scholar] [CrossRef]
  54. Wolf, B.; van de Wolfshaar, J.; van Netten, S. Three-dimensional multi-source localization of underwater objects using convolutional neural networks for artificial lateral lines. J. R. Soc. Interface 2020, 17, 20190616. [Google Scholar] [CrossRef] [PubMed]
  55. Xu, D.; Lv, Z.; Zeng, H.; Bessaih, H.; Sun, B. Sensor placement optimization in the artificial lateral line using optimal weight analysis combining feature distance and variance evaluation. ISA Trans. 2019, 86, 110–121. [Google Scholar] [CrossRef] [PubMed]
  56. Verma, S.; Papadimitriou, C.; Lüthen, N.; Arampatzis, G.; Koumoutsakos, P. Optimal sensor placement for artificial swimmers. J. Fluid Mech. 2020, 884, A24. [Google Scholar] [CrossRef] [Green Version]
  57. Kern, S.; Koumoutsakos, P. Simulations of optimized anguilliform swimming. J. Exp. Biol. 2006, 209, 4841–4857. [Google Scholar] [CrossRef] [Green Version]
  58. Gazzola, M.; Chatelain, P.; van Rees, W.M.; Koumoutsakos, P. Simulations of single and multiple swimmers with non-divergence free deforming geometries. J. Comput. Phys. 2011, 230, 7093–7114. [Google Scholar] [CrossRef]
  59. Kern, S.; Chatelain, P.; Koumoutsakos, P. Modeling, Simulation and Optimization of Anguilliform Swimmers. In Bio-Mechanisms of Swimming and Flying: Fluid Dynamics, Biomimetic Robots, and Sports Science; Springer: Berlin/Heidelberg, Germany, 2008; p. 167. [Google Scholar]
  60. Carling, J.; Williams, T.L.; Bowtell, G. Self-propelled anguilliform swimming: Simultaneous solution of the two-dimensional Navier–Stokes equations and Newton’s laws of motion. J. Exp. Biol. 1998, 201, 3143–3166. [Google Scholar]
  61. Angot, P.; Bruneau, C.H.; Fabrie, P. A penalization method to take into account obstacles in incompressible viscous flows. Numer. Math. 1999, 81, 497–520. [Google Scholar] [CrossRef]
  62. Coquerelle, M.; Cottet, G.H. A vortex level set method for the two-way coupling of an incompressible fluid with colliding rigid bodies. J. Comput. Phys. 2008, 227, 9121–9137. [Google Scholar] [CrossRef] [Green Version]
  63. Towers, J.D. Finite difference methods for approximating Heaviside functions. J. Comput. Phys. 2009, 228, 3478–3489. [Google Scholar] [CrossRef]
  64. Chorin, A.J. Numerical Solution of the Navier–Stokes Equations. Math. Comput. 1968, 22, 745–762. [Google Scholar] [CrossRef]
  65. Novati, G.; Verma, S.; Alexeev, D.; Rossinelli, D.; van Rees, W.M.; Koumoutsakos, P. Synchronised Swimming of Two Fish. arXiv 2016, arXiv:1610.04248. [Google Scholar]
  66. Kroese, A.B.; Schellart, N.A. Velocity- and acceleration-sensitive units in the trunk lateral line of the trout. J. Neurophysiol. 1992, 68, 2212–2221. [Google Scholar] [CrossRef] [PubMed]
  67. Bleckmann, H.; Zelick, R. Lateral line system of fish. Integr. Zool. 2009, 4, 13–25. [Google Scholar] [CrossRef] [PubMed]
  68. Virtanen, P.; Gommers, R.; Oliphant, T.E.; Haberl, M.; Reddy, T.; Cournapeau, D.; Burovski, E.; Peterson, P.; Weckesser, W.; Bright, J.; et al. SciPy 1.0: Fundamental Algorithms for Scientific Computing in Python. Nat. Methods 2020. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  69. Dierckx, P. An algorithm for smoothing, differentiation and integration of experimental data using spline functions. J. Comput. Appl. Math. 1975, 1, 165–184. [Google Scholar] [CrossRef] [Green Version]
  70. Huan, X.; Marzouk, Y.M. Simulation-based optimal Bayesian experimental design for nonlinear systems. J. Comput. Phys. 2013, 232, 288–317. [Google Scholar] [CrossRef] [Green Version]
  71. Papadimitriou, C.; Lombaert, G. The effect of prediction error correlation on optimal sensor placement in structural dynamics. Mech. Syst. Signal Process. 2012, 28, 105–127. [Google Scholar] [CrossRef]
  72. Simoen, E.; Papadimitriou, C.; Lombaert, G. On prediction error correlation in Bayesian model updating. J. Sound Vib. 2013, 332, 4136–4152. [Google Scholar] [CrossRef]
  73. Ryan, K.J. Estimating Expected Information Gains for Experimental Designs With Application to the Random Fatigue-Limit Model. J. Comput. Graph. Stat. 2003, 12, 585–603. [Google Scholar] [CrossRef]
  74. Papadimitriou, D.I.; Papadimitriou, C. Optimal sensor placement for the estimation of turbulence model parameters in CFD. Int. J. Uncertain. Quant. 2015, 5, 545–568. [Google Scholar] [CrossRef] [Green Version]
  75. Papadimitriou, C. Optimal sensor placement methodology for parametric identification of structural systems. J. Sound Vib. 2004, 278, 923–947. [Google Scholar] [CrossRef]
Figure 1. Parametrization of the swimmer surface as described in Equation (1).
Figure 1. Parametrization of the swimmer surface as described in Equation (1).
Biomimetics 05 00010 g001
Figure 2. Simulation setup used for determining the optimal sensor distribution on a fish-like body. The follower is initially located inside the rectangular area. The number of swimmers in the leading group is varied between one and eight. The sensor-placement algorithm attempts to find the arrangement of sensors s that allows the follower to determine with lowest uncertainty the relative position r and the number of swimmers n f in the leading group of swimmers. For each sensor s k the swimmer collects measurements y k 1 and y k 2 at locations x 1 ( s k ) and x 2 ( s k ) on the skin, respectively.
Figure 2. Simulation setup used for determining the optimal sensor distribution on a fish-like body. The follower is initially located inside the rectangular area. The number of swimmers in the leading group is varied between one and eight. The sensor-placement algorithm attempts to find the arrangement of sensors s that allows the follower to determine with lowest uncertainty the relative position r and the number of swimmers n f in the leading group of swimmers. For each sensor s k the swimmer collects measurements y k 1 and y k 2 at locations x 1 ( s k ) and x 2 ( s k ) on the skin, respectively.
Biomimetics 05 00010 g002
Figure 3. Snapshots of the pressure field in the environment of the follower swimmer generated by one (a), four (b) and seven (c) schooling swimmers. The snapshots are taken at the moment the measurement was performed for one particular location of the follower in the prior region. High pressure is shown in red and low pressure in blue.
Figure 3. Snapshots of the pressure field in the environment of the follower swimmer generated by one (a), four (b) and seven (c) schooling swimmers. The snapshots are taken at the moment the measurement was performed for one particular location of the follower in the prior region. High pressure is shown in red and low pressure in blue.
Biomimetics 05 00010 g003
Figure 4. Utility curves for the first sensor using pressure measurements. In (a) the utility estimator for the “size of the leading school” experiment is presented. (b) corresponds to the utility estimator for the “relative position” experiment. We show the resulting curves for one, three and seven swimmer in the leading group and the total expected utility. We observe that although the form does not drastically change, the total utility increases with increasing size of the leading group.
Figure 4. Utility curves for the first sensor using pressure measurements. In (a) the utility estimator for the “size of the leading school” experiment is presented. (b) corresponds to the utility estimator for the “relative position” experiment. We show the resulting curves for one, three and seven swimmer in the leading group and the total expected utility. We observe that although the form does not drastically change, the total utility increases with increasing size of the leading group.
Biomimetics 05 00010 g004
Figure 5. Optimal sensor placement for the pressure sensors and the “size of the leading school” experiment. In (a) the utility estimator for the first five sensors and in (b) the value of the utility estimator at the optimal sensor location for the first 20 sensors are presented. In (c), the distribution of the sensors on the swimmer surface is presented. Here, the numbers associated to each sensor indicate that this location is the i-th sensor location chosen according to Equation (21).
Figure 5. Optimal sensor placement for the pressure sensors and the “size of the leading school” experiment. In (a) the utility estimator for the first five sensors and in (b) the value of the utility estimator at the optimal sensor location for the first 20 sensors are presented. In (c), the distribution of the sensors on the swimmer surface is presented. Here, the numbers associated to each sensor indicate that this location is the i-th sensor location chosen according to Equation (21).
Biomimetics 05 00010 g005
Figure 6. Optimal sensor placement for the pressure gradient sensors for the “relative position” experiment. In (a), the utility estimator for the first five sensors and in (b) the value of the utility estimator at the optimal sensor location for the first 20 sensors are presented. In (c), the distribution of the sensors on the swimmer surface is presented. Here, the numbers associated to each sensor indicate that this location is the i-th sensor location chosen according to Equation (21).
Figure 6. Optimal sensor placement for the pressure gradient sensors for the “relative position” experiment. In (a), the utility estimator for the first five sensors and in (b) the value of the utility estimator at the optimal sensor location for the first 20 sensors are presented. In (c), the distribution of the sensors on the swimmer surface is presented. Here, the numbers associated to each sensor indicate that this location is the i-th sensor location chosen according to Equation (21).
Biomimetics 05 00010 g006
Figure 7. (a) Estimated posterior probability for a single sensor optimally placed and a single configuration per group size. The posterior shows clear peaks at the correct number of swimmer for all cases, leading to perfect inference of the parameter of interest. The posterior probability for (b) optimal and (c) worst sensor location for multiple configurations per group size. Here, for the optimal sensor location and one, two, three and five swimmer we see a clear peak for the true size of the group. For the worst sensor location the posterior is almost uniform and does not allow to extract any information about the size of group.
Figure 7. (a) Estimated posterior probability for a single sensor optimally placed and a single configuration per group size. The posterior shows clear peaks at the correct number of swimmer for all cases, leading to perfect inference of the parameter of interest. The posterior probability for (b) optimal and (c) worst sensor location for multiple configurations per group size. Here, for the optimal sensor location and one, two, three and five swimmer we see a clear peak for the true size of the group. For the worst sensor location the posterior is almost uniform and does not allow to extract any information about the size of group.
Biomimetics 05 00010 g007
Figure 8. Estimated posterior for the final location for the best (left column) and worst (right column) sensor-location for one (upper row) and three sensors (lower row). Light colors correspond to high probability density values. We marked the actual location with a black circle.
Figure 8. Estimated posterior for the final location for the best (left column) and worst (right column) sensor-location for one (upper row) and three sensors (lower row). Light colors correspond to high probability density values. We marked the actual location with a black circle.
Biomimetics 05 00010 g008
Figure 9. Optimal sensor locations for the shear stress measurements for the “size of the leading school” in (a) experiment and “relative position” experiment in (b).
Figure 9. Optimal sensor locations for the shear stress measurements for the “size of the leading school” in (a) experiment and “relative position” experiment in (b).
Biomimetics 05 00010 g009

Share and Cite

MDPI and ACS Style

Weber, P.; Arampatzis, G.; Novati, G.; Verma, S.; Papadimitriou, C.; Koumoutsakos, P. Optimal Flow Sensing for Schooling Swimmers. Biomimetics 2020, 5, 10. https://doi.org/10.3390/biomimetics5010010

AMA Style

Weber P, Arampatzis G, Novati G, Verma S, Papadimitriou C, Koumoutsakos P. Optimal Flow Sensing for Schooling Swimmers. Biomimetics. 2020; 5(1):10. https://doi.org/10.3390/biomimetics5010010

Chicago/Turabian Style

Weber, Pascal, Georgios Arampatzis, Guido Novati, Siddhartha Verma, Costas Papadimitriou, and Petros Koumoutsakos. 2020. "Optimal Flow Sensing for Schooling Swimmers" Biomimetics 5, no. 1: 10. https://doi.org/10.3390/biomimetics5010010

Article Metrics

Back to TopTop