Pleobot: a modular robotic solution for metachronal swimming

Metachronal propulsion is widespread in aquatic swarming organisms to achieve performance and maneuverability at intermediate Reynolds numbers. Studying only live organisms limits our understanding of the mechanisms driving these abilities. Thus, we present the design, manufacture, and validation of the Pleobot—a unique krill-inspired robotic swimming appendage constituting the first platform to study metachronal propulsion comprehensively. We combine a multi-link 3D printed mechanism with active and passive actuation of the joints to generate natural kinematics. Using force and fluid flow measurements in parallel with biological data, we show the link between the flow around the appendage and thrust. Further, we provide the first account of a leading-edge suction effect contributing to lift during the power stroke. The repeatability and modularity of the Pleobot enable the independent manipulation of particular motions and traits to test hypotheses central to understanding the relationship between form and function. Lastly, we outline future directions for the Pleobot, including adapting morphological features. We foresee a broad appeal to a wide array of scientific disciplines, from fundamental studies in ecology, biology, and engineering, to developing new bio-inspired platforms for studying oceans across the solar system.


Introduction
Nature-inspired robotics have become valuable tools in biology and ecology as they allow better-controlled comparative studies than traditional organismal methods 1 .Independent manipulation of a particular trait or behavior of interest enables testing hypotheses central to exploring the underlying mechanisms driving maneuverability and performance.While several test platforms have shed light on important aspects of biological underwater propulsion [2][3][4][5] , they have also produced practical engineering solutions for unmanned underwater vehicles (UUV).Previous robotic solutions for engineering problems, including tunas 6 , scallop 7 , and dolphins 8 demonstrate improved propulsive efficiency and optimal operation in a range of environments.However, because of experimental or technical challenges, most advances have been for relatively high Reynolds numbers (Re).Consequently, despite the tremendous diversity of species inhabiting the moderate realm, we still need a comprehensive understanding of aquatic propulsion and its derived engineering applications at low-to-intermediate Re.
Metachronal swimming is a ubiquitous swimming mode among many of the most abundant invertebrate species inhabiting the transitional flow regime (Re < 1000), such as copepods and krill.This propulsive gait is characterized by the sequential beating of several closely spaced swimming appendages, phase-shifted in time, in a tail-to-head traveling wave 9 .Metachrony enables large-scale behaviors in krill, particularly the ability to perform diel vertical migrations (DVM) from the sea surface to oxygen minimum zones located up to 1 km at depth 10,11 .Efficient swimming is achieved by leveraging drag during properly tuned power-and-recovery strokes 9,12 .During the power stroke, the swimming appendages (pleopods) move opposite the swimming direction while maximizing their surface area.This generates drag in each pleopod, producing thrust to propel the organism.Conversely, during the recovery stroke, the profile area decreases, thus reducing the drag on the appendages and creating a net thrust force sufficient to overcome the drag on the body.Previous works showed drag-based metachronal swimming is more effective than lift-based propulsion at intermediate Re for accelerating, braking, and turning because it generates significant thrust over short periods 9,[13][14][15] .
Fluid flow experiments on live krill have established the foundations of our understanding of the effect of fluid-structure interactions on the far-field flow 16,17 .Murphy et al. characterized the swimming kinematics of live krill and linked the success of the propulsion system to appendage morphology, stroke kinematics, and the resulting hydrodynamic effects 16 .Tip vortices forming on the pleopods of krill 17 and mantis shrimp 18 were shown to be central in the production of thrust 19,20 .Kim and Gharib extended these results using idealized pleopod shapes with simplified kinematics.They demonstrated that the area over the surface of the appendage enclosed by the tip vortex, specifically, modulates thrust production 19 .However, challenges in characterizing the near-field flow and measuring the forces generated during swimming have hindered efforts to establish the relationship between fluid dynamics and force distribution (i.e., thrust and lift) in intermediate Re metachronal swimmers.Simplified robotic models and numerical simulations have complemented these efforts, shedding light on the role of varying Re, phase lag, and appendage spacing on the hydrodynamics 18,[21][22][23][24][25] .For instance, Ford et al. showed that at low Re, the individual jets produced by adjacent pleopods do not interact due to viscous dissipation, but at Re around 800, they form a near-steady jet 22 .The outcome of this difference is an increase in vertical and horizontal momentum that enables krill to generate thrust and lift forces required for locomotion and hovering.Phase lag also emerges as an important factor contributing to metachronal swimming efficiency, as it yields near-maximal efficiency and thrust, and achieves higher average body velocity when compared to synchronous motions 21,22,24,26 .Likewise, appendage spacing is a critical morphological factor linked to stroke kinematics.Previous works demonstrated that greater swimming performance is achieved when the ratio of appendage spacing to appendage length is less than one 27 .This was hypothesized to be the direct consequence of inter-pleopod fluid flow interactions enhancing the vortex strength and circulation during the power stroke 18 .
While the insights gathered through simplified systems are invaluable, we still need a unifying theory of drag-based metachronal swimming to explain the relationship between appendage kinematics, fluid flows, and force production.For instance, even the production of lift needed for hovering still needs to be explored.Adding complexity to better match the natural biological and hydrodynamic phenomena is necessary to elucidate the effects of morphological and kinematic characteristics on the near-and far-field flow.We should then adopt an integrative framework allowing for the comprehensive study of the link between locomotor kinematics, form, and function.
Motivated by a robotics-inspired biology approach to address this problem, we present the first fully articulated, multi-link, metachronal robotic appendage reproducing the swimming kinematics of krill.Our novel design achieves active control of both the proximal and distal appendages and the passive out-of-plane actuation of the two rami of the distal appendages through the interaction between the pleopod and the surrounding fluid.In addition, we implemented the main morphological features of krill pleopods through the assembly of 3D-printed modules.Using force and fluid flow measurements, we quantify the effects of particular flow features on thrust and lift production and disentangle their role during the power and recovery phases of a stroke.We use biological data to validate our model as a suitable test platform.In this study, we leverage the Pleobot to study the near-field flow of a metachronal pleopod for the first time and provide evidence for a previously unreported leading-edge suction effect that contributes to lift during the power stroke.The repeatability and modularity of the Pleobot enables a range of configurations and motion programs to investigate metachrony across Re, taxa, and morphological attributes.As such, we anticipate the Pleobot will broadly appeal to a wide array of scientific disciplines, such as ecology, biology, and engineering.

