Hydrodynamic model of fish orientation in a channel flow

  1. Maurizio Porfiri  Is a corresponding author
  2. Peng Zhang
  3. Sean D Peterson  Is a corresponding author
  1. Department of Biomedical Engineering, New York University, United States
  2. Department of Mechanical and Aerospace Engineering, New York University, United States
  3. Center for Urban Science and Progress, New York University Tandon School of Engineering, United States
  4. Mechanical and Mechatronics Engineering Department, University of Waterloo, Canada
8 figures, 4 tables and 1 additional file

Figures

Fish rheotaxis.

(a) Illustration of the problem with notation, showing a fish swimming in a background flow described by Equation 4. (b) Schematic of the cross-stream sweeping movement of some fish species swimming without visual cues; snapshots of fish at earlier time instants are illustrated by lighter shading.

Method of images.

Schematic of the fish (black) in the channel (thick lines) and the set of images (gray) needed to generate the channel. The streamlines generated by the fish in an otherwise quiescent fluid are shown in the channel colored by local velocity magnitude (red: high; blue: low). Dashed and solid lines are mirroring planes for the method of images, the pattern for which continues ad infinitum.

Velocity field around a swimming fish from computational fluid dynamics.

(a) Mean velocity field around the steady swimming giant danio relative to the background flow. (b) Velocity field predicted by a dipole with θ=π located at 0.315l from the fish head along its centerline relative to the background flow. The selection of the dipole location and strength is detailed in Appendix 2.

Qualitative dynamics of equation set Equation 10a.

(a) Cross-stream equilibria for upstream swimming as a function of β. (b,c) Phase plot for downstream and upstream swimming in the case α=0.1, ρ=0.1, and κ=1, so that β=20. In all panels, red refers to unstable equilibria and green to stable equilibria.

Visual illustration of the process of determining the roots of Equation 14.

(a) Plot of the function π332cot(π(ξ+12))csc2(π(ξ+12)) (black), superimposed with three lines of different slope: 200 (red), -200 (dashed blue), and -2 (solid blue). (b) Zoomed-in view of the curves in (a) showing that the blue line can only intersect the black curve at the origin.

Appendix 2—figure 1
Details about the implementation of the computational fluid dynamics simulations.

(a) Mesh implemented in the simulations, with definitions of coordinate systems and a zoomed-in view of the refined mesh around the fish. (b) Mesh convergence analysis, showing the mean drag as a function of number of elements in the simulation.

Appendix 2—figure 2
Instantaneous velocity fields around a swimming fish relative to the background flow from computational fluid dynamics.

(a) – (e) correspond to t=0, T/8, T/4, 3T/8, and T/2, respectively, where T is the period of a tail beat. White curves are streamlines with arrows indicating flow directions.

Appendix 2—figure 3
Analysis of the boundary layer thickness along a swimming fish.

x-component of the flow velocity, ut, extracted across the half-length of the fish body. The values of ut are measured in a coordinate system moving at the speed of v0 along the fish swimming direction.

Tables

Table 1
Estimation of model parameters from data in the literature.
Referenceρϵακβ
Bak-Coleman et al., 20130.05[10-2,10-1][0,0.17]
Bak-Coleman and Coombs, 20140.04[10-2,10-1][0,0.16]0[0,100]
Baker and Montgomery, 1999 and Montgomery et al., 19970.1[10-2,10-1][0,0.32][2,7][0,256]
Elder and Coombs, 20150.066[10-2,10-1][0,0.24]0[0,55]
Kulpa et al., 20150.041 near center of jet1.3
Oteiza et al., 2017[0.018,0.066][0.20,0.82]
Peimani et al., 20170.0441
Suli et al., 20120.018[0.1,1]
Van Trump and McHenry, 2013[0.055,0.127][10-2,10-1][0,0.32]0[0,106]
  1. LL+ cavefish swimming speed v05cm/s in zero background flow in Bak-Coleman and Coombs, 2014 is used to estimate α.

