Proudman–Johnson Equation

• The Proudman–Johnson equation is a family of nonlinear 1-D evolution equations characterized by transport–stretching dynamics, flow-map degeneration, and finite-time singularity formation.
• It offers multiple parameterizations that connect it to classical models like Burgers, Hunter–Saxton, and even non-right-invariant geometric flows.
• Modern analysis employs Lagrangian reduction and flow-map formulations to study phenomena such as blow-up behavior, global existence, and stability under various boundary conditions.

The Proudman–Johnson equation denotes a family of nonlinear one-dimensional evolution equations that appear in several closely related normal forms, and whose modern analysis centers on transport–stretching dynamics, flow-map degeneration, and finite-time singularity formation. In contemporary PDE literature, the term usually refers to the generalized inviscid family

$u_{txx}+u\,u_{xxx}-a\,u_xu_{xx}=0,$

or equivalent parameterizations such as

$u_{txx}+(1+2\lambda)u_xu_{xx}+u\,u_{xxx}=0,$

together with periodic, Dirichlet, or non-periodic boundary conditions; the name also persists in the older reversed-stagnation-point literature for a distinct similarity-reduced equation of boundary-layer type (Kogelbauer, 2019, Bauer et al., 2021, Chio, 2013).

1. Nomenclature and normal forms

A standard modern formulation is the generalized Proudman–Johnson equation

$u_{txx}+u\,u_{xxx}-a\,u_xu_{xx}=0,$

studied on either $[0,1]$ with Dirichlet boundary conditions or on the periodic domain with mean-free periodic boundary conditions. In this convention, the classical Proudman–Johnson equation is the case $a=1$; the same paper records $a=-3$ as Burgers and $a=-2$ as Hunter–Saxton (Kogelbauer, 2019).

A second common convention writes the generalized inviscid family on $\mathbb R$ or $S^1$ as

$u_{txx}+(1+2\lambda)u_xu_{xx}+u\,u_{xxx}=0.$

In the non-periodic theory, this is paired with the identification $u_{txx}+(1+2\lambda)u_xu_{xx}+u\,u_{xxx}=0,$0, which links the equation to the $u_{txx}+(1+2\lambda)u_xu_{xx}+u\,u_{xxx}=0,$1-Hunter–Saxton equation (Bauer et al., 2021). A third convention, used in an information-geometric formulation, is

$u_{txx}+(1+2\lambda)u_xu_{xx}+u\,u_{xxx}=0,$2

where $u_{txx}+(1+2\lambda)u_xu_{xx}+u\,u_{xxx}=0,$3 gives the original Proudman–Johnson equation, $u_{txx}+(1+2\lambda)u_xu_{xx}+u\,u_{xxx}=0,$4 gives the Hunter–Saxton equation, and $u_{txx}+(1+2\lambda)u_xu_{xx}+u\,u_{xxx}=0,$5 gives the $u_{txx}+(1+2\lambda)u_xu_{xx}+u\,u_{xxx}=0,$6-Burgers equation (Bauer et al., 1 Aug 2025).

For periodic problems, an integrated form is often more convenient: $u_{txx}+(1+2\lambda)u_xu_{xx}+u\,u_{xxx}=0,$7 with periodic boundary conditions on $u_{txx}+(1+2\lambda)u_xu_{xx}+u\,u_{xxx}=0,$8 and $u_{txx}+(1+2\lambda)u_xu_{xx}+u\,u_{xxx}=0,$9 (Sarria et al., 2013). In this formulation, the nonlocal term $u_{txx}+u\,u_{xxx}-a\,u_xu_{xx}=0,$0 enforces periodic compatibility.

The older aerodynamic usage is different. In the reversed stagnation-point setting, the similarity-reduced equation associated with Proudman and Johnson’s 1962 analysis is

$u_{txx}+u\,u_{xxx}-a\,u_xu_{xx}=0,$1

with

$u_{txx}+u\,u_{xxx}-a\,u_xu_{xx}=0,$2

That equation arises from a stream-function ansatz for two-dimensional reversed stagnation-point flow near a flat wall (Chio, 2013). This suggests that “Proudman–Johnson equation” functions partly as a historical label spanning more than one reduction framework.

2. Characteristic and flow-map formulations

The modern theory is largely driven by Lagrangian reduction. For periodic generalized inviscid Proudman–Johnson,

$u_{txx}+u\,u_{xxx}-a\,u_xu_{xx}=0,$3

and

$u_{txx}+u\,u_{xxx}-a\,u_xu_{xx}=0,$4

Introducing

$u_{txx}+u\,u_{xxx}-a\,u_xu_{xx}=0,$5

together with

$u_{txx}+u\,u_{xxx}-a\,u_xu_{xx}=0,$6

one obtains

