🔢 mathematics

A parallel space-time $p$ -adaptive discontinuous Galerkin method for nonlinear acoustics

This paper introduces and analyzes a parallel space-time $p$ -adaptive discontinuous Galerkin method for nonlinear acoustics, establishing its well-posedness and a priori error estimates while demonstrating through numerical experiments that the approach effectively captures nonlinear phenomena like harmonic generation and reduces computational cost via adaptive refinement.

Original authors: Daniele Corallo, Pascal Lehner, Christian Wieners

Published 2026-02-16

📖 4 min read🧠 Deep dive

CC BY 4.0

Original authors: Daniele Corallo, Pascal Lehner, Christian Wieners

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 complex sound wave travels through a jelly-like substance, like the human body or the ocean. This isn't just a simple ripple; it's a nonlinear wave, meaning the wave changes shape as it moves, creating sharp peaks and strange harmonics (like a musical note suddenly turning into a chord).

This paper presents a new, super-smart way to simulate these waves on a computer. Here is the breakdown using everyday analogies:

1. The Problem: The "Too-Complicated" Equation

Traditionally, scientists used equations that were like trying to drive a car while looking in the rearview mirror and the side mirrors simultaneously. They involved "second-order" time derivatives, which made the math incredibly messy and hard to solve on a computer without the simulation crashing or giving wrong answers.

The Solution: The authors found a "shortcut." They reformulated the physics so that the equations only look at the current state and the immediate change (first-order derivatives). It's like switching from a complex, winding mountain road to a straight, flat highway. It's much easier to drive (solve) on this new road.

2. The Method: The "Space-Time LEGO" Approach

Most computer simulations work like a movie reel: they calculate one frame (a slice of time), then move to the next. This paper uses a Space-Time Discontinuous Galerkin (DG) method.

The Analogy: Imagine building a 3D model of a wave not by stacking flat 2D sheets (time steps) on top of each other, but by building a giant, solid block of LEGO bricks that represents both space and time all at once.
The "Discontinuous" Part: These LEGO bricks don't have to fit perfectly together at the edges. They can be slightly misaligned. The math has special "glue" (called interior penalty) that holds them together even if they don't match perfectly. This gives the computer the freedom to use different sized bricks in different places.

3. The "Smart" Feature: P-Adaptivity (The Chameleon)

This is the coolest part of the paper. The method is p-adaptive.

The Analogy: Imagine you are painting a landscape. In the sky, you use big, broad brushstrokes because it's smooth. But near a tree with tiny leaves, you switch to a tiny, fine brush.
How it works: The computer automatically looks at the wave. Where the wave is calm and smooth, it uses simple, low-degree math (big brushstrokes). Where the wave gets crazy, steep, or forms a shock (like a sonic boom), it automatically switches to high-degree, complex math (tiny brushstrokes) to capture the details.
The Benefit: This saves massive amounts of computing power. Instead of using a super-fine brush for the entire ocean, you only use it where the action is.

4. Solving the Puzzle: Newton's Method

Because the waves are "nonlinear" (they change based on how loud they are), the math is a giant, tangled knot. To untie it, the authors use Newton's Method.

The Analogy: Imagine you are trying to find the bottom of a dark valley. You take a step, check the slope, and adjust your path. You do this over and over, getting closer and closer to the bottom. The computer does this "guess-and-check" loop incredibly fast until it finds the perfect solution.

5. The Results: Why It Matters

The authors tested their method with two types of experiments:

The "Known Answer" Test: They created a fake wave where they already knew the answer. Their method found the answer with high precision, proving the math works.
The "Real World" Test: They simulated sound traveling through tissue. The simulation successfully reproduced real-world phenomena, like harmonic generation (where a single tone splits into multiple frequencies) and shock waves.

The Big Win: They showed that by using their "smart brush" (p-adaptivity), they could get the same accuracy as a standard method but using far fewer computer resources (degrees of freedom).

Summary

This paper introduces a new, highly efficient way to simulate complex sound waves. It treats space and time as a single 3D puzzle, uses "smart" math that automatically focuses its power where it's needed most, and solves the equations with a reliable iterative process. It's a major step forward for medical imaging, underwater sonar, and understanding how sound behaves in the real world.

1. Problem Statement

The paper addresses the numerical simulation of nonlinear acoustic wave propagation, a phenomenon critical in medical applications (e.g., tissue ablation, lithotripsy, imaging) and engineering (e.g., nondestructive testing, underwater acoustics).

Challenges: Classical models like the Westervelt or Kuznetsov equations involve second-order time derivatives acting on nonlinear terms. This structure complicates the design of stable and efficient numerical methods, particularly when resolving steep gradients or shock-like structures.
Goal: To develop a robust, high-order numerical method capable of handling the interplay between high-frequency wave propagation and nonlinear effects, while supporting adaptive refinement and parallel computation.

2. Mathematical Model

