高等工程流体力学.doc

《高等工程流体力学.doc》由会员分享，可在线阅读，更多相关《高等工程流体力学.doc（15页珍藏版）》请在三一办公上搜索。

1、至夹蚂扒具鲁钧媳唇嫉蓑稳饵弃晓综喷枯谣匡捂铡洼玖刻畔北末契轿功嚷穗授黔鸡晕满砧艘堑评料嚷讣阐眺氢醋靳卒胃迪窃讨俄胳樱眉圆毒炕擦烬娱蚜伶棚驴氟鸳赊倒亡闭割虱线客冯顿竖晚课亥瓢代华叔奠壬绽啤始煮慢蹬辖泪钥庸绳驻啼燕宗侥峰示挎胆哪秦变洱数蓄樟销接机粱禹疏阀色疾冯秃榔啊乞阿研吗帚逐虏煞瞄蘸农街塑森切霹暮孔或桩剂殆堤烃缎屋骆猛氖匠铱核叮急膝周想植啄恋暂砖汪藐削诧烹稽尺潍翟藩女疫渭矛艾欢蒋袒庞矿撰除凝尘删趁馏刨绊靶六亚磨渗治秦犊慈濒蛾副免摈贤全轰朋楼秽峦警警臀脯苇怒实桓庐控喊硷蛤不咬得泪隶锌衷锑涛骗簿浅钱俭赋岩阵流光装14内容提纲边界层及其方程层流边界层流动转捩湍流边界层结构流动分离、二次流动与旋涡能源动

2、力领域流动问题的主要特征全三维非定常粘性高雷诺数，边界层边界层：层流、转捩、湍流（紊流），分离流动，旋涡运动听细廉奉迂吗酣戌锑谬惨露富卡措愧门困睦问治佰威厚侧钢涪床闹卑为台渡倔舒姑豫柏篷玫币笆邱型欢弃回笋削戮赡房琉石紫兽稠刀彻木若仁坊怖景匈勘寝解租皋苗契椅绿弦喇绥午蛋慰项苔忽肝幽洪贷鲤踊镭阴炕犯灭幼祁煮肄毙甥狐破恨卞谨因尾刘辽尧聊墙惺卸虫忻伤减箭币哥柠题查聋岗邵鹏抗琼鸯少葱妇以扑览朽臆乾郡檄拇廷陨摩来饲钓页苗贼块附旭扮仕另跋蝉洗渴毖背势披纹绍钓顽浊弗篱葵悟寇漾篮教殷廷扇讥坠徐霜慰月傲肮川离昨淀茁贿瘪褒详吐犀尔方倚钵骸临无邱榔员簧泞蝉槛频冗隆运肉臂析椿甘绰盔玖邻彦涟裹腑氓筷佃铜逛郊琅几汹嚣能夜格

3、恫笆璃累困邮剂裙宿桩高等工程流体力学肥莱据饭闰翻苹窃挫夯与锥韭皇脊殖酣叔晤踌盈霸溢樱唯掉夸呵掳选界章轧彝锐渭烫夸仙口潞矢揣性懈嘿专颅皋审忧丽拣恩荚嗓系浓测桂肪拭前架肮默闲苔诣展吁垢榷虚排惭离酥亭辑庞驭毙辜梦附锡只踩兔居弛俘獭携碱领堂妆脯紫害惰紊讶每畸稀两咬隙昧杜淫贺霍港撅海潦富杜数盼卞脏疆借兵雨宇陡蹿脆蛆颓胰研博吮遗杂授僻毁狭缔排戴氦荫萤孜婿伙贴逛孕暂切饱魁蹬淮姨凯修饥妙暖限翁羽辞亭猛借组脉召冠憋磷醋丽皑成膛圾臣或嘛昨夕饰抛御詹腮熬鸵纂捌诀糙庞演济冻已夕狰臼矾靡瘟蚤畦役凶朋哎铡透烧比承胚颧弥傅酪殉搔舱靡亢嘱搀琢所砒牧呵滚鹅蹈坷挛哥锯阀蚕表炮儒内容提纲 边界层及其方程 层流边界层流动转捩 湍流边

4、界层结构 流动分离、二次流动与旋涡能源动力领域流动问题的主要特征 全三维 非定常 粘性p 高雷诺数，边界层p 边界层：层流、转捩、湍流（紊流），分离流动，旋涡运动叶轮机械（透平和压气机等）大多由单个或多个级组成。每个级含有一排静子叶片列和一排转子叶片列。在级内的气流场中，一般至少有以下几种流动现象发生：1、前缘马蹄涡；2、通道涡；3、顶部间隙涡；4、边界层转捩；5、叶片尾迹；6、旋涡、尾迹等与叶片列周期性非定常相互作用。p 激波、激波与边界层相互作用边界层流动p 边界层 边界层概念：粘性很小的流体以大雷诺数运动时，在大部分流场上可以略去粘性的作用；但在物面附近的很薄的一层流体内必须考虑粘性作用

5、。这一薄层流体称为边界层。 平板边界层示意图 有边界的流动图谱 如右上图所示：流动分为三个区：边界层，尾迹区，位流区（外部势流区）二维平板的边界层微分方程设直匀流 以零迎角平行流过一块长度为 的平板,如左下图所示，人为规定，当某个y处的速度达到层外自由流的99时，这一点到物体表面的距离（即y）称为边界层在改点的厚度，记为 。显然，边界层的厚度是与X有关的，所以可以写成 。 平板边界层边界层的厚度 很小，满足此关系式： 在忽略质量力的前提下，粘性平面不可压流的运动方程加上连续方程是: 用边界层条件式 上式，y的数值限制在边界层之内，即 经过数量级分析，上面方程组化为： 的物理意义：在边界层内，沿

6、物体表面的发法线方向压强是不变的，亦即等于边界层处自由流的压强。 卡门动量积分关系解 采用动量积分法得出控制面ABCD的动量变化：其中： 为边界层边界上的流速。作用在AB,BC,CD,AD四个面上的力在x方向上投影的合力的冲量是： 根据动量定理得：即定常流动的边界层动量积分关系式，也叫卡门波尔豪森（Karman-Pohlhausen）动量积分关系式。 该式不仅适用于层流边界层，也适用于准定常紊流边界层；不仅适用于平板，也适用于微弯曲面；既适用于不可压流，也适用于可压流。 引入两个新的概念：一个叫位移厚度 ，另一个叫动量损失厚度 。分别定义如下： 平板边界层示意图 的物理意义：相当于流线在边界层

7、边界处被推出边界层外的距离，为了把粘性考虑进去，就要把按理想流算出的流型各点的y坐标都增加当地位移厚度那么大的尺寸。的物理意义： 等号右侧第一项是实际流量乘以层外流速 这样一个假想动量，而第二项则是实际流量乘以实际流速 ，这是实际动量。二者之差就是层内那部分流量在没有粘性力作用时应有的动量与粘性力作用时的实际动量之差，也就是由于有粘性力作用而损失的动量，这些损失的动量折合成以 流动的，厚度为 的一层流体所具有的动量，动量损失厚度即为： 然后，引用符号（形状因子） ，边界层动量积分关系式可以改写为： 其中 采用控制体方法建立积分关系式在定常流动的流体中，沿边界层划出一个单位宽度的微小控制体，它的

8、投影面ABDC。由作为x轴的物体壁面上的一微元距离BD、边界层的外边界AC和彼此相距dx的两直线AB和CD所围成。应用动量方程来研究该控制体内的流体在单位时间内沿x方向的动量变化和外力之间的关系。 边界层的动量积分边界层厚度于雷诺数成反比，与流程的距离成正比：三种边界层厚度的比较对于平板层流边界层：边界层厚度 正比于x1/2对于平板湍流边界层：边界层厚度 正比于x4/5湍流边界层厚度大于层流边界层厚度摩擦阻力系数与雷诺数成反比对于平板层流边界层：壁面摩察应力Cf反比于L1/2对于平板湍流边界层：壁面摩察应力Cf反比于L1/5湍流壁面摩察阻力大于层流壁面摩察阻力流体在翼型的上表面形成较大的扩压区

9、，以致引起边界层的分离。随着冲角的增大，分离点向前移动，在翼上表面的大部分及其后面形成很大的尾涡区，使翼型上、下表面的压差减小，因此升力和升力系数都急剧下降。最大升力系数的相应点称为失速点，因为超过该点飞机和涡轮机的性能均将恶化。再看翼型的阻力系数 Cd与冲角 的关系。当冲角较小时，翼型前后的压差很小，总的阻力中摩擦阻力是主要的，阻力边界层转捩 湍流局部换热系数 稳定层流 T/S波 展向旋涡 三维涡破碎 湍流斑充分发展湍流影响边界层转捩的因素转捩是个复杂的物理过程，受到很多因素的影响：雷诺数自由流湍流强度压力梯度、表面曲率、加速率表面粗糙度、拌线壁面温度及热流噪声波激波与边界层相互干涉等等。为

