The role of mechanical resonance in the neural control of swimming in fishes

Mapping spatial patterns to energetic benefits in groups of flow-coupled swimmers

The coordinated motion of animal groups through fluids is thought to reduce the cost of locomotion to individuals in the group. However, the connection between the spatial patterns observed in collectively moving animals and the energetic benefits at each position within the group remains unclear. To address this knowledge gap, we study the spontaneous emergence of cohesive formations in groups of fish, modeled as flapping foils, all heading in the same direction. We show in pairwise formations and with increasing group size that (1) in side-by-side arrangements, the reciprocal nature of flow coupling results in an equal distribution of energy requirements among all members, with reduction in cost of locomotion for swimmers flapping inphase but an increase in cost for swimmers flapping antiphase, and (2) in inline arrangements, flow coupling is non-reciprocal for all flapping phase, with energetic savings in favor of trailing swimmers, but only up to a finite number of swimmers, beyond which school cohesion and energetic benefits are lost at once. We explain these findings mechanistically and we provide efficient diagnostic tools for identifying locations in the wake of single and multiple swimmers that offer opportunities for hydrodynamic benefits to aspiring followers. Our results imply a connection between the resources generated by flow physics and social traits that influence greedy and cooperative group behavior.

Benchmarking the Immersed Boundary Method for Viscoelastic Flows

We present and analyze a series of benchmark tests regarding the application of the immersed boundary (IB) method to viscoelastic flows through and around non-trivial, stationary geometries. The IB method is widely used to simulate biological fluid dynamics and other modeling scenarios in which a structure is immersed in a fluid. Although the IB method has been most commonly used to model systems involving viscous incompressible fluids, it also can be applied to visoelastic fluids, and has enabled the study of a wide variety of dynamical problems including the settling of vesicles and the swimming of elastic filaments in fluids modeled by the Oldroyd-B constitutive equation. In the viscoelastic context, however, relatively little work has explored the accuracy or convergence properties of this numerical scheme. Herein, we present benchmarking results for an IB solver applied to viscoelastic flows in and around non-trivial geometries using either the idealized Oldroyd-B constitutive model or the more physcially realistic, polymer-entanglementbased Rolie-Poly constitutive equations. We use two-dimensional numerical test cases along with results from rheology experiments to benchmark the IB method and compare it to more complex finite element and finite volume viscoelastic flow solvers. Additionally, we analyze different choices of regularized delta function and relative Lagrangian grid spacings which allow us to identify and recommend the key choices of these numerical parameters depending on the present flow regime.

Pose analysis in free-swimming adult zebrafish, Danio rerio : "fishy" origins of movement design

Movement requires maneuvers that generate thrust to either make turns or move the body forward in physical space. The computational space for perpetually controlling the relative position of every point on the body surface can be vast. We hypothesize the evolution of efficient design for movement that minimizes active (neural) control by leveraging the passive (reactive) forces between the body and the surrounding medium at play. To test our hypothesis, we investigate the presence of stereotypical postures during free-swimming in adult zebrafish, Danio rerio . We perform markerless tracking using DeepLabCut, a deep learning pose estimation toolkit, to track geometric relationships between body parts. To identify putative clusters of postural configurations obtained from twelve freely behaving zebrafish, we use unsupervised multivariate time-series analysis (B-SOiD machine learning software). When applied to single individuals, this method reveals a best-fit for 36 to 50 clusters in contrast 86 clusters for data pooled from all 12 animals. The centroids of each cluster obtained over 14,000 sequential frames recorded for a single fish represent an apriori classification into relatively stable "target body postures" and inter-pose "transitional postures" that lead to and away from a target pose. We use multidimensional scaling of mean parameter values for each cluster to map cluster-centroids within two dimensions of postural space. From a post-priori visual analysis, we condense neighboring postural variants into 15 superclusters or core body configurations. We develop a nomenclature specifying the anteroposterior level/s (upper, mid and lower) and degree of bending. Our results suggest that constraining bends to mainly three levels in adult zebrafish preempts the neck, fore- and hindlimb design for maneuverability in land vertebrates.

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.

Linking neural circuits to the mechanics of animal behavior in Drosophila larval locomotion

