Amneet Pal Singh Bhalla

An immersed boundary method for rigid bodies

We develop an immersed boundary (IB) method for modeling flows around fixed or moving rigid bodies that is suitable for a broad range of Reynolds numbers, including steady Stokes flow. The spatio-temporal discretization of the fluid equations is based on a standard staggered-grid approach. Fluid-body interaction is handled using Peskins IB method; however, unlike existing IB approaches to such problems, we do not rely on penalty or fractional-step formulations. Instead, we use an unsplit scheme that ensures the no-slip constraint is enforced exactly in terms of the Lagrangian velocity field evaluated at the IB markers. Fractionalstep approaches, by contrast, can impose such constraints only approximately, which can lead to penetration of the flow into the body, and are inconsistent for steady Stokes flow. Imposing no-slip constraints exactly requires the solution of a large linear system that includes the fluid velocity and pressure as well as Lagrange multiplier forces that impose the motion of the body. The principal contribution of this paper is that it develops an efficient preconditioner for this exactly constrained IB formulation which is based on an analytical approximation to the Schur complement. This approach is enabled by the near translational and rotational invariance of Peskins IB method. We demonstrate that only a few cycles of a geometric multigrid method for the fluid equations are required in each application of the preconditioner, and we demonstrate robust convergence of the overall Krylov solver despite the approximations made in the preconditioner. We empirically observe that to control the condition number of the coupled linear system while also keeping the rigid structure impermeable to fluid, we need to place the immersed boundary markers at a distance of about two grid spacings, which is significantly larger from what has been recommended in the literature for elastic bodies. We demonstrate the advantage of our monolithic solver over split solvers by computing the steady state flow through a two-dimensional nozzle at several Reynolds numbers. We apply the method to a number of benchmark problems at zero and finite Reynolds numbers, and we demonstrate first-order convergence of the method to several analytical solutions and benchmark computations.

Undulating fins produce off-axis thrust and flow structures

While wake structures of many forms of swimming and flying are well characterized, the wake generated by a freely swimming undulating fin has not yet been analyzed. These elongated fins allow fish to achieve enhanced agility exemplified by the forward, backward and vertical swimming capabilities of knifefish, and also have potential applications in the design of more maneuverable underwater vehicles. We present the flow structure of an undulating robotic fin model using particle image velocimetry to measure fluid velocity fields in the wake. We supplement the experimental robotic work with high-fidelity computational fluid dynamics, simulating the hydrodynamics of both a virtual fish, whose fin kinematics and fin plus body morphology are measured from a freely swimming knifefish, and a virtual rendering of our robot. Our results indicate that a series of linked vortex tubes is shed off the long edge of the fin as the undulatory wave travels lengthwise along the fin. A jet at an oblique angle to the fin is associated with the successive vortex tubes, propelling the fish forward. The vortex structure bears similarity to the linked vortex ring structure trailing the oscillating caudal fin of a carangiform swimmer, though the vortex rings are distorted because of the undulatory kinematics of the elongated fin.

A Forced Damped Oscillation Framework for Undulatory Swimming Provides New Insights into How Propulsion Arises in Active and Passive Swimming

A fundamental issue in locomotion is to understand how muscle forcing produces apparently complex deformation kinematics leading to movement of animals like undulatory swimmers. The question of whether complicated muscle forcing is required to create the observed deformation kinematics is central to the understanding of how animals control movement. In this work, a forced damped oscillation framework is applied to a chain-link model for undulatory swimming to understand how forcing leads to deformation and movement. A unified understanding of swimming, caused by muscle contractions (“active” swimming) or by forces imparted by the surrounding fluid (“passive” swimming), is obtained. We show that the forcing triggers the first few deformation modes of the body, which in turn cause the translational motion. We show that relatively simple forcing patterns can trigger seemingly complex deformation kinematics that lead to movement. For given muscle activation, the forcing frequency relative to the natural frequency of the damped oscillator is important for the emergent deformation characteristics of the body. The proposed approach also leads to a qualitative understanding of optimal deformation kinematics for fast swimming. These results, based on a chain-link model of swimming, are confirmed by fully resolved computational fluid dynamics (CFD) simulations. Prior results from the literature on the optimal value of stiffness for maximum speed are explained.

Energy efficiency and allometry of movement of swimming and flying animals

Significance Is a whale or a tuna more efficient? The answer depends on the definition of efficiency. If one uses cost of transport, defined as the energy required to move unit distance, then the tuna wins. If energy of transporting unit mass over unit distance is used then the whale wins. We address this ambiguity by rationally deriving a new efficiency measure called the energy-consumption coefficient (CE), which is a nondimensional measure of fuel consumption. CE is a fundamental metric to quantify efficiency of self-propelled bodies analogous to the drag coefficient (Cd) to quantify aerodynamic shapes of vehicles. The analysis also leads to allometric scalings of frequency and velocity of swimming and flying organisms over more than 20 orders of magnitude of mass. Which animals use their energy better during movement? One metric to answer this question is the energy cost per unit distance per unit weight. Prior data show that this metric decreases with mass, which is considered to imply that massive animals are more efficient. Although useful, this metric also implies that two dynamically equivalent animals of different sizes will not be considered equally efficient. We resolve this longstanding issue by first determining the scaling of energy cost per unit distance traveled. The scale is found to be M2/3 or M1/2, where M is the animal mass. Second, we introduce an energy-consumption coefficient (CE) defined as energy per unit distance traveled divided by this scale. CE is a measure of efficiency of swimming and flying, analogous to how drag coefficient quantifies aerodynamic drag on vehicles. Derivation of the energy-cost scale reveals that the assumption that undulatory swimmers spend energy to overcome drag in the direction of swimming is inappropriate. We derive allometric scalings that capture trends in data of swimming and flying animals over 10–20 orders of magnitude by mass. The energy-consumption coefficient reveals that swimmers beyond a critical mass, and most fliers are almost equally efficient as if they are dynamically equivalent; increasingly massive animals are not more efficient according to the proposed metric. Distinct allometric scalings are discovered for large and small swimmers. Flying animals are found to require relatively more energy compared with swimmers.

