[论文解读] Research of convergence of the iterative method for solution of the Cauchy problem for the Navier - Stokes equations based on estimated formula

A. Tsionskiy, Mikhail Tsionskiy|arXiv (Cornell University)|Jan 10, 2011

Iterative Methods for Nonlinear Equations被引用 6

一句话总结

本文研究了在初始条件（柯西问题）下，利用傅里叶变换和拉普拉斯变换求解二维和三维纳维-斯托克斯方程的迭代方法的收敛性。建立了估计收敛边界公式 $ F_n/( u u^4) < 1 $，在多个周向模态上进行了数值验证，并通过速度幅值比较确认了快速收敛，所有测试参数下均满足 $ |υ_2^*| \ll |υ_1| $。

ABSTRACT

Solution of the Navier-Stokes equations with initial conditions (Cauchy problem) for 2D and 3D cases was obtained in the convergence series form by iterative method using Fourier and Laplace transforms in paper $\cite{TT02}$. For several combinations of problem parameters numerical results were obtained and presented as graphs. Estimated formula for the border of the parameter area of convergence of the iterative method was obtained in paper $\cite{TT03}$. This paper describes numerical proof of convergence of the iterative method for solution of the Cauchy problem for the Navier - Stokes equations. Usage of estimated formula for the border of the parameter area of convergence of the iterative method is shown for wide ranges of the problem's parameters.

研究动机与目标

• 分析二维和三维纳维-斯托克斯方程柯西问题中迭代方法的收敛行为。
• 验证先前工作中推导出的迭代方法收敛边界的估计公式。
• 在不同粘性系数、力幅值和径向模态参数下，数值展示迭代方法的有效性和鲁棒性。
• 通过第一、二次迭代间速度幅值比较，提供收敛性的可视化与定量证据。
• 确认在固定 $ \nu $ 和 $ F_n $ 时，随着 $ \mu $ 值增大，收敛性得到改善。

提出的方法

• 在极坐标系下，对初始速度为零且施加径向力的纳维-斯托克斯方程应用迭代方法。
• 利用傅里叶变换和拉普拉斯变换，推导出二维和三维情况下柯西问题的级数形式解析解。
• 推导出估计收敛边界 $ F_n/( u \mu^4) < 1 $，适用于所有周向模态 $ n = 1,2,3,\dots $。
• 对 $ n = 1 $ 至 $ 5 $ 进行数值模拟，变化 $ F_n $、$ \mu $ 和 $ \nu $，其中 $ F_n $ 设为 $ 10^k/n $，$ k = 0,1,2,3 $。
• 在 $ \nu = 0.01 $ 及其他取值下，以对数尺度图形比较第一轮迭代 ($ |\vec{u}_1| $) 和第二轮迭代 ($ |\vec{u}_2^*| $) 的速度幅值。
• 在 $ F_n $-$ \mu $、$ F_n $-速度和 $ \nu $-速度平面上绘制图表，收敛区域定义为曲线 $ F_n/( u \mu^4) = 1 $ 之上的区域。

实验结果

研究问题

• RQ1 估计收敛边界 $ F_n/( u \mu^4) < 1 $ 是否能准确预测纳维-斯托克斯方程柯西问题中迭代方法的收敛区域？
• RQ2 在不同周向模态和参数组合下，该方法的第一轮与第二轮迭代中速度幅值如何演化？
• RQ3 该迭代方法在宽范围的粘性系数 $ \nu $、力幅值 $ F_n $ 和径向模态参数 $ \mu $ 下是否具有鲁棒性？
• RQ4 $ |\vec{u}_1| $ 与 $ |\vec{u}_2^*| $ 之间是否存在定量关系？在所有测试配置中是否均满足 $ |\vec{u}_2^*| \ll |\vec{u}_1| $？
• RQ5 增大 $ \mu $ 对收敛速率有何影响？在固定 $ F_n $ 和 $ \nu $ 时，更高 $ \mu $ 值是否能改善收敛性？

主要发现

• 该迭代方法在所有测试的周向模态 $ n = 1 $ 至 $ 5 $ 中均表现出强收敛性，所有情况下均满足 $ |\vec{u}_2^*| \ll |\vec{u}_1| $。
• 所有测试的参数组合（$ F_n $、$ \mu $、$ \nu $）中，当 $ 0.2 \leq F_n \leq 1000 $ 且 $ 2.0 \leq \mu \leq 18.0 $ 时，均落在由 $ F_n/( u \mu^4) < 1 $ 定义的收敛区域内。
• 对于 $ \nu = 0.01 $，所有数值结果均以黑点形式绘制在曲线 $ F_n/( u \mu^4) = 1 $ 之上，证实了估计边界的正确性。
• 第二轮迭代的速度幅值 $ |\vec{u}_2^*| $ 显著小于第一轮迭代的 $ |\vec{u}_1| $，表明其快速衰减，收敛性显著。
• 随着 $ \mu $ 增大，收敛性得到改善，表明在固定 $ F_n $ 和 $ \nu $ 时，更高的 $ \mu $ 值可增强收敛速率。
• 在 $ \varphi $-平面上 $ 0 \leq r \leq 3 $ 的图形比较结果表明，所有半径和角度下均满足 $ |\vec{u}_2^*| < |\vec{u}_1| $，放大图突出显示了第二步的低幅值行为。