Navier-Stokes精确解：Poiseuille流、Couette流与Stokes扩散

纳维-斯托克斯方程以非线性著称。那么怎么可能精确求解它们呢?

对称性。这就是全部的诀窍。

当流动的几何形状足够简单(直管、两个平板、无限平面)时,速度只能指向一个方向,并沿一个或两个坐标变化。在许多这些对称设置中,非线性项要么消失,要么大幅简化,使得方程坍缩为线性偏微分方程。在稳态情况下,你通常只剩下一个可以用纸笔求解的常微分方程。

这是直觉理解。纳维-斯托克斯方程描述所有可能的流体运动。但如果你强迫流体进入一个非常有序的情况,没有旋涡,没有混沌,一切都朝一个方向前进,方程的大部分复杂性就变得无关紧要了。纳维-斯托克斯的难点在于流体自我推动的反馈回路。在这些对称流动中,没有什么可以相互推动。非线性自平流项就直接消失了。

这些精确解不是好奇的产物。它们是流体力学教育的基础,是数值代码的基准,也是理解真实流动何时以及如何失控的起点。

纳维-斯托克斯方程的精确解是满足以下条件的速度场 $u$ 和压力 $p$:

$$\partial_t u + (u \cdot \nabla)u = -\nabla p + \nu \Delta u, \qquad \nabla \cdot u = 0$$

以闭形式,无需数值近似。

关键机制是非线性平流项 $(u \cdot \nabla)u$ 的消失或简化。对于 平行流(在笛卡尔坐标系中 $u = (u(y,t),\, 0,\, 0)$ 或在柱坐标系中 $u = (0,\, 0,\, u(r,t))$ 的流动),速度场自动满足无散度条件,平流项完全消失,因为速度在流向方向上没有分量,也没有沿流向方向的变化。动量方程随后简化为线性偏微分方程,或对于稳态流动,简化为常微分方程。

当流动也是 稳态($\partial_t u = 0$)时,剩下的是

$$\nu \Delta u = \nabla p,$$

一个泊松方程,其解是粘性流动的经典精确解:泊肃叶流、库埃特流及其相关流动。

对于非稳态平行流($\partial_t u \neq 0$ 但平流仍然消失),方程变为扩散方程 $\partial_t u = \nu\, \partial_{yy} u$,这产生了斯托克斯第一和第二问题。

如果你在寻找真正的闭式Navier-Stokes解，下面这些经典解是最好的起点。以下每个公式都在特定几何和边界假设下求解了不可压缩方程的一个简化版本。

重要的警告是:这些并不是一般解。它们之所以精确,是因为几何设置消除了使三维Navier-Stokes问题困难的大部分非线性行为。

下表列出了不可压缩牛顿Navier-Stokes方程在无滑移边界条件下的标准闭式解。这里$\mu$是动力黏度,$\nu=\mu/\rho$是运动黏度,$R$是管道半径,$h$是平板间距,$dp/dx$是常数流向压力梯度。

这些公式是研究生流体力学教材中的标准内容,例如Batchelor的An Introduction to Fluid Dynamics、Landau和Lifshitz的Fluid Mechanics以及Schlichting的Boundary-Layer Theory。下面的推导直接展示这些简化过程,而不是把公式当作引用结果。

泊肃叶流(也称为 哈根-泊肃叶流)是纳维-斯托克斯方程最重要的精确解,也是大多数工程师首先遇到的解。

想象水稳定地流过一根长直管道。两端之间的恒定压差驱动流动。管壁不动,因此接触管壁的流体停滞在零速度(这是无滑移条件)。离管壁越远,流体速度越快。中心移动最快。

速度分布?一条 抛物线。在管壁处为零。在中心处最大。中间是平滑的曲线。剖开管道看横截面:它是一个倒置的碗。

关键的定量事实是:总流量与管道半径的 四次方 成比例。四次方。半径翻倍,你得到的不是两倍的流量,甚至不是四倍。你得到十六倍的流量。这就是哈根-泊肃叶定律,它解释了为什么即使是微小的动脉狭窄也能切断血液供应。

假设: 泊肃叶流假设流体是不可压缩的牛顿流体(恒定粘度),流动是稳态的且充分发展的(不再从入口处加速),并且流动是层流的。平滑。有序。实际上,管道流在雷诺数约为2,300时转变为湍流,这是一个经验观察,没有人成功从底层理论推导出来。

