Blow-Up Equations: Analysis & Criteria

• Blow-up equations are differential equations—including ODEs, PDEs, SPDEs, and integro-differential equations—that exhibit finite-time singularities through unbounded growth or loss of regularity.
• They are analyzed using comparison principles, scaling laws, critical exponents, and invariant quantities, which help distinguish regimes of global existence versus blow-up.
• Extensions to stochastic, fractional, and geometric settings demonstrate how memory effects, noise, and compactification techniques contribute to singularity formation, supported by both theoretical and numerical studies.

Blow-up equations are ordinary differential equations, partial differential equations, stochastic partial differential equations, and integro-differential equations whose maximal solutions can cease to exist in finite time through unbounded growth, loss of regularity, or breakdown of the underlying geometric flow. The literature represented here treats blow-up as a structural phenomenon rather than a single pathology: in some settings every nontrivial positive solution becomes unbounded in finite time, in others only derivatives or moments diverge, and in geometric formulations the decisive event is degeneration of a Jacobian or loss of invertibility of the flow map. Across these formulations, blow-up is typically organized by comparison principles, scaling laws, critical integrals, and invariant quantities (Laister et al., 2019).

1. Definitions and singularity notions

For autonomous ODEs $y'=f(y)$, a solution is a blow-up solution when its maximal existence time $t_{\max}$ is finite; this is distinct from grow-up, where divergence occurs only as $t\to\infty$ (Takayasu et al., 2016). In semilinear fractional heat equations on $\mathbb R^n$, the phrase “blow-up property” is used for the stronger statement that every nontrivial nonnegative solution becomes unbounded in finite time (Laister et al., 2019). For nonlinear Volterra integro-differential equations, the same finite-time/non-finite-time dichotomy is retained, but the growth law is mediated by the full past history through a convolution kernel (Appleby et al., 2017).

Blow-up does not always mean pointwise divergence of the solution itself. For the periodic $b$-family,

$u_t-u_{txx}+(b+1)uu_x=b\,u_xu_{xx}+uu_{xxx},$

finite-time blow-up occurs if and only if the slope becomes unbounded below,

$$\lim_{t\to T^*}\Big((2b-1)\inf_{x\in\mathbb T}u_x(t,x)\Big)=-\infty,$$

so the singularity mechanism is wave breaking: the solution remains bounded while $u_x\to-\infty$ (Estrella, 2014). In the Lagrangian formulation of Euler–Arnold equations, blow-up is detected by the vanishing of the radial Jacobian factor $\gamma_r$; once $\gamma_r$ reaches zero, the flow ceases to be locally invertible, and this forces $t_{\max}$0-blow-up of the velocity field through $t_{\max}$1 (Bauer et al., 2023).

Several papers use probabilistic notions of blow-up. For semilinear parabolic equations driven by Lévy noise on bounded domains, the relevant event is finite-time blow-up in mean square,

$t_{\max}$2

for some finite $t_{\max}$3 (Mohan et al., 2024). For space-time fractional stochastic reaction-diffusion equations, blow-up means non-existence of a global random field solution because moments such as $t_{\max}$4 or $t_{\max}$5 become infinite in finite time (Asogwa et al., 2018). Other settings distinguish finer singularity classes: the critical wave equation on curved backgrounds admits finite-energy type II blow-up (Nahas et al., 2013), singular Liouville equations distinguish simple from non-simple blow-up through the spherical Harnack inequality (Wu, 2023), and the fourth-order parabolic equation $t_{\max}$6 exhibits complete finite-time and complete infinite-time blow-up, meaning that every spatial point diverges, possibly with different signs in different regions (Escudero, 2021).

2. Thresholds, critical exponents, and dichotomies

A central theme is that blow-up is often governed by a sharp threshold. For the fractional semilinear heat equation

$t_{\max}$7

with $t_{\max}$8 locally Lipschitz, non-decreasing, convex, positive for $t_{\max}$9, and satisfying the near-zero scaling condition (S), there is an exact equivalence between three statements: every nontrivial nonnegative PDE solution blows up in finite time, every positive solution of

$t\to\infty$0

blows up in finite time, and the critical near-zero divergence condition

$t\to\infty$1

holds (Laister et al., 2019). The complementary integrability condition yields positive global solutions. For pure powers $t\to\infty$2, this reproduces the fractional Fujita exponent

