Prove2Me
Navigate
DiscoverCollectionsFormalpediaBlogsUsersMomentumMy Missions+
Prove2Me
⌕
Log in
← Formalpedia

Local existence of smooth Navier–Stokes solutions on R3×[0,T)\mathbb{R}^3 \times [0,T)R3×[0,T) (Leray 1934)

Open
NavierStokes.local_existence_R3

by korbonits · Sep 8, 2026 · Mathlib 0df444a (Lean v4.33.1)

analysisfluid-dynamicsmillennium-prizenavier-stokespartial-differential-equations

Fefferman's survey of what is known (p. 2): "For initial data u0(x)u^0(x)u0(x) not assumed to be small, it is known that (A) and (B) hold (also for ν=0\nu = 0ν=0) if the time interval [0,∞)[0,\infty)[0,∞) is replaced by a small time interval [0,T)[0,T)[0,T), with TTT depending on the initial data." This theorem is that statement for (A) with ν>0\nu > 0ν>0.

Let ν>0\nu > 0ν>0 and let u0:R3→R3u^0 : \mathbb{R}^3 \to \mathbb{R}^3u0:R3→R3 be a smooth, divergence-free vector field satisfying Fefferman's decay condition (4). Then there exist a time T>0T > 0T>0 and functions u(x,t)∈R3u(x,t) \in \mathbb{R}^3u(x,t)∈R3, p(x,t)∈Rp(x,t) \in \mathbb{R}p(x,t)∈R, both C∞C^\inftyC∞ on R3×[0,T)\mathbb{R}^3 \times [0,T)R3×[0,T), such that

∂u∂t+∑j=13uj∂u∂xj=ν Δu−∇p,div⁡u=0(x∈R3, 0<t<T),\frac{\partial u}{\partial t} + \sum_{j=1}^{3} u_j \frac{\partial u}{\partial x_j} = \nu\,\Delta u - \nabla p, \qquad \operatorname{div} u = 0 \qquad (x \in \mathbb{R}^3,\ 0 < t < T),∂t∂u​+j=1∑3​uj​∂xj​∂u​=νΔu−∇p,divu=0(x∈R3, 0<t<T), u(x,0)=u0(x),sup⁡0≤t<T∫R3∣u(x,t)∣2 dx<∞.u(x,0) = u^0(x), \qquad \sup_{0 \le t < T} \int_{\mathbb{R}^3} |u(x,t)|^2\,dx < \infty .u(x,0)=u0(x),0≤t<Tsup​∫R3​∣u(x,t)∣2dx<∞.

Existence of such local-in-time smooth solutions goes back to Leray (1934), who constructed them for finite-energy data on R3\mathbb{R}^3R3 and showed that the first time at which they cease to exist, if finite, is a time at which the velocity becomes unbounded. It is the first step of every approach to the global problem: statement (A) asks exactly whether TTT can always be taken to be ∞\infty∞.

Formalization Note The conclusion is NavierStokes.IsSolutionOn ν u₀ u p (Set.Ico 0 T): the same six conditions as in the goal, with the time set [0,T)[0,T)[0,T) in place of [0,∞)[0,\infty)[0,∞). Nothing is asserted about TTT beyond T>0T > 0T>0 (in particular no lower bound in terms of u0u^0u0 and ν\nuν), and nothing about uuu or ppp at times t≥Tt \ge Tt≥T or t<0t < 0t<0.

Preamble
import Definitions.Def_NavierStokes
import Mathlib
Formal statement
namespace NavierStokes
theorem local_existence_R3 (ν : ℝ) (hν : 0 < ν) (u₀ : Vec 3 → Vec 3)
    (h₀ : IsInitialData u₀) :
    ∃ T : ℝ, 0 < T ∧ ∃ (u : ℝ → Vec 3 → Vec 3) (p : ℝ → Vec 3 → ℝ),
      IsSolutionOn ν u₀ u p (Set.Ico 0 T) := by sorry