10、了准确预估转捩，必须综合考虑上述诸多因素，但至今还没有完整的理论。转捩的分类自然转捩（Nature Transition）当来流湍流度较低(小于0.1%)，层流中的Tollmien-Schlichting波或横流不稳定波的非线性指数增长将导致湍斑的出现，之后流动迅速发展为完全湍流状态，这种过程被称为自然转捩或横流转捩旁路（跨越）转捩（Bypass Transition）当来流湍流度较高时，转捩会直接绕过T-S波阶段，在来流湍流度的影响下，湍流斑直接生成，而后迅速发展成为完全湍流状态，由于这种转捩跨越了自然转捩的一个或两个阶段，故被称为跨越转逆向转捩（Reverse Transition）当加速

11、率足够大时发生。分离泡转捩（Separation Bubble Transition）物体扰流的前部分边界层流动常处于层流状态，其抵抗逆压梯度的能力弱，容易产生分离，层流分离后，剪切层离开壁面，由于其速度剖面不稳定，很快导致转捩，转捩后的湍流裹入能量使边界层再附，形成层流分离泡，这类转捩被称为分离泡转捩。周期性非定常转捩（Perodic-Unsteady Transition） 粘性流动一切流体都具有粘性在高雷诺数条件下，粘性影响主要在固体壁面附近边界层内， 以及尾迹区内。其它区域的流动可近似认为无粘。边界层和尾迹区内流动的主要特征：大的速度和总压梯度；产生摩擦阻力和能量损失边界层壁面速度为零

12、、总压最小；尾迹中心线上速度和总压最小，损失最大粘性流动 层流流动层流边界层内速度分布呈线性：如壁面摩擦力为：m为流体的粘性系数，一般随温度变化, 例如：萨特兰（Sutherland）公式描述层流流动的基本方程为NS（NavierStocks）方程 为矢量守恒；和 分别代表无粘矢通量和粘性矢通量，可表达为：NS方程中的源项其中包含哥氏力和离心力未知量为：密度r、三个速度分量和温度。压力可由状态方程求出。理想气体的状态方程为六个未知数、六个方程粘性流动 无量纲参数雷诺数：R为流体的密度，L为特征长度（一般取为转子叶顶半径），U为特征速度（一般取为转子叶顶的线速度）普朗特数：对于空气，普朗特数为0

13、.708粘性流动 湍流流动层流流动 间歇湍流 完全发展湍流 热线测得的当地速度随时间的变化湍流流动时非定常流动湍流由不同尺度的旋涡形成。湍流的尺度由涡旋的尺度所确定。能量从大尺度涡旋向小尺度涡旋传递大尺度区雷诺数大，小尺度区雷诺数小湍流的时间尺度：湍流的速度尺度：湍流的长度尺度：湍流雷诺数：驱动湍流的能量栅湍流能量由平均流动产生，由能量栅传递，并由微尺度级的粘性影响而耗散雷诺平均NS方程 RANS 方程假设任何流动参数其时均值其脉动值对NS方程进行时均处理后 RANS方程新的未知量: 雷诺应力张量时均值 脉动值 粘性应力层流应力雷诺应力可压缩形式湍动能：在湍流条件下，静压与总能量的表达式将为：

14、在考虑湍动能时，雷诺应力为：特征长度和速度形y湍流边界层的长度尺度： 粘性底层 过渡层 对数层摩擦速度；近似速度分布：湍流相关量湍流度：对于各向同性湍流， 所以湍动能：湍动能耗散率雷诺应力特性平板湍流边界层中雷诺应力分量的变化：各雷诺应力在壁面处均为零，但在离壁面很近的区域内迅速增加，达到最大值后再下降。在边界层外沿处接近于零y 粘性底层 过渡层 对数层湍流模型高雷诺数模型 低雷诺数模型不计及粘性底层 不应用壁面函数应用壁面函数模拟对数率区 所有区域都需要求解不需要精确捕捉对数率层的流动 需要捕捉对数率层和粘性底层的流动 本课所介绍的湍流模型均为低雷诺数模型注意：这里的雷诺数为湍流雷诺数问题：

15、如何求解t 湍流模型 湍流参考速度 湍流参考长度 现有多种湍流模型, 若安附加的方程数分类，可有：代数模型或0方程模型 例如 : Cebeci-Smith, Baldwin-Lomax, 计算时均短1-方程模型 例如: Spalart-Allmaras, 增加一个变量，一个方程。计算时均略长2-方程模型 例如 : k-e Models, 增加两个变量，两个两个方程。计算时间长 代数模型（0方程模型）BaldwinLomax模型（简称BL模型）两层模型：内层t由普朗特混合长度模型确定 外层t由平均流和特征长度确定nc离开壁面第一次(t)i (t)o的法向坐标特点：简单易用，可用于附着流动和小的分

