🌀 nonlinear sciences

Shock waves of spherical/cylindrical KdV-B: Asymptotic, stability, superposition

This paper provides a detailed asymptotic description, stability analysis via conservation laws, and effective superposition rules for the diverging shock wave solutions of spherical and cylindrical KdV-B equations.

Original authors: Alexey Samokhin

Published 2026-04-27

📖 4 min read☕ Coffee break read

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 standing on a vast, calm lake. If you throw a single stone into the water, you see ripples spreading outward in perfect circles. Now, imagine instead of a lake, you are looking at a massive, expanding explosion or a ripple moving through a thick, jelly-like substance.

This paper, written by Alexey Samokhin, is about a specific type of "mathematical ripple" called a KdV-B shock wave. Specifically, it looks at how these ripples behave when they aren't just moving in a straight line (like a wave in a narrow canal), but are instead expanding outward like a growing bubble (spherical) or a widening ring (cylindrical).

Here is the breakdown of the paper using everyday concepts:

1. The "Shape-Shifter" Waves (The Problem)

In a simple world, waves are predictable. But in the world of the KdV-B equation, two forces are constantly fighting:

Dispersion: This is the "spreading" force. It tries to pull the wave apart, making it wider and thinner.
Dissipation: This is the "friction" force. It tries to smooth the wave out, like honey slowing down a moving object.

When these two forces fight, they create "shock waves." Depending on who is winning the fight, the wave might look like a smooth, gentle hill (monotonic) or a choppy, vibrating series of mini-waves (oscillatory). The author notes that while mathematicians have found some "perfect" formulas for these waves, those formulas often don't make sense in the real world—they describe waves that behave in ways physics wouldn't allow.

2. The "Patchwork Quilt" Solution (The Asymptotics)

Since we can't find one single, perfect formula to describe the whole wave, Samokhin uses a clever "patchwork" approach.

Think of a wave like a comet.

The Tail: Near the center (where the explosion started), the wave behaves in a very specific, curved way because of the expanding space.
The Head: Further out, the wave starts to act like a standard, flat wave moving through a straight line.

The author creates a mathematical "patchwork quilt" by taking the formula for the "tail" and stitching it to the formula for the "head." He even adds a tiny bit of "overlap" (called truncation and shifting) to make sure the seam between the two parts is invisible, creating a smooth, realistic model that matches what happens in computer simulations.

3. The "Self-Correcting" Wave (Stability)

One of the most important parts of the paper is proving that these waves are stable.

Imagine you are driving a car on a highway. If a gust of wind hits you and nudges the car slightly to the left, you steer back to the center. The car is "stable."

Samokhin uses a concept called "momentum decay" to prove that these shock waves are like that car. If a little "bump" or error hits the wave, the natural friction (dissipation) in the equation acts like a corrective steering wheel. The error eventually "fades away," and the wave returns to its original, intended shape.

4. The "Social Rules" of Waves (Superposition)

Finally, the paper asks: What happens when two waves crash into each other?

In many complex systems, when two things collide, they create a chaotic mess. But Samokhin discovers that these specific shock waves follow a surprisingly simple "social rule" called superposition.

If you have two waves, Wave A and Wave B, and they collide, they don't explode into chaos. Instead, they essentially "merge" into a new, single wave that has the combined height and speed of the original two. It’s like two small puddles merging into one larger puddle—the math stays clean and predictable.

Summary

In short, Samokhin has provided a "user manual" for expanding shock waves. He has shown how to build a realistic model of them, proved that they won't fall apart if they get bumped, and explained how they behave when they meet their neighbors.

This technical summary outlines the research presented by Alexey Samokhin regarding the asymptotic behavior, stability, and superposition of shock waves in spherical and cylindrical Korteweg-de Vries-Burgers (KdV-B) equations.

1. Problem Statement

The paper addresses the spherical and cylindrical KdV-B equations, defined by:
$u_t = -2uu_x - \frac{n}{t}u + \epsilon u_{xx} + \lambda u_{xxx}$
where $n=1$ for spherical and $n=1/2$ for cylindrical cases.

While the flat KdV-B equation ( $n=0$ ) is well-understood and possesses known exact traveling wave solutions, the spherical and cylindrical variants lack physically interpretable exact solutions. Existing exact solutions for these equations often result in velocity profiles that increase linearly or become negative as $x \to \infty$ , making them physically unrealistic. The research aims to provide a reliable asymptotic description of these waves, prove their stability, and establish rules for their interaction (superposition).

2. Methodology

The author employs a multi-faceted approach combining analytical asymptotics, dynamical systems theory, and numerical verification:

Asymptotic Analysis: The author decomposes the shock wave into two distinct regions:
1. The Near-Field (Convex Part): Near the origin ( $x \approx 0$ ), the term $\frac{n}{t}u$ dominates. The author derives analytical forms for this region using a similarity transformation $u(x,t) = y(\frac{x}{Vt})$ .
2. The Far-Field (Shock Front): As $x$ increases, the $\frac{n}{t}u$ term becomes negligible, and the solution behaves like a standard flat KdV-B traveling wave.
Correction Modeling: To prevent mathematical "breaks" or discontinuities at the junction of these two regions, the author introduces a truncation parameter ( $r^*$ ) and a small shift ( $s^*$ ). The truncation is estimated using the analytical solution of the Burgers equation.
Stability Analysis: The author utilizes the selective decay of the momentum conservation law. By analyzing the time derivative of the momentum integral, the author demonstrates that perturbations in the KdV-B system inevitably decay over time.
Numerical Simulation: The proposed asymptotic formulas are validated against direct numerical simulations to ensure visual and mathematical indistinguishability.

3. Key Contributions and Results

A. Asymptotic Approximations

The paper provides explicit piecewise formulas for the shock waves:

Cylindrical Shocks: A combination of a specific algebraic decay function for the near-field and the flat KdV-B traveling wave solution for the far-field.
Spherical Shocks: A combination involving the Lambert W function for the near-field and the flat KdV-B solution for the far-field.
These formulas successfully bridge the gap between the divergent nature of the geometry and the standard traveling wave behavior.

B. Stability Proof

The author proves that these shock waves are stable in the class of rapidly decreasing solutions. Using the momentum decay argument, it is shown that any initial perturbation $\delta(x,t)$ with zero boundary conditions will satisfy $\lim_{t \to \infty} \delta(x,t) = 0$ . This confirms that the shock wave profile is a robust attractor.

C. Superposition Rules

Based on the stability results, the author derives two practical rules for the interaction of these shocks:

Additive Superposition: The sum of two shock waves $F_{H1} + F_{H2}$ asymptotically evolves into a single shock wave with the combined height $H_1 + H_2$ .
Velocity Additivity: Because the velocity $V$ is proportional to the height $H$ in these geometries, the velocity of the resulting superposed shock is the sum of the individual velocities ( $V_{total} = V_1 + V_2$ ).

4. Significance

This work is significant because it provides a physically meaningful analytical framework for modeling diverging shock waves in non-flat geometries. By reducing the complex spherical/cylindrical problem to a corrected version of the well-studied flat KdV-B problem, the author offers a highly efficient tool for researchers. The ability to predict the behavior of superposed shocks without intensive numerical computation is particularly valuable for studying complex wave interactions in fluid dynamics and plasma physics.