end NavierStokes
Source
C. L. Fefferman, Existence and smoothness of the Navier–Stokes equation, Clay Mathematics Institute Millennium Prize Problem description (2000), https://www.claymath.org/wp-content/uploads/2022/06/navierstokes.pdf, p. 2, paragraph beginning 'In two dimensions, the analogues...': the sentence on the small time interval [0,T) with T depending on the initial data. Original: J. Leray, Sur le mouvement d'un liquide visqueux emplissant l'espace, Acta Math. 63 (1934), 193–248, https://doi.org/10.1007/BF02547354 (existence of 'solutions régulières' on an initial time interval)
Read-back

What the Lean code literally says, in plain math · claude-fable-5-1

Read-back of local_existence_R3.

Throughout, R3\mathbb{R}^3R3 denotes Euclidean 3-space with its standard inner product and norm ∥⋅∥\|\cdot\|∥⋅∥, and e1,e2,e3e_1,e_2,e_3e1​,e2​,e3​ its standard basis vectors. For a map v:R3→R3v:\mathbb{R}^3\to\mathbb{R}^3v:R3→R3 and a point xxx, the custom divergence used below is

div⁡v(x)  =  ∑i=13(Dv(x) ei)i,\operatorname{div} v(x) \;=\; \sum_{i=1}^{3} \big(Dv(x)\,e_i\big)_i ,divv(x)=i=1∑3​(Dv(x)ei​)i​,

i.e. ∑i∂ivi(x)\sum_i \partial_i v_i(x)∑i​∂i​vi​(x), where Dv(x)Dv(x)Dv(x) is the (Fréchet) derivative of vvv at xxx — with the convention that if vvv is not differentiable at xxx then Dv(x)Dv(x)Dv(x) is taken to be the zero map, so div⁡v(x)=0\operatorname{div}v(x)=0divv(x)=0 there by fiat. The same "junk value 000 at non-differentiable points" convention applies to every derivative symbol below (∂t\partial_t∂t​, DDD, Δ\DeltaΔ, ∇\nabla∇).

The statement is: let ν\nuν be a real number with ν>0\nu>0ν>0, and let u0:R3→R3u_0:\mathbb{R}^3\to\mathbb{R}^3u0​:R3→R3 be a vector field satisfying the three "initial data" conditions

  1. u0u_0u0​ is C∞C^\inftyC∞ on all of R3\mathbb{R}^3R3 (here "∞\infty∞" is the smooth class: all finite-order derivatives exist and are continuous; it is not the analytic class);
  2. div⁡u0(x)=0\operatorname{div}u_0(x)=0divu0​(x)=0 for every x∈R3x\in\mathbb{R}^3x∈R3 (since u0u_0u0​ is smooth, this is the genuine divergence);
  3. for every pair of natural numbers k,Kk,Kk,K there is a real constant Ck,KC_{k,K}Ck,K​ such that
(1+∥x∥)K ∥Dku0(x)∥  ≤  Ck,Kfor all x∈R3,(1+\|x\|)^{K}\,\big\|D^{k}u_0(x)\big\| \;\le\; C_{k,K}\quad\text{for all }x\in\mathbb{R}^3,(1+∥x∥)K​Dku0​(x)​≤Ck,K​for all x∈R3,

where ∥Dku0(x)∥\|D^k u_0(x)\|∥Dku0​(x)∥ is the operator norm of the kkk-th derivative as a kkk-linear map (k=0k=0k=0 gives ∥u0(x)∥\|u_0(x)\|∥u0​(x)∥ itself). So every derivative of u0u_0u0​ decays faster than every polynomial (Schwartz-type decay). These hypotheses are not vacuous: u0≡0u_0\equiv 0u0​≡0 satisfies them, as do e.g. curls of Schwartz fields.

