N-S纳维斯托克斯方程推导过程

1、很多人一听到N-S方程就有点头皮发麻，因为涉及到流体力学的知识比较多，如果没有一个完整有逻辑的思路，理解N-S方程是有点困难。其中涉及到欧拉法，场论，随体导数，流体力学连续性方程(即质量守恒方程)，流体力学N-S方程(即动量方程)，动量方程在流体力学中有两种，一种是理想流体动量方程，一种是粘性流体动量方程，粘性流体的动量方程也叫纳维-斯托克斯方程，也简称N-S方程。我试图想把N-S方程弄清楚点，所以写了一点东西，分享一下。首先要讲一下流体力学的欧拉法，在课本中还讲了拉格朗斯法，因为连续性方程和N-S方程是用欧拉法得出的，和拉格朗日法没什么关系。我就不讲拉格朗日法，以免产生混乱。欧拉方法的着眼点

2、不是流体质点而是空间点。设法在空间中的每一点上描述出流体运动随时间的变化状况。如果每一点的流体运动都已知道，则整个流体的运动状况也就清楚了。欧拉方法中流体质点的运动规律数学上可表示为下列矢量形式：假设空间一点的坐标(x,y,z,t),其中x,y,z是该空间的坐标，t是此刻时间。u,v,w是这一空间点的三个方向速度。p,P,T是这一空间点的压力，密度和温度。这样就有了每一个点的速度，压力，密度，温度，就可以描述运动流体的状态。这里需要强调一点的是下面这六个式子，可以换一个角度把他们看成方程，对后面理解连续性方程和N-S方程有帮助，比如u=x+2y+3zu=u(x,y,z,t);v=v(x,y,z

3、,t);w=w(x,y,z,t);p=p(x,y,z,t);：=：(x,y,z,t);T=T(x,y,z,t)因为后面需要随体导数的概念，还需要把速度函数表示成矢量的形式。前面u,v,w是标量，是v在(x,y,z,t)直角坐标系三个方向的速度。、='.(r,t)随体导数表示流体质点在欧拉场内(见流体运动学)运动时所具有的物理量对时间的全、,上+、r口d.一导数。上面定义了空间一点速度，那么加速度就是。设有一流体质点在运动，t时刻在dtM点(x,y,z,t),速度为v(M,t),过了事之后，在M'点，速度为S(M：t+&)。根据定义，加速度表达式可以写成如下：d、.(M；

4、tt)-、(M,t)二limdt'Qt二lim、(M,t：t)i(M；t),ljm-(M,t)-(M,t)-tQt-t0t基于时间A的变化基于空间位置mtm的变化、(M,tt)-、(M,t)MM、(M,t)(M,t)=limlimlim.t-0;t.J。:tMM0MM=(M，。，lim'(M(M,t)4MMQMM.I,M而M靠近，MTM的变化会引起下个方向速度的变化用M点速度可以分解成这两部分，是因为从MtM'点，一方面有时间&t的变化，一方面有空间dt位置的变化，分解成这两部分，正是基于这两个原因。写成直角坐标系，用u,v,w三个方向速度表示成如下：u=u(x

5、,y,z,t);v=v(x,y,z,t);w=w(x,y,z,t);代入上面加速度公式，得到chc(M,t).(M；t)->(M,t)=、.lim=dt尊叫沸，+MM-&_,M?DM靠近，MT的变化会引起妻个万向速度的变化Mmh速度dufu(x,y,z,t)u(x,y,z,t)u(x,y,z,t)u(x,y,z,t)，v(x,y,z,t)dt二t；x二yu(x,y,z,t)w(x,y,z,t)：zFuFuFu；:u二u-：-v-：-wft；xFy；zdu；:u：:u：:u；:uu一vwdtft；:xfy.zdv.:v：:v：:v2Vuvw.;:zdtft三x：:ydw；:w：:w

6、：:w：:wu-v-wdtct二xcy二z至此已经用欧拉法推到出了流体速度和加速度(即随体导数)的公式。随体导数也可以用复合函数求导的方法得到。用复合函数链导法则会更容易理解一些。后面接下来要推导的是流体力学连续方程。连续性方程是质量守恒定律在流体力学中的具体表述形式。它的前提是对流：:udxu-fx2体采用连续介质模型，速度和密度都是空间坐标及时间的连续、可微函数。；:udxu-；x2假设有一个微体积正六面体，正六面体的中心三个方向的速度是u,v,w。一、士-:udx左表面的流速uM=u：x2右表面的流速uN=udx；x2单位时间内x方向流出和流进的质量流量差：M出一M进=iPu+1&quo

7、t;m)dxbydz,Pu&)dx'dydz=u)dxdydz2Fx2三xJEx同理y方向和z方向的质量流量差：(%)dxdydz：:y：(、w)dxdydz:z在dt时间内因为密度变化而减少的质量为：:dxdydz-(：二)dxdydz=-dxdydz:tFt由质量守恒，单位时间内流出与流入六面体的流体质量差综合应等于六面体因密度变化而减少的质量。I：，(：u)：:(Pv)2(5)-dxdydz=dxdydz二|x二y二z二t:uPvFPw八=0.:t：:xZ；z以上就已经得到了连续性方程。对不可压缩流体，连续性方程可以简化，可以得到以下简化的连续性方程：-：u：vWNc一二