The motions that make up animal behavior arise from the interplay between neural circuits and the mechanical parts of the body. Therefore, in order to comprehend the operational mechanisms governing behavior, it is essential to examine not only the underlying neural network but also the mechanical characteristics of the animal's body. The locomotor system of fly larvae serves as an ideal model for pursuing this integrative approach. By virtue of diverse investigation methods encompassing connectomics analysis and quantification of locomotion kinematics, research on larval locomotion has shed light on the underlying mechanisms of animal behavior. These studies have elucidated the roles of interneurons in coordinating muscle activities within and between segments, as well as the neural circuits responsible for exploration. This review aims to provide an overview of recent research on the neuromechanics of animal locomotion in fly larvae. We also briefly review interspecific diversity in fly larval locomotion and explore the latest advancements in soft robots inspired by larval locomotion. The integrative analysis of animal behavior using fly larvae could establish a practical framework for scrutinizing the behavior of other animal species.

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.

Role of the caudal peduncle in a fish-inspired robotic model: how changing stiffness and angle of attack affects swimming performance

In fish, the tail is a key element of propulsive anatomy that contributes to thrust during swimming. Fish possess the ability to alter tail stiffness, surface area and conformation. Specifically, the region at the base of the tail, the caudal peduncle, is proposed to be a key location of fish stiffness modulation during locomotion. Most previous analyses have focused on the overall body or tail stiffness, and not on the effects of changing stiffness specifically at the base of the tail in fish and robotic models. We used both computational fluid dynamics analysis and experimental measurements of propulsive forces in physical models with different peduncle stiffnesses to analyze the effect of altering stiffness on the tail angle of attack and propulsive force and efficiency. By changing the motion program input to the tail, we were able to alter the phase relationship between the front and back tail sections between 0° and 330°. Computational simulations showed that power consumption was nearly minimized and thrust production was nearly maximized at the kinematic pattern whereφ= 270°, the approximate phase lag observed in the experimental foils and in free swimming tuna. We observed reduced thrust and efficiency at high angles of attack, suggesting that the tail driven during these motion programs experiences stalling and loss of lift. However, there is no single peduncle stiffness that consistently maximizes performance, particularly in physical models. This result highlights the fact that the optimal caudal peduncle stiffness is highly context dependent. Therefore, incorporating the ability to control peduncle stiffness in future robotic models of fish propulsion promises to increase the ability of robots to approach the performance of fish.

Planktos: An Agent-Based Modeling Framework for Small Organism Movement and Dispersal in a Fluid Environment with Immersed Structures

Multiscale modeling of marine and aerial plankton has traditionally been difficult to address holistically due to the challenge of resolving individual locomotion dynamics while being carried with larger-scale flows. However, such problems are of paramount importance, e.g., dispersal of marine larval plankton is critical for the health of coral reefs, and aerial plankton (tiny arthropods) can be used as effective agricultural biocontrol agents. Here we introduce the open-source, agent-based modeling software Planktos targeted at 2D and 3D fluid environments in Python. Agents in this modeling framework are relatively tiny organisms in sufficiently low densities that their effect on the surrounding fluid motion can be considered negligible. This library can be used for scientific exploration and quantification of collective and emergent behavior, including interaction with immersed structures. In this paper, we detail the implementation and functionality of the library along with some illustrative examples. Functionality includes arbitrary agent behavior obeying either ordinary differential equations, stochastic differential equations, or coded movement algorithms, all under the influence of time-dependent fluid velocity fields generated by computational fluid dynamics, experiments, or analytical models in domains with static immersed mesh structures with sliding or sticky collisions. In addition, data visualization tools provide images or animations with kernel density estimation and velocity field analysis with respect to deterministic agent behavior via the finite-time Lyapunov exponent.

On the Lagrangian-Eulerian Coupling in the Immersed Finite Element/Difference Method

