The underwater Brachistochrone

This paper formulates and solves the brachistochrone problem for bodies moving through dense fluids by incorporating buoyancy, viscous drag, and added mass effects, revealing that the classical cycloid becomes suboptimal or unreachable as body density approaches fluid density and establishing a critical framework for planning short-range trajectories of buoyancy-driven underwater vehicles.

The Gist

Imagine you are trying to slide a toy car from the top of a hill to a specific spot at the bottom as fast as possible.

In a perfect, frictionless world (like a cartoon), the fastest path isn't a straight line down, nor is it a gentle curve. It's a specific, swooping curve called a cycloid. This is the "Brachistochrone" problem, a famous puzzle solved by math geniuses like Newton and Bernoulli over 300 years ago.

But what happens if you try to do this underwater?

This paper, written by Mohammad-Reza Alam, asks that exact question. It turns out that underwater, the rules of the game change completely because of three new "villains" that don't exist in the air: Buoyancy (water pushing up), Drag (water pushing back), and Added Mass (the water you have to push out of the way).

Here is the breakdown of the paper's findings using simple analogies:

1. The "Heavy" vs. The "Light" Swimmer

The most important factor is how heavy your object is compared to the water.

• The Heavy Rock (High Density): If you drop a heavy rock, it barely notices the water. It plummets fast. In this case, the classic "cycloid" curve is still the best path. The water is just a minor annoyance.
• The Light Sponge (Low Density): If you drop a sponge that is almost as light as water, the water fights back hard. The "cycloid" curve becomes a terrible idea. Why? Because the cycloid dips deep to gain speed, but then has to climb back up. In water, climbing up is a nightmare because the water drags on you, stealing all your speed.
  • The Result: For light objects, the fastest path is actually much flatter and straighter. You don't want to dive deep and fight your way back up; you want to skim along a shallow path to avoid losing energy.

2. The "Added Mass" Ghost

This is the paper's biggest new discovery. When you accelerate underwater, you aren't just moving your body; you are also dragging a cloud of water along with you.

• The Analogy: Imagine trying to run through a crowded hallway. You aren't just moving your own legs; you have to push the people around you out of the way. That crowd has "weight."
• The Effect: This "crowd" (added mass) makes the object feel much heavier than it actually is. If you ignore this, you will think your underwater glider will arrive in 10 seconds, but in reality, it will take 12 or 13 seconds. The paper shows that ignoring this "ghost weight" leads to a 20% error in your time predictions.

3. The "Drag Crisis" (The Magic Switch)

Water resistance isn't constant. It behaves like a light switch that flips at a specific speed.

• The Analogy: Imagine driving a car. At low speeds, the air feels thick and sticky. But if you hit a specific speed (around 200,000 Reynolds number for a sphere), the air suddenly "slips" around the car, and drag drops by half!
• The Paper's Finding: For objects of a certain size and weight, they might hit this "magic speed" right in the middle of their trip.
  • If they hit it, they suddenly zoom forward.
  • If they miss it (because they are slightly heavier or lighter), they stay stuck in the "sticky" slow zone.
  • The Danger: If you plan a path assuming the drag is constant, you might design a route that looks perfect on paper but fails in real life because the object never hits that "magic switch" speed.

4. The "Dead End" Zone

In the old math problem, you could go from Point A to Point B no matter where B was. Underwater, that's not true.

• The Analogy: Think of a hiker with a limited battery. If they walk too far uphill, they run out of battery and get stuck.
• The Finding: If your underwater glider is too light, there are some destinations it simply cannot reach. If it dives too deep to gain speed, the water drag will steal all its energy before it can climb back up to the target. The paper maps out a "reachable zone"—a safe area where the glider can actually get to. Outside that zone, the mission is impossible without extra power.

5. The "Three-Point" Challenge

The paper also looked at a harder version: "Go from A, touch a specific point M (to avoid a rock), and then go to B."

