🌀 nonlinear sciences

Superposition of shock waves of the generalized BBM equation

This paper investigates the generalized BBM equation with an additional dissipative term, deriving a stable two-parameter family of travelling shock wave solutions and establishing effective superposition rules that apply to both continuous and discontinuous waves.

Original authors: Alexey Samokhin

Published 2026-03-16

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Alexey Samokhin

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 watching a river. In some theoretical models of water waves, if you create a perfect, solitary wave (a "soliton"), it travels forever without changing shape, like a ghost ship that never fades. This is the classic BBM equation model.

But in the real world, water isn't perfect. It has friction, it loses energy to the air, and it gets "sticky." This paper introduces a Generalized BBM equation that adds a "dissipation" term—essentially, a mathematical way of saying, "Hey, there's friction here."

Because of this friction, the perfect ghost ships disappear. Instead, we get Shock Waves. Think of these not as smooth, solitary humps, but as traffic jams or a sudden change in water level that moves down the river. They have a high side and a low side, and they eventually settle into a steady shape.

Here is the breakdown of what the author, Alexey Samokin, discovered about these waves, explained simply:

1. The Waves Have a "Personality" (Stability)

The paper first asks: If I disturb a shock wave, will it fall apart, or will it fix itself?

The answer is it fixes itself. The author proves that these shock waves are incredibly stable. Imagine you have a perfectly formed traffic jam moving down a highway. If you suddenly add a few extra cars or take a few away (a "perturbation"), the traffic jam might get messy for a moment. But because of the "friction" in the system, the chaos eventually smooths out, and the jam returns to its original, steady shape.

The math behind this is like a "conservation law" that acts as a self-correcting mechanism. The system naturally "dissipates" the extra energy until it finds its comfortable, steady state again.

2. The "Magic" of Adding Waves (Superposition)

This is the most exciting part of the paper. In normal physics, if you add two waves together, they usually just pass through each other or create a messy interference pattern. But here, the author found a surprisingly simple rule for what happens when two shock waves crash into each other.

The Analogy: The Vector Sum
Imagine each shock wave is an arrow pointing in a specific direction with a specific length.

Wave A has a "high side" of 5 and a "low side" of 2.
Wave B has a "high side" of 3 and a "low side" of 1.

When these two waves collide and merge, they don't just bounce off. They add up like vectors.

The new "high side" becomes $5 + 3 = 8$ .
The new "low side" becomes $2 + 1 = 3$ .

The two waves merge into a single, new shock wave that has these new boundary values. It's like two teams of people merging into one bigger team; the new team's size is just the sum of the old teams.

3. The Speed of the New Wave

The paper also gives a simple rule for how fast this new, merged wave will travel.
If Wave 1 travels at speed $V_1$ and Wave 2 travels at speed $V_2$ , the new merged wave travels at:
$V_{new} = V_1 + V_2 - 1$
(Note: The "-1" is just a mathematical quirk of how this specific equation is written, like a tax you have to pay when two things combine.)

4. Smooth vs. Bumpy (Oscillations)

Not all merged waves look the same.

Monotonic (Smooth): Sometimes, when two waves merge, the result is a smooth, S-shaped curve. It's like a gentle ramp.
Oscillatory (Bumpy): Sometimes, the result is a wave that wiggles back and forth (like a spring) before settling down. The paper provides a way to predict whether the result will be smooth or bumpy based on the height of the waves and the "friction" of the medium.

5. Real-World Examples (The Computer Simulations)

The author ran computer simulations to show this in action:

The "Turbulent Union": Two waves crash, create a messy, wiggly mess, and then slowly calm down into a single, stable wave.
The "Peaceful Fusion": Two waves merge so smoothly that they just glide into a new shape without much fuss.
The "Standstill": Two waves moving in opposite directions crash and cancel each other out perfectly, leaving a stationary wave that doesn't move at all.

Why Does This Matter?

This isn't just about math equations. The rules discovered here apply to any shock wave, even ones that are jagged or discontinuous (like a sudden cliff in a graph), not just the smooth ones.

This is useful for understanding:

Traffic flow: How two jams merge into one.
Fluid dynamics: How waves in oil pipelines or water channels interact.
Plasma physics: How energy moves through hot gases in stars or fusion reactors.

In a nutshell: The paper tells us that when shock waves collide in a "friction-filled" world, they don't destroy each other. Instead, they merge into a new, stable wave whose size and speed can be predicted by simply adding their parts together. It's a beautiful example of how chaos (the collision) naturally organizes itself into order (the new wave).

1. Problem Statement

The paper investigates the Generalized Benjamin-Bona-Mahony (gBBM) equation, which extends the classic BBM equation by incorporating a dissipative term. The classic BBM equation ( $u_t + u_x + uu_x - \lambda u_{xxt} = 0$ ) is known for supporting solitary wave (soliton) solutions. However, the inclusion of the dissipative term ( $-\varepsilon u_{xx}$ ) in the generalized form:
$u_t + u_x + uu_x - \lambda u_{xxt} - \varepsilon u_{xx} = 0$
fundamentally alters the nature of the solutions. Instead of solitons, the system generates traveling shock waves (both continuous and discontinuous).