Appendix 2—table 1
Parameters employed in (27) to describe the locomotory pattern of a giant danio.
Parameterslc1c2kLf
Values7.3cm0.004-2.33m-178.5m-13Hz
Appendix 3—table 1
Relevant publications on fish rheotaxis in the absence of visual cues, identified through literature review.
ReferenceFishSwimming domainFlow properties*Sensory cuesRheotaxis threshold speed
SpeciesLengthFlow speedFlow gradient
Bak-Coleman et al., 2013Giant danio (Devario aequipinnatus)6.0 –7.3 cmFlow tank of 25×25×25cm (L×h×W)0, 3, and 7cm/sRe7500 at LL+ threshold speed; flow gradient expected to be small near center of tankLL+/LL-3cm/s
Bak-Coleman and Coombs, 2014blind cavefish (Astyanax mexicanus)4.2 –5.0 cmFlow tank of 25×25×25cm (L×h×W)0, 1, 2, 3, 4, 7 and8cm/sRe2000 at LL+ threshold speed; flow gradient expected to be small near center of tankLL+/LL-; fish made transient contacts with substrateLL+: 0.90cm/s; LL-:0.54cm/s
Baker and Montgomery, 1999blind cavefish (Astyanax fasciatus)4 –7 cmFlow tank of 51×9×20cm (L×h×W)0, 2, 3, 5, 9 and 16cm/sRe2000 at LL+ threshold speed; flow gradient expected to be small near center of tankLL+/LL-; tactile sensesLL+: 2–3cm/s; LL-: 9–16cm/s
Elder and Coombs, 2015Mexican tetras (Astyanax mexicanus)8.3cmFlow tank of 25×25×25cm (L×h×W)0, 1, 2, 4, 7, and 12cm/sRe5000 at threshold speed; flow gradient expected to be small near center of tankLL+/LL-2cm/s for LL+ and LL-
Kulpa et al., 2015blind cavefish (Astyanax mexicanus)4.4 –5.3 cmFlow tank of 25×25×10cm (L×h×W)Maximum speed of 8cm/sJet flow across center of tank; flow gradient expected to be largeLL+/LL-8cm/s
Lyon, 1904blind FundulusunspecifiedTrough with unspecified dimensions; tideway leading to pond“not too strong current” in trough and current with “more or less eddy and irregularity” in tidewayFlow gradient expected to be smallLL+; some fish gained tactile sensesNot measured; rheotaxis elicited only by tactile cues
Lyon, 1904blind FundulusunspecifiedTrough with unspecified dimensionsflow “gushing rather violently”Jet flow; flow gradient expected to be largeLL+Not measured; rheotaxis elicited by flow
Montgomery et al., 1997blind cavefish (Astyanax fasciatus)4 –7 cm§Flow tank of 51×9×20cm (L×h×W)0, 2, 3, 5, 9, and 16cm/sRe2000 at LL+ threshold speed; flow gradient expected to be small near center of tankLL+/LL-; tactile sensesLL+: 2–3 cm/s; LL-: 9–16 cm/s
Oteiza et al., 2017zebrafish (Danio rerio) larva 5–7 days post fertilization (dpf)unspecified13 cm-long circular tube with diameter 1.27– 4.76cm0.2–0.8 cm/sLow to high flow gradients identified through particle image velocimetryLL+/LL-LL+: rheotaxis observed as low as 0.2cm/s
Peimani et al., 2017zebrafish (Danio rerio) larva 5–7 dpfestimated 
0.35cm
Flow channel of 63.3×1.6×0.55mm (L×h×W)0.95–3.8 cm/sRe10 at threshold speed; flow gradient expected to be largeLL+0.95cm/s
Suli et al., 2012zebrafish (Danio rerio) larva 5 dpf0.33cmFlume of 110×3.7×2.8cm (L×h×W)0.075, 0.15, 0.2cm/sRe<75; flow gradient expected to be largeLL+/LL-Not quantified
Van Trump and McHenry, 2013blind Mexican cavefish (Astyanax fasciatus)3 –7 cmCylindrical channel of 150×11cm (L×D)0, 1, 2, 4, 6, 8, 10, 13,16cm/sRe>2000 at threshold speed; flow gradient expected to be small near center of tankLL+/LL-2–4 cm/s
  1. *

    LL+: lateral line enabled; LL−: lateral line disabled

  2. Data are extracted from the same set of experiments

  3. Two experiments are considered from the same paper

  4. §
