Optimal Flow Sensing for Schooling Swimmers

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-widthw(s) along the midline [57,58]

w(s)=2_whs−^s2,0≤s<_sb,_wh−(_wh−_wt)s−_sbst−_sb,_sb≤s<_st,_wtL−sL−_st,_st≤s≤L.

A sketch of the parametrization is presented in Figure 1. Following [59], we use_wh=_sb=0.04L,_st=0.95Land_wt=0.01L. 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)sin2πt_Tp−2πsL.

Here_Tp=1is the tail-beat period and A is the undulation amplitude which linearly increases fromA(0)=0.82/LtoA(L)=5.7/L to replicate the anguilliform swimming motion described by [60]. Given the curvature along s and a center of mass, the coordinatesr(s,t) of the swimmer’s midline can be computed by integrating the Frenet–Serret formulas [59]. In turn, the half-widthw(s)and the coordinatesr(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 velocityuinside the swimmers’ bodies and by including in the NSE a penalization term to enforce no-slip and no-through boundary conditions,

^uk+1−^ukδt=−∇^pk−(^uk·∇)^uk+νΔ^uk+∑i=1_Nsλ_χi(_us,ik−^uk).

Here,νis the kinematic viscosity,λ=1/δtis the penalization coefficient,_Nsis the number of swimmers,_us,ikis the velocity field imposed by swimmer i (composed of translational, rotational and undulatory motions), and_χiis 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],

Δ^pk=1δt∇·^u˜k−1δt∑i=1_Nsχi∇·^u˜k,

where^u˜k=^uk−(^uk·∇)^uk+νΔ^uk . 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_Tp,iof each swimmer i is increased or decreased if it lags behind or surpasses respectively a desired positionΔ_xiin the direction of the school’s motion,

_Tp,i=_Tp(1−Δ_xi).

The mean school trajectory is adjusted by imposing an additional uniform curvature_kC,ialong each swimmer’s midline in order to minimize its lateral deviationΔ_yiand its angular deflectionΔ_θi,

_kC,i=_{[Δyi,⟨Δθi⟩]–}+_{[⟨Δyi⟩,Δθi]–}+_{[⟨Δyi⟩,⟨Δθi⟩]–}.

Here,⟨·⟩defines an exponential moving average with weightδt/_Tp, which approximates the integral term found in PI controllers and

_[a,b]–=|a|b,ifab<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Δ_yiandΔ_θihave 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 velocityu=(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 velocityu 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 S_CIP_Y [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 locationrof the center of mass and (ii) the population_nfof the leading group. We denote withϑ=rorϑ=_nfthese unknown quantities and allow the follower to update its prior beliefp(ϑ)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_x1(_sk)and_x2(_sk), respectively. The corresponding measurements are denoted by_yk1and_yk2 , respectively (see Figure 2 for a sketch of the setup).

We denote byF(ϑ;s)∈^R2nthe 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)∈^R2n×2n. In this case the likelihood of a measurement is given by,

py|ϑ,s=1^(2π)2ndet(Σ(s))exp−12^{y−F(ϑ;s)⊤ Σ−1}(s)y−F(ϑ;s).

The covariance matrix depends on the sensor positionssand 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,

_Σij(s)=^σ2exp−∥_x1(_si)−_x1(_sj)∥ℓ,if1≤i,j≤n,^σ2exp−∥_x2(_si−n)−_x2(_sj−n)∥ℓ,ifn<i,j≤2n,0,otherwise,

whereℓ>0is 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.1L. The correlation strength is set to be two times the average of the signals coming from the simulations,

σ=1n_Nϑ∑j=12n∑i=1_Nϑ|F(^ϑ(i);_sj)|,

where^ϑ(i)are samples from the distributionp(ϑ). 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_si=_sjfori≠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¹⁰⁻⁷when_si=_sj . 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 locationssyielding 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 functionsp(ϑ)andp(ϑ|y,s), respectively. We denote byTthe support ofp(ϑ). The two densities are connected through Bayes’ theorem,

p(ϑ|y,s)=p(y|ϑ,s)p(ϑ|s)p(y|s),

wherep(y|ϑ,s) is the likelihood function defined in Equation (9) andp(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):=_DKL(p(ϑ|y,s)||p(ϑ))=_∫Tlnp(ϑ|y,s)p(ϑ)p(ϑ|y,s)ϑ.

The expected utility is defined as the average value over all possible measurements,

U(s):=_Ey|su(s,y)=_∫Yu(s,y)p(y|s)dy=_{∫Y ∫T}lnp(ϑ|y,s)p(ϑ)p(ϑ|y,s)dϑp(y|s)dy,

whereY is the domain of all possible measurements. Using Equation (12) the expected utility can be expressed as,

U(s)=_{∫Y ∫T}lnp(y|ϑ,s)p(y|s)p(y|ϑ,s)p(ϑ)dϑdy.

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

Whenϑ=ris a continuous random variable andϑ∈Ω⊂^R2. The estimator for the expected utility in this case can be obtained by approximating the two integrals by Monte Carlo integration using_Nϑsamples fromp(ϑ)and_Nysamples fromp(y|ϑ,s) [70]. The resulting estimator is given by,

U(s)≈U^(s)=1_{Nϑ Ny}∑j=1_Ny∑i=1_Nϑlnp(^y(i,j)|^ϑ(i),s)−ln1_Nϑ∑k=1_Nϑp(^y(i,j)|^ϑ(k),s),

where^ϑ(i)∼_pϑ(·)fori=1,…,_Nϑand^y(i,j)∼_py(·|^ϑ(i),s)forj=1,…,_Ny. 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_Ny 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)_∫Ylnp(y|_ϑi,s)p(y|s)p(y|_ϑi,s)dy.

