The Extended KdV Equation: Augmented Lagrangian and Variational Solitary Waves with Applications to Dispersive Hydrodynamics

This paper extends the averaged Lagrangian method to derive explicit variational formulas for solitary wave solutions of the extended Korteweg–de Vries (eKdV) equation and validates their accuracy in predicting the leading edge of dispersive shock waves through comparison with numerical simulations.

The Gist

Imagine the ocean as a giant, restless drum. When you hit it, waves ripple out. For a long time, scientists used a simple, classic recipe (called the KdV equation) to predict how these waves move. It works great for small, gentle ripples. But what happens when the waves get bigger, or when the water is shallow and the waves start to interact in more complex, "bumpy" ways? The old recipe starts to lose its flavor.

This paper is about creating a new, upgraded recipe (called the Extended KdV or eKdV equation) that can handle these bigger, messier waves. However, this new recipe is so complicated that solving it directly is like trying to untangle a knot while wearing oven mitts.

Here is how the authors, Saleh Baqer and Hamid Said, cracked the code, explained simply:

1. The Problem: A Recipe Too Complex to Cook

The new equation describes waves that have "tails" and "bumps" that the old equation ignores. Because it's so complex, finding a perfect, exact solution (a specific wave shape that fits the math perfectly) is nearly impossible for most real-world scenarios. It's like trying to predict the exact path of a leaf in a hurricane using only a ruler.

2. The Solution: The "Augmented" Lagrangian (The Magic Constraint)

The authors used a clever mathematical trick called the Averaged Lagrangian method. Think of this as a way to find the "average" behavior of a system without getting lost in every tiny detail.