Pleobot kinematics
The swimming kinematics of the Pleobot are characterized by three main angles: 1) α, the angle between the body axis and the protopodite; 2) β , the angle between the protopodite and the biramous distal appendage (including the endopodite and exopodite); and 3) γ, the angle between the endopodite and exopodite (Fig. 1A-D).The cupping angle, ζ , is also implemented in the Pleobot as seen in Euphausia superba (Fig. 1E,F).The resulting krill-inspired model (Fig. 1G) uses a train of gears to actuate both α and β kinematics (Fig. 1B).γ is actuated passively due to hydrodynamic forces during the power and recovery strokes (Fig. 1D), and the ζ angle of the exopodite is fixed (Fig. 1F).
The kinematics of the Pleobot were validated against organismal data 16 by tracking the motion of the robotic system over several consecutive cycles (n = 5).The cycle-averaged α and β angles over time were in line with kinematics data for E. superba (Fig. 2A,B) 16 .The slight deviation of the β angle from krill data was associated with the tolerance between gears, which is directly proportional to the number of units in the train.The motion of the passively actuated exopodite ranged from 14 • to 95 • , with abduction occurring at the end of the return stroke and adduction coinciding with the start of the recovery strokes (Fig. 1).Comparative values over a complete beat cycle have not been reported for E. superba, but Murphy and collaborators 16 described γ oscillating between 0 • and 77 • .The kinematic analysis of the Pleobot demonstrates it can perform the prescribed motion program with minimal error and agrees with the values reported for live krill.As such, Pleobot is an accurate and highly repeatable robotic analog of the swimming appendages of live metachronal organisms 16 .

