Khokhlov–Zabolotskaya–Kuznetsov Equation

Updated 17 December 2025

KZK equation is a fundamental nonlinear PDE that models unidirectional sound wave propagation by incorporating quadratic nonlinearity, transverse diffraction, and dissipative effects.
It is derived via paraxial reduction of the isentropic compressible Navier–Stokes or Euler equations using a multiple-scale expansion that balances weak nonlinearity with viscosity and dispersion.
The equation underpins advanced modeling in high-intensity focused ultrasound, laser pulses, and nonlinear elastic wave phenomena while exhibiting rich symmetry and conservation properties.

The Khokhlov–Zabolotskaya–Kuznetsov (KZK) equation is a fundamental model in nonlinear acoustics and wave propagation theory, describing the unidirectional evolution of weakly nonlinear, weakly diffracting, and absorbing sound beams in isotropic, nearly homogeneous media. It arises as a leading-order paraxial approximation of the isentropic compressible Navier–Stokes or Euler equations, combining quadratic nonlinearity, transverse diffraction, and dissipative or dispersive effects in a single, tractable partial differential equation. Mathematically, the KZK equation encapsulates a canonical class of (2+1)-dimensional evolution equations with nontrivial symmetry, reduction, and conservation properties, and serves as both a paradigm for theoretical analysis (integrability, symmetry classification) and as the standard tool for modeling intense ultrasound beams, laser pulses, or finite-amplitude shear waves in nonlinearly elastic solids.

1. Canonical Formulation and Physical Derivation

In Cartesian coordinates, for acoustic pressure $p(x,y,z,t)$ or a velocity-related variable $I(\tau, z, y)$ (with $\tau$ retarded time), the KZK equation typically takes the form: $c\,\frac{\partial^2 I}{\partial \tau \,\partial z} \;-\; \frac{\gamma+1}{4\rho_0} \frac{\partial^2}{\partial \tau^2}(I^2) \;-\;\frac{\nu}{2c^2\rho_0} \frac{\partial^3 I}{\partial \tau^3} \;-\;\frac{c^2}{2} \Delta_y I = 0,$ where $c$ is the sound speed, $\gamma$ is the adiabatic index, $\nu$ is the kinematic viscosity, and $\rho_0$ is the equilibrium density. The equation is derived via the following hierarchy (Dekkers et al., 2020, Dekkers et al., 2018, Destrade et al., 2013):

Navier–Stokes or Euler base: Start with the isentropic compressible equations, apply irrotationality, and expand about a quiescent state.
Small-amplitude (Mach-number) expansion: Weak nonlinearity parameter $\varepsilon \ll 1$ organizes amplitude and constitutive expansions.
Paraxial change of variables: Introduce retarded time and slow spatial variables aligned with beam direction, enabling reduction to a first-order evolution PDE in $z$ (propagation axis).
Order matching: Balance propagation, nonlinearity, viscosity, and diffraction at the same asymptotic order ( $O(\varepsilon)$ ), yielding the paraxial KZK equation.

For elastic solids, an analogous hierarchy holds starting from nonlinear elasticity with specialized constitutive assumptions (generalized neo-Hookean, etc.), leading to the scalar KZK equation for anti-plane shear waves (Destrade et al., 2013).

2. Structure, Symmetries, and Generalizations

The general KZK class (gKZK) admits further generalization with variable coefficients: $\bigl(u_t + p(t)\,u\,u_x + q(t)\,u_{xx}\bigr)_x + \sigma(t)\,u_{yy} = 0,$ where $p(t), q(t), \sigma(t)$ are nonzero smooth functions (Gungor et al., 2014). Via equivalence transformations, two of these coefficients can be normalized, yielding the canonical representative: $\bigl(u_t + u\,u_x + u_{xx}\bigr)_x + \sigma(t)\,u_{yy} = 0.$ The KZK equation admits an infinite-dimensional nilpotent Lie symmetry algebra generated by vector fields

$X(f) = f(t)\,\partial_x + \dot f(t)\,\partial_u,\quad Y(g) = g(t)\,\partial_y - \frac{\dot g(t)}{2\sigma(t)} y \partial_x - \frac{\ddot g(t)}{2\sigma(t)} y \partial_u,$

where $f, g$ are arbitrary smooth functions, and additional (finite-dimensional) symmetries arise under special choices of $\sigma(t)$ related to Möbius invariants of $t$ . For subclasses such as the generalized dispersionless KP (gdKP) or generalized KP (gKP), integrability and the existence of a full Kac–Moody–Virasoro (KMV) algebra are characterized by explicit coefficient relations, e.g., $p(t)=(c_1\int^t \sigma(s)\,ds+c_2)^{-3/2}\sigma(t)$ for gdKP (Gungor et al., 2014).

3. Paraxial Reduction and Relationship to Underlying Physics

The KZK equation systematically arises as a reduced model capturing the essential balance of acoustic nonlinearities, transverse diffraction, and small absorption in the so-called paraxial regime. The reduction process involves the following (Dekkers et al., 2020, Dekkers et al., 2018):

