[论文解读] Implementation for blow up of tornado-type solutions for complex version of 3D Navier-Stokes system

Максим Дмитриевич Арнольд, А. В. Хохлов|ArXiv.org|Jun 26, 2008

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

一句话总结

本文通过使用纯虚数初始数据的傅里叶变换方程，数值实现了复3D纳维-斯托克斯系统中龙卷风型解的自相似爆破情景。通过在(0,0,R)处构造紧支撑的缩放厄米函数线性组合作为初始条件，作者模拟了在临界时间 $T_{\text{cr}}$ 附近呈现多项式爆破行为的能量增长，观测到的增速为15–20阶，接近西蒙和李所预测的理论渐近行为 $1/(T_{\text{cr}}-t)^5$。

ABSTRACT

We consider Cauchy problem for Fourier transformation of 3-dimensional Navier-Stokes system with zero external force. Using initial data purposed by Dong Li and Ya.G.Sinai we implement self-similar regime producing fast growing behavior of the energy of solution while time tends to critical value.

研究动机与目标

• 数值实现西蒙和李提出的复3D纳维-斯托克斯系统中龙卷风型解的理论爆破情景。
• 通过在(0,0,R)处使用至多 $D=3$ 阶的缩放厄米函数线性组合，构造傅里叶空间中紧支撑的初始数据，以实现自相似爆破动力学的数值模拟。
• 研究解在趋近有限临界时间 $T_{\text{cr}}$ 时的能量增长速率，检验理论预测的 $1/(T_{\text{cr}}-t)^5$ 渐近行为。
• 验证爆破行为在离散化参数和初始条件形式的小扰动下的鲁棒性。

提出的方法

• 研究采用零外力、粘性系数设为1的傅里叶变换3D纳维-斯托克斯系统，通过幅度参数 $A$ 的幂级数展开求解积分方程 (1)。
• 初始数据构造为 $v_A(k,0) = A c_0(k)$，其中 $c_0(k)$ 是以 (0,0,R) 为中心、半径为 $r$ 的球体内紧支撑的函数，使用至多 $D=3$ 阶的缩放厄米函数线性组合。
• 在 $nX=nY=160$、$nZ=1440$ 的3D结构化网格上，以 $dt=0.001$ 的时间步长进行数值演化，使用 Intel MKL 进行3D卷积计算。
• 能量 $\mathcal{E}(t) = \int |v(k,t)|^2 dk$ 作为 $L^2$ 范数的积分平方根，于连续时间步长中进行监测。
• 通过 $g_p(y,s) = h_p((0,0,R)+\sqrt{R}y,s)$ 的重整化方法分析自相似行为，关注递推关系谱分析中的不稳定模态，极限为 $p \to \infty$。
• 通过将时间步长、空间分辨率和初始条件形式最多扰动25%，测试爆破的稳定性，观察到在自然缩放下一致的爆破行为。

实验结果

研究问题

• RQ1 能否通过傅里叶变换方程和紧支撑初始数据，实现3D纳维-斯托克斯系统中龙卷风型解的自相似爆破情景的数值实现？
• RQ2 在该类模拟中，当 $t$ 趋近临界时间 $T_{\text{cr}}$ 时，能量增长速率 $\mathcal{E}(t)$ 的观测值是多少？
• RQ3 爆破行为在离散化参数和初始条件形式的小扰动下是否具有鲁棒性？
• RQ4 数值模拟是否表现出与西蒙和李预测的理论 $1/(T_{\text{cr}}-t)^5$ 渐近行为一致的能量增长？
• RQ5 递推关系 $h_p(k,s)$ 中的不稳定特征值模态在爆破动力学出现过程中起什么作用？

主要发现

• 解的能量（以 $\mathcal{M}(t) = \sqrt{\mathcal{E}(t)}$ 衡量）从 $t=0.000$ 时的1.00000 初始下降至 $t=0.042$ 时的0.3664，随后急剧上升至 $t=0.045$ 时的19.972，$t=0.046$ 时的101.20，以及 $t=0.047$ 时的698.57，表明在 $T_{\text{cr}}$ 附近存在快速增长。
• 在另一组离散化参数（空间步长0.20328，$dt=0.0008$）下，能量从 $t=0.0424$ 时的0.3047 上升至 $t=0.0480$ 时的1027.7，表明在参数变化下仍保持一致的爆破行为。
• 观测到的能量增长速率估计在15至20之间，显著高于理论的 $1/(T_{\text{cr}}-t)^5$ 速率，但在计算限制下仍处于预期范围内。
• 在时间步长、空间分辨率和初始条件形式最多25%的扰动下，爆破行为保持稳定，且解结构在自然缩放下发生变换。
• 数值解表现出清晰的空间结构：在傅里叶空间中呈现明显的“云团”——第一团位于 $z \sim 20$，第二团位于 $z \sim 40$，第三团位于 $z \sim 60$，与 $h_p$ 支撑在 $O(Rp)$ 处的自相似缩放一致。
• 结果支持存在一个10参数族的初始数据可导致有限时间爆破，能量在 $t \to T_{\text{cr}}$ 时发散，尽管由于分辨率限制，尚未完全实现精确的 $1/(T_{\text{cr}}-t)^5$ 速率。