Interactions between internal forces, body stiffness, and fluid environment in a neuromechanical model of lamprey swimming

Form and function of anguilliform swimming

Anguilliform swimmers are long and narrow animals that propel themselves by undulating their bodies. Observations in nature and recent investigations suggest that anguilliform swimming is highly efficient. However, understanding the underlying reasons for the efficiency of this type of locomotion requires interdisciplinary studies spanning from biology to hydrodynamics. Regrettably, these different fields are rarely discussed together, which hinders our ability to understand the repeated evolution of this swimming mode in vertebrates. This review compiles the current knowledge of the anatomical features that drive anguilliform swimming, compares the resulting kinematics across a wide range of anguilliform swimmers, and describes the resulting hydrodynamic interactions using data from both in vivo experiments and computational studies.

Robot motor learning shows emergence of frequency-modulated, robust swimming with an invariant Strouhal number

Fish locomotion emerges from diverse interactions among deformable structures, surrounding fluids and neuromuscular activations, i.e. fluid-structure interactions (FSI) controlled by fish's motor systems. Previous studies suggested that such motor-controlled FSI may possess embodied traits. However, their implications in motor learning, neuromuscular control, gait generation, and swimming performance remain to be uncovered. Using robot models, we studied the embodied traits in fish-inspired swimming. We developed modular robots with various designs and used central pattern generators (CPGs) to control the torque acting on robot body. We used reinforcement learning to learn CPG parameters for maximizing the swimming speed. The results showed that motor frequency converged faster than other parameters, and the emergent swimming gaits were robust against disruptions applied to motor control. For all robots and frequencies tested, swimming speed was proportional to the mean undulation velocity of body and caudal-fin combined, yielding an invariant, undulation-based Strouhal number. The Strouhal number also revealed two fundamental classes of undulatory swimming in both biological and robotic fishes. The robot actuators were also demonstrated to function as motors, virtual springs and virtual masses. These results provide novel insights in understanding fish-inspired locomotion.

Complex Modal Characteristic Analysis of a Tensegrity Robotic Fish's Body Waves

A bionic robotic fish based on compliant structure can excite the natural modes of vibration, thereby mimicking the body waves of real fish to generate thrust and realize undulate propulsion. The fish body wave is a result of the fish body's mechanical characteristics interacting with the surrounding fluid. Thoroughly analyzing the complex modal characteristics in such robotic fish contributes to a better understanding of the locomotion behavior, consequently enhancing the swimming performance. Therefore, the complex orthogonal decomposition (COD) method is used in this article. The traveling index is used to quantitatively describe the difference between the real and imaginary modes of the fish body wave. It is defined as the reciprocal of the condition number between the real and imaginary components. After introducing the BCF (body and/or caudal fin) the fish's body wave curves and the COD method, the structural design and parameter configuration of the tensegrity robotic fish are introduced. The complex modal characteristics of the tensegrity robotic fish and real fish are analyzed. The results show that their traveling indexes are close, with two similar complex mode shapes. Subsequently, the relationship between the traveling index and swimming performance is expressed using indicators reflecting linear correlation (correlation coefficient (Rc) and p value). Based on this correlation, a preliminary optimization strategy for the traveling index is proposed, with the potential to improve the swimming performance of the robotic fish.

Control of movement of underwater swimmers: Animals, simulated animates and swimming robots

The control of movement in living organisms represents a fundamental task that the brain has evolved to solve. One crucial aspect is how the nervous system organizes the transformation of sensory information into motor commands. These commands lead to muscle activation and subsequent animal movement, which can exhibit complex patterns. One example of such movement is locomotion, which involves the translation of the entire body through space. Central Pattern Generators (CPGs) are neuronal circuits that provide control signals for these movements. Compared to the intricate circuits found in the brain, CPGs can be simplified into networks of neurons that generate rhythmic activation, coordinating muscle movements. Since the 1990s, researchers have developed numerous models of locomotive circuits to simulate different types of animal movement, including walking, flying, and swimming. Initially, the primary goal of these studies was to construct biomimetic robots. However, it became apparent that simplified CPGs alone were not sufficient to replicate the diverse range of adaptive locomotive movements observed in living organisms. Factors such as sensory modulation, higher-level control, and cognitive components related to learning and memory needed to be considered. This necessitated the use of more complex, high-dimensional circuits, as well as novel materials and hardware, in both modeling and robotics. With advancements in high-power computing, artificial intelligence, big data processing, smart materials, and electronics, the possibility of designing a new generation of true bio-mimetic robots has emerged. These robots have the capability to imitate not only simple locomotion but also exhibit adaptive motor behavior and decision-making. This motivation serves as the foundation for the current review, which aims to analyze existing concepts and models of movement control systems. As an illustrative example, we focus on underwater movement and explore the fundamental biological concepts, as well as the mathematical and physical models that underlie locomotion and its various modulations.

Mechanosensory Feedback in Lamprey Swimming Models and Applications in the Field of Spinal Cord Regeneration

The central pattern generator (CPG) in anguilliform swimming has served as a model for examining the neural basis of locomotion. This system has been particularly valuable for the development of mathematical models. As our biological understanding of the neural basis of locomotion has expanded, so too have these models. Recently, there have been significant advancements in our understanding of the critical role that mechanosensory feedback plays in robust locomotion. This work has led to a push in the field of mathematical modeling to incorporate mechanosensory feedback into CPG models. In this perspective piece, we review advances in the development of these models and discuss how newer complex models can support biological investigation. We highlight lamprey spinal cord regeneration as an area that can both inform these models and benefit from them.

Computational Study of Stiffness-Tuning Strategies in Anguilliform Fish

Biological evidence demonstrates that fish can tune their body stiffness to improve thrust and efficiency during swimming locomotion. However, the stiffness-tuning strategies that maximize swimming speed or efficiency are still unclear. In the present study, a musculo-skeletal model of anguilliform fish is developed to study the properties of variable stiffness, in which the planar serial-parallel mechanism is used to model the body structure. The calcium ion model is adopted to simulate muscular activities and generate muscle force. Further, the relations among the forward speed, the swimming efficiency, and Young's modulus of the fish body are investigated. The results show that for certain body stiffness, the swimming speed and efficiency are increased with the tail-beat frequency until reaching the maximum value and then decreased. The peak speed and efficiency are also increased with the amplitude of muscle actuation. Anguilliform fish tend to vary their body stiffness to improve the swimming speed and efficiency at a high tail-beat frequency or small amplitude of muscle actuation. Furthermore, the midline motions of anguilliform fish are analyzed by the complex orthogonal decomposition (COD) method, and the discussions of fish motions associated with the variable body stiffness and the tail-beat frequency are also presented. Overall, the optimal swimming performance of anguilliform fish benefits from the matching relationships among the muscle actuation, the body stiffness, and the tail-beat frequency.

Proprioceptive feedback amplification restores effective locomotion in a neuromechanical model of lampreys with spinal injuries

Spinal injuries in many vertebrates can result in partial or complete loss of locomotor ability. While mammals often experience permanent loss, some nonmammals, such as lampreys, can regain swimming function, though the exact mechanism is not well understood. One hypothesis is that amplified proprioceptive (body-sensing) feedback can allow an injured lamprey to regain functional swimming even if the descending signal is lost. This study employs a multiscale, integrative, computational model of an anguilliform swimmer fully coupled to a viscous, incompressible fluid and examines the effects of amplified feedback on swimming behavior. This represents a model that analyzes spinal injury recovery by combining a closed-loop neuromechanical model with sensory feedback coupled to a full Navier-Stokes model. Our results show that in some cases, feedback amplification below a spinal lesion is sufficient to partially or entirely restore effective swimming behavior.

