首页/文章/ 详情

说说纳维尔-斯托克斯方程(Navier-Stokes Equations)

1年前浏览2131

      流体力学(Fluid Mechanics)有着非常漫长的历史,最古老的涉及流体力学的人物可能就是古希腊的阿基米德了(Archimedes),大家都知道的阿基米德浮力定理,当时阿基米德在羊皮纸上用希腊语写了论浮力On Floating Bodies)的文章并流传至今。但是在很长时间里流体力学并不被当作一门独立的学科,直到1687年牛顿(Isaac Newton)在其划时代的巨著《自然哲学的数学原理》(Philosophiae Naturalis Principia Mathematica)之后流体力学才真正走上历史舞台。在书中除了著名的牛顿三大定理以外,牛顿还提到流体力学相关的内容,他写到:对于平直的均匀流体,流体层与层之间的剪应力正比于垂直流体方向上的速度梯度。

阿基米德的手稿

牛顿的《自然哲学的数学原理》

        除了阿基米德和牛顿以外,在流体力学的发展历程中,做出贡献的人物可以说是阵容豪华,比如非粘性流体(Inviscid flow)数学分析涉及的人物有欧拉(Leonhard Euler), 达朗贝尔(Jean le Rond d'Alembert),拉格朗日(Joseph Louis Lagrange), 拉普拉斯(Pierre-Simon Laplace), 泊松(Siméon Denis Poisson);非粘性流体的数学分析涉及的人物有高斯(Guass),泊松,圣维南(Saint-Venant) 等等,其它在流体力学中涉及湍流/紊流,粘性/非粘性的推动者有雷诺(Osborne Reynolds),泰勒(Geoffrey Ingram Taylor)等等。而对这些理论归纳整理的两位科学家就是纳维尔(Claude-Louis Navier)和斯托克斯(George Gabriel Stokes),当今有名的纳维尔-斯托克斯方程组(Navier-Stokes Equations)就是以他们俩的名字命名的,也是今天要讲的主题


       克劳德-路易·纳维尔Claude-Louis Navier,1785年2月10日-1836年8月21日)是一位法国工程师与物理学家,主要贡献在力学领域。

       纳维尔的父亲在1793年过世,纳维尔的母亲让她的叔叔Emiland Gauthey来教育纳维尔。1824年纳维尔进入法国科学院,1830年时担任国立桥路学校的教授,次年接替奥古斯丁·路易·柯西,在巴黎综合理工学院担任微积分及力学的教授。

       纳维尔在1821年将弹性理论以数学公式的形式表示,使这个领域第一次可以计算有足够精确度的结果。1826年纳维尔确认弹性模量是材料的一个基本属性,和物体的截面二次轴矩无关,因此纳维尔也是结构分析的创始者之一。

       当然,纳维尔主要的贡献还是纳维尔-斯托克斯方程,是流体力学的理论中心。由于其卓越的贡献,纳维尔也是法国艾菲尔铁塔上所刻的72人名字之一(还有拉格朗日,拉普拉斯,傅里叶等大家)。

纳维尔


埃菲尔铁塔上刻有象征法兰西科学精神的72个人物名字的位置


       乔治·加布里埃尔·斯托克斯爵士eorge Gabriel Stokes,1819年8月13日-1903年2月1日),爱尔兰数学家物理学家,就读和任教于剑桥大学,主要贡献在流体动力学、光学数学物理学(如大家熟知的斯托克斯公式)。他曾任英国皇家学会秘书和会长。

       他出生在新基 督家庭,父亲是爱尔兰斯莱戈郡教区牧师。先后在都柏林布里斯托尔就读,1837年考入剑桥大学彭布罗克学院四年后以最高分毕业,并获得史密斯奖。1849年获卢卡斯数学教授席位,1854年出任皇家学会秘书,1885–1890年期间出任会长。1889年被封为从男爵。1899年6月1日,他任卢卡斯教授50周年,剑桥大学举行了盛大庆祝会,校监向他颁发金牌。

       1842至1843年,他刊出第一批论文,关于不可压缩流体的稳定流动。其后在1845年,他就流体流动的摩擦力弹性固体的平衡和运动发表论文,在1850年探讨流体的内部摩擦力对运动的影响。他也曾就声音的理论作出贡献,如对声音强度的影响。他的研究标志着流体动力学的新里程,不但有助解释自然现象(如空中的运动、水中的运动等),更有助解决技术问题,如水在河和管道中的流动以及船只的表面阻力等。


        前面我们曾说到,在CFD流体力学中,最重要的方程就是纳维尔-斯托克斯方程(Navier-Stokes equations),这并不是一个单一的方程,而是一个方程组,是在众多科学家和工程师的推动下产生的。在流体力学中,有很多方程,但很多方程都和纳维尔-斯托克斯方程有着联系,可以说,该方程在流体力学中起着基础性的作用,但也起着决定性的作用,套用指环王的故事,Navier-Stokes equations可以比做至尊魔戒,其它大大小小的方程比做各个流体力学国度国王拥有的普通戒指,那么可谓是一戒统领众戒(one ring rules all rings)。  

       数学家、科学家伊恩·斯图尔特(Ian Stewart)出过一本书,名叫《17 Equations That Changed The World(改变世界的17个方程)》。其中大多数公式我们都见过,里面就包括纳维尔-斯托克斯方程。

