Conservation laws and exact solutions of a nonlinear acoustics equation by classical symmetry reduction

This paper investigates the point symmetries, conservation laws, and exact solutions of the generalized Westervelt equation—a nonlinear model for sound wave propagation in compressible media—utilizing modern symmetry reduction methods to derive potential symmetries and analyze shock wave formation.

The Gist

Imagine you are trying to understand how a sound wave travels through a crowd of people. If the crowd is sparse and calm, the wave moves smoothly, like a gentle ripple in a pond. But if the crowd is dense and the sound is loud (like a shout in a stadium), the people bump into each other, the wave gets squashed, and it behaves in a messy, unpredictable way. This is nonlinear acoustics.

The paper you're asking about is a mathematical detective story. The authors are trying to solve a very complicated equation (the Westervelt equation) that describes exactly this messy behavior of sound in things like human tissue (for ultrasound) or underwater.

Here is the breakdown of their investigation, explained simply:

1. The Problem: A Messy Equation

The equation they are studying is like a recipe for sound waves, but it has a "secret ingredient" (nonlinearity) that makes it very hard to solve. Usually, when you try to predict how sound moves, you get a smooth curve. But in real life, sound can suddenly "break" or form a shock wave (like a sonic boom or a sudden snap). The authors wanted to find the hidden rules that govern this equation so they could predict exactly what happens.

2. The Detective Tool: "Symmetry"

Think of Symmetry as a magic mirror. If you look at a snowflake in the mirror, it looks the same. In math, if you shift a sound wave forward in time or slide it to the left, and the equation still works perfectly, that's a symmetry.

The authors used a special kind of math (Lie Group theory) to find all the "magic mirrors" for their sound equation.

• Time Translation: The rules of sound don't change just because it's Tuesday instead of Monday.
• Space Translation: The rules don't change if you move the experiment from Cadiz to New York.
• Scaling: They found that for certain types of sound, if you stretch the wave or shrink it, the math still holds up.

By finding these symmetries, they could simplify the giant, scary equation into much smaller, manageable pieces.

3. The Hidden Treasure: "Conservation Laws"

In physics, Conservation Laws are like an unbreakable bank account. No matter what happens, the total amount of "stuff" (like energy or mass) stays the same; it just moves around.

The authors discovered that this sound equation has a hidden bank account: The Net Mass of the Sound Wave.

• Imagine the sound wave is a cloud of air. Even as the cloud squishes and stretches, the total "weight" of that cloud (adjusted for the weird physics of the equation) never disappears.
• They used a modern method (the "multiplier method") to find these rules. It's like finding the receipt that proves you didn't lose any money, even though you spent it in weird ways.

4. The "Ghost" Variables: Potential Systems

Sometimes, the rules of the game are hidden. To see them, the authors invented "ghost variables" (called potentials).

• Imagine you are looking at a shadow on a wall. You can't see the object casting the shadow, but you can guess its shape.
• The authors created a "shadow world" (a potential system) where the equations look different. In this shadow world, they found new symmetries and conservation laws that didn't exist in the original "real world" equation.
• These are called nonlocal laws. Think of it like knowing that if you push the left side of a giant rubber sheet, the right side moves instantly, even though they aren't touching. It's a connection that spans across space.

5. The Grand Finale: The Shock Wave

After all this detective work, they used their findings to solve the equation for a specific, real-world scenario: Shock Waves.

• The Analogy: Imagine a line of cars on a highway. If everyone drives at the same speed, traffic flows smoothly. But if the car in front slams on the brakes, the cars behind pile up instantly. That pile-up is a "shock."
• The Result: The authors found an exact formula for how this "traffic jam" of sound looks. They showed that the wave starts smooth, then suddenly steepens into a sharp cliff (the shock), and then levels out again.
• They even drew pictures (Figures 1 and 2 in the paper) showing what this looks like in 3D space and time. It looks like a sudden, sharp wall of sound moving through the air.

Why Does This Matter?

You might ask, "Who cares about a math equation for sound?"

• Medical Imaging: When doctors use ultrasound to look inside your body, the sound waves can get distorted by your tissues. Understanding these "shock waves" helps doctors get clearer pictures and avoid hurting patients.
• Underwater Tech: Sonar systems need to know how sound behaves in the ocean, where pressure changes everything.
• Music: Even the sound of a trumpet involves these nonlinear effects!

Summary

The authors took a messy, complicated equation describing loud, distorted sound. They used symmetry (finding patterns) and conservation laws (finding what stays the same) to simplify it. They even looked at the problem through a "ghost lens" (potentials) to find hidden rules. Finally, they used these tools to write down the exact formula for a shock wave, helping us understand how sound behaves when it gets too loud to be smooth.

It's a bit like taking a chaotic jazz improvisation, finding the underlying rhythm, and writing down the sheet music so anyone can play it perfectly.

Technical Summary

1. Problem Statement

The paper investigates a generalized Westervelt equation, a nonlinear partial differential equation (PDE) modeling sound wave propagation in a compressible medium. This equation is crucial for applications in biomedical ultrasound, underwater imaging, and sonochemistry.

The specific equation studied is:
$$f(p)_{tt} - \beta p_{xxt} = c^2 p_{xx}$$
where:

