Suitable Weak Solutions in Nonlinear PDEs

• Suitable weak solutions are defined as weak solutions to nonlinear PDEs that satisfy an additional local energy inequality reflecting physical dissipation.
• They play a critical role in partial regularity theory by quantifying singular sets and underpinning ε-regularity criteria in systems like Navier–Stokes.
• This framework has been adapted to various complex systems, including hyperdissipative fluids, non-Newtonian flows, and multi-physics PDEs, ensuring rigorous control over local energy behavior.

A suitable weak solution is a distinguished class of weak solution to nonlinear PDEs for which a local energy inequality reflecting physical dissipation (or entropy) holds in addition to the usual distributional and global energy properties. This notion was first introduced for the incompressible Navier–Stokes equations by Scheffer and Caffarelli–Kohn–Nirenberg (CKN), and is central to the modern partial regularity theory of singularities in critical nonlinear evolution and stationary systems. The definition and analysis of suitable weak solutions have been extended to a wide range of systems (including hyperdissipative fluids, non-Newtonian flows, boundary and high-dimensional settings, as well as stochastic and coupled multi-physics PDEs), often providing the essential framework for establishing $\varepsilon$-regularity criteria and quantifying the size and structure of singular sets.

1. Core Definition and Analytical Structure

For the (possibly fractional) Navier–Stokes equations, a suitable weak solution imposes not only the usual functional and weak distributional requirements but also an explicit local energy-type inequality. In the hyperdissipative Navier–Stokes case, with dissipation $(-\Delta)^{\alpha}$, $\alpha\in(1,5/4]$, and for divergence-free initial data $u_0\in L^2(\mathbb{R}^3)$, a pair $(u,p)$ is suitable if:

• $$u\in L^{\infty}(0,T;L^2(\mathbb{R}^3))\cap L^2(0,T;\dot{H}^{\alpha}(\mathbb{R}^3))$$, $p\in L^{3/2}_{loc}(\mathbb{R}^3\times(0,T))$;
• the pair solves the system in the sense of distributions;
• the global energy inequality holds for almost every $0

$$\frac12\|u(t)\|_{L^2}^2 + \int_s^t\|(-\Delta)^{\alpha/2}u(\tau)\|_{L^2}^2\,d\tau \le \frac12\|u(s)\|_{L^2}^2;$$

• and for all nonnegative $$\varphi\in C_c^{\infty}(\mathbb{R}^3_x\times[0,\infty)_y\times(0,T))$$ with suitable boundary behavior in the Caffarelli–Silvestre–Yang extension variable $(-\Delta)^{\alpha}$0, the local energy inequality (with the Dirichlet form for the extension) is satisfied:

$(-\Delta)^{\alpha}$1

$(-\Delta)^{\alpha}$2

In classical (α=1) or stationary contexts, the inequality simplifies accordingly; the unifying feature is the inclusion of a localized dissipation term, integration against arbitrary nonnegative space–time test functions, and the control of flux terms—see (Colombo et al., 2017, Huang et al., 4 Nov 2025, Donatelli et al., 2010).

2. Distinction from Generic Weak Solutions

Whereas Leray–Hopf weak solutions (“energy class”) are constructed via global energy bounds and compactness, they generally lack enough local structure to control fine properties of singularities. The additional local energy inequality required for suitability ensures that potential “dissipation defects” do not concentrate on sets of positive parabolic measure. In particular, the Caffarelli–Kohn–Nirenberg partial regularity theory applies only to suitable weak solutions. For instance, even Leray solutions constructed via mollification or artificial compressibility require delicate pressure and time-derivative bounds in order to pass to the suitable (localized) energy regime (Donatelli et al., 2010).

3. Suitability in Generalized and Complex Fluid Systems

The notion of suitability extends to a variety of nonlinear evolution systems beyond Newtonian incompressible fluid mechanics. In each such context, the technical form of the local energy inequality is adapted to the nonlinearities of the system. Examples include:

• Magnetohydrodynamics (MHD): The velocity $(-\Delta)^{\alpha}$3, magnetic field $(-\Delta)^{\alpha}$4, and pressure $(-\Delta)^{\alpha}$5 satisfy a coupled system, and the suitable weak solution class incorporates additional local energy inequalities involving $(-\Delta)^{\alpha}$6, dissipation $(-\Delta)^{\alpha}$7, and cross terms reflecting Lorentz-force and induction effects (Kang et al., 2012).
• Non-Newtonian (Power-law) Flows: Suitability consists in the distributional formulation and a generalized local energy inequality for the power-law stress $(-\Delta)^{\alpha}$8; the dissipative term is $(-\Delta)^{\alpha}$9, and the local energy inequality reflects the higher integrability and the appropriate monotonicity structure (Kim, 7 Feb 2026).
• Multi-physics Systems: For coupled PDEs—e.g., Navier–Stokes–Planck–Nernst–Poisson (Gong et al., 2019), Beris–Edwards Q-tensor models (Du et al., 2019), and chemotaxis–Navier–Stokes (Chen et al., 2023)—the suitable weak solution requires additional tradeoffs between energy, entropy, and higher-order dissipations, all in localized (test function) form, ensuring that cross-effects (chemical, electrostatic, or orientational) respect dissipation and regularity at the localization level.

4. Partial Regularity and $\alpha\in(1,5/4]$0-Regularity Criteria

The essential motivation for introducing the class of suitable weak solutions is to facilitate partial regularity theory. The Caffarelli–Kohn–Nirenberg theorem asserts: for any suitable weak solution (Navier–Stokes or its generalizations), the set of singular points (where the solution is not locally Hölder continuous) has parabolic Hausdorff dimension at most $\alpha\in(1,5/4]$1 (or lower depending on the dissipation). For the hyperdissipative case, this bound improves to $\alpha\in(1,5/4]$2 (Colombo et al., 2017).

Analytically, this result is underpinned by an $\alpha\in(1,5/4]$3-regularity principle: if the “excess” (a scaling-invariant quantity measuring deviation from regularity, e.g., $\alpha\in(1,5/4]$4-norm of velocity, $\alpha\in(1,5/4]$5-norm of pressure) is sufficiently small on some parabolic cylinder, then the solution is smooth on a smaller subcylinder. The strategy involves:

• Scaling reduction to minimal blow-up;
• A compactness and Liouville-type classification of possible blow-up limits (self-similar, ancient, or homogeneous solutions, depending on the PDE);
• Correction and absorption of nonlocal and nonlinear error terms (including fractional extension energies, pressure projections, or subgrid stress);
• Iteration of excess decay to propagate regularity (Colombo et al., 2017, Huang et al., 4 Nov 2025, Seregin, 2019, Jiu et al., 2018).

This structure extends to various inhomogeneous and stochastic models, provided the suitable local energy balance is maintained (Chen et al., 2024).

5. Extensions to Numerical Schemes and Approximation

The suitability property has been established for several classes of approximate or numerical solutions to Navier–Stokes and its generalizations. Fully discretized schemes in space–time, finite element or variational multiscale (VMS) with pressure stabilization and subgrid modeling, can be shown to converge to suitable weak solutions provided that:

• The schemes are consistent and stable in energy;
• Discrete commutator properties and inf–sup stability are available or compensated by orthogonal subgrid scales;
• Testing discrete solutions by localized energy multipliers (mollified or projected) passes to the limit without creating nonphysical artifacts.

This ensures the reliability of such schemes for capturing physically meaningful dissipation and the structure of possible singular sets (Berselli et al., 2017, Badia et al., 2016).

6. Impact, Uniqueness, and Further Developments

The framework of suitable weak solutions has deep implications for uniqueness, regularity, and the structure of turbulence. Weak–strong uniqueness results state that if a suitable weak solution coincides (initially) with a classical (strong or mild) solution, then the two must agree as long as the strong solution exists. Such results have been strengthened to include solution classes in weighted Lebesgue and Besov spaces (Lemarié-Rieusset, 2021).

Moreover, the quantification of “gaps” between energy inequalities and equalities in the suitable class connects to nonuniqueness scenarios (convex integration, concentration–dissipation anomalies) (Maremonti et al., 2022, Crispo et al., 2019). In stochastic PDEs, martingale suitable weak solutions are developed to provide a probabilistic analogue of local energy decay, with stochastic correction and martingale supermartingale structure, preserving dissipative regularity mechanisms (Chen et al., 2024).

The concept continues to evolve, with explicit regularity criteria and localized versions (e.g., weighted and exterior-Serrin-type conditions), as well as extensions to high dimensions, flows with complex boundaries, and coupled field interactions (Huang et al., 4 Nov 2025, Neustupa, 2013).

7. Summary Table: Canonical Criteria for Suitability

System/Class	Key Integrability for Suitability	Local Energy Inequality
3D Navier–Stokes	$\alpha\in(1,5/4]$6	For all $\alpha\in(1,5/4]$7, $\alpha\in(1,5/4]$8 and $\alpha\in(1,5/4]$9-level with pressure term
Hyperdissipative Navier–Stokes	$u_0\in L^2(\mathbb{R}^3)$0	With fractional extension variable, bulk Dirichlet energy (Colombo et al., 2017)
Power-law Non-Newtonian	$u_0\in L^2(\mathbb{R}^3)$1	Generalized, with $u_0\in L^2(\mathbb{R}^3)$2 dissipation (Kim, 7 Feb 2026)
MHD, Chemotaxis models, Nernst–Planck	Various: velocity, fields, densities	Summed or coupled local energy/entropy balancing each physical effect
Stochastic Navier–Stokes	$u_0\in L^2(\mathbb{R}^3)$3	Martingale local energy inequality in expectation (Chen et al., 2024)

All variants maintain the central principle: suitable weak solutions not only satisfy the distributional and global energy inequalities, but also a strong enough local energy property to ensure that the singular set is small and regularity can be propagated from smallness in critical quantities, often via $u_0\in L^2(\mathbb{R}^3)$4-regularity. This framework is foundational for the analysis of singularities, regularity, and dissipative mechanisms in nonlinear PDEs.