[论文解读] On Euler's equation and `EPDiff'

David Mumford, Peter W. Michor|arXiv (Cornell University)|Sep 28, 2012

Nonlinear Waves and Solitons参考文献 6被引用 3

一句话总结

本文证明了不可压缩、无粘性流体流动的欧拉方程解是通过正则化算子 $ L_{\varepsilon,\eta} $ 在微分同胚群上对 EPDiff（Euler-Poincaré 方程）族测地线方程的解的极限，其中参数 $ \varepsilon, \eta \geq 0 $。关键结果是在时间局部范围内，EPDiff 解与欧拉方程之间存在收敛性界，且流场中的动量守恒确保了稳定性与紧支集的传播。

ABSTRACT

We study a family of approximations to Euler's equation depending on two parameters $\varepsilon,η\ge 0$. When $\varepsilon=η=0$ we have Euler's equation and when both are positive we have instances of the class of integro-differential equations called EPDiff in imaging science. These are all geodesic equations on either the full diffeomorphism group $\operatorname{Diff}_{H^\infty}(\mathbb R^n)$ or, if $\varepsilon = 0$, its volume preserving subgroup. They are defined by the right invariant metric induced by the norm on vector fields given by $$ \|v\|_{\varepsilon,η} = \int_{\mathbb R^n} dx $$ where $L_{\varepsilon,η} = (I- frac{η^2}{p} riangle)^p \circ (I- frac1{\varepsilon^2} abla \circ ÷)$. All geodesic equations are locally well-posed, and the $L_{\varepsilon,η}$-equation admits solutions for all time if $η>0$ and $p\ge (n+3)/2$. We tie together solutions of all these equations by estimates which, however, are only local in time. This approach leads to a new notion of momentum which is transported by the flow and serves as a generalization of vorticity. We also discuss how delta distribution momenta lead to "vortex-solitons", also called "landmarks" in imaging science, and to new numeric approximations to fluids.

研究动机与目标

• 在微分同胚群上建立欧拉方程与 EPDiff 方程族测地线方程之间的严格联系。
• 分析算子 $ L_{\varepsilon,\eta} = (I - \frac{\eta^2}{p}\Delta)^p \circ (I - \frac{1}{\varepsilon^2}\nabla\circ\text{div}) $ 对流体方程的正则化效应。
• 证明当 $ \eta > 0 $ 且 $ p \geq (n+3)/2 $ 时，EPDiff 方程的解全局适定，从而可逼近欧拉方程。
• 引入广义动量概念，该动量由流场传输，扩展了涡度的概念。
• 探索从 $ \delta $-函数动量中涌现出的涡旋孤子（特征点）及其在数值流体逼近中的作用。

提出的方法

• 在 $ \text{Diff}_{H^\infty}(\mathbb{R}^n) $ 上使用右不变度量，其范数定义为 $ \|v\|_{\varepsilon,\eta}^2 = \int \langle L_{\varepsilon,\eta}v, v \rangle \, dx $，其中 $ L_{\varepsilon,\eta} $ 结合了平滑与散度校正功能。
• 推导 EPDiff 方程 $ \partial_t m = -(v\cdot\nabla)m - \text{div}(v)m - m\cdot(Dv)^T $，其中 $ m = L_{\varepsilon,\eta}v $，该方程为微分同胚群上的测地线方程。
• 证明动量 1-形式 $ \widetilde{m} $ 在流场 $ \varphi $ 下保持不变，即 $ \widetilde{m}(\cdot,t) = \varphi(\cdot,t)_* \widetilde{m}(\cdot,0) $，意味着支集与衰减性质的保持。
• 建立局部时间估计，将 $ \varepsilon, \eta > 0 $ 时 EPDiff 方程的解与欧拉方程解在 $ \varepsilon, \eta \to 0 $ 时的联系。
• 利用拉格朗日形式 $ \partial_t \varphi = K^{\varphi} \ast (\varphi_* \widetilde{m}(\cdot,0)) $，其中 $ K $ 是 $ L_{\varepsilon,\eta} $ 的格林函数，以初始动量表示流场。
• 通过 $ \delta $-函数动量对特征点测地线进行数值探索，显示其根据角动量与能量表现出排斥或吸引动力学。

实验结果

研究问题

• RQ1 当 $ \varepsilon, \eta \to 0 $ 时，带有正则化参数 $ \varepsilon, \eta > 0 $ 的 EPDiff 方程解如何收敛到欧拉方程的解？
• RQ2 广义动量 $ m = L_{\varepsilon,\eta}v $ 在流场中保持结构的作用是什么？它如何推广涡度概念？
• RQ3 在何种条件下 EPDiff 方程存在全局解？这与欧拉方程的适定性有何关联？
• RQ4 如何从 $ \delta $-分布动量中产生涡旋孤子（特征点）？它们表现出何种动力学行为？
• RQ5 EPDiff 框架能否作为欧拉方程的稳定数值逼近，特别是在 $ \eta \to 0 $ 的极限下？

主要发现

• 当 $ \varepsilon, \eta \to 0 $ 时，欧拉方程的解是 EPDiff 方程解的极限，且已建立时间局部范围内的收敛性界。
• 当 $ \eta > 0 $ 且 $ p \geq (n+3)/2 $ 时，EPDiff 方程全局适定，对所有时间存在唯一解，确保了正则化逼近的稳定性。
• 动量 1-形式 $ \widetilde{m} $ 在流场 $ \varphi $ 下保持不变，因此紧支集与快速衰减性质得以保持，表明结构稳定性。
• 涡旋孤子（特征点）由 $ \delta $-分布动量产生，其动力学行为取决于能量与角动量，表现为散射或捕获，临界能量阈值为 $ E = \frac{8}{5}|\omega|^2 $。
• 由核 $ K_{0,1} $ 在原点附近的导数导出的向量场在 $ (x_1,x_2) $-平面上表现出螺旋流线，并沿 $ x_3 $-轴向外流动，表明具有局域涡旋环状行为。
• 对特征点测地线的数值积分显示，改变背景动量 $ \overline{m} $ 会导致排斥与吸引相互作用，表明其在模拟流体不稳定性中的潜在应用价值。