The immersed boundary (IB) method is a non-body conforming approach to fluid-structure interaction (FSI) that uses an Eulerian description of the momentum, viscosity, and incompressibility of a coupled fluid-structure system and a Lagrangian description of the deformations, stresses, and resultant forces of the immersed structure. Integral transforms with Dirac delta function kernels couple the Eulerian and Lagrangian variables, and in practice, discretizations of these integral transforms use regularized delta function kernels. Many different kernel functions have been proposed, but prior numerical work investigating the impact of the choice of kernel function on the accuracy of the methodology has often been limited to simplified test cases or Stokes flow conditions that may not reflect the method's performance in applications, particularly at intermediate-to-high Reynolds numbers, or under different loading conditions. This work systematically studies the effect of the choice of regularized delta function in several fluid-structure interaction benchmark tests using the immersed finite element/difference (IFED) method, which is an extension of the IB method that uses a finite element structural discretization combined with a Cartesian grid finite difference method for the incompressible Navier-Stokes equations. Whereas the conventional IB method spreads forces from the nodes of the structural mesh and interpolates velocities to those nodes, the IFED formulation evaluates the regularized delta function on a collection of interaction points that can be chosen to be denser than the nodes of the Lagrangian mesh. This opens the possibility of using structural discretizations with wide node spacings that would produce gaps in the Eulerian force in nodally coupled schemes (e.g., if the node spacing is comparable to or broader than the support of the regularized delta functions). Earlier work with this methodology suggested that such coarse structural meshes can yield improved accuracy for shear-dominated cases and, further, found that accuracy improves when the structural mesh spacing is increased. However, these results were limited to simple test cases that did not include substantial pressure loading on the structure. This study investigates the effect of varying the relative mesh widths of the Lagrangian and Eulerian discretizations in a broader range of tests. Our results indicate that kernels satisfying a commonly imposed even-odd condition require higher resolution to achieve similar accuracy as kernels that do not satisfy this condition. We also find that narrower kernels are more robust, in the sense that they yield results that are less sensitive to relative changes in the Eulerian and Lagrangian mesh spacings, and that structural meshes that are substantially coarser than the Cartesian grid can yield high accuracy for shear-dominated cases but not for cases with large normal forces. We verify our results in a large-scale FSI model of a bovine pericardial bioprosthetic heart valve in a pulse duplicator.

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.

Exploring the sensitivity in jellyfish locomotion under variations in scale, frequency, and duty cycle

Jellyfish have been called one of the most energy-efficient animals in the world due to the ease in which they move through their fluid environment, by product of their bell kinematics coupled with their morphological, muscular, material properties. We investigated jellyfish locomotion by conducting in silico comparative studies and explored swimming performance across different fluid scales (i.e., Reynolds Number), bell contraction frequencies, and contraction phase kinematics (duty cycle) for a jellyfish with a fineness ratio of 1 (ratio of bell height to bell diameter). To study these relationships, an open source implementation of the immersed boundary method was used (IB2d) to solve the fully coupled fluid-structure interaction problem of a flexible jellyfish bell in a viscous fluid. Thorough 2D parameter subspace explorations illustrated optimal parameter combinations in which give rise to enhanced swimming performance. All performance metrics indicated a higher sensitivity to bell actuation frequency than fluid scale or duty cycle, via Sobol sensitivity analysis, on a higher performance parameter subspace. Moreover, Pareto-like fronts were identified in the overall performance space involving the cost of transport and forward swimming speed. Patterns emerged within these performance spaces when highlighting different parameter regions, which complemented the global sensitivity results. Lastly, an open source computational model for jellyfish locomotion is offered to the science community that can be used as a starting place for future numerical experimentation.

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.

Elastohydrodynamic Scaling Law for Heart Rates

Animal hearts are soft shells that actively pump blood to oxygenate tissues. Here, we propose an allometric scaling law for the heart rate based on the idea of elastohydrodynamic resonance of a fluid-loaded soft active elastic shell that buckles and contracts axially when twisted periodically. We show that this picture is consistent with numerical simulations of soft cylindrical shells that twist-buckle while pumping a viscous fluid, yielding optimum ejection fractions of 35%-40% when driven resonantly. Our scaling law is consistent with experimental measurements of heart rates over 2 orders of magnitude, and provides a mechanistic basis for how metabolism scales with organism size. In addition to providing a physical rationale for the heart rate and metabolism of an organism, our results suggest a simple design principle for soft fluidic pumps.

Collective Pulsing in Xeniid Corals: Part II-Using Computational Fluid Dynamics to Determine if There are Benefits to Coordinated Pulsing

