粘性流体力学(4)

理论流体力学的计算结果与实验结果不一致的原因是在运动方程中没有将流 体的粘性考虑进去。计及流体粘性而由Navier与Stokes分别于1821及1845年 建立起来的Navier-Stokes方程，又因过于复杂而不可能求解。所以，流体力 学的发展遇到了困难。一方面是无粘流理论解决不了实际问题，另一方面是 粘流理论无法求解（蠕流有解，但没有实用价值）。边界层概念提出来以后， 把这两方面的困难都解决了。
第三章 粘性流体力学基本方程组 1、质量守恒定律：连续方程（质量方程）

按取和约定，上式可表示为

对于定常流，此式成为

对于不可压流，上式成为
第三章 粘性流体力学基本方程组 2、粘性流体的运动方程一动量守恒定律




为建立流体方程，须分析流 体微团受力。 微元体受到两种力：一为彻 体力，它是作用于微元体内 所有质量上的力，如重力、 惯性力、电磁力等；另一类 为表面力，它是作用在微团 界面上的力，如压力、摩擦 力等。 本节只研究表面力，用P表示 作用于单位容积的表面力。 在垂直于X轴的两个外表面上， 分别作用有合应力 :




将（3.9）带入（3.10）,有


（3.11）
(3.10)(3.11)就是著名的广义牛顿粘性应力公式。此式的适用范围很宽。

粘性流体的运动方程是动量守恒定律的数学表述。 对于流体运动应考虑彻体力和表面力两类外力：若F表示作用在单位质 量上的彻体力，P表示作用在单位容积上的表面力，由牛顿第二定律， 则运动方程可写为如下向量形式： （3.12）
常采用完全气体假设，即满足：
热流密度矢量q,以及对热传导系数，粘性系数 可采用公式：
利用以上的假设和关系，就可使可压缩粘性气体动力学方程组封闭
第三章 粘性流体力学基本方程组 5. NS方程的适用性





NS方程推导过程中引人了一些假设。方程组的适用性问题实质上就是 这些假设是否正确的问题。 所引入的假设：完全气体的状态方程、广义牛顿粘性应力公式、连续 介质假设。 连续介质的假设：只要最小涡的尺度大大超过分子平均自由行程，这 个假设就是可用的。气体分子的平均自由行程为0.0001mm,液体为 0.0000001mm。 NS方程是玻耳兹曼（Boltzmann）方程的第一次近似，即分子运动达到 近似平衡状态的玻耳兹曼方程。 总之，除了稀薄气体和个别特殊情况（研究激波厚度内的结构等特殊 问题），对于水和空气一类的流体，NS方程总是可用的。
层流底层
Rex=ux/
边界层的发展
x
流体流过光滑平板时，边界层由层流转变为湍流发生在 Rec=21053106
边界层
umax
充分发展的流动
u Le
进口段长度
(a)层流 u Le
边界层
充分发展的流动
进口段长度
(b)湍流 圆管内边界层的发展
层流时 Le/d=0.05Re； 湍流时 Le=4050d 这里，雷诺数 Re=ubd/。
第三章 粘性流体力学基本方程组
连续介质假设
第三章 粘性流体力学基本方程组



流体运动所遵循的规律是由物理学三大守恒定律规定的， 即质量守恒定律，动量守恒定律和能量守恒定律。这三大 定律对流体运动的数学描写就是流体动力学基本方程组。. 现在还求不出这个方程组的解析解，但研究这个方程组具 有极其重要的意义，因为千变万化的流动现象都是由这个 方程组所规定的。 本课程的全部内容，实质上就是在各种具体条件下、用各 种不同的方法、以不同的近似程度求解这个方程组，研究 解的性质，及其流动表现。

下标X表示应力向量作用在与x轴垂直的微元面上。

பைடு நூலகம்
可得作用在垂直于X轴的微元面上的表面力的合力为

同样可得作用在垂直于y轴和Z轴的微元面上的表面力的合力分别为

于是可得作用于单位容积的表面力的合力为

(3.1) 式中Px、Py、Pz都是向量，可以把它们沿三个坐标方向分解，即分解成正 应力和平行于各微元面的切应力 。
第三章 粘性流体力学基本方程组 1、质量守恒定律：连续方程（质量方程）

对于任一流体微元，根据质量守恒定律可得连续方程 （也称质量方程）：

此即三维可压缩流的连续方程，它表示：
净输出控制体 控制体内的质量 的质量流量 随时间的变化率 0
第四章 边界层理论 1、概述


粘性很小的流体以大 雷诺数运动时，在大 部分流场上可以略去 粘性作用，但在物面 附近的很薄一层的流 体内必须考虑粘性的 作用。这一个薄层流 体，称为边界层。 Prandtl (1904) 首先 提出了边界层概念， 边界层被称为:Mile stone of the Modern Fluid Mechanics。
(3.5)