Irrotational, isentropic flows: Start from the full set of compressible equations with weak viscosity.
Multiple-scale ansatz: Seek solutions concentrated around a dominant direction (quasi-plane wave) with slow evolution in both axial and transverse coordinates.
Nonlinear, diffractive, dissipative effects: Retain propagation, quadratic nonlinearity, and transverse Laplacian at leading $\varepsilon$ order.
Hierarchical validity: The KZK equation is justified for times/distances on $O(1)$ scales in the slow (axial) variable, with the solution remaining $O(\varepsilon)$ close (in $L^2$ norm) to that of the parent Kuznetsov or Navier–Stokes/Euler system.

These features underpin KZK's role in modeling HIFU, lithotripsy, nonlinear ocean acoustics, and related beam phenomena (Dekkers et al., 2020, Dekkers et al., 2018).

4. Lie Symmetry Classification and Integrable Subclasses

Comprehensive Lie symmetry analysis reveals both an infinite-dimensional nilpotent structure and finite algebraic extensions for special coefficient cases. Key findings (Gungor et al., 2014):

Infinite-dimensional ideal: For arbitrary $\sigma(t)$ , symmetry generators $X(f), Y(g)$ form an ideal via commutator $[Y(g_1),Y(g_2)]=X(g_1\dot g_2-g_2\dot g_1)$ .
Finite-dimensional enlargements: For $\sigma(t)=\sigma_0$ , additional generators $\partial_t$ and $2t\partial_t + x\partial_x + y\partial_y - u\partial_u$ appear, completing a solvable algebra.
Reduction to (1+1) Burgers-type PDE: Application of symmetry reduction yields an ansatz reducing gKZK to a Burgers equation with variable coefficients; further change of variables removes residual time dependence, mapping to the classical Burgers equation.
Integrability and Virasoro symmetry: For the gdKP equation, integrability (and KMV symmetry) holds if and only if the coefficients satisfy explicit algebraic relations, mirroring the transformation theory to constant-coefficient dKP or KP equations.

This formalism underpins the mathematical structure and links to integrable systems within the KZK class.

5. Conservation Laws in Standard Geometries

The KZK equation admits an infinite family of local conservation laws in Cartesian, cylindrical, and spherical geometries (Hereman et al., 15 Dec 2025). Using the multiplier (characteristic) method, conserved currents are constructed by multipliers of the form $\Lambda=\phi(x,y,z)+t\psi(x,y,z)$ (in Cartesian coordinates), where $\psi$ is harmonic and $\phi$ solves a Poisson equation driven by $\psi_z$ . The corresponding conserved densities and fluxes are: $T^t = (\phi + t \psi)\big(\delta\,p_{tt} + 2\,\tilde\beta\,p\,p_t - 2c_0^3\,p_z\big) - \psi\big(\delta\,p_t + \tilde\beta\,p^2\big),$ with fluxes in $x$ , $y$ , and $z$ similarly constructed. In cylindrical or spherical symmetry, the transverse Laplacian and associated equations are modified accordingly.

Physical interpretation: Each conservation law encodes a generalized energy-momentum invariant, with the freedom to select harmonic $\psi$ yielding infinitely many moment-type invariants. These are instrumental for analytic estimates, existence and uniqueness proofs, and as diagnostics in numerical computation.

6. Well-posedness, Error Estimates, and Range of Applicability

The KZK initial-boundary value problem (e.g., in a half-space $z\ge0$ with periodic-in-time, small-amplitude data of zero mean) admits global well-posedness for $I_0$ in appropriate Sobolev spaces, provided the viscosity parameter is positive (Dekkers et al., 2020, Dekkers et al., 2018). Key mathematical results:

For sufficiently smooth, small initial data, there exists a unique global solution $I$ that remains periodic in retarded time $\tau$ and of zero mean for each $z$ .
An $L^2$ -uniform error estimate between the KZK solution and the solution to the underlying Kuznetsov (or Navier–Stokes/Euler) system holds for $z=O(1)$ , with $O(\varepsilon)$ accuracy, provided initial and boundary data are appropriately matched.
In the absence of viscosity, shock formation is possible in finite propagation distance, with analyticity breaking at the (nonlinear) gradient blow-up.

Boundary conditions and initial data are prescribed based on the physical directionality of the beam (inflow at $x_1=0$ ), enforced according to the KZK-generated velocity field.

The KZK equation is linked to a hierarchy of related nonlinear evolution equations:

Kuznetsov equation: The originating model for weakly nonlinear acoustics from which the KZK arises via paraxial reduction (Dekkers et al., 2020, Dekkers et al., 2018).
Nonlinear Progressive wave Equation (NPE): Alternative paraxial scaling leads to NPE, differing from the KZK principally in the sign of the diffraction term; the analysis of well-posedness and decay carries over via a bijective change of variables.
Burgers, Zabolotskaya–Khokhlov (ZK), KP, and related limits: Certain asymptotic limits (elimination of dispersion, diffraction, or nonlinearity) or further dimensional reduction yield well-known one-dimensional or lower-dimensional model equations, such as Burgers or the Zakharov–Kuznetsov (ZK) equation (Destrade et al., 2013).

These connections elucidate the regime of KZK validity and its interpretative placement within the broader context of nonlinear wave propagation.

References: