流体力学-N-S方程.pptVIP

实际流体的运动微分方程;以应力表示的黏性流体运动微分方程式; 粘性流体的法向应力和切向应力都必须同时考虑。在粘性流体表面上任取一点N，过N作微元面积ΔA。其外法线方向矢量为 ，切线方向为 ，N点的表面应力分为法向应力pn和切向应力τ，pn和τ随微元面积ΔA在空间的位置而变化。在直角坐标系中将pn和τ沿x，y，z三个坐标轴分解成9个应力分量，即 ; 下面来分析控制体所受表面力的合力。为了简单起见，以y方向为例。设流体微元的密度为ρ，则微元质量为ρdxdydz，其表面应力在y轴上分量如图1所示。 y轴上合力为：;对于x、z轴同理有 （3） 方程（3）就是以应力表示的粘性流体运动微分方程，通常X、Y、Z作为已知量，不可压缩流体 已知，方程应包含六个应力及三个速度分量，共9个未知数。而方程（3）加上连续性方程也只有4个方程，无法求解，必须找出新的补充关系式。 ;应力和变形速度的关系;联系到流体微元运动分析，如图2所示，则流体微元绕z轴的剪切角速度为; 流体运动时在内部产生的切应力是阻碍流体剪切变形的。其实流体的粘性不只对剪切变形有阻碍作用。下面我们分析一下法向应力与直线变形速度之间的关系。 如图3所示，ABCD经过dt时间变成ABCD时，微元体在x方向的粘性阻力阻止其拉长。微元体在y方向缩短，亦必遇到y方向的粘性阻力阻止其缩短。 这种表面法向作用的应力的大小与各该方向的直线应变速度有关，方向与直线变形的方向相反。; 由式（6）可以看出，由于各个方向的直线应变速度不见得相等，因而这种由于粘性阻碍作用所产生的法向应力也是各向不等的，pxxpyypzz统称为一点上的各项异性压强。 于是在实际流体运动时，一点上的法向应力除了由于分子运动统计平均的各向同性压强p之外，还需加上由于粘性影响而与直线变形有关的各向异性压强，最后可以得到法向应力与直线应变速度之间的关系为;将（5）和（7）联合在一起，写成矩阵形式，则为; 对于不可压缩流体来说，将（7）式中的三式相加，因??;p的含义; 纳维－斯托克斯方程; 为简化计算，我们再引入一个数学符号拉普拉斯（Laplace）算子，或标志性微分算子，为;式中，;理想流体运动微分方程及其积分;这就是流体的平衡微分方程，也叫欧拉平衡方程。 所以，N-S方程是不可压缩均质流体运动的普遍方程。如果将牛顿第二定律理解为力的平衡关系式，则从上面的推导不难看出N-S方程也是作用在流体上的力的平衡关系式。 ; 该方程组有四个未知数p,vx,vy,vz，它和连续性方程共有四个方程式，从理论上讲，在一定初始条件和边界条件下，任何一个不可压缩均质粘性流体的运动问题，是可以求解的。但是由于实际流体中的粘性影响非常复杂，单纯用求解N-S方程的方法去解决各种实际问题是有困难的。 而且N-S方程式二阶非线性非齐次的偏微分方程组，除针对具体情况用数值计算方法外，还不能积分求普遍解，只有在某些简单的或特殊的条件下，才能求得精确解。 ;非线性偏微分方程的常用解法;非线性偏微分方程的常用解法