边界层分析求解

合集下载

工程流体力学中的管流与边界层区分析

工程流体力学中的管流与边界层区分析

工程流体力学中的管流与边界层区分析在工程流体力学中,管流和边界层是两个重要的概念,它们在流体力学的研究和工程应用中具有重要作用。

管流指的是在管道内部通过液体或气体流动的现象,边界层则是指流体在固体壁面附近的一层流动区域。

管流是流体力学中常见的一种流动形式,它是在管道内部流动的一种现象。

在管流中,流体以一定的速度经过管道内壁面,流动过程中会产生一些特殊的现象和参数。

首先,流体在管流中会产生阻力。

阻力是指流体流动时受到的阻碍力,阻力大小与管道内的摩擦相互作用有关。

在管流中,阻力对于流体能量的转化和流速的损失是至关重要的。

其次,管流中的流速分布是不均匀的。

由于管道内壁面的摩擦阻力,流体的速度在不同位置存在差异。

一般情况下,流速在管道中心最大,在靠近壁面时逐渐减小,这种速度分布又被称为速度剖面。

了解流体在管流中的速度分布,对于管道设计、流量计算和其他工程应用都具有重要意义。

此外,管流中的流量计算也是一个重要的问题。

流量是指单位时间内流体通过管道的体积。

在工程实践中,准确计算管流的流量对于保证工程的安全和稳定运行至关重要。

相比之下,边界层是指流体在靠近固体壁面附近的一层流动区域。

边界层具有以下特点:首先,边界层内的流速分布是不均匀的。

由于壁面摩擦阻力的作用,流体的速度在边界层内存在梯度变化,即流速随离壁面距离的增加而增加。

其次,边界层内的流体粘性较大。

由于流体分子的力学相互作用力,流体在边界层内流动时会受到粘性阻尼的影响,流体粘性较大。

再次,边界层对壁面的摩擦力很重要。

流体在边界层内产生的摩擦力对于固体壁面的运动和稳定性具有重要影响。

最后,边界层的厚度会随着流体的速度变化而变化。

一般情况下,流体速度越大,边界层越厚。

边界层的厚度对于壁面附近的流动状态和特性具有重要影响。

在工程流体力学中,管流和边界层的分析和研究有助于了解和预测流体在管道中的流动行为。

例如,在管道系统的设计和优化中,需要对管流的速度分布、流量计算和阻力特性进行准确的评估和计算。

哈特曼边界层的初级稳定性分析

哈特曼边界层的初级稳定性分析