• The Result: In the air, you can always do this. Underwater, the "Dead End" zone gets bigger. If you have to stop at a specific point, you might not have enough energy left to finish the trip. The path you take to get to the middle point might leave you stranded before you reach the end.

Why Does This Matter?

This isn't just a math puzzle; it's a manual for underwater robots (gliders).

• These robots are used to map the ocean, track oil spills, and monitor climate change.
• They run on batteries, so every second counts.
• If engineers use the old "cycloid" math, they might program the robot to take a path that takes 20% longer or, worse, gets the robot stuck in the middle of the ocean.
• By using this new "underwater Brachistochrone" math, engineers can plot paths that are faster, safer, and more energy-efficient, especially for robots that are designed to be very light and buoyant.

In short: The fastest way down in a vacuum is a swooping curve. The fastest way down underwater is a smarter, flatter path that respects the water's weight, its stickiness, and the fact that you have to push a cloud of water along with you.

Technical Summary

1. Problem Statement

The paper addresses the underwater brachistochrone problem: determining the trajectory $y(x)$ that minimizes the transit time for a body moving from a starting point $A$ to a target point $B$ in a dense fluid.

Unlike the classical vacuum brachistochrone (which yields a cycloid), the underwater problem is fundamentally altered by three physical effects:

1. Buoyancy: Reduces the effective gravitational driving force.
2. Viscous Drag: Dissipates kinetic energy, preventing the body from accelerating indefinitely.
3. Added Mass: The inertial effect of the fluid entrained by the moving body, which increases the system's effective inertia.

The study specifically focuses on buoyancy-driven underwater gliders (e.g., autonomous vehicles) where the body density ($\rho_b$) is comparable to the fluid density ($\rho_f$), making the added mass effect significant. The problem is formulated for a sphere moving in a quiescent fluid, accounting for a Reynolds-number-dependent drag coefficient ($C_d$) to capture the drag crisis phenomenon.

2. Methodology

Governing Equations

The dynamics are derived from the Basset-Boussinesq-Oseen (BBO) equation, simplified for a body moving monotonically forward. The authors retain:

• Gravity and Buoyancy: Net driving force $(\rho_b - \rho_f)Vg \sin\theta$.
• Viscous Drag: $F_d = \frac{1}{2} C_d \rho_f A v |v|$, where $C_d$ is a function of the Reynolds number ($Re$).
• Added Mass: $F_a = c_m \rho_f V \dot{v}$, where $c_m$ is the added mass coefficient (0.5 for a sphere).
• Note: The Basset history force (memory term) and normal added mass components are neglected as they are negligible for trajectory optimization at moderate-to-high Reynolds numbers or do not affect the tangential equation of motion.

The tangential equation of motion is derived as:
$$\dot{v} (\gamma + c_m) = (\gamma - 1) g \sin \theta - \frac{1}{2} \frac{C_d}{L} v |v|$$
where $\gamma = \rho_b/\rho_f$ is the density ratio, $L = V/A$ is the volume-to-area ratio, and $\theta$ is the path angle.

Dimensionless Formulation

The equations are non-dimensionalized using characteristic length $L$ and time $\sqrt{L/g}$. The problem is governed by:

• Density ratio ($\gamma$)
• Added mass coefficient ($c_m$)
• A dimensionless parameter $G$ related to the Reynolds number and object size.
• Endpoint coordinates $(x_e, y_e)$.

The drag coefficient $C_d$ is modeled using a phenomenological correlation for a smooth sphere that captures the sharp drop in drag (drag crisis) around $Re \approx 2 \times 10^5$.

Numerical Solution

• Optimization: The trajectory $y(x)$ is represented as a cubic spline with control points. The transit time is minimized using the Nelder-Mead simplex method (derivative-free).
• Integration: Equations of motion are integrated forward using ode45 with high tolerance ($10^{-13}$).
• Extensions: The framework is extended to a three-point brachistochrone (passing through an intermediate waypoint $M$) to analyze reachability constraints.

