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$,这产生了斯托克斯第一和第二问题。

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

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

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

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

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

考虑半径为 $R$ 的圆管中稳态、充分发展、轴对称的流动,由轴向 $x$ 上的均匀压力梯度 $dp/dx < 0$ 驱动。在柱坐标系 $(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)$。

简单库埃特流。 顶部平板在 $x$ 方向以速度 $U$ 移动;底部平板静止;$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}.$

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

相似变量 $\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代码的标准验证案例。
• 杰弗里-哈梅尔流。 半角为 $\alpha$ 的楔形中的稳态纯径向流 $u_r(r,\theta)$。流函数 $\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))$ 中。精确解没有解决这个问题,因为它们居住在解空间的子空间中,在那里完全非线性动力学从未参与。

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

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

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

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

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

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

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