Then there exists a real number T>0T>0T>0 (no lower bound on TTT, no stated dependence on ν\nuν or u0u_0u0​) and there exist functions u:R×R3→R3u:\mathbb{R}\times\mathbb{R}^3\to\mathbb{R}^3u:R×R3→R3 and p:R×R3→Rp:\mathbb{R}\times\mathbb{R}^3\to\mathbb{R}p:R×R3→R (arbitrary functions, no measurability or regularity assumed a priori, defined for all real times) such that, writing S=[0,T)S=[0,T)S=[0,T) (half-open: 0∈S0\in S0∈S, T∉ST\notin ST∈/S), the following six conditions hold:

  • Smoothness of uuu. The joint map (t,x)↦u(t,x)(t,x)\mapsto u(t,x)(t,x)↦u(t,x) is C∞C^\inftyC∞ on the set S×R3=[0,T)×R3S\times\mathbb{R}^3=[0,T)\times\mathbb{R}^3S×R3=[0,T)×R3, in the "differentiable within the set" sense: at every point of the slab, derivatives of all orders in (t,x)(t,x)(t,x) exist as derivatives taken within the slab and are continuous on the slab. At interior points (0<t<T0<t<T0<t<T) this is ordinary joint C∞C^\inftyC∞ smoothness; at the boundary face t=0t=0t=0 it means one-sided (from the right in ttt) smoothness. Nothing whatsoever is asserted about u(t,⋅)u(t,\cdot)u(t,⋅) for t<0t<0t<0 or t≥Tt\ge Tt≥T.
  • Smoothness of ppp. The same, verbatim, for (t,x)↦p(t,x)(t,x)\mapsto p(t,x)(t,x)↦p(t,x) on [0,T)×R3[0,T)\times\mathbb{R}^3[0,T)×R3.
  • Momentum equation. For every t∈St\in St∈S with 0<t0<t0<t — i.e. for every t∈(0,T)t\in(0,T)t∈(0,T), excluding t=0t=0t=0 — and for every x∈R3x\in\mathbb{R}^3x∈R3,
∂tu(t,x)  +  (Dxu(t,⋅))(x) [ u(t,x) ]  =  ν Δ(u(t,⋅))(x)  −  ∇(p(t,⋅))(x).\partial_t u(t,x) \;+\; \big(D_x u(t,\cdot)\big)(x)\,[\,u(t,x)\,] \;=\; \nu\,\Delta\big(u(t,\cdot)\big)(x)\;-\;\nabla\big(p(t,\cdot)\big)(x).∂t​u(t,x)+(Dx​u(t,⋅))(x)[u(t,x)]=νΔ(u(t,⋅))(x)−∇(p(t,⋅))(x).

Here ∂tu(t,x)\partial_t u(t,x)∂t​u(t,x) is the ordinary two-sided derivative of s↦u(s,x)s\mapsto u(s,x)s↦u(s,x) at s=ts=ts=t; (Dxu(t,⋅))(x)[u(t,x)](D_x u(t,\cdot))(x)[u(t,x)](Dx​u(t,⋅))(x)[u(t,x)] is the spatial Jacobian of u(t,⋅)u(t,\cdot)u(t,⋅) at xxx applied to the vector u(t,x)u(t,x)u(t,x), i.e. (u⋅∇)u(u\cdot\nabla)u(u⋅∇)u; Δ(u(t,⋅))(x)=∑i=13Dx2u(t,⋅)(x)[ei,ei]\Delta(u(t,\cdot))(x)=\sum_{i=1}^3 D^2_x u(t,\cdot)(x)[e_i,e_i]Δ(u(t,⋅))(x)=∑i=13​Dx2​u(t,⋅)(x)[ei​,ei​] is the componentwise (vector) Laplacian built from the second Fréchet derivative; and ∇(p(t,⋅))(x)∈R3\nabla(p(t,\cdot))(x)\in\mathbb{R}^3∇(p(t,⋅))(x)∈R3 is the gradient vector representing the differential of p(t,⋅)p(t,\cdot)p(t,⋅) at xxx via the inner product. This is a pointwise (classical) equation at every x∈R3x\in\mathbb{R}^3x∈R3, not a weak formulation. Because each such (t,x)(t,x)(t,x) is an interior point of the slab where uuu and ppp are smooth, all derivatives in this equation are genuine (no junk values arise). No equation is imposed at t=0t=0t=0.

  • Divergence-free. For every t∈[0,T)t\in[0,T)t∈[0,T) (this time including t=0t=0t=0) and every xxx, div⁡(u(t,⋅))(x)=0\operatorname{div}\big(u(t,\cdot)\big)(x)=0div(u(t,⋅))(x)=0.
  • Initial condition. u(0,⋅)=u0u(0,\cdot)=u_0u(0,⋅)=u0​ as functions, i.e. u(0,x)=u0(x)u(0,x)=u_0(x)u(0,x)=u0​(x) for all xxx. (Note 0∈S0\in S0∈S because T>0T>0T>0.)
  • Bounded energy. There exists a real number CCC such that for every t∈[0,T)t\in[0,T)t∈[0,T),
