🔢 mathematics

On the Jordan-Moore-Gibson-Thompson equation of nonlinear acoustics

This paper provides a systematic overview of the Jordan-Moore-Gibson-Thompson (JMGT) equation in nonlinear acoustics, focusing on its well-posedness, memory and fractional attenuation effects, singular limits, and illustrative examples of control and inverse problems.

Original authors: Barbara Kaltenbacher

Published 2026-04-09

📖 6 min read🧠 Deep dive

CC BY 4.0

Original authors: Barbara Kaltenbacher

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

The Big Picture: Fixing a "Ghost" in the Machine

Imagine you are trying to predict how sound waves travel through water or air, especially when those waves are very loud (like a sonic boom or a medical ultrasound).

For a long time, scientists used "classic" math models to do this. But these models had a weird, impossible flaw: they predicted that sound travels infinitely fast. It's like saying if you snap your fingers in New York, someone in London hears it instantly, before the sound wave even leaves your finger. In the real world, nothing travels faster than light, and sound is definitely slower.

This paper is about a "fix" for that problem. The author, Barbara Kaltenbacher, reviews a specific mathematical equation called the Jordan-Moore-Gibson-Thompson (JMGT) equation. Think of this equation as a new, more realistic rulebook for how sound behaves when it gets loud and fast.

The Main Characters: The Three Equations

The paper discusses a few versions of this rulebook, depending on how complicated the sound is:

The "Westervelt" Version: Good for tracking the pressure of the sound (how hard the air is being squeezed).
The "Kuznetsov" Version: Good for tracking the movement of the air particles themselves.
The "Linear" Version: A simplified version used when the sound is quiet and behaves nicely.

The paper focuses mostly on the first two because they handle the "messy" parts of loud sound, like when waves crash into each other and create new, weird frequencies.

The Secret Ingredient: The "Time Lag" ( $\tau$ )

The magic ingredient in this new equation is a tiny number called $\tau$ (tau).

The Old Way: Imagine a relay race where the baton is passed instantly. If Runner A moves, Runner B moves immediately. This is the old model, and it causes the "infinite speed" paradox.
The JMGT Way: Imagine Runner A has to take a tiny breath before passing the baton to Runner B. There is a tiny delay or lag.

In physics, this lag represents thermal relaxation. When sound waves heat up the air, the air doesn't cool down instantly; it takes a split second. The JMGT equation includes this split-second delay. This tiny delay fixes the "infinite speed" problem and makes the math match reality.

What the Paper Actually Does

Kaltenbacher isn't just saying "this equation is cool." She is acting like a rigorous engineer checking if the blueprint works. She tackles four main questions:

1. Does the Math Make Sense? (Well-Posedness)

Before you build a bridge, you need to know if the math says it will stand up.

The Question: If we start with a specific sound wave, does the equation give us one clear, stable answer? Or does it explode into nonsense?
The Finding: Yes! As long as the sound isn't too loud (small data), the equation works perfectly. It predicts a stable wave. However, if the sound is incredibly loud (large data), the wave can "blow up" (mathematically speaking, the pressure goes to infinity in a finite time), which in real life would mean the wave shatters or creates a shockwave.

2. Does the Sound Fade Away? (Decay)

The Question: If you stop making noise, does the sound eventually die out, or does it bounce around forever?
The Finding: Because of the "friction" in the air (viscosity) and the thermal lag, the sound energy dissipates. The paper proves that under the right conditions, the sound waves will eventually calm down and disappear, just like ripples in a pond.

3. What Happens When the Delay Vanishes? (Singular Limits)

The Question: What if we set that tiny delay ( $\tau$ ) to zero? Do we get the old, flawed models back?
The Finding: Yes. The paper shows that as you shrink the delay to zero, the JMGT equation smoothly transforms back into the older, classical equations. This proves that the new equation is a "parent" that contains the old ones as special cases. It's like showing that a high-definition video is just a regular video with more pixels.

4. Can We Control or "See" Inside? (Control & Inverse Problems)

This is where it gets practical for medicine and engineering.