Flow visualization
PIV measurements were performed in four different planes to evaluate the flow field (Fig. 3A).At the beginning of the power stroke (t/T, the non-dimensional stroke time, is 0), a vortex forms at the tip of the endopodite accompanied by a bound vortex uniformly distributed along most of the length of this segment (Fig. 3E).This pattern agrees with biological data acquired using live ghost shrimp (Fig. 3F).Meanwhile, the exopodite abducts as α increases and β decreases (Fig. 2).During the first half of the power stroke (before reaching t/T = 0.25), the exopodite maintains a low angle of attack (AoA, measured relative to the body axis, which is horizontal, in a clockwise direction), promoting the formation of an attached leading-edge vortex (LEV) along the anterior face (in the swimming direction) of the exopodite (Fig. 3B).The LEV remains attached until the exopodite reaches AoA = 73 o ± 1.9 o at t/T = 0.3.As α angle increases, the AoA of the exopodite also increases to become vertical relative to the swimming direction (AoA = 90 o ) during the second half of the power stroke (Fig. 4I).Given that the pleopod decelerates during this phase, the LEV separates from the exopodite (Figs 3C, 4I) and the bound vortex of the endopodite sheds at the tip ( Figs 3E,F, 4E).The combined effect of the flow entrained by the two counter-rotating vortices results in the production of a large posteriorly oriented downward jet.
Upon initiating the recovery stroke, the exopodite adducts and overlaps with the endopodite to effectively reduce the overall profile area of the pleopod subject to the flow (Figs 2C,E, 3D, 4C).Contrary to the power stroke, this configuration generates a pair of counter-rotating side edge vortices that induces spanwise flow in the swimming direction (Figs 3D, 4F,G).

Hydrodynamic forces
Force measurements were performed to evaluate the contribution of lift and thrust forces during steady forward swimming (Fig. 4).The lift (C L ) and thrust (C T ) coefficients were computed for a 20x-scaled endopodite-and-exopodite model.At the beginning of the power stroke (t/T = 0), there is a significant amount of thrust force, likely because of the incoming spanwise flow produced behind the pleopod at the end of the previous beat (Fig. 4G).During the power stroke, thrust and lift strongly correlate with the abduction of the exopodite.The forces on the endopodite decompose in the thrust and lift directions, and the LEV forming along the exopodite drives this association.Throughout the first third of the power stroke (t/T ∼ 0.16), the lift coefficient increases sharply to reach its global maximum (C L = 0.44) as the exopodite is abducting and produces an attached LEV (Fig. 4B,C,H).In contrast, while thrust is being produced, it remains mostly constant during this phase, and only increases significantly after lift peaks and decreases (Fig. 4A).This inverse relationship coincides with the AoA of the exopodite shifting to a more vertical orientation that causes the LEV to shed.It also results from the increased contribution of the endopodite to thrust, given the overall vertical orientation of the entire pleopod (see Fig. 4A,B).During the second half of the power stroke, thrust decreases sharply as the bound vortex of the endopodite moves toward the tip and forms a stopping vortex (Fig. 4E).Given the increase of the effective angle of the pleopod relative to the flow (angle Ψ, the overall angle between the body axis and the frontal face of the endopodite; see Fig. 4), the endopodite and the exopodite produce a downward force inducing minimum lift.The subsequent increase in lift force is likely due to the occurrence of spanwise flow entrained by the stopping vortex that pushes against the anterior face of the decelerating endopodite and exopodite (Fig. 4E).
During the first 2/3 of the recovery stroke, thrust is negative, suggesting this phase of a beat is mostly dominated by drag forces.The thrust coefficient is lowest around the mid-recovery stroke (Fig. 4A) when the pleopod accelerates.Thrust increases gradually to become positive, thus indicating some thrust is produced even during the recovery stroke.After the start of the recovery stroke, lift drops sharply to another minimum with decreasing ψ angle (Fig. 4B).Lift rises again at the end of the recovery stroke when the appendage decelerates.Increasing thrust and lift coincide with the occurrence of spanwise flow behind the overlapping endopodite and exopodite (Fig. 4G).