考虑半径为 $R$ 的圆管中稳态、充分发展、轴对称的流动,由均匀压力梯度 $dp/dx < 0$ 沿轴向 $x$ 驱动。在柱坐标系 $(r, \theta, x)$ 中,速度场的形式为 $u = (0,\, 0,\, u(r))$。

不可压缩性条件 $\nabla \cdot u = 0$ 自动满足。平流项消失,因为 $u$ 不依赖于 $x$。轴向动量方程简化为

$$0 = -\frac{dp}{dx} + \mu\left(\frac{d^2 u}{dr^2} + \frac{1}{r}\frac{du}{dr}\right),$$

或等价地

$$\frac{1}{r}\frac{d}{dr}\left(r\frac{du}{dr}\right) = \frac{1}{\mu}\frac{dp}{dx}.$$

由于 $dp/dx$ 是常数,这是关于 $r$ 的常微分方程。用边界条件 $u(R) = 0$(无滑移)和 $du/dr|_{r=0} = 0$(对称性)积分两次,得到 抛物线速度分布:

$$u(r) = \frac{1}{4\mu}\left(-\frac{dp}{dx}\right)(R^2 - r^2).$$

中心线处的最大速度为 $u_{\max} = \frac{R^2}{4\mu}\left(-\frac{dp}{dx}\right)$,平均速度为 $\bar{u} = u_{\max}/2$。

对横截面积分得到体积流量的 哈根-泊肃叶定律:

$$Q = \frac{\pi R^4 \Delta p}{8 \mu L},$$

其中 $\Delta p > 0$ 是管道长度 $L$ 上的压降。

$R^4$ 依赖关系是决定性特征:半径的小变化产生流量的大变化。在血液动力学中,这就是为什么血流对动脉狭窄如此敏感。

有效性: 该解适用于不可压缩、牛顿、稳态、充分发展的层流。实验上,湍流转变发生在 $\text{Re} = \bar{u} D / \nu \approx 2{,}300$,其中 $D = 2R$ 是管道直径。这个阈值是经验观察;没有定理可以从纳维-斯托克斯方程预测它。关于雷诺数和层流-湍流转变的讨论,见 雷诺数、湍流以及为何小尺度重要。

库埃特流 是当一块平板移动而另一块保持静止时,被困在两块平行平板之间的流体的精确解。它是流体力学中最简单的精确解之一。

想象一副扑克牌平放在桌子上。拖动顶部的牌侧向移动,下面的牌也会移动,每张牌的移动都比上面的牌少一点。就是这样。速度从底部平板的零线性变化到顶部平板的速度,仅此而已。

直线速度分布。底部平板:静止。顶部平板:移动。中间的一切只是线性插值,不需要压力梯度,不需要复杂的设置,运动纯粹由移动边界拖动流体产生。

当你还沿着通道施加压力梯度时,事情变得更有趣,因为这时你在同一间隙中结合了剪切驱动流和压力驱动流。速度分布变形为叠加在线性分布上的抛物线,有时称为 平面泊肃叶-库埃特流,根据压力梯度相对于平板速度的强度,你甚至可以在一个壁面附近得到回流。

完全去掉移动平板,保持两壁静止,让压差完成所有工作。这就是 平面泊肃叶流,管道流的平板类比。抛物线。中间最快。两壁为零。

考虑两块相距 $h$ 的无限平行平板之间的稳态流动。设 $y$ 为垂直于平板的坐标,底部平板在 $y = 0$,顶部平板在 $y = h$。流动是平行的:$u = (u(y),\, 0,\, 0)$。

简单库埃特流。 顶部平板以速度 $U$ 沿 $x$ 方向移动;底部平板静止;$dp/dx = 0$。动量方程简化为 $d^2u/dy^2 = 0$,得到

$$u(y) = U\frac{y}{h}.$$

这是纳维-斯托克斯方程最简单的非平凡精确解:完全由边界条件驱动的线性速度分布。

平面泊肃叶流。 两块平板都静止;恒定压力梯度 $dp/dx < 0$ 驱动流动。动量方程变为

$$\mu \frac{d^2 u}{dy^2} = \frac{dp}{dx},$$

边界条件为 $u(0) = u(h) = 0$。解为

$$u(y) = \frac{1}{2\mu}\left(-\frac{dp}{dx}\right)y(h - y),$$

关于中心线 $y = h/2$ 对称的抛物线分布。

一般库埃特-泊肃叶流。 结合移动顶板和压力梯度给出叠加

$$u(y) = U\frac{y}{h} + \frac{1}{2\mu}\left(-\frac{dp}{dx}\right)y(h - y).$$

第一项是剪切驱动的贡献;第二项是压力驱动的贡献。根据 $dp/dx$ 相对于 $\mu U / h^2$ 的符号和大小,分布可以是单调的,具有内部极大值,甚至在一个壁面附近表现出回流。

泊肃叶流和库埃特流是稳态的。没有随时间变化。斯托克斯问题 是最简单的非稳态精确解,它们揭示了一些美妙的东西:粘度使动量在流体中 扩散,就像热量在固体中扩散一样。

斯托克斯第一问题(也称为 瑞利问题):想象一个平板上方的巨大静止流体池。在时间零时,平板突然开始以恒定速度侧向滑动。紧邻平板的流体立即被拖动,但更远处的流体需要时间才能注意到。一个平滑的边界层从平板向外生长,随着时间推移变得越来越厚。

平板上方任何高度的速度取决于该高度与特征扩散长度 $\sqrt{\nu t}$ 的比值,其中 $\nu$ 是粘度,$t$ 是经过的时间。更粘的流体?运动向上传播得更快。

斯托克斯第二问题:相同的设置,但现在平板正弦振荡而不是以恒定速度移动。振荡只渗透到流体中的有限距离。更远的地方,流体几乎注意不到。运动的幅度随高度呈指数衰减,产生一个薄的振荡边界层。这是振荡边界层背后的机制:振荡平板使附近的流体运动,但扰动随离平板的距离呈指数衰减。

两个斯托克斯问题都涉及半无限流体($y > 0$),它位于 $y = 0$ 处的平板上方。流动是平行的:$u = (u(y,t),\, 0,\, 0)$。平流项消失,控制方程是一维扩散方程:

$$\frac{\partial u}{\partial t} = \nu \frac{\partial^2 u}{\partial y^2}.$$

斯托克斯第一问题(瑞利问题)。 初始条件:$u(y,0) = 0$,其中 $y > 0$。边界条件:$u(0,t) = U$,其中 $t > 0$。远场:$u \to 0$(当 $y \to \infty$)。

相似变量 $\eta = y / (2\sqrt{\nu t})$ 将偏微分方程简化为常微分方程。解为

$$u(y,t) = U\,\operatorname{erfc}\!\left(\frac{y}{2\sqrt{\nu t}}\right),$$

其中 $\operatorname{erfc}$ 是余误差函数。边界层厚度以 $\delta \sim \sqrt{\nu t}$ 增长,这是扩散传播的标志。

斯托克斯第二问题。 平板振荡:$u(0,t) = U\cos(\omega t)$。寻求形式为 $u(y,t) = \operatorname{Re}[\hat{u}(y)\, e^{i\omega t}]$ 的解得到

$$\hat{u}(y) = U\exp\!\left(-(1+i)\frac{y}{\delta_s}\right), \qquad \delta_s = \sqrt{\frac{2\nu}{\omega}},$$

因此

$$u(y,t) = U\exp\!\left(-\frac{y}{\delta_s}\right)\cos\!\left(\omega t - \frac{y}{\delta_s}\right).$$

幅度随穿透深度 $\delta_s$ 呈指数衰减,相位线性滞后:一个横波,其能量完全被粘度耗散。更高的频率穿透深度更小。

泊肃叶流、库埃特流和斯托克斯流获得了最多的关注。但它们并不是唯一的精确解。远不止这些。

• 泰勒-格林涡:二维中具有真正涡旋结构的衰减旋涡模式。测试过CFD代码?你可能已经用它运行过了。这是每个人首先想到的基准,几十年来一直如此。
• 杰弗里-哈梅尔流:在楔形通道中收敛或发散的流动。它捕捉了流体如何加速进入变窄的间隙或减速进入扩大的间隙。
• 希门茨驻点流:流体正面撞击平壁,在表面减速到零并向侧向转向。风击中建筑物。射流撞击板。

对称性。这些中的每一个都利用特定的几何对称性使方程易于处理,它们在专门的情况下很重要,但对于日常的流体力学入门,泊肃叶流和库埃特流仍然承担所有的重任。

除了平行流,其他几个精确解族利用特定的对称性或自相似结构:

