Scaling Laws for Caudal Fin Swimmers Incorporating Hydrodynamics, Kinematics, Morphology, and Scale Effects

This study derives and validates physics-based scaling laws for thrust, power, and efficiency in caudal fin swimmers by utilizing high-fidelity simulations and a leading-edge vortex model to analyze the interplay between hydrodynamics, kinematics, morphology, and scale effects.

The Gist

Imagine you are watching a school of mackerel dart through the ocean. They move with a fluid, wave-like motion that starts at their nose and ripples all the way to their tail, which whips back and forth like a propeller. For decades, scientists have tried to figure out the "secret sauce" of how these fish swim so efficiently. How big should the tail be? How fast should it wiggle? Does a tiny baby fish swim differently than a giant whale shark?

This paper is like a master blueprint for underwater swimming. The authors, researchers from Johns Hopkins University, used powerful computer simulations to build a digital mackerel and test it in a virtual ocean. They didn't just watch the fish; they looked inside the water to see the invisible forces at play.

Here is the story of their discovery, broken down into simple concepts:

1. The Invisible "Grip" (The Leading-Edge Vortex)

Imagine you are running with an umbrella. If you hold it just right, the wind creates a low-pressure pocket that pulls the umbrella forward. Fish do something similar with their tails.

When the fish's tail whips through the water, it creates a swirling tornado of water right at the front edge of the tail. The researchers call this a Leading-Edge Vortex (LEV). Think of this vortex as a magnetic grip. The water swirls around the tail, creating a suction effect that pulls the fish forward. The paper proves that this "swirling grip" is the main engine of the fish's speed. Without it, the tail would just push water away without much forward motion.

2. The "Goldilocks" Zone (Scaling Laws)

The researchers wanted to know: Does the size of the fish change the rules of swimming?

They found that swimming isn't just about being big or small; it's about the relationship between three things:

• How fast the tail beats (Frequency).
• How wide the tail swings (Amplitude).
• How fast the fish is moving (Speed).

They discovered a "Goldilocks" zone. If the tail swings too slowly or too wildly, the fish wastes energy. If it swings just right, the fish hits a sweet spot of efficiency. They created a mathematical "recipe" (scaling laws) that predicts exactly how fast a fish can go based on its shape and how it moves its tail.

3. The "Slip" Factor (The Wave vs. The Fish)

Imagine a surfer riding a wave. The wave moves forward, and the surfer moves forward.

• The Wave Speed: This is how fast the ripple travels down the fish's body.
• The Fish Speed: This is how fast the fish actually moves through the water.

The paper introduces a concept called "Slip." It's the difference between the wave speed and the fish speed.

• If the fish swims too fast (faster than the wave), the tail starts acting like a brake, slowing the fish down.
• If the fish swims too slow, it's just flailing uselessly.
• The Sweet Spot: The most efficient fish swim at about 80-90% of the wave speed. This is the "slip ratio" that nature has optimized over millions of years.

4. The Shape of the Wave (The "A-star" Parameter)

Here is the most surprising finding. It's not just about how much the tail moves, but how the wave grows along the body.

Imagine a wave traveling down a rope.

• Scenario A: The wave gets bigger very smoothly and gradually.
• Scenario B: The wave gets bigger, but then suddenly jerks or changes angle right at the tip of the tail.

The researchers found that Scenario A is much more efficient. They identified a specific number (which they call A') that measures this "smoothness" or "phase mismatch" at the tail.

• The Biological Twist: Real fish (like mackerel) don't have the perfect smooth wave (A' = 0). They have a slight "jerk" at the tail (A' ≈ 0.4). Why? Because fish are made of muscle and bone, not computer code. Their tails are passive and floppy at the end. They can't perfectly control the angle of the tail tip.
• The Lesson for Robots: If we want to build the perfect underwater robot, we shouldn't just copy a fish's shape exactly. We should design the robot's tail to have that "perfectly smooth" wave (A' = 0), even if real fish can't do it. This could make our robots much faster and more energy-efficient.

5. The Size Matters (Reynolds Number)

Finally, the paper explains why a tiny fish (like a baby zebrafish) swims differently than a giant tuna.