Dipole- and vortex sheet-based models of fish swimming

Elucidating the hydrodynamics of fish swimming is critical to identifying the processes underlying fish orientation and schooling. Due to their mathematical tractability, models based on potential flow are preferred in the study of bidirectional interactions of fish with their surroundings. Dipole-based models that assimilate fish to pairs of vortices are particularly enticing, but yet to be thoroughly validated. Here, we embark on a computational fluid dynamics (CFD) campaign informed by experimental data to validate the accuracy of dipole-based models. The locomotory patterns of a fish undergoing carangiform swimming are reconstructed from existing experimental data, which are used as inputs to CFD simulations of a fish swimming in a channel flow. We demonstrate that dipole-based models are accurate in capturing key features of the fluid flow, but cannot predict the elongated flow streamlines around the fish that are evident in CFD. To address this issue, we propose an alternative model that replaces each vortex in the pair with a sheet along the fish length. Using a pair of vortex sheets that span approximately 80% of the fish body length with a separation distance of approximately 50% of the body width, the model is successful in predicting the fluid flow around the swimming fish for a range of background flow speeds and channel widths. The proposed model shows improved accuracy at the cost of a mildly increased computational effort, thereby constituting an ideal basis for research on fish hydrodynamics.

The Core Concepts, Competencies, and Grand Challenges of Comparative Vertebrate Anatomy and Morphology

Core concepts offer coherence to the discourse of a scientific discipline and facilitate teaching by identifying large unifying themes that can be tailored to the level of the class and expertise of the instructor. This approach to teaching has been shown to encourage deeper learning that can be integrated across subdisciplines of biology and has been adopted by several other biology subdisciplines. However, Comparative Vertebrate Anatomy, although one of the oldest biological areas of study, has not had its core concepts identified. Here, we present five core concepts and seven competencies (skills) for Comparative Vertebrate Anatomy that came out of an iterative process of engagement with the broader community of vertebrate morphologists over a 3-year period. The core concepts are (A) evolution, (B) structure and function, (C) morphological development, (D) integration, and (E) human anatomy is the result of vertebrate evolution. The core competencies students should gain from the study of comparative vertebrate anatomy are (F) tree thinking, (G) observation, (H) dissection of specimens, (I) depiction of anatomy, (J) appreciation of the importance of natural history collections, (K) science communication, and (L) data integration. We offer a succinct description of each core concept and competency, examples of learning outcomes that could be used to assess teaching effectiveness, and examples of relevant resources for both instructors and students. Additionally, we pose a grand challenge to the community, arguing that the field of Comparative Vertebrate Anatomy needs to acknowledge racism, androcentrism, homophobia, genocide, slavery, and other influences in its history and address their lingering effects in order to move forward as a thriving discipline that is inclusive of all students and scientists and continues to generate unbiased knowledge for the betterment of humanity. Despite the rigorous process used to compile these core concepts and competencies, we anticipate that they will serve as a framework for an ongoing conversation that ensures Comparative Vertebrate Anatomy remains a relevant field in discovery, innovation, and training of future generations of scientists.

Hydrodynamic Interaction of Two Self-Propelled Fish Swimming in a Tandem Arrangement

Collective locomotion in biological systems is ubiquitous and attracts much attention, and there are complex hydrodynamics involved. The hydrodynamic interaction for fish schooling is examined using two-dimensional numerical simulations of a pair of self-propelled swimming fish in this paper. The effects of different parameters on swimming speed gain and energy-saving efficiency are investigated by adjusting swimming parameters (initial separation distance d0, tail beat amplitude A, body wavelength λ, and period of oscillation T) at different phase difference δϕ between two fish. The hydrodynamic interaction performance of fish swimming in a tandem arrangement is analyzed with the help of the instantaneous vorticity contours, pressure contours, and mean work done. Using elementary hydrodynamic arguments, a unifying mechanistic principle, which characterizes the fish locomotion by deriving a scaling relation that links swimming speed u to body kinematics (A, T, and λ), arrangement of formation (d0), and fluid properties (kinematic viscosity ν), is revealed. It is shown that there are some certain scaling laws between similarity criterion number (Reynolds number (Re) and Strouhal number (St)) and energy-consuming coefficient (CE) under different parameters (Δ). In particular, a generality in the relationships of St–Re and CE–(Re ·Δ) can emerge despite significant disparities in locomotory performance.

Fish-inspired segment models for undulatory steady swimming

Many aquatic animals swim by undulatory body movements and understanding the diversity of these movements could unlock the potential for designing better underwater robots. Here, we analyzed the steady swimming kinematics of a diverse group of fish species to investigate whether their undulatory movements can be represented using a series of interconnected multi-segment models, and if so, to identify the key factors driving the segment configuration of the models. Our results show that the steady swimming kinematics of fishes can be described successfully using parsimonious models, 83% of which had fewer than five segments. In these models, the anterior segments were significantly longer than the posterior segments, and there was a direct link between segment configuration and swimming kinematics, body shape, and Reynolds number. The models representing eel-like fishes with elongated bodies and fishes swimming at high Reynolds numbers had more segments and less segment length variability along the body than the models representing other fishes. These fishes recruited their anterior bodies to a greater extent, initiating the undulatory wave more anteriorly. Two shape parameters, related to axial and overall body thickness, predicted segment configuration with moderate to high success rate. We found that head morphology was a good predictor of its segment length. While there was a large variation in head segments, the length of tail segments was similar across all models. Given that fishes exhibited variable caudal fin shapes, the consistency of tail segments could be a result of an evolutionary constraint tuned for high propulsive efficiency. The bio-inspired multi-segment models presented in this study highlight the key bending points along the body and can be used to decide on the placement of actuators in fish-inspired robots, to model hydrodynamic forces in theoretical and computational studies, or for predicting muscle activation patterns during swimming.

Robotics as a Comparative Method in Ecology and Evolutionary Biology

Comparative biologists have typically used one or more of the following methods to assist in evaluating the proposed functional and performance significance of individual traits: comparative phylogenetic analysis, direct interspecific comparison among species, genetic modification, experimental alteration of morphology (for example by surgically modifying traits), and ecological manipulation where individual organisms are transplanted to a different environment. But comparing organisms as the endpoints of an evolutionary process involves the ceteris paribus assumption: that all traits other than the one(s) of interest are held constant. In a properly controlled experimental study, only the variable of interest changes among the groups being compared. The theme of this paper is that the use of robotic or mechanical models offers an additional tool in comparative biology that helps to minimize the effect of uncontrolled variables by allowing direct manipulation of the trait of interest against a constant background. The structure and movement pattern of mechanical devices can be altered in ways not possible in studies of living animals, facilitating testing hypotheses of the functional and performance significant of individual traits. Robotic models of organismal design are particularly useful in three arenas: (1) controlling variation to allow modification only of the trait of interest, (2) the direct measurement of energetic costs of individual traits, and (3) quantification of the performance landscape. Obtaining data in these three areas is extremely difficult through the study of living organisms alone, and the use of robotic models can reveal unexpected effects. Controlling for all variables except for the length of a swimming flexible object reveals substantial non-linear effects that vary with stiffness. Quantification of the swimming performance surface reveals that there are two peaks with comparable efficiency, greatly complicating the inference of performance from morphology alone. Organisms and their ecological interactions are complex, and dissecting this complexity to understand the effects of individual traits is a grand challenge in ecology and evolutionary biology. Robotics has great promise as a "comparative method," allowing better-controlled comparative studies to analyze the many interacting elements that make up complex behaviors, ecological interactions, and evolutionary histories.