Coordinated movements have been shown to enhance the speed or efficiency of swimming, flying, and pumping in many organisms. Coordinated pulsing has not been observed in many cnidarians (jellyfish, anemones, corals), as is the case for the xeniid corals considered in our corresponding paper. This observation opens the question as to whether xeniid corals, and cnidarians in general, do not coordinate their pulsing behavior for lack of a hydrodynamic advantage or for other reasons. For example, a diffuse nervous system with lack of substantial sensory input may not be capable of such coordination. Similarly, grouping may serve a defensive role rather than a fluid dynamic role. In this paper, the immersed boundary method is used to quantify the volumetric flux of fluid generated by an individual xeniid coral polyp in comparison with a pair of polyps. Both the distances between the polyps and the phase difference between each polyp are considered. More specifically, the fully coupled fluid-structure interaction problem of a coral polyp driving fluid flow is solved using a hybrid version of the immersed boundary method where the Navier-Stokes equations are solved using a finite differences and the elasticity equations describing the coral are solved using finite elements. We explore three possible hypotheses: (1) pulsing in pairs increases upward flow above the polyps and is thus beneficial, (2) these benefits vary with the polyps' pulsing phase difference, and (3) these benefits vary with the distance between the polyps. We find that there is no substantial hydrodynamic advantage to pulsing in a pair for any phase difference. The volumetric flux of fluid generated by each coral also decreases as the distance between polyps is decreased. This surprising result is consistent with measurements taken from another cnidarian with similar behavior, the upside down jellyfish, in which each medusa drives less flow when in a group.

An Immersed Interface Method for Discrete Surfaces

Fluid-structure systems occur in a range of scientific and engineering applications. The immersed boundary (IB) method is a widely recognized and effective modeling paradigm for simulating fluid-structure interaction (FSI) in such systems, but a difficulty of the IB formulation of these problems is that the pressure and viscous stress are generally discontinuous at fluid-solid interfaces. The conventional IB method regularizes these discontinuities, which typically yields low-order accuracy at these interfaces. The immersed interface method (IIM) is an IB-like approach to FSI that sharply imposes stress jump conditions, enabling higher-order accuracy, but prior applications of the IIM have been largely restricted to numerical methods that rely on smooth representations of the interface geometry. This paper introduces an immersed interface formulation that uses only a C 0 representation of the immersed interface, such as those provided by standard nodal Lagrangian finite element methods. Verification examples for models with prescribed interface motion demonstrate that the method sharply resolves stress discontinuities along immersed boundaries while avoiding the need for analytic information about the interface geometry. Our results also demonstrate that only the lowest-order jump conditions for the pressure and velocity gradient are required to realize global second-order accuracy. Specifically, we demonstrate second-order global convergence rates along with nearly second-order local convergence in the Eulerian velocity field, and between first- and second-order global convergence rates along with approximately first-order local convergence for the Eulerian pressure field. We also demonstrate approximately second-order local convergence in the interfacial displacement and velocity along with first-order local convergence in the fluid traction along the interface. As a demonstration of the method's ability to tackle more complex geometries, the present approach is also used to simulate flow in a patient-averaged anatomical model of the inferior vena cava, which is the large vein that carries deoxygenated blood from the lower extremities back to the heart. Comparisons of the general hemodynamics and wall shear stress obtained by the present IIM and a body-fitted discretization approach show that the present method yields results that are in good agreement with those obtained by the body-fitted approach.

A novel mechanism of mixing by pulsing corals

The dynamic pulsation of xeniid corals is one of the most fascinating phenomena observed in coral reefs. We quantify for the first time the flow near the tentacles of these soft corals, the active pulsations of which are thought to enhance their symbionts' photosynthetic rates by up to an order of magnitude. These polyps are approximately 1 cm in diameter and pulse at frequencies between approximately 0.5 and 1 Hz. As a result, the frequency-based Reynolds number calculated using the tentacle length and pulse frequency is on the order of 10 and rapidly decays as with distance from the polyp. This introduces the question of how these corals minimize the reversibility of the flow and bring in new volumes of fluid during each pulse. We estimate the Péclet number of the bulk flow generated by the coral as being on the order of 100-1000 whereas the flow between the bristles of the tentacles is on the order of 10. This illustrates the importance of advective transport in removing oxygen waste. Flow measurements using particle image velocimetry reveal that the individual polyps generate a jet of water with positive vertical velocities that do not go below 0.1 cm s^-1 and with average volumetric flow rates of approximately 0.71 cm³ s^-1 Our results show that there is nearly continual flow in the radial direction towards the polyp with only approximately 3.3% back flow. 3D numerical simulations uncover a region of slow mixing between the tentacles during expansion. We estimate that the average flow that moves through the bristles of the tentacles is approximately 0.03 cm s^-1 The combination of nearly continual flow towards the polyp, slow mixing between the bristles, and the subsequent ejection of this fluid volume into an upward jet ensures the polyp continually samples new water with sufficient time for exchange to occur.