3. Key Contributions

1. First Inclusion of Added Mass: This is the first formulation of the brachistochrone problem that explicitly incorporates the added mass effect. The authors demonstrate that neglecting added mass leads to significant errors in predicted transit times, especially when $\rho_b \approx \rho_f$.
2. Reynolds-Number Dependent Drag: Unlike previous studies assuming constant drag coefficients, this work utilizes a variable $C_d(Re)$, revealing the critical impact of the drag crisis on optimal trajectories.
3. Decomposition of Physical Effects: The paper quantifies the individual contributions of drag and added mass, showing that neglecting both yields transit time predictions roughly 50% lower than reality, while neglecting added mass alone underestimates time by ~20%.
4. Finite Reachable Domain: In the three-point problem, the authors identify a finite reachable domain for the endpoint. Unlike the classical vacuum problem where any point is reachable, underwater drag creates energy constraints where certain endpoints become inaccessible after passing a waypoint.

4. Key Results

Trajectory Shapes and Density Ratios ($\gamma$)

• High Density ($\gamma \gg 1$): Drag is negligible; the optimal path is virtually identical to the classical cycloid.
• Low Density ($\gamma \to 1$): The net driving force weakens. The optimal path becomes shallow and nearly straight to minimize drag losses. The classical cycloid becomes suboptimal (or even fails to reach the target) because its deep excursion incurs excessive drag on the ascent.
• Critical Density Ratio ($\gamma \approx 1.4$): Trajectories are acutely sensitive to $\gamma$ due to the drag crisis. Small changes in density shift the trajectory such that it crosses the critical Reynolds number ($Re \approx 2 \times 10^5$), causing a sudden drop in $C_d$ and a significant reduction in transit time.

Quantitative Findings

• Improvement over References: For $\gamma \approx 1.1$, the optimal path is 24% faster than the classical cycloid and 7% faster than a straight line. For $\gamma \approx 1.5$, the improvement over a straight line peaks at ~27%.
• Drag Crisis Sensitivity: For objects operating near the drag crisis, a constant-$C_d$ approximation produces qualitatively incorrect paths. A trajectory planned with pre-crisis drag coefficients may fail to reach the destination entirely.
• Size Dependence: The transit time is not a universal function of $\gamma$ alone; it depends on the physical size of the object (via the Reynolds number). Intermediate-sized spheres exhibit a "kink" in transit time vs. density curves due to the transition through the drag crisis.

Three-Point Problem

• The inclusion of an intermediate waypoint introduces a reachability constraint.
• There exists a "black region" in the endpoint plane where targets cannot be reached because the object dissipates all kinetic energy via drag before arriving.
• This provides a practical tool for mission planning, delineating feasible target zones for underwater gliders based on vehicle parameters.

5. Significance and Applications

• Underwater Vehicle Planning: The results provide a simple yet rigorous planning tool for short-range trajectories of buoyancy-driven gliders. It highlights that standard vacuum-based or constant-drag models can lead to mission failure or significant time overestimation.
• Design Implications: The study suggests that for vehicles operating near neutral buoyancy, optimizing the trajectory is as critical as optimizing the vehicle's shape. It also motivates the use of variable density (ballast pumping) to overcome reachability limits.
• Future Directions: The framework sets the stage for more complex optimizations, including:
  • Variable density schedules (active buoyancy control).
  • Morphing vehicle shapes (changing $C_d$ and $c_m$ along the path).
  • Extraterrestrial applications (e.g., probes in the subsurface oceans of Europa), where gravity and fluid density vary significantly.

In conclusion, the paper establishes that the underwater brachistochrone is a complex, multi-parameter optimization problem where added mass and variable drag are not minor corrections but dominant factors that fundamentally alter the optimal path and the feasibility of reaching specific targets.