$u_{txx}+u\,u_{xxx}-a\,u_xu_{xx}=0,$7

and the representation formula

$u_{txx}+u\,u_{xxx}-a\,u_xu_{xx}=0,$8

The same framework gives

$u_{txx}+u\,u_{xxx}-a\,u_xu_{xx}=0,$9

so the sign pattern of $[0,1]$0 is transported by the flow (Sarria et al., 2013).

A complementary flow-map formulation on $[0,1]$1 or the periodic domain starts from

$[0,1]$2

with $[0,1]$3. Along trajectories,

$[0,1]$4

For $[0,1]$5, setting $[0,1]$6 yields the inhomogeneous Liouville equation

$[0,1]$7

This can be solved explicitly up to a scalar auxiliary function $[0,1]$8, leading to

$[0,1]$9

In this formulation, continuation is equivalent to preserving the diffeomorphism property $a=1$0; breakdown occurs when

$a=1$1

for some $a=1$2 (Kogelbauer, 2019).

3. Blow-up, global existence, and the role of initial curvature

For smooth periodic data with quadratic behavior near extrema, the generalized inviscid Proudman–Johnson equation admits a sharp parameter-dependent classification. In the $a=1$3-formulation, solutions are global for $a=1$4. For $a=1$5, finite-time blow-up is two-sided and everywhere in the sense that

$a=1$6

and for every non-extremal label $a=1$7, $a=1$8. For $a=1$9, blow-up is one-sided and discrete: only the minimum diverges. For $a=-3$0, blow-up is again two-sided and everywhere, but now non-minimizing points diverge to $a=-3$1 while minimizing points diverge to $a=-3$2 (Sarria et al., 2013).

Norm growth is more delicate than $a=-3$3-blow-up. For $a=-3$4, $a=-3$5 for $a=-3$6, and also for $a=-3$7. The energy

$a=-3$8

remains finite for $a=-3$9 and diverges for $a=-2$0, so one-sided blow-up on $a=-2$1 can remain compatible with finite $a=-2$2-energy up to blow-up time (Sarria et al., 2013).

The curvature of the initial slope near its extrema refines these thresholds. If, near each maximizing point,

$a=-2$3

then the decisive positive-$a=-2$4 threshold is

$a=-2$5

Solutions are global for $a=-2$6; if $a=-2$7, then $a=-2$8, whereas $a=-2$9 gives convergence to a nontrivial steady state. Finite-time two-sided, everywhere blow-up occurs for $\mathbb R$0. For negative $\mathbb R$1, $\mathbb R$2 produces one-sided discrete blow-up, while for $\mathbb R$3 and

$\mathbb R$4

the blow-up becomes two-sided and everywhere (Sarria et al., 2013). For smooth data with $\mathbb R$5 vanishing to order $\mathbb R$6 at the relevant extremum, $\mathbb R$7, so the positive threshold becomes $\mathbb R$8 (Sarria et al., 2013).

An alternative interval/periodic flow-map theory in the $\mathbb R$9-parameterization gives global existence for $S^1$0, global existence for $S^1$1 provided

$S^1$2

and global existence for $S^1$3 under the smallness condition

$S^1$4

The same theory gives a blow-up criterion for $S^1$5 in terms of square-integrability of

$S^1$6

on

$S^1$7

(Kogelbauer, 2019).

4. Damping and nonstandard boundary conditions

A damped version of the generalized inviscid Proudman–Johnson equation has been studied on $S^1$8 under the homogeneous three-point boundary condition

$S^1$9

The equation is

$u_{txx}+(1+2\lambda)u_xu_{xx}+u\,u_{xxx}=0.$0

with

$u_{txx}+(1+2\lambda)u_xu_{xx}+u\,u_{xxx}=0.$1

Writing $u_{txx}+(1+2\lambda)u_xu_{xx}+u\,u_{xxx}=0.$2, one obtains the identity

$u_{txx}+(1+2\lambda)u_xu_{xx}+u\,u_{xxx}=0.$3

and therefore the Riccati-type inequality

$u_{txx}+(1+2\lambda)u_xu_{xx}+u\,u_{xxx}=0.$4

This scalar reduction is the core of the blow-up mechanism (Obi-Okoye et al., 2021).

For bounded smooth damping, if

$u_{txx}+(1+2\lambda)u_xu_{xx}+u\,u_{xxx}=0.$5

then the solution blows up in finite time at $u_{txx}+(1+2\lambda)u_xu_{xx}+u\,u_{xxx}=0.$6,

$u_{txx}+(1+2\lambda)u_xu_{xx}+u\,u_{xxx}=0.$7

with explicit upper bound

$u_{txx}+(1+2\lambda)u_xu_{xx}+u\,u_{xxx}=0.$8

The result shows that bounded damping shifts the threshold quantitatively but does not prevent blow-up for sufficiently negative boundary data (Obi-Okoye et al., 2021).

The same paper proves that even unbounded smooth damping may fail to arrest singularity formation. For the explicit choice

$u_{txx}+(1+2\lambda)u_xu_{xx}+u\,u_{xxx}=0.$9

if

$u_{txx}+(1+2\lambda)u_xu_{xx}+u\,u_{xxx}=0,$00

where $u_{txx}+(1+2\lambda)u_xu_{xx}+u\,u_{xxx}=0,$01 is the exponential integral, then there exists a finite $u_{txx}+(1+2\lambda)u_xu_{xx}+u\,u_{xxx}=0,$02 such that

$u_{txx}+(1+2\lambda)u_xu_{xx}+u\,u_{xxx}=0,$03

The proved singularity is boundary-value blow-up of $u_{txx}+(1+2\lambda)u_xu_{xx}+u\,u_{xxx}=0,$04 itself, rather than a derivative-only blow-up scenario (Obi-Okoye et al., 2021).

5. Geometric and variational interpretations

On the real line, the non-periodic generalized inviscid Proudman–Johnson equation is formally equivalent to the $u_{txx}+(1+2\lambda)u_xu_{xx}+u\,u_{xxx}=0,$05-Hunter–Saxton equation under

$u_{txx}+(1+2\lambda)u_xu_{xx}+u\,u_{xxx}=0,$06

In this regime, it is the Eulerian geodesic equation of a right-invariant homogeneous $u_{txx}+(1+2\lambda)u_xu_{xx}+u\,u_{xxx}=0,$07-Finsler metric on

$u_{txx}+(1+2\lambda)u_xu_{xx}+u\,u_{xxx}=0,$08

The decisive linearizing map is

$u_{txx}+(1+2\lambda)u_xu_{xx}+u\,u_{xxx}=0,$09

which is an isometric embedding into an open convex subset of $u_{txx}+(1+2\lambda)u_xu_{xx}+u\,u_{xxx}=0,$10. Geodesics therefore become straight lines, and the explicit solution formula is

$u_{txx}+(1+2\lambda)u_xu_{xx}+u\,u_{xxx}=0,$11

The maximal positive existence time is

$u_{txx}+(1+2\lambda)u_xu_{xx}+u\,u_{xxx}=0,$12

whenever some $u_{txx}+(1+2\lambda)u_xu_{xx}+u\,u_{xxx}=0,$13; otherwise the solution exists for all $u_{txx}+(1+2\lambda)u_xu_{xx}+u\,u_{xxx}=0,$14. The same work emphasizes that the equivalence with $u_{txx}+(1+2\lambda)u_xu_{xx}+u\,u_{xxx}=0,$15-Hunter–Saxton fails on the circle (Bauer et al., 2021).

A different geometric interpretation uses information geometry. On the space of positive densities, the Amari–Čencov $u_{txx}+(1+2\lambda)u_xu_{xx}+u\,u_{xxx}=0,$16-connections $u_{txx}+(1+2\lambda)u_xu_{xx}+u\,u_{xxx}=0,$17 are realized as Levi-Civita connections of explicit metrics $u_{txx}+(1+2\lambda)u_xu_{xx}+u\,u_{xxx}=0,$18. Pulling $u_{txx}+(1+2\lambda)u_xu_{xx}+u\,u_{xxx}=0,$19 back by

$u_{txx}+(1+2\lambda)u_xu_{xx}+u\,u_{xxx}=0,$20

gives the metric

$u_{txx}+(1+2\lambda)u_xu_{xx}+u\,u_{xxx}=0,$21

and the corresponding Euler–Arnold equation is

$u_{txx}+(1+2\lambda)u_xu_{xx}+u\,u_{xxx}=0,$22

Here the metric is right-invariant iff $u_{txx}+(1+2\lambda)u_xu_{xx}+u\,u_{xxx}=0,$23; for $u_{txx}+(1+2\lambda)u_xu_{xx}+u\,u_{xxx}=0,$24, the generalized Proudman–Johnson equation is the Euler–Arnold equation of a non-right-invariant metric. The same theory states that these generalized Proudman–Johnson equations are globally well-posed on $u_{txx}+(1+2\lambda)u_xu_{xx}+u\,u_{xxx}=0,$25 iff $u_{txx}+(1+2\lambda)u_xu_{xx}+u\,u_{xxx}=0,$26 (Bauer et al., 1 Aug 2025).

6. Self-similar regimes and recent developments

Recent periodic work near the classical threshold $u_{txx}+(1+2\lambda)u_xu_{xx}+u\,u_{xxx}=0,$27 studies

$u_{txx}+(1+2\lambda)u_xu_{xx}+u\,u_{xxx}=0,$28

under odd symmetry. In this convention, $u_{txx}+(1+2\lambda)u_xu_{xx}+u\,u_{xxx}=0,$29 is the classical Proudman–Johnson equation. There exists $u_{txx}+(1+2\lambda)u_xu_{xx}+u\,u_{xxx}=0,$30 such that for some $u_{txx}+(1+2\lambda)u_xu_{xx}+u\,u_{xxx}=0,$31 initial data the inviscid equation $u_{txx}+(1+2\lambda)u_xu_{xx}+u\,u_{xxx}=0,$32 blows up in finite time when $u_{txx}+(1+2\lambda)u_xu_{xx}+u\,u_{xxx}=0,$33. The same paper constructs self-similar solutions

$u_{txx}+(1+2\lambda)u_xu_{xx}+u\,u_{xxx}=0,$34

for $u_{txx}+(1+2\lambda)u_xu_{xx}+u\,u_{xxx}=0,$35, with $u_{txx}+(1+2\lambda)u_xu_{xx}+u\,u_{xxx}=0,$36 in the supercritical range $u_{txx}+(1+2\lambda)u_xu_{xx}+u\,u_{xxx}=0,$37, $u_{txx}+(1+2\lambda)u_xu_{xx}+u\,u_{xxx}=0,$38 at $u_{txx}+(1+2\lambda)u_xu_{xx}+u\,u_{xxx}=0,$39, and $u_{txx}+(1+2\lambda)u_xu_{xx}+u\,u_{xxx}=0,$40 for $u_{txx}+(1+2\lambda)u_xu_{xx}+u\,u_{xxx}=0,$41, yielding $u_{txx}+(1+2\lambda)u_xu_{xx}+u\,u_{xxx}=0,$42 decay in the subcritical regime. The same analysis also proves finite-time self-similar blow-up for some $u_{txx}+(1+2\lambda)u_xu_{xx}+u\,u_{xxx}=0,$43 data at the critical case $u_{txx}+(1+2\lambda)u_xu_{xx}+u\,u_{xxx}=0,$44, and finite-time blow-up in the viscous case $u_{txx}+(1+2\lambda)u_xu_{xx}+u\,u_{xxx}=0,$45 for $u_{txx}+(1+2\lambda)u_xu_{xx}+u\,u_{xxx}=0,$46 sufficiently close to $u_{txx}+(1+2\lambda)u_xu_{xx}+u\,u_{xxx}=0,$47 (Guo et al., 19 Nov 2025).

A separate recent development uses the periodic mean-free equation

$u_{txx}+(1+2\lambda)u_xu_{xx}+u\,u_{xxx}=0,$48

as the renormalized form of a blow-up problem for the incompressible porous medium equation. Under the corresponding change of variables, the explicit blow-up profile becomes the steady state

$u_{txx}+(1+2\lambda)u_xu_{xx}+u\,u_{xxx}=0,$49

The stationary $u_{txx}+(1+2\lambda)u_xu_{xx}+u\,u_{xxx}=0,$50 solutions are completely classified: they are exactly

$u_{txx}+(1+2\lambda)u_xu_{xx}+u\,u_{xxx}=0,$51

For every $u_{txx}+(1+2\lambda)u_xu_{xx}+u\,u_{xxx}=0,$52, small $u_{txx}+(1+2\lambda)u_xu_{xx}+u\,u_{xxx}=0,$53 perturbations of $u_{txx}+(1+2\lambda)u_xu_{xx}+u\,u_{xxx}=0,$54 are globally stable and satisfy

$u_{txx}+(1+2\lambda)u_xu_{xx}+u\,u_{xxx}=0,$55

where

$u_{txx}+(1+2\lambda)u_xu_{xx}+u\,u_{xxx}=0,$56

is selected by the initial curvature at the point of maximum. By contrast, for every $u_{txx}+(1+2\lambda)u_xu_{xx}+u\,u_{xxx}=0,$57 there exist perturbations arbitrarily small in $u_{txx}+(1+2\lambda)u_xu_{xx}+u\,u_{xxx}=0,$58 for which the solution cannot converge in $u_{txx}+(1+2\lambda)u_xu_{xx}+u\,u_{xxx}=0,$59 to any $u_{txx}+(1+2\lambda)u_xu_{xx}+u\,u_{xxx}=0,$60. This identifies a sharp regularity threshold between $u_{txx}+(1+2\lambda)u_xu_{xx}+u\,u_{xxx}=0,$61-type stability and $u_{txx}+(1+2\lambda)u_xu_{xx}+u\,u_{xxx}=0,$62-instability, and it simultaneously yields corresponding stability statements for associated special classes of solutions of 2D Euler and inviscid primitive or hydrostatic Euler (Collot et al., 23 Jul 2025).