[论文解读] Classical analytic solutions of the non-stationary Navier-Stokes equation in two, three and higher dimensions

R. K. M. Thambynayagam|arXiv (Cornell University)|Jul 29, 2013

Fluid Dynamics and Vibration Analysis参考文献 1被引用 13

一句话总结

本文提出了一种方法，通过构造初始速度场和外力，使得非定常纳维-斯托克斯方程中的非线性'脐带力'消失，从而推导出二维、三维及更高维度下的经典解析解。关键贡献在于：当初始数据和外力被选择为使脐带力为零并保持不可压缩性时，热方程的光滑解也能满足纳维-斯托克斯方程，从而在 $\mathbb{R}^n$ 中得到精确的、$C^{\infty}$ 解。该方法明确构造了 $n=2,3$ 的解，并验证了其满足能量边界和无散条件。

ABSTRACT

In this paper we present a method to derive classical solutions of the Navier-Stokes equations for non-stationary initial value problems in domain $\mathbb{R}^n$ ($n=2,3$ or higher). Exact solutions in $\mathbb{R}^2$ and $\mathbb{R}^3$ in the presence of an externally specified force are obtained by imposing certain restrictions on the initial data.

研究动机与目标

• 为 $\mathbb{R}^n$ 中 $n=2,3,$ 及更高维的非定常纳维-斯托克斯方程开发精确的经典解析解。
• 通过引入一种消除非线性'脐带力' $\mathcal{U}_i(v)$ 的条件，解决纳维-斯托克斯方程中的非线性挑战。
• 确保速度场在所有时刻保持散度为零（无散），满足不可压缩性约束。
• 证明在特定初始和外力条件下，热方程的光滑解也能满足完整的纳维-斯托克斯系统。
• 为 $\mathbb{R}^2$ 和 $\mathbb{R}^3$ 中的速度场和压力场提供显式、封闭形式的解，使其满足能量边界和正则性条件。

提出的方法

• 该方法引入了一个非线性项，即'脐带力' $\mathcal{U}_i(v)$，其源于压力梯度与速度非线性性；通过构造使其消失。
• 压力场通过简化的压力泊松方程（PPE）求解：$\Delta p = \rho \sum_{i=1}^n \frac{\partial (f_i - g_i)}{\partial x_i}$，利用伽马函数和距离核的基本解求解。
• 通过假设解具有可分形式：$v_i(x,t) = \phi_i(x) \Omega(t) e^{-\pi^2 \kappa t}$，其中 $\phi_i(x)$ 为满足不可压缩性的三角函数。
• 初始速度场被选为不同空间变量中正弦与余弦项的组合，以确保无散结构。
• 引入外力，使得非线性惯性项 $g_i$ 与压力梯度项相互平衡，从而实现 $\mathcal{U}_i(v) \equiv 0$，使系统简化为线性热方程。
• 振幅的时间演化由 $\Omega(t) = 1 + \int_0^t \mathcal{G}(\tau) e^{\pi^2 \kappa \tau} d\tau$ 控制，该式在特定外力情况下被显式求解。

实验结果

研究问题

• RQ1 在特定初始和外力条件下，是否可以在 $\mathbb{R}^n$（$n \geq 2$）中为非定常纳维-斯托克斯方程构造出精确的解析解？
• RQ2 在何种条件下，非线性'脐带力' $\mathcal{U}_i(v)$ 会消失，从而使纳维-斯托克斯系统退化为线性热方程？
• RQ3 如何设计初始速度场，以确保在 $\mathbb{R}^n$ 中所有时间下流动始终保持无散？
• RQ4 当非线性项被消除时，压力泊松方程在构造精确解中起到什么作用？
• RQ5 能否为二维和三维流动显式推导出满足能量边界和无散约束的光滑、$C^\infty$ 解？

主要发现

• 通过适当选择初始和外力条件，使脐带力 $\mathcal{U}_i(v) \equiv 0$，在 $\mathbb{R}^2$ 和 $\mathbb{R}^3$ 中导出了速度场和压力场的精确 $C^{\infty}$ 解。
• 对于 $f_1 = f_I e^{-\lambda t}$，$f_2 = f_3 = 0$ 的情形，$v_1$ 的解为 $v_1 = \{a\sin(\pi x_3) - c\cos(\pi x_2)\} e^{-\pi^2 \kappa t} + \begin{cases} f_I t, & \lambda = 0 \\ \frac{f_I}{\lambda}(1 - e^{-\lambda t}), & \lambda > 0 \end{cases}$，$v_2$、$v_3$ 和 $p$ 具有类似形式。
• 压力场显式给出为 $p = -\rho \Omega^2(t) e^{-2\pi^2 \kappa t} \left[ bc\cos(\pi x_1)\sin(\pi x_2) + ab\cos(\pi x_3)\sin(\pi x_1) + ac\cos(\pi x_2)\sin(\pi x_3) \right]$，满足泊松方程和光滑性。
• 由于初始数据被构造为满足 $\sum_{i=1}^n \partial v_i^0 / \partial x_i = 0$，且时间演化保持该性质，因此速度场在所有 $t \geq 0$ 下保持无散。
• 能量耗散受初始动能有界，且解满足 $\int_{\mathbb{R}^n} |v(x,t)|^2 dx < \mathcal{C}$ 对所有 $t \geq 0$ 成立，确认了物理一致性。
• 涡量定义良好，可通过速度场的旋度计算，确认了解在空间和时间上的光滑性与正则性。