在这种层状运动中v/x=0.所以此式右端的速度梯度实为应变变化率张量的 分量Syx的二倍, 即有：
(3.6)




式(3.5)(3.6)称为牛顿粘性应力公式， 斯托克斯将牛顿的这个公式推广到粘性流体的任意流动情形中去，其假设为 1．流体是牛顿流体，并是连续的；2．流体是各向同性的，即其性质与方向 无关；3．所建立的关系应适合运动和静止的情况。 考虑（3.3）式，假设更一般关系式：



将式（3.2)(3.3)带入（3.1）可得出单位容积的表面力公式：

(3.4)

定义一个反映流体运动与变形的对称张量，称为应变变化率张量

展开写为：

它是由性质不同的两组分量组成：主对角线上的分量描写线应变变化率i= ui/xi ，反映拉压；其余分量，描写角变形率，反映剪切。
第三章 粘性流体力学基本方程组
V 2 d 2 e d W1 W2 Q dt

其中左边项中分别表示系统内流体的动能和内能，W1是单位时间内边界 上应力对流体所作的功，W2是单位时间内体积力所作的功，Q是单位时间 内从边界面传入流体团的热量。
第三章 粘性流体力学基本方程组
其中




称为物质导数或随体导数，代表微团的某性质对时间的变化率。 将(3.3)（3.4)(3.11）带入(3.12)有:

（3.13）
这就是粘性流体的运动方程——著名的纳维一斯托克斯（NavierStokes）方程，简称NS方程。 一般可认与空间位置无关，整理上式，如第一式（x方向）,可得：



(3.14)

或写为矢量式：

对于不可压流体，也可写为：

由上诸式可见：尽管流体都有粘性，但只在速度梯度变化剧烈的 地方粘性应力才起重要作用。这一点很重要。

粘性流体和理想流体的差别， 以图为例：对于理想流体， 通过界面F，微元体A只对微 元体B作用了压力p；而对于 粘性流体，除正应力 y 外, 微元体A还对微元体B作用了 粘性切应力 力y 的大小也不等于压力p.
第三章 粘性流体力学基本方程组
连续介质假设




研究对象是流体的宏观运动，即大量分子的平均行为， 而不是单个分子的个别行为，因而可以不去考虑物质的分 子结构和单个分子的运动细节。 物质的分子结构和分子的热运动只对宏观运动存在间接的 影响，即只能通过影响物质的热力学特性来影响物体的运 动。 因此，当研究物体的变形、流动等宏观运动特性时，就可 以将物体作为一种连续体对待，而无须计及它的微观分子 结构。 分子动力论是用质点力学和统计学相结合的方法来研究物 质宏观力学和热力学性质的科学。这一理论取得了很大成 就，但它目前也只能应用于某些简单的气体，远不能解决 范围十分宽广的流体力学的大量问题。



物体的应力与运动学参数之间存在着一定的关系，在弹性 力学中这种关系是由胡克定律表示的，即弹性固体中应力 与应变成正比。 不同的流体有不同的性质，这种关系有不同的类型。对于 大多数流体, 应力与应变变化率成正比，或者说，应力与 应变变化率之间存在着线性关系，服从这种关系的流体称 为牛顿流体 牛顿根据实验最早提出了切应力与层间速度梯度成正比：
边界层概念的历史原因

十九世纪末，流体力学的研究工作有两个互不沟通的方向：一是理论流体力 学（亦称水动力学），用数学方法研究流体对固体物的绕流，当时已达到较 高的水平。但计算结果往往与实验结果不一致。例如圆柱体绕流，计算结果 是没有阻力，但实验表明有阻力。于是，在流体力学史上留下了达朗倍尔疑 题(d’Alembert’s Paradox)；二是水力学，主要是用实验方法进行研究，将实 验结果归纳成经验公式或半经验公式应用于工程实际，但缺乏理论基础。
yx
,而且正应
航空粘性流体 力学问题示例
航空粘性流体力学问题示例
航天粘性流体力学问题示例
航天粘性流体力学问题示例
第三章 粘性流体力学基本方程组 3、粘性流体的能量方程一能量守恒定律

能量方程反映了粘性流体在流动过程中满足的能量守恒定律。 当流体流动时，能量守恒定律可叙述为：封闭系统内流体能量随时间的 变化率等于单位时间内作用在该系统上所有外力所作的功和由边界传入 的热量，相应的数学表达式为
ij= 2sij
该式对于切应力的表示是符合牛顿粘性应力公式的. 但静止时，sij ＝0，按 上式， ij 也应为零，但流体静压p0 总是存在，所以上式需要修改

应在上式中再加一项，此项应只影响主对角线上的元素（因为静压是ij 对角线上正应力），由此，假设：
(3.7)

其中b为待定标量，其选定方式有很多种。 根据流体各向同性假设，所寻求的关系应与坐标选择无关，而由张量理论 可知，二阶张量中主对角线上三个分量之和是不变量，即不随坐标旋转而 变化，所以假设：