07粘性流体动力学基础.docxVIP

第七章粘性流体动力学基础 第一节粘性流体运动的基本方程 采用流体力学微元体平衡分析方法可以推导出粘性流体运动的基本方程组，该方法可参 考本书的第二章和第三章。本节将直接由两大守恒定律(质量守恒定律和动量守恒定律)来 建立控制流体运动的基本方程组。首先需要给出空间某点物理量的随体时间导数表达式、雷 诺输运方程以及本构关系。 一、随体导数 描述流体运动规律有拉格朗日和欧拉两种基本方法。 拉格朗日法着眼于确定的流体质点, 观察它的位置随时间的变化规律。欧拉法着眼于从空间坐标去研究流体流动，它的描述对象 是流场。随体导数的物理意义是：将流体质点物理量 q的拉格朗日变化率以欧拉导数的形式 表示出来。随体时间导数的数学表达式为： 四V、q (7-1) dt ft 式中右边第一项代表由时间的变化所引起的变化率，也就是由于场的时间不定性所造成 的变化率，叫做当地导数。第二项代表假定时间不变时，流体质点在流场中的位置变化所引 起的变化率。这是由于场的不均匀性造成的，叫做迁移导数。 二、雷诺输运方程 雷诺输运方程描述了积分形式的拉格朗日法和欧拉法的时间导数的变换关系。设封闭系 统在t时刻占有体积 Q (t )，如图7-1所示。其中关于物理量 q的总量的随体时间导数有 图7-1封闭系统输运示意图 d ；:a qd半崇」qV ndS （7-2） dL】t it :t st 其中s（t沪封闭体积的曲面， n为曲面的法向向量。上式表明：封闭系统中，某物理量总和 的随体导数等于该瞬间与该系统重合的控制域中该物理量总和的当地时间导数 （非定常效应） 和通过控制面流出的该物理量的流量（对流效应）之和，此即为流体的雷诺输运方程。用广 义的高斯公式将面积分转换成体积分，上式也可以写成 § qd” =彳qV 心 （7-3） dt“t：i 二F M 三、连续方程 连续性方程反映了流体在运动过程中必须满足质量守恒定律。其中拉格朗日法的研究对 象是流体中一个确定质量的流体物质团（称为封闭系统），随着流体的运动，封闭系统的表 面的位置会不断随时间而变化，但没有流体穿过它的边界。质量守恒定律可表述为：封闭系 统内流体的质量在流体运动的过程中不发生变化。而欧拉法的研究对象则是流场空间中一个 固定的区域（称为控制域），控制域表面的位置不随时间而变化，由于流体的运动，控制域 的表面通常会有流体通过。质量守恒定律可表述为：控制域内流体质量随时间的增加与流体 经控制体表面流入的质量相等。 在式（7-3）中令q = P,可得连续方程 (7-4)皿 一+^ 何）dQ =。 (7-4) 考虑到积分体积的任意性并假定被积函数连续，上式可以写成 （7-5） 这是基于欧拉观点的微分形式的连续方程。它表明控制体中流体质量在单位时间内的增加来 自流体质量经控制体表面的流入速率。将随体时间导数表达式代入上式，便得到基于拉格朗 日观点的微分形式的连续方程。 （7-6） 对于不可压缩流动，恒有 成立，此时连续方程简化为 （7-7） 连续方程仅反映了流体的运动学特性，与流体的本构关系无关。动量方程反映了流体的 动力学特性，因此需要先介绍本构方程。 四、本构方程 本构方程反应了应力和应变率之间存在的制约关系，这是建立流体动力学方程的基础。 真实流体的力学性质是很复杂的，不同种类的流体可能表现出完全不同的力学特性，即便是 同一种流体在不同的外部条件下，比如温度不同时，力学特性也会有很大的差异。因此要建 立一个普适的本构方程几乎是不可能的。 Stokes提出了适用于牛顿流体的如下三条假设： 流体是各向同性的，也就是说流体的物理性质与方向无关，只是坐标位置的函数； 应力张量 bj是应变率张量 罚的线性函数，与旋度无关。 静止流体中，切应力为零，正应力的值为流体的静压。 根据以上假设，考虑到应力张量和应变率张量的对称性，由张量理论便可以推导出应力 和应变率间的关系如下： Jj=-P、ij?，ekk、?j 2忸 (7-8) 其中 为动力粘性系数,丸为第二粘度。静压 是一个热力学状态参数 其中 为动力粘性系数, 丸为第二粘度。静压 是一个热力学状态参数 O在热力 学平衡态下，它总是等于三个相互垂直方向上正应力的平均值 的问题中，有相当一部分是接近平衡态的非平衡体系，这时 定的差别的。将上式的下标缩并后两边除以 -3后得到(力学压强)。在流体力学研究 学平衡态下，它总是等于三个相互垂直方向上正应力的平均值 的问题中，有相当一部分是接近平衡态的非平衡体系，这时 定的差别的。将上式的下标缩并后两边除以 -3后得到 (力学压强)。在流体力学研究 与一点处的平均压强 是有 P = - 3。kk = P - § 穴昼=P - 穴 Bekk (7-9) 其中 (7-10) 称为体积粘性系数。 这表明热力学平衡压强或静水压强 p与力学压