Appendix 3—table 2
Results of the bibliographical research on fish rheotaxis in the absence of visual cues, used to validate the proposed model.
ReferenceFish species*EvidenceComparison with model
SupportiveInconclusive
Within studies
Effect of lateral line
Bak-Coleman et al., 2013Giant danio (Devario aequipinnatus)No significant difference in fish heading angle against current was detected between LL+ and LL-×
Bak-Coleman and Coombs, 2014blind cavefish (Astyanax mexicanus)Rheotaxis threshold speed was slightly (but not significantly) lower in LL- condition×
Baker and Montgomery, 1999 and Montgomery et al., 1997blind cavefish (Astyanax fasciatus)Rheotaxis threshold speed was significantly higher in LL- condition; fish received intermittent tactile senses×
Elder and Coombs, 2015Mexican tetras (Astyanax mexicanus)No significant influence of LL condition was detected on rheotactic performance×
Kulpa et al., 2015blind cavefish (Astyanax mexicanus)Significantly higher rheotaxis index in LL+ fish than LL- fish in jet stream×
Oteiza et al., 2017zebrafish (Danio rerio) larva 5–7 dpfPosterior lateral line ablation or chemical neuromast ablation severely reduced rheotaxis×
Suli et al., 2012zebrafish (Danio rerio) larva 5 dpfLL hair cell damage led to a significant decrease in rheotaxis; regeneration of LL hair cells restored rheotaxis×
Van Trump and McHenry, 2013blind Mexican cavefish (Astyanax fasciatus)In LL+ and LL-, fish exhibited statistically indistinguishable rheotaxis behavior×
Effect of flow gradient
Lyon, 1904blind FundulusIn a flow with small gradient, rheotaxis was elicited only when fish received tactile cues; in jet flow with large gradient, rheotaxis was elicited by flow without tactile cues. Lack of data on statistical significance×
Oteiza et al., 2017zebrafish (Danio rerio) larva 5–7 dpfRheotaxis of fish improved with increasing gradient magnitudes×
Across studies
Effect of channel width
Bak-Coleman and Coombs, 2014; Baker and Montgomery, 1999blind cavefish (Astyanax mexicanus); blind cavefish (Astyanax fasciatus)Significantly different threshold speed for LL+ fish: 0.90±0.137cm/s (mean ±s.e.m.) in 25cm wide tunnel; between 2cm/s and 3cm/s in 9cm wide tunnel. Tactile cues available to fish in Bak-Coleman and Coombs, 2014×
Bak-Coleman and Coombs, 2014; Van Trump and McHenry, 2013blind cavefish (Astyanax mexicanus); blind cavefish (Astyanax fasciatus)Significantly different threshold speed for LL+ fish: 0.90±0.137cm/s (mean ±s.e.m.) in 25cm wide tunnel; between 2cm/s and 4cm/s in 11cm diameter tunnel. Tactile cues available to fish in Bak-Coleman and Coombs, 2014×
Oteiza et al., 2017; Peimani et al., 2017zebrafish (Danio rerio) larva 5–7 dpfOnset of rheotaxis in LL+ fish observed at flow speed 0.95cm/s in 1.6mm wide tunnel; rheotaxis observed in LL+ fish at flow speed 0.2cm/s in 2.22cm diameter tunnel×
Effect of body length
Bak-Coleman and Coombs, 2014; Elder and Coombs, 2015blind cavefish (Astyanax mexicanus); Mexican tetras (Astyanax mexicanus)Significantly different threshold speed for LL+ fish: 0.90±0.137cm/s (mean ±s.e.m.) for 4.2–5.0cm long fish; 1.96±0.350cm/s (mean ±s.e.m.) for 8.3cm long fish. Tactile cues available to fish in Bak-Coleman and Coombs, 2014×
Total59
  1. *

    LL+: lateral line enabled; LL−: lateral line disabled

Additional files

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Maurizio Porfiri
  2. Peng Zhang
  3. Sean D Peterson
(2022)
Hydrodynamic model of fish orientation in a channel flow
eLife 11:e75225.
https://doi.org/10.7554/eLife.75225