Third-order Conservation Law

• Third-order conservation laws are invariance properties that maintain a constant third-order functional—often involving higher derivatives or differences—across evolving systems.
• They are derived using methods such as variational principles, symmetry analysis, and Fourier transforms, yielding explicit forms like mixed space-time energy functionals.
• These laws are pivotal for ensuring a priori bounds, enhancing global existence theories, and improving numerical scheme accuracy in applications from turbulence to discrete integrable systems.

A third-order conservation law refers to an invariance property or structure that guarantees the constancy of a third-order functional (often involving third derivatives or differences, or third-order correlation functions) under evolution governed by a dynamical equation, discrete system, or turbulent regime. Such laws extend classical conservation principles (mass, energy, momentum) to higher-order functionals and emerge in diverse contexts including nonlinear PDEs, discrete integrable systems, turbulent transport, numerical schemes, and generalized evolution equations.

1. Foundational Examples of Third-order Conservation Laws

Third-order conservation laws have been rigorously constructed and analyzed in several canonical equations and systems. The Kirchhoff–Pokhozhaev equation in the work of Boiti–Manfrin establishes an explicit third-order conserved functional $I_3$ for the nonlinear wave equation $u_{tt} - (a\|\nabla u\|^2 + b)^{-2} \Delta u = 0$, defined for $u$ on a domain $\Omega \subset \mathbb{R}^n$ with Dirichlet boundary conditions. Under regularity and nondegeneracy conditions, $I_3$ incorporates mixed space-time derivatives, cross terms, and higher moments, ensuring $dI_3/dt = 0$ (Boiti et al., 13 Dec 2025). For difference equations, first integrals (discrete conservation laws) are constructed via characteristic invariants, as in $u_{n+3} = w(n, u_n, u_{n+1}, u_{n+2})$ where the conserved function $\Phi(n, u_n, u_{n+1}, u_{n+2})$ satisfies a shift-invariance condition (Mamba et al., 2017). In turbulence, the third-order moments of MHD Elsässer fields satisfy an exact vectorial divergence relation for energy flux across scales (Galtier, 2010).

2. Mathematical Construction and Derivation Methods

Third-order conservation laws are established through variational or symmetry methods, Fourier/Liouville transform techniques, or combinatorial invariants. For PDEs, approaches include:

• Construction of quadratic forms in Fourier space and solution of associated linear ODEs for coefficients, as in the Kirchhoff–Pokhozhaev case.
• Self-adjointness analysis enables the application of Ibragimov’s theorem, producing conserved vectors for self-adjoint, quasilinear evolution equations of up to third order, even in the absence of classical Lagrangian structures (Freire, 2010).
• For discrete equations, Lie symmetry generators and characteristic flows yield invariants; determining equations for multipliers $P_2$ are solved to construct conserved functionals (Mamba et al., 2017).

In MHD turbulence, third-order divergence relations are derived from statistical assumptions of homogeneity, incompressibility, and stationarity, supplemented by foliation ansatz and critical balance scaling (yielding, for axisymmetric cases, a unique exact third-order law) (Galtier, 2010).

3. Explicit Formulations and Functional Structure

The explicit structure of third-order conservation laws varies by context:

• In nonlinear wave equations, the third-order conserved functional has the form