The core problem addressed is the interaction and superposition of these traveling shock waves. Unlike linear systems where superposition is additive, or integrable nonlinear systems (like KdV) where solitons pass through each other preserving their shape, the behavior of shock waves in dissipative systems is complex. The paper seeks to define the rules governing how two such waves merge and what the resulting wave profile and velocity will be.

2. Methodology

The author employs a combination of analytical derivation, dynamical systems theory, and conservation laws:

Traveling Wave Reduction: The partial differential equation (PDE) is reduced to an Ordinary Differential Equation (ODE) by assuming a traveling wave solution of the form $u(x,t) = u(x - Vt)$.
Phase Plane Analysis: The resulting ODE is analyzed as a dynamical system. The author identifies fixed points corresponding to the boundary conditions at infinity ( $H$ and $h$ ) and classifies them (saddle vs. focus) to determine whether the shock wave is monotonic or oscillatory.
Stability Analysis via Conservation Laws: The paper utilizes the momentum conservation law ( $\int u^2 dx$ ). While the classic BBM conserves momentum, the gBBM equation exhibits selective decay of momentum due to the dissipative term. The author proves that the momentum of any perturbation decays monotonically to zero, implying asymptotic stability.
Numerical Simulation: The theoretical rules are validated using computer simulations (Maple 21) to visualize the evolution of wave superposition, including cases involving discontinuous initial data (Riemann breaks).

3. Key Contributions and Results

A. Characterization of Traveling Shock Waves (TWS)

The paper derives a two-parameter family of traveling wave solutions defined by their boundary conditions at infinity, $u(-\infty) = H$ and $u(+\infty) = h$ .

Velocity Relation: The wave velocity $V$ is explicitly determined by the boundary values:
$V = 1 + \frac{H + h}{2}$
Monotonic vs. Oscillatory: A criterion is established to distinguish between monotonic and oscillatory shock waves based on the parameters $\lambda, \varepsilon, H,$ $λ, ε, H,$ and $h$ $h$ .
- If the discriminant of the characteristic equation is positive, the solution is monotonic.
- If negative, the solution is oscillatory (exhibiting damped oscillations near the boundaries).
Explicit Solutions: An explicit formula (involving hyperbolic tangents) is provided for a specific subset of these waves, though the paper notes that the general family is defined by the boundary conditions rather than a single closed-form expression.

B. Stability Theorem

Theorem 1: Every solution to the gBBM equation with specific boundary conditions $H$ and $h$ (and sufficiently smooth initial data) asymptotically converges to the unique traveling shock wave solution $T_{H,h}$ as $t \to \infty$ .

Mechanism: The proof relies on the fact that the "excess" momentum (the difference between the initial state and the target shock wave) dissipates over time due to the $\varepsilon u_{xx}$ term.

C. Superposition Rules

Based on the stability theorem, the paper derives effective rules for the superposition of two shock waves ( $T_{H_1, h_1}$ and $T_{H_2, h_2}$ ):

Additivity of Boundary Conditions: The superposition of two waves results in a new wave whose boundary conditions are the sum of the individual boundary conditions:
$T_{H_1, h_1} + T_{H_2, h_2} \to T_{H_1+H_2, h_1+h_2}$
Note: This is an asymptotic result; the intermediate state is turbulent, but the final state is a single shock wave with summed boundaries.
Velocity of the Resulting Wave: The velocity of the resulting wave is given by:
$V_{new} = V_1 + V_2 - 1$
This rule holds regardless of whether the initial waves are continuous or discontinuous (e.g., Riemann shocks).

D. General Applicability

A significant finding is that these superposition rules are not limited to exact analytical traveling wave solutions. They apply to a wider class of shock waves, including discontinuous initial data and perturbations that are not exact solutions to the equation, provided they satisfy the boundary conditions.

4. Illustrative Examples

The paper presents five numerical examples demonstrating these rules:

Turbulent Union: Two continuous waves merging into an oscillatory shock.
Discontinuous Union: The sum of two Riemann breaks (step functions) evolving into a smooth traveling shock.
Peaceful Fusion: Merging waves resulting in a monotonic (non-oscillatory) shock.
Collision: Waves moving in opposite directions colliding and forming a new wave moving in the direction of the net momentum.
Standstill: A collision resulting in a stationary shock wave ( $V=0$ ).

5. Significance

Theoretical Insight: The work bridges the gap between soliton theory (conservative systems) and shock wave theory in dissipative systems. It demonstrates that even in non-integrable, dissipative systems, predictable "superposition" rules exist based on boundary conditions.
Practical Utility: The derived rules ( $V_{new} = V_1 + V_2 - 1$ and additive boundary conditions) provide a simple, effective tool for predicting the outcome of wave interactions in physical systems modeled by the gBBM equation (e.g., long waves in liquids, acoustic-gravity waves, plasma physics).
Methodological Validation: The paper reinforces the utility of the selective decay approach (using decaying conservation laws) for analyzing the stability and long-term behavior of dissipative partial differential equations.

In conclusion, Samokhin's paper establishes that while the gBBM equation destroys the soliton property, it replaces it with a robust class of traveling shock waves whose interactions are governed by simple, additive rules derived from the system's dissipative nature.