Closed-Loop Control of Active Sensing Movements Regulates Sensory Slip

Active sensing involves the production of motor signals for the purpose of acquiring sensory information [1-3]. The most common form of active sensing, found across animal taxa and behaviors, involves the generation of movements-e.g., whisking [4-6], touching [7, 8], sniffing [9, 10], and eye movements [11]. Active sensing movements profoundly affect the information carried by sensory feedback pathways [12-15] and are modulated by both top-down goals (e.g., measuring weight versus texture [1, 16]) and bottom-up stimuli (e.g., lights on or off [12]), but it remains unclear whether and how these movements are controlled in relation to the ongoing feedback they generate. To investigate the control of movements for active sensing, we created an experimental apparatus for freely swimming weakly electric fish, Eigenmannia virescens, that modulates the gain of reafferent feedback by adjusting the position of a refuge based on real-time videographic measurements of fish position. We discovered that fish robustly regulate sensory slip via closed-loop control of active sensing movements. Specifically, as fish performed the task of maintaining position inside the refuge [17-22], they dramatically up- or downregulated fore-aft active sensing movements in relation to a 4-fold change of experimentally modulated reafferent gain. These changes in swimming movements served to maintain a constant magnitude of sensory slip. The magnitude of sensory slip depended on the presence or absence of visual cues. These results indicate that fish use two controllers: one that controls the acquisition of information by regulating feedback from active sensing movements and another that maintains position in the refuge, a control structure that may be ubiquitous in animals [23, 24].

The role of curvature feedback in the energetics and dynamics of lamprey swimming: A closed-loop model

Like other animals, lampreys have a central pattern generator (CPG) circuit that activates muscles for locomotion and also adjusts the activity to respond to sensory inputs from the environment. Such a feedback system is crucial for responding appropriately to unexpected perturbations, but it is also active during normal unperturbed steady swimming and influences the baseline swimming pattern. In this study, we investigate different functional forms of body curvature-based sensory feedback and evaluate their effects on steady swimming energetics and kinematics, since little is known experimentally about the functional form of curvature feedback. The distributed CPG is modeled as chains of coupled oscillators. Pairs of phase oscillators represent the left and right sides of segments along the lamprey body. These activate muscles that flex the body and move the lamprey through a fluid environment, which is simulated using a full Navier-Stokes model. The emergent curvature of the body then serves as an input to the CPG oscillators, closing the loop. We consider two forms of feedback, each consistent with experimental results on lamprey proprioceptive sensory receptors. The first, referred to as directional feedback, excites or inhibits the oscillators on the same side, depending on the sign of a chosen gain parameter, and has the opposite effect on oscillators on the opposite side. We find that directional feedback does not affect beat frequency, but does change the duration of muscle activity. The second feedback model, referred to as magnitude feedback, provides a symmetric excitatory or inhibitory effect to oscillators on both sides. This model tends to increase beat frequency and reduces the energetic cost to the lamprey when the gain is high and positive. With both types of feedback, the body curvature has a similar magnitude. Thus, these results indicate that the same magnitude of curvature-based feedback on the CPG with different functional forms can cause distinct differences in swimming performance.

When animals move, they receive sensory inputs, which in turn are used to modulate the movement. Relatively little is known about how these inputs affect performance during steady locomotion. Using a computational model of a swimming lamprey, we investigated two different types of feedback, both consistent with experimental data. Both have strong, but different, effects on swimming speed and energy consumption, suggesting that sensory feedback is crucial not just for responding to perturbations, but also for high performance steady locomotion.

Soft Biomimetic Fish Robot Made of Dielectric Elastomer Actuators

