90弯管内流动的理论模型及流动特性的数值研究110精.docx
《90弯管内流动的理论模型及流动特性的数值研究110精.docx》由会员分享,可在线阅读,更多相关《90弯管内流动的理论模型及流动特性的数值研究110精.docx(22页珍藏版)》请在冰豆网上搜索。
90弯管内流动的理论模型及流动特性的数值研究110精
第21卷第3期
2004年6月 计算力学学报
ChineseJournalofComputationalMechanics
Vol.21,No.3
June2004
文章编号:
100724708(2004)0320314208
90°弯管内流动的理论模型及流动特性的数值研究
丁 珏, 翁培奋3
(上海大学应用数学和力学研究所,上海200072)
摘 要:
从三维不可压缩雷诺时均Navier2Stokes方程出发,对90°弯曲管道内湍流流动进行数值模拟。
网格划分采用六面体网格,湍流模型为RNGk2Ε模型,在近壁区采用两层壁面模型进行修正,流场的计算结果与实验数据吻合较好。
在此基础上,本文数值研究了来流方向对流场结构和流动特性的影响。
得出在弯管流场中发生了分离现象,且随着来流侧滑角的增大,分离区范围增大。
此外,随着来流从同一侧滑角变换至同一攻角时,横截面的二次流图像中也从具有两个对称主涡变成只具有一个主涡的现象。
关键词:
弯曲管道;Navier2Stokes方程;湍流模型;流动特性中图分类号:
O357.5 文献标识码:
A
1 引 言
弯曲管道广泛地应用于工业、农业等机械设备上,如压缩机、泵及各种类型的热交换器,甚至核动力的管道系统中。
这些形形色色的弯管实现了流体输运和热量、质量交换等功能。
一般情况下,影响流体流动性质的因素有很多,诸如弯管的弯曲程度,流体的来流马赫数Ma,流体运动方向等。
在纵多因素的影响下,弯曲管道内的流场呈现出十分复杂的流动特性。
如一定来流条件下,在管壁附近形成分离区,管道横截面上产生二次流动,这些现象不仅造成流体总压和能量的损失,而且形成的局部障碍区域也使流动系统的阻力增大,降低了热量、质量的交换效率。
因此,弯管内的流动一直受到内流研究者的关注
[123]
较大成绩。
史峰[4]等人利用标准k2Ε湍流模型研究了Rc弯道内的湍流流动。
计算发D=1.69的90°现利用k2Ε模型得到的流场时均速度分布与实验数据相比,在弯段Η=45°截面前,二者符合较好,即计算误差较小;当弯段Η>45°,存在较大的误差。
文献[5]采用湍流大涡模型,对90°弯管进行了数值研究,拓宽了湍流大涡模型在具有强曲率弯曲管道内流体流动问题的应用。
基于以往研究的基础,本文引进Yahhot和Orszag
[7]
新近应用重整化群方法提出的RNGk2Ε
湍流模型,建立90°大曲率弯管内部流场的Navier2Stokes方程,应用交错网格系统下的SIMPLE算
法求解,并将计算结果与文献[1]中的实验数据进行对比。
此外,针对以往研究较少的流体运动方向对流场结构及流动特性影响,本文着重进行了分析。
给出了几种来流侧滑角和攻角下流场内重要参量的分布情况,揭示出流场中出现的物理现象和规律。
。
早期的实验研究以Taylor为代表,他采用激光多谱勒测速仪对90°方形截面弯管内雷诺数分别为790和40000的来流条件进行了系统的研究,给出了层流和湍流两种流态下的时均速度分布,以及弯曲段壁面的压力分布等重要实验结果[1]。
文献[2]应用LDV测量了弯曲度为RcD=1.69(Rc,D分别代表曲率的平均半径与弯管的水动力直径)的矩形截面弯道内的湍流。
而文献[3]则通过热线测速系统对三维90°弯管内流场的气流特性进行了实验研究。
在理论研究方面[426],国内研究者已取得
收稿日期:
2002209229;修改稿收到日期:
20032032311
基金项目:
教育部青年教师教学科研奖励基金;上海市曙光
计划项目及上海市重点学科资助项目1
作者简介:
丁 珏(19732),女,副研究员;
翁培奋3(19652),男,教授,博士生导师1
2 数学模型
2.1 控制方程
通过对连续方程和瞬间Navier2Stokes方程时均化
得到直角坐标系下定常条件、不可压缩流体流动遵循的方程。
连续方程
=0 (i=1,2,3;j=1,2,3)xi
-
(1)
N2S方程
©1995-2006TsinghuaTongfangOpticalDiscCo.,Ltd.Allrightsreserved.
第3期
--
丁 珏,等:
90°弯管内流动的理论模型及流动特性的数值研究
-
31
5
--′(Θuj=-+-Θu′gΦΛΘiuj)-i2
xjxixjxj
-
5=u,v,w,k,Ε
其中S5,#5分别代表通用变量5的源项和有效扩散系数。
J是Jacobi变换矩阵。
u,v和w分别代表曲线坐标下的速度分量。
Φi2=
--
1,0,
i=2i≠2
(2)
其中 ui,p分别代表流体的平均流速和平均压力。
雷诺应力张量-Θuiuj使用由Boussinesgue等提出的涡粘性模型来计算。
2.2 RNGk-Ε湍流模式[7]
RNGk2Ε湍流模型既适应于高雷诺数情况,也
-′-′
3 数值方法
3.1 离散格式
采用有限体积法对控制方程进行数值离散,并采用SIMPLE算法来实施计算。
在计算域内使用交错网格[8]。
离散后方程的一般形式为
Ap5p=
适应于低雷诺数下的湍流流动,即提供了一个微分形式的有效粘性系数表示式,来说明低雷诺数流动效应。
湍动能方程 Θuj
-
∑A
nb
5nb+S5p
其中5表示某个待求物理量,s5p为5参量方程的源项,下标nb代表控制容积p点周围的各网格点。
3.2 边界条件
()
=ΕΑ+Gk+Gb-ΘkΛeff
5xj5xj5xj
(3)
进口条件:
流体均匀地进入管道中,初始运动
速度为1.0ms(雷诺数40000),温度为293K。
出口条件:
压力出口条件。
壁面条件:
固体壁面上采用无滑移条件;温度和压力分别采用绝热壁和等压梯度条件。
湍能耗散率方程
()
Θuj=Α+ΕΛeff
xjxjxj
-
C1Ε
k
(
Gk+C3ΕGb)-
3
C2Ε
Θ
2k
(4)
RNGk2Ε湍流模式只适用于离开固体壁面边
界一定距离的区域,因此对近壁区的湍流采用两层壁面模型的处理技巧,即将求解区域分成:
壁面流动区和核心湍流区。
壁面流动区主要包括粘性子层和部分的完全发展区。
利用湍流Reynolds数Rey=(
y表示由计算域某一点到壁面的法向距Λ
离)来区分这两个区域。
当Rey≥200时,采用上述的湍流模式中的k和Ε方程来求解;当Rey<200时,使用Wolfshtein一方程模型求解近壁区的湍流流动,即仍使用上述的动量方程和湍动能k方程,但是湍流粘性系数为Λt=ΘC对湍流的影响
lΛ=cly[1-lΕ=cly[1-klΛ。
湍流的耗散
其中Α。
k,ΑΕ分别代表k,Ε的反向有效Prandtl数
对于低雷诺数、近壁面流动情况,可以通过消除变量,导出湍流粘性系数的微分方程:
()=1.72d
ΕΛ
3
^Μ-1+CΜ
d^Μ
^其中 Μ=
CΜ=100Λ
3, Γ=SkΕ,3
1+ΒΓ
--() S=(2SijSij)2, Sij=+2xjxi
Γ0=4.38,Β=0.012,C1Ε=1.42,C2Ε=1.68
3C2Ε
=C2Ε+
率可用Ε=k1.5lΕ计算。
长度尺度lΛ,lΕ反映了壁面
A
Λ
为了适应大曲率复杂弯管内对计算网格生成的要求,将上述直角坐标系下的方程变换到任意曲
),通用的输运方程为线坐标系下(Ν,Γ,Φ
++=(5)(5)#+#+5ΝJ5Ν5ΓJ5Γ
exp(-exp(-
)],)]
A
Ε
其中cl=ϑCΛ
-0.75
AΛ=70,AΕ=2cl。
4 算例及结果分析
为了方便分析问题,将弯管分成上游直线段、弯曲段和下游直线段三个部分。
90°弯管的方形截面的
)
(#5+JS5
J©1995-2006TsinghuaTongfangOpticalDiscCo.,Ltd.Allrightsreserved.
316计算力学学报 第21
卷 边长为0.04m。
弯管半径比RcD=2.3,Rc为曲率的平均半径,Rc=0.5(ri+ro),D为水动力直径,本文中取0.04m。
弯曲段内侧壁面的曲率半径ri=0.072m,外侧壁面的曲率半径r0=0.112m。
上下则指均匀来流速度和它在x2z平面投影的夹角。
首先,我们利用上述的模型和数值方法来计算0°侧滑角和0°攻角下的流场中重要物理参量的分布,并与文献[1]中的实验数据比较,来验证模型和算
法的正确性。
图2为计算所得的在不同横截面和不同位置
处,沿主流方向的时均速度分布曲线及其与实验结果的对比。
可以看到,两者符合较好
。
直线段l的长度均为0.3m。
弯管的几何形状如图1所示。
坐标系的原点O位于管道入口截面的中心。
流体运动方向的变化会对流场结构与流动性质产生影响。
文中,流体的侧滑角是均匀来流速度
与其在x-y平面投影的夹角,用Β来表示;攻角Α
©1995-2006TsinghuaTongfangOpticalDiscCo.,Ltd.Allrightsreserved.
第3期
丁 珏,等:
90°
弯管内流动的理论模型及流动特性的数值研究
317
较好。
通过上面几组数据的比较,可以说明本文建立的数学模型能够真实反映流场中流体的运动状况,且计算结果是可靠的。
4.1 不同来流侧滑角下流场的流动特性
现考察来流在0°攻角,侧滑角分别为0°,40°,80°条件下,管道内三维流场出现的分离现象、横截面上的二次流图像及沿主流方向的壁面压力分布情况。
4.1.1 二次流图像
0°侧滑角、0°攻角来流表示均匀流体垂直进入
图中xH表示沿直线段的轴向距离与水动力学直径D之比,D取0.04m。
Vc为体积平均速度。
r表示无量纲的径向位置,rr0)。
z
3
3
3
管道。
图4(a)为该条件下弯管的主流入口截面Η=0°的速度矢量图。
由于此位置受管道曲率的影响很
≡(r-
(ri-r0)
为无量纲的展向位置,z3≡zz12;z12是
小,因此没有出现明显的回流现象,二次流速度也很
小。
图4(b)为弯曲段出口截面Η=90°的速度矢量图。
此时截面上二次流速度较大,已出现一对涡流,且两个主涡的涡心靠近内侧壁面,说明外侧壁面附近的压力大于内侧壁面,推挤流体向内侧壁面流动。
随着流体向弯管下游流去,压力梯度逐渐减小,因此出口截面位置L=0.776m(L表示沿弯管壁面中心线的沿程长度)的二次流图像中一对主涡的涡心远离内侧壁面,接近外侧壁面,如图4(c)所示。
此外,在截面的4个顶角处,还出现了速度值较小、强度较弱的小涡,该涡的放大图如图4(d)所示
。
弯管宽度的12。
图3显示了计算所得的沿弯曲段外侧壁面中心线位置的压力系数cp随极角Η的变化关系。
cp=2(p-pref)(ΘVc),pref是Η=0°,r3=0,z3=0
2
位置处的参考压力值。
观察图3可以看到,由于流
体运动受到弯管曲率和离心力的影响,沿弯曲段外侧壁面的压力值呈现出先增大后减小的变化趋势。
此外,图中计算的cp~Η曲线与实验数据[1]一致性
©1995-2006TsinghuaTongfangOpticalDiscCo.,Ltd.Allrightsreserved.
318
计算力学学报
第21
卷
图5展示了侧滑角为40°情况下沿主流方向上各横截面的速度分布情况。
由于受到来流方向和流场产生分离现象的共同影响,管道上游直线段L=0.06m位置处,横截面上的二次流运动速度已产
“增加”,极易导致流场中小涡消失。
然而,本文计算得出的小涡清晰可见,这也说明了本文采用的数值方法是正确的。
4.1.2 分离现象
0°侧滑角下,沿主流方向的流场分离现象不明
生。
图5(a)显示了在靠近后侧壁面与外侧壁面之间拐角处出现的回流区的放大图。
弯管出口截面的速度矢量如图5(b)所示。
可以看到,图中出现与图4相似的二次流图像,即存在一对大小相近的对称
显。
当侧滑角增大到40°时,情况已发生变化。
在靠近管道入口附近处,分离现象发生,流场中出现一个较大的回流区域,如图7(a)所示,这主要是受到流体粘性和来流方向的影响所致。
此条件下回流区的面积占弯管上游直线段部分的15左右。
而且,在来流的迎风壁面侧,流体速度较大,而背风侧附近即分离现象发生的区域,流体运动速度最小。
此外,在背风侧(外侧壁面)的附面层内,可清晰看到速度分布的特点。
图7(b)是侧滑角为80°时纵剖面(z3=0)上的流体速度分布图。
可以看出,分离现
主涡。
而且,在该截面的4个拐角处,同样也出现了强度较弱的小涡。
当侧滑角进一步增大到80°时,如图6所示,由于二次流造成流体间动量和能量的交换,使得截面上速度分布向外侧壁面和前后侧壁面推移,因此在距离外侧壁面较近的区域又生成一对反向的次涡,如图6(c)所示。
可见,此条件下管道出口截面呈现了多个大小不一的回流区。
通过以上几组图,可以得出:
由于受到弯管曲率及流场逆压梯度的影响,管道截面上产生了二次旋流。
该现象的出现给流场带来了横截面方向上的动能和能量交换,既改变了主流方向的平均速度,又减少了流场的总压。
此外,随着侧滑角的增大,管道出口截面出现了一对主涡、一对次涡,以及两对小涡的图像。
主涡流的强度较大,在管道存在曲率的情况下就会出现;次涡、小涡的强度依次减弱。
而且,随着侧滑角从0°增加到80°,管道出口截面上主涡流的运动速度逐渐减小。
另一方面,由于数值计算中会出现寄生的“数值湍流”,而使湍流粘性
象十分严重。
即相比于40°侧滑角情况,该条件下回流区域的范围明显增大。
在发生分离现象的流场中,逆着主流方向流动的流体在来流介质的冲击、带动下,又重新进入回流区。
这一现象的发生,会使流场产生局部流动障碍,增加了输运系统的阻力。
同时,使流动总能量部分以热量形式耗散掉,降低了管道传输流体的效率。
如果涉及的是热交换器的弯管,就有可能造成管内外热量交换效率降低的后果。
因此,针对本文所讨论的情况,为了实施流动控制,应尽量减小来流速度的攻角,即均匀来流垂直弯管的入口截面进入管道
。
©1995-2006TsinghuaTongfangOpticalDiscCo.,Ltd.Allrightsreserved.
第3期丁 珏,等:
90°
弯管内流动的理论模型及流动特性的数值研究319
©1995-2006TsinghuaTongfangOpticalDiscCo.,Ltd.Allrightsreserved.
320
计算力学学报
第21
卷
4.1.3 内外侧壁面的压力系数分布
图8分别显示了0°,40°和80°侧滑角情况下,沿弯管壁面中心线长度的静压分布。
这里的参考压力等物理量均取弯管入口处来流的参量值(即L=0.0m位置)。
图中可以看到:
对于0°侧滑角的来流
图8(a),弯曲段区域的内侧壁面(弯曲段曲率半径较小的壁面)所受压力值明显低于外侧壁面(曲率半径较大的壁面)的压力值,这主要是流体介质沿弯管流动时受惯性离心力作用所致。
此外,由于流场发生分离现象,加上流动过程中附面层变化引起的损失,使得弯管出口截面内、外侧壁面上的压力小于进口处的压力值。
图8(b)和图8(c)的曲线中,在管道入口附近的内、外侧壁面上均出现压力平台,这主要是流场分离现象发生的区域。
4.2 不同来流攻角下流场的流动特性
上节已对0°侧滑角、0°攻角的情况进行了描述。
下面,着重讨论0°侧滑角下,攻角分别为40°,80°条件下弯管内流体的流动特性。
4.2.1 二次流图像
图9为来流攻角40°条件下弯管出口截面上流体运动速度的矢量图。
可以看到,不同攻角引起截面上二次旋流的图像不同于来流侧滑角的情况。
即二次流图像在管道曲率和逆压梯度的影响下,已不再生成对称的涡流,而是以一个强度较大的固体涡形式存在。
涡流运动方向仍是从外侧壁面流向内侧壁面。
同时,在壁面顶角处均有4个流动方向与固体涡相反的次涡。
当来流攻角增大到80°时,在管道入口附近(
L
=0.1m),来流介质强烈冲击迎风侧壁面(前侧壁
面),造成此时前侧壁面所受的压力值增大,流体运动速度也随之增大。
而且,在迎风侧壁面的阻挡下,流体在靠近两侧壁面的拐角区域急剧转弯,并沿着两侧壁面向背风侧壁面(后侧壁面)流去,如图10(a)所示。
此时,二次流图像中靠近后侧壁面的
拐角处分别出现一个较对称的回流区。
图10(b)显示了弯管出口截面的速度矢量图,仍呈现出单涡现象,但此条件下流体的二次流运动速度小于40°攻角情况。
4.2.2 内外侧壁面的压力系数分布
图11显示了0°侧滑角,40°、80°攻角条件下,沿弯管曲壁面中心线的压力分布曲线。
可以看到,在弯曲段Η=0°~90°区域,曲率半径较大的外侧壁面所受的压力值仍大于内侧壁面。
而且,曲线中出现一段压力平台,这主要是流场发生分离现象所致。
©1995-2006TsinghuaTongfangOpticalDiscCo.,Ltd.Allrightsreserved.
第3期
丁 珏,等:
90°
弯管内流动的理论模型及流动特性的数值研究
321
5 结 论
本文利用数值模拟的方法研究了90°弯曲管道内流体的流动特性。
其中,采用雷诺时均Navier2Stokes方程和RNGk2Ε湍流模型来描述流
.ExperiCheng,GuoRongweimentalinvestigationforseparatedflowinthree2dimensionalbendduct[J].ActaAeronauticaetAstronauticaSinica,1994,15(9):
109521099.(inChinese))
[4] 史 峰,徐 忠,马材芬.大曲率弯道内湍流数值计算
体的湍流流动,同时利用有限体积法和SIMPLE算法实施计算。
计算结果与实验数据吻合得很好,说明建立的数学模型和数值方法是正确的。
此外,文章还着重讨论了以往研究较少的不同来流攻角和来流侧滑角下,流体介质的流动特性。
如横截面上产生的二次旋流、流场中出现的分离现象、迎风及背风壁面上的压力分布等,这些物理图像是符合实际流动的基本特性。
而且,本文的研究结果也为揭示弯管内复杂流动现象的作用机理及实施流动控制提供了参考。
与测量[J].空气动力学学报,1990,8(4):
4232429.
(ShiFeng,XuZhong,MaCaifen.Turbulentflowincurvedductwithlargecurvaturenumericalcomputationandmeasurement[J].
ActaAerod2
ynamicaSinica,1990,8(4):
4232429.(inChinese))
[5] 樊洪明,何钟怡,王小华.弯曲管道内流动的大涡模
拟[J].水动力学研究与进展,2001,16
(1):
78283.
(FanHongming,HeZhongyi,WangXiaohua.Largeeddysimulationofacurvedductflow[J].JournalofHydrodynamics,2001,16
(1):
78283.(in
Chinese))
[6] 许 承,郭荣伟.一种分离流数值求解方法[J].空气
参考文献(References):
[1] TaylorAMKP,WhitelawJH,YianneskisM.
Curved
ductswith
strong
secondarymotion:
动力学学报,1994,12(3):
3632366.(XuCheng,
.AneffectivemethodfornumericalGuoRongwei
simulationofseparatedflows[J].ActaAerodyna2micaSinica,1994,12(3):
3632366.(inChinese))
velocitymeasurementsofdevelopinglaminarandturbulentflow[J].JournalofFluidsEngineering,1982,104:
3502359.
[2] 马材芬,史 峰,徐 忠.90°弯道内紊流的试验研究
[J].工程热物理学报,1990,11
(1):
37240.(MaCaifen,ShiFeng,XuZhong.Developingturbulentflow
in
a
90°curved
duct[J].
Journal
of
EngineeringThermophysics,1990,11
(1):
37240.(in
[7] YakhotV,OrszagSA.Renormalizationgroup
analysisofturbulence1.Basictheory[J].JournalofScientificComputing,1986,1
(1):
39251.
[8] 范维澄,万跃鹏.流动及燃烧的模型与计算[M].
合肥:
中国科学技术大学出版社,1992.(FanWei2
cheng,WanYuepeng.TheModelandComputationofFlowandCombustion[M].Hefei:
Universityof
Chinese))
[3] 许 承,郭荣伟.三维弯管内分离流场的试验研究
[J].
ScienceandTechnologyofChinaPress,1992.(inChinese))
航空学报,1994,15(9):
109521099.(Xu
(下转第329页)
©1995-2006TsinghuaTongfangOpticalDiscCo.,Ltd.Allrightsreserved.
第3期荣见华,等:
一种基于应力的双方向结构拓扑优化算法329
.TheESOalgorithmcatersfortopologyoptiinvolvetheremovalofmaterialmizationbyallowingthe
removalofmaterialfromallpartsofthedesignspace.Withappropriatechequer2boardcontrolsandcontrolsonthenumberofcavitiesformed,
.constraintsthemethodcanreproducetraditionalfullystressedtopologies,andhasbeenappliedintotheproblemswithstaticstress,stiffness,displacementetc.Ifthealgorithmwasrestrictedtotheremovalofsurface2onlymaterial,thenashape
Recentresearch(Q.M.Quern)haspresentedabi2directionaloptimizationproblemissolved.
evolutionarystructuraloptimization(BESO)methodwherebymaterialcanbeaddedtoandremovedoff.ButtherearemuchmoreoscillationstatesinoptimizingiterationprocessesofthisBESOmethod,itleadstolongcalculationtimeforanoptim