最怪异的流体行为,任意小的粘性,就会产生巨大的能量耗散

康托的天堂 2023-05-28 14:01:56

纳维-斯托克斯方程(Navier-Stokes Equations)是流体力学中的基本方程,描述了粘性流体的运动。这组方程由法国科学家克劳德·路易·马里·纳维埃和英国科学家乔治·斯托克斯在19世纪初独立推导出来,可以看作是欧拉方程在考虑粘性效应后的延伸。

欧拉方程是描述理想流体(即没有粘性的流体)运动的基本的方程,而纳维-斯托克斯方程则是描述有粘性的流体运动的基本方程,以下简记为N-S方程。这些方程在科学和工程上都是重要的,然而人们对它们的了解还很不够,它们是对数学的重大挑战。

为了描述这些方程,我们在欧几里得空间R^d中进行研究,这里d表示维数,是2或者3。设在点x=(x₁,…,xd)∈R^d处,以及时间t∈R时刻,流体具有速度向量

流体的压强则是

欧拉方程就是对于所有(x,t),

而N-S方程则是对于所有(x,t),

这里ν>0是一个摩擦系数,称为流体的“粘性”。

本文中限于讨论不可压缩流体,就是说,除了要求适合方程(1)或(2)以外,还要求对于一切(x,t)适合

欧拉方程和N-S方程只不过就是把牛顿定律F= ma应用于流体的无穷小部分而已。事实上,很容易看到向量

就是流体的微元在位置x时刻t所经历的加速度。

在欧拉方程中,牛顿运动定律F= ma里的力F完全来自压强梯度,举例来说,如果压强随高度增加,则负梯度

指向下方,所以有一个净力驱使流体向下运动),但在N-S方程(2)中还有一个附加项

来自摩擦力。

N-S方程在许多各种各样的情况下都与对实际流体进行的实验十分符合。因为流体是很重要的,所以N-S方程也是很重要的。

欧拉方程只不过是N-S方程在v=0时的极限情况,但是我们将会看到,二者的解的性态很不相同,哪怕是ν很小时也是。

我们想要了解在初始条件,即对于所有x∈R^d,

下,欧拉方程和(3)的解或者N-S方程和(3)的解,这里u⁰(x)是已给的初速度。为

了与(3)式相容,假设对于所有x∈R^4,

还有为了避免在物理上不合理的情况,例如无穷能量,还假设当

u⁰(x)以及具有固定的t的u(r,t)都“足够快”地趋于零。这里不解释“足够快”的确切意义是什么,但是从现在起就假设只讨论这种急速衰减的速度。

物理学家和工程师想要知道怎样有效而又快速地计算N-S方程(2)和(4)的解,并理解这些解的性态。数学家首先要问的是是否有解存在,如果有,又是否只有一个解。虽然欧拉方程已经有了300年的历史,N-S方程的历史也接近了200年,专家们对于N-S方程或者欧拉方程的解是否对于所有时间都存在,或者会发生“破裂",仍然没有共识。有严格证明支持的肯定的答案似乎还很遥远。

让我们把关于欧拉方程和N-S方程的“破裂问题”讲得更明确一点。方程(1)-(3)涉及了u(x,t)的一阶和二阶导数。假设(4)中的初始速度u⁰(x)具有一切阶数的导数

而且当x趋于无穷时,它们都“足够快”地趋于零,这些都是很自然的。于是我们要问,是否N-S方程(2)-(4)或者欧拉方程(1),(3)和(4)对于所有的x∈R^d和t>0都有解u(x,t)和p(x,t),使得对于所有的x∈R^d和t>0,导数

都存在(而且当x趋于无穷时,都“足够快”地趋于零)。具有这种性质的一对u和p称为欧拉方程或N-S方程的“光滑解”。谁也不知道(在3维情况下)这种解是否存在。我们知道有某个依赖于初始速度(4)的正的时刻T=T(u⁰)>0,使得欧拉方程或N-S方程有“光滑解”在x∈R^d,t∈[0,T)中存在。