8、0.x二y二z这个不可压缩流体的连续性方程很重要，下面推导N-S方程的时候要用到。接下来要推导出流体力学的N-S方程。在推导N-S方程之前，有很多人都在这里有困惑。这里有两个概念要搞清楚，那就是什么是理想流体和粘性流体。我们很多课本在讲流体力学的时候是先讲了理想流体的动量方程，之后又没有接着讲粘性流体的动量方程，所以有些人到后面再讲N-S方程就混淆了。另外就是很多人一听到N-S方程就心里有点害怕，畏惧了，还没来得及去仔细研究就放弃了，如果仔细研究一下，其实也不难，很多流体力学的书是用场论的知识去推导出N-S方程的，我们工科学校对场论没有接触，最好还是用正六面体的方法来推导N-S方程。哈工大陈卓

9、如和王洪杰老师的工程流体力学对N-S方程的推导用的是正六面体法，很容易看懂。清华大学的书就比较难，可以参考。在这里得先推到一下理想流体的动量方程，后面再推导粘性流体的动量方程。理想流体是一种不可压缩、不计粘性（粘这里必须先分清理想流体和粘性流体的概念。度为零）的流体。欧拉在忽略粘性的假定下，建立了描述理想流体运动的基本方程。实际上,理想流体在自然界中是不存在的，它只是真实流体的一种近似模型。但实际上由于流体中存在着粘性，流体的一部分机械能将不可逆地转化为热能，并使流体流动出现许多复杂现象，例如边界层效应、摩阻效应、非牛顿流动效应等。自然界中各种真实流体都是粘性流体。下ABFpdxp；x2y在理

10、想流体内部取一微体积正六面体。中心点压力p(x,y,z),受力分析沿x轴方向:面推导理想流体动量方程。理想流体和粘性流体的区别在于是否有粘性力，即切应力。；:pdxp-fx2T=0。1 .表面力，因为是理想流体，没有切应力,pdx左表面Pm=pmA=（p）dydz：x2右表面Pn=pnA=（pdx）dydz二x22 .质量力，单位质量力在三个坐标轴上分量是fx,fy,fz。x轴方向质量力fxPdxdydz。在x轴方向由牛顿第二定律：''二F'ma二/：:pdxpdxdu(p)dydz-(p)dydzfxdxdydz:-dxdydz一='fx2x2dtr1pdu由

11、前面速度随体导数fx-二xdt11甲二u二u二u二ufx-二一u一v一w一人:：x：x:v:z同样y轴z轴列方程，可以得到理想流体动量方程组。,1；:p；:u；:u；:u；:ufx-+u+V+W,:，：x：:t：:x；:y；:z工1jp.:V'VjV'Vfy-二一uVW；二yct二x二y二z12PFw；:wjwjwfz二uVW.-z二t二x二y二z卜面要推导粘性流体动量方程，也就是纳维斯托克斯方程，也叫做N-S方程。同样取一个微六面体，但粘性流体有切应力，分别对六个面做受力分析如图所示：x方向的受力，质量力，左右方向压力，前后面切力，上下面切力。由牛顿第二定律列方程:Afx：&

12、quot;dydzKydyxxYdx)dydz二yx-yxdxdz-(yxyydy)dxdzj%dydx-3x+Udz)dydx一:dxdydz-du对粘性流体，切应力由广义牛顿内摩擦定律确定:i/：uyx，zy，xy=:(一)二；y；xi/fv；z.-yz二()二；zzw；:u_-zx_（）一-xz：xjz粘性流体中某一点三个方向的压力是不相等的，任意点的压力与三个方向的正应力有以下关系式：1，、P=-(-xx；-yy,；-zz)；3cryy-zzu=p-2;:x二p”;yw=P-2J-；z以上广义牛顿内摩擦定律以及压力与正应力的关系可以找陈卓如老师的工程流体力学,关的解释。代入化简上面由牛

13、顿第二定律得到的方程：有相fx；dxdydz'卜'一，：yx一-yxdxdz-(yxILyxxxxdydz-仁dx)dydzexdy)dxdz-|izxdydx一(Ezx+-zxdz)dydx'1czJ-PdxdydzCdfx：.xdu一dt-,二yx,:；zx.:y：:z二P曲(1)dt把%=P-2J；yxxxyx二n.y二x=N（私十生）代入上式（1）中:二x二zfx-xx-yx-zx.:x:y.:zfx«-2噌）/：uN、口二（）-（：:yfx_2；:w；x:y:z色）ddtfx2-2mujx.2.x-2m二u.2-7上.x；y2:xFz-2m二ufx-2m二ujx-2二ufx-2m二ujx一2.x一2-y.2二u-2m:u-2m二u-2:x22vw,J_2_J.:xjy.:x：:z二工dt1fp+v?ufz_2-7j2u2:y-2mu一2.-z.:u+v?uftjvjw+一十L-L.二y二z由不可压缩流体连续性可以得到这部分为0.:u.:xdu）二代入1dt；:y代入刖随体导数1一P.x；：2u一2一x+v；：2u2-y+v；：2u一2:zFu一次u2221二p+v一一v+v二v+v二v_二vuP-:y-2一x-2-y-2;zFtfy-2二w-2二w-2二w1Fp-