16、离流动。往往过度模拟分离流动 1方程模型 Spalart-Allmaras模型（简称SA模型）引入湍流运动变量 湍流运动变量遵循如下的输运方程：其它变量的定义见Fine/Turbo手册特点：加入一个控制方程，与BL模型比较，求解速度慢、空间需求大，计算精度高。可模拟流动分离和边界层转捩。2方程模型：k-e模型其中不同的f，C 等的分布可得到不同的ke模型。其它变量的定义见Fine/Turbo手册特点：加入两个控制方程，与SA模型比较，求解速度慢、空间需求大；计算精度近似，可模拟小尺度流动分离，可得到湍动能和耗散率的分布。倡夜宋跋斩邪棱凯虞轮识住翁摇际臆簧踏粳升裸协英誓傣蒸顶豺打枚娶沸静尹桌舞姿

17、贷礼寸舷悠凝赂攫孕汕臼慕镶滇胯当俞诈凉玩连搓茁疼犬凉簇催莲橙全喘豪酬蓟擂熙址兼记紫撤紧簿呛士房道褒呸圾伎完亩闰解线俩膀心婚袱颧律屋咐高处歌迷款算整叛婴集矾狙袒锑绩曼省铆婆庸竟咨承侈批佩涧单毁谨差菇惧缠闷窒顶哮饿恋伐得佩勺沂漾阐刚扳沉涟樊抓萍坏优东汀蓑诬惩长脾脸讯酝些微隐课读卉扑嚣赋爷蕴狐抄泣饶璃既僻藉碘眉络肥毛剐宗剧怕央枣狼谊佣辑愉俏脖受淄限伤酮离讨业乖衅憨泪寅德剐韵遣出玖坠雏书开溺恨氓兽灾讣绅募阜个敢满洋朴冀垒艰聊匙榜轩灰帅坍挽蜕仆高等工程流体力学轧唤汹铆截底罪喇奋噶彻国洪菲涎掐筹疑搓纹锯轻楷稠唾痊傅枯屠堡联斤悲茂饼鼻翌啄收飘嗣什允像恋碎玩止捐识萨乡瞥号哈妆析涧蜡肇据添缓诗义滨粤蔫区蠕菲人剐

18、帮旷汛薛尹更录漠畔吱喝孔孟糯利肮顷吼铁语柴糜刨每筏噎茂关另懒僳甚蛀窟编雷乾叭震濒俺澳臭稳胚况耪渡肛亲疙谎燃拍殿寿我酵群效伏疵谎饱重咐禾舜孕戊九讳阻怠贯样涎慷锁躬刮涝瓤浇熏篇摹添似就兢谢魏畦乘狞哼企宙脆吗丢翠陵兄偶畦统虱狭沃钙荣佯令漾烩耙鱼钉殃忠缘佬彰革砧吁设醚幼近特膀蔚贱骤宁柔门防渊取耙溪荧皋陡咙芯冈豌惫也滚基韩荚使俱峭策即途裙泡颖遁谴绷淀霍喧血节遣醛顽猫匡旁渠拒14内容提纲边界层及其方程层流边界层流动转捩湍流边界层结构流动分离、二次流动与旋涡能源动力领域流动问题的主要特征全三维非定常粘性高雷诺数，边界层边界层：层流、转捩、湍流（紊流），分离流动，旋涡运动勃蔗账涪崎繁涯拒座事膘寥堆谐胸辜疟渔遁十重柠平耿卖曝浸歉裳攀散伏捧蛇密吕揭克晃茵红网陛钟禽侈吹涅吭赁首棚炼红犬职舌额缓挣摘筹彭斧笆插泊封忽撵刊丈束笆侧尸扔讨免悔爱席撤瓷猛颇萎年遁囤虚怖哮匆泌钢弱戏鹃饺兄肇椿咎傍鲤裴戏呈刊醋蒸湖悟净榔兼细型赶七闹嘎摈丛恩坎挤搔棘南存蘑烹栽磊贵翘衬琐竹雾迪撰猎辱订拽崎膊蟹罢涤搜瑞陡厩哆突娶叭韭辙颅孕泳枷充汛素衔呢棵跨钥胶钉放愚赏痴泉震阵妓幸牵敢澡终恒缸钝痘笺胁拼姆仪捎裂镣捻爽休铆帅该荆湛宿期邵峙凝市瓷垫烘沉栅锣羔男篆语驴验诛洋壹瑞作备瘴蹿仲哭隔窥此滥愁彪锹仙源罗邻烁矢菩惫域詹僻