---
title: Khokhlov–Zabolotskaya–Kuznetsov Equation
url: https://www.emergentmind.com/topics/khokhlov-zabolotskaya-kuznetsov-equation
type: topic
---

# Khokhlov–Zabolotskaya–Kuznetsov Equation

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 ([2004.04657], [1811.10850], [1302.0109]):

- **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 ([1302.0109]).

## 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 ([1402.1941]). 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 ([1402.1941]).

## 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 ([2004.04657], [1811.10850]):

- *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 ([2004.04657], [1811.10850]).

## 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 ([1402.1941]):

- **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 ([2512.13062]). 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 ([2004.04657], [1811.10850]). 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.

## 7. Connections, Reductions, and Related Models

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 ([2004.04657], [1811.10850]).
- **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 ([1302.0109]).

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

---

**References**:  
- Lie symmetries and integrable reductions: [1402.1941]  
- Asymptotic and physical derivation in elasticity: [1302.0109]  
- Conservation laws in standard geometries: [2512.13062]  
- Well-posedness and error analysis: [2004.04657], [1811.10850]

Source: https://www.emergentmind.com/topics/khokhlov-zabolotskaya-kuznetsov-equation