Hydrodynamical Fingerprint of a Neighbour in a Fish Lateral Line

For fish, swimming in group may be favorable to individuals. Several works reported that in a fish school, individuals sense and adjust their relative position to prevent collisions and maintain the group formation. Also, from a hydrodynamic perspective, relative-position and kinematic synchronisation between adjacent fish may considerably influence their swimming performance. Fish may sense the relative-position and tail-beat phase difference with their neighbors using both vision and the lateral-line system, however, when swimming in dark or turbid environments, visual information may become unavailable. To understand how lateral-line sensing can enable fish to judge the relative-position and phase-difference with their neighbors, in this study, based on a verified three-dimensional computational fluid dynamics approach, we simulated two fish swimming adjacently with various configurations. The lateral-line signal was obtained by sampling the surface hydrodynamic stress. The sensed signal was processed by Fast Fourier Transform (FFT), which is robust to turbulence and environmental flow. By examining the lateral-line pressure and shear-stress signals in the frequency domain, various states of the neighboring fish were parametrically identified. Our results reveal that the FFT-processed lateral-line signals in one fish may potentially reflect the relative-position, phase-differences, and the tail-beat frequency of its neighbor. Our results shed light on the fluid dynamical aspects of the lateral-line sensing mechanism used by fish. Furthermore, the presented approach based on FFT is especially suitable for applications in bioinspired swimming robotics. We provide suggestions for the design of artificial systems consisting of multiple stress sensors for robotic fish to improve their performance in collective operation.

Eulerian simulation of complex suspensions and biolocomotion in three dimensions

Fluid–structure interactions are ubiquitous in many natural and man-made environments. They are difficult to study analytically, and therefore accurate and flexible computational methods are an indispensable tool in the field. Typically, fluids are simulated with a fixed background computational mesh, whereas a solid is simulated with a mesh that moves with it, making it challenging to couple the two. Here we develop a three-dimensional computational method where both fluid and solid can be represented on a fixed computational grid, which simplifies the coupling between the two phases considerably. Our method can simulate scenarios that are difficult to do with other methods, such as complex suspensions containing soft particles that are both heavier and lighter than the fluid.

We present a numerical method specifically designed for simulating three-dimensional fluid–structure interaction (FSI) problems based on the reference map technique (RMT). The RMT is a fully Eulerian FSI numerical method that allows fluids and large-deformation elastic solids to be represented on a single fixed computational grid. This eliminates the need for meshing complex geometries typical in other FSI approaches and greatly simplifies the coupling between fluid and solids. We develop a three-dimensional implementation of the RMT, parallelized using the distributed memory paradigm, to simulate incompressible FSI with neo-Hookean solids. As part of our method, we develop a field extrapolation scheme that works efficiently in parallel. Through representative examples, we demonstrate the method’s suitability in investigating many-body and active systems, as well as its accuracy and convergence. The examples include settling of a mixture of heavy and buoyant soft ellipsoids, lid-driven cavity flow containing a soft sphere, and swimmers actuated via active stress.

Tunable stiffness in fish robotics: mechanisms and advantages

One of the emerging themes of fish-inspired robotics is flexibility. Adding flexibility to the body, joints, or fins of fish-inspired robots can significantly improve thrust and/or efficiency during locomotion. However, the optimal stiffness depends on variables such as swimming speed, so there is no one 'best' stiffness that maximizes efficiency in all conditions. Fish are thought to solve this problem by using muscular activity to tune their body and fin stiffness in real-time. Inspired by fish, some recent robots sport polymer actuators, adjustable leaf springs, or artificial tendons that tune stiffness mechanically. Models and water channel tests are providing a theoretical framework for stiffness-tuning strategies that devices can implement. The strategies can be thought of as analogous to car transmissions, which allow users to improve efficiency by tuning gear ratio with driving speed. We provide an overview of the latest discoveries about (1) the propulsive benefits of flexibility, particularlytunableflexibility, and (2) the mechanisms and strategies that fish and fish-inspired robots use to tune stiffness while swimming.

Innovative Energy-Saving Propulsion System for Low-Speed Biomimetic Underwater Vehicles

This article covers research on an innovative propulsion system design for a Biomimetic Unmanned Underwater Vehicle (BUUV) operating at low speeds. The experiment was conducted on a laboratory test water tunnel equipped with specialised sensor equipment to assess the Fluid-Structure Interaction (FSI) and energy consumption of two different types of propulsion systems. The experimental data contrast the undulating with the drag-based propulsion system. The additional joint in the drag-based propulsion system is intended to increase thrust and decrease energy input. The tests were conducted at a variety of fins oscillation frequencies and fluid velocities. The experiments demonstrate that, in the region of low-speed forward movement, the efficiency of the propulsion system with the additional joint is greater.

Swimming kinematics and performance of spinal transected lampreys with different levels of axon regeneration

Axon regeneration is critical for restoring neural function after spinal cord injury. This has prompted a series of studies on the neural and functional recovery of lampreys after spinal cord transection. Despite this, there are still many basic questions remaining about how much functional recovery depends on axon regeneration. Our goal was to examine how swimming performance is related to degree of axon regeneration in lampreys recovering from spinal cord transection by quantifying the relationship between swimming performance and percent axon regeneration of transected lampreys after 11 weeks of recovery. We found that while swimming speeds varied, they did not relate to percent axon regeneration. In fact, swimming speeds were highly variable within individuals, meaning that most individuals could swim at both moderate and slow speeds, regardless of percent axon regeneration. However, none of the transected individuals were able to swim as fast as the control lampreys. To swim fast, control lampreys generated high amplitude body waves with long wavelengths. Transected lampreys generated body waves with lower amplitude and shorter wavelengths than controls, and to compensate, transected lampreys increased their wave frequencies to swim faster. As a result, transected lampreys had significantly higher frequencies than control lampreys at comparable swimming velocities. These data suggest that the control lampreys swam more efficiently than transected lampreys. In conclusion, there appears to be a minimal recovery threshold in terms of percent axon regeneration required for lampreys to be capable of swimming; however, there also seems to be a limit to how much they can behaviorally recover.

Summary: Lampreys that have recovered from having their spinal cords transected do not fully regain swimming abilities and are not able to swim as efficiently as non-transected lampreys.

Increasing Viscosity Helps Explain Locomotor Control in Swimming Polypterus senegalus