Control: Imagine you are a doctor using ultrasound to heat a tumor. You want the sound to hit exactly the tumor and nowhere else. The paper discusses how to mathematically calculate the perfect "push" (control) on the speakers to get the sound to the right spot.
Inverse Problems (The Superpower): This is the coolest part. Usually, if you listen to a sound coming out of a box, it's hard to guess what's inside. But because the JMGT equation handles nonlinearity (loud sounds creating new frequencies), it's easier to figure out what's inside.
- Analogy: Imagine tapping a drum. If the drum is made of wood, it makes a "thud." If it's metal, it makes a "ping." If you hit it really hard, the wood might make a weird "crack" sound that the metal doesn't.
- The paper shows that by listening to these "weird sounds" (higher harmonics) created by the nonlinearity, we can uniquely identify the properties of the tissue (like how fast sound travels through it) much better than with the old, linear models.

The Takeaway

This paper is a massive "user manual" for the Jordan-Moore-Gibson-Thompson equation.

It tells us:

Why we need it: To fix the "infinite speed" error in old models.
That it works: The math is stable and predictable for realistic scenarios.
How it connects: It links perfectly to the old models when the delay is tiny.
Why it matters: It helps us design better ultrasound machines, control sound waves more precisely, and "see" inside the human body with greater clarity by using the unique properties of loud sound.

In short, Kaltenbacher has taken a complex, third-order differential equation and proven that it is the robust, reliable, and physically accurate tool we need to understand the loud, messy, and fascinating world of nonlinear acoustics.

1. Problem Statement and Context

The paper addresses the Jordan-Moore-Gibson-Thompson (JMGT) equation, a third-order in time quasilinear partial differential equation (PDE) used to model finite-amplitude sound propagation in gases and liquids.

Physical Motivation: Classical models of nonlinear acoustics, such as the Westervelt and Kuznetsov equations, are second-order in time and rely on the Fourier heat flux law. This leads to the "paradox of infinite speed of propagation." The JMGT equation resolves this by incorporating the Maxwell–Cattaneo law for heat flux, introducing a time lag parameter $\tau > 0$ (thermal relaxation). This restores the hyperbolic nature of the model, ensuring finite propagation speeds.
Mathematical Challenge: The JMGT equation is a third-order hyperbolic PDE. Its analysis is more complex than classical second-order models due to the presence of the third time derivative ( $\tau u_{ttt}$ $τ u_{ttt}$ ). The paper focuses on two primary nonlinear variants:
1. JMGT-Westervelt: For acoustic pressure $u$ , with nonlinearity $(u^2)_{tt}$ .
2. JMGT-Kuznetsov: For acoustic velocity potential $u$ , with nonlinearity $(\tilde{\eta}u_t^2 + |\nabla u|^2)_t$ .
Scope: The paper provides a systematic overview of the mathematical literature, focusing on well-posedness (local and global), singular limits ( $\tau \to 0$ ), memory/fractional attenuation, periodic solutions, and applications in control and inverse problems.

2. Methodology

The author employs a unified analytical framework based on energy estimates and fixed-point arguments.

Energy Identity Construction: A central technique involves defining a combined quantity $z = \tau u_t + u$ . This transforms the third-order equation into a perturbed wave equation. By testing the equation with a specific function close to $-\Delta(\tau u_{tt} + \sigma u_t + \rho u)$ , the author derives a comprehensive energy identity (Equation 20).
Stability Criterion: The analysis relies heavily on the condition $\delta = b - \tau c^2 > 0$ , where $b$ is the damping coefficient, $c$ is the sound speed, and $\tau$ is the relaxation time. This condition ensures exponential stability and is linked to the Routh-Hurwitz criterion for the linearized system.
Fixed-Point Theorems: For the nonlinear problems, the author uses Banach's Contraction Principle. The nonlinearity is treated entirely as a source term on the right-hand side of the linearized equation, simplifying the exposition compared to previous literature that incorporated parts of the nonlinearity into the highest-order derivative.
Singular Limit Analysis: To connect JMGT to classical models, the paper derives $\tau$ -uniform energy estimates. This allows for the rigorous proof of convergence as $\tau \to 0$ , where the hyperbolic JMGT equation converges to the parabolic Westervelt/Kuznetsov equations.
Inverse Problem Linearization: For imaging applications, the author linearizes the operator equation relating coefficients to observations and utilizes the Inverse Function Theorem, leveraging the frequency-domain structure of the solution.