• Tiny Fish: They swim in "thick" water (like swimming in honey). The water sticks to them, creating a lot of friction. They have to wiggle their tails very fast and wide to get moving.
• Big Fish: They swim in "thin" water (like flying through air). The water slides off them easily. They can swim faster with less frantic movement.

The researchers' new formulas connect these two worlds. They show that as fish get bigger, the "rules" of swimming shift, but they all follow the same underlying physics of that swirling vortex grip.

Why Does This Matter?

This isn't just about fish. This is a design manual for the future.

• For Engineers: If you are building a bio-robotic underwater vehicle (like a robot fish for ocean exploration), you can use these formulas to know exactly how big to make the tail, how fast to make it beat, and what shape the body should be to save battery life.
• For Biologists: It helps us understand why fish evolved the way they did. They aren't just "randomly" shaped; their bodies are perfectly tuned to the physics of water, constrained only by their muscle limitations.

In a nutshell: The fish are nature's masters of hydrodynamics, using swirling water tornadoes to pull themselves forward. By decoding their "secret recipe," we can build better underwater robots and understand the physics of movement in a whole new way.

Technical Summary

1. Problem Statement

While body-caudal fin (BCF) propulsion is a dominant mode of locomotion for many fish and bio-robotic underwater vehicles (BUVs), a comprehensive understanding of how scale (Reynolds number), kinematics, and morphology interact to determine swimming performance remains incomplete.

• The Gap: Previous studies have proposed scaling laws for swimming speed and efficiency (e.g., relating Strouhal and Reynolds numbers), but these often lack a direct mechanistic link to the specific hydrodynamic forces generated by the caudal fin.
• The Challenge: Existing models often rely on idealized flow assumptions or dimensional analysis without accounting for the complex vortex dynamics (specifically Leading-Edge Vortices) that dominate thrust generation in BCF swimmers. Furthermore, the influence of specific kinematic parameters, such as the amplitude envelope slope at the tail, on optimal performance has not been fully quantified.

2. Methodology

The authors employed High-Fidelity Direct Numerical Simulations (DNS) to analyze a mackerel-inspired carangiform swimmer across a wide range of conditions.

• Computational Setup:

  • Solver: Vicar3D, a sharp-interface immersed boundary method solver for incompressible Navier-Stokes equations.
  • Geometry: A 3D model of a mackerel (Scomber scombrus) with a zero-thickness membrane caudal fin.
  • Kinematics: Prescribed undulatory motion where lateral displacement increases quadratically from head to tail.
  • Range: Simulations covered Reynolds numbers ($Re_L$) from 500 to 50,000 and various Strouhal numbers ($St_A$).
  • Validation: Results were validated against experimental data (e.g., Videler & Hess) and prior computational studies (e.g., Borazjani & Sotiropoulos).

• Analytical Framework:

  • Force Partitioning Method (FPM): Used to decompose forces into pressure and viscous components, identifying that thrust is primarily generated by pressure forces on the caudal fin.
  • LEV-Based Model: The study extends a Leading-Edge Vortex (LEV) model (previously applied to flapping foils) to the caudal fin. This model posits that thrust is proportional to the circulation generated by the LEV, which depends on the effective angle of attack and pitching/heaving motions.
  • Derivation: The authors derived analytical scaling laws for thrust, power, efficiency, and cost-of-transport (COT) by integrating the LEV model with the specific kinematics of BCF swimming.

3. Key Contributions

The paper introduces several novel parameters and scaling relationships that bridge the gap between fluid physics and biological/robotic performance:

• Identification of $A'^*$: The authors define a new non-dimensional kinematic parameter, $A'^*$, which represents the slope of the amplitude envelope at the tail ($dA/dx$) normalized by wavelength. This parameter quantifies the phase mismatch between the heaving velocity and the pitching angle of the caudal fin.
• Mechanistic Scaling Laws: Derivation of scaling laws for:
  • Thrust ($C_T$): Linked to the product of $\sin(\alpha_{eff})$ and $\sin(\theta)$.
  • Power ($C_W$): Linked to the work done against pressure loads.
  • Efficiency ($\eta$): Shown to be a function of the slip ratio ($U/U_c$) and $A'^*$.
  • Drag ($C_D$): Scaled as $Re_U^{-2/3}$ for terminal swimming, differing from the classic laminar $Re^{-1/2}$ scaling due to the moving boundary layer.