$t\to\infty$3

with universal blow-up for $t\to\infty$4 and small-data global existence for $t\to\infty$5 (Laister et al., 2019).

Nonlocal diffusion with convolution kernels exhibits an analogous Fujita-type structure. For

$t\to\infty$6

where $t\to\infty$7 is radially symmetric, $t\to\infty$8, and $t\to\infty$9 near $\mathbb R^n$0, the backward-kernel functional

$\mathbb R^n$1

satisfies a scalar inequality implying blow-up whenever $\mathbb R^n$2, with $\mathbb R^n$3 (Biler, 2018). For $\mathbb R^n$4, the same critical exponent appears,

$\mathbb R^n$5

and the supercritical regime $\mathbb R^n$6 admits a partial dichotomy: small critical Morrey norm gives global-in-time smooth solutions with decay $\mathbb R^n$7, while large critical Morrey norm gives finite-time blow-up (Biler, 2018).

Memory terms can alter the mechanism without destroying sharp thresholds. For the Volterra equation

$\mathbb R^n$8

finite-time blow-up occurs if and only if

$\mathbb R^n$9

for some $b$0; if the integral diverges for all $b$1, solutions are global (Appleby et al., 2017). The same paper derives sharp asymptotic laws

$b$2

showing that the local kernel mass at the origin, rather than the detailed global shape of the kernel, controls both explosive and nonexplosive superlinear growth (Appleby et al., 2017).

3. Local breakdown and continuation criteria

Many blow-up theories are formulated as continuation criteria: a strong solution persists as long as a critical quantity stays integrable or bounded. For the periodic $b$3-family with $b$4, the main result is local in space. If there exists $b$5 such that

$b$6

then the corresponding strong solution blows up in finite time, with the explicit upper bound

$b$7

Here $b$8 is defined through the variational quantity $b$9, and the proof proceeds along characteristics through the combinations $u_t-u_{txx}+(b+1)uu_x=b\,u_xu_{xx}+uu_{xxx},$0 and Riccati-type differential inequalities (Estrella, 2014).

For the Thermal Quasi-Geostrophic system on $u_t-u_{txx}+(b+1)uu_x=b\,u_xu_{xx}+uu_{xxx},$1, the continuation mechanism is of Beale–Kato–Majda type. If

$u_t-u_{txx}+(b+1)uu_x=b\,u_xu_{xx}+uu_{xxx},$2

then the strong solution extends beyond $u_t-u_{txx}+(b+1)uu_x=b\,u_xu_{xx}+uu_{xxx},$3. Consequently, if the maximal time is finite, then necessarily

$u_t-u_{txx}+(b+1)uu_x=b\,u_xu_{xx}+uu_{xxx},$4

and in particular $u_t-u_{txx}+(b+1)uu_x=b\,u_xu_{xx}+uu_{xxx},$5 (Crisan et al., 2022). The key analytic input is a logarithmic estimate for the velocity arising from the modified Helmholtz relation $u_t-u_{txx}+(b+1)uu_x=b\,u_xu_{xx}+uu_{xxx},$6, which plays the role that the Biot–Savart law plays in the classical Euler BKM criterion (Crisan et al., 2022).

The three-dimensional compressible Navier–Stokes equations admit analogous Serrin-type criteria tailored to compressible flow. In the barotropic case, finite-time breakdown of a strong solution implies

$u_t-u_{txx}+(b+1)uu_x=b\,u_xu_{xx}+uu_{xxx},$7

for some $u_t-u_{txx}+(b+1)uu_x=b\,u_xu_{xx}+uu_{xxx},$8 depending only on $u_t-u_{txx}+(b+1)uu_x=b\,u_xu_{xx}+uu_{xxx},$9; in the heat-conducting case the criterion becomes

$$\lim_{t\to T^*}\Big((2b-1)\inf_{x\in\mathbb T}u_x(t,x)\Big)=-\infty,$$0

for some $$\lim_{t\to T^*}\Big((2b-1)\inf_{x\in\mathbb T}u_x(t,x)\Big)=-\infty,$$1 depending on the physical parameters (Choe et al., 2017). A notable feature is that these criteria are proved under only the physical viscosity conditions $$\lim_{t\to T^*}\Big((2b-1)\inf_{x\in\mathbb T}u_x(t,x)\Big)=-\infty,$$2 and $$\lim_{t\to T^*}\Big((2b-1)\inf_{x\in\mathbb T}u_x(t,x)\Big)=-\infty,$$3, without additional restrictions on the Lamé coefficients (Choe et al., 2017).