3. Key Contributions and Results

A. Well-Posedness Analysis

Local Well-Posedness: Proven for both bounded domains and $\mathbb{R}^d$ with small initial data. The proof relies on the energy estimates derived from the unified identity.
Global Well-Posedness & Exponential Decay: Established in the dissipative case ( $\delta > 0$ ) for small initial data. The author proves that solutions exist globally and decay exponentially.
Critical Case ( $\delta = 0$ ): In the inviscid limit, the system is marginally stable. Global well-posedness can still be achieved via absorbing boundary conditions (feedback control).
Blow-up: For large initial data, the paper discusses results showing finite-time blow-up (non-existence of global solutions), specifically when the nonlinearity grows sufficiently fast.

B. Memory and Fractional Attenuation

Memory Terms: The paper reviews models where the standard damping is replaced by convolution integrals (memory kernels). It highlights that Type III memory (involving $z = \tau u_t + u$ ) can induce exponential decay even in the critical case $\delta = 0$ , whereas Type I memory generally requires $\delta > 0$ .
Fractional Attenuation: The paper discusses extensions using fractional time derivatives (Djirbashian-Caputo) to model frequency-dependent attenuation more accurately. These models arise from fractional generalizations of the Maxwell-Cattaneo law and are motivated by thermodynamic consistency.

C. Time Periodic Solutions

The paper establishes the existence of time-periodic solutions (relevant for continuous wave ultrasound).
It demonstrates that the nonlinearity generates higher harmonics (multiples of the fundamental frequency).
The problem is reformulated in the frequency domain as a coupled system of Helmholtz-type equations, facilitating the analysis of harmonic generation.

D. Singular Limits ( $\tau \to 0$ )

The author provides a rigorous convergence analysis showing that as $\tau \to 0$ , the solutions of the JMGT equation converge to those of the classical Westervelt and Kuznetsov equations.
A key technical achievement is the derivation of $\tau$ -independent energy estimates, which are necessary to handle the change in the equation's type from hyperbolic to parabolic.

E. Control and Inverse Problems

Optimal Control: The paper reviews a control problem where boundary controls (Neumann conditions) are used to track a desired pressure profile. It highlights the challenge of unbounded control operators in hyperbolic systems and the use of Riccati equations for feedback synthesis.
Inverse Problems (Imaging): The paper argues that nonlinearity aids uniqueness in reconstructing spatially varying coefficients (sound speed $c$ $c$ , attenuation $b$ $b$ , nonlinearity $\eta$ $η$ ).
- Unlike linear models which often require infinite sources for uniqueness, the nonlinear JMGT model allows for the unique local reconstruction of the nonlinearity coefficient $\eta(x)$ using a single source and time-domain observations.
- The relaxation time $\tau$ is crucial here; it restores finite propagation speed, preventing the "smearing" of information caused by strong attenuation in parabolic models, thereby improving the stability of the inverse problem.

4. Significance

This paper serves as a comprehensive bridge between the physical modeling of nonlinear acoustics and rigorous mathematical analysis. Its significance lies in:

Unification: It offers a unified energy-based framework that handles initial value problems, periodic problems, and singular limits simultaneously.
Physical Realism: It validates the JMGT equation as a superior physical model over classical second-order equations by resolving the infinite speed paradox and providing a mathematically sound basis for finite propagation speeds.
Inverse Problem Advancement: It demonstrates that the hyperbolic nature of the JMGT equation, combined with nonlinear effects, significantly enhances the solvability and uniqueness of inverse problems in medical ultrasound imaging, potentially leading to better diagnostic tools.
Literature Synthesis: It consolidates a vast body of recent mathematical work (including memory effects, fractional derivatives, and blow-up analysis) into a coherent narrative, highlighting open questions and the interplay between physical parameters and mathematical stability.