• 泰勒-格林涡。 在二维中,$u = (A\cos(ax)\sin(by),\, -B\sin(ax)\cos(by))$ 满足 $aA = bB$(不可压缩性)和指数时间衰减 $e^{-\nu(a^2+b^2)t}$。这是完整二维纳维-斯托克斯方程(包括非线性项,它变成梯度并被吸收到压力中)的精确解。它是DNS代码的标准验证案例。
• 杰弗里-哈梅尔流。 稳态纯径向流 $u_r(r,\theta)$,位于半角为 $\alpha$ 的楔形中。流函数 $\psi = \nu f(\theta)$ 满足关于 $\theta$ 的三阶非线性常微分方程。收敛和发散通道都存在解,在较高雷诺数下具有丰富的分岔结构。
• 希门茨驻点流。 流体冲击平壁的二维相似解。流函数 $\psi = \sqrt{a\nu}\, x\, f(\eta)$,$\eta = y\sqrt{a/\nu}$,将纳维-斯托克斯简化为希门茨常微分方程:$f''' + ff'' - f'^2 + 1 = 0$,边界条件为 $f(0)=f'(0)=0$,$f'(\infty)=1$。

我们可以在所有这些情况下精确求解纳维-斯托克斯方程。那么为什么仍然存在百万美元的开放问题?

技巧。每一个都依赖于技巧。几何形状被如此仔细地选择,使得方程中最困难的部分(非线性项)要么完全消失,要么简化为可管理的东西,问题之所以可解,正是因为它已被剥离了使纳维-斯托克斯困难的一切。管道流?一维的。库埃特?一条线。泰勒-格林涡把它的非线性隐藏在压力里面。

千禧年大奖问题 询问的是 一般 三维流动。没有对称性技巧。没有简化几何。光滑无散度初始数据,所有尺度上的完全非线性相互作用。在那种情况下,没有人证明解总是保持光滑,也没有人证明它们会爆破。我们真的不知道。

所以这些精确解告诉我们一些东西,但远远不够。它们证明了当你给方程强对称性作为支撑时,方程可以产生显式光滑解。百万美元的问题是光滑性是否 总是 成立,对于标准三维形式中的任意光滑无散度初始数据,还是在湍流的完全暴力中的某个地方,某些东西会灾难性地且不可逆转地出错。

上面讨论的每个精确解都通过消除或平凡化非线性平流项 $(u \cdot \nabla)u$ 来实现可处理性。在平行流中它完全消失。在泰勒-格林涡中它是被吸收到压力中的梯度。在相似解中它通过变量变换简化为低维问题。

克雷千禧年大奖问题涉及具有任意光滑、快速衰减初始数据的三维不可压缩纳维-斯托克斯方程的 初值问题,正是这些简化都不适用的情况。

问题是以下方程的解

$$\partial_t u + (u \cdot \nabla)u = -\nabla p + \nu \Delta u, \quad \nabla \cdot u = 0, \quad u(x,0)=u_0(x)$$

在 $\mathbb{R}^3$ 上,给定光滑、无散度、快速衰减的初始数据,是否对所有时间保持在 $C^\infty(\mathbb{R}^3 \times [0,\infty))$ 中。精确解没有解决这个问题,因为它们居住在解空间的子空间中,在那里完全非线性动力学从未参与。

简而言之:精确解在高度对称的情况下提供显式光滑解。开放问题是适定性是否扩展到无限制的三维流动。关于精确的表述,见 纳维-斯托克斯存在性与光滑性。

要理解方程本身以及每个项的含义,从 什么是纳维-斯托克斯方程? 开始

要了解这些方程如何从第一原理建立,阅读 纳维-斯托克斯方程如何推导。

要理解像泊肃叶流和库埃特流这样的层流何时分解为湍流,阅读 雷诺数、湍流以及为何小尺度重要。

要理解精确解无法回答的百万美元问题,阅读 千禧年大奖问题:存在性与光滑性。

推荐在本站阅读的下一页:

• 什么是纳维-斯托克斯方程?(控制系统及其项)
• 纳维-斯托克斯方程如何推导(从动量平衡和牛顿本构定律到偏微分方程)
• 雷诺数、湍流以及为何小尺度重要(从这里描述的层流状态过渡到完全湍流)
• 千禧年大奖问题:存在性与光滑性(精确解无法解决的开放问题的精确表述)