• $p(t, x)$ is the pressure fluctuation.
• $\beta > 0$ is the damping coefficient.
• $c$ is the speed of sound.
• $f(p)$ is an arbitrary function representing the nonlinearity. The physical case of interest is when $f(p) = p + h(p)$, where $h(p)$ is a nonlinear function (e.g., quadratic for the original Westervelt equation).

Key Challenges:

• The equation lacks a local Lagrangian formulation in terms of $p$, rendering the classical Noether's theorem inapplicable for deriving conservation laws directly.
• There is a need to classify symmetries and conservation laws for the general function $f(p)$, not just specific cases, to understand the equation's integrability and solution structure.
• Finding exact solutions, particularly those describing physical phenomena like shock waves, is essential for modeling real-world acoustic propagation.

2. Methodology

The authors employ Lie symmetry analysis and the multiplier method within the framework of jet space and variational calculus.

• Point Symmetry Classification: They determine the infinitesimal generators of Lie point symmetries by solving the symmetry determining equations. This involves classifying the admitted symmetries based on the form of the arbitrary function $f(p)$.
• Multiplier Method: To find conservation laws without a Lagrangian, they use the modern generalization of Noether's theorem. They identify multipliers $Q$ such that the product of the equation and $Q$ yields a total divergence ($D_t T + D_x \Phi = 0$).
• Potential Systems: To uncover nonlocal symmetries and conservation laws not present in the original equation, they introduce potential variables. They construct two distinct potential systems:
  1. A first-layer system involving a potential $u$.
  2. A second-layer system involving potentials $v$ and $w$.
     They analyze the symmetries and conservation laws of these extended systems and project them back to the original PDE.
• Symmetry Reduction: They apply the identified point symmetries to reduce the PDE to Ordinary Differential Equations (ODEs). Specifically, they focus on translation symmetries to derive travelling wave solutions.

3. Key Contributions and Results

A. Classification of Point Symmetries

The paper provides a complete classification of point symmetries for the generalized equation with $f''(p) \neq 0$ and $\beta \neq 0$:

• General Case: Admits only time-translation ($X_1 = \partial_t$) and space-translation ($X_2 = \partial_x$).
• Specific Nonlinearities:
  • If $f(p) = k(p+p_0)^{1+q}$, an additional scaling symmetry exists ($X_3$).
  • If $f(p) = k(p+p_0)^{-3}$, an additional non-rigid scaling symmetry exists ($X_4$).
• The authors explicitly derive the commutator structure of these symmetry generators, showing how they form a Lie algebra.

B. Local Conservation Laws (Low-Order)

Using the multiplier method, four low-order multipliers ($Q_1=1, Q_2=t, Q_3=x, Q_4=tx$) are identified, leading to four conservation laws.

• Physical Interpretation: The conserved quantities derived from these laws are related to the net mass of the sound waves.
• By relating pressure $p$ to density $\rho$ via the equation of state, the conserved integrals are shown to represent the generalized net mass and the generalized $x$-weighted net mass.

C. Potential Systems and Nonlocal Laws

The authors derive two potential systems to find symmetries and conservation laws that are nonlocal (not inherited from the original PDE).

• First Potential System: Introduces potential $u$.
  • Yields a new potential symmetry (a $u$-translation) not present in the original PDE.
  • Yields nonlocal conservation laws involving integrals of $f(p)$.
• Second Potential System: Introduces potentials $v$ and $w$.
  • Yields additional potential symmetries (translations and shifts in $v$ and $w$).
  • Identifies a specific nonlocal conservation law involving the potential $v$.
• These results demonstrate that the equation admits structures (symmetries and conservation laws) that are invisible when analyzing the PDE in its original form.

D. Exact Travelling Wave and Shock Solutions

Applying the translation symmetry ($p = U(x - \nu t)$), the PDE is reduced to a third-order ODE, which integrates to a first-order ODE.

• General Solution: For a general nonlinear function $h(p)$, the solution is expressed as a quadrature (an integral form).
• Shock Wave Solutions: By assuming a power-law nonlinearity $h(U) = \kappa U^n$ (with $n > 1$), the authors derive explicit shock wave solutions.
  • For the specific physical case $n=2$ (original Westervelt), the solution is:
    $$U(\xi) = \frac{c^2 - \nu^2}{U_0(c^2 - \nu^2)e^{-(c^2-\nu^2)\xi/\beta\nu} + \kappa\nu^2}$$
  • Behavior: These solutions represent shock waves characterized by a sharp transition between two asymptotic values. Unlike solitons, these waves dissipate energy. The paper confirms that for $c^2 > \nu^2$, the solution is non-singular and exhibits the expected shock profile.

4. Significance

• Mathematical Rigor: The paper fills a gap in the literature by providing a complete symmetry and conservation law classification for the generalized Westervelt equation, rather than just specific cases.
• Methodological Demonstration: It effectively demonstrates the utility of the multiplier method for non-variational PDEs and the power of potential systems in uncovering hidden (nonlocal) structures.
• Physical Relevance: The derivation of explicit shock wave solutions provides a theoretical basis for understanding acoustic shock formation in biomedical and industrial applications (e.g., ultrasound imaging and sonochemistry).
• Foundation for Future Work: The classification of symmetries and the identification of conservation laws provide a robust framework for future studies, including numerical accuracy checks, integrability analysis, and extensions to 2D/3D geometries.

In summary, this work establishes the fundamental mathematical structure of the generalized Westervelt equation, linking abstract symmetry analysis to concrete physical phenomena like mass conservation and shock wave propagation.