哈特曼边界层的初级稳定性分析董帅;林殿吉;吕玉坤【摘要】导电流体在法向外置磁场的作用下,在贴近壁面处会形成哈特曼边界层.哈特曼边界层的稳定性研究对电磁冶金过程和热核聚变反应冷却系统等相关设备的设计和运行都有着十分重要的意义.本文采用非正则模态稳定性分析方法,对两无限大绝缘平行平板内导电流体流动的稳定性进行了研究.通过在时间上迭代求解扰动变量的控制方程组和伴随控制方程组,获得了在磁场作用下初级扰动的增长情况及其空间分布形式,分析了磁场强度对最优扰动增长倍数Gmax、最优展向波数βopt和最优时刻topt的影响,并考察了上下两个哈特曼边界层之间的相互作用.结果表明,最优初始扰动的空间分布形式为沿着流场方向的漩涡,关于法向方向对称或者反对称.当哈特曼数Ha较大时(Ha>10),对称漩涡和反对称漩涡形式的初始扰动增长倍数基本相等;上下两个哈特曼边界层可以认为是彼此独立的,不会相互影响,此时最优扰动增长倍数Gmax与局部雷诺数R的平方成正比,相应的最优展向波数βopt和最优时刻topt均正比于哈特曼数Ha.当哈特曼数Ha较小时(Ha<10),反对称漩涡形式的初始扰动更为不稳定,其增长倍数大于对称漩涡的增长倍数,且上下两个边界层之间存在着一定的相互作用,并对整个流场的稳定性产生一定的影响.【期刊名称】《力学学报》【年(卷),期】2016(048)002【总页数】9页(P327-335)【关键词】哈特曼边界层;初级扰动;稳定性分析;非正则模态;数值模拟【作者】董帅;林殿吉;吕玉坤【作者单位】华北电力大学能源动力与机械工程学院,河北保定071003;华北电力大学能源动力与机械工程学院,河北保定071003;华北电力大学能源动力与机械工程学院,河北保定071003【正文语种】中文【中图分类】O361.52015–05–18收稿,2015–11–27录用,2016–01–05网络版发表.1)国家自然科学基金青年科学基金(11302076),河北省自然科学基金(A2014502047)和中央高校基本科研业务费(2014MS111)资助项目.2)董帅,讲师,主要研究方向:磁流体力学、流动稳定性分析.E-mail:********************.cn,*****************引用格式:董帅,林殿吉,吕玉坤.哈特曼边界层的初级稳定性分析.力学学报,2016,48(2):327-335Dong Shuai,Lin Dianji,Lu¨ Yukun.The primary stability analysis of Hartmann boundary layer.Chinese Journal of Theoretical and AppliedMechanics,2016,48(2):327-335冶金行业中双金属复层合金的制备、高分子聚合物的挤压成型、气动挤压塑料片、吹制玻璃等加工工艺中均会涉及到磁流体边界层内的流动问题[1-3].在这些领域中,导电流体被置于电磁场环境下,通过施加电磁力可以控制导电流体的运动状态[4-5].此外,在机翼上合理布置电磁场及相关设备,可以有效改变绕流边界层的结构,从而达到减震、减阻等效果[6-7].导电流体在法向外置磁场的作用下,在贴近壁面处会形成速度梯度很大、厚度很薄的边界层,即哈特曼边界层[8-9].哈特曼边界层的稳定性研究对电磁制动[10]、电磁搅拌[11]、晶体增长[12]、热核聚变反应冷却系统[13]等相关设备的设计和制造都有着十分重要的意义.Lundquist[14]最早提出了判断哈特曼边界层是否稳定(层流还是湍流)的参数——局部雷诺数R式中,Re为主流雷诺数,Ha为哈特曼数.Murgatroyd[15]对哈特曼边界层的稳定性进行了实验研究,发现局部临界雷诺数Rc位于150~250的范围内,哈特曼边界层内流体开始从层流向紊流转捩.Moresco和Alboussi´ere[16]实验得到的临界雷诺数Rc约为380,但是Lingwood和Alboussi`ere[17]采用正则模态线性稳定性分析得到的临界雷诺数Rc均明显高于实验结果.最近,Krasnov等[18]采用直接数值模拟的方法,通过在流场中施加随机背景噪声扰动来观测流场变化,得到了与实验值相近的结果,临界雷诺数Rc约为350.基于直接数值模拟结果,他们提出了如下转捩过程假说:(1)流向漩涡和哈特曼层流流动相互作用形成流向条纹结构;(2)流向条纹结构失稳,进而过渡到湍流.Gerard-Varet[19]通过非正则模态稳定性的分析方法对单一无限大平板上的哈特曼边界层进行研究,认为在壁面绝缘情况下,其最优扰动为流向漩涡,最优扰动增长倍数Gmax和局部雷诺数R呈平方关系,即Gmax≈5.65×10−4R2.Airiau和Castets[20]获得的结果与之类似.Krasnov等[21]随后对大哈特曼数Ha下,导电流体在方形截面管道内的流动进行了直接数值模拟.当哈特曼数Ha从0变化到400时,流动从湍流逐渐过渡到层流,哈特曼边界层厚度不断减小,上下两个哈特曼边界层之间的相互作用也逐渐减弱. 目前,对哈特曼边界层稳定性的研究主要针对较大哈特曼数Ha下的磁流体流动,对小哈特曼数Ha (Ha < 10)下流动的研究还不充分.本文对两块平行的无限大绝缘平板内导电流体流动的稳定性进行分析,重点考察初级扰动的增长情况,以及哈特曼数Ha和雷诺数Re的影响.平板流内上下两个哈特曼边界层之间可能存在着一定的相互作用,从而对整个流场的稳定性产生影响,这在以前的研究中并未引起足够的重视,本文将就此展开研究和讨论.1.1物理模型本研究考察的物理模型为不可压缩导电流体在两块无限大的绝缘平行平板之间的流动,在垂直流场的方向施加一个均匀恒定的静态磁场.其中,两平板的间距为d,磁场强度为B,层流流场中心线上的速度为U.平板内的流体在压力作用下流动,其体积流量保持恒定.1.2控制方程组本文是对低磁雷诺数下导电流体展开的研究,即忽略了流场对外界磁场的影响,只考虑外界磁场对流场的影响[22].其无量纲控制方程组和边界条件为其中,e = (0,0,1)为单位矢量,x,y,z分别代表了流向,展向和法向方向.u为速度矢量,p为压力,φ为电势,u,v和w分别代表了流向、展向、法向速度.进行无量纲化时,特征速度取哈特曼层流流动中心线上的速度U,特征长度L取平板间距的一半,即L = d/2;磁场强度特征值取B,相应的电势特征值为σUB,时间和压力的特征值分别为L/U和ρU2.式(2)中出现了两个重要的无量纲参数,哈特曼数和主流雷诺数本文研究的基本流是哈特曼层流流动,流场中的任一物理量ψ可以分解为其中,Ψ(z,t)表示基本流场的物理量,ψp(x,y,z,t)表示扰动变量.将式(3)代入式(2)中,减去基本流的控制方程组,并忽略关于扰动的高阶项,即可获得关于扰动的线性方程组.本研究主要采用非正则模态稳定性分析方法[23-27],考察扰动的增长模式及其获得的最大峰值,为了对扰动的增长情况进行分析,将扰动项展开成如下形式其中,α为流向波数,β为展向波数.为了定量的描述扰动的增长,选取变量E(t)代表扰动动能随时间的变化情况其中,上标“*”代表变量的共轭复数.初始扰动的增长或者衰减可由下式表示式中,T是所要考察的时刻,0为初始时刻,G称为扰动增长倍数.当哈特曼数Ha、局部雷诺数R、流向波数α、展向波数β和时刻T分别取某一定值时,扰动的分布形式仍有若干种可能,这些扰动都是控制方程组的解,目的是要从中找到某一特定的扰动分布,在T时刻获得扰动增长倍数G的最大值.这一特定的扰动分布被称之为最优扰动.1.3求解方法要获得G因子的最大值,数学上可以采用拉格朗日乘子法求解[28-29].即引入扰动变量的伴随场,通过耦合求解关于扰动的原始控制方程组和相应的关于伴随场变量的控制方程组,来获取最优扰动的空间分布形式和增长倍数.为了求解T时刻对应的最优初始扰动ˆu(z,0),可以采用循环迭代的数值算法[23],如图1所示.图1中顶标“~”代表变量的伴随场,顶标“ˆ”代表变量的原始场.首先由初始时刻的扰动分布ˆu(z,0),在时间轴上正向迭代求解控制方程组,即可获得在T时刻的扰动分布ˆu(z,T);该扰动分布经由适当变换可作为伴随场变量的初始值˜u(z,T),然后在时间轴上逆向迭代求解伴随场变量满足的控制方程组,可得˜u(z,0),并用该扰动分布去更新初始扰动ˆu(z,0).如此循环迭代下去,直至˜u(z,0) 和ˆu(z,0)满足一定的收敛标准,最后需要再次求解控制方程组,计算最优扰动ˆu(z,0)在T时刻的扰动动能,并代入式(6)中求解扰动增长倍数G.本文采用了高精度的伪谱方法求解无量纲方程组,在法向方向上采用切比雪夫配置点布置方式.这种网格节点布置方式的特点是中间稀疏,两边密集,从而可以较好地求解近壁面边界层内的流动.当哈特曼数较小时,磁场效果不明显,哈特曼边界层厚度较大,采用较少的网格节点(N = 64)可以精确求解控制方程组.当哈特曼数较大时,哈特曼边界层厚度变薄,需要相应增加节点个数(N = 128)以满足计算精度要求.对于一组确定的局部雷诺数R和哈特曼数Ha,某一初始扰动能够在某一流向波数α、展向波数β和时刻T下,得到所有扰动增长倍数G的最大值,称为最优扰动增长倍数Gmax,对应的流向波数α、展向波数β和时刻T分别称为最优流向波数αopt、最优展向波数βopt和最优时刻topt.对剪切流的研究发现,无论是平行平板流还是圆管中的流动,最优扰动的初始分布形式均为沿着流场方向的漩涡(α= 0),并且通常为关于法向反对称或者对称的漩涡[2631].所以本文选取初始扰动为沿着流场方向的对称漩涡或者反对称漩涡两种情况进行分析,并设定流向波数α= 0,只考虑了展向波数β的变化情况.图2和图3分别表示局部雷诺数R = 300,哈特曼数Ha = 5和15时,不同展向波数β下、不同扰动增长倍数G和时刻t之间的关系.当哈特曼数Ha = 5时,相同的展向波数β下,反对称漩涡在相应时刻的扰动增长倍数G均大于对称漩涡的扰动增长倍数G,即反对称漩涡获得了更大的扰动能量,但反对称漩涡达到最优扰动增长倍数Gmax的时刻大于对称漩涡.在哈特曼数Ha = 15时,相同的展向波数β下,初始对称漩涡和反对称漩涡取得了相等的最优扰动增长倍数Gmax,说明在较大哈特曼数Ha时,两种不同结构的初始扰动其增长情况基本相同.图4和图5分别表示局部雷诺数R = 300,哈特曼数Ha = 5和15时,最优扰动增长倍数Gmax和展向波数β随时间的变化情况.当哈特曼数Ha = 5时,反对称漩涡扰动动能达到最大值的时刻要大于对称漩涡,展向波数β要小于对称漩涡,这说明磁场对反对称漩涡的抑制效果更加明显.但是必须明确的是,反对称漩涡的最大扰动动能要大于对称漩涡.汇总结果如表1所示.当哈特曼数Ha = 15时,磁场对两种漩涡的影响相同,其最优扰动增长倍数Gmax和展向波数β随时间的变化基本相同,并在相同的时刻和展向波数β下取得最优扰动.两种情况下,随着展向波数β和时间t的增加,扰动增长倍数G均逐渐增大,直到G因子在某一时刻t 和某一展向波数β下获取峰值,即Gmax,相应的取得最优展向波数βopt和最优时刻topt.图6和图7分别显示局部雷诺数R = 300,哈特曼数Ha = 5和15时,最优初级扰动在初始时刻的空间分布.哈特曼数Ha = 5时,对称漩涡和反对称漩涡扰动的幅值是不同的;而哈特曼数Ha = 15时,对称漩涡和反对称漩涡的幅值相同,并且无论流向、展向还是法向速度,在流场中心区域均保持在零值上下,说明最优扰动主要分布在边界层附近,当哈特曼数越大(Ha = 15)时,这种趋势更加明显.从图中还可以看出,当哈特曼数Ha较小时,上下两块平板上的边界层之间有着一定程度的相互作用;而当哈特曼数较大时,两个边界层之间基本上不会互相影响.图8显示了初始扰动为对称漩涡和反对称漩涡两种情况下,局部雷诺数R = 300,500,700,1000时,最优扰动增长倍数Gmax随哈特曼数Ha的变化规律.在局部雷诺数R取值不同时,最优扰动增长倍数Gmax随哈特曼数Ha的主要变化趋势都是相同的,即哈特曼数Ha < 8时,最优扰动增长倍数Gmax随着哈特曼数Ha的增大而增大,当哈特曼数Ha > 8时,最优扰动增长倍数Gmax基本不变. 在局部雷诺数R相同,哈特曼数Ha < 8情况下,反对称漩涡的最优扰动增长倍数Gmax均大于对称漩涡.说明这种形式的扰动更为不稳定,而上下两个边界层之间的影响也是不可以忽略的.在哈特曼数Ha = 5附近,反对称漩涡的最优.最动扰动增长倍数Gmax出现了先增大后减小的特征.这表明哈特曼数Ha < 5时,磁场作用较弱,扰动增长主要受主流雷诺数的制约,即Gmax∝Re2;当哈特曼数Ha = 5时,最优扰动增长倍数Gmax达到最大值;当哈特曼数Ha > 5时,磁场强度的增加,扰动增长主要受局部雷诺数的制约,即Gmax∝R2.图9所示为不同哈特曼数Ha、不同漩涡结构时,最优扰动增长倍数Gmax和局部雷诺数R之间的关系,并与Gerard-Varet[19]的结果进行了对比.在哈特曼数Ha = 1,3,5,9时,对称和反对称漩涡的Gmax-R曲线出现了明显的差别,反对称漩涡的最优扰动增长倍数Gmax均大于对称漩涡的值.哈特曼数Ha = 10,15,20时,对称和反对称漩涡的Gmax-R曲线全部重合,即不仅哈特曼数Ha的大小不影响最优扰动增长倍数Gmax的大小,漩涡结构对最优扰动增长倍数Gmax的影响也非常有限.本文研究表明,在壁面绝缘情况下,最优扰动增长倍数Gmax和局部雷诺数R呈平方关系,Gmax≈6.12×10−4R2,该结果与Gerard-Varet[19]的结果Gmax≈5.65×10−4R2一致.因此,在哈特曼数Ha较大时,哈特曼边界层很薄,平行平板间的上、下两个边界层可以认为是相互独立的,不会互相影响.但是对于小哈特曼数,上下两平板间的哈特曼边界层距离较近,会互相影响.图10显示了不同漩涡结构时,最优展向波数βopt和哈特曼数Ha之间的关系.哈特曼数Ha≥10时,最优展向波数βopt和哈特曼数Ha呈线性增长关系;哈特曼数Ha < 10时,最优展向波数βopt和哈特曼数Ha呈指数增长关系.哈特曼数Ha < 10时,在相同的哈特曼数Ha下,反对称漩涡在更小的最优展向波数βopt下达到了最优扰动增长倍数Gmax.图11表明不同局部雷诺数R、不同漩涡结构时,最优时刻topt和哈特曼数Ha之间的关系.哈特曼数Ha < 10时,最优时刻topt和哈特曼数Ha二者增长关系不明显;哈特曼数Ha≥10时,最优时刻topt和哈特曼数Ha呈线性增长关系,且局部雷诺数R越大,斜率越大.哈特曼数Ha < 10时,反对称漩涡对磁场更加敏感,其最优时刻topt大于对称漩涡的最优时刻topt.哈特曼数Ha≥10时,磁场对两种漩涡的抑制效果相同.(1)当哈特曼数Ha较大时,磁场可以有效抑制扰动的增长,无论是对称漩涡还是反对称漩涡,其效果相同.在哈特曼数Ha较小的时候,对称漩涡和反对称漩涡对磁场的敏感性是不同的,反对称漩涡需更长的时间取得最优扰动增长倍数,说明磁场抑制了其增长速度,但是其取得了比对称漩涡更大的扰动增长倍数,即这种形式的扰动更为不稳定.(2)哈特曼数Ha > 10时,最优扰动增长倍数Gmax与局部雷诺数R的平方成正比,与哈特曼数Ha大小无关,对应的最优展向波数βopt和最优时刻topt均正比于哈特曼数Ha的一次方.(3)哈特曼数Ha > 10时,上下两边界层之间无相互影响;哈特曼数Ha < 10时,上下两边界层存在一定的相互作用,并且对整个流场的稳定性产生一定的影响.【相关文献】1 Mukhopadhyay S,Layek GC,Samad SA.Study of MHD boundary layer fl w over a heated stretching sheet with variable viscosity.International Journal of Heat and Mass Transfer,2005,48(21):4460-44662 Sharma R,Bhargava R,Bhargava P.A numerical solution of unsteady MHD convection heat and mass transfer past a semi-infinit vertical porous moving plate using element free Galerkin putational Materials Science,2010,48(3):537-5433 Rao JS,Sankar H.Numerical simulation of MHD effects on convective heat transfer characteristics of fl w of liquid metal in annular tube.Fusion Engineering andDesign,2011,86(2):183-1914 Molokov S,Moreau R,Moffatt HK.Magnetohydrodynamics-Historical Evolution and Trends.New Yorker:Springer,2007:231-2335 Moreau R.Magnetohydrodynamics.New York:Springer,1990:1-46 周本谋,范宝春,陈志华等.流体边界层上电磁力的控制效应研究.力学学报,2004,36(4):472-478 (Zhou Benmou,Fan Baochun,Chen Zhihua,et al.Flow control effects of electromagnetic force in the boundary layer.Chinese Journal of Theoretical and AppliedMechanics,2004,36(4):472-478 (in Chinese))7陈耀慧,范宝春,周本谋等.翼型绕流的电磁力控制.力学学报,2008,40(1):121-127 (Chen Yaohui,Fan Baochun,Zhou Benmou,et al.Electro-magnetic control of hydrofoilwake.Chinese Journal of Theoretical and Applied Mechanics,2008,40(1):121-127 (in Chinese))8 Hartmann J.Hg-Dynamics I,Theory of the laminar fl w of an electrically conductive liquid in a homogeneous magnetic field Kongelige Danske Videnskabernes Selskab,Matematisk-Fysiske Meddelelser,1937,15(6):1-289 Hartmann J,Lazarus F.Hg-dynamics II:Experimental investigation on the fl w of mercury in a homogeneous magnetic field Kongelige Danske Videnskabernes Selskab,Matematisk-Fysiske Meddelelser,1937,15(7):1-4510 Ha MY,Lee HG,Seong SH.Numerical simulation of threedimensional fl w,heat transfer,and solidificatio of steel in continuous casting mold with electromagnetic brake.Journal of Materials Processing Technology,2003,133(3):322-33911 Fujisaki K.In-Mold electromagnetic stirring in continuous casting.IEEE Transactions onIndustry Applications,2001,37(4):1098-110412 Dold P,Benz KW.Rotating magnetic fields Fluid fl w and crystal growth applications.Progress in Crystal Growth and Characterization of Materials,1999,38(1):39-58 13 Smolentsev S,Moreau R,Buhler L,et al.MHD thermoflui issues of liquid-metal blankets:Phenomena and advances.Fusion Engineering and Design,2010,85(7):1196-120514 Lundquist S.Studies in magnetohydr odynamics.Arkiv f¨ur Fysik,1952,5(4):297-34715 Murgatroyd W.Experiments on magnetohydrodynamic channel fl w.Philosophical Magazine,1953,44(359):1348-135416 Moresco P,Alboussi´ere T.Experimental study of the instability of the Hartmann layer.Journal of Fluid Mechanics,2004,504(10):167-18117 Lingwood RJ,Alboussi´ere T.On the stability of the Hartmann layer.Physics ofFluids,1999,11(8):2058-206818 Krasnov D,Zienicke E,Zikanov O,et al.Numerical study of the instability of the Hartmann layer.Journal of Fluid Mechanics,2004,504(10):183-21119 Gerard-Varet D.Amplificatio of small perturbations in a Hartmann layer.Physics of Fluids,2002,14(4):1458-146720 Airiau C,Castets M.On the amplificatio of small disturbances in a channel fl w with a normal magnetic field Physics of Fluids,2004,16(8):2991-300521 Krasnov D,Zikanov O,Boeck T.Numerical studyof magnetohydrodynamic duct fl w at high Reynolds and Hartmann numbers.Journal of Fluid Mechanics,2012,704(2):421-446 22 Davidson PA.An Introduction to Magnetohydrodynamics.Cambridge:Cambridge University Press,2001:117-11923 Dong S,Krasnov D,Boeck T.Secondary energy growth and turbulence suppression in conducting channel fl w with streamwise magnetic field Physics ofFluids,2012,24(7):07410124 李春曦,陈朋强,叶学民.含活性剂液滴在倾斜粗糙壁面上的铺展稳定性.物理学报,2015,64(1):014702 (Li Chunxi,Chen Pengqiang,Ye Xuemin.Stability of surfactant-laden droplet spreading over an inclined heterogeneous substrate.Acta PhysicaSinica,2015,64(1):014702 (in Chinese))25 Luchini P.Reynolds-number-independent instability of the boundary layer over a fla surface:optimal perturbations.Journal of Fluid Mechanics,2000,404(2):298-30926 Butler KM,Farrell BF.Three-dimensional optimal perturbations in viscousshearflw.Physicsof Fluids A:Fluid Dynamics,1992,4(8):1637-165027 Andersson P,Berggren M,Henningson DS.Optimal disturbances and bypass transition in boundary layers.Physics of Fluids,1999,11(1):134-15028 Schmid PJ,Henningson DS.Stability and Transition in Shear Flows.NewYorker:Springer,2001:515-51929 Reddy SC,Schmid PJ,Baggett JS,et al.On stability of streamwise streaks and transition thresholds in plane channel fl ws.Journal of Fluid Mechanics,1998,365:269-30330 Zikanov O.On the instability of pipe Poiseuille fl w.Physics of Fluids,1996,8(11):2923-293231 Gustavsson LH.Energy growth of three-dimensional disturbances in plane Poiseuille fl w.Journal of Fluid Mechanics,1991,224:241-260。

边界层流动特性分析

边界层流动特性分析
▪ 边界层方程
1.边界层方程是描述边界层内流体运动规律的基本方程,主要包括Navier-Stokes 方程和连续方程。 2.边界层方程的求解通常需要采用数值模拟或者近似解析方法,如普适函数法和相 似理论法。 3.边界层方程的研究对于揭示边界层流动的内在机制和预测流动行为具有关键作用 。
边界层概念与定义
▪ 边界层厚度测量方法
1.热膜风速计法:通过测量热膜上的热量传递来推算流体的速 度分布,从而得到边界层的厚度。 2.皮托管法:利用皮托管测量总压和静压差,计算出平均速度 ,再根据速度分布推导出边界层厚度。 3.激光多普勒测速技术(LDV):通过发射激光束并接收反射 光的多普勒频移信号,精确测量流场速度,进而确定边界层厚 度。
边界层分离
1.边界层分离是指当流体流过曲率半径较小的固体表面时,边 界层内的流体由于离心力的作用而从固体表面分离的现象。 2.边界层分离会导致流体在分离点后方形成涡旋,从而增加流 体与固体表面的摩擦阻力并影响流体的整体流动性能。 3.边界层分离的研究对于理解和控制流体流动中的能量损失、 噪声辐射以及流体机械的性能具有重要的实际意义。
边界层的分类
1.根据流体运动的特征,边界层可以分为层流边界层和湍流边 界层。层流边界层是指流体流动呈现有序、稳定的流动状态, 而湍流边界层则表现为无序、随机的流动状态。 2.根据流体与固体表面的相对运动关系,边界层还可以分为静 止边界层和动边界层。静止边界层是指固体表面静止不动时形 成的边界层,而动边界层则是指固体表面运动时形成的边界层 。 3.根据流体与固体表面的接触方式,边界层可以进一步细分为 光滑表面边界层和粗糙表面边界层。
边界层控制技术
1.边界层控制技术是通过改变边界层的流动特性来提高流体机 械效率、降低能耗和减少环境污染的一类技术。 2.常见的边界层控制技术包括流动诱导分离控制、湍流减阻技 术和热边界层控制等。 3.边界层控制技术在航空航天、能源、交通等领域具有广泛的 应用前景,对于推动相关行业的技术进步和可持续发展具有重 要作用。

大气边界层风场特征模拟与分析研究

大气边界层风场特征模拟与分析研究

大气边界层风场特征模拟与分析研究在大气科学领域中,大气边界层是指地球表面与上层大气之间的过渡层,对于气象预测、空气质量评估等领域具有重要意义。

而其中的风场特征是其研究的重点之一。

本文将探讨大气边界层风场特征的模拟与分析方法。

一、背景介绍大气边界层的风场特征是受多种因素综合影响的,在不同的地理环境中表现出差异性。

由于大气边界层风场的复杂性,为了更好地理解和模拟其特征,研究者们提出了不同的方法和技术。

二、数值模拟方法数值模拟方法是研究大气边界层风场特征的常用手段之一。

该方法通过数学模型和计算机算法来模拟大气边界层风场的演变过程。

常见的数值模拟方法包括气象模式、雷达回波模拟等。

1. 气象模式气象模式是通过一系列参数和方程来描述大气运动、能量传递和物质运输的数值模型。

通过运行气象模式,可以模拟大气边界层的风场特征。

常用的气象模式包括欧洲中期天气预报中心的欧洲中期天气预报模式(ECMWF)、美国国家环境预报中心的全球预报系统(GFS)等。

2. 雷达回波模拟雷达回波模拟是利用雷达回波数据和数学模型,通过计算和推算得到大气边界层的风场信息。

通过分析雷达回波的特征,可以获取大气边界层风场的分布和运动情况。

三、实测资料分析方法除了数值模拟方法外,实测资料的分析也是研究大气边界层风场特征的重要手段之一。

通过各种地面、航空、卫星观测站点所获取的实测数据,可以对大气边界层的风场特征进行分析。

1. 地面观测站点地面观测站点是通过建立气象观测站网络,采集并记录大气各种要素的实测资料。

通过对地面观测站点资料的分析,可以得到不同地理环境中大气边界层风场的特征。

2. 航空观测资料航空观测资料是通过飞机或无人机等航空平台所采集的数据。

通过对航空观测资料的分析,可以获取大气边界层风场在垂直方向上的变化情况,进而揭示其垂直结构特征。

3. 卫星观测资料卫星观测资料是通过卫星对地球表面进行遥感探测所获取的数据。

卫星观测资料具有广覆盖区域、高时空分辨率的特点,通过对卫星观测资料的分析,可以更全面地认识大气边界层风场特征。

边界层的形成与流动特性分析

边界层的形成与流动特性分析

边界层的形成与流动特性分析边界层是指在固体物体表面和流体之间的一个细小区域,这个区域内由于粘性效应的存在,流体流动速度逐渐从静止状态递增,直到达到与远离固体的自由流动速度相同的状态。

边界层形成与流动特性的分析对于许多领域具有重要意义,包括航空航天、工程设计、地质地球物理等。

边界层形成的过程可以通过物理原理和数学模型进行解释。

当流体在静止的固体表面上流动时,由于粘性作用,流体分子与静止物体表面接触后减速,形成运动速度减小的速度梯度。

这种速度梯度会逐渐向上游传播,形成一个层状结构,即边界层。

边界层的厚度取决于流体的速度、密度、粘性以及固体表面的粗糙度等因素。

边界层的流动特性与其形成过程紧密相关。

边界层的流动可以分为层流和湍流两种形式。

在边界层的初始部分,流体分子按层状结构有序运动,形成层流流动。

然而,在远离边界层的区域,由于速度梯度的变化,流体分子开始混乱运动,形成湍流流动。

层流和湍流的比例可以通过雷诺数来描述。

当雷诺数较小时,层流占主导地位;而当雷诺数较大时,湍流占主导地位。

边界层的流动特性也会受到影响因素的改变而发生变化。

例如,当固体表面的粗糙度增加时,边界层的湍流程度也会增加。

此外,边界层也受到来流速度的影响。

当来流速度增加时,边界层的厚度会减小,流动的剪切力也会增加。

这对于工程设计和流体力学的分析非常重要,因为它可以影响到一些重要的参数,例如风的压力、阻力、换热和质量传输等。

边界层的形成与流动特性分析对于实际问题的研究具有重要意义。

例如,在航空航天工程中,了解边界层的形成与流动特性可以帮助设计更加优化的机翼和机身,减小空气阻力,提高飞行性能。

在工程设计中,通过分析边界层的形成与流动特性可以改善传热和传质过程,提高设备的效率。

在地质地球物理研究中,边界层的分析可以帮助解释地下流体运动和岩石物理现象。

综上所述,边界层的形成与流动特性分析对于许多领域的研究和应用具有重要意义。

通过物理原理和数学模型的分析,我们可以更好地理解流体与固体表面的相互作用,并优化相关系统的设计与运行。

边界层理论(Boundary layer theory)--西安交大

边界层理论(Boundary layer theory)--西安交大

)之外的流体速度就形成:润湿→附着→内摩擦力→减速→梯度
边界层内:沿板面法向的速度梯度很边界层外:不存在速度梯度或速度梯度
流体在平板上流动时的边界层:
流动边界层:存在着较大速度梯度的流体层区域,即流速降为主体流速的99%以内的区域。

边界层厚度:边界层外缘与壁面间的垂直距离。

层流边界层:在平板的前段,边界层内的流型为层流。

湍流边界层:离平板前沿一段距离后,边界层内的流
直管内:流体须经一定的距离才能形成稳定的边界层。

由于总流量不变,中心流速增加。

边界层占据整个管截面。

与物体的长度相比,边界层的厚度很小;边界层内沿边界层厚度的速度变化非常急边界层沿着流体流动的方向逐渐增厚;
边界层中各截面上的压强等于同一截面上在边界层内粘滞力和惯性力是同一数量级边界层内流体的流动存在层流和紊流两种
圆柱后部:猫眼
扩张管(上壁有抽吸)
B
S′
A
涡,这种旋涡具有一定的脱落频率,称为卡门涡街.
湍流产生的原因:
湍动强度
在模型实验中,模拟湍流,要求雷诺数和湍动强边界层的转变、分离以及热量和质量传递系数等
依微分方程的个数:零方程模型、一方
FLUENT软件在化学处理领域主要可应用 于:
燃烧 干燥 过滤 传热和传质 材料处理 混合 反应 分离 蒸馏 喷射控制 成型 焚化 测量/控制 聚合 沉淀 通风



















大气边界层特征的观测与分析

大气边界层特征的观测与分析

大气边界层特征的观测与分析大气边界层是地球大气层中与我们的日常生活和环境息息相关的重要部分。

它是靠近地球表面的一层大气,厚度通常在几百米到几千米之间,其特征对天气、气候、污染物扩散等都有着显著的影响。

因此,对大气边界层特征的观测与分析具有重要的科学意义和实际应用价值。

大气边界层的特征主要包括温度、湿度、风速、风向等气象要素的垂直分布,以及湍流运动等。

为了观测这些特征,科学家们采用了多种手段和方法。

一种常见的观测方法是使用气象塔。

气象塔通常高达几十米甚至上百米,在不同高度上安装了各种气象传感器,如温度传感器、湿度传感器、风速仪等,可以实时获取不同高度处的气象数据。

通过对这些数据的分析,我们能够了解大气边界层中气象要素随高度的变化情况。

例如,在白天,由于太阳辐射的加热作用,地面温度升高,空气受热上升,形成对流,导致温度在垂直方向上的分布呈现出明显的梯度。

除了气象塔,飞机观测也是一种重要的手段。

飞机可以在大气边界层中飞行,并携带各种测量仪器,获取大范围、高空间分辨率的气象数据。

然而,飞机观测的成本较高,且受飞行条件和航线的限制。

近年来,随着卫星遥感技术的发展,为大气边界层的观测提供了新的视角。

卫星可以通过测量大气中的水汽、温度等参数,反演得到大气边界层的特征信息。

但卫星观测也存在一定的局限性,比如分辨率相对较低,对某些细节特征的捕捉能力有限。

在获取了大量的观测数据之后,接下来就是对这些数据进行分析。

数据分析的方法多种多样,其中常用的有统计分析和数值模拟。

统计分析是通过对观测数据进行整理、计算,得出各种统计量,如均值、方差、相关性等,从而揭示大气边界层特征的一般规律。

例如,通过对多年的温度观测数据进行统计分析,可以发现大气边界层温度的季节变化和年际变化特征。

数值模拟则是利用计算机模型来模拟大气边界层的物理过程。

模型中考虑了大气的热力学、动力学方程以及各种物理过程,如辐射、湍流交换等。

通过输入观测数据和边界条件,模型可以预测大气边界层的演化和变化。

大气边界层高度测量方法研究及比较分析

大气边界层高度测量方法研究及比较分析

大气边界层高度测量方法研究及比较分析在大气科学研究中,大气边界层高度的测量是至关重要的。

边界层高度是指大气中温度、湿度、风速等物理特性发生显著变化的区域,对于气象预报、空气污染监测、气候变化研究等都具有重要意义。

本文将研究和比较几种常用的大气边界层高度测量方法。

一、大气边界层高度的重要性大气边界层是指气象学上特指地球大气直接受大气地面相互作用影响的底部大气层。

边界层高度的变化与大气中的温度、湿度、气压和风速等物理参数的垂直分布密切相关。

准确测量大气边界层高度有助于了解大气层的垂直结构和变化,进而提高天气预报的准确度。

二、常用的大气边界层高度测量方法1. 气象雷达法气象雷达法是通过雷达测量附近的云层、积雨云或大气湿度变化来确定边界层高度。

该方法具有测量范围广、实时性强的优点,适用于大面积的边界层高度探测。

2. 气象探空法气象探空法是通过使用气象探空仪,将探空仪悬挂在气球上并释放到大气中,记录高度和大气参数的变化。

这种方法能够提供边界层高度的垂直剖面,并且在可用性和准确性方面具有优势。

3. 激光测量法激光测量法是通过激光束向大气中发射,通过测量激光的散射和反射来确定大气边界层高度。

该方法具有非接触式测量和高精度测量的优点,但受到天气和环境条件的限制。

4. 遥感卫星法遥感卫星法通过使用卫星传感器来观测大气的物理特性,如温度、湿度和云量等,从而确定边界层高度。

这种方法具有全球范围内的监测能力,但分辨率较低,无法提供较高精度的边界层高度数据。

三、不同方法的比较分析以上介绍了几种常用的大气边界层高度测量方法,它们各自具有优缺点。

选择合适的方法取决于应用的具体需求。

例如,气象雷达法适用于大范围的边界层高度探测,而遥感卫星法则更适合全球范围的监测。

此外,测量方法的实用性和准确性也需要考虑。

气象探空法具有较高的精度和垂直分辨率,但需要人力物力投入较大。

激光测量法在非接触式测量和高精度测量方面具有优势,但受到天气和环境条件的限制。

高等热质交换技术 边界层的分析

高等热质交换技术  边界层的分析

1.速度边界层(流动边界层) 速度发生明显变化 的 流体薄层。
流动边界层厚度 :
u 0.99u
流场的划分: 主流区:y< 理想流体
紊流核 心
u 边界层区: y 速度梯度存在与粘性力的作用区。
边界层的流态: 层流边界层、过渡区、紊流边界层
2
临界距离xc : 边界层从层流开始向紊流过渡的距 离。其大小取决于流体的物性、固体壁面的粗糙度等 几何因素以及来流的稳定度,由实验确定的临界雷诺 数Rc给定。
对流项 导热项
对流项进行比较分析
1.假定Δ/Δ ≥ 1,在这种情形下,导热项的系数就应 该等于 1 1 t P r Re 2 t
由于
1 Re
t
1 Re P r
Pr t
这种情形下成立的条件只能是Pr ≥ 1,意味着这样的流 体,流动边界层厚度大于或等于温度边界层厚度.
温度的无量纲参数采用: T Tw T Tw

变化范围0-1
热边界层厚度δ t度不一定等于流动边界层厚度δ, δ=LΔt , Δt是某个百分数。
无量纲参数:
x X L
y Y L
1 2 2 2 U V 2 2 2 X t Y c p Lu t X Y
导热项是偏微分方程中最高阶项,必须考虑,但 2 2 2 远小于 可以忽略不计。 2 2 X Y
无量纲能量方程
1 1 2 U V 2 X t Y Pr Re 2 Y t
有量纲边界层能量 方程
T T 2T u v x y cp p y 2
u v 0 x y
边界层对流 换热微分方 程组:

传热学(chapter7)-边界层 数量级分析

传热学(chapter7)-边界层 数量级分析

t − tw = 0.99 t f − tw
1 >> δ
∂ 2u ∂u ∂p ∂u ρ = − ∂x + µ u ∂x + v ∂y ∂y 2
∂ 2u ∂u ∂p ∂u ρ = − ∂x + µ u ∂x + v ∂y ∂y 2
∂ 2u ∂u ∂p ∂u ρ = − ∂x + µ u ∂x + v ∂y ∂y 2
∂ 2u ∂u ∂u u +v =ν 2 ∂y ∂y ∂x
∂ 2t ∂t ∂t ρc p u ∂x + v ∂y = λ ∂y 2 Fra bibliotek物理现象
(1)当粘性流体在壁面上流动时,由于粘性的作用, 在贴附于壁面的流体速度实际上等于零,即在 y=0,u=0处;。 (2)此后随 y 增加 ,u增加 。
(3)壁面法向(y 向)的速度分布,如上图所示。 经过一个薄层后u接近主流速度。 (4)
空气 u∞ = 10m / s 平壁
2mm
µ =0
第六章对流换热的基本方程1什么是边界层流动和热2为什么研究边界层流动和热变化的分布区域边界层外流动和热稳定物理现象1当粘性流体在壁面上流动时由于粘性的作用在贴附于壁面的流体速度实际上等于零即在增加u增加向的速度分布如上图所示
第六章
对流换热的基本方程
1、什么是边界层(流动和热)
2、为什么研究边界层(流动和热变化的 分布区域,边界层外流动和热稳定)

大气边界层特征分析

大气边界层特征分析

大气边界层特征分析1.1大气稳定度特征利用A地区2019年地面常规观测资料,对A地区2019年夏、冬两季每天7:00和19:00的大气稳定度级别进行频数统计。

结果如下:表1.1 A地区2019年夏冬两季大气稳定度频数统计表夏冬7:00 19:00 7:00 19:00A 0 0 0 0B 65 0 0 0C 9 0 59 0D 16 31 30 20D~E 0 22 0 22E 0 32 0 22F 0 7 0 25由表1.1可以看出,夏冬两季每天7:00的状态均为不稳定状态(B、C)和中性状态(D);19:00的状态均为中性状态(D)和稳定状态(D~E、E、F),中性状态(D)在早晚都会出现。

但是在7:00并未出现极不稳定状态(A),这是因为早晨太阳刚刚出来,辐射非常弱,到达地面的辐射能力很少,使得湍流运动达不到极不稳定状态的量值。

夏季7:00出现中性不稳定(B)的天数最多,为65天,冬季为0天,出现这种现象是因为夏季太阳直射点在北半球,昼长夜短,因此日出时间早,到7:00时的太阳高度角最高达到27º57´,最低也可以达到21º34´,此时太阳辐射到地面的能量加热地面使得地面增温,湍流运动较强,因此造成大气层结达到中性不稳定(B)的程度。

冬季19:00出现中等稳定(F)的天数最多,为25天,其主要是因为冬季昼短夜长,太阳在19:00时早已落山,此时地面由于接收不到太阳辐射能量的补给,温度开始降低,地面的冷却通过湍流交换向上传播而形成逆温层,此外,大气辐射冷却同样加快了逆温层的形成。

此时,边界层底层为稳定层结,当以上所述的湍流交换过程以及大气辐射冷却过程越强时,层结越稳定。

1.2边界层高度特征对夏冬两季每天7:00和19:00的不同稳定度级别下的边界层高度进行频数统计,如表1.2~1.3。

表1.2夏季不同稳定度下边界层高度频数统计表7:00 19:00不稳定(B、C)中性(D)中性(D)稳定(D~E、E、F)0-500(米)32 10 1 5500-1000(米)38 6 13 151000-1500(米) 4 0 4 191500-2000(米)0 0 7 82000以上(米)0 0 6 14 总数74 16 31 61平均582 494 1336 1542表1.3 冬季不同稳定度下边界层高度频数统计表7:00 19:00不稳定(C)不稳定(D)中性(D)稳定(D~E、E、F)0-500(米)57 30 7 19500-1000(米) 2 0 8 421000-1500(米)0 0 5 61500-2000(米)0 0 0 22000以上(米)0 0 0 0 总数59 30 20 69平均244 203 795 727对比表1.2和表1.3,夏季边界层的月平均高度高于冬季,一年中夏季的太阳辐射能量最强,对流能量大,造成边界层高度较高。

平板边界层实验报告

平板边界层实验报告

平板边界层实验报告引言平板边界层实验是一种常见的流体力学实验方法,用于研究在流体与固体界面发生的各种现象和特性。

通过实验可以获取边界层厚度、速度剖面、摩擦系数等参数,对于理解流体边界层的特性具有重要意义。

本实验报告将详细介绍平板边界层实验的原理、实验装置、实验过程和实验结果,并对实验结果进行分析和讨论。

实验原理在实验中,我们使用平板边界层实验装置对流体的边界层进行研究。

其原理基于以下几点:1.边界层理论:边界层是指流体流动过程中处于流体与固体物体之间的一层流动区域,其特点是速度梯度较大、流动剪切应力较高。

边界层的特性对于流体的运动、传热和传质等过程具有重要影响。

2.平板边界层:平板边界层是指位于平板表面附近的边界层,它是边界层研究中最常见的情况之一。

通过对平板边界层的研究,可以深入理解边界层的结构、特性及其对流体流动的影响。

3.流动速度剖面:边界层中流体的速度随距离平板表面的距离而变化,一般呈现一定的速度剖面形态。

通过测量流体速度剖面,可以确定边界层的厚度和速度分布特性。

实验装置实验装置由以下几个主要部分组成:1.平板:平板用于产生平板边界层。

通常采用光滑的表面,材质多为金属或塑料。

2.流体:实验中常使用空气或水作为流体介质。

流体通过输送装置注入到实验装置中。

3.流量计:流量计用于精确测量流体的流量,以保证实验条件的准确性。

4.速度测量装置:速度测量装置用于测量流体在平板边界层中的速度。

常见的测量方法包括热线法、激光多普勒测速法等。

5.数据记录系统:数据记录系统用于记录实验过程中获得的各项数据,包括流体流量、速度剖面等。

实验步骤本实验的具体步骤如下:1.准备工作:清洁实验装置,确保平板表面光滑且无杂质。

2.实验装置搭建:按照实验要求搭建实验装置,包括安装平板、连接流体输送装置和速度测量装置。

3.流体注入:启动流体输送装置,将流体注入实验装置中,并调节流量控制阀以控制流体的流量。

4.测速:使用速度测量装置对流体在平板边界层中的速度进行测量。

k-ω低雷诺数

k-ω低雷诺数

k-ω低雷诺数
k-ω模型是一种流体动力学模型,用于描述湍流流动。

它包含两个方程:k方程和ω方程,其中k表示湍动动能,ω表示湍动耗散率。

低雷诺数(Low Reynolds Number)指的是在较小的雷诺数下进行流动分析,即在这种情况下湍流效应较小,流动主要由粘性效应决定。

在k-ω模型中,低雷诺数下的流动可以通过以下步骤进行分析:
1.确定边界条件:在进行低雷诺数分析之前,需要确定流体
流动的几何形状和边界条件。

这通常包括流动域的形状、
入口速度、壁面条件等。

2.解析求解:使用k-ω模型的方程进行数值求解,得到k和
ω的分布。

这可以通过使用计算流体力学(CFD)软件或
其他数值方法来实现。

3.边界层分析:对于低雷诺数流动,边界层的发展和特性对
整个流动的影响很大。

因此,在低雷诺数分析中,通常会
对边界层进行详细的分析,以了解壁面摩擦、边界层厚度
等参数的变化。

4.结果解释:根据计算结果,进行相应的结果解释和分析。

可以关注转捩点位置、流动分离和再附着等方面的特征。

需要注意的是,低雷诺数下的湍流流动具有一些特点,例如湍流的产生和维持主要依赖粘性效应,湍流的能量耗散较小等。

因此,在低雷诺数下进行流动分析时,需要特别关注流动特性,
并结合实际问题和结果来解释流动现象。

大气边界层结构及其影响因素分析

大气边界层结构及其影响因素分析

大气边界层结构及其影响因素分析大气边界层是介于地球表面和大气上层之间的一个区域,是大气和地面之间的交界处,气象学上也称为对流层。

边界层的高度通常在几百米到几千米之间,其中包含了地面的摩擦层和大气上部的惯性层。

边界层的结构和层高是大气研究和气象预报中的重要参数,影响着大气的动力学和化学过程。

大气边界层结构主要由气温、风速、湿度、气压和大气组分等因素控制,其中,气温和风速是影响边界层结构最重要的因素。

气温是影响边界层结构最重要的因素之一,它对于边界层的发展和演化有着重要的影响。

在白天,地面受太阳辐射的热量加热,导致地面温度升高,热空气向上流,形成了较强的对流运动。

在这种情况下,边界层往往呈现出温度递减、风速递增、湿度递减等明显特征。

然而,夜间地面会散放出热量,导致温度逐渐降低,空气不再产生对流运动,边界层开始稳定,呈现出温度递增、风速递减、湿度递增等特征,这种现象我们通常称之为夜间稳定层。

风速是另一个影响边界层结构的重要因素。

由于边界层受到地面摩擦的影响,近地表层的风速会逐渐减小,形成风速递减层,这种现象通常表现为物体在不同高度风速不一致,例如风筝在地面上飘荡,而在较高的高度上则会遇到更强的风速。

在边界层中,风速的垂直梯度也很重要,垂直风速梯度的大小决定了边界层的动力学特征和边界层混合的程度。

湿度也是影响边界层结构的重要因素之一。

湿度的变化会影响边界层的稳定度和边界层内部的热和水汽交换,同时也影响着大气的化学反应过程。

在潮湿的条件下,水汽会加强大气中的凝结和降水过程,从而影响边界层的垂直分布和动力学特征。

气压和大气组分也对边界层结构产生了一定的影响。

气压的变化会影响着空气的压强梯度力、重力和惯性力等,进而对流场和温度场发生改变。

大气组分中水汽、二氧化碳、氧气、氮气等组分在边界层中的扩散和混合也会影响边界层结构的演化和化学反应过程。

总之,大气边界层的结构和演化是由多种因素决定的。

气温和风速是边界层结构最为重要的两个因素,湿度、气压和大气组分等因素也都对边界层结构和化学反应过程产生了重要的影响。

边界层理论

边界层理论
19世纪中,随着航海、水利工程等的迅速发展,流体力学的另一个重要分支,研究不可压缩粘性流体流动的 水力学得到很大的发展。它是建立在大量实验测量的基础上。当时如哈根、泊肃叶、雷诺等用实验研究水和其他 粘性流体在管道和槽渠中流动时的阻力和压强损失问题、得到的有关粘性流体的实验研究成果,有助于解决某些 工程实际问题。但由于水力学在理论指导上的不足,由实验成果得出的经验公式和半经验理论公式有一定的局限 性。于是在19世纪中叶产生了粘性流体运动的理论,1827年,纳维尔在欧拉运动微分方程中加上粘性项,第一个 得到粘性流体运动微分方程。1846年,斯托克斯严格地导出了这个方程,称为纳维尔-斯托克斯方程,简称N-S方 程。虽然N-S方程对粘性流体流动问题的研究分析有所帮助,但对这个方程数学上的求解是十分复杂和困难的。 1851年,斯托克斯对N-S方程作了某些简化,略去方程中的惯性项,也就是在非常缓慢的流体流动条件下,计算 出球体在流动的粘性流体中所受到的阻力。
边界层方程组
边界层方程组
不可压缩流体在大雷诺数的层流情况下绕过平滑壁面的情况。在此考虑二维定常不可压缩流动。规定沿物体 壁面的方向为x轴,垂直于壁面的方向为y轴。由于边界层厚度δ比物面特征尺寸L小得多,因此对二维的忽略重 力的纳维-斯托克斯方程逐项进行数量级分析,在忽略数量级小的各项后,可近似认为边界层垂直方向的压力不 变,从而得到层流边界层方程组为:
发展
1907年,布拉修斯成功地应用边界层理论计算在流体中运动物体的摩擦阻力。1921年,卡门和波耳豪森提 出了边界层动能积分方程,以计算边界层问题,这个方程经霍尔斯坦-博伦(1940)和瓦茨进行简化和改进,到 现在还被广泛应用。另外边界层动能积分方程和热能积分方程分别由莱本森和弗兰克尔提出。这三个边界层的近 似计算方法使边界层理论在工程界中很快地推广开来。1925年,普朗特提出的混合长度理论和1930年卡门提出的 相似性理论,将边界层理论推广到紊流边界层、射流和物体后的尾迹流中去。从层流向紊流的转捩现象是流体动 力学中的基本现象。早在19世纪末,雷诺就首先对转捩现象进行了研究。1914年,普朗特做了著名的圆球实验, 正确地指出:边界层中的流动可以是层流的,也可以是紊流的,还指出边界层分离的问题,因此计算阻力的问题 是受这种转捩支配的。从层流向紊流的转捩过程的理论研究,是以雷诺的假设为基础的,即承认紊流是由于层流 边界层产生不稳定性的结果。1921年,普朗特开始进行转捩的理论研究,1929年获得成功。当时托尔明从理论上 算出零冲角平板转捩的临界雷诺数,后被别人所进行非常仔细的实验所证实。稳定性理论能够考虑到对转捩有影 响的压强梯度、抽吸、马赫数和传热等许多因素。这个理论已得到很多重要的应用,如设计阻力非常小的层流翼 型。

大气边界层案例分析

大气边界层案例分析

⼤⽓边界层案例分析⼤⽓边界层案例分析1. 由下图分析晴天⽩天和夜间典型的风温垂直分布。

分析:⼤⽓边界层中温度层级起着重要作⽤,层结的稳定与否决定了湍流的强弱,也就决定了边界层中⽓象要素的垂直分布(廓线)。

图1.3.1是晴天⽩天和夜间典型的理想的风温垂直分布。

在贴近地⾯的薄⽓层内(近地层SL),⽩天由于地⾯强烈受热,形成贴近地⾯⼤⽓中超绝热温度递减率,⽽反映在位温上,即是/0z θ??<,风速则随⾼度递增。

再向上,在边界层的⼤部分范围内θ有⼀个不随⾼度变化的⽓层,风速也是如此,相应温度呈绝热下降,我们称之为混合层(ML)。

其原因是强烈的湍流混合使风、位温等垂直梯度减⼩,造成均匀分布。

在边界层以上的⾃由⼤⽓(FA)中,温度恢复为⾃由⼤⽓的递减率,位温则随⾼度⽽增,风则接近地转风速。

在⾃由⼤⽓与边界层间有⼀个过渡区域,其中各⽓象要素由边界层值逐渐过渡到⾃由⼤⽓。

此层称为夹卷层(EZ),在夹卷层中,发⽣着复杂的物理过程,从边界层中受热上升的⽓块可以穿透边界层与⾃由⼤⽓间的逆温⽽进⼊⾃由⼤⽓。

同样,湍流、重⼒波等亦可使⾃由⼤⽓中具有较⾼位温的⽓块进⼊边界层,这种过程称为夹卷,在夹卷层中即进⾏着边界层与⾃由⼤⽓间的各种交换。

典型夜间的风温廓线从图 1.3.1可看出在地⾯附近有⼀个逆温层,亦即稳定边界层(SBL),在T 和θ上均体现出来,这是由于地⾯强烈冷却造成地⾯温度低于⼤⽓造成,在其上则是⼀个θ随⾼度变化很⼩的“残留层”(RL),从成因来说,⽩天的对流边界层在夜间由于地⾯降温⽽在近地⾯形成逆温,但上部⼀段却保持着⽩天混合层的特征,使θ近于随⾼度不变,并且在残留层与⾃由⼤⽓间仍有顶盖逆温(CI),但残留层由于逆温层的存在已与地⾯脱离关系,其中湍流得不到发展的动⼒⽽逐渐衰减。

夜间边界层的风场由于夜间湍流弱,湍流摩擦⼒减⼩,风速与⽩天⽐得到加强,因⽽呈现出有最⼤风在某⾼度出现。

2. ⼤⽓边界层是与⼈类活动关系最为密切的⼀层,⼤⽓边界层具有哪些基本特点?分析:⼤⽓边界层的基本特点有:(1)运动的湍流性⼤⽓边界层有别于其上的⾃由⼤⽓的基本特点就是其运动的湍流性。

边界层厚度计算方法详述

边界层厚度计算方法详述

边界层厚度计算方法详述
1.拉格朗日方法:
拉格朗日方法是一种基于粒子运动的方法,其基本假设是跟踪粒子的运动轨迹,通过观察粒子运动状态的变化来计算边界层厚度。

具体步骤如下:
-使用流场模拟软件进行计算,获取流体在不同位置上的速度分布。

-设定一个壁面附近的初始点P,在该点上标记一粒子,并设定初始速度。

-模拟粒子在流体中的运动轨迹直到离开边界层或到达边界层厚度要求的位置。

-记录每个位置上粒子的速度和距离
-根据速度和距离的变化情况,推断出边界层厚度。

2.欧拉方法:
欧拉方法是一种基于流场的全局分析方法,通过分析流体的整体特性来计算边界层厚度。

具体步骤如下:
-建立流体的速度剖面模型。

常见的模型包括线性剖面模型、幂函数剖面模型和指数剖面模型等。

-使用流场模拟软件进行数值计算,得到固定位置上的速度分布。

-分析速度分布和壁面摩擦阻力等特性,进而计算边界层的厚度。

值得注意的是,在具体计算边界层厚度时
-不同流体属性的影响:边界层厚度与流体的粘性相关,较高的粘性
会导致边界层厚度增加。

-壁面状况的影响:粗糙的壁面会增加摩擦阻力,进而使边界层厚度
减小。

-局部流动特性的影响:如流动存在涡旋、偏转或分离现象,边界层
厚度的变化可能会出现非线性的特点。

总结起来,边界层厚度的计算方法主要包括拉格朗日方法和欧拉方法。

拉格朗日方法基于粒子运动的追踪来计算边界层厚度,而欧拉方法则是通
过整体流场分析来获取边界层厚度。

在具体计算中,还应考虑流体属性、
壁面状况和局部流动特性等因素的影响。

边界层分析

边界层分析

流体外掠平板时的层流边界层与紊流边界层
临界距离:由层流边界层开始 xc 向紊流边界层过度的距离:
临界雷诺数: Rec
惯性力 Rec 粘性力 u xc u xc v
平板: Rec 2 105 ~ 3106 ; 取Rec 5 105
紊流边界层: 为何是一个范围? 层流底层:紧靠壁面处,粘滞力会占绝对优势,使粘附于壁的一 层薄层仍然会保持层流特征,具有最大的速度梯度。
流动边界层的厚度与流速、流体的运动粘度 和离平板前缘的距离x的关系
δcm
平板长度l (cm) 空气沿平板流动时边界层增厚的情况
x, 空气速度 u 10m / s :
x100mm 1.8mm; x200mm 2.5mm
由牛顿粘性定律:
u y
速度梯度越大,粘滞应力越大。
流动边界层的几个重要特征: 1. 边界层厚度 与壁的定型尺寸L相比极小
L
2. 边界层内存在较大的速度梯度。
3. 边界层流态分层流和紊流,紊流边界层紧靠壁面处仍有 层流特征,层流底层。
4. 流场可以分为边界层区和主流区。
边界层区:由粘性流体运动微分方程组描述。 主流区:由理想流体运动微分方程—欧拉方程描述。 边界层理论的基本论点: 边界层概念也可以用来分析其他情况下的流动和换热。 如流体在管内的受迫流动,流体外掠圆管流动,流体在竖直
层流靠流体导热换热,紊流 依靠流体微团脉动对流换热

故:紊流换热比层流换热强。
与t 的关系:分别反映流体分子和流体微团的动量和热
扩散的深度
t

Pr
1 3
0.6 Pr 50 层流:
壁面加热流体时热边界层的形成和发展
  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

5
对于管内的流动运 动,取临界雷诺数 2300
粘性底层:在紊流边界层内,由于紧贴壁面处那一层薄层内
粘滞力甚大,流体仍具有层流的特征。 紊流支层:粘性底层上方称为紊流支层,在该层内粘滞力较 小,流体具有紊流的特点。 边界层厚度=粘性底层+紊流支层
底 =29.4 x w
9
9 1 10

m
t∞ u
流体流过固体壁面的流场就 人为地分成两个不同的区域。
δ 0
t
δ
tw x
其一是边界层流动区,这里流体的黏性力与流体的惯性力共 同作用,引起流体速度发生显著变化;其二是势流区,这里 流体黏性力的作用非常微弱,可视为无黏性的理想流体流动, 也就是势流流动。
2)边界层的厚度
当速度变化达到 u u 0.99 时的空间位置为速度边界层的 外边缘,那么从这一点到壁面的距离就是边界层的厚度 x
0
x
x x 5.0 w x
1
2
5.0 Re
1
2
要使边界层的厚度远小于流动方向上的尺度(即 x x 1 ), 也就是所说的边界层是一个薄层,这就要求雷诺数必须足够 Re 1 的大,即
因此,对于流体流过平板,满足边界层假设的条件就是雷 诺数足够大。由此也就知道,当速度很小、黏性很大时或 在平板的前沿,边界层是难以满足薄层性条件。
3) 临界雷诺数
随着x的增大,δ(x)也逐步增大,同时黏性力对流 场的控制作用也逐步减弱,从而使边界层内的流动变得紊乱。 把边界层从层流过渡到紊流的x值称为临界值,记为xc, 其所对应的雷诺数称为临界雷诺数,即 Re c u xc
流体平行流过平板的 临界雷诺数大约是
Re c 5 10
• 过度边界段: Re数处于之间,粘性力和惯性力相当 • 紊流边界段:Re数很大,惯性力主导,忽略粘性力
惯性力
指当物体加速时,惯性会使物体有保持原有运动状态的
倾向,若是以该物体为参照物,看起来就仿佛有一股方向相 反的力作用在该物体上,因此称之为惯性力。因为惯性力实 际上并不存在,实际存在的只有原本将该物体加速的力,因 此惯性力又称为假想力。 惯性系:相对于地球静止或作匀速直线运动的物体. 非惯性系:相对地面惯性系做加速运动的物体.
当Pr>1时,Pr=υ/a,υ>a,粘性扩散 >热量扩散,速度边 界层厚度>温度边界层厚度。 当Pr<1时,Pr=υ/a,υ<a,粘性扩散 <热量扩散,速度边 界层厚度<温度边界层厚度。 也可以从公式得出 T∞ δ δt
1 1 t x x Pr 3 1.026
T∞ u∞
δt
3)速度边界层厚度与热边界层厚度的关系 在同一位置上热边界层厚度与速度边界层厚度的相对大小与 流体的普朗特数Pr有关,也就是与流体的热扩散特性和动量 扩散特性的相对大小有关。
1 1 t x x Pr 3 1.026
由此式可以看出,热边界层是否满足薄层性的条件,除了 Re×足够大之外还取决于普朗特数的大小,当普朗特数非常 小时(Pr<<1),热边界层相对于速度边界层就很厚,反之 则很薄。 普朗特数Pr 的物理意义: 表征流体的热扩散特性和动量扩散 特性的相对大小
对于导热系数不高的流体,紊流换热热阻主要取决于粘性底
层的导热过程,边界层的温度梯度在粘性底层最大,而在紊流 支层变化平缓。
流型:层流和紊流
3.引入边界层概念的意义
缩小计算区域。对对流换热问题的研究可集中在边界 层区域内 边界层内的流动与换热可以利用边界层的特点进一步 简化
二、数量级分析与边界层微分方程
u u p 2u u v ) 2 ( x y x y
dp du 由上式: u dx dx p ~ 0( ) 可视为边界层的又一特性 y
u v 0 x y
u u 1 dp 2u u v 2 x y dx y
§5-4 边界层积分方程组及其求解 一.用边界层积分方程求解对流换热问题的基本思想
(1) (2) (3) 建立边界层积分方程 针对包括固体边界及边界层 外边界在内的有限大小的控制容积; 对边界层内的速度和温度分布作出假设,常用的函 数形式为多项式; 利用边界条件确定速度和温度分布中的常数,然后 将速度分布和温度分布带入积分方程,解出 和 t 的计算式; 根据求得的速度分布和温度分布计算固体边界上的
2. 热(温度)边界层(Thermal boundary layer)
1)定义 当流体流过平板而平板的温度tw与来流流体的温度t∞不相等 时,在壁面上方也能形成温度发生显著变化的薄层,常称为 热边界层。
Tw
2)热边界层厚度 当壁面与流体之间的温差达到壁面与来流流体之间的温差的 0.99倍时,即 (t w t ) /(t w t ) 0.99 ,此位置就是边界层 的外边缘,而该点到壁面之间的距离则是热边界层的厚度, 记为 t x 层流:温度呈抛物线 分布 湍流:温度呈幂函数 分布 湍流边界层贴壁处的 温度梯度明显大于层 流 故:湍流换热比层流换热强!
2t 2t t t c p u v 2 2 x x y y
a.定义5个基本量的数量级 主流速度:u ~ 0(1); 温度:
t ~ 0(1);
壁面特征长度: l ~ 0(1);
边界层厚度:
~ 0( ); t ~ 0( )
u∞ δ x
0
0
x
(a)Pr<1
(b)Pr>1
要点:
热边界层的边界线将流体的温度场划分为两个区域,只有在
热边界层中才有温度变化,而在热边界层以外可以认为温度梯 度为零,当做等温流动区。
对于层流,壁面法线方向热量传递靠导热方式,边界层内温
度分布为抛物线;对于紊流,粘性底层的热量传递靠导热,而 在底层以外的紊流支层,除导热外,主要靠速度脉动引起的对 流混合作用。
1 1 u u p 2u 2u (u v ) ( 2 2 ) x y x x y
u v 0 x y
(a)
(b)
v v p 2v 2v (u v ) ( 2 2 ) x y y x y 2 1 1 ( ) ( 2 ) 2 1 1
§15-3 边界层微分方程组的解 边界层的概念是1904年德国科学 家普朗特提出的。 (Boundary layer) 一. 边界层 1.流动边界层 1)定义 垂直于壁面的方向上流体 流速发生显著变化的流体 薄层定义为流动边界层。
在离平壁前端x处用热线风速仪,测得沿壁
面方向方向上各点的流速,这一分布呈现类似 抛物线型。在u=0.99u∞ 处以外的流体,可以认 为不受流体粘性的影响,称其为主流区。而 u=0.99u∞ 以内的区域,存在明显速度梯度,称 为边界层区。
紊 =0.37 w x 4 5 m
15
紊 x = 0.37 Re1x 5
层 =5.0 xv w
底 =29.4 x w
9
12
9 1 10

m
由以上两式可以发现,流体的主流速度w∞越大,层流边 界层厚度δ层以及粘性底层的厚度δ底越薄;x增大时,层流边
界层厚度δ层随x0.5成正比增加,而粘性底层则随x0.1成正比增
热线风速仪
流体流过固体壁面时,由于壁面层流体分子的不滑移特 性,在流体黏性力的作用下,近壁流体流速在垂直于壁面的 方向上会从壁面处的零速度逐步变化到来流速度。
t∞ u
δ 0
t
δ
tw x
普朗特通过观察发现,对于低黏度的流体,如水和空气 等,在以较大的流速流过固体壁面时,在壁面上流体速度发 生显著变化的流体层是非常薄的。
t t 2t u v a 2 x y y
du dp u dx dx

层流边界层对流 换热微分方程组: 3个方程、3个未 知量:u、v、t,
方程封闭
如果配上相应的 定解条件,则可 以求解
du dp 若 0,则 0 dx dx
三、外掠平板层流换热边界层微分方程式分析解简述
小:空气外掠平板, u=10m/s: x 100mm
理论关系式为 :
1.8mm ; x200mm 2.5mm
t∞ u
层 =5.0 xv w
12
层 x =5.0 w x
1
2
5.0 u x
δ
1 2
t
δ
tw
=5.0 Re 1 2 x
结果为:

x C f , x 0.664 Re 1/ 2 x hx 0.332
5.0 Re 1/ 2 x

x Nu x 0.332 Re 1/ 2 Pr1/ 3 x
Re 1/ 2 Pr1/ 3 x
h 0.332

l
Re 1/ 2 Pr1/ 3
式中,Nux=αx/λ为局部努赛尔数,无量纲,其大小反映了局部 换热的程度。

1

) (
1 1
2
1

2
2

t t t c(u v ) 2 p x y y
p p ~ 0( ) ~ 0(1) y x 表明:边界层内的压力梯度仅沿 x 方向变化,而边界层 内法向的压力梯度极小。
边界层内任一截面压力与 y 无关而等于主流压力

p dp x dx
微分方程:
u v 0 x y
u u 2u (u v ) 2 x y y
2t t t c p u x v y y 2
边界条件:
u | y 0 0, v | y 0 0, t | y 0 t w u | y u , t | y t u const
相关文档
最新文档