改变世界的17个公式

       Navier-Stokes equations已经渗透到描述流体的方方面面,也在众多实际工程中得到了应用。但意外的是,这个方程的数学特性——解的存在性和光滑性至今没有得到证明

       纳维-斯托克斯存在性与光滑性是有关纳维-斯托克斯方程其解的数学性质有关的数学问题。方程可以描述空间中流体液体气体)的运动。纳维尔-斯托克斯方程式的解可以用到许多实际应用的领域中。不过对于纳维-斯托克斯方程式解的理论研究仍然不足,尤其纳维-斯托克斯方程式的解常会包括紊流。虽然紊流在科学及工程中非常的重要,不过紊流仍是未解决的物理学问题之一。

      许多纳维尔-斯托克斯方程式解的基本性质都尚未被证明。例如数学家就尚未证明在三维座标,特定的初始条件下,纳维尔-斯托克斯方程式是否有符合光滑性的解。也尚未证明若这样的解存在时,其动能有其上下界,这就是“纳维-斯托克斯存在性与光滑性”问题

       由于了解纳维尔-斯托克斯方程式被视为是了解难以捉摸的紊流现象的第一步,克雷数学研究所在2000年5月提供了一百万美元的奖金给第一个证明该方程的人。这也是美国克雷数学研究所Clay Mathematics Institute)在2000年提出的7个千禧年大奖难题中的问题之一克雷数学研究所设定了该问题具体的数学描述:

       证明或反证下的叙述:在三维的空间及时间下,给定一起始的速度场,存在一向量的速度场及纯量的压强场,为纳维-斯托克斯方程式的解,其中速度场及压强场需满足光滑及全局定义的特性。


我们不妨看看该问题具体的描述(以下方程的描述来自维 基百科):

        以数学的观点来看,纳维尔-斯托克斯方程是一个针对任意维度向量场的非线性偏微分方程。从物理及工程的观点来看,纳维-斯托克斯方程是一个用连续介质力学描述液体或非稀疏气体运动的方程式组。此方程式是以牛顿第二运动定律为基础,考虑一黏滞性牛顿流体的所有受力,包括压强、黏滞力及外界的体积力。

       由于克雷数学研究所提出的问题是在三维空间下,不可压缩的匀质流体为准,以下也只考虑此条件下的纳维-斯托克斯方程。

       令为描述流体速度的三维向量场,且为流体压强。纳维尔-斯托克斯方程为:

其中

  • 为动黏滞度

  • 为外力

  • 梯度运算子

  • 拉普拉斯算子,也可写为

       上述方程是向量方程,可以分解为三个纯量的方程,将速度及外力分解为三个座标下的分量:

       则纳维-斯托克斯方程可写成以下的形式,

       其中的未知数有速度及压强。由于只考虑三维空间,因此有三个方程及四个未知数,分别是速度的三个分量及压强,还需要一个方程才能解出所有的未知数。这个新增的方程是描述流体不可压缩性连续性方程式:

       所以纳维-斯托克斯方程解的速度会是无散度的向量函数。对于在均匀介质中的无散度流,其密度及动黏滞度为定值。具体方程的展开这里不再说明。