This article presents the design, fabrication, and characterization of a soft biomimetic robotic fish based on dielectric elastomer actuators (DEAs) that swims by body and/or caudal fin (BCF) propulsion. BCF is a promising locomotion mechanism that potentially offers swimming at higher speeds and acceleration rates, and efficient locomotion. The robot consists of laminated silicone layers wherein two DEAs are used in an antagonistic configuration, generating undulating fish-like motion. The design of the robot is guided by a mathematical model based on the Euler–Bernoulli beam theory and takes account of the nonuniform geometry of the robot and of the hydrodynamic effect of water. The modeling results were compared with the experimental results obtained from the fish robot with a total length of 150 mm, a thickness of 0.75 mm, and weight of 4.4 g. We observed that the frequency peaks in the measured thrust force produced by the robot are similar to the natural frequencies computed by the model. The peak swimming speed of the robot was 37.2 mm/s (0.25 body length/s) at 0.75 Hz. We also observed that the modal shape of the robot at this frequency corresponds to the first natural mode. The swimming of the robot resembles real fish and displays a Strouhal number very close to those of living fish. These results suggest the high potential of DEA-based underwater robots relying on BCF propulsion, and applicability of our design and fabrication methods.

Automatic control: the vertebral column of dogfish sharks behaves as a continuously variable transmission with smoothly shifting functions

During swimming in dogfish sharks, Squalus acanthias, both the intervertebral joints and the vertebral centra undergo significant strain. To investigate this system, unique among vertebrates, we cyclically bent isolated segments of 10 vertebrae and nine joints. For the first time in the biomechanics of fish vertebral columns, we simultaneously characterized non-linear elasticity and viscosity throughout the bending oscillation, extending recently proposed techniques for large-amplitude oscillatory shear (LAOS) characterization to large-amplitude oscillatory bending (LAOB). The vertebral column segments behave as non-linear viscoelastic springs. Elastic properties dominate for all frequencies and curvatures tested, increasing as either variable increases. Non-linearities within a bending cycle are most in evidence at the highest frequency, 2.0 Hz, and curvature, 5 m−1. Viscous bending properties are greatest at low frequencies and high curvatures, with non-linear effects occurring at all frequencies and curvatures. The range of mechanical behaviors includes that of springs and brakes, with smooth transitions between them that allow for continuously variable power transmission by the vertebral column to assist in the mechanics of undulatory propulsion.

Batoid fish locomotion: effects of speed on pectoral fin deformation in the little skate Leucoraja erinacea

Most batoid fishes have a unique swimming mode in which thrust is generated by either oscillating or undulating expanded pectoral fins that form a disc. Only one previous study of the freshwater stingray has quantified three-dimensional motions of the wing, and no comparable data are available for marine batoid species that may differ considerably in their mode of locomotion. Here we investigate three-dimensional kinematics of the pectoral wing of the little skate, Leucoraja erinacea, swimming steadily at two speeds (1 and 2 body lengths per second, BL×s−1). We measured the motion of nine points in three dimensions during wing oscillation and determined that there are significant differences in movement amplitude among wing locations, as well as significant differences as speed increases in body angle, wing beat frequency, and speed of the traveling wave on the wing. In addition, we analyzed differences in wing curvature with swimming speed. At 1 BL×s−1, the pectoral wing is convex in shape during the downstroke along the medio-lateral fin midline, but at 2 BL×s−1 the pectoral fin at this location cups into the flow indicating active curvature control and fin stiffening. Wing kinematics of the little skate differed considerably from previous work on the freshwater stingray, which does not show active cupping of the whole fin on the downstroke.

Axial systems and their actuation: new twists on the ancient body of craniates

Craniate animals--vertebrates and their jawless sister taxa--have evolved a body axis with powerful muscles, a distributed nervous system to control those muscles, and an endoskeleton that starts at the head and ends at the caudal fin. The body axis undulates, bends, twists, or holds firm, depending on the behavior. In this introduction to the special issue on axial systems and their actuation, we provide an overview of the latest research on how the body axis functions, develops, and evolves. Based on this research, we hypothesize that the body axis of craniates has three primary, post-cranial modules: precaudal, caudal, and tail. The term "module" means a portion of the body axis that functions, develops, and evolves in relative independence from other modules; "relative independence" means that structures and processes within a module are more tightly correlated in function, development, and behavior than the same processes are among modules.