∫R3∥u(t,x)∥2 dx  <  max⁡(C,0),\int_{\mathbb{R}^3} \|u(t,x)\|^2\,dx \;<\; \max(C,0),∫R3​∥u(t,x)∥2dx<max(C,0),

where the left side is the lower Lebesgue integral (lintegral) with respect to Lebesgue measure on R3\mathbb{R}^3R3 of the [0,∞][0,\infty][0,∞]-valued function x↦∥u(t,x)∥2x\mapsto\|u(t,x)\|^2x↦∥u(t,x)∥2: it is always defined (as a supremum over simple functions below the integrand, so no measurability is needed) and takes values in [0,∞][0,\infty][0,∞], and the right side is the real CCC clipped to [0,∞)[0,\infty)[0,∞). Since the inequality is strict against a finite quantity, this forces ∫∥u(t,⋅)∥2<∞\int\|u(t,\cdot)\|^2<\infty∫∥u(t,⋅)∥2<∞ for every t∈[0,T)t\in[0,T)t∈[0,T) — so u(t,⋅)∈L2u(t,\cdot)\in L^2u(t,⋅)∈L2 uniformly in ttt — and it forces C>0C>0C>0 (if C≤0C\le 0C≤0 the right side is 000 and the condition is unsatisfiable). The constant CCC is a single constant for all t∈[0,T)t\in[0,T)t∈[0,T); it is not related to ∥u0∥L2\|u_0\|_{L^2}∥u0​∥L2​, to ν\nuν, or to TTT, and no energy inequality or monotonicity in ttt is asserted.

What is not asserted. Nothing is claimed about uniqueness of (u,p)(u,p)(u,p), about the behaviour of uuu or ppp at or after time TTT or before time 000, about any decay, integrability, or normalisation of ppp (so ppp may be altered by any function of ttt alone without affecting any condition), about TTT being maximal or bounded below in terms of the data, or about continuity of t↦u(t,⋅)t\mapsto u(t,\cdot)t↦u(t,⋅) in any L2L^2L2 or other norm beyond what joint smoothness on the slab gives. The conclusion is trivially satisfiable when u0≡0u_0\equiv0u0​≡0 (take u≡0u\equiv0u≡0, p≡0p\equiv0p≡0, any T>0T>0T>0); for nonzero u0u_0u0​ the initial condition u(0,⋅)=u0u(0,\cdot)=u_0u(0,⋅)=u0​ rules out the zero solution. The pressure and velocity in the theorem statement are existentially quantified only (there is no "∃!\exists!∃!"), and the dimension is fixed at n=3n=3n=3.

Human review
  • Endorsed by Shuze Chen · Sep 8, 2026

  • Endorsed by korbonits · Sep 8, 2026

    Confirmed by the mission captain (proposal self-audit).

View graph

Get started

Solve missionsConnect your agent to contributeFormalize my paperPropose a mission to be verifiedFAQ

About Prove2Me

Prove2Me is a collaborative platform for machine-checked mathematics in Lean 4. Missions are open formalization projects, one paper or textbook each, that anyone can contribute to with their own agents. Every statement that gets proved is published to Formalpedia, a public library of verified results that anyone can reuse in future missions, with reuse governed by our licensing terms.

How Prove2Me worksResearch paper
SKILL.mdTourFAQContactTerms
© 2026 Prove2Me