Locomotion relies on the successful integration of sensory information to adjust brain commands and basic motor rhythms created by central pattern generators. It is not clearly understood how altering the sensory environment impacts control of locomotion. In an aquatic environment, mechanical sensory feedback to the animal can be readily altered by adjusting water viscosity. Computer modeling of fish swimming systems shows that, without sensory feedback, high viscosity systems dampen kinematic output despite similar motor control input. We recorded muscle activity and kinematics of six Polypterus senegalus in four different viscosities of water from 1 cP (normal water) to 40 cP. In high viscosity, P. senegalus exhibit increased body curvature, body wave speed, and body and pectoral fin frequency during swimming. These changes are the result of increased muscle activation intensity and maintain voluntary swimming speed. Unlike the sensory-deprived model, intact sensory feedback allows fish to adjust swimming motor control and kinematic output in high viscous water but maintain typical swimming coordination.

Body and Tail Coordination in the Bluespot Salamander (Ambystoma laterale) During Limb Regeneration

Animals are incredibly good at adapting to changes in their environment, a trait envied by most roboticists. Many animals use different gaits to seamlessly transition between land and water and move through non-uniform terrains. In addition to adjusting to changes in their environment, animals can adjust their locomotion to deal with missing or regenerating limbs. Salamanders are an amphibious group of animals that can regenerate limbs, tails, and even parts of the spinal cord in some species. After the loss of a limb, the salamander successfully adjusts to constantly changing morphology as it regenerates the missing part. This quality is of particular interest to roboticists looking to design devices that can adapt to missing or malfunctioning components. While walking, an intact salamander uses its limbs, body, and tail to propel itself along the ground. Its body and tail are coordinated in a distinctive wave-like pattern. Understanding how their bending kinematics change as they regrow lost limbs would provide important information to roboticists designing amphibious machines meant to navigate through unpredictable and diverse terrain. We amputated both hindlimbs of blue-spotted salamanders (Ambystoma laterale) and measured their body and tail kinematics as the limbs regenerated. We quantified the change in the body wave over time and compared them to an amphibious fish species, Polypterus senegalus. We found that salamanders in the early stages of regeneration shift their kinematics, mostly around their pectoral girdle, where there is a local increase in undulation frequency. Amputated salamanders also show a reduced range of preferred walking speeds and an increase in the number of bending waves along the body. This work could assist roboticists working on terrestrial locomotion and water to land transitions.

Emergent metachronal waves using tension-driven, fluid-structure interaction models of tomopterid parapodia

Metachronal waves are ubiquitous in propulsive and fluid transport systems across many different scales and morphologies in the biological world. Tomopterids are a soft-bodied, holopelagic polychaete that use metachrony with their flexible, gelatinous parapodia to deftly navigate the midwater ocean column that they inhabit. In the following study, we develop a three-dimensional, fluid-structure interaction model of a tomopterid parapodium to explore the emergent metachronal waves formed from the interplay of passive body elasticity, active muscular tension, and hydrodynamic forces. After introducing our model, we examine the effects that varying material properties have on the stroke of an individual parapodium. We then explore the temporal dynamics when multiple parapodia are placed sequentially and how differences in the phase can alter the collective kinematics and resulting flow field.

Automatic segmentation of fish midlines for optimizing robot design

While fish use continuous and flexible bodies to propel themselves, fish robots are often made from interconnected segments. How many segments do robots need to represent fish movements accurately? We propose a new method to automatically determine parsimonious robot models from actual fish data. We first identify key bending points (i.e., joint positions) along the body and then study the concerted movement of the segments so that the difference between actual fish and modelled bending kinematics is minimized. To demonstrate the utility of our method, we analyse the steady swimming kinematics of 10 morphologically distinct fish species. Broadly classified as sub-carangiform (e.g., rainbow trout) and carangiform (e.g., crevalle jack) swimmers, these species exhibit variations in the way they undulate when traditional parameters (including head and tail beat amplitudes, body wavelength and maximum curvature along the body) are considered. We show that five segments are sufficient to describe the kinematics with at least 99% accuracy. For optimal performance, segments should progressively get shorter towards the tail. We also show that locations where bending moments are applied vary among species, possibly because of differences in morphology. More specifically, we find that wider fish have shorter head segments. We discover that once bending points are factored in, the kinematics differences observed in these species collapse into a single undulatory pattern. The amplitude and timing of how body segments move entirely depend on their respective joint positions along the body. Head and body segments are also coupled in a timely manner, which depends on the position of the most anterior joint. Our findings provide a mechanistic understanding of how morphology relates to kinematics and highlight the importance of head control, which is often overlooked in current robot designs.

A Tail of Four Fishes: An Analysis of Kinematics and Material Properties of Elongate Fishes

The elongate body plan is present in many groups of fishes, and this morphology dictates functional consequences seen in swimming behavior. Previous work has shown that increasing the number of vertebrae, or decreasing the intervertebral joint length, in a fixed length artificial system increases stiffness. Tails with increased stiffness can generate more power from tail beats, resulting in an increased mean swimming speed. This demonstrates the impacts of morphology on both material properties and kinematics, establishing mechanisms for form contributing to function. Here, we wanted to investigate relationships between form and ecological function, such as differences in dietary strategies and habitat preferences among fish species. This study aims to characterize and compare the kinematics, material properties, and vertebral morphology of four species of elongate fishes: Anoplarchus insignis, Anoplarchus purpurescens, Xiphister atropurpureus, and Xiphister mucosus. We hypothesized that these properties would differ among the four species due to their differential ecological niches. To calculate kinematic variables, we filmed these fishes swimming volitionally. We also measured body stiffness by bending the abdominal and tail regions of sacrificed individuals in different stages of dissection (whole body, removed skin, and removed muscle). Finally, we counted the number of vertebrae from CT scans of each species to quantify vertebral morphology. Principal component and linear discriminant analyses suggested that the elongate fish species can be distinguished from one another by their material properties, morphology, and swimming kinematics. With this information combined, we can draw connections between the physical properties of the fishes and their ecological niches.

Neuromechanical wave resonance in jellyfish swimming

For organisms to have robust locomotion, their neuromuscular organization must adapt to constantly changing environments. In jellyfish, swimming robustness emerges when marginal pacemakers fire action potentials throughout the bell's motor nerve net, which signals the musculature to contract. The speed of the muscle activation wave is dictated by the passage times of the action potentials. However, passive elastic material properties also influence the emergent kinematics, with time scales independent of neuromuscular organization. In this multimodal study, we examine the interplay between these two time scales during turning. A three-dimensional computational fluid-structure interaction model of a jellyfish was developed to determine the resulting emergent kinematics, using bidirectional muscular activation waves to actuate the bell rim. Activation wave speeds near the material wave speed yielded successful turns, with a 76-fold difference in turning rate between the best and worst performers. Hyperextension of the margin occurred only at activation wave speeds near the material wave speed, suggesting resonance. This hyperextension resulted in a 34-fold asymmetry in the circulation of the vortex ring between the inside and outside of the turn. Experimental recording of the activation speed confirmed that jellyfish actuate within this range, and flow visualization using particle image velocimetry validated the corresponding fluid dynamics of the numerical model. This suggests that neuromechanical wave resonance plays an important role in the robustness of an organism's locomotory system and presents an undiscovered constraint on the evolution of flexible organisms. Understanding these dynamics is essential for developing actuators in soft body robotics and bioengineered pumps.

Tunabot Flex: a tuna-inspired robot with body flexibility improves high-performance swimming

