第五章 边界层理论
- 格式:pdf
- 大小:1.14 MB
- 文档页数:5
1Transport Phenomena, Xu Jian, 2009第五章边界层理论边界层概念 边界层方程 边界层分离2Transport Phenomena, Xu Jian, 20095.1 边界层概念在上述层流动量传递的若干实例的分析中,(1)形状简单;(2)引入了假设:管道无限长、忽略进口段影响。
实际问题要复杂得多。
边界层理论,粘滞力对动量传递影响的一般理论,是粘性流体力学的基础,也与热量传递过程和质量传递过程有着密切的关系。
3Transport Phenomena, Xu Jian, 20095.1 边界层概念Prandtl(1904)提出边界层概念,把统一的流场,划分成两个区域,边界层和外流区;其流体流动(沿流动方向和沿与流动方向垂直的方向)有不同的特点。
边界层:流体速度分布明显受到固体壁面影响的区域。
边界层的形成:¾壁面处流体的“不滑脱”no-slip ¾流体的“内摩擦”作用 边界层厚度δ¾U =0Æ0.99 U 04Transport Phenomena, Xu Jian, 20095Transport Phenomena, Xu Jian, 20095.1 边界层概念流过一物体壁面的流体分成两部分¾边界层,粘性流体,不能忽略粘滞力¾外流区,理想流体,可以忽略粘滞力6Transport Phenomena, Xu Jian, 2009边界层理论的要点边界层厚度δ的变化¾前缘处,δ=0¾x ↑, δ↑;沿壁面的法向将有更多的流体被阻滞¾δ<<x边界层内,δ<<x (距离很小);0Æ0.99 U 0(速度变化大)¾速度梯度很大,剪切力很大¾流体速度减慢Æ惯性力<<层外,惯性力与粘性力数量级相当7Transport Phenomena, Xu Jian, 2009边界层流动的转变x<x c (临界距离)层流边界层 过渡区 湍流边界层转变判据:¾临界值:5×105;¾特征长度:距前缘的距离;¾特征速度:来流速度0Re xU ρμ=8Transport Phenomena, Xu Jian, 2009圆管进口段效应靠近管壁部分:边界层,速度减慢;厚度不断增大,进口段长度之后,汇交在管中心处;充分发展段的流动状态取决于交汇处边界层的流动状态;进口段的中心部分:无粘性流动区,速度均匀,区域不断缩小,在边界层汇交时消失;沿程速度不断增大Î压降增大(附加压降);9Transport Phenomena, Xu Jian, 20095.2 边界层方程普兰德边界层方程:量级比较 边界层积分动量方程:动量衡算沿平壁层流边界层的计算:动量积分方程的应用10Transport Phenomena, Xu Jian, 20095.2.1 普兰德边界层方程2222222211x x x x xy y y y y x y u u u u P u u x y x x y u u u u P u u x y y x y μρρμρρ⎛⎞∂∂∂∂∂+=−++⎜⎟∂∂∂∂∂⎝⎠⎛⎞∂∂∂∂∂+=−++⎜⎟⎜⎟∂∂∂∂∂⎝⎠讨论不可压缩流体在平板壁面上的稳态二维层流2222221x x x x Du u u u PDt x x y z υρ⎛⎞∂∂∂∂=−+++⎜⎟∂∂∂∂⎝⎠2222221y y y yDu u u uPDtyx y z υρ⎛⎞∂∂∂∂=−+++⎜⎟⎜⎟∂∂∂∂⎝⎠不可压缩流体的Navier-Stocks 方程不可压缩流体在边界层中作稳态二维流动,方程简化为:y0x u u x y∂∂+=∂∂连续性方程:11Transport Phenomena, Xu Jian, 20095.2.1 普兰德边界层方程普兰德首先发现可以通过比较数量级简化方程:¾Re 较大时,边界层的厚度δ<<x¾边界层内的惯性力和粘性力数量级相当 标准数量级:¾x 为距离的标准数量级,记为x=O(1)¾u 0为速度的标准数量级,记为u 0=O(1)¾边界层厚度δ的数量级记为δ= O(δ),远远小于O(1) 其他物理量的数量级:¾u x 与u 0是一个数量级,记为u x =O(1)¾y 与u 0是一个数量级,记为u x =O(1)12Transport Phenomena, Xu Jian, 2009其他物理量的数量级(1)(1)(1)x x u u O O x x O ∂Δ≈==∂Δ()222(1)(1)(1)(1)x x u u O O x O O x ∂Δ≈==∂Δyx u u x y ∂∂+=∂∂(1)x u O x∂=∂+(1)y u O y∂=∂()y u O δ=(1)1()()x x u u O O y y O δδ∂Δ≈==∂Δ()22222(1)1()()x x u u O O y O y δδ∂Δ≈==∂Δ22221x x x x xy u u u u P u u x y x x y μρρ⎛⎞∂∂∂∂∂+=−++⎜⎟∂∂∂∂∂⎝⎠数量级(1)(1)×1()()δδ×(1)21()δ13Transport Phenomena, Xu Jian, 2009其他物理量的数量级22221x x x x xy u u u u P u u x y x x y μρρ⎛⎞∂∂∂∂∂+=−++⎜⎟∂∂∂∂∂⎝⎠(1)(1)×1()()δδ×(1)21()δInertial Force=Viscous Force:2()O μδρ=1(1)PO xρ∂≤∂22221y yy yx y u u u u P u u x y y xy μρρ⎛⎞∂∂∂∂∂+=−++⎜⎟⎜⎟∂∂∂∂∂⎝⎠(1)()δ×()(1)δ×()δ1()δ2()δ()δ≤()δ(1)14Transport Phenomena, Xu Jian, 2009普兰德边界层方程2210x x xxy yx u u u dP u u x ydx y u u x yμρρ∂∂∂+=−+∂∂∂∂∂+=∂∂000x y x y u u y u u ====∞=时,时,普兰德边界层方程B.C.通过数量级比较得到的简化方程:应用条件:不可压缩流体在边界层中作稳态二维流动,而且Re 比较大15Transport Phenomena, Xu Jian, 20095.2.2 边界层积分动量方程卡门避开使用N-S 方程,直接对边界层进行衡算x 方向质量衡算:¾左侧进入:¾右侧流出:¾上部外流区进入yxz dxdy 1个单位距离δlyu 0, ρμlx u dy ρ∫()00ll x xu dy u dy dxxρρ∂+∂∫∫()lx u dy dxxρ∂∂∫()()2220000000u (-u )ll l l x x x xlx x u dy u dy dx u dy u dy dxx xdx u u dyx ρρρρρ∂∂+−−∂∂∂=∂∫∫∫∫∫x 向净动量变化率:不可压缩流体沿平板壁面的稳态二维流动16Transport Phenomena, Xu Jian, 2009边界层积分动量方程作用于控制体的x 向外力¾壁面剪切力:¾作用在左右侧面的压力差:1s dx τ−⋅⋅1Pdx l x∂−⋅⋅∂00(u )l x x s Pu u dy l x xρτ∂∂−=+∂∂∫0[,]x y l u u δ∈=00(u )x x sP u u dy x xδρδτ∂∂−=+∂∂∫只考虑x 方向的流动00(u )x x s d dPu u dy dx dxδρδτ−=+∫边界层内外压力近似相等00(u )x x sd u u dy dx δρτ−=∫卡门边界层积分动量方程17Transport Phenomena, Xu Jian, 2009边界层积分动量方程可以求出边界层厚度、流体阻力、曳力系数等;方程有u x ,τw ,δ三个变量,需要补充u x =f 1(y),τw =f 2(δ)的关系;需要预先假定一个速度分布方程才能求解,故只能算是一种近似的方法。
第五章边界层理论及层流边界层中的传递现象5.1 边界层理论的要点5.1.1 问题的提出前述,Re∝惯性力/粘性力当Re<1时,惯性力<<粘性力,可用“爬流”模型,略去惯性力项,N-S方程==>爬流方程(stokes近似),解决一些实际问题(沉降、润滑、渗流等),获得比较满意的结果。
但工程流动问题,绝大多数的Re很大。
这时,是否可以完全略去粘性力,使Navier-Stokes方程==>Euler方程(理想流体)。
但是,这样的结果与实际情况相差很大。
突出的一例即“达朗倍尔佯谬(paradox)——在流体中作等速运动的物体不受阻力”。
究竟应当怎样才能正确地处理大Re数的流动呢?这个矛盾一直到1904年,德国流体学家普兰德(Prandtl)提出了著名的边界层理论(大Re数的流动中,大部分区域的惯性力>>粘性力,但在紧靠固——流边界的极薄流层中,惯性力≈粘性力),才令人满意地解决了大Re数的流动的阻力问题。
后人把Prandtl 提出的流动边界层概念,推广到流动系统的传热边界层和传质边界层,从而确定传热、传质的速率以及了解有关的影响因素。
还有人研究了边界层中的化学反应,解决了一些实际问题。
因此,边界层理论被认为是近代流体力学的奠基石。
5.1.2 流动边界层(速度边界层)以平板流动为例,x方向一维稳态流动,在垂直壁面的y方向上,流动可划分为性质不同的两个区域:(1)y<δ(边界层):受壁面影响,法向速度变化急剧,du/dy很大,粘性力大(与惯性同阶),不能忽略。
(2)y>δ(层外主流层):壁面影响很弱,法向速度基本不变,du/dy≈0。
所以可忽略粘性力(即忽略法向动量传递)。
按理想流体处理,Euler方程适用。
这两个区域在边界层的外缘衔接起来,由于层内的流动趋近于外流是渐进的,不是突变的,因此,通常约定:在流动边界层的外缘处(即y=δ处),u x=0.99u∞,δ——流动边界层厚度,δ=δ(x)。
1Transport Phenomena, Xu Jian, 2009第五章边界层理论边界层概念 边界层方程 边界层分离2Transport Phenomena, Xu Jian, 20095.1 边界层概念在上述层流动量传递的若干实例的分析中,(1)形状简单;(2)引入了假设:管道无限长、忽略进口段影响。
实际问题要复杂得多。
边界层理论,粘滞力对动量传递影响的一般理论,是粘性流体力学的基础,也与热量传递过程和质量传递过程有着密切的关系。
3Transport Phenomena, Xu Jian, 20095.1 边界层概念Prandtl(1904)提出边界层概念,把统一的流场,划分成两个区域,边界层和外流区;其流体流动(沿流动方向和沿与流动方向垂直的方向)有不同的特点。
边界层:流体速度分布明显受到固体壁面影响的区域。
边界层的形成:¾壁面处流体的“不滑脱”no-slip ¾流体的“内摩擦”作用 边界层厚度δ¾U =0Æ0.99 U 04Transport Phenomena, Xu Jian, 20095Transport Phenomena, Xu Jian, 20095.1 边界层概念流过一物体壁面的流体分成两部分¾边界层,粘性流体,不能忽略粘滞力¾外流区,理想流体,可以忽略粘滞力6Transport Phenomena, Xu Jian, 2009边界层理论的要点边界层厚度δ的变化¾前缘处,δ=0¾x ↑, δ↑;沿壁面的法向将有更多的流体被阻滞¾δ<<x边界层内,δ<<x (距离很小);0Æ0.99 U 0(速度变化大)¾速度梯度很大,剪切力很大¾流体速度减慢Æ惯性力<<层外,惯性力与粘性力数量级相当7Transport Phenomena, Xu Jian, 2009边界层流动的转变x<x c (临界距离)层流边界层 过渡区 湍流边界层转变判据:¾临界值:5×105;¾特征长度:距前缘的距离;¾特征速度:来流速度0Re xU ρμ=8Transport Phenomena, Xu Jian, 2009圆管进口段效应靠近管壁部分:边界层,速度减慢;厚度不断增大,进口段长度之后,汇交在管中心处;充分发展段的流动状态取决于交汇处边界层的流动状态;进口段的中心部分:无粘性流动区,速度均匀,区域不断缩小,在边界层汇交时消失;沿程速度不断增大Î压降增大(附加压降);9Transport Phenomena, Xu Jian, 20095.2 边界层方程普兰德边界层方程:量级比较 边界层积分动量方程:动量衡算沿平壁层流边界层的计算:动量积分方程的应用10Transport Phenomena, Xu Jian, 20095.2.1 普兰德边界层方程2222222211x x x x xy y y y y x y u u u u P u u x y x x y u u u u P u u x y y x y μρρμρρ⎛⎞∂∂∂∂∂+=−++⎜⎟∂∂∂∂∂⎝⎠⎛⎞∂∂∂∂∂+=−++⎜⎟⎜⎟∂∂∂∂∂⎝⎠讨论不可压缩流体在平板壁面上的稳态二维层流2222221x x x x Du u u u PDt x x y z υρ⎛⎞∂∂∂∂=−+++⎜⎟∂∂∂∂⎝⎠2222221y y y yDu u u uPDtyx y z υρ⎛⎞∂∂∂∂=−+++⎜⎟⎜⎟∂∂∂∂⎝⎠不可压缩流体的Navier-Stocks 方程不可压缩流体在边界层中作稳态二维流动,方程简化为:y0x u u x y∂∂+=∂∂连续性方程:11Transport Phenomena, Xu Jian, 20095.2.1 普兰德边界层方程普兰德首先发现可以通过比较数量级简化方程:¾Re 较大时,边界层的厚度δ<<x¾边界层内的惯性力和粘性力数量级相当 标准数量级:¾x 为距离的标准数量级,记为x=O(1)¾u 0为速度的标准数量级,记为u 0=O(1)¾边界层厚度δ的数量级记为δ= O(δ),远远小于O(1) 其他物理量的数量级:¾u x 与u 0是一个数量级,记为u x =O(1)¾y 与u 0是一个数量级,记为u x =O(1)12Transport Phenomena, Xu Jian, 2009其他物理量的数量级(1)(1)(1)x x u u O O x x O ∂Δ≈==∂Δ()222(1)(1)(1)(1)x x u u O O x O O x ∂Δ≈==∂Δyx u u x y ∂∂+=∂∂(1)x u O x∂=∂+(1)y u O y∂=∂()y u O δ=(1)1()()x x u u O O y y O δδ∂Δ≈==∂Δ()22222(1)1()()x x u u O O y O y δδ∂Δ≈==∂Δ22221x x x x xy u u u u P u u x y x x y μρρ⎛⎞∂∂∂∂∂+=−++⎜⎟∂∂∂∂∂⎝⎠数量级(1)(1)×1()()δδ×(1)21()δ13Transport Phenomena, Xu Jian, 2009其他物理量的数量级22221x x x x xy u u u u P u u x y x x y μρρ⎛⎞∂∂∂∂∂+=−++⎜⎟∂∂∂∂∂⎝⎠(1)(1)×1()()δδ×(1)21()δInertial Force=Viscous Force:2()O μδρ=1(1)PO xρ∂≤∂22221y yy yx y u u u u P u u x y y xy μρρ⎛⎞∂∂∂∂∂+=−++⎜⎟⎜⎟∂∂∂∂∂⎝⎠(1)()δ×()(1)δ×()δ1()δ2()δ()δ≤()δ(1)14Transport Phenomena, Xu Jian, 2009普兰德边界层方程2210x x xxy yx u u u dP u u x ydx y u u x yμρρ∂∂∂+=−+∂∂∂∂∂+=∂∂000x y x y u u y u u ====∞=时,时,普兰德边界层方程B.C.通过数量级比较得到的简化方程:应用条件:不可压缩流体在边界层中作稳态二维流动,而且Re 比较大15Transport Phenomena, Xu Jian, 20095.2.2 边界层积分动量方程卡门避开使用N-S 方程,直接对边界层进行衡算x 方向质量衡算:¾左侧进入:¾右侧流出:¾上部外流区进入yxz dxdy 1个单位距离δlyu 0, ρμlx u dy ρ∫()00ll x xu dy u dy dxxρρ∂+∂∫∫()lx u dy dxxρ∂∂∫()()2220000000u (-u )ll l l x x x xlx x u dy u dy dx u dy u dy dxx xdx u u dyx ρρρρρ∂∂+−−∂∂∂=∂∫∫∫∫∫x 向净动量变化率:不可压缩流体沿平板壁面的稳态二维流动16Transport Phenomena, Xu Jian, 2009边界层积分动量方程作用于控制体的x 向外力¾壁面剪切力:¾作用在左右侧面的压力差:1s dx τ−⋅⋅1Pdx l x∂−⋅⋅∂00(u )l x x s Pu u dy l x xρτ∂∂−=+∂∂∫0[,]x y l u u δ∈=00(u )x x sP u u dy x xδρδτ∂∂−=+∂∂∫只考虑x 方向的流动00(u )x x s d dPu u dy dx dxδρδτ−=+∫边界层内外压力近似相等00(u )x x sd u u dy dx δρτ−=∫卡门边界层积分动量方程17Transport Phenomena, Xu Jian, 2009边界层积分动量方程可以求出边界层厚度、流体阻力、曳力系数等;方程有u x ,τw ,δ三个变量,需要补充u x =f 1(y),τw =f 2(δ)的关系;需要预先假定一个速度分布方程才能求解,故只能算是一种近似的方法。
00(u )x x sd u u dy dx δρτ−=∫18Transport Phenomena, Xu Jian, 20095.2.3 沿平壁层流边界层的计算待定系数法:0nix i i u a y ==∑需要n +1个边界条件以确定n +1个待定系数23x u a by cy dy =+++边界条件:i. y=0, u x =0 (不滑脱)ii. y=0,220d 0x y u dy =⎛⎞=⎜⎟⎝⎠iii. y =δ, u x =u 0iv. y =δ, d 0xudy =3031u 22x u y y δδ⎛⎞=−⎜⎟⎝⎠19Transport Phenomena, Xu Jian, 2009沿平壁层流边界层的计算()300022000000011220031(u )223939280280314012134.641,0,0,10 4.64 4.64Re x x x x s x x sx s y xu d y y u u dy dx u d u u u dy u u dx du u d dxdy u x c x c u x u x δδρτδδδδρτμτμμδδδρμδμδδρρ=−⎛⎞⎛⎞−==−⎜⎟⎜⎟⎝⎠⎝⎠⇒−=⇒=⎛⎞==⇒=⎜⎟⎝⎠⎛⎞⇒=+===⇒==⎜⎟⎝⎠∫∫++处的剪切应力:1220300200.323Re 0.646,,1.2921Re 2sx xL sx D Lu b dx b Lu b L C u bL τρτμρρ−=====∫d d 总曳力F 宽度长度为F 20Transport Phenomena, Xu Jian, 2009【例】沿平壁层流边界层的计算温度为20℃的空气在常压下以5m/s 的速度流过一块宽1 m 的平板壁面。
试计算距平板前缘0.5m 处的边界层厚度及进入边界层内的质量流率,并计算这一段平板壁面的曳力系数与承受的摩擦曳力。
假设临界雷诺数Re xc =5×105。
()53550.5120.530.50000120.51.8110. 1.205/Re 1.664105100.52 4.64Re 0.00569315(3)0.0214/228(4) 1.292.Re 0.x D C Pa s kg m x my y u bdy u bdy u b kg sC δδμρδωρρρδδδ−−−=×==×<×==⎡⎤⎛⎞==−==⎢⎥⎜⎟⎝⎠⎢⎥⎣⎦==∫∫o 解:(1)判断边界层流型:20空气,处的边界层为层流边界层20003170.02382Du C bL Nρ==d F 21Transport Phenomena, Xu Jian, 20095.3 边界层分离边界层分离的必要条件¾流体具有粘性¾流动过程中存在逆压梯度。