在2维空间情况,可以取T=+∞,就是说对于2D欧拉或2D N-S不会发生“破裂”。在3维空间的情况,谁也不能排除有这样的可能:对于某个如上所述的T=T(u°),可以在区域

中找到欧拉方程或N-S方程的解u(r,t),p(x,t),它在此区域中光滑,但是有一个导数使得

这将意味着在越过时刻T以后就不会有光滑解了,我们就说3D欧拉或3D N-S在时刻T破裂。说不定对于3D欧拉以及/或者3D N-S真就发生了这样的事情,谁也不知道相信什么才好。

对于3D N-S以及3D欧拉,人们已经做出过不少计算机仿真。纳维-斯托克斯的仿真并没有显示出有破裂的证据,但是这可能只意味着会产生破裂的初速度u⁰极为罕见。3D欧拉的形态可能极为狂野,所以很难判断一项给定的数值研究是否表示发生了破裂。事实上,要完成一项可靠的3D欧拉的数值仿真,是难得出了名的事。

假设会发生破裂,并且研究这时N-S方程或欧拉方程解的性态,这是很有用的。例如设在时刻T<∞对于3D欧拉发生了破裂,Beale, Kato和Majda有一个定理断言,当t→T时“涡度”

会增长得这样快,使得积分

发散。这件事被用来指明某些声称证明了3D欧拉发生了破裂的计算机仿真其实是无效的。现在也知道当t接近一个有限的破裂时刻T时,涡度向量的方向会随着x的变化而狂野地变动。

(5)式的向量ω有很自然的物理意义:它表示在时刻t流体怎样绕着x点旋转。如果在x点处,在时刻t,放一个小涡轮,使它的轴的定向平行于w(x,t),则这个涡轮会被流体以角速度ω旋转起来。

对于3D N-S,Sverak最近有一个结果说,如果发生了破裂,则压强p(x,t)上下方均为无界。

1930年代,由勒雷首创地提出了一个非常有前途的思想:研究N-S方程的“弱解”。他的思想如下:初看起来,N-S方程(2)和(3)只有当u(x,t)和p(x,t)为充分光滑时才有意义,例如,我们可能希望u对x_j具有2阶导数,因为这种导数出现在(2)式的右方。然而,如果做形式的运算就会发现,(2)和(3)等价于下面的(2')和(3'),而在那里,即令u(x,t),p(x,t)很粗糙,也是有意义的。让我们先来看一下(2')和(3')是怎样导出来的,然后再来看它们的用处。

起点是以下的观察:R^n上的一个函数F恒等于零当且仅当对于每一个光滑的且在某一紧集合外恒为0的函数θ都有

把这一点说明用于方程(2)和(3),再作一点形式运算(即分部积分)就知道它们等价于

以及