Geometric Euler–Arnold theory replaces pointwise comparison by comparison in an infinite-dimensional function space. In radial variables, the transported momentum identity

$$\lim_{t\to T^*}\Big((2b-1)\inf_{x\in\mathbb T}u_x(t,x)\Big)=-\infty,$$4

leads to a Liouville-type comparison principle for a weighted Jacobian quantity $$\lim_{t\to T^*}\Big((2b-1)\inf_{x\in\mathbb T}u_x(t,x)\Big)=-\infty,$$5. Under positivity and log-supermodularity hypotheses on the Green kernel and a sign condition $$\lim_{t\to T^*}\Big((2b-1)\inf_{x\in\mathbb T}u_x(t,x)\Big)=-\infty,$$6, the comparison theorem forces $$\lim_{t\to T^*}\Big((2b-1)\inf_{x\in\mathbb T}u_x(t,x)\Big)=-\infty,$$7 and hence $$\lim_{t\to T^*}\Big((2b-1)\inf_{x\in\mathbb T}u_x(t,x)\Big)=-\infty,$$8 to vanish in finite time, yielding $$\lim_{t\to T^*}\Big((2b-1)\inf_{x\in\mathbb T}u_x(t,x)\Big)=-\infty,$$9-breakdown (Bauer et al., 2023). In the radial Hunter–Saxton case this reduces to the exact criterion

$u_x\to-\infty$0

which is equivalent to blow-up in any dimension (Bauer et al., 2023).

4. Blow-up profiles, rates, and local structure

Besides criteria, several papers construct explicit or asymptotic blow-up profiles. For the critical focusing wave equation on a curved three-dimensional background,

$u_x\to-\infty$1

there exists, for every $u_x\to-\infty$2, a finite-energy type II blow-up solution of the form

$u_x\to-\infty$3

The parameter $u_x\to-\infty$4 yields a continuum of blow-up rates. The construction extends the slow blow-up theory of Krieger, Schlag, and Tataru to curved backgrounds, but the curvature-induced perturbation $u_x\to-\infty$5 breaks exact scale invariance and limits the method to $u_x\to-\infty$6 (Nahas et al., 2013).

The local structure of blow-up can be highly sensitive to oscillation properties. For singular Liouville equations

$u_x\to-\infty$7

Wei–Zhang proved vanishing theorems for non-simple blow-up sequences. Wu showed that these vanishing conclusions are genuinely tied to the non-simple regime. In the non-quantized case $u_x\to-\infty$8, there exist blow-up sequences satisfying spherical Harnack near the origin while both

$u_x\to-\infty$9

remain uniformly nonzero. In the quantized case $\gamma_r$0, there exist simple blow-up sequences with

$\gamma_r$1

showing that the non-simple assumption in Wei–Zhang’s Laplacian vanishing theorem is essential (Wu, 2023).

The fourth-order parabolic equation

$\gamma_r$2

provides explicit examples in which blow-up is complete and sign-sensitive. On the square, the disc, and $\gamma_r$3, the constructed solutions exhibit both finite-time and infinite-time blow-up; in several families every spatial point diverges, but different regions may approach $\gamma_r$4 and $\gamma_r$5 in finite time (Escudero, 2021). The same paper refines an earlier criterion by showing that, in the initial-Dirichlet weak-solution framework, $\gamma_r$6-blow-up implies blow-up in the stronger norm $\gamma_r$7, via a Gagliardo–Nirenberg interpolation argument (Escudero, 2021). A plausible implication is that a priori criteria based only on $\gamma_r$8-control can miss the actual geometric sharpness of Hessian-driven singularity formation.

5. Stochastic, fractional, and coupled hyperbolic extensions

In stochastic parabolic problems, blow-up theory must absorb Itô corrections and jump compensators. For semilinear stochastic parabolic equations on a bounded smooth domain with Dirichlet boundary conditions,

$\gamma_r$9

finite-time mean-square blow-up is proved for both additive Lévy noise and linear multiplicative Lévy noise (Mohan et al., 2024). The proof adapts the classical concavity method by introducing

$\gamma_r$0

and deriving

$\gamma_r$1