$$I_3 = q \|\Delta u_t\|^2 + \frac{1}{q}\|\nabla \Delta u\|^2 - q' \int \Delta u\, \Delta u_t\, dx + \frac{1}{8}(q')^2 [q \|\nabla u_t\|^2 + \frac{1}{q} \|\Delta u\|^2] - \frac{a}{16}\left(\frac{a^2 (s')^4}{4} + q^2 (s'')^2\right)$$

where $q(t)$, $q'(t)$, $s'(t)$, $s''(t)$ are as defined in the underlying equation; the functional combines space-time energies and cross terms (Boiti et al., 13 Dec 2025).

• In discrete third-order difference equations, the conserved function has the generic form

$\Phi(n, u_n, u_{n+1}, u_{n+2}) = \text{const.}$

with explicit formulas derived from recurrence relations for multipliers $B_i(n)$ (Mamba et al., 2017).

• In axisymmetric MHD turbulence, the third-order moment (Kolmogorov–Yaglom analog) reads

$$\langle |\delta \mathbf{z}^\pm(\mathbf{r})|^2 \, \delta \mathbf{z}^\mp(\mathbf{r}) \rangle = -g(\theta)\, \varepsilon^\pm\, r\, \mathbf{e}_T(\theta)$$

where $g(\theta)$ and $\mathbf{e}_T(\theta)$ encode anisotropy and geometric projection due to a strong mean field (Galtier, 2010).

4. A Priori Bounds and Significance for Existence Theory

Third-order conservation laws are instrumental in providing a priori estimates on higher derivatives or moments, crucial for global existence and regularity results. For the Kirchhoff–Pokhozhaev equation, under "small energy" assumptions ($I_1 \leq 1/(6|ab|)$), it is proven that the $L^2$ norms of derivatives up to third order are uniformly bounded in time, via coercivity properties of $I_3$ (Boiti et al., 13 Dec 2025). This strengthens the global-existence theory for solutions in low-regularity regimes and demonstrates the practical functional analytic utility of higher-order conservation.

5. Applications in Discrete Systems and Numerical Schemes

Conservation laws of third order are crucial for characterizing integrable discrete systems and ensuring accuracy in numerical simulation of nonlinear PDEs:

• In discrete third-order difference equations, Lie symmetry–based methods facilitate reduction of order and construction of first integrals, determining solvability and invariant structure (Mamba et al., 2017).
• The truncation error analysis for third-order MUSCL finite-volume schemes proves that, with specific reconstruction parameter ($\kappa=1/3$ in Van Leer’s scheme), exact cubic accuracy and true third-order convergence are attained for cell–averaged nonlinear conservation laws; improper coupling with inconsistent diffusion operators may degrade this to second order (Nishikawa, 2020).
• In the relative-velocity Lattice Boltzmann algorithm, the third-order equivalent conservation law explicitly captures the leading dispersive and diffusive corrections, with terms systematically derived via Taylor and moment expansions and parameterized by relaxation and moment matrix structures (Graille et al., 2015).

6. Generalizations, Limitations, and Open Directions

Third-order conservation laws exhibit specificity to structural features–for example, the Kirchhoff–Pokhozhaev third-order law is unique to the profile $m(s)=(a\,s+b)^{-2}$. For generic Kirchhoff-type equations, only first-order energy is typically conserved (Boiti et al., 13 Dec 2025). Open problems include the extension to weak solutions, alternative boundary conditions, degenerate parameter cases (e.g., $b=0$ in Kirchhoff–Pokhozhaev), and the construction of higher-order conservation laws (fourth, fifth, etc.) via generalized quadratic forms or symmetry flows.

In numerical analysis, the compatibility of discretization schemes with underlying conservation law order remains an area of ongoing refinement, especially for complex unstructured or higher-dimensional grids.

7. Interdisciplinary Connections and Physical Relevance

Third-order conservation laws have foundational implications for fluid dynamics, plasma physics, nonlinear wave propagation, and statistical turbulence. In MHD turbulence, such laws underpin explicit expressions for energy transfer and anisotropic scaling, generalizing classical Kolmogorov hypotheses (Galtier, 2010). In continuum mechanics and evolutionary biology, self-adjointness and derived conservation vectors relate to invariant measures and pattern formation (Freire, 2010). The concept of third-order invariance provides a unifying theme linking analytic, algebraic, and computational perspectives across mathematical physics.

---

Third-order law context	Conserved quantity structure	Significance
Kirchhoff–Pokhozhaev PDE (Boiti et al., 13 Dec 2025)	$I_3$: mixed space-time energy functional	A priori bounds, global existence
Discrete difference equations (Mamba et al., 2017)	$\Phi$: shift-invariant functional	Integrability, symmetry reductions
MHD turbulence (Galtier, 2010)	Third-order moment of Elsässer fields	Energy flux, anisotropy scaling
Self-adjoint evolution PDE (Freire, 2010)	Conserved vector $(C^0,C^1)$ (via symmetries)	Non-Lagrangian invariants
Numerical schemes (Nishikawa, 2020, Graille et al., 2015)	Discrete third-order truncation/conservation	High-order accuracy, algorithm design

A plausible implication is that as analytic, discrete, and computational models become more sophisticated, third-order conservation laws provide essential foundation for advanced control of solution regularity, structure-preserving algorithms, and deeper statistical insight into multi-scale phenomena.