Discussion
In recent years, the growing need for maneuverable AUVs for underwater exploration has galvanized approaches inspired by nature [28][29][30][31][32] .In particular, the discovery of extraterrestrial oceans motivates the development of novel robotic platforms that will likely require the efficiency, versatility, and maneuverability of metachronal swimming.As such, the Pleobot lays down the foundations for the upcoming work on metachronal swimming to inform the design of underwater explorers.
Pleobot represents the first fully articulated robot to study metachronal swimming.The design constitutes a modular, 3D-printed platform incorporating the kinematic and morphological characteristics of krill -a marine organism observed to form large aggregations and migrate hundreds of meters into the ocean.While krill kinematics have been explored in great detail, a quantitative analysis of the fluid flow near individual appendages is needed to link kinematics to hydrodynamics and understand force production.Pleobot is the first system to disentangle the roles of kinematic parameters, morphology, and near-field hydrodynamics on thrust and lift generation at the scale of a single appendage.
Our unique study shows that Pleobot serves as an analog system to investigate a wide range of parameters and characteristics that collectively contribute to the efficiency and flow characteristics of metachronal propulsion.Although krill has been used as a model organism in biological, computational, and robotic studies, the focus has mostly been on the flow field of several appendages rather than the singular mechanisms for the generation of thrust and lift.Conclusions regarding fluid-structure interactions are thus solely possible by measuring the near field flow and by introducing complexities in the model that enable independent testing of hypotheses relevant to the ecology and behavior of the natural systems (i.e., maneuverability, migrations, fast swimming).
In our study, the analysis of force and PIV data allows us to evaluate the contribution of specific flow features to the generation of lift and thrust in Pleobot.At the beginning of the power stroke, the endopodite contributes to thrust and lift through the formation of tip and bound vortices.Similarly, the exopodite also contributes to the generation of both forces.Creation of lift during steady forward swimming is necessary for krill because they are negatively buoyant and must create lift to maintain their position in the water column 17 .As the stroke progresses, the exopodite becomes more vertically oriented and its contribution to lift decreases, while its contribution to thrust increases (Fig. 4A,B).This represents a trade-off between the generation of lift and thrust production.While the primary purpose of thrust in steady swimming is to propel the organism forward, lift is needed to guarantee the vertical position in the water column.Our data suggest that rather than being determined by behavioral processes, thrust and lift production are constrained to a relatively narrow range of kinematics and morphological criteria.Behavioral changes in the gait and beat frequency invariably modulate these forces.However, the dependency on the ψ angle of the passive abduction of the exopodite and its orientation in the flow demonstrates that a limited set of system requirements are fundamentally important for thrust and lift production in metachronal propulsion.
Notably, we reveal the existence of a LEV along the exopodite during the power stroke.Force measurements show that lift is enhanced by the formation of this attached LEV on the exopodite during the power stroke.Furthermore, this vortex delays the loss in lift forces slightly to maintain positive lift throughout part of the power stroke.This lift-based mechanism likely supplements the effects of the drag-based mechanism stemming from the pleopod pushing on water.These results in the near field illustrate the importance of the Pleobot, not only as a platform to engineer solutions for ocean exploration in the intermediate Reynolds number regime but also as a framework to understand nature.Force measurements on a biologicallyinspired appendage highlight the importance of understanding the vortex dynamics around the pleopod and reveal that both lift and thrust are generated simultaneously.
The incorporation of natural kinematics and morphological traits into the Pleobot provides promising results and incentives for increasing the complexity of our robotic model and elucidate additional mechanisms contributing to efficiency in krill and other metachronal swimmers.Markedly, the presence of hair-like setae around the endopodite and exopodite, and flexibility of the propulsors have the most potential for altering the flow field around the pleopods.Setae are defining features in metachronal organisms in the low-to-intermediate Re that increase the surface area in contact with water to enhance thrust 9,33,34 .Flexibility has been shown to be an important characteristic of propulsors in smoothing out thrust peaks to generate nearly constant thrust during the power stroke 19 .Pleobot provides a platform to evaluate the importance of these additional morphological features.Moreover, its modular design can be easily replicated and grouped in several units to investigate mechanisms such as constructive pleopod interactions (Fig. 5).For example, the proximity of several appendages beating at a phase lag promotes fluid flow interactions that likely contribute to enhanced thrust 18,35 .The Pleobot emerges as an ideal platform to quantify these effects.
Finally, while organismal investigations form the core of our knowledge about metachronal swimming 13,[16][17][18] , working with live animals provides only a partial understanding of the biomechanical and hydrodynamic mechanisms leading to efficiency and performance.Recent mechanical and robotics approaches allowing repeatable measurements highlight the valuable contribution of low-cost, controllable test platforms for biological studies 22,27,35 .Such a comparative approach will be central to developing a unifying theory of metachronal swimming.Manufacturing of the Pleobot is accessible to both institutions and individuals who want to experiment (see oliveira:zenodo:2022)). Its modular design allows for the quick prototyping of different appendage shapes and sizes.Furthermore, the control of both its proximal and distal segments facilitates the study of a wide range of systems, from small copepods, to krill, and larger species like mantis shrimp and lobsters.The Pleobot can thus be employed to facilitate the comprehensive characterization of metachronal, drag-based propulsion and help establish a unifying theory for this locomotor mechanism.