for suitable $\gamma_r$2. In the multiplicative case the noise strength enters through

$\gamma_r$3

and the restriction $\gamma_r$4 ensures that the destabilizing noise can be absorbed by Laplacian dissipation via Poincaré’s inequality (Mohan et al., 2024).

Space-time fractional stochastic reaction-diffusion equations introduce both temporal memory and anomalous diffusion,

$\gamma_r$5

Here blow-up means non-existence of a global random field solution because suitable moments become infinite in finite time (Asogwa et al., 2018). Under superlinear lower bounds such as $\gamma_r$6 or $\gamma_r$7, Walsh-isometry and Jensen-type arguments lead to nonlinear renewal inequalities of the form

$\gamma_r$8

which force finite-time blow-up under the stated kernel and initial-data hypotheses (Asogwa et al., 2018). The resulting theorems cover white noise, spatially colored noise, Riesz kernel covariances, whole-space and bounded-domain geometries, and deterministic drift-plus-noise models (Asogwa et al., 2018).

Weakly coupled Euler–Poisson–Darboux–Tricomi systems show that time-dependent propagation speeds and derivative nonlinearities preserve a sharp blow-up geometry when the two components share the same speed $\gamma_r$9. For

$t_{\max}$00

the blow-up region is

$t_{\max}$01

and the maximal existence time satisfies polynomial or exponential upper bounds depending on whether the problem is subcritical, critical, or double-critical (Hassen et al., 2023). The same-speed assumption is used repeatedly in the support property, in the construction of the adjoint test functions, and in closing the coupled ODE comparison argument. The mass terms $t_{\max}$02 do not appear in the blow-up region or the lifespan estimate (Hassen et al., 2023).

6. Compactification, dynamical systems at infinity, and numerical validation

A major line of work interprets blow-up as asymptotic dynamics on a compactified phase space. For asymptotically quasi-homogeneous ODEs, quasi-Poincaré compactification sends infinity to a horizon $t_{\max}$03, and a time desingularization produces a vector field that extends continuously to that horizon (Matsue, 2016). Hyperbolic equilibria and periodic orbits on $t_{\max}$04 then generate finite-time blow-up along their stable manifolds. If $t_{\max}$05 is a hyperbolic equilibrium at infinity, the corresponding original solution satisfies

$t_{\max}$06

and, for nonzero components,

$t_{\max}$07

The qualitative dynamics at infinity are proved to be topologically equivalent across quasi-Poincaré, directional, and intermediate compactifications (Matsue, 2016).

Validated numerics make this dynamical-systems picture constructive. For polynomial ODEs, admissible compactifications together with a normalized vector field and a quadratic Lyapunov function around a critical point at infinity allow one to prove that a computed orbit truly blows up and to enclose the blow-up time rigorously (Takayasu et al., 2016). The method uses interval or affine arithmetic to validate a Lyapunov domain and then estimates

$t_{\max}$08

by explicit upper and lower bounds. Poincaré compactification is used when the normalized field is regular at the boundary; parabolic compactification is used when the Poincaré Jacobian becomes singular near infinity (Takayasu et al., 2016).

For unresolved hydrodynamic PDEs, the numerical problem is different: the issue is not formal validation of a known trajectory but distinguishing genuine singularity formation from discretization artifacts. In a complexified Euler–EPDiff model on $t_{\max}$09, a geometrically consistent Zeitlin-type discretization is used to identify a numerical signature of blow-up: the supremum norm of vorticity should show a resolution-dependent fingerprint as the truncation level $t_{\max}$10 increases (Jansson et al., 2022). The proposed diagnostic is not merely large $t_{\max}$11, but persistent and stronger growth of $t_{\max}$12 under refinement, together with computational stability checks based on the structure-preserving discretization. This suggests a practical criterion for future numerical blow-up studies in equations whose analytical status remains open (Jansson et al., 2022).

Taken together, these results show that blow-up equations are unified less by a single model than by a repertoire of analytic and geometric mechanisms: Fujita-type threshold integrals, Riccati inequalities along characteristics, Liouville comparison in Banach spaces, stochastic concavity, renewal inequalities, explicit profile construction, and compactification at infinity. The common structure is that finite-time breakdown becomes detectable once the correct critical quantity has been identified—whether it is a solution norm, a derivative, a weighted moment, a Jacobian, or a trajectory on the horizon at infinity.