Tunas are flexible, high-performance open ocean swimmers that operate at high frequencies to achieve high swimming speeds. Most fish-like robotic systems operate at low frequencies (≤3 Hz) resulting in low swim speeds (≤1.5 body lengths per second), and the cost of transport (COT) is often one to four orders of magnitude higher than that of tunas. Furthermore, the impact of body flexibility on high-performance fish swimming remains unknown. Here we design and test a research platform based on yellowfin tuna (Thunnus albacares) to investigate the role of body flexibility and to close the performance gap between robotic and biological systems. This single-motor platform, termed Tunabot Flex, measures 25.5 cm in length. Flexibility is varied through joints in the tail to produce three tested configurations. We find that increasing body flexibility improves self-propelled swimming speeds on average by 0.5 body lengths per second while reducing the minimum COT by 53%. The most flexible configuration swims 4.60 body lengths per second with a tail beat frequency of 8.0 Hz and a COT measuring 18.4 J kg^-1m^-1. We then compare these results in addition to the midline kinematics, stride length, and Strouhal number with yellowfin tuna data. The COT of Tunabot Flex's most flexible configuration is less than a half-order of magnitude greater than that of yellowfin tuna across all tested speeds. Tunabot Flex provides a new baseline for the development of future bio-inspired underwater vehicles that aim to explore a fish-like, high-performance space and close the gap between engineered robotic systems and fish swimming ability.

Fish without Tail Fins-Exploring the Function of Tail Morphology of the First Vertebrates

We use a series of hydrodynamic experiments on abstracted models to explore whether primitive vertebrates may have swum under various conditions without a clearly-differentiated tail fin. Cambrian vertebrates had post-anal stubby tails, some had single dorsal and ventral fins, but none had yet evolved a clearly differentiated caudal fin typical of post-Cambrian fishes, and must have relied on their long and flexible laterally-compressed bodies for locomotion, i.e., by bending their bodies side-to-side in order to propagate waves from head to tail. We approach this problem experimentally based on an abstracted model of Metaspriggina walcotti from the 506-million-year old Burgess Shale by using oscillating thin flexible plates while varying the tail fin geometry from rectangular to uniform, and finally to a no tail-fin condition. Despite a missing tail fin, this study supports the observation that the abstracted Metaspriggina model can generate a strong propulsive force in cruise conditions, both away from, and near the sea bed (in ground effect). However, when the abstracted Metaspriggina model moves in ground effect, a weaker performance is observed, indicating that Metaspriggina may not necessarily have been optimized for swimming near the sea bed. When considering acceleration from rest, we find that the Metaspriggina model's performance is not significantly different from other morphological models (abstracted truncate tail and abstracted heterocercal tail). Statistical analysis shows that morphological parameters, swimming modes, and ground effect all play significant roles in thrust performance. While the exact relationships of Cambrian vertebrates are still debated, as agnathans, they share some general characteristics with modern cyclostomes, in particular an elongate body akin to lampreys. Lampreys, as anguilliform swimmers, are considered to be some of the most efficient swimmers using a particular type of suction thrust induced by the traveling body wave as it travels from head to tail. Our current experiments suggest that Metaspriggina's ability in acceleration from rest, through possibly a similar type of suction thrust, which is defined as the ability to generate low pressure on upstream facing sections of the body, might have evolved early in response to increasing predator pressure during the Cambrian Explosion.

Reproducing Five Motor Behaviors in a Salamander Robot With Virtual Muscles and a Distributed CPG Controller Regulated by Drive Signals and Proprioceptive Feedback

Diverse locomotor behaviors emerge from the interactions between the spinal central pattern generator (CPG), descending brain signals and sensory feedback. Salamander motor behaviors include swimming, struggling, forward underwater stepping, and forward and backward terrestrial stepping. Electromyographic and kinematic recordings of the trunk show that each of these five behaviors is characterized by specific patterns of muscle activation and body curvature. Electrophysiological recordings in isolated spinal cords show even more diverse patterns of activity. Using numerical modeling and robotics, we explored the mechanisms through which descending brain signals and proprioceptive feedback could take advantage of the flexibility of the spinal CPG to generate different motor patterns. Adapting a previous CPG model based on abstract oscillators, we propose a model that reproduces the features of spinal cord recordings: the diversity of motor patterns, the correlation between phase lags and cycle frequencies, and the spontaneous switches between slow and fast rhythms. The five salamander behaviors were reproduced by connecting the CPG model to a mechanical simulation of the salamander with virtual muscles and local proprioceptive feedback. The main results were validated on a robot. A distributed controller was used to obtain the fast control loops necessary for implementing the virtual muscles. The distributed control is demonstrated in an experiment where the robot splits into multiple functional parts. The five salamander behaviors were emulated by regulating the CPG with two descending drives. Reproducing the kinematics of backward stepping and struggling however required stronger muscle contractions. The passive oscillations observed in the salamander's tail during forward underwater stepping could be reproduced using a third descending drive of zero to the tail oscillators. This reduced the drag on the body in our hydrodynamic simulation. We explored the effect of local proprioceptive feedback during swimming and forward terrestrial stepping. We found that feedback could replace or reduce the need for different drives in both cases. It also reduced the variability of intersegmental phase lags toward values appropriate for locomotion. Our work suggests that different motor behaviors do not require different CPG circuits: a single circuit can produce various behaviors when modulated by descending drive and sensory feedback.

A semi-automated finite difference mesh creation method for use with immersed boundary software IB2d and IBAMR

Numerous fluid-structure interaction problems in biology have been investigated using the immersed boundary method. The advantage of this method is that complex geometries, e.g., internal or external morphology, can easily be handled without the need to generate matching grids for both the fluid and the structure. Consequently, the difficulty of modeling the structure lies often in discretizing the boundary of the complex geometry (morphology). Both commercial and open source mesh generators for finite element methods have long been established; however, the traditional immersed boundary method is based on a finite difference discretization of the structure. Here we present a software library for obtaining finite difference discretizations of boundaries for direct use in the 2D immersed boundary method. This library provides tools for extracting such boundaries as discrete mesh points from digital images. We give several examples of how the method can be applied that include passing flow through the veins of insect wings, within lymphatic capillaries, and around starfish using open-source immersed boundary software.

Decoding the Relationships between Body Shape, Tail Beat Frequency, and Stability for Swimming Fish

As fish swim through a fluid environment, they must actively use their fins in concert to stabilize their motion and have a robust form of locomotion. However, there is little knowledge of how these forces act on the fish body. In this study, we employ a 3D immersed boundary model to decode the relationship between roll, pitch, and yaw of the fish body and the driving forces acting on flexible fish bodies. Using bluegill sunfish as our representative geometry, we first examine the role of an actuating torque on the stability of the fish model, with a torque applied at the head of the unconstrained fish body. The resulting kinematics is a product of the passive elasticity, fluid forces, and driving torque. We then examine a constrained model to understand the role that fin geometry, body elasticity, and frequency play on the range of corrective forces acting on the fish. We find non-monotonic behavior with respect to frequency, suggesting that the effective flexibility of the fins play an important role in the swimming performance.

Edible additive effects on zebrafish cardiovascular functionality with hydrodynamic assessment

