流体力学课件--.pptVIP

* C.H.5 粘性流体动力学 5.1 应力分析 1.应力 与理想流体一样，考虑粘性后，流体所承受的力也可归结为两类：即质量力与表面力。我们称单位面积上的表面力为应力，用 表示， 为该面的外法线方向。一般说来， 可分解为法向应力与切向应力，即： 对于理想流体因无摩擦存在，切应力为零，只存在法应力，且总是指向内法线方向，即 。而考虑粘性后， 不再为零。 在直角坐标系中，来研究分别作用于垂直于x, y, z轴各平面上的应力 。 设 为作用于正方向上的表面力(yz平面，左边流体作用于x面上的力)， 在x,y,z方向的投影为： 同理可作出y,z面上的表面力 在x,y,z方向的投影： 法向应力： 切向应力： 二阶应力张量： y z x o 2. 应力的性质 1)切应力之间具有对称性，即： 2)某点的应力不仅与该点的位置、作用的时间有关，而且与作用面的方向有关。即同一点取不同的面， 是不同的。但可以证明，任意三个相互垂直的表面力的法向应力之和相等，即 三个主应力值之和与作用面的方向无关，因此，仅是作用点的坐标和时间函数。把主应力的算术平均值定义为粘性流体的压力，即： 上式的负号是由于规定法向应力以拉力为正，压力为负. 3)应力张量与变形速度的关系 考虑粘性后，引入6个未知项，为使方程组封闭，在一定的假设下建立应力张量与变形速度间的关系： 若流体不可压： ，上式可简化为： 对于理想流体， 4) 关于粘性系数 当 时，理想流体;关于 的项是由于考虑了粘性而引起的，称 为粘性系数。 在水平方向（x方向）以力F拉面积为A的板，板的速度为U，在速度不大和板矩h都很小的情况下，试验证明下式成立： 取极限时有： 称作牛顿切应力公式或牛顿定律，满足上式牛顿流体；不满足上式非牛顿流体。 粘性系数是反应流体粘性大小的一种度量，它的大小与流体的种类有关，而且随温度与压强而变。 的量纲可据牛顿切应力公式导出 N·s/m2=牛顿·秒/平方米,动力学粘性系数 ; m2/s= 平方米/秒 ,运动学粘性系数 当温度T=15时，水和空气的年性系数分别为： 5) 粘性流体压力与理想流体压力之间区别 对于理想流体，不存在切向应力。两个特性: ①理想流体压力P的大小与它的作用面无关，即： ②理想流体的压力就是作用在物体表面上沿内法线方向上的单位面积上的表面应力，即： 对于粘性流体，存在切应力，特性: ①作用于物体表面沿法线方向的表面压力是法向应力 ，而不是粘性流体压力P。 ② 不仅与点的位置、时间有关，而且与作用面的方向有关 ③定义粘性流体的压力P为三个主应力的算术平均值，即： 5.2 Naiver-Stokes方程 1. 粘性流体的运动方程 B A C E D F H G z y x dx dz dy D’ 1)表面力：作用于六面体上的表面力在x,y,z方向的分量为： 2)质量力：单位质量流体所受质量力： 据牛顿第二定律，便可写出粘性流体的运动方程如下： 把应力和变形速度之间的关系代入粘性流体的运动方程中，得到Ｎ－Ｓ方程 Naiver-Stokes运动方程 若流体不可压，则上式可简化成： 未知数共有五个，方程只有四个, 方程组不封闭，但对于均质不可压流体或正压流体, 方程组封闭 不可压流体在直角坐标系中的运动方程及应力关系总结为： 5.3 N-S方程的几个解析解 均质不可压或正压流体，N-S方程与连续方程构成了一完整方程组。给出定解条件便可进行求解。但由于该方程组是一个二阶非线性偏微分方程组，求解困难。 理想均质不可压或正压流体的欧拉方程中也存在非线性项，但部分实际问题可以看作是无旋，于是问题归结为求解势函数满足的拉氏方程。 粘性流体，运动有旋，必须去解原始的二阶偏微分方程组。数学家们还没有提供求解非线性偏微分方程的普遍有效方法，就是对于粘性不可压均质流体方程组解的存在性与唯一性问题迄今还没有全面解决，这是对于某些简单情形才得到证明。 所以为求出N-S方程的解。我们不得不据力学考虑，作一些近似或假定，对方程组进行简化，以便找出一定准确度的解。 一是找准确解:对于简单问题，忽略非线性项，方程组简化，求方程的解析解。 二是求近似解：据问题的特点，抓大放小，得到近似方程，然后求解，具体处理上还可分两种情况: （1）小Re情形： （2）大Re情形： (3)中等Re情形：惯性力与粘性力同样重要，数值解。 忽略惯性力，方程线性化；