更准确地说,给定任意的光滑的函数u(x,t)和p(x,t),方程(2)和(3)成立的充分必要条件是(2'),(3)对于任意光滑而且在R^3×(0,∞)的某个紧集合外恒为0的函数

O₁(x,t),θ₂(x,t),θ₃(x,t)和φ(x,t)

均成立。

称函数θ₁(x,t),θ₂(x,t),θ₃(x,t)和φ(x,t)为试验函数,而说u和p是3D N-S的弱解。因为在(2'),(3')中,所有的导数都是对于光滑的试验函数取的,所以它们对于很粗糙的函数u(x,t)和p(x,t)仍是有意义的。总结起来,我们有如下的结论:

一对光滑的(u,p)解出了3D N-S当且仅当它们是弱解,然而弱解的思想即令对于粗糙的(u,p)也是有意义的。

我们希望能按以下的步骤来应用弱解:

证明3D N-S在整个R^3×(0,∞)上有适当的弱解存在。

证明3D N-S的任意适当的弱解必定是光滑的。

由此得到结论,在步骤(1)中做出的适当弱解就是3D N-S在整个R^3×(0,∞)上的光滑解。

在这里“适当的”就是指“不太大的”,它的精确定义我们就略过不说了。

上面的计划对于某些有趣的偏微分方程是成功的,但是对于3D N-S,这个计划只部分地实现了,很久以来人们就知道如何构造出3D N-S的适当的弱解,但是其唯一性一直没有得到证明。感谢Sheffer,Lin, Caffarelli, Kohn和Nirenberg等人的工作,现在已经知道了,在一个具有小的分形维数的集合E之外,3D N-S的任一适当的弱解都是光滑的,即必定具有任意阶的导数。特别是E中不会包含曲线。要想排除破裂,就必须证明E是空集合。

对于欧拉方程,弱解仍是有意义的,但是Sheffer和Schnirelman举出的例子表明,它们的性态可能很奇怪。一个2维的流体在开始可能是平静的,后来不受任何外力的作用却突然在一个有界的空间区域里运动起来,然后又恢复平静。对于2D欧拉的弱解,这种行为可能会发生的。

除了上面讲过了的破裂问题以外,纳维-斯托克斯和欧拉方程还会引起一些其他的基本问题。设固定了3D N-S或3D欧拉的初始值u⁰(x),则在时刻t=0的能量是

对于v≥0,用

来表示以u⁰(x)为初速度,粘性为v的3D N-S的解。假设

在一切时间都存在(至少当v>0时如此)。则在t≥0时,

的能量是

以(1)-(3)为基础作一些初等的运算即得

特别是对于欧拉方程有v=0,因此(6)式给出,能量等于一个常数E_0,而只要解是存在的,它就与时间无关。

现在假设v很小,但不是零。由(6)式自然会猜想,当v很小时,

所以在很长的时间段里几乎是常数。然而数值试验和物理实验都强烈地指出并不如此。相反,似乎存在一个依赖于u⁰但是与v无关的T_0>0,使得流体的初始能量到了时刻T_0至少损失一半,而不论v多么小(当然要大于0)。

如果能够证明或否定这个结论是很重要的。我们需要理解,为什么这么小的粘性会产生如此大的能量耗散。

37 阅读:4594
评论列表
  • 2023-11-25 22:58

    兄弟们能看懂吗?我是有点看不懂

    全能 回复:
    nvxhjhdrdcv 回复:
    没什么问题,主要说了以下3点...
  • 2024-03-19 19:56

    天书

  • 2023-08-13 21:04

    字母我都没认全[笑着哭][笑着哭][笑着哭]

  • 2024-02-21 20:41

    你跟我聊这个,那我得睡一会[笑着哭]

  • 2024-03-20 08:34

    这个方程可能只是一个近似方程,自然界中不存在严格解。

  • 2023-06-06 21:45

    很好

  • 2023-06-14 09:35

    就跟三体运动微分方程一下解不出来通解的,所以飞机只能靠风洞吹和软件模拟

    黄恒 回复:
    风洞吹的是战机风阻系数
  • 2023-11-26 04:12

    不知这篇文章是抄来的,还是自己写的,抄来就没意思了。有这个兴趣和爱好自己写,就值得鼓励,等你把这个问题解决了,你就是大师级的人物了,加油!

  • 2024-03-19 21:47

    脑瓜子嗡嗡的[不开心]

  • yrnd 2
    2024-04-11 08:27

    韦东奕的研究方向。

  • 2024-04-15 23:50

    99.9999999%都看不懂

  • 2023-11-28 20:21

    是韦神写的么?

    Dirichlet 回复:
    韦神写不出这种
  • 2024-03-31 02:16

    数学太高深,勉强大学读完高数,再深点反正我是整不懂了。

  • 2024-04-17 20:46

    元代码0和1可以解释一切科学公式[点赞][点赞]

康托的天堂

简介:科学如此美妙,我想让你知道