The authors utilize a first-order in time formulation derived in previous work [22], which avoids the complications of second-order time derivatives found in classical models.

Variables: $p$ (pressure fluctuation) and $v$ (velocity field fluctuation).
Governing Equations: The system is a set of coupled partial differential equations involving:
- First-order time derivatives ( $\partial_t p, \partial_t v$ ).
- Spatial Laplacians ( $\Delta p, \Delta v$ ) representing damping/viscosity.
- Nonlinear terms involving only spatial derivatives (e.g., $p\nabla \cdot v$ , $\nabla p \cdot v$ , $\nabla(v^2)$ ).
Structure: The system is treated as a symmetric hyperbolic system (Friedrichs system) with added damping and nonlinear source terms.

3. Methodology: Space-Time Discontinuous Galerkin (DG)

The core contribution is the application of a space-time Discontinuous Galerkin (DG) method to this nonlinear model.

Discretization Strategy

Space-Time Mesh: The domain is discretized using tensor-product space-time elements ( $I_j \times K$ ), allowing for independent polynomial degrees in time ( $p$ ) and space ( $q$ ).
Symmetric Friedrichs System: The hyperbolic part ( $M\partial_t + A$ ) is discretized using a symmetric Friedrichs system approach with upwind fluxes based on local Riemann problems.
Interior Penalty (IP): The Laplacian (damping) terms are handled via an interior penalty method to ensure stability and consistency across element boundaries.
Nonlinear Terms: The nonlinear volume terms are treated explicitly within the variational formulation.
Solver: The resulting nonlinear algebraic system is solved using Newton's method (Newton-Raphson), requiring the solution of a linearized system at each iteration.

Theoretical Analysis

The paper provides a rigorous mathematical analysis of the discrete problem:

Well-posedness:
- A stability analysis is first performed on the linearized system.
- Using a fixed-point argument and the Newton-Kantorovich theorem, these results are extended to the fully nonlinear discrete system.
- The analysis proves the existence of a unique solution for sufficiently small data and time steps.
Error Estimates:
- A priori error estimates are derived in a natural DG norm.
- The convergence rate depends on the polynomial degrees in space ( $q$ ) and time ( $p$ ), as well as the regularity of the exact solution.
- The analysis accounts for projection errors and the consistency of the discrete residual.

Adaptivity

p-adaptivity: The method employs a residual-based error indicator to dynamically adjust the polynomial degrees ( $p$ and $q$ ) of space-time cells.
Strategy: Cells with high error indicators are refined (degree increased), while those with low indicators are derefined (degree decreased). This aims to minimize degrees of freedom (DoFs) while maintaining accuracy for specific goal functionals.

4. Key Contributions

First Application: To the authors' knowledge, this is the first application of a space-time DG method to nonlinear acoustics.
First-Order Formulation: The work leverages a first-order time formulation that simplifies the numerical treatment of nonlinearities compared to second-order models.
Rigorous Convergence Proof: The paper establishes a complete well-posedness and convergence analysis for the nonlinear space-time DG scheme, including a priori error estimates.
Parallel Implementation: The method is implemented in the M++ framework, demonstrating parallel solvability and efficiency.
p-Adaptive Framework: The introduction of a residual-based $p$ -adaptive strategy specifically tailored for space-time nonlinear acoustic problems.

5. Numerical Results

The authors conducted numerical experiments in 2D to validate the theory:

Known Solution Test: Using an analytical solution, the method demonstrated convergence rates that generally exceeded theoretical predictions in the DG norm.
Nonlinear Phenomena: The model successfully reproduced characteristic nonlinear acoustic effects, specifically:
- Harmonic Generation: The emergence of frequencies at integer multiples of the base frequency.
- Waveform Distortion: The shift in peak values (steepening) due to local changes in sound speed, distinguishing it from linear wave propagation.
Adaptivity Performance:
- The error estimator effectively guided refinement.
- Goal Functionals: Adaptive refinement significantly reduced the DoFs required to accurately approximate specific functionals (e.g., total pressure or energy) compared to uniform refinement, particularly at high resolutions (around $10^7$ DoFs).
- The adaptive method showed better convergence for energy functionals ( $E$ ) than for integrated pressure ( $P$ ).

6. Significance and Conclusion

Robustness: The proposed method offers a stable and flexible framework for simulating complex nonlinear acoustic waves, handling shocks and high frequencies effectively.
Efficiency: The space-time formulation combined with $p$ -adaptivity allows for efficient resource allocation, reducing the computational cost for high-accuracy simulations.
Future Work: The authors note that a rigorous convergence analysis for the adaptive algorithm itself remains an open question. Future directions include exploring $hp$-adaptivity and hybrid DG methods to further enhance efficiency.

In summary, this paper bridges the gap between advanced space-time numerical methods and the complex physics of nonlinear acoustics, providing both a theoretical foundation and a practical, parallelizable solver for high-fidelity simulations.