skip to main content


Title: Hydrodynamics of Metachronal Rowing at Intermediate Reynolds Numbers
Abstract

Metachronous rowing is a swimming mechanism widely adopted by small marine invertebrate like comb jellies, in which rows of appendages perform propulsive strokes sequentially in a coordinated manner with a fixed phase difference. To simulate metachronous rowing at intermediate Reynolds number, in this paper, a row of flexible cilia models was placed inside the flow field, with their roots stroke at a sinusoidal function of time and a fixed phase difference. A fully coupled two-way numerical solver was developed, which solves the Navier-Stokes equations for the fluid field coupled with the differential equation for the flexible cilia model. This numerical solver is applied to investigate how the row of cilia models are deformed by the hydrodynamic forces (pressure and shear) and momentum and thus impact hydrodynamic performance. Results show that the passive deformation of cilia potentially improve the hydrodynamic performance compared to the rigid cilia. With the metachronous rowing mechanism, the cilia generate the thrust to move forward. The approach used in this study presents a general way to explore the fluid dynamics of complex fluid-structure interaction problems.

 
more » « less
Award ID(s):
2120505
PAR ID:
10499346
Author(s) / Creator(s):
; ; ; ;
Publisher / Repository:
American Society of Mechanical Engineers
Date Published:
Journal Name:
ASME 2023 International Mechanical Engineering Congress and Exposition
ISBN:
978-0-7918-8766-0
Format(s):
Medium: X
Location:
New Orleans, Louisiana, USA
Sponsoring Org:
National Science Foundation
More Like this
  1. Metachronal rowing is a biological propulsion mechanism employed by many swimming invertebrates (e.g. copepods, ctenophores, krill and shrimp). Animals that swim using this mechanism feature rows of appendages that oscillate in a coordinated wave. In this study, we used observations of a swimming ctenophore (comb jelly) to examine the hydrodynamic performance and vortex dynamics associated with metachronal rowing. We first reconstructed the beating kinematics of ctenophore appendages based on a high-speed video of a metachronally coordinated row. Following the reconstruction, two numerical models were developed and simulated using an in-house immersed-boundary-method-based computational fluid dynamics solver. The two models included the original geometry (16 appendages in a row) and a sparse geometry (8 appendages, formed by removing every other appendage along the row). We found that appendage tip vortex interactions contribute to hydrodynamic performance via a vortex-weakening mechanism. Through this mechanism, appendage tip vortices are significantly weakened during the drag-producing recovery stroke. As a result, the swimming ctenophore produces less overall drag, and its thrust-to-power ratio is significantly improved (up to 55.0 % compared with the sparse model). Our parametric study indicated that such a propulsion enhancement mechanism is less effective at higher Reynolds numbers. Simulations were also used to investigate the effects of substrate curvature on the unsteady hydrodynamics. Our results illustrated that, compared with a flat substrate, arranging appendages on a curved substrate can boost the overall thrust generation by up to 29.5 %. These findings provide new insights into the fluid dynamic principles of propulsion enhancement underlying metachronal rowing.

     
    more » « less
  2. Discher, Dennis (Ed.)
    Hydrodynamic flow produced by multiciliated cells is critical for fluid circulation and cell motility. Hundreds of cilia beat with metachronal synchrony for fluid flow. Cilia-driven fluid flow produces extracellular hydrodynamic forces that cause neighboring cilia to beat in a synchronized manner. However, hydrodynamic coupling between neighboring cilia is not the sole mechanism that drives cilia synchrony. Cilia are nucleated by basal bodies (BBs) that link to each other and to the cell’s cortex via BB-associated appendages. The intracellular BB and cortical network is hypothesized to synchronize ciliary beating by transmitting cilia coordination cues. The extent of intracellular ciliary connections and the nature of these stimuli remain unclear. Moreover, how BB connections influence the dynamics of individual cilia has not been established. We show by focused ion beam scanning electron microscopy imaging that cilia are coupled both longitudinally and laterally in the ciliate Tetrahymena thermophila by the underlying BB and cortical cytoskeletal network. To visualize the behavior of individual cilia in live, immobilized Tetrahymena cells, we developed Delivered Iron Particle Ubiety Live Light (DIPULL) microscopy. Quantitative and computer analyses of ciliary dynamics reveal that BB connections control ciliary waveform and coordinate ciliary beating. Loss of BB connections reduces cilia-dependent fluid flow forces. 
    more » « less
  3. Summary

    In this paper, a three‐dimensional numerical solver is developed for suspensions of rigid and soft particles and droplets in viscoelastic and elastoviscoplastic (EVP) fluids. The presented algorithm is designed to allow for the first time three‐dimensional simulations of inertial and turbulent EVP fluids with a large number particles and droplets. This is achieved by combining fast and highly scalable methods such as an FFT‐based pressure solver, with the evolution equation for non‐Newtonian (including EVP) stresses. In this flexible computational framework, the fluid can be modeled by either Oldroyd‐B, neo‐Hookean, FENE‐P, or Saramito EVP models, and the additional equations for the non‐Newtonian stresses are fully coupled with the flow. The rigid particles are discretized on a moving Lagrangian grid, whereas the flow equations are solved on a fixed Eulerian grid. The solid particles are represented by an immersed boundary method with a computationally efficient direct forcing method, allowing simulations of a large numbers of particles. The immersed boundary force is computed at the particle surface and then included in the momentum equations as a body force. The droplets and soft particles on the other hand are simulated in a fully Eulerian framework, the former with a level‐set method to capture the moving interface and the latter with an indicator function. The solver is first validated for various benchmark single‐phase and two‐phase EVP flow problems through comparison with data from the literature. Finally, we present new results on the dynamics of a buoyancy‐driven drop in an EVP fluid.

     
    more » « less
  4. ABSTRACT

    Tidal interactions in coalescing binary neutron stars modify the dynamics of the inspiral and hence imprint a signature on their gravitational wave (GW) signals in the form of an extra phase shift. We need accurate models for the tidal phase shift in order to constrain the supranuclear equation of state from observations. In previous studies, GW waveform models were typically constructed by treating the tide as a linear response to a perturbing tidal field. In this work, we incorporate non-linear corrections due to hydrodynamic three- and four-mode interactions and show how they can improve the accuracy and explanatory power of waveform models. We set up and numerically solve the coupled differential equations for the orbit and the modes and analytically derive solutions of the system’s equilibrium configuration. Our analytical solutions agree well with the numerical ones up to the merger and involve only algebraic relations, allowing for fast phase shift and waveform evaluations for different equations of state over a large parameter space. We find that, at Newtonian order, non-linear fluid effects can enhance the tidal phase shift by $\gtrsim 1\, {\rm radian}$ at a GW frequency of 1000 Hz, corresponding to a $10{{\%}}-20{{\%}}$ correction to the linear theory. The scale of the additional phase shift near the merger is consistent with the difference between numerical relativity and theoretical predictions that account only for the linear tide. Non-linear fluid effects are thus important when interpreting the results of numerical relativity and in the construction of waveform models for current and future GW detectors.

     
    more » « less
  5. In this paper, we introduce a new nonlocal modal hydrodynamic theory for fluid–structure interactions (FSI) of light, flexible cantilever beams and plates undergoing small amplitude vibrations in Newtonian, incompressible, viscous, heavy fluids otherwise at rest. For low aspect ratio flexible structures and high mode numbers, three dimensional (3D) and nonlocal fluid effects become prominent drivers of the coupled dynamics, to the point that existing local hydrodynamic theories based on two dimensional (2D) fluid approximations become inadequate to predict the system response. On the other hand, our approach is based on a rigorous, yet efficient, 3D treatment of the hydrodynamic loading on cantilevered thin structures. The off-line solution of the FSI problem results in the so-called nonlocal modal hydrodynamic function matrix, that is, the representation of the nonlocal hydrodynamic load operator on a basis formed by the structural modes. Our theory then integrates the nonlocal hydrodynamics within a fully coupled structural modal model in the frequency domain. We compare and discuss our theory predictions in terms of frequency response functions, mode shapes, hydrodynamic loads, quality factors, added mass ratios with the predictions of the classical local approaches, for different actuation scenarios, identifying the limitations of the hypotheses underlying existing treatments. Importantly, we also validate our new model with experiments conducted on flexible square plates. While computationally efficient, our fully coupled theory is exact up to numerical truncation and can bridge knowledge gaps in the design and analysis of FSI systems based on low aspect ratio flexible beams and plates. 
    more » « less