Food coloring is often used as a coloring agent in foods, medicines and cosmetics, and it was reported to have certain carcinogenic and mutagenic effects in living organisms. Investigation of physiological parameters using zebrafish is a promising methodology to understand disease biology and drug toxicity for various drug discovery on humans. Zebrafish (Danio rerio) is a well-acknowledged model organism with combining assets such as body transparency, small size, low cost of cultivation, and high genetic homology with humans and is used as a specimen tool for the in-vivo throughput screening approach. In addition, recent advances in microfluidics show a promising alternative for zebrafish manipulation in terms of drug administration and extensive imaging capability. This pilot work highlighted the design and development of a microfluidic detection platform for zebrafish larvae through investigating the effects of food coloring on cardiovascular functionality and pectoral fin swing ability. The zebrafish embryos were exposed to the Cochineal Red and Brilliant Blue FCF pigment solution in a concentration of (0.02‰, 0.2‰) cultured in the laboratory from the embryo stage to hatching and development until 9 days post fertilization (d.p.f.). In addition, zebrafish swimming behaviors in terms of pectoral fin beating towards the toxicity screening were further studied by visualizing the induced flow field. It was evidenced that Cochineal Red pigment at a concentration of 0.2‰ not only significantly affected the zebrafish pectoral fin swing behavior, but also significantly increased the heart rate of juvenile fish. The higher concentration of Brilliant Blue FCF pigment (0.2%) increased heart rate during early embryonic stages of zebrafish. However, zebrafish exposed to food coloring did not show any significant changes in cardiac output. The applications of this proposed platform can be further extended towards observing the neurobiological/hydrodynamic behaviors of zebrafish larvae for practical applications in drug tests.

Diving into a Simple Anguilliform Swimmer’s Sensitivity

Computational models of aquatic locomotion range from modest individual simple swimmers in 2D to sophisticated 3D multi-swimmer models that attempt to parse collective behavioral dynamics. Each of these models contain a multitude of model input parameters to which its outputs are inherently dependent, that is, various performance metrics. In this work, the swimming performance’s sensitivity to parameters is investigated for an idealized, simple anguilliform swimming model in 2D. The swimmer considered here propagates forward by dynamically varying its body curvature, similar to motion of a Caenorhabditis elegans. The parameter sensitivities were explored with respect to the fluid scale (Reynolds number), stroke (undulation) frequency, as well as a kinematic parameter controlling the velocity and acceleration of each upstroke and downstroke. The input Reynolds number and stroke frequencies sampled were from [450, 2200] and [1, 3] Hz, respectively. In total, 5000 fluid–structure interaction simulations were performed, each with a unique parameter combination selected via a Sobol sequence, in order to conduct global sensitivity analysis. Results indicate that the swimmer’s performance is most sensitive to variations in its stroke frequency. Trends in swimming performance were discovered by projecting the performance data onto particular 2D subspaces. Pareto-like optimal fronts were identified. This work is a natural extension of the parameter explorations of the same model from Battista in 2020.

How zebrafish turn: analysis of pressure force dynamics and mechanical work

Whereas many fishes swim steadily, zebrafish regularly exhibit unsteady burst-and-coast swimming, which is characterized by repeated sequences of turns followed by gliding periods. Such a behavior offers the opportunity to investigate the hypothesis that negative mechanical work occurs in posterior regions of the body during early phases of the turn near the time of maximal body curvature. Here, we used a modified particle image velocimetry (PIV) technique to obtain high-resolution flow fields around the zebrafish body during turns. Using detailed swimming kinematics coupled with body surface pressure computations, we estimated fluid-structure interaction forces and the pattern of forces and torques along the body during turning. We then calculated the mechanical work done by each body segment. We used estimated patterns of positive and negative work along the body to evaluate the hypothesis (based on fish midline kinematics) that the posterior body region would experience predominantly negative work. Between 10% and 20% of the total mechanical work was done by the fluid on the body (negative work), and negative work was concentrated in the anterior and middle areas of the body, not along the caudal region. Energetic costs of turning were calculated by considering the sum of positive and negative work and were compared with previous metabolic estimates of turning energetics in fishes. The analytical workflow presented here provides a rigorous way to quantify hydrodynamic mechanisms of fish locomotion and facilitates the understanding of how body kinematics generate locomotor forces in freely swimming fishes.

Swimming and defence: competing needs across ontogeny in armoured fishes (Agonidae)

Biological armours are potent model systems for understanding the complex series of competing demands on protective exoskeletons; after all, armoured organisms are the product of millions of years of refined engineering under the harshest conditions. Fishes are no strangers to armour, with various types of armour plating common to the 400-500 Myr of evolution in both jawed and jawless fishes. Here, we focus on the poachers (Agonidae), a family of armoured fishes native to temperate waters of the Pacific rim. We examined armour morphology, body stiffness and swimming performance in the northern spearnose poacher (Agonopsis vulsa) over ontogeny. As juveniles, these fishes make frequent nocturnal forays into the water column in search of food, while heavily armoured adults are bound to the benthos. Most armour dimensions and density increase with body length, as does body stiffness. Juvenile poachers have enlarged spines on their armour whereas adults invest more mineral in armour plate bases. Adults are stiffer and accelerate faster than juveniles with an anguilliform swimming mode. Subadults more closely approximate adults more than smaller juveniles, with regards to both swimming and armour mechanics. Poacher armour serves multiple functions over ontogeny, from facilitating locomotion, slowing sinking and providing defence.

Experimental-numerical method for calculating bending moments in swimming fish shows that fish larvae control undulatory swimming with simple actuation

Most fish swim with body undulations that result from fluid-structure interactions between the fish's internal tissues and the surrounding water. Gaining insight into these complex fluid-structure interactions is essential to understand how fish swim. To this end, we developed a dedicated experimental-numerical inverse dynamics approach to calculate the lateral bending moment distributions for a large-amplitude undulatory swimmer that moves freely in three-dimensional space. We combined automated motion tracking from multiple synchronised high-speed video sequences, computation of fluid dynamic stresses on the swimmer's body from computational fluid dynamics, and bending moment calculations using these stresses as input for a novel beam model of the body. The bending moment, which represent the system's net actuation, varies over time and along the fish's central axis due to muscle actions, passive tissues, inertia, and fluid dynamics. Our three-dimensional analysis of 113 swimming events of zebrafish larvae ranging in age from 3 to 12 days after fertilisation shows that these bending moment patterns are not only relatively simple but also strikingly similar throughout early development and from fast starts to periodic swimming. This suggests that fish larvae may produce and adjust swimming movements relatively simply, yet effectively, while restructuring their neuromuscular control system throughout their rapid development.

Hydrodynamics of Biomimetic Marine Propulsion and Trends in Computational Simulations

The aim of the present paper is to provide the state of the works in the field of hydrodynamics and computational simulations to analyze biomimetic marine propulsors. Over the last years, many researchers postulated that some fish movements are more efficient and maneuverable than traditional rotary propellers, and the most relevant marine propulsors which mimic fishes are shown in the present work. Taking into account the complexity and cost of some experimental setups, numerical models offer an efficient, cheap, and fast alternative tool to analyze biomimetic marine propulsors. Besides, numerical models provide information that cannot be obtained using experimental techniques. Since the literature about trends in computational simulations is still scarce, this paper also recalls the hydrodynamics of the swimming modes occurring in fish and summarizes the more relevant lines of investigation of computational models.

Fish-like aquatic propulsion studied using a pneumatically-actuated soft-robotic model