Modeling of krill appendage kinematics
To reproduce the kinematics of Euphausia superba, as reported by Murphy et al. 16 , we use values for α and β to find kinematic relationships for Pleobot.The Pleobot is dynamically scaled based on the Reynolds number.This dimensionless parameter is defined using the velocity and the length scale of the appendages and is approximately 600 for live krill: where L is the length of the pleopod, U tip the average velocity at its tip, and ν the kinematic viscosity of the fluid.U tip is calculated using the stroke amplitude θ , the frequency n, and the length L as 2θ nL.
The locomotor system is implemented by a multi-link mechanism using a transmission gear box actuating both the proximal and distal segments (corresponding to α and β ) with independently controlled servo motors.Angle γ represents out-of-plane motion, posing challenges for active actuation and is thus passively actuated via the hydrodynamic interaction between the structure (pleopods) and the fluid.

4/13
The kinematic relationships governing the positions of the gears are 36,37 : This equation gives the relationship between the driving and driven gears, ∆φ k+1 and ∆φ k−1 , respectively, and the link connecting both, ∆φ k (Supplementary figure 1), where N and r are the number of teeth and radius of the corresponding gear.The position of the pleopod, R pi , is governed for a given appendage i as (Supplementary figure 2,B): where The local reference frame vectors for the protopodite and exopodite are r 1i and r 2i .Similarly, R(θ 1i ) and R(θ 2i ) are the global rotation matrices and x 1i and x 2i the lengths (Supplementary figure 2,A).
Finally, the angular displacements of the gears are obtained using the relationships described in Equation 2: where ψ is the rotation of the gear along link 1 (endopodite), θ is the angle of the first link measured from the global reference frame (Supplementary figure 2,A), and r e is the gear ratio.Equations ( 3) and ( 5) allow for the solution of direct kinematics of the mechanism by using α i and β i , to calculate ψ 1i , which moves the second bar (endopodite)(Supplementary figure 2,C).The first bar (protopodite) is moved by angle α 1i (Supplementary figure 2,A).The metachronal trajectory at the end of pleopod 1, is given by R p1 = [x p1 , y p1 ] T .Manufacturing was completed by 3D printing, keeping small tolerances to reduce vibrations and loss of movement of the gear train.The supports house the servos and gears, and were designed to be an above-water structure.Bearings are used to reduce friction between gears and the axes.The transmission has an amplification of 2.5 to achieve the desired angular speed for the links.

Exopodite and endopodite
As E. superba propels forward, the endopodite and exopodite abduct and adduct to effectively change the profile area of the appendages to generate thrust and reduce drag.This motion is characterized by γ, the angle between the exopodite and the endopodite.In the Pleobot, this process is actuated on the horizontal plane, lateral to the protopodite.
In forward-swimming krill, exopodite abduction occurs at the beginning of the power stroke, and reaches a maximum of 77 •16 .This motion induces cupping of the appendages creating a V-shaped structure 17 , reminiscent of those observed in swimming fish that have been shown to produce greater thrust compared to flat fins (e.g., see 2 ).The exopodite and endopodite adduct during the recovery stroke.
Appendage cupping forms angle ζ between the endopodite and the exopodite.Photographic evidence was used to quantify ζ by measuring the angle between the midplanes of the endopodite and the exopodite (Fig. 1E) .The Pleobot was set to match the mean quantified value of 37 • .