• Morphological Parameter ($K_{morph}$): Introduction of a parameter combining body drag, fin shape, and relative sizes ($K_{morph} = 2(C_f/\beta_T)(S_b/S_f)$) to predict the relationship between Strouhal and Reynolds numbers.
• Optimal Kinematics: Demonstration that the optimal Strouhal number ($St_{opt}$) and maximum efficiency are strictly determined by $A'^*$.

4. Key Results

A. Wake Topology and Strouhal Number

• The wake structure transitions from a "double vortex row" at low Reynolds numbers (high $St_A$) to a "single vortex row" (or narrow wake) at high Reynolds numbers (low $St_A$).
• The wake spreading angle ($\theta_w$) is directly proportional to the Strouhal number: $\theta_w \approx \tan^{-1}(St_A/2)$.
• At high $Re$, the free-swimming $St_A$ approaches a minimum value ($St_{min}$), resulting in a very narrow wake (e.g., $\sim 8^\circ$ at $Re_U = 36,000$).

B. Thrust and Power Scaling

• Thrust: The thrust coefficient scales with $St_A$ but deviates from the simple $St_A^2$ scaling at low Strouhal numbers. A minimum Strouhal number ($St_{min}$) is required to generate positive thrust; below this, the fin produces drag.
• Power: The power coefficient scales roughly with $St_A^3$ but includes corrections based on the amplitude envelope slope.
• Efficiency: The Froude efficiency peaks at an optimal Strouhal number ($St_{opt}$). Crucially, efficiency decreases as $A'^*$ increases.
  • For the mackerel model ($A'^* \approx 0.38$), peak efficiency occurs at $St \approx 0.31$.
  • If $A'^*$ could be reduced to 0 (perfect phase alignment), efficiency would increase significantly, and the optimal $St$ would decrease.

C. The $St_A$ vs. $Re_U$ Relationship

The study derives a unified scaling law connecting Strouhal and Reynolds numbers for terminal swimming:
$$\frac{St_A}{St_{min}} = \frac{1}{2} + \frac{1}{2} \sqrt{1 + \left(\frac{Re_{cr}}{Re_U}\right)^{2/3}}$$

• **Low $Re$($Re_U \ll Re_{cr}$):**$St_A \sim Re_U^{-1/3}$.
• **High $Re$($Re_U \gg Re_{cr}$):**$St_A \to St_{min}$ (constant).
• Critical Reynolds Number ($Re_{cr}$): A threshold (approx. 42,600 for mackerel) where viscous effects become negligible. This value depends on morphology ($K_{morph}$) and kinematics.

D. Biological and Robotic Implications

• Biological Observation: Real fish typically operate with $A'^* \approx 0.4$, which is sub-optimal for pure hydrodynamic efficiency. The authors suggest this is a trade-off: biological swimmers use passive elastic tail motion (which naturally creates a non-zero $A'^*$) to reduce muscular effort, sacrificing some peak efficiency for control and energy economy.
• BUV Design: To achieve optimal performance, BUVs should not simply scale geometry linearly.
  • Smaller vehicles should use relatively larger caudal fins.
  • Larger vehicles should use relatively smaller fins (relative to body size) to maintain optimal $K_{morph}$.
  • Kinematics should be tuned to minimize $A'^*$ (phase mismatch) if maximum efficiency is the goal, though this may require active control of the tail tip.

5. Significance

This work provides a mechanistic framework for understanding aquatic locomotion that moves beyond empirical correlations.

1. Predictive Power: The derived scaling laws allow researchers to predict swimming speed, thrust, and efficiency for any BCF swimmer given its morphology and kinematics, without needing full CFD simulations.
2. Design Guidelines: It offers specific design rules for bio-inspired underwater vehicles, emphasizing the coordination of body shape, fin size, and kinematic phase to optimize performance across different scales.
3. Biological Insight: It resolves the paradox of why fish do not swim at the theoretically "most efficient" kinematic configuration, attributing it to the constraints of passive elastic mechanics and neuromuscular control.
4. Unified Theory: It successfully unifies low-Re (viscous-dominated) and high-Re (inertial-dominated) regimes into a single scaling law, bridging a gap in the literature regarding how swimming performance transitions with size.