🔢 mathematics

A variational approach to the derivation of reduced models for bubbly flows

This paper utilizes Hamilton's least action principle to derive reduced models for bubbly flows by constraining bubble interfaces to evolve within finite-dimensional parameterized families, such as spheres, and establishes the corresponding interface conditions and well-posedness for curl-free liquid flows with homogeneous bubble pressure.

Original authors: Cosmin Burtea, David Gérard-Varet

Published 2026-05-13

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Cosmin Burtea, David Gérard-Varet

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine you are trying to predict how a swarm of bubbles moves through a glass of water. In the real world, these bubbles are messy. They wobble, stretch, squash, and change shape as they float. To describe this perfectly using standard physics, you would need to track every single point on the surface of every bubble as it deforms. It's like trying to write a script for a dance where every dancer can change their costume, height, and shape at every single second. It is mathematically possible, but it is incredibly heavy, complex, and difficult to solve.

This paper proposes a clever shortcut. Instead of tracking every wobble, the authors ask: "What if we pretend the bubbles are perfect spheres (or perfect ellipsoids) that can only grow, shrink, or move?"

Here is how they do it, explained simply:

1. The "Least Action" Rule (The Universe's GPS)

The authors use a famous principle from physics called Hamilton's Least Action Principle. Think of this as the universe's GPS. Nature always tries to take the "easiest" or most efficient path between two points. In physics, this "path" is calculated using a formula called a "Lagrangian" (a fancy way of measuring energy).

The Standard Way: If you let the bubbles change shape however they want, the GPS calculates a path where the pressure is perfectly equal on both sides of the bubble skin. This gives you the "perfect" but very complicated model.
The New Way: The authors say, "Let's tell the GPS: 'You can only choose paths where the bubbles stay round.'"

2. The Problem with the Shortcut

If you force a bubble to stay round, you run into a mathematical snag. In the real world, the pressure inside the bubble must exactly match the pressure outside. But if you force the bubble to stay round, it might want to squish into an oval shape to balance the pressure, but you are telling it, "No, stay round!"

So, the pressure on the inside and outside won't match perfectly at every single point on the surface. If you just ignore this, the math breaks.

3. The Creative Solution: The "Average" Compromise

This is the paper's main breakthrough. Instead of demanding the pressure match at every single point on the bubble's skin, the authors say: "Let's demand that the pressure matches on average for the specific ways the bubble is allowed to move."

They use a creative analogy of "degrees of freedom" (ways the bubble can move):

Moving the bubble: The bubble can move left, right, up, down.
Growing the bubble: The bubble can get bigger or smaller.

The authors prove that if you force the bubble to stay round, you don't need the pressure to be equal everywhere. You only need the total push from the inside to equal the total push from the outside, specifically for the directions the bubble is allowed to move.

Analogy: Imagine a balloon. If you squeeze it, it wants to bulge out the sides. If you force it to stay round, the pressure inside might be higher on the left and lower on the right. But, if you add up all those pushes, they must balance out the push of the water trying to squeeze it, specifically for the "move" and "grow" buttons you gave the balloon.

4. What They Found

By using this "average" rule, they derived new, simpler equations (models) for:

Spherical bubbles: Bubbles that stay perfectly round.
Ellipsoidal bubbles: Bubbles that stay perfectly egg-shaped.

These new models are much easier to solve on a computer because they don't have to track a wobbly, changing surface. Instead, they just track a few numbers: the center of the bubble, its size, and (for egg-shapes) how stretched it is.

5. Proving It Works

Finally, the authors checked if these new, simpler equations actually make sense mathematically. They looked at a specific, simplified scenario:

The water isn't swirling (it's "curl-free").
The air inside the bubbles is the same pressure everywhere.

In this case, they proved that their new equations turn into a set of standard motion equations (like the ones used to calculate how a car accelerates). They showed that if you start with a specific setup, there is exactly one correct answer for how the bubbles will move, and the math doesn't fall apart.

Summary

The paper is a guide on how to simplify a very messy physics problem (wobbly bubbles) by imposing a strict rule (bubbles must stay round or egg-shaped). The authors figured out the specific mathematical "rules of the road" (interface conditions) that allow this simplification to work without breaking the laws of physics. They proved that these simplified rules lead to a solvable, predictable system, making it much easier to simulate bubbly flows on computers.

Technical Summary: A Variational Approach to the Derivation of Reduced Models for Bubbly Flows

Problem Statement
The paper addresses the mathematical modeling of gas bubbles moving within an ambient inviscid liquid. While the full dynamics of such two-phase flows are governed by coupled Euler (or Navier-Stokes) equations with free boundaries, solving these systems is analytically and numerically demanding. The authors focus on deriving "reduced models" where the bubble interfaces are constrained to evolve within a specific subclass of hypersurfaces described by a finite number of parameters (e.g., maintaining sphericity or ellipsoidal shapes).

