[论文解读] The Euler Equations and Non-Local Conservative Riccati Equations

Peter Constantin|ArXiv.org|Apr 11, 2000

Navier-Stokes equation solutions参考文献 4被引用 6

一句话总结

本文提出了一类无穷维的精确、光滑解，适用于不可压缩3D欧拉方程，其动能无限，且在有限时间内发生爆破。通过垂直涡度分量的非局部保守型Riccati方程，作者利用特征线方法推导出解的显式公式，证明在具有负均值涡度的一般初始条件下，解的最大值与最小值会同时发散——排除了一侧爆破的可能性。

ABSTRACT

We present an infinite dimensional family of of exact solutions of the incompressible three-dimensional Euler equations. These solutions, proposed by Gibbon and Ohkitani, have infinite kinetic energy and blow up in finite time.

研究动机与目标

• 构造3D不可压缩欧拉方程的显式解，其动能无限且在有限时间内爆破。
• 分析具有非平凡涡度结构的解中有限时间奇点形成的机制。
• 确立爆破不可能为单侧（即仅最大值或仅最小值发散）的结论，解决数值模拟提出的问题。
• 利用特征线方法与非局部Riccati方程，对解的动力学行为提供完整的解析描述。
• 证明速度场的旋度在爆破机制中仅起次要作用，因为即使初始涡度为零，相同的奇点行为依然出现。

提出的方法

• 通过假设构造解：u(x,y,z,t) = (u(x,y,t), zγ(x,y,t))，其中γ为周期函数且u为无散度。
• γ的演化由一个非局部保守型Riccati方程控制：∂γ/∂t + u·∇γ = -γ² + I(t)，其中I(t) = (2/|Q|)∫γ²dx。
• 涡度ω满足∂ω/∂t + u·∇ω = γω，且通过泊松方程与u通过流函数和势函数相关联。
• 采用特征线方法求解系统，其中从标签到位置的映射A(x,t)与时间重参数化τ(t)由初始数据导出。
• 爆破时间T*通过初始数据的二重积分计算得出：T* = (1/|Q|²)∫∫ [1/(γ₀(a)−γ₀(b))] log[(γ₀(a)−m₀)/(γ₀(b)−m₀)] dada。
• 特征线上的解表示为γ(X(a,t),t) = α(τ) [γ₀(a)/(1+τγ₀(a)) − φ̄(τ)]，其中α(τ)与φ̄(τ)通过初始数据的积分定义。

实验结果

研究问题

• RQ1 能否构造出具有无限动能的3D不可压缩欧拉方程的精确、光滑解，使其在有限时间内爆破？
• RQ2 此类解中奇点形成的精确机制是什么，特别是关于涡度与速度分量的关系？
• RQ3 在一般初始条件下，是否可能发生单侧爆破（即仅最大值或仅最小值发散）？
• RQ4 γ方程的非局部Riccati结构如何影响全局动力学与爆破时间？
• RQ5 涡度ω在爆破过程中起多大作用？其是否可被解耦？

主要发现

• 对于光滑、均值为零的初始数据γ₀(x) ∈ C∞(𝕋²)且满足min(γ₀) = m₀ < 0的解，其在有限时间T*内爆破。
• 爆破时间由下式显式给出：T* = (1/|Q|²)∫∫ [1/(γ₀(a)−γ₀(b))] log[(γ₀(a)−m₀)/(γ₀(b)−m₀)] dada。
• γ的最大值与最小值同时发散：在T*时刻，max(γ) → +∞且min(γ) → −∞，从而排除了单侧爆破的可能性。
• 时间上的爆破速率渐近为T* − t ∼ (τ* − τ)(1 + log(1/(τ* − τ)))²，其中τ* = −1/m₀。
• 特征线上的解为γ(X(a,t),t) = α(τ) [γ₀(a)/(1+τγ₀(a)) − φ̄(τ)]，其中当τ → τ*时，α(τ) ∼ (log(τ* − τ))⁻²且φ̄(τ) ∼ −C/(τ* − τ) log(1/(τ* − τ))。
• 流映射的雅可比行列式为J(a,t) = [1/(1+τγ₀(a))] / (1/|Q|∫ da/(1+τγ₀(a)))，确保爆破前体积保持不变且映射可逆。