Convergent evolution of mechanically optimal locomotion in aquatic invertebrates and vertebrates.

Examples of animals evolving similar traits despite the absence of that trait in the last common ancestor, such as the wing and camera-type lens eye in vertebrates and invertebrates, are called cases of convergent evolution. Instances of convergent evolution of locomotory patterns that quantitatively agree with the mechanically optimal solution are very rare. Here, we show that, with respect to a very diverse group of aquatic animals, a mechanically optimal method of swimming with elongated fins has evolved independently at least eight times in both vertebrate and invertebrate swimmers across three different phyla. Specifically, if we take the length of an undulation along an animal’s fin during swimming and divide it by the mean amplitude of undulations along the fin length, the result is consistently around twenty. We call this value the optimal specific wavelength (OSW). We show that the OSW maximizes the force generated by the body, which also maximizes swimming speed. We hypothesize a mechanical basis for this optimality and suggest reasons for its repeated emergence through evolution.

Hydrodynamics of suspensions of passive and active rigid particles: a rigid multiblob approach

Author(s): Usabiaga, FB; Kallemov, B; Delmotte, B; Bhalla, APS; Griffith, BE; Donev, A | Abstract:

A fully resolved active musculo-mechanical model for esophageal transport

Esophageal transport is a physiological process that mechanically transports an ingested food bolus from the pharynx to the stomach via the esophagus, a multilayered muscular tube. This process involves interactions between the bolus, the esophagus, and the neurally coordinated activation of the esophageal muscles. In this work, we use an immersed boundary (IB) approach to simulate peristaltic transport in the esophagus. The bolus is treated as a viscous fluid that is actively transported by the muscular esophagus, and the esophagus is modeled as an actively contracting, fiber-reinforced tube. Before considering the full model of the esophagus, however, we first consider a standard benchmark problem of flow past a cylinder. Next a simplified version of our model is verified by comparison to an analytic solution to the tube dilation problem. Finally, three different complex models of the multi-layered esophagus, which differ in their activation patterns and the layouts of the mucosal layers, are extensively tested. To our knowledge, these simulations are the first of their kind to incorporate the bolus, the multi-layered esophagus tube, and muscle activation into an integrated model. Consistent with experimental observations, our simulations capture the pressure peak generated by the muscle activation pulse that travels along the bolus tail. These fully resolved simulations provide new insights into roles of the mucosal layers during bolus transport. In addition, the information on pressure and the kinematics of the esophageal wall resulting from the coordination of muscle activation is provided, which may help relate clinical data from manometry and ultrasound images to the underlying esophageal motor function.

Separability of drag and thrust in undulatory animals and machines

For nearly a century, researchers have tried to understand the swimming of aquatic animals in terms of a balance between the forward thrust from swimming movements and drag on the body. Prior approaches have failed to provide a separation of these two forces for undulatory swimmers such as lamprey and eels, where most parts of the body are simultaneously generating drag and thrust. We nonetheless show that this separation is possible, and delineate its fundamental basis in undulatory swimmers. Our approach unifies a vast diversity of undulatory aquatic animals (anguilliform, sub-carangiform, gymnotiform, bal-istiform, rajiform) and provides design principles for highly agile bioinspired underwater vehicles. This approach has practical utility within biology as well as engineering. It is a predictive tool for use in understanding the role of the mechanics of movement in the evolutionary emergence of morphological features relating to locomotion. For example, we demonstrate that the drag-thrust separation framework helps to predict the observed height of the ribbon fin of electric knifefish, a diverse group of neotropical fish which are an important model system in sensory neurobiology. We also show how drag-thrust separation leads to models that can predict the swimming velocity of an organism or a robotic vehicle.

Gray's paradox: A fluid mechanical perspective

Nearly eighty years ago, Gray reported that the drag power experienced by a dolphin was larger than the estimated muscle power – this is termed as Grays paradox. We provide a fluid mechanical perspective of this paradox. The viewpoint that swimmers necessarily spend muscle energy to overcome drag in the direction of swimming needs revision. For example, in undulatory swimming most of the muscle energy is directly expended to generate lateral undulations of the body, and the drag power is balanced not by the muscle power but by the thrust power. Depending on drag model utilized, the drag power may be greater than muscle power without being paradoxical.

A moving control volume approach to computing hydrodynamic forces and torques on immersed bodies

Abstract We present a moving control volume (CV) approach to computing hydrodynamic forces and torques on complex geometries. The method requires surface and volumetric integrals over a simple and regular Cartesian box that moves with an arbitrary velocity to enclose the body at all times. The moving box is aligned with Cartesian grid faces, which makes the integral evaluation straightforward in an immersed boundary (IB) framework. Discontinuous and noisy derivatives of velocity and pressure at the fluid–structure interface are avoided and far-field (smooth) velocity and pressure information is used. We re-visit the approach to compute hydrodynamic forces and torques through force/torque balance equations in a Lagrangian frame that some of us took in a prior work (Bhalla et al., 2013 [13] ). We prove the equivalence of the two approaches for IB methods, thanks to the use of Peskins delta functions. Both approaches are able to suppress spurious force oscillations and are in excellent agreement, as expected theoretically. Test cases ranging from Stokes to high Reynolds number regimes are considered. We discuss regridding issues for the moving CV method in an adaptive mesh refinement (AMR) context. The proposed moving CV method is not limited to a specific IB method and can also be used, for example, with embedded boundary methods.