Robotic design
Modular CAD designs were printed with polylactic acid (PLA) using a Prusa i3 MKS3+ 3D printer (Prusa Research, Prague, Czech Republic) for fast prototyping.Each pair of appendages is actuated by two servos (HS-5087MH, Hitec RCD, San Diego, CA, USA), controlled by a microcontroller (ELEGOO Mega 2560, Elegoo Industries, Shenzhen, China) programmed using Simulink (MathWorks, Natick, MA, USA) via two repeating sequence interpolated blocks, one for α and one for β , prescribing the angles adapted from 16 .The full CAD library and assembly as well as the list of purchased components can be accessed in the open-access repository by Oliveira Santos et al. (see oliveira:zenodo:2022).
The swimming kinematics of the Pleobot were analyzed and compared to those reported for live krill 16 .One robotic appendage was tethered to a traverse beam and submerged in a glycerin-water mixture at room temperature (60% glycerin and 40% water).Appendage motion was recorded at 125 fps using a scientific camera (FASTCAM MINI WX, Photron, 2048 pixels x 2048 pixels).Black markers on the surface of the robotic pleopod were tracked, both on the protopodite and the exopodite by digitizing the video recordings via DLTdv8 for MATLAB using automatic point tracking 38 .

Flow field measurements
The flow field around the beating appendage was measured using 2D PIV.The experiments were carried out in a cubic tank of 30 cm by side using a high-speed camera (FASTCAM MINI WX, Photron, 2048 pixels x 2048 pixels) at 125 frames per second, a Nikon lens (Nikon AF-S VR Micro-NIKKOR 50 mm), and a continuous laser (Laserglow, 1 W at 532 nm) with a cylindrical lens to create a laser sheet.The flow was seeded with 10 µm particles (Dantec Dynamics, Skovlunde, Denmark) (Supplementary Figure 3).The field of view is 174 mm x 174 mm, resulting in a spatial resolution of 0.085 mm per pixel.The velocity fields were calculated using DaVis 10 (LaVision) with decreasing interrogation window sizes (initial 48 pixels x 48 pixels, 50% overlap, 1 pass, and final 32 pixels x 32 pixels, 50%, 3 passes).Standard vector post-processing was performed to remove outliers in the flow field.The vorticity field was calculated from the velocity field.

Force measurements
The force measurements were conducted on a scaled-up (20x) model of the distal appendage by using a 6-axis force transducer (Nano 17 F/T transducer, ATI) set up in a cubic tank of 30 cm by side (Supplementary Figure 4).The Reynolds number was matched to ensure dynamically scaled experiments.To this end, the beating frequency and the fluid viscosity were modified from the baseline.A glycerin-water mixture (60% glycerin and 40% water) was used with a dynamic viscosity of 9 cSt, measured with a standard rheometer (Ares-G2, TA Instruments).The appendage was designed to be neutrally buoyant in the water-glycerin mixture.
The force transducer was mounted in the axis of rotation of the bi-ramous pleopod, such that the x-axis was aligned with the pleopod, and the y-axis was perpendicular to it.The forces measured, parallel (F ) and perpendicular (F ⊥ ) to the pleopod, are decomposed into lift and thrust by taking into account the orientation of this section in the body via the angle Ψ, the overall angle between the body axis and the endopodite, defined as α + 180 o − β .The lift and thrust coefficients (C L and C T , respectively) are calculated using the lift and thrust forces measured from the decomposed components (Equation 6).Using the instantaneous force along a cycle, averaged over five cycles, we calculate lift and thrust as: and C L and C T as where ρ is the density of the glycerin-water mixture, U is the appendage tip velocity, and A the area of the abducted appendage.

Organism experiments
Live ghost shrimp Palaemonetes paludosus were used to validate the hydrodynamics and γ angle kinematics of the Pleobot.Ghost shrimp are not obligate swimmers like krill, but they are appropriate analogs for several reasons.First, they have a ubiquitous reliance on metachronal swimming.Second, they have comparable morphology and swimming kinematics (i.e., parameters related to their exopodites, endopodites, and setae) to krill.Finally, they are similar in body size and anatomy.
Specimens (n = 4) were tethered to a vertical wire using cyanoacrylate applied to the dorsal surface of their carapace and hung in a filming vessel (15 × 10 × 5 cm 3 ) brought to 100% oxygen saturation.We performed bright-field PIV by seeding the water with 10 µm particles (Dantec Dynamics, Skovlunde, Denmark) and back-lighting with a fiber optic illuminator (Fiber-Lite MI-152, Dolan-Jenner Industries, Boxborough, MA, USA) coupled with a collimating lens (N-BK7 Plano-Convex Lens, Thorlabs, Newton, NJ, USA) 39 .We recorded the lateral and posterior views using a high-speed digital video camera (Fastcam Nova R2, Photron, Tokyo, Japan) at 2000 fps and a resolution of 2048 × 1472 squared pixels.
Velocity vectors were calculated using the DaVis 10 software package (LaVision, Göttingen, Germany).Image pairs were analyzed with three passes of overlapping interrogation windows (75%) with decreasing size (96 × 96 to 64 × 64 squared pixels).All frames were used for analysis, yielding a time separation between frames of 0.5 ms.Masking of the body and pleopods before image interrogation confirmed the absence of surface artifacts in the PIV measurements.