However, the standard trick didn't work for this new, complex equation because the equation doesn't follow the usual rules of energy conservation (it's "non-conservative").

To fix this, the authors invented an "Augmented Lagrangian."

• The Analogy: Imagine you are trying to balance a stack of plates on a moving truck. The "Lagrangian" is the physics of the truck. The "constraints" are the rules you must follow (like "don't drop the plates").
• Usually, you just follow the rules. But here, the rules were hidden inside the math. The authors realized they had to add the rules explicitly into their equation using "Lagrange multipliers" (think of these as invisible strings holding the rules in place).
• By tying these invisible strings to their math, they created a "Master Equation" that forced the complex wave to behave in a predictable way.

3. The Result: A Perfect "Solitary Wave"

Using this new method, they found a specific type of wave called a solitary wave (or soliton).

• What is it? Imagine a single, perfect hump of water traveling down a river without changing its shape. It's like a surfer who never falls off the board, no matter how far they go.
• The Shape: The authors proved this wave has a specific shape called a sech² profile. In plain English, it's a smooth, bell-shaped curve that looks like a perfect hill.
• The Formulas: They didn't just find the shape; they wrote down the exact math for:
  • How tall the wave is.
  • How fast it travels.
  • How wide it is.
  • Crucially, they showed that if you turn off the "extra" complex parts of the equation, their new formulas magically shrink down to become the old, classic formulas. This proves their new method is a true upgrade, not a replacement that breaks the old rules.

4. Testing the Theory: The "Shock Wave" Test

To see if their math actually works in the real world, they applied it to a Dispersive Shock Wave (DSW).

• The Analogy: Imagine a dam breaking. A wall of water rushes out. In a dispersive fluid, this wall doesn't stay a solid block; it breaks apart into a train of waves, with a big, fast wave leading the pack.
• The authors used their new formulas to predict the height and speed of that leading wave.
• The Verdict: They compared their math predictions against powerful computer simulations (which act as a "digital laboratory"). The results matched almost perfectly. Even when the waves got very large and started to develop "resonant radiation" (faint, wiggly tails that the simple math ignores), their prediction for the main wave's height and speed remained incredibly accurate.

Summary

In short, the authors took a very difficult, complex equation for water waves that was previously too messy to solve easily. They built a new mathematical "scaffolding" (the Augmented Lagrangian) to hold the complex parts in place. This allowed them to derive simple, clear formulas for the height, speed, and width of a single traveling wave. They proved these formulas work by showing they match computer simulations perfectly, even for large, complex waves.

What they did NOT do:

• They did not apply this to clinical medicine or human biology.
• They did not claim this solves the problem of all ocean waves (like tsunamis hitting a city).
• They did not invent a new physical device.
• They strictly focused on the math of fluid dynamics and verifying it against computer simulations.

Technical Summary

Technical Summary: The Extended KdV Equation: Augmented Lagrangian and Variational Solitary Waves with Applications to Dispersive Hydrodynamics

Problem Statement
The extended Korteweg–de Vries (eKdV) equation is a fifth-order nonlinear dispersive partial differential equation used to model a wide spectrum of hydrodynamic phenomena, including shallow water waves, plasma physics, and nonlinear optics. Unlike the classical KdV equation, the eKdV equation includes second-order nonlinear and dispersive terms characterized by coefficients $B_1, B_2, B_3,$ and $B_4$. A significant challenge in the study of the eKdV equation is that, for general coefficients, the system is non-integrable and does not conserve energy. Consequently, the existence of exact, closed-form solitary wave solutions is not guaranteed. While algebraic methods have been employed to find solutions for specific coefficient sets (e.g., shallow water waves), these approaches often fail to yield expressions that asymptotically reduce to the classical KdV results when higher-order terms are neglected, limiting their utility in broader dispersive hydrodynamic applications such as Dispersive Shock Waves (DSWs).

Methodology
The authors extend the method of averaged Lagrangians to the general, non-conservative eKdV equation. The core of their approach involves constructing a "master" or "augmented" Lagrangian, modeled on Luke's variational principle for shallow water waves.

1. Augmented Lagrangian Construction: The authors identify that deriving the eKdV equation from Luke's Lagrangian requires specific constraints relating the velocity potential to the wave displacement. These constraints, which arise at specific asymptotic orders (specifically involving the first-order Boussinesq equation), are incorporated into the variational formulation using the method of Lagrange multipliers. This results in an augmented Lagrangian ($L_{aug}$) that includes the governing constraints via multipliers $\lambda$ and $\sigma$.
2. Variational Solitary Wave Ansatz: A single solitary wave solution with a $\text{sech}^2$ profile is assumed, defined by an amplitude parameter $a_0$, an inverse width $w_s$, and a velocity $V_s$. The augmented Lagrangian is averaged over the traveling wave phase $\theta = x - V_s t$.
3. Derivation of Parameters: By applying the Euler-Lagrange equations to the averaged augmented Lagrangian, the authors derive explicit algebraic relations for the wave velocity $V_s$ and inverse width $w_s$ in terms of the fixed amplitude $a_0$ and the equation coefficients, retaining terms up to second-order in the small parameters $\alpha$ and $\beta$.
4. Energy Method for Height: Since the eKdV equation is non-conservative, a direct energy integral is not exact. The authors employ an asymptotic energy estimate, utilizing the classical KdV equation to approximate higher-order terms, to derive a second-order expression for the actual solitonic height $A_s$.
5. Application to DSWs: The derived variational formulas are applied to the Riemann problem for the eKdV equation. Using the "DSW equal amplitude approximation" method, the authors estimate the height and velocity of the leading solitary wave edge of a dispersive shock wave.

Key Contributions and Results

• Variational Formulation: The paper successfully constructs an augmented Lagrangian for the general non-conservative eKdV equation, overcoming the difficulty that standard Lagrangian methods are typically inapplicable to non-integrable, non-conservative systems.
• Explicit Second-Order Formulas: The authors obtain explicit formulas for the solitary wave velocity ($V_s$), inverse width ($w_s$), and height ($A_s$) up to second-order in $\alpha$ and $\beta$.
• Asymptotic Reduction: A critical feature of the derived expressions is their asymptotic reduction. When the second-order coefficients ($B_1, B_2, B_3, B_4$) are set to zero, the formulas reduce directly to the classical KdV results. This property distinguishes the current work from previous algebraic approaches (e.g., Karczewska et al.) which did not preserve this reduction.
• Numerical Validation: The theoretical predictions for solitary wave parameters and DSW leading-edge properties are compared against direct numerical simulations (using pseudo-spectral methods and Runge-Kutta integration). The results show strong agreement across a wide range of parameters. The study notes that while resonant radiation (oscillatory tails) emerges at higher parameter values, the variational $\text{sech}^2$ profile accurately captures the main wave structure and the decreasing trend in wave height caused by energy loss to radiation.
• DSW Application: The derived formulas successfully predict the leading edge velocity and height of eKdV dispersive shock waves, demonstrating the method's utility for non-integrable models where exact Whitham modulation theory is unavailable.

Significance and Claims
The authors claim that their variational approach provides a robust and general machinery for constructing single solitary wave solutions for the eKdV equation without imposing restrictions on the second-order coefficients. The primary significance of this work lies in the ability to generate analytical solutions that are both mathematically tractable (due to the $\text{sech}^2$ form) and physically consistent (due to the asymptotic reduction to KdV). This makes the derived solutions particularly valuable for studying complex dispersive hydrodynamic regimes, such as classical and non-classical DSWs, where the dependence of wave velocity on amplitude and parameters is crucial. The paper concludes that the method effectively bridges the gap between variational principles and non-integrable higher-order dispersive models, offering a practical tool for analyzing wave propagation in various physical contexts.