目前纳维尔-斯托克斯方程的证明进展如下:

  1. 二维空间下的纳维尔-斯托克斯问题已在1960年代已经被证明:存在光滑及全局定义解的解。

  2. 在初速相当小时此问题也已得证:存在光滑及全局定义解的解。

  3. 若给定一初速,且存在一有限、依而变动的时间T,使得在的范围内,纳维-斯托克斯方程有平滑的解,还无法确定在时间超过T后,是否仍存在平滑的解。

  4. 数学家让·勒雷在1934年时证明了所谓纳维-斯托克斯问题弱解的存在,此解在平均值上满足纳维-斯托克斯问题,但无法在每一点上满足。


       当然,虽然说纳维尔-斯托克斯方程描述了流体领域的大部分条件,不过也有其适用范围,该方程只适用于牛顿流体,什么是牛顿流体呢?简单说就是:任一点上的剪应力都同剪切变形速率呈线性函数关系的流体。一般高黏度的流体是不满足这种关系的,说明牛顿流体和非牛顿流体有个简单的例子就是大家熟知的虹吸现象。在低黏度下,虹吸要进行下去,吸取口必须在页面以下,但非牛顿流体的高黏度流体下,吸取口哪怕高于液面,其虹吸依然能够进行,因为黏度太大了,这个还是很容易想象到的。

       那么,可能有人要问了,居然有不适用纳维尔-斯托克斯方程的流体,那么该方程是不是就不完备了。是的,方程是不完备,不能一统流体世界的天下,不过对于工程应用来说,大部分情况还是处理牛顿流体,或者可以近似为牛顿流体,这样纳维尔-斯托克斯方程基本就是流体世界的王者了。

低黏度下的虹吸现象(吸取口低于液面)