Figure 1 .
Figure 1.Morphology and kinematic parameters of the pleopod.Panels A (Euphausia superba), C (Palaemonetes paludosus) and E show the kinematic parameters of free-swimming shrimp incorporated in the Pleobot presented in panel B, D and F: α is defined as the angle between the axis of the body and the proximal segment (protopodite), β is the angle between the protopodite and the distal biramous segment (formed by the endopodite and the exopodite), γ appears during the power stroke as the exopodite and endopodite separate, and ζ characterizes the cupping formed between the exopodite and the endopodite.The Pleobot (panel G) is engineered based on a mechanical gear train to actively control α and β , while passively integrating γ.Note that panel C was captured from the back of the organism, while panel E shows a close up view of the pleopods from the side (as in panel A).

Figure 2 .
Figure 2. Kinematic measurements.Validation of the Pleobot was achieved by tracking its motion and plotting each kinematic parameter alongside biological data 16 .Panels A-C show the evolution of kinematic angles α, β , and γ.Note that the standard deviation over 5 cycles is indicated as shaded regions.Panels D and E show the evolution of these angles from the lateral (panel D) and frontal (panel E) views.

Figure 3 .
Figure 3. Flow field measurements.The planes at which PIV data were acquired are presented in panel A. Representative instantaneous velocity and vorticity fields are shown in panels: B for the vertical plane I, C for the frontal plane II, D for the horizontal plane III, and E for the vertical plane IV.Flow field measurements acquired from a tethered live shrimp are presented in panel F for reference.The black regions in panels B, and D-E are the cross-sectional profiles of the Pleobot corresponding to the planes in A and the white area in panel E represents a shadow area of the laser.In C, the Pleobot is out-of-plane and its position is superimposed for clarity.Note that the size and direction of the arrows scale with the value of the velocity at a given grid point, while the colors scale with the value of the vorticity according to the colormap.The red scale arrow represents 20 cm s −1 .The scale bar represents 1 cm in (B-E) and 1 mm in (F).

Figure 4 .
Figure 4. Force measurements.The thrust (A) and lift (B) coefficients are averaged over five consecutive cycles.The orange dashed lines show the kinematics of the appendage, characterized by the angle Ψ. Corresponding pleopod profiles in panel C for the anterior view of the γ angle emphasize the change in pleopod surface area over a beat cycle .Panels D-I show PIV measurements at different times of a cycle, corresponding to the black dashed lines in panels A and B. Panels D and E correspond to the lateral view along the endopodite during the power stroke, showing a tip and bound vortex (D) that are shed at the end of the power stroke (E).Panels F and G show a bottom view of the exopodite and endopodite during the recovery stroke, with strengthening spanwise flow behind the appendages between two side edge vortices.Panels H and I correspond to a lateral view along the midplane of the exopodite during the power stroke, where we identify a leading-edge vortex (LEV) that initially contributes primarily to lift and then thrust as the AoA changes.The scale bar and the red scale arrow represent 1 cm and 20 cm s −1 , respectively for all frames.

Figure 5 .
Figure 5. Scaled KRill Inspired Modular Pleobots (SKRIMP).This sketch illustrates how the modularity and rapid prototyping of the Pleobot enables the investigation of metachronal swimming by employing several units.This includes fluid flows and forces -in several taxa, body plans, and appendage configurations.SKRIMP constitutes the baseline to engineer a new generation of AUVs operating in complex marine environments leveraging the swimming characteristics of metachronal swimmers.