07粘性流体动力学基础.docx
- 文档编号:10918170
- 上传时间:2023-02-23
- 格式:DOCX
- 页数:26
- 大小:267.89KB
07粘性流体动力学基础.docx
《07粘性流体动力学基础.docx》由会员分享,可在线阅读,更多相关《07粘性流体动力学基础.docx(26页珍藏版)》请在冰豆网上搜索。
07粘性流体动力学基础
第七章粘性流体动力学基础
第一节粘性流体运动的基本方程
采用流体力学微元体平衡分析方法可以推导出粘性流体运动的基本方程组,该方法可参考本书的第二章和第三章。
本节将直接由两大守恒定律(质量守恒定律和动量守恒定律)来建立控制流体运动的基本方程组。
首先需要给出空间某点物理量的随体时间导数表达式、雷诺输运方程以及本构关系。
一、随体导数
描述流体运动规律有拉格朗日和欧拉两种基本方法。
拉格朗日法着眼于确定的流体质点,
观察它的位置随时间的变化规律。
欧拉法着眼于从空间坐标去研究流体流动,它的描述对象是流场。
随体导数的物理意义是:
将流体质点物理量q的拉格朗日变化率以欧拉导数的形式
表示出来。
随体时间导数的数学表达式为:
四V'、q(7-1)
dtft
式中右边第一项代表由时间的变化所引起的变化率,也就是由于场的时间不定性所造成的变化率,叫做当地导数。
第二项代表假定时间不变时,流体质点在流场中的位置变化所引起的变化率。
这是由于场的不均匀性造成的,叫做迁移导数。
二、雷诺输运方程
雷诺输运方程描述了积分形式的拉格朗日法和欧拉法的时间导数的变换关系。
设封闭系
统在t时刻占有体积Q(t),如图7-1所示。
其中关于物理量q的总量的随体时间导数有
图7-1封闭系统输运示意图
d;:
a
qd半崇」qVndS(7-2)
dL】t tst 其中s(t沪封闭体积的曲面,n为曲面的法向向量。 上式表明: 封闭系统中,某物理量总和 的随体导数等于该瞬间与该系统重合的控制域中该物理量总和的当地时间导数(非定常效应) 和通过控制面流出的该物理量的流量(对流效应)之和,此即为流体的雷诺输运方程。 用广义的高斯公式将面积分转换成体积分,上式也可以写成 §qd”=彳"qV心(7-3) dt“t: i二FM 三、连续方程 连续性方程反映了流体在运动过程中必须满足质量守恒定律。 其中拉格朗日法的研究对象是流体中一个确定质量的流体物质团(称为封闭系统),随着流体的运动,封闭系统的表面的位置会不断随时间而变化,但没有流体穿过它的边界。 质量守恒定律可表述为: 封闭系统内流体的质量在流体运动的过程中不发生变化。 而欧拉法的研究对象则是流场空间中一个固定的区域(称为控制域),控制域表面的位置不随时间而变化,由于流体的运动,控制域的表面通常会有流体通过。 质量守恒定律可表述为: 控制域内流体质量随时间的增加与流体经控制体表面流入的质量相等。 在式(7-3)中令q=P,可得连续方程 (7-4) 皿一+^何)dQ=。 考虑到积分体积的任意性并假定被积函数连续,上式可以写成 (7-5) 这是基于欧拉观点的微分形式的连续方程。 它表明控制体中流体质量在单位时间内的增加来自流体质量经控制体表面的流入速率。 将随体时间导数表达式代入上式,便得到基于拉格朗日观点的微分形式的连续方程。 (7-6) 对于不可压缩流动,恒有成立,此时连续方程简化为 (7-7) 连续方程仅反映了流体的运动学特性,与流体的本构关系无关。 动量方程反映了流体的动力学特性,因此需要先介绍本构方程。 四、本构方程 本构方程反应了应力和应变率之间存在的制约关系,这是建立流体动力学方程的基础。 真实流体的力学性质是很复杂的,不同种类的流体可能表现出完全不同的力学特性,即便是 同一种流体在不同的外部条件下,比如温度不同时,力学特性也会有很大的差异。 因此要建 立一个普适的本构方程几乎是不可能的。 Stokes提出了适用于牛顿流体的如下三条假设: (1)流体是各向同性的,也就是说流体的物理性质与方向无关,只是坐标位置的函数; (2)应力张量bj是应变率张量罚的线性函数,与旋度无关。 (3)静止流体中,切应力为零,正应力的值为流体的静压。 根据以上假设,考虑到应力张量和应变率张量的对称性,由张量理论便可以推导出应力 和应变率间的关系如下: Jj=-P、'ij・,ekk、・j2忸(7-8) 其中 (7-10) 称为体积粘性系数。 这表明热力学平衡压强或静水压强p与力学压强相差rBekk。 式(7-8) 也可写成: 1 气=(_p+%ekk勤+2Peij—gekkS|(7-11) 对于单原子气体=,PB=0。 对于多原子牛顿流体,根据Stokes假设,通常满足体积粘 性系数Pb为零的条件,不必区分力学压强与热力学压强,本构方程简化为 sj=-p§j+2卜ej-―ekk^ij|(7-12) <3) 其中只含有动力粘性系数,该本构关系样适用于静止流体、理想流体()。 五、动量方程 动量方程在物理上反映了流体在流动过程中满足的动量守恒定律。 基于拉格朗日观点,动量守恒定律可叙述为: 封闭系统内流体动量随时间的变化率等于作用在该系统上所有外力之和。 其数学表达式可以写成 d一— VdC;.: iiEfd",iinES(7-13) dt"•”S 在雷诺输运方程中(7-3)式中,令q=pV并代入上式,可得到基于欧拉观点的积分形 式的动量方程 --N'、: VVd,)=..f"nES(7-14) 利用广义高斯公式将上式中的面积分项改写成体积分,考虑到积分体积的任意性并假定被积函数连续,则有 (7-15) 这是基于欧拉观点的微分形式的动量方程。 以连续方程(7-5)代入上式,得到动量方程的另 一种常见的形式 (7-16) 将牛顿流体的本构方程式(7-12)代入式(7-16)后,得到牛顿流体的动量方程(或称为 Navier-Stokes方程) (7-17) 式中V2即为Laplace算子,学为运动粘滞系数。 在不可压缩流动中,有 (7-18) 对于理想流体的假设,则可简化简化为欧拉方程 (7-19) 第二节边界层的概念 由于方程的非线性特征,使得问题的求解非常困难。 在许多情况下,需要根据流 动的特点对方程进行不同程度的简化。 在低雷诺数流动中,由于粘性力远大于惯性力的特点, Stokes近似将方程的惯性力项略去,使基本方程得以线化,得到了具有一定精度的小 球阻力公式。 在Oseen近似中,在方程中保留了线化的惯性力项,使小球绕流的远场特性得到了改善。 大雷诺数流动的情况相反,惯性力项远大于粘性力项。 作为近似将粘性力项略去后, 方程化为无粘流体的欧拉方程。 若使用与它相匹配的无粘流的可滑移边界条件,对固 体的绕流问题会出现零阻力的非物理解(达朗贝尔佯谬);若使用无滑移的粘性固壁条件会导 致数学模型在边界条件上的过约束。 为了解决大雷诺数情况下欧拉方程和粘性边界条件间的矛盾,普朗特(1904)引入了边 界层的概念。 对绕流问题,他认为在固壁附近的很薄的一层区域内,沿固壁切向的速度由外部势流的值迅速下降为零,以满足粘性流体的固壁边界条件。 如图7-3所示,边界层形成的 原因也可通过从涡旋传输的观点来解释。 流动中的任何固体边界层都相当于连续分布的涡源, 它不断的在流动中产生涡旋。 紧靠表面附近的涡旋,一方面向外扩散,另一方面随着流体向下游流动。 涡旋扩散的速度取决于流体的粘性系数,粘度越大,扩散得越快,而涡旋向下游流动的速度取决于来流速度。 当雷诺数足够大时,平板表面附近的涡旋向下游流动的速度比向垂直于流动方向的速度大得多,以致包含这些涡旋的流动仅仅限于贴近表面的一个向下游伸展的薄层,这个薄层就是边界层。 在边界层内,流动是有旋的;而边界层以外的流动则可视为无旋的。 目前边界层理论已成为近代流体力学的重要基石,它澄清了大雷诺数流动问题中粘性对流动的影响。 在许多情况下,大雷诺数与湍流相互关联,本章将分节讨论低速层流边界层和湍流边界层。 边界层理论基于大雷诺数流动的近似,首先需在近似中保留部分粘性项而建立Prandtl边界层方程。 为了说明边界层的基本特征,本章将先引出描述边界层的数学方程式,接着讨论一个最典型的边界层流动(平板边界层),然后再介绍边界层分离现象。 二「昼--一—三 二门成3膏9」,」*nA— U————一.卫 图7-3边界层内的涡旋和速度分布示意图 第三节边界层的微分方程式 由粘性流体力学的基本方程,采用量级分析方法和普朗特展开方法都可以推导出边界层 的微分方程式。 本节将介绍第一种方法。 考虑大雷诺数的二维绕流问题,假定固壁是平直的(平板或楔)。 设轴与壁面垂直,轴与壁面平行且指向下游,坐标原点和顶点重合,如图7-3所示。 连续性方程和动量方程的两个投影分别为 代入式(7-20)中并将各项的量级标注如下: .: t *Vx .: VxVL x——V U : x ■: Vy Rect Vx ■: Vy -Vy 当Re>>1时,在式( y : x Re -y UcVx _+V 印y * Lex a * 1 1 Po0 * 叩 / 1 + PU2 * ex Re 1 Re Peg _PU2 * £P_ * 剑 工1 *Re2 * * * * (i ■2X/: -Vx 一*2 : x (-2 3Vy 一*2 ;x Re2 =0 2*2 、•2-y Re) Re -2: Vy *2 (7-21a) (7-21b) (1 Re) (7-21c) 7-21)中略去高阶小量,并恢复为有量纲的形式可得 -: Vx3y XyM L'\ x: y (7-22a) .: Vx: : Vx: : vx1;: p ——Vx——Vv tx-y -t: x (7-22b) 式(7-22c)表明边界层内压力的法向梯度近似为零, 只是 分布来描述。 考虑到边界层外的势流区的速度只有 简化为 £Ve ft 将(7-22b)式中的压力项用外流速度表示,得到 -: vx3x3x —-Vx—-Vv— ■xy txy (7-22c) 的函数,它可由外部势流区的压力 个方向的分量 dve1dp Ve一=—二一 dx|dx vx=ve(x),则欧拉方程 (7-23) (7-24) 对于曲率不大的二维曲面壁而言,分析表明,只要将x取成沿壁面的流线坐标,边界层 方程的形式与上式完全相同。 在边界层方程中,保留了惯性力项,部分保留了粘性力项,压力项由外部势流解给出定。 与Navier-Stokes方程组相比,边界层方程组是大大的简化了,方 程由椭圆型方程变为了抛物型方程,使问题的求解由二维无穷域变为一个半无限的长条域。 对于前者必须在封闭边界上给出边界条件,而对于后者则下游边界条件无需给定,只需给出: 边界层方程仍然是非线性的。 边界层内的解与外部势流区的解在边界层的边缘上衔接,在给定边界层方程外部边界条件后,对边界层方程的求解时,则需要对边界层厚度的定义加以说明。 第四节边界层厚度 边界层是在大雷诺数流动中近壁处的涡量集中区。 由于全流场中从粘性区向无粘区的过渡是逐渐进行的,不存在一个非此即彼的明确界限,因此边界层的边缘并不非常清晰。 为了实际应用的方便,边界层厚度有着如下几种较为严格的定义。 即边界层的位移厚度&,边界层的动量损失厚度&,以及边界层的能量损失厚度63, 上节中提及的边界层厚度6即为边界层与外部势流的边界,亦称为名义厚度。 一、边界层的位移厚度-1 由于壁面摩擦的影响,与理想流体相比,边界层内实际流过的体积流量会有所减少。 为了使基于理想流体理论计算得到的流量与粘性流的实际情况一致,需要把原来的固壁向外推一个距离51,该距离被称为边界层的位移厚度。 如图7-4所示,矩形OACE勺面积与相当与减 y>: ■-一 少流量的面积OD函相等,对于不可压流动Pve81=fP(ve—vxdy,即: ■0 (7-26) vx 1dy 在实际问题中,往往应该考虑边界层的存在对外部势流场的影响。 例如溢流坝面流动中, 下泄流量不变,但随着边界层发展,必然迫使自由水面抬高一个位移厚度。 又例如,对于低速风洞的试验段,不能设计成一个平直段,通常有一个约0.5°的扩散角,以补偿边界层增厚 的影响。 8。 在定常流 式(7-26)的积分上限为无穷,在实际计算中,通常取为边界层名义厚度 度总小于边界层厚度。 二、边界层的动量损失厚度2 边界层厚度的另一种定义是基于考虑边界层存在所导致的无粘流中流体动量的损失。 在 通常将52称为动量损失厚度或简称动量厚度。 三、动能损失厚度3 与动量损失类似,边界层由于粘性阻滞而造成的动能流通量损失为: dy。 这一动能通量损失,如果用物体表面的一层无粘性流体流动的动能通 上式就是动能损失厚度53的定义。 从理论上说,粘性的影响可以由壁面一直延伸至无穷 远处。 本节在动量厚度和能量厚度的表达式中,已将积分上限由*换成8;前者基于渐近边 界层概念,后者基于有限边界层概念。 第五节边界层动虽积分关系式 边界层微分方程的求解通常可分为近似解和相似解(精确解)。 在工程计算中,求解边界层问题可以使用各种近似解法,以期较为迅速地得到具有一定精度的计算结果。 这种方法是针对边界层微分方程积分后得到的动量方程进行求解的,故本节先介绍卡门动量积分关系式,下节将介绍它在平板问题中的运用。 首先推导动量积分关系式,从二维层流边界层的微分方程出发,即: (7-29a) -Vx: Vx: Vx: VedVe1: /、 +V^—+Vy=+Ver+;丁(7-29b) t: x: y: tdx-: y 式中t=卜史兰。 以Ve乘以式(7-29a)两边,考虑到Ve=Ve(x,t),得到: y (7-30a) -: VeVx: 'VeVy;: Ve 一、"一、=Vx一x: y: x Vx乘以式(7-29a)再加上式(7-29b)而得到: (7-30b) 式(7-30a)与式(7-30b)逐项相减后: -一-Ve: 1C,、 了(Ve—Vx)+=(VeVx—VxVx)+(Ve—Vx)=+=(VeVy—VxVy)=—;=(7-31a)-t: x: x: y|: y 固壁边界条件如下: y=°°: Vx=Ve,7=0 利用边界条件式(7-31b)对式(7-31a)沿边界层厚度方向积分;由于式(7-31a)左边第4项 积分后自动为零,故得到: 二〃w 0Ve-Vxdy=一 .、.: Ve Ve-Vxdy一IVxVe-Vxdy : X0fx 根据边界层位移厚度和动量损失厚度定义式,上式可以写成: 《三(Ve&1)+^2+仲1+2^2)1冬=^2=色(7-32) Ve一t;xVe一x: ? Ve2 A&.... 令H=—为边界层的形状因子,对于层流边界层其变化范围大约为H财2.0~3.5,对于 2 湍流边界层其变化范围大约为H充1.3~2.5。 可以得出定常不可压层流边界层的积分关系 式: d、2-,,2dVeCf ——+(2+H)=——(7-33) dxvedx2 该式即为卡门动量积分关系。 在动量积分关系中,未知量是壁面摩擦系数Cf、动量厚度^2和形状参数H(或位移厚 度每),它们都可由速度分布Vx确定。 作为近似解,如果能给出一个恰当的取决于单参数的 速度分布Vx,边界层问题便归结为求解单参数的卡门动量积分关系。 以下平板边界层的例子来说明卡门动量积分关系的应用。 第六节平板层流边界层的计算 d、2 卡门动量积分关系式简化成 (7-34) d(Pv32)=iwdx表明壁面阻力所做的功等于边界层动量损失的增量。 动量积分法求解,需 要先假设一个满足边界条件的速度分布,但并不要求与边界层内每点的流速都符合。 设工=f0),其中n=y,边界条件(7-3ib)需要作进一步的分析。 为了使近似的速度分布 Ve 与实际情况相近,速度在物面和边界层外缘必须满足边界条件 为使边界层外缘的速度过渡平滑,速度的一阶导数满足边界条件(零切应力) 速度二阶导数满足的边界条件是 其中上面左式由边界层方程式(7-24)使用壁面边界条件得到的,反映了壁面速度剖面曲率 和压力梯度间的关系,称为“壁面第一兼容条件”;右式代表边界层外缘速度的过渡平滑。 速度的三阶导数满足 =0, yz0 其中第一式由边界层方程对 y求导后在物面取值得到,第二式仍代表边界层外缘速度的过渡 平滑。 对平板边界层即有: f(0)=f''(0)=f"'(0)=0,f (1)=1和f' (1)=f" (1)=f"' (1)=0。 如果必要可以使用速度更高阶导数的边界条件。 (7-35) \17 3 7- ndf - 1f 1! o - 2 a 中 其 1 、2=、f1-fdy=a2、 式(7-34)可写为: d、.f0 v=: dx上2ve 可以认为平板边界层是从前缘开始的,当x=0时&=0,积分上式得: (7-38) 、x=.2f'。 X ■,: •2\Ve 由上式可以得到边界层的各种参数,例如: 壁面摩擦阻力系数: 平均壁面摩擦阻力系数: (7-39) 动量损失厚度: (7-41) 位移厚度: (7-42) 1 其中%=0(1-fdn f=a0a1a22a33 一2 3a〔 =0.646 Ve 、x =1.740 表7-1所列为选用各种函数的计算结果比较。 最简单的情况是线性分布fg)=n,边界条 件自动满足,由于没有待定常数,关于速度导数的边界条件不必考虑。 读者可进一步用其它形式的函数进行计算,并将它们与Bulasius的精确解(相似解)进行比较。 表7-1平板边界层近似解 速度剖面 Ve «1 «2 1一 f(0: %产 Vvx Vvx 我: VexPve2VV ML CdJe \V f(nZ 1 2 1 6 1 1.732 0.577 0.289 1.155 ")=2—芋 1 3 2_ 15 2 1.825 0.730 0.365 1.460 fL)=着_Ln3 22 3 8 39 280 3 2 1.740 0.646 0.323 1.292 f「)=sinW <2) Jl—2 31 4* 2兀 兀 2 1.741 0.655 0.327 1.310 Blasius解 1.721 0.664 0.332 1.328 第七节平板湍流边界层 实验证明,当雷诺数ReL=veL/v较小时,平板上的边界层流动可能全部都是层流边界 层。 当ReL足够大时,就可能同时存在三种流态: 在靠近平板的前缘一段是层流,在后缘一段是湍流边界层,在此两段之间是过渡段。 通常边界层的概念出现在大雷诺数的流动中,但是大雷诺数是指外流而言的。 在边界层 内的流动,特征尺度是边界层的厚度 当地雷诺数应表示成ReB=ve6/v。 由于边界层的 厚度远小于全流场的特征尺度,其中的扰动并不能马上发展起来。 若在平板前缘施加一个随 机的小扰动,当雷诺数小于某一临界值时,随机脉动在向下游的过程中会很快消失,边 界层保持稳定的层流状态。 当一旦达到某一临界雷诺数的值(根据稳定性理论来分析) 便进入向湍流转振的过渡区域。 在向下游继续发展的过程中,最后将使边界层迅速进入完全的湍流状态。 概括起来,湍流边界层的有如下的分层结构: 粘性底层(0主y<5Y) v* : : : y•壬30—) v*v* 海迹律层(0.28& 外层 [粘性顶层(0.485 (7-43) 其中 v・="w 具有速度的量纲,与壁面切应力有关,称为摩擦力速度(或剪切流速)。 湍流边界层的速度分布可使用Prandtl混合长度理论来分析。 考虑无限长平板在自身平面 内运动引起的二维定常平行剪切流,湍流的平均速度只有沿轴方向的一个分量,且只与坐 标有关,这时可得到 ,粘性力项与雷诺应力项比较 考虑到湍流的大雷诺数特性,在壁面之外的大部分区域内 是个高阶小量,上式可进一步近似为 Pl”dvx/dyf=Tw,这表明 (7-45) 称为平板湍流的速度对数分布律。 由于该处雷诺应力变小,粘性力增大, 以上分析在非常接近壁面处是不成立的, 上式的对数项在该处已经发散。 在非常接近壁面的区域,考虑到粘性力占主导地位的事实,有 dVx •-=-w dy 积分得到满足边界条件的解为 2 V.. Vx==—yV 上式表明流体的速度呈线性分布,由于该层中粘性力占主导,通常称为粘性底层。 粘性底层 比边界层薄得多,其外满足速度对数分布律的区域称为对数律层,两者间存在一个过渡层。 而工程上通常采用如下的近似方法求解平板湍流边界层。 在没有压力梯度的情况下, 事实上, (7-47) (7-48) 在105 当Rex变小的话,n将比7小,随着Rex的增大,n会变 大,在有压力梯度的情况下,指数定律的精度会降低。 采用上述指数函数形式求解平板湍流边界层,并借鉴圆管湍流关于壁面阻力的研究成果,假设 17 与平板层流边界层的计算大似,有&2=8gf(1—fdy=,2&=—6,利用动量积分式(7-34) 有 -/-1/4 7d5v =0.0225、(7-50)72dx) 可以认为平板边界层是从前缘开始的,当x=0时6=0,积分上式得: 6^>0.3^^^(7-51) Re; vex,,„,,工一,一、,,,,,, 其中Rex=^,上式即为湍流平板边界层厚度沿x方向的计算公式,与式(7-38)相比,湍 v 流边界层厚度是x的0.8次方,而层流边界层厚度是x的0.5次方。 可见湍流边界层厚度增 长比层流边界层要迅速。 与第六节类似,可求解如下物理量 边界层动量厚度: 壁面摩擦阻力系数: 平均壁面摩擦阻力系数: CDf=-[Cfdxn.072’^i(7-54) l0 式中的系数可根据实验数据稍加修正,得出Prandtl湍流边界层阻力计算公式: CDf=0.074’担),适用范围为105 相应的Blasius层流边界层阻力计算)* 公式见第六节,图7-5给出它们随雷诺数变化的曲线表示形式(曲线①和②所示)。 针对实 验数据的其它拟合曲线形式本节不作详细介绍。 s 为00001 图7-5平板阻
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 07 粘性 流体动力学 基础