Here,ϑ=_nfrepresents the number of swimmers in the leading group. An estimator of the given utility can be obtained by Monte Carlo integration using_Nysamples from the likelihood distributionp(y|_ϑi,s). The estimator is given by

U(s)≈U^(s)=1_Ny∑j=1_Ny∑i=1_Nϑp(_ϑi)lnp(^y(i,j)|_ϑi,s)−ln∑k=1_Nϑp(_ϑk)p(^y(i,j)|_ϑk,s).

where^y(i,j)∼_py(·|^ϑ(i),s)forj=1,…,_Ny. Letφbe the random variable representing one of the group configurations. Each group configuration is associated with a unique number_φi,ℓforℓ=1,…,_ni, where_niis the total number of configurations containing i swimmers. With this notation,φtakes values in the set{_φi,ℓ|i=1,…,8,ℓ=1,…,_ni} . For examples of different configurations see Appendix A.

Using the fact that fori=1,…,_Nϑ,

p(y,φ=_φk,ℓ|ϑ=_ϑi,s)=0,fork≠i,

and

p(y|ϑ=_ϑi,φ=_φi,ℓ,s)=p(y|φ=_φi,ℓ,s),forℓ=1,…,_ni,

and the assumption

p(φ=_φi,ℓ|ϑ=_ϑi,s)=1_ni,forℓ=1,…,_ni,

the likelihood function can be written as,

p(y|ϑ=_ϑi,s)=∑k=1_Nϑ∑ℓ=1_nip(y,φ=_φk,ℓ|ϑ=_ϑi,s)=∑ℓ=1_nip(y,φ=_φi,ℓ|ϑ=_ϑi,s)=∑ℓ=1_nip(y|ϑ=_ϑi,φ=_φi,ℓ,s)p(φ=_φi,ℓ|ϑ=_ϑi)=1_ni∑ℓ=1_nip(y|φ=_φi,ℓ,s).