高黏度下的虹吸现象(吸取口高于液面)


        说到这里,或许可以告一段落了,但或许有人也有点困惑,有些人会想:虽然我数学没那么厉害,但感觉纳维尔-斯托克斯方程的证明对于数学家来说应该不难吧,这么几百年过去了难道都没有一个数学家站出来证明他吗。我们已经征服了描述时空的相对论,让我们翱翔到宇宙深处;我们也已经很大程度上征服了量子力学所描述的微观世界,让我们体验到神奇的量子世界,那里有宇宙大爆炸,有浪漫的平行宇宙。可是,对于一个函数来说,其解的存在性和光滑性(可以看做高中数学中函数的连续性)都不能被证明,实在匪夷所思啊。

       有些数学问题看似简单,但对数学家来说,要证明其绝对的正确性,却并不一定是个易事,比如前面提到的美国克雷数学研究所提出的7个千禧年难题,其起草者之一安德鲁·怀尔斯(Andrew Wiles),他在1993证明了难倒数学界300多年的费马大定理(Fermat's Last Theorem),而费马大定理本身看似也是非常的简单,其问题描述如下:

当整数n >2时,关于x, y, z的方程 x^n + y^n = z^n 没有正整数解

       就是这样一个人人都能看懂的问题,却耗费了数学界300多年才得以解决,可想而知,对于数学证明来说,远没有我们想象的简单。


       希望,在我们这个时代,纳维尔-斯托克斯方程(Navier-Stokes Equations)能够得到证明,为我们打开自然界的奥秘!!!



Navier-Stokes的背景知识

简介 NS方程,全称:纳维-斯托克斯(Navier-Stokes)方程 ,2000年5月24日,美国克莱数学研究所的科学顾问委员会把NS方程列为七个“千禧难题”(又称世界七大数学难题)之一,这七道问题被研究所认为是“重要的经典问题,经许多年仍未解决。”克莱数学研究所的董事会决定建立七百万美元的大奖基金,每个“千年大奖问题”的解决都可获得百万美元的奖励。另外六个“千年大奖问题”分别是: NP完全问题, 霍奇猜想(Hodge),黎曼假设(Riemann),杨-米尔斯理论(Yang-Mills),庞加莱猜想和BSD猜想(Birch and Swinnerton-Dyer)。


NS方程的存在性与光滑性
  起伏的波浪跟随着我们的正在湖中蜿蜒穿梭的小船,湍急的气流跟随着我们的现代喷气式飞机的飞行。数学家和物理学家深信,无论是微风还是湍流,都可以通过理解NS方程的解,来对它们进行解释和预言。虽然这些方程是19世纪写下的,我们对它们的理解仍然极少。挑战在于对数学理论作出实质性的进展,使我们能解开隐藏在NS方程中的奥秘。
深度描述
  描述粘性不可压缩流体动量守恒的运动方程。简称N-S方程。因1821年由C.-L.-M.-H.纳维和1845年由G.G.斯托克斯分别导出而得名。在直角坐标系中,其矢量形式为=-Ñp+ρF+μΔv,式中ρ为流体密度,p为压强,u(u,v,w)为速度矢量,F(X,Y,Z)为作用于单位质量流体的彻体力,Ñ为哈密顿算子 ,Δ为拉普拉斯算子。后人在此基础上又导出适用于可压缩流体的N-S方程。N-S方程反映了粘性流体(又称真实流体)流动的基本力学规律,在流体力学中有十分重要的意义。它是一个非线性偏微分方程,求解非常困难和复杂,目前只有在某些十分简单的流动问题上能求得精确解;但在有些情况下,可以简化方程而得到近似解。例如当雷诺数Re1时,绕流物体边界层外,粘性力远小于惯性力,方程中粘性项可以忽略,N-S方程简化为理想流动中的欧拉方程(=-Ñp+ρF);而在边界层内,N-S方程又可简化为边界层方程,等等。在计算机问世和迅速发展以后,N-S方程的数值求解才有了很大的发展。
  在解释纳维-斯托克斯方程的细节之前,首先,必须对流体作几个假设。第一个是流体是连续的。这强调它不包含形成内部的空隙,例如,溶解的气体的气泡,而且它不包含雾状粒子的聚合。另一个必要的假设是所有涉及到的场,全部是可微的,例如压强,速度,密度,温度,等等。该方程从质量,动量,和能量的守恒的基本原理导出。对此,有时必须考虑一个有限的任意体积,称为控制体积,在其上这些原理很容易应用。该有限体积记为Omega,而其表面记为partialOmega。该控制体积可以在空间中固定,也可能随着流体运动。
其他数学难题的解决情况
  ①庞加莱猜想:2006年6月3日,中国数学家曹怀东和邱成桐已补足了佩雷尔曼证明的漏洞,将庞加莱猜想完全证明!
  ②黎曼假设:很多人攻关,没看到希望
  ③霍奇猜想:进展不大
  ④杨—米尔理论:太难,几乎没人做
  ⑤P与NP问题:没什么进展
  ⑥BSD猜想:有希望破解
相关人物
  斯托克斯,G。G(George Gabriel stokes1819~1903)英国力学家、数学家。1819年8月13日生于斯克林,1903年2月1日卒于剑桥。
  斯托克斯1849年起在剑桥大学任卢卡斯座教授,1851年当选皇家学会会员,1854年起任学会书记,30年后被选为皇家学会会长。斯托克斯为继I.牛顿之后任卢卡斯座教授、皇家学会书记、皇家学会会长这三项职务的第二个人。
  斯托克斯的主要贡献是对粘性流体运动规律的研究。C.-L.-M.-H.纳维从分子假设出发,将L.欧拉关于流体运动方程推广,1821年获得带有一个反映粘性的常数的运动方程。1845年斯托克斯从改用连续系统的力学模型和牛顿关于粘性流体的物理规律出发,在《论运动中流体的内摩擦理论和弹性体平衡和运动的理论》中给出粘性流体运动的基本方程组,其中含有两个常数,这组方程后称纳维-斯托克斯方程,它是流体力学中最基本的方程组。1851年,斯托克斯在《流体内摩擦对摆运动的影响》的研究报告中提出球体在粘性流体中作较慢运动时受到的阻力的计算公式,指明阻力与流速和粘滞系数成比例,这是关于阻力的斯托斯公式。斯托克斯发现流体表面波的非线性特征,其波速依赖于波幅,并首次用摄动方法处理了非线性波问题(1847)。
  斯托克斯对弹性力学也有研究,他指出各向同性弹性体中存在两种基本抗力,即体积压缩的抗力和对剪切的抗力,明确引入压缩刚度的剪切刚度(1845),证明弹性纵波是无旋容胀波,弹性横波是等容畸变波(1849)。
  斯托克斯在数学方面以场论中关于线积分和面积分之间的一个转换公式(斯托克斯公式)而闻名。
来源:WELSIM
非线性湍流光学理论爆炸材料Mathematica控制
著作权归作者所有,欢迎分享,未经许可,不得转载
首次发布时间:2023-06-24
最近编辑:1年前
WELSIM
一枚搞仿真的老员工
获赞 23粉丝 63文章 253课程 0
点赞
收藏
未登录
还没有评论
课程
培训
服务
行家
VIP会员 学习 福利任务 兑换礼品
下载APP
联系我们
帮助与反馈