The central difficulty identified is that imposing such geometric constraints on the interface renders the classical boundary condition of pressure continuity ( $p_+ = p_-$ ) incompatible with the reduced dynamics. The paper seeks to determine the appropriate substitute interface conditions that allow for a well-posed system of partial differential equations (PDEs) or ordinary differential equations (ODEs) while respecting the shape constraints.

Methodology
The authors employ Hamilton's least action principle as the primary tool for derivation. The methodology proceeds in two main stages:

Recovery of the Classical Model: In Section 2, the authors formally derive the standard sharp-interface two-phase flow model (Euler equations for bulk phases coupled by continuity of normal velocities and pressures) using a Lagrangian framework. They define an action functional representing the difference between kinetic and potential energies. By considering variations of the flow maps that preserve a common interface, they recover the standard interface conditions, specifically the equality of normal stresses (pressure continuity).
Derivation of Reduced Models: In Section 3, the authors modify the variational principle by restricting the class of admissible diffeomorphisms. Instead of allowing arbitrary interface deformations, they constrain the interface $I(t)$ to be a union of surfaces drawn from a finite-dimensional manifold of shapes (e.g., spheres $S(c, r)$ or ellipsoids $E(c, S)$ ).
- Relaxation of Boundary Conditions: Because the variations are restricted to those preserving the shape class, the standard condition $p_+ = p_-$ can no longer be enforced pointwise. Instead, the authors show that the stationarity of the action leads to a set of integral conditions. Specifically, the pressure jump across the interface must be orthogonal to the space of admissible normal velocity fields associated with the shape deformations.
- Mathematical Formulation: If $V_m(x)$ represents the normal velocity induced by a tangent vector $\dot{m}$ in the parameter space $M$ , the new interface condition becomes:
  $\int_{I_k(t)} (p_- - p_+) V_{m_k}(x) \, d\sigma = 0$
  for each bubble $k$ . This replaces the pointwise pressure equality with a finite number of integral constraints corresponding to the degrees of freedom of the bubble shape.

Key Contributions and Results

General Formalism: The paper establishes a general variational framework for deriving reduced models for any shape class defined by a finite-dimensional manifold $M$ . This includes spherical and ellipsoidal bubbles as specific examples.
Explicit Reduced Systems:
- Spherical Bubbles: For $N$ $N$ spherical bubbles with centers $c_k(t)$ $c_{k} (t)$ and radii $r_k(t)$ $r_{k} (t)$ , the authors derive a system where the interface conditions are replaced by:
  1. $\int_{I_k} (p_- - p_+) \mathbf{n} \, d\sigma = 0$ (force balance)
  2. $\int_{I_k} (p_- - p_+) \, d\sigma = 0$ (volume/pressure balance)
    These conditions replace the pointwise pressure continuity.
- Ellipsoidal Bubbles: A similar derivation is provided for ellipsoidal bubbles, yielding integral conditions involving the normal vector and the shape tensor.
Well-Posedness Analysis (Section 4): The authors address the well-posedness of the derived reduced models under specific simplifying assumptions:
- The liquid flow is irrotational ( $\text{curl } \mathbf{u} = 0$ ) and fills the whole space ( $\Omega = \mathbb{R}^3$ ).
- The pressure (or density) inside each bubble is spatially homogeneous.
- Under these hypotheses, the PDE system reduces to a system of second-order ODEs governing the evolution of the shape parameters $m(t)$ .
- The authors prove the local existence and uniqueness of solutions to this ODE system by constructing a "shape-dependent" gradient field and utilizing a non-degeneracy condition on the normal velocity map.
Specific Example: For a single spherical bubble, the derived ODEs recover the classical Rayleigh-Plesset equation (when translation is neglected) and a coupled system for translation and expansion, consistent with previous works like Hermans [1973] and Wang and Smereka [2003].

Significance and Claims
The paper claims to provide a rigorous variational justification for reduced models of bubbly flows that go beyond the classical Rayleigh-Plesset framework. By utilizing Hamilton's principle, the authors ensure that the resulting reduced models inherently respect conservation laws (via Noether's theorem) and possess a clear physical interpretation.

The authors modestly note that the derivation is a matter of mathematical modeling; the existence of critical points for the action functional is a separate problem. They explicitly state that the well-posedness analysis provided is limited to the case of irrotational flows and homogeneous bubble pressures. They acknowledge that the general well-posedness of the system (3.1.7) for arbitrary flows remains an open question.

The significance of the work lies in offering a systematic method to derive interface conditions for constrained bubble shapes, which the authors suggest could be computationally advantageous for numerical simulations where interfaces are naturally approximated by a finite number of degrees of freedom. The paper does not propose new experimental applications but situates its theoretical results within the context of existing applications in acoustics, oceanography, and industry mentioned in the introduction.