N-S(纳维斯托克斯)方程推导过程.docxVIP

很多人一听到N-S 方程就有点头皮发麻，因为涉及到流体力学的知识比较多，如果没有一个完整有逻辑的思路，理解N-S 方程是有点困难。其中涉及到欧拉法，场论，随体导数， 流体力学连续性方程（即质量守恒方程），流体力学N-S 方程（即动量方程），动量方程在流体力学中有两种，一种是理想流体动量方程，一种是粘性流体动量方程，粘性流体的动量 方程也叫纳维-斯托克斯方程，也简称N-S 方程。我试图想把N-S 方程弄清楚点，所以写了一点东西，分享一下。 首先要讲一下流体力学的欧拉法，在课本中还讲了拉格朗斯法，因为连续性方程和N-S 方程是用欧拉法得出的，和拉格朗日法没什么关系。我就不讲拉格朗日法，以免产生混乱。欧拉方法的着眼点不是流体质点而是空间点。设法在空间中的每一点上描述出流体运动随时 间的变化状况。如果每一点的流体运动都已知道，则整个流体的运动状况也就清楚了。欧拉 方法中流体质点的运动规律数学上可表示为下列矢量形式：假设空间一点的坐标(x,y,z,t)， 其中 x,y,z 是该空间的坐标，t 是此刻时间。u,v,w 是这一空间点的三个方向速度。p, ? ,T 是这一空间点的压力，密度和温度。这样就有了每一个点的速度，压力，密度，温度，就可以描述运动流体的状态。这里需要强调一点的是下面这六个式子，可以换一个角度把他们看 成方程，对后面理解连续性方程和N-S 方程有帮助，比如u=x+2y+3z u ? u(x, y, z,t); v ? v(x, y, z, 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 是标量， 是?? 在(x,y,z,t)直角坐标系三个方向的速度。 ?? ???(r?, t) 随体导数表示流体质点在欧拉场内 随体导数表示流体质点在欧拉场内（见流体运动学）运动时所具有的物理量对时间的全 导数。上面定义了空间一点速度，那么加速度就是 d?? 。设有一流体质点在运动，t 时刻在 dt M 点（x,y,z,t），速度为??(M , t) ，过了?t 之后，在M ? 点，速度为??(M ?, t ? ?t) 。根据 定义，加速度表达式可以写成如下： 定义，加速度表达式可以写成如下： d?? ? lim ??(M ?, t ? ?t) ???(M , t) ? lim ??(M ?, t ? ?t) ???(M ?, t) ? lim ??(M ?, t) ???(M , t) dt ?t ?0 ?t ?t ?0 ?t ?t ?0 ?t ? lim ?? (M ?, t ? ?t) ?? (M , t) ? lim MM ? lim ? (M ?, t) ???(M , t) 基于时间?t的变化? ? 基于空间位置M ?M??的变化 ?t ?0 ?t ?t ?0 ?t MM ??0 MM ? ? ???(M , t) ?t ? ?? lim ??(M ?, t) ???(M , t) MM ??0 MM ? d?? dt ?t ?0,M和M ?靠近， M ?M ?的变化会引起??三个方向速度的变化用M点速度 可以分解成这两部分，是因为从M ? M ? 点，一方面有时间?t 的变化，一方面有空间 位置的变化，分解成这两部分，正是基于这两个原因。 位置的变化，分解成这两部分，正是基于这两个原因。 写成直角坐标系，用u,v,w 三个方向速度表示成如下： u u ? u(x, y, z,t); v ? v(x, y, z, t); w ? w(x, y, z, t); d? d?? ? dt ???(M , t) ?t ? ?? lim ??(M ?, t) ???(M , t) ? MM ??0 MM ? ?t ?0,M和M ?靠近， M ?M ?的变化会引起??三个方向速度的变化 用M点速度 du ? ?u(x, y, z, t) ? u(x, y, z, t) ?u(x, y, z, t) ? v(x, y, z, t) ?u(x, y, z, t) dt ?t ?x ?y ? w(x, y, z, t) ?u(x, y, z, t) ?z ? ?u ? u ?u ? v ?u ? w ?u ?t ?x ?y ?z du ? ?u ? u ?u ? v ?u ? w ?u ; dt ?t ?x ?y ?z dv ? ?v ? u ?v ? v ?v ? w ?v ; dt ?t ?x ?y ?z dw ? ?w ? u ?w ? v ?w ? w ?w dt ?t ?x ?y ?z 力学连续方程。连续性方程是质量守恒定律在流体力学中的具体表