Fish locomotion is characterized by waves of muscle electrical activity that proceed from head to tail, and result in an undulatory pattern of body bending that generates thrust during locomotion. Isolating the effects of parameters like body stiffness, co-activation between the right and left sides of the body, and frequency on thrust generation has proven to be difficult in live fishes. We use a pneumatically-actuated fish-like model to investigate how these parameters affect locomotor force generation. We measure thrust as well as side forces and torques generated during propulsion. Using a statistical linear model we examine the effects of input parameter combinations on thrust generation. We show that both stiffness and frequency substantially affect swimming kinematics, and that there are complex interactive effects of these two parameters on thrust. The stiffer the backbone the more impact that increasing frequency has on thrust production. For stiffer models, increasing frequency resulted in higher values for both thrust and lateral forces. Large side forces reduce swimming efficiency but this effect could be mitigated by decreasing undulatory wavelength and allowing appropriate phasing of left and right body movements to reduce amplitudes of side force.

In vitro virtual reality: an anatomically explicit musculoskeletal simulation powered by in vitro muscle using closed-loop tissue-software interaction

Muscle force-length dynamics are governed by intrinsic contractile properties, motor stimulation and mechanical load. Although intrinsic properties are well characterised, physiologists lack in vitro instrumentation to account for combined effects of limb inertia, musculoskeletal architecture and contractile dynamics. We introduce in vitro virtual reality (in vitro-VR) which enables in vitro muscle tissue to drive a musculoskeletal jumping simulation. In hardware, muscle force from a frog plantaris was transmitted to a software model where joint torques, inertia and ground reaction forces were computed to advance the simulation at 1 kHz. To close the loop, simulated muscle strain was returned to update in vitro length. We manipulated (1) stimulation timing and (2) the virtual muscle's anatomical origin. This influenced interactions among muscular, inertial, gravitational and contact forces dictating limb kinematics and jump performance. We propose that in vitro-VR can be used to illustrate how neuromuscular control and musculoskeletal anatomy influence muscle dynamics and biomechanical performance.

Surprising simplicities and syntheses in limbless self-propulsion in sand

Animals moving on and in fluids and solids move their bodies in diverse ways to generate propulsion and lift forces. In fluids, animals can wiggle, stroke, paddle or slap, whereas on hard frictional terrain, animals largely engage their appendages with the substrate to avoid slip. Granular substrates, such as desert sand, can display complex responses to animal interactions. This complexity has led to locomotor strategies that make use of fluid-like or solid-like features of this substrate, or combinations of the two. Here, we use examples from our work to demonstrate the diverse array of methods used and insights gained in the study of both surface and subsurface limbless locomotion in these habitats. Counterintuitively, these seemingly complex granular environments offer certain experimental, theoretical, robotic and computational advantages for studying terrestrial movement, with the potential for providing broad insights into morphology and locomotor control in fluids and solids, including neuromechanical control templates and morphological and behavioral evolution. In particular, granular media provide an excellent testbed for a locomotion framework called geometric mechanics, which was introduced by particle physicists and control engineers in the last century, and which allows quantitative analysis of alternative locomotor patterns and morphology to test for control templates, optimality and evolutionary alternatives. Thus, we posit that insights gained from movement in granular environments can be translated into principles that have broader applications across taxa, habitats and movement patterns, including those at microscopic scales.

Quantification of the influence of drugs on zebrafish larvae swimming kinematics and energetics

The use of zebrafish larvae has aroused wide interest in the medical field for its potential role in the development of new therapies. The larvae grow extremely quickly and the embryos are nearly transparent which allows easy examination of its internal structures using fluorescent imaging techniques. Medical treatment of zebrafish larvae can directly influence its swimming behaviours. These behaviour changes are related to functional changes of central nervous system and transformations of the zebrafish body such as muscle mechanical power and force variation, which cannot be measured directly by pure experiment observation. To quantify the influence of drugs on zebrafish larvae swimming behaviours and energetics, we have developed a novel methodology to exploit intravital changes based on observed zebrafish locomotion. Specifically, by using an in-house MATLAB code to process the recorded live zebrafish swimming video, the kinematic locomotion equation of a 3D zebrafish larvae was obtained, and a customised Computational Fluid Dynamics tool was used to solve the fluid flow around the fish model which was geometrically the same as experimentally tested zebrafish. The developed methodology was firstly verified against experiment, and further applied to quantify the fish internal body force, torque and power consumption associated with a group of normal zebrafish larvae vs. those immersed in acetic acid and two neuroactive drugs. As indicated by our results, zebrafish larvae immersed in 0.01% acetic acid display approximately 30% higher hydrodynamic power and 10% higher cost of transport than control group. In addition, 500 μM diphenylhydantoin significantly decreases the locomotion activity for approximately 50% lower hydrodynamic power, whereas 100 mg/L yohimbine has not caused any significant influences on 5 dpf zebrafish larvae locomotion. The approach has potential to evaluate the influence of drugs on the aquatic animal’s behaviour changes and thus support the development of new analgesic and neuroactive drugs.

From single neurons to behavior in the jellyfish Aurelia aurita

Jellyfish nerve nets provide insight into the origins of nervous systems, as both their taxonomic position and their evolutionary age imply that jellyfish resemble some of the earliest neuron-bearing, actively-swimming animals. Here, we develop the first neuronal network model for the nerve nets of jellyfish. Specifically, we focus on the moon jelly Aurelia aurita and the control of its energy-efficient swimming motion. The proposed single neuron model disentangles the contributions of different currents to a spike. The network model identifies factors ensuring non-pathological activity and suggests an optimization for the transmission of signals. After modeling the jellyfish's muscle system and its bell in a hydrodynamic environment, we explore the swimming elicited by neural activity. We find that different delays between nerve net activations lead to well-controlled, differently directed movements. Our model bridges the scales from single neurons to behavior, allowing for a comprehensive understanding of jellyfish neural control of locomotion.

Energetics of optimal undulatory swimming organisms

Energy consumption is one of the primary considerations in animal locomotion. In swimming locomotion, a number of questions related to swimming energetics of an organism and how the energetic quantities scale with body size remain open, largely due to the difficulties with modeling and measuring the power production and consumption. Based on a comprehensive theoretical framework that incorporates cyclic muscle behavior, structural dynamics and swimming hydrodynamics, we perform extensive computational simulations and show that many of the outstanding problems in swimming energetics can be explained by considering the coupling between hydrodynamics and muscle contraction characteristics, as well as the trade-offs between the conflicting performance goals of sustained swimming speed U and cost of transport COT. Our results lead to three main conclusions: (1) in contrast to previous hypotheses, achieving optimal values of U and COT is independent of producing maximal power or efficiency; (2) muscle efficiency in swimming, in contrast to that in flying or running, decreases with increasing body size, consistent with muscle contraction characteristics; (3) the long-standing problem of two disparate patterns of longitudinal power output distributions in swimming fish can be reconciled by relating the two patterns to U-optimal or COT-optimal swimmers, respectively. We also provide further evidence that the use of tendons in caudal regions is beneficial from an energetic perspective. Our conclusions explain and unify many existing observations and are supported by computational data covering nine orders of magnitude in body size.

Tuna robotics: A high-frequency experimental platform exploring the performance space of swimming fishes

