《2022年环境水力学》由会员分享,可在线阅读,更多相关《2022年环境水力学(10页珍藏版)》请在金锄头文库上搜索。
1、学习必备欢迎下载第一章液体流动的基本概念和基本方程(4 学时)1.1 基本概念:一研究对象:连续介质假定,使物理量为空间坐标和时间的函数。描述流体运动特性的物理量v,p,T,C。基本特征参量。 lagrange Method(拉格朗日)Euler Method (欧拉)两种方法研究对象不同:流体质点空间点流体微团微团控制体流体系统控制体二.基本参量表示法:用两种方法表示的基本参量方法不同。Lagrange 法:标量( p,T,C)p=(a,b,c,t)质点迹线= (a,b, c,t)矢量( , , , )da b c tdttdadtEuler 法:( , , )x y z tx,y,z 变,
2、()dadtt附体性dxudtdyudtdzudtijkxyzuiv jwk所以iiiijuuautu(张量形式)三.迹线和流线。迹线:dxudtdyudtdzudt流线:0ddxdydzuvw(恒定流时重合)四.质点导数。液体质点的流动参数B 随时间的变化律的欧拉法表示。也称为随体导数。DBBDtt迁移变率当地变率算子精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 1 页,共 10 页学习必备欢迎下载恒定流:B0t均匀流:()0B不可压:0DDt五.任意度量中系统体积分的随体导数。00BdB 为0系统体内积分。例:00dm0DmDt(连续性方程积
3、分形式)(一般将其变为欧拉法形式)00()dM动量DMFDt(动量方程)2()2e020()2edEDEWDt(外力所做功)e为内能(随温度、压力变化的能量)单位质量流体所具有的内能,状态函数。输运方程(transport Equation) 、 (L-E)00()EularBDdddDtt:0在0t时刻所占领的控制体。物理意义:00DdDt系统体积分的随体导数。00dt控制体内物理量体积分( B)的当地随时变化率。()AddAn(高斯定理)An dA从封闭面A 流出的的体积分, 也就是系统中一个位置移动到另一个位置,由于流场的不均匀性而改变引起的的体积分的迁移变化率。1.2 运动液体的应力和
4、应变关系本构方程一.流体微团运动的分析。1.微团运动 =平移 +变形(线变形、角变形)+转动精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 2 页,共 10 页学习必备欢迎下载若 A 点流速为:Auvw,则距其距离为(,)d r dx dy dz处点流速可表示为:xxxyxzyxyyyzAzxzyzzSSSuudxvvSSSdyw drwwdzSSS转动引起的流速增量平移变形引起的流速增量其中xyzijkw drwwwdxdydz2.变形率张量:12jiijjiuuSxx1 2 3i、,1 2 3j、ij线变形,ij角变形变形率张量具有对称性:ij
5、jiSS(6 个独立)ijS称为体积膨胀率,ijxxyyzzuvwSSSSvxyz(速度的散度)角转速分量,定义为12rotv(12速度的旋度)1122xyzijkvijkxyzuvw若0,无旋、有势流动,速度有势v0lvdl(环流量 =0)存在势函数,称为流速势满足20(拉普拉斯方程)区别,力有势f,力势函数,0lFdl 二、运动流体中的应力静止:000000ppp,只有压应力(负号表示与作用面外法线方向相反)精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 3 页,共 10 页学习必备欢迎下载运动、理想液体(不存在粘性切应力)000000xxyy
6、zzppp,xxyyzzmpppp粘性实际运动:xxxyxzyxyyyzzxzyzzppp,zx:x 表示作用面法线方向,z 表示力方向。由于存在切应力,所以,法向应力xxyyzzppp,但xxyyzzpppconst(不因坐标变化而变),所以,引入动水压强:13mxxyyzzpppp与作用面无关,各方上压力应力被认为是mp加上一个附加压应力,如xxmxxppp。切应力具有对称性:yzzyA:对于层流000000mxxxyxzmyxyyyzLmzxzyzzpppppp其中xxxyxzyxyyyzzxzyzzppp粘性附加应力张量B:对于紊流222xxyxzyxzyztLzxzyzuu uu u
7、u uuu uu uu uu其中222xxyxzyxzyzzxzyzuu uu uu uuu uu uu uu为紊流附加切应力项(混渗应力or Re 应力)简写为ijijijtL三、应力应变的关系(本构方程)1、牛顿流体ni ji jsn=1。本构关系符合内摩擦定律,应力和应变线性关系。非牛顿流体1n2、本构方程:各向同性、牛顿流体、层流,应力应变关系:牛顿内摩擦定律,Stokes(斯托克斯)推广。精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 4 页,共 10 页学习必备欢迎下载正应力:2iijmjuppxij(不可压)223iijmjuppvx
8、ij(可压)切压力:jiijjijiuuxxij1.3 连续性方程质量守恒0DMDt(对理想、实际流体都适用)一、积分形式:0dv ndAtv ndAv d:控制体一维恒定:1 112 22v Av A二、微分形式0vtv:流出控制体的质量(单位时间、单位体积)0DvDt不可压0DDt0v(恒定、非恒定流均可)恒定0DDt0v(不限制压缩否)1.4 运动方程一、微分形式:Fma1、ideal:Eular 方程:1dvvfpvvdtt2、实际牛顿流体、各项同性、不可压缩、层流N-S 方程21vvvFpvtvt:当地变化率,vv:迁移变化率(非线性项),vvvt:惯性项。F:质量力项,1p:压强梯
9、度项,2v:粘性项。3、对紊流,将上述方程取时均,得到(雷诺方程)精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 5 页,共 10 页学习必备欢迎下载一、上述方程可以简化为:I II III A:恒定I0B:绝热 II=I0C:理想流体I I InAAp vdApv ndA(压力)D:只有重力作用F则 IIAv dzv n dAE:不可压缩、绝热、理想、恒定、内能不变Aev n dA则推出202Apvv n dA一维问题的伯努利方程,(22pv为机械能)。二、微分形式:化面积分为体积分:nAAp vdAn vdAv d引入 Forier 定律:Tq
10、Tnn-热传导系数AAq dATndAT d将0去掉积分量得(1.40)22RDveF vvqTDt微分形式在直角坐标系下、不可压缩,上述微分形式可简化如下:公式(1.43)(1 6)(79)RDeqDt项项DeDt-内能变化,(16)项应力做功(粘性耗散) ,(79)项热传导,Rq辐射。对大多数液体:vec Tvc- 定容比热引入热扩散率(导温系数)pcpc- 定压比热一般取vpccc1-6=项- 耗散函数则方程为:2RqDTTDtcc去掉(忽略)、Rq影响,即得热传导方程2DTTDt(有关方程柱坐标下表示,参见书,不详述) 。1.6 基本方程组的封闭问题精选学习资料 - - - - - -
11、 - - - 名师归纳总结 - - - - - - -第 6 页,共 10 页学习必备欢迎下载对于不可压缩、各项同性、牛顿流体、层流。方程组有5 个,未知数v,p,T(6 个未知数)补充状态方程pRT。对于不可压缩液体const,所以,直接用连续运动方程求p,v 后求温度场。1.71.81.9 有势流动有涡流动边界层概念1、有势流动,实际流动满足何条件可简化为有势流动求解。0或者AudlndA单连域-涡通量,(2v)根据 kelvin s(开尔文)定理:在有势质量力作用下,理想流体,正压条件下,(流体质点组成的封闭曲线的速度环量不随时间变化)即0t。对于高 Re 流动粘性项作用可忽略,认为理想
12、流体, 由静止状态开始的流动,如重力坝泄流,圆柱绕流,闸孔出流。波浪等,可近似认为是势流。2、势流特征:存在势函数, , ,x y z t,v图不可压势流20基本方程20,22pvzcg(Eular 方程积分)理、恒、不可压、有势(力)势流。单独变量, p。平面、不可压势流:存在流函数( , , )x y z,量纲2ms图1满足20求解方法:势流叠加法、复变函数法、流网法、实验法。3、涡量满足涡量方程(推求不详细讲)。4、边界层概念:如上所述,对高Re,N-S 忽略得到 Euler(物面为流线) ,但用此法求解平板绕流时,阻力=0,不符合实际,物面应给粘滑条件。引入边界层概念:在固体壁面附近的
13、一个薄层中,流速梯度变大, 粘性作用必须考虑称为边界层。流动分区:边界层内,粘性流。边界层外,理想流。分界:引入边界层厚度概念(又分为排挤厚度或流量损失厚度1动量损失厚度2能量损失厚度3) 。层外势流解1构成的线为流线,解外层势流时边界条件。精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 7 页,共 10 页学习必备欢迎下载恒定二维、不可压缩、质量力不计(曲率小),假定L,yx,vu量纲比较,忽略量,简化为边界层方程(1.97) 、 (1.801.81) 。由边界层方程知0py即,边界层内法向压强基本不变,即可用边界层边界上势流解的p,所以边界层方
14、程中p 已知,只有u、v 两个变量,一般数值解。引入边界层概念目的:解决沿程水头损失计算绕流阻力计算,紊动射流(边界层流动)边界层动量积分方程(1.85)未知数u、0三个未知数,一般补充ufy,00ydudy(层流边界层,紊流边界层不同) 。1.10 小密度差流动,Bousinesq s 近似一、 Bousinesq s 近似:密度变化的影响仅仅在运动方程的重力项中考虑对其他各项影响忽略不计。二、小密度差流动的基本方程:许多环境问题, 由于水中加入污染物质,使密度发生变化,但这种密度差很小,严格按考虑密度变化的可压缩流体分析比较困难,常常采用Bousinesq s 假定、作近似处理。对于质量力
15、只有重力、不可压流体、其运动方程为:2DvpgvDt若,g以静止时值为参考状态,即0,0ppp代入上式,并知静止时00g所以,2000011DvpgvDt由于01所以,惯性项中可忽略,令*0gg(折减重力加速度)*2001DvpgvDt连续性方程不变0v以上为基本方程1.11 旋转流体流动,科里奥利力效应(Coriolis effect )一、对于较大的江河、湖泊、海湾等水域的流动问题,将地球看成惯性坐标系统会引起一定的误差,此时要考虑地球转动对水域流动的影响。根据理论力学知:精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 8 页,共 10 页学习
16、必备欢迎下载绝对运动(定坐标系)=相对运动(动坐标系)+牵连运动(动坐标系对定坐标系)对于动坐标系为定轴转动的情况:IRekaaaaIa绝对加速度,Ra相对加速度,ea牵连加速度,kaCoriolis 加速度。earkazv相对运动和牵连运动相互作用结果。由于牛顿第二定律对于定坐标系进行的,所以,对于流体质量力为重力,N-S 方程为:21IIDvuapguDt所以,对于旋转动坐标系的相对运动微分方程为:21IekRDvuaaapgurzvDtr可表示为2212r与p合并考虑。二科化力效应判别数。有两个:(Rossby)罗斯比数0RuL惯性力科化力(Ekman)埃克曼数2kEL粘性力科化力由公式
17、可知:0R减小,kE减小,科化力大,不能忽略。对环境问题:地球自转不变, L 大时,0R减小,kE减小。所以,对大面积水体流动问题,如海湾、湖泊、大气环流、科化力要计入。第三章稳流模型(一)方程的封闭问题及稳流模型的概念。1.瞬时不可压缩流动N-S 方程:连续:0iiux运动:21iiijijijjuuuuFtxxxx方程组 4 个,未知数4 个,理论上可直接求出解,也就是说明N-S 方程包含了所有的湍动信息,可直接求解,也就是进行直接模拟(DNS ) 。直接法要求:t,iu都很小与脉动量同一量级。计算量很大。精选学习资料 - - - - - - - - - 名师归纳总结 - - - - -
18、- -第 9 页,共 10 页学习必备欢迎下载边界条件给定难度大。目前只有对Re 小流动,简单边界条件进行计算的算例无法用于工程实际。对工程问题, 往往不要求物理量的瞬时值,只要知道时均值随时间变化就可以了,因此转向解时均流动。2.时均不可压缩紊流(Reynolds 方程)011iiiiijijijijjuxuuuuu uFtxxxx方程 4 个,未知数,iu, 6个独立的ijuu( i,j=1,2,3)总共 10 个未知数。所以方程不封闭。3.紊流模型:为了使上述方程封闭,根据紊流特性需附加的条件。精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 10 页,共 10 页