We present a novel approach to the Liouville problem for the stationary Navier-Stokes equations. As an application of our method, we prove conditional Liouville theorems with assumptions on the antiderivative of the velocity that represent substantial improvements on what was heretofore known.
论文ID : 2506.14533标题 : Conditional Liouville Theorems for the Navier-Stokes Equations作者 : Matei P. Coiculescu (Princeton University), Jincheng Yang (Johns Hopkins University)分类 : math.AP (Analysis of PDEs)发表时间 : November 25, 2025 (arXiv v2)论文链接 : https://arxiv.org/abs/2506.14533v2 本文提出了一种研究稳态Navier-Stokes方程Liouville问题的新方法。作为该方法的应用，作者证明了在速度场反导数（antiderivative）假设下的条件Liouville定理，这些结果相比已有研究有实质性改进。
Liouville问题是稳态Navier-Stokes方程研究中的核心问题：在三维空间R³中，是否存在非平凡的D-解（具有有限Dirichlet能量且在无穷远处趋于零的解）？
具体而言，考虑稳态Navier-Stokes方程：
− Δ u + u ⋅ ∇ u + ∇ p = 0 -\Delta u + u \cdot \nabla u + \nabla p = 0 − Δ u + u ⋅ ∇ u + ∇ p = 0 div u = 0 \text{div } u = 0 div u = 0 lim ∣ x ∣ → ∞ u ( x ) = 0 \lim_{|x|\to\infty} u(x) = 0 lim ∣ x ∣ → ∞ u ( x ) = 0 ∫ R 3 ∣ ∇ u ∣ 2 d x < ∞ \int_{R^3} |\nabla u|^2 dx < \infty ∫ R 3 ∣∇ u ∣ 2 d x < ∞ 物理意义 ：从物理角度，粘性流体的能量耗散应使其趋向静止状态（u ≡ 0 u \equiv 0 u ≡ 0 数学基础 ：这是理解流体力学偏微分方程长时间行为的关键问题尺度不变性 ：该问题具有尺度不变性，需要尺度不变的假设才能获得深刻结果Galdi的结果 ：若u ∈ L 9 / 2 ( R 3 ) u \in L^{9/2}(R^3) u ∈ L 9/2 ( R 3 ) u ≡ 0 u \equiv 0 u ≡ 0 u ∈ L 6 ( R 3 ) u \in L^6(R^3) u ∈ L 6 ( R 3 ) Seregin的BMO⁻¹方法 ：若u ∈ L 6 ∩ BMO − 1 u \in L^6 \cap \text{BMO}^{-1} u ∈ L 6 ∩ BMO − 1 u ≡ 0 u \equiv 0 u ≡ 0 Seregin-Chae-Wolf的量化版本 ：对流函数ψ \psi ψ u = curl ψ u = \text{curl }\psi u = curl ψ s > 3 s > 3 s > 3 现有方法都在固定尺度 的球上进行局部估计，本文提出变尺度的胶囊（capsule）方法 ：
在每个点周围构造适应性的胶囊区域 胶囊的尺寸根据局部Dirichlet能量和平均速度自适应调整 核心思想：如果速度场"拉伸"不太严重，则可证明Liouville定理 提出了胶囊（capsule）方法 ：一种新的变尺度局部估计框架，胶囊的半径和长度根据局部信息自适应选择改进了流函数平均振荡条件 （定理1.1）：扩展到s ≥ 1 s \geq 1 s ≥ 1 s > 3 s > 3 s > 3 对3 < s < 7 3 < s < 7 3 < s < 7 对3 < s < 9 / 2 3 < s < 9/2 3 < s < 9/2 阈值α < 1 / 9 \alpha < 1/9 α < 1/9 s s s 建立了线积分控制条件 （定理1.2）：首次使用速度场线积分的增长控制：∫ x 0 x u ⋅ d ℓ ≤ C ∣ x − x 0 ∣ β \int_{x_0}^x u \cdot d\ell \leq C|x-x_0|^\beta ∫ x 0 x u ⋅ d ℓ ≤ C ∣ x − x 0 ∣ β 当β < 29 / 193 \beta < 29/193 β < 29/193 u ≡ 0 u \equiv 0 u ≡ 0 提供了统一的条件Liouville定理框架 ：通过控制"拉伸比"L / R L/R L / R 证明在特定条件下，稳态Navier-Stokes方程的D-解必然是平凡解（u ≡ 0 u \equiv 0 u ≡ 0
定义 ：胶囊C R , L , e ( x ) C_{R,L,e}(x) C R , L , e ( x ) R R R L ≥ R L \geq R L ≥ R e ∈ S 2 e \in S^2 e ∈ S 2 x x x C R , L , e ( x ) : = ⋃ t ∈ [ − L + R , L − R ] ( t e + B R ( x ) ) C_{R,L,e}(x) := \bigcup_{t\in[-L+R,L-R]} (te + B_R(x)) C R , L , e ( x ) := ⋃ t ∈ [ − L + R , L − R ] ( t e + B R ( x ))
几何上，这是一个"胶囊形"区域，由沿方向e e e B R B_R B R
在胶囊C C C
\fint B R u d x = U e 1 \fint_{B_R} u dx = Ue_1 \fint B R u d x = U e 1 \fint 2 C ∣ ∇ u ∣ 2 d x = Ξ 2 \fint_{2C} |\nabla u|^2 dx = \Xi^2 \fint 2 C ∣∇ u ∣ 2 d x = Ξ 2 关键引理2.1 ：在ansatz条件L Ξ ≤ ε 0 R ( L / R ) δ L\Xi \leq \varepsilon_0 \frac{R}{(L/R)^\delta} L Ξ ≤ ε 0 ( L / R ) δ R ∥ u − U e 1 ∥ L ∞ ( C ) ≲ ε 0 ( U R L + 1 R ) \|u - Ue_1\|_{L^\infty(C)} \lesssim \varepsilon_0\left(\frac{UR}{L} + \frac{1}{R}\right) ∥ u − U e 1 ∥ L ∞ ( C ) ≲ ε 0 ( L U R + R 1 )
证明策略 ：
控制平均速度 （2.1节）：利用不可压缩性，证明沿流线方向的平均速度振荡有界：
∣ u ˉ ( t ) − U e 1 ∣ ≲ L Ξ |\bar{u}(t) - Ue_1| \lesssim L\Xi ∣ u ˉ ( t ) − U e 1 ∣ ≲ L Ξ 速度分解 （2.2节）：将u u u u = U e 1 + v + h u = Ue_1 + v + h u = U e 1 + v + h v = curl ( − Δ ) − 1 ( ϕ ω ) v = \text{curl}(-\Delta)^{-1}(\phi\omega) v = curl ( − Δ ) − 1 ( ϕ ω ) h h h 3 2 C \frac{3}{2}C 2 3 C 涡度方程估计 （2.3节）：对涡度方程
b ⋅ ∇ ω + curl ( ω × ζ ) = Δ ω b \cdot \nabla\omega + \text{curl}(\omega \times \zeta) = \Delta\omega b ⋅ ∇ ω + curl ( ω × ζ ) = Δ ω ω ∈ L 6 \omega \in L^6 ω ∈ L 6 u ∈ L ∞ u \in L^\infty u ∈ L ∞ 流线极大函数 （定义3.1）：对流映射Φ s \Phi_s Φ s ∂ s Φ s = u ( Φ s ) \partial_s\Phi_s = u(\Phi_s) ∂ s Φ s = u ( Φ s ) ( M Φ f ) ( x ) : = sup s > 0 1 2 s ∫ − s s ∣ f ∘ Φ τ ( x ) ∣ d τ (M_\Phi f)(x) := \sup_{s>0} \frac{1}{2s}\int_{-s}^s |f \circ \Phi_\tau(x)| d\tau ( M Φ f ) ( x ) := sup s > 0 2 s 1 ∫ − s s ∣ f ∘ Φ τ ( x ) ∣ d τ
引理3.2 ：M Φ M_\Phi M Φ ( p , p ) (p,p) ( p , p ) 1 < p < ∞ 1 < p < \infty 1 < p < ∞
胶囊构造算法 （命题3.4）：对每个点x ∈ R 3 x \in R^3 x ∈ R 3 C x = C R ( x ) , L ( x ) , e ( x ) ( x ) C_x = C_{R(x),L(x),e(x)}(x) C x = C R ( x ) , L ( x ) , e ( x ) ( x )
计算平均速度：\fint B R ( x ) u d x = U ( x ) e ( x ) \fint_{B_R(x)} u dx = U(x)e(x) \fint B R ( x ) u d x = U ( x ) e ( x ) 定义长度：L ( x ) = max { U ( x ) R , 1 } 1 1 + σ R L(x) = \max\{U(x)R, 1\}^{\frac{1}{1+\sigma}}R L ( x ) = max { U ( x ) R , 1 } 1 + σ 1 R 选择R ( x ) R(x) R ( x ) Ξ ~ ( x ) L 1 − δ R 1 + δ = ε 0 \tilde{\Xi}(x)L^{1-\delta}R^{1+\delta} = \varepsilon_0 Ξ ~ ( x ) L 1 − δ R 1 + δ = ε 0 Ξ ~ 2 = \fint C x M Φ [ M ( ∣ ∇ u ∣ 2 ) ] d x \tilde{\Xi}^2 = \fint_{C_x} M_\Phi[M(|\nabla u|^2)] dx Ξ ~ 2 = \fint C x M Φ [ M ( ∣∇ u ∣ 2 )] d x 点的分类 ：
圆点 （round point）：L ( x ) = R ( x ) L(x) = R(x) L ( x ) = R ( x ) 长点 （long point）：L ( x ) > R ( x ) L(x) > R(x) L ( x ) > R ( x ) 引理3.6 （可比性）：若x , z x, z x , z C x ∩ C z ≠ ∅ C_x \cap C_z \neq \emptyset C x ∩ C z = ∅ R ( z ) ≤ 2 R ( x ) R(z) \leq 2R(x) R ( z ) ≤ 2 R ( x ) C z ⊂ K C x C_z \subset KC_x C z ⊂ K C x
引理3.7 （Vitali型覆盖）：对长点集合A ⊂ L A \subset L A ⊂ L { C x i } \{C_{x_i}\} { C x i } ∣ ⋃ x ∈ A C x ∣ ≤ ∑ i ∣ K C x i ∣ \left|\bigcup_{x\in A} C_x\right| \leq \sum_i |KC_{x_i}| ⋃ x ∈ A C x ≤ ∑ i ∣ K C x i ∣
推论3.8 ：Ξ ~ ∈ L 2 + ε , ∞ ( R 3 ) \tilde{\Xi} \in L^{2+\varepsilon,\infty}(R^3) Ξ ~ ∈ L 2 + ε , ∞ ( R 3 ) L 2 + ε L^{2+\varepsilon} L 2 + ε
自适应尺度选择 ：不同于传统固定尺度方法，胶囊尺寸根据局部能量和速度自适应，能更精细地捕捉解的结构流线极大函数 ：引入M Φ M_\Phi M Φ 拉伸比控制 ：通过控制L / R L/R L / R 统一的阈值 ：定理1.1中α < 1 / 9 \alpha < 1/9 α < 1/9 s ≥ 1 s \geq 1 s ≥ 1 s s s 假设 ：u = curl ψ u = \text{curl }\psi u = curl ψ ( \fint B R ( x 0 ) ∣ ψ − ( ψ ) B R ( x 0 ) ∣ s d x ) 1 / s ≤ C R α , ∀ R > 1 , ∀ x 0 ∈ R 3 \left(\fint_{B_R(x_0)} |\psi - (\psi)_{B_R(x_0)}|^s dx\right)^{1/s} \leq CR^\alpha, \quad \forall R > 1, \forall x_0 \in R^3 ( \fint B R ( x 0 ) ∣ ψ − ( ψ ) B R ( x 0 ) ∣ s d x ) 1/ s ≤ C R α , ∀ R > 1 , ∀ x 0 ∈ R 3 s ≥ 1 s \geq 1 s ≥ 1 C > 0 C > 0 C > 0 0 < α < 1 0 < \alpha < 1 0 < α < 1
结论 ：
u ∈ L p , ∞ u \in L^{p,\infty} u ∈ L p , ∞ p > 4 1 − α p > \frac{4}{1-\alpha} p > 1 − α 4 若α < 1 9 \alpha < \frac{1}{9} α < 9 1 u ≡ 0 u \equiv 0 u ≡ 0 证明思路 （Section 4.2）：
利用引理4.2，将流函数振荡与U R UR U R ∫ B R ( ψ − ψ ˉ ) ⋅ ( e × y ) d y ≈ U R 5 \int_{B_R} (\psi - \bar{\psi}) \cdot (e \times y) dy \approx UR^5 ∫ B R ( ψ − ψ ˉ ) ⋅ ( e × y ) d y ≈ U R 5 结合假设得到：( L / R ) 1 + σ ≲ R α (L/R)^{1+\sigma} \lesssim R^\alpha ( L / R ) 1 + σ ≲ R α 通过选择最优参数σ → ∞ \sigma \to \infty σ → ∞ ∣ u ∣ p / Ξ ~ 2 |u|^p/\tilde{\Xi}^2 ∣ u ∣ p / Ξ ~ 2 u ∈ L p , ∞ u \in L^{p,\infty} u ∈ L p , ∞ 当α < 1 / 9 \alpha < 1/9 α < 1/9 p < 9 / 2 p < 9/2 p < 9/2 u ≡ 0 u \equiv 0 u ≡ 0 假设 ：
∫ x 0 x u ⋅ d ℓ ≤ C ∣ x − x 0 ∣ β , ∀ x 0 , x ∈ R 3 \int_{x_0}^x u \cdot d\ell \leq C|x-x_0|^\beta, \quad \forall x_0, x \in R^3 ∫ x 0 x u ⋅ d ℓ ≤ C ∣ x − x 0 ∣ β , ∀ x 0 , x ∈ R 3 0 < β < 1 0 < \beta < 1 0 < β < 1
结论 ：
u ∈ L p , ∞ u \in L^{p,\infty} u ∈ L p , ∞ p > 4 − 34 β / 29 1 − β p > \frac{4-34\beta/29}{1-\beta} p > 1 − β 4 − 34 β /29 若β < 29 193 \beta < \frac{29}{193} β < 193 29 u ≡ 0 u \equiv 0 u ≡ 0 证明思路 （Section 4.1）：
引理4.1建立：对长点x x x ∫ x − L e x + L e u ⋅ d ℓ ≈ L U \int_{x-Le}^{x+Le} u \cdot d\ell \approx LU ∫ x − L e x + L e u ⋅ d ℓ ≈ LU 结合假设：L U ≈ ( L / R ) 2 + σ ≲ L β LU \approx (L/R)^{2+\sigma} \lesssim L^\beta LU ≈ ( L / R ) 2 + σ ≲ L β 选择最优σ = 5 / 12 \sigma = 5/12 σ = 5/12 δ = 5 / 12 \delta = 5/12 δ = 5/12 p p p 当β < 29 / 193 \beta < 29/193 β < 29/193 p < 9 / 2 p < 9/2 p < 9/2 考虑方程：b ⋅ ∇ θ − Δ θ = f + div g b \cdot \nabla\theta - \Delta\theta = f + \text{div } g b ⋅ ∇ θ − Δ θ = f + div g
基本解 （引理A.1）：
Γ ( x ) = 1 4 π r e − λ ( r − x 1 ) , λ = U 2 \Gamma(x) = \frac{1}{4\pi r}e^{-\lambda(r-x_1)}, \quad \lambda = \frac{U}{2} Γ ( x ) = 4 π r 1 e − λ ( r − x 1 ) , λ = 2 U
关键估计 （引理A.2）：∇ Γ ∈ L x 1 3 / 2 , ∞ L x 2 , x 3 3 / 2 \nabla\Gamma \in L^{3/2,\infty}_{x_1}L^{3/2}_{x_2,x_3} ∇Γ ∈ L x 1 3/2 , ∞ L x 2 , x 3 3/2 λ \lambda λ
局部L p L^p L p （引理A.3）：若θ \theta θ C C C ∥ θ ∥ L r ( C / 2 ) ≲ R ∥ f ∥ L q ( C ) + ∥ g ∥ L q ( C ) + ( U R L + 1 R ) ∥ θ ∥ L q ( C ) \|\theta\|_{L^r(C/2)} \lesssim R\|f\|_{L^q(C)} + \|g\|_{L^q(C)} + \left(\frac{UR}{L} + \frac{1}{R}\right)\|\theta\|_{L^q(C)} ∥ θ ∥ L r ( C /2 ) ≲ R ∥ f ∥ L q ( C ) + ∥ g ∥ L q ( C ) + ( L U R + R 1 ) ∥ θ ∥ L q ( C ) 1 / r = 1 / q − 1 / 3 1/r = 1/q - 1/3 1/ r = 1/ q − 1/3
核心假设（方程22）：
L Ξ ≤ ε 0 R ( L / R ) δ ≲ ε 0 ( U R L + 1 R ) L\Xi \leq \varepsilon_0 \frac{R}{(L/R)^\delta} \lesssim \varepsilon_0\left(\frac{UR}{L} + \frac{1}{R}\right) L Ξ ≤ ε 0 ( L / R ) δ R ≲ ε 0 ( L U R + R 1 )
这确保了局部估计中的非线性项和漂移项可控，是整个方法的关键。
结果 假设条件 适用范围 阈值 Seregin 7 ( \fint B R ∥ ψ − ψ ˉ ∥ s ) 1 / s ≤ C R α \left(\fint_{B_R}\|\psi-\bar{\psi}\|^s\right)^{1/s} \leq CR^\alpha ( \fint B R ∥ ψ − ψ ˉ ∥ s ) 1/ s ≤ C R α s > 3 s > 3 s > 3 α = s − 3 6 ( s − 1 ) \alpha = \frac{s-3}{6(s-1)} α = 6 ( s − 1 ) s − 3 Chae-Wolf 3 ( \fint B R ∥ T − T ˉ ∥ s ) 1 / s ≤ C R α \left(\fint_{B_R}\|T-\bar{T}\|^s\right)^{1/s} \leq CR^\alpha ( \fint B R ∥ T − T ˉ ∥ s ) 1/ s ≤ C R α s > 3 s > 3 s > 3 α = min { 1 3 − 1 s , 1 6 } \alpha = \min\{\frac{1}{3}-\frac{1}{s}, \frac{1}{6}\} α = min { 3 1 − s 1 , 6 1 } 本文定理1.1 ( \fint B R ∥ ψ − ψ ˉ ∥ s ) 1 / s ≤ C R α \left(\fint_{B_R}\|\psi-\bar{\psi}\|^s\right)^{1/s} \leq CR^\alpha ( \fint B R ∥ ψ − ψ ˉ ∥ s ) 1/ s ≤ C R α s ≥ 1 s \geq 1 s ≥ 1 α < 1 9 \alpha < \frac{1}{9} α < 9 1
改进之处 ：
扩展到1 ≤ s ≤ 3 1 \leq s \leq 3 1 ≤ s ≤ 3 ：首次处理低指数情况3 < s < 7 3 < s < 7 3 < s < 7 3 < s < 9 / 2 3 < s < 9/2 3 < s < 9/2 统一阈值 ：α < 1 / 9 \alpha < 1/9 α < 1/9 s s s s s s 衰减假设 ：Galdi 4 ：u ∈ L 9 / 2 ( R 3 ) ⇒ u ≡ 0 u \in L^{9/2}(R^3) \Rightarrow u \equiv 0 u ∈ L 9/2 ( R 3 ) ⇒ u ≡ 0 Chae-Wolf 2 ：对数改进版本 BMO⁻¹方法 ：Seregin-Silvestre-Šverák-Zlatoš 8 ：带散度自由漂移的Laplace方程 Seregin 6 ：u ∈ L 6 ∩ BMO − 1 ⇒ u ≡ 0 u \in L^6 \cap \text{BMO}^{-1} \Rightarrow u \equiv 0 u ∈ L 6 ∩ BMO − 1 ⇒ u ≡ 0 量化BMO⁻¹条件 ：Seregin 7 ：流函数振荡控制（s > 3 s > 3 s > 3 Chae-Wolf 3 ：张量T T T Bang-Yang 1 ：对数改进 本文属于量化BMO⁻¹条件方向，但通过胶囊方法 实现了：
更广的指数范围（s ≥ 1 s \geq 1 s ≥ 1 统一的阈值 新的线积分控制条件 方法论贡献 ：胶囊方法提供了一个灵活的框架，可根据局部信息自适应选择估计尺度定理1.1 ：在流函数平均振荡假设下，对所有s ≥ 1 s \geq 1 s ≥ 1 α < 1 / 9 \alpha < 1/9 α < 1/9 u ≡ 0 u \equiv 0 u ≡ 0 定理1.2 ：首次建立了基于速度场线积分增长的Liouville定理核心洞察 ：Liouville性质的本质是控制速度场的"拉伸"——当L / R L/R L / R u ≡ 0 u \equiv 0 u ≡ 0 阈值gap ：定理1.1：α < 1 / 9 \alpha < 1/9 α < 1/9 α < 1 \alpha < 1 α < 1 定理1.2：β < 29 / 193 ≈ 0.15 \beta < 29/193 \approx 0.15 β < 29/193 ≈ 0.15 β < 1 \beta < 1 β < 1 D-解假设 ：仍需要有限Dirichlet能量假设，未完全解决一般Liouville问题技术复杂性 ：胶囊构造和覆盖引理的证明较为技术性最优性未知 ：阈值1 / 9 1/9 1/9 29 / 193 29/193 29/193 改进阈值 ：寻找更优的参数选择（σ , δ \sigma, \delta σ , δ 推广到其他方程 ：胶囊方法可能适用于其他带漂移的椭圆/抛物方程去除能量假设 ：探索是否可以放松有限Dirichlet能量条件数值验证 ：构造数值算例验证阈值的锐性方法创新性强 ：胶囊方法是genuinely new的思路，突破了固定尺度估计的局限 流线极大函数M Φ M_\Phi M Φ 技术严谨 ：证明完整细致，从局部估计到全局结论逻辑清晰 带漂移Poisson方程的处理（附录A）具有独立价值 结果实质改进 ：扩展到s ≥ 1 s \geq 1 s ≥ 1 统一阈值优于之前依赖s s s 线积分条件是全新视角 写作清晰 ：引言充分阐述了动机和背景 技术路线层次分明 符号记号使用规范 阈值保守 ：α < 1 / 9 \alpha < 1/9 α < 1/9 β < 29 / 193 \beta < 29/193 β < 29/193 α , β < 1 \alpha, \beta < 1 α , β < 1 参数选择（σ = 5 / 12 \sigma = 5/12 σ = 5/12 几何直观不足 ：胶囊构造的几何意义可以更直观地解释 为何L / R L/R L / R 数值验证缺失 ：推广性讨论不足 ：方法是否适用于其他维度或其他方程讨论较少 与时间演化问题的联系未涉及 理论贡献 ：为Liouville问题提供了新的技术工具 胶囊方法可能启发其他问题的研究 预期将被后续工作引用和改进 实用价值 ：主要是理论结果，直接应用价值有限 但对理解流体稳定性有概念意义 可复现性 ：理论研究 ：为研究稳态Navier-Stokes方程提供新工具相关问题 ：可能适用于：其他带漂移的椭圆方程 抛物方程的长时间行为 其他流体力学模型 不适用 ：不适用于时间依赖问题 不适用于可压缩流体 不直接处理湍流等复杂现象 4 G. Galdi. An introduction to the mathematical theory of the Navier-Stokes equations: steady-state problems . Springer, 2011.
6 G. Seregin. Liouville type theorem for stationary Navier-Stokes equations. Nonlinearity , 29(8):2191, 2016.
7 G. Seregin. Remarks on Liouville type theorems for steady-state Navier-Stokes equations. St. Petersburg Mathematical Journal , 30(2):321-328, 2019.
8 G. Seregin, L. Silvestre, V. Šverák, A. Zlatoš. On divergence-free drifts. Journal of Differential Equations , 252(1):505-540, 2012.
本文通过创新的胶囊方法 ，在稳态Navier-Stokes方程的Liouville问题上取得了实质性进展。核心思想是用自适应尺度的胶囊区域进行局部估计，通过控制"拉伸比"L / R L/R L / R s ≥ 1 s \geq 1 s ≥ 1 α < 1 / 9 \alpha < 1/9 α < 1/9