Tuna and related scombrid fishes are high-performance swimmers that often operate at high frequencies, especially during behaviors such as escaping from predators or catching prey. This contrasts with most fish-like robotic systems that typically operate at low frequencies (< 2 hertz). To explore the high-frequency fish swimming performance space, we designed and tested a new platform based on yellowfin tuna (Thunnus albacares) and Atlantic mackerel (Scomber scombrus). Body kinematics, speed, and power were measured at increasing tail beat frequencies to quantify swimming performance and to study flow fields generated by the tail. Experimental analyses of freely swimming tuna and mackerel allow comparison with the tuna-like robotic system. The Tunabot (255 millimeters long) can achieve a maximum tail beat frequency of 15 hertz, which corresponds to a swimming speed of 4.0 body lengths per second. Comparison of midline kinematics between scombrid fish and the Tunabot shows good agreement over a wide range of frequencies, with the biggest discrepancy occurring at the caudal fin, primarily due to the rigid propulsor used in the robotic model. As frequency increases, cost of transport (COT) follows a fish-like U-shaped response with a minimum at ~1.6 body lengths per second. The Tunabot has a range of ~9.1 kilometers if it swims at 0.4 meter per second or ~4.2 kilometers at 1.0 meter per second, assuming a 10-watt-hour battery pack. These results highlight the capabilities of high-frequency biological swimming and lay the foundation to explore a fish-like performance space for bio-inspired underwater vehicles.

3D computational models explain muscle activation patterns and energetic functions of internal structures in fish swimming

How muscles are used is a key to understanding the internal driving of fish swimming. However, the underlying mechanisms of some features of the muscle activation patterns and their differential appearance in different species are still obscure. In this study, we explain the muscle activation patterns by using 3D computational fluid dynamics models coupled to the motion of fish with prescribed deformation and examining the torque and power required along the fish body with two primary swimming modes. We find that the torque required by the hydrodynamic forces and body inertia exhibits a wave pattern that travels faster than the curvature wave in both anguilliform and carangiform swimmers, which can explain the traveling wave speeds of the muscle activations. Notably, intermittent negative power (i.e., power delivered by the fluid to the body) on the posterior part, along with a timely transfer of torque and energy by tendons, explains the decrease in the duration of muscle activation towards the tail. The torque contribution from the body elasticity further clarifies the wave speed increase or the reverse of the wave direction of the muscle activation on the posterior part of a carangiform swimmer. For anguilliform swimmers, the absence of the aforementioned changes in the muscle activation on the posterior part is consistent with our torque prediction and the absence of long tendons from experimental observations. These results provide novel insights into the functions of muscles and tendons as an integral part of the internal driving system, especially from an energy perspective, and they highlight the differences in the internal driving systems between the two primary swimming modes.

For undulatory swimming, fish form posteriorly traveling waves of body bending by activating their muscles sequentially along the body. However, experimental observations have shown that the muscle activation wave does not simply match the bending wave. Researchers have previously computed the torque required for muscles along the body based on classic hydrodynamic theories and explained the higher wave speed of the muscle activation compared to the curvature wave. However, the origins of other features of the muscle activation pattern and their variation among different species are still obscure after decades of research. In this study, we use 3D computational fluid dynamics models to compute the spatiotemporal distributions of both the torque and power required for eel-like and mackerel-like swimming. By examining both the torque and power patterns and considering the energy transfer, storage, and release by tendons and body viscoelasticity, we can explain not only the features and variations in the muscle activation patterns as observed from fish experiments but also how tendons and body elasticity save energy. We provide a mechanical picture in which the body shape, body movement, muscles, tendons, and body elasticity of a mackerel (or similar) orchestrate to make swimming efficient.

Immersed Methods for Fluid-Structure Interaction

Fluid-structure interaction is ubiquitous in nature and occurs at all biological scales. Immersed methods provide mathematical and computational frameworks for modeling fluid-structure systems. These methods, which typically use an Eulerian description of the fluid and a Lagrangian description of the structure, can treat thin immersed boundaries and volumetric bodies, and they can model structures that are flexible or rigid or that move with prescribed deformational kinematics. Immersed formulations do not require body-fitted discretizations and thereby avoid the frequent grid regeneration that can otherwise be required for models involving large deformations and displacements. This article reviews immersed methods for both elastic structures and structures with prescribed kinematics. It considers formulations using integral operators to connect the Eulerian and Lagrangian frames and methods that directly apply jump conditions along fluid-structure interfaces. Benchmark problems demonstrate the effectiveness of these methods, and selected applications at Reynolds numbers up to approximately 20,000 highlight their impact in biological and biomedical modeling and simulation.

Lift and Drag Acting on the Shell of the American Horseshoe Crab (Limulus polyphemus)

The intertidal zone is a turbulent landscape where organisms face numerous mechanical challenges from powerful waves. A model for understanding the solutions to these physical problems, the American horseshoe crab (Limulus polyphemus), is a marine arthropod that mates in the intertidal zone, where it must contend with strong ambient flows to maintain its orientation during locomotion and reproduction. Possible strategies to maintain position include either negative lift generation or the minimization of positive lift in flow. To quantify flow over the shell and the forces generated, we laser-scanned the 3D shape of a horseshoe crab, and the resulting digital reconstruction was used to 3D-print a physical model. We then recorded the movement of tracking particles around the shell model with high-speed video and analyzed the time-lapse series using particle image velocimetry (PIV). The velocity vector fields from PIV were used to validate numerical simulations performed with the immersed boundary (IB) method. IB simulations allowed us to resolve the forces acting on the shell, as well as the local three-dimensional flow velocities and pressures. Both IB simulations and PIV analysis of vorticity and velocity at a flow speed of 13 cm/s show negative lift for negative and zero angles of attack, and positive lift for positive angles of attack in a free-stream environment. In shear flow simulations, we found near-zero lift for all orientations tested. Because horseshoe crabs are likely to be found primarily at near-zero angles of attack, we suggest that this negative lift helps maintain the orientation of the crab during locomotion and mating. This study provides a preliminary foundation for assessing the relationship between documented morphological variation and potential environmental variation for distinct populations of horseshoe crabs along the Atlantic Coast. It also motivates future studies which could consider the stability of the horseshoe crab in unsteady, oscillating flows.

Red muscle activity in bluegill sunfish Lepomis macrochirus during forward accelerations

Fishes generate force to swim by activating muscles on either side of their flexible bodies. To accelerate, they must produce higher muscle forces, which leads to higher reaction forces back on their bodies from the environment. If their bodies are too flexible, the forces during acceleration could not be transmitted effectively to the environment, but fish can potentially use their muscles to increase the effective stiffness of their body. Here, we quantified red muscle activity during acceleration and steady swimming, looking for patterns that would be consistent with the hypothesis of body stiffening. We used high-speed video, electromyographic recordings, and a new digital inertial measurement unit to quantify body kinematics, red muscle activity, and 3D orientation and centre of mass acceleration during forward accelerations and steady swimming over several speeds. During acceleration, fish co-activated anterior muscle on the left and right side, and activated all muscle sooner and kept it active for a larger fraction of the tail beat cycle. These activity patterns are both known to increase effective stiffness for muscle tissue in vitro, which is consistent with our hypothesis that fish use their red muscle to stiffen their bodies during acceleration. We suggest that during impulsive movements, flexible organisms like fishes can use their muscles not only to generate propulsive power but to tune the effective mechanical properties of their bodies, increasing performance during rapid movements and maintaining flexibility for slow, steady movements.