Notice that the likelihood function for fixed_ϑi, is a mixture of Gaussian distributions with equal weights and thatp(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_niand then drawy∼_py(·|_φi,^ℓ*,s).

The final form of the estimator is given by

U^(s)=1_Ny∑j=1_Ny∑i=1_Nϑp(_ϑi)ln1_ni∑ℓ=1_nip(^y(i,j)|_φi,ℓ,s)−ln1_ni∑k=1_Nϑp(^ϑ(k))∑ℓ=1_nip(^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 estimatorU^(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_s1☆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 settings=(_s1☆,s)and repeating the grid search for the new optimal location_s2☆. This procedure is then continued by defining

_si☆=argmaxsU^(s)wheres=(_s1☆,…,_si−1☆,s).

We note that the scalar variable s denotes the mid-line coordinate of a single sensor-pair, whereas the vectorsholds 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 numberRe=^L2ν=2000. In all experiments, we use 4096 points to discretize the horizontal directionx∈[0,1]and all artificial swimmers have a length ofL=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_nidifferent initial configurations. In each configuration we assign a number_φi,ℓfori=1,⋯,8andℓ=1,⋯,_ni. In total, we consider_Ntot=_{∑i ni}=61distinct configurations each having the same prior probability1/_Ntot . In Appendix A we present the initial condition for all configurations. The center of mass of the school is located atx=0.3and 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 with21×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 forT=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 timeT+2. The resulting sensor-signal associated to the midline coordinatessfor a given configurationϑis denotedF(ϑ;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 locationss . 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 gradientsF(ϑ;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 ats=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 sety=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×¹⁰⁻⁴. For the “relative position” experiment we find the optimal location^s*=3.84×¹⁰⁻⁴. 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.

Figure 1. Parametrization of the swimmer surface as described in Equation (1).

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 sensorssthat allows the follower to determine with lowest uncertainty the relative positionrand the number of swimmersnfin the leading group of swimmers. For each sensorskthe swimmer collects measurementsyk1andyk2at locationsx1(sk)andx2(sk)on the skin, respectively.

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 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 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 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 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 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 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).

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 extenty∈[0,0.5]was discretized using 2048 gridpoints, for the ones with four rows it was extended toy∈[0,0.75]and discretized using 3072 gridpoints.

Biomimetics 05 00010 g0a1 550

Figure A1.Configurations for two leading swimmers.

Figure A1.Configurations for two leading swimmers.

Biomimetics 05 00010 g0a2 550

Figure A2.Configurations for three leading swimmers.

Figure A2.Configurations for three leading swimmers.

Biomimetics 05 00010 g0a3 550

Figure A3.Configurations for four leading swimmers.

Figure A3.Configurations for four leading swimmers.

Biomimetics 05 00010 g0a4 550

Figure A4.Configurations for five leading swimmers.

Figure A4.Configurations for five leading swimmers.

Biomimetics 05 00010 g0a5 550

Figure A5.Configurations for six leading swimmers.

Figure A5.Configurations for six leading swimmers.

Biomimetics 05 00010 g0a6 550

Figure A6.Configurations for seven leading swimmers.

Figure A6.Configurations for seven leading swimmers.

Biomimetics 05 00010 g0a7 550

Figure A7.Configurations for eight leading swimmers.

Figure A7.Configurations for eight leading swimmers.

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 locations.

In Section 2.2.3 we showed that the distribution of_ϑiconditioned on measurementsy, under the assumption of uniform prior, is proportional to

p(_ϑi|y)∝p(y|_ϑi)p(_ϑi)=181_ni∑ℓ=1_ni12π^σ2exp-^{(y-F(_φi,ℓ))2}2^σ2.

We want to show that for anyi≠j,

p(_ϑi|_ϑj)=p(_ϑi|y=F(_φj,k))≠p(_ϑj|y=F(_φi,ℓ))=p(_ϑj|_ϑi),

for any configurations_φj,kand_φi,ℓcorresponding to a school of size_ϑjand_ϑi , respectively. From Equation (A1) it is easy to see that (A2) is true due to the fact that

1_ni∑^ℓ'=1_niexp-^{(F(_φj,k)-F(_φi,ℓ'))2}2^σ2≠1_nj∑^k'=1_njexp-^{(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.,_ni=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.

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.

3. Triantafyllou, M.S.; Weymouth, G.D.; Miao, J. Biomimetic Survival Hydrodynamics and Flow Sensing. Annu. Rev. Fluid Mech. 2016, 48, 1-24.

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.

5. Puckett, J.G.; Pokhrel, A.R.; Giannini, J.A. Collective gradient sensing in fish schools. Sci. Rep. 2018, 8, 7587.

6. Dykgraaf, S. Untersuchungen über die Funktion der Seitenorgane an Fischen. Zeitschrift für Vergleichende Physiologie 1933, 20, 162-214.

7. Dykgraaf, S. The functioning and significance of the lateral-line organs. Biol. Rev. Camb. Philos. Soc. 1963, 38, 51-105.

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.

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.

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.

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.

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.

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.

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.

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.

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.

17. Pitcher, T.; Partridge, B.; Wardle, C. A blind fish can school. Science 1976, 194, 963-965.

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.

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.

20. Montgomery, J.C.; Baker, C.F.; Carton, A.G. The lateral line can mediate rheotaxis in fish. Nature 1997, 389, 960-963.

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.

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.

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.

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.

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.

26. Bleckmann, H. Peripheral and central processing of lateral line information. J. Comp. Physiol. A 2008, 194, 145-158.

27. Jiang, Y.; Ma, Z.; Zhang, D. Flow field perception based on the fish lateral line system. Bioinspir. Biomim. 2019, 14, 041001.

28. Engelmann, J.; Hanke, W.; Mogdans, J.; Bleckmann, H. Hydrodynamic stimuli and the fish lateral line. Nature 2000, 408, 1476-4687.

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.

30. Tao, J.; Yu, X. Hair flow sensors: From bio-inspiration to bio-mimicking-A review. Smart Mater. Struct. 2012, 21, 113001.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

52. Colvert, B.; Alsalman, M.; Kanso, E. Classifying vortex wakes using neural networks. Bioinspir. Biomim. 2018, 13, 025003.

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.

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.

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.

56. Verma, S.; Papadimitriou, C.; Lüthen, N.; Arampatzis, G.; Koumoutsakos, P. Optimal sensor placement for artificial swimmers. J. Fluid Mech. 2020, 884, A24.

57. Kern, S.; Koumoutsakos, P. Simulations of optimized anguilliform swimming. J. Exp. Biol. 2006, 209, 4841-4857.

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.

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.

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.

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.

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.

63. Towers, J.D. Finite difference methods for approximating Heaviside functions. J. Comput. Phys. 2009, 228, 3478-3489.

64. Chorin, A.J. Numerical Solution of the Navier-Stokes Equations. Math. Comput. 1968, 22, 745-762.

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.

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.

67. Bleckmann, H.; Zelick, R. Lateral line system of fish. Integr. Zool. 2009, 4, 13-25.

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.

69. Dierckx, P. An algorithm for smoothing, differentiation and integration of experimental data using spline functions. J. Comput. Appl. Math. 1975, 1, 165-184.

70. Huan, X.; Marzouk, Y.M. Simulation-based optimal Bayesian experimental design for nonlinear systems. J. Comput. Phys. 2013, 232, 288-317.

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.

72. Simoen, E.; Papadimitriou, C.; Lombaert, G. On prediction error correlation in Bayesian model updating. J. Sound Vib. 2013, 332, 4136-4152.

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.

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.

75. Papadimitriou, C. Optimal sensor placement methodology for parametric identification of structural systems. J. Sound Vib. 2004, 278, 923-947.