收藏本站
收藏 | 手机打开
二维码
手机客户端打开本文

气固两相平面湍流射流的直接数值模拟

夏钧  
【摘要】: 自由剪切流是一种典型的包含转捩和湍流的流动现象,而平面射流是一种比平面混合层复杂的自由剪切流。迄今为止,关于圆射流和矩形射流的研究较多,但平面射流的数值成果,尤其是对湍流研究具有指导意义的DNS数据,仍相对较少。据作者所知,本文是国内较先开展的关于气固两相平面湍流射流的DNS研究。 本文的求解对象为有弱可压缩性的流场控制方程组,并采用单向耦合的方法模拟各种典型尺寸的颗粒在射流流场中的扩散行为和浓度分布规律。DNS的数值方法基于有限差分法。文中采用具有高分辨率的4阶精度的紧致差分格式离散空间导数,同时,时间步进采用低存储量要求的4阶显式Runge-Kutta格式。非物理的出口区(PML缓冲区)方法结合Thompson无反射边界条件被用于边界条件的设置。入口边界采用特征入流边界条件,并在出口边界作压力修正。数值结果揭示出的细致的流动结构证实该套数值方法是切实可行的。 文中首先用DNS方法模拟了二维弱剪切自由射流中随时间发展的流场结构及典型尺寸的颗粒在流场中的运动规律。计算结果得到了以完全对称模式发展的涡量场,同时精确地捕捉到了展向涡的卷起、两涡配对及在弱剪切自由射流中出现的特殊的三涡合并现象。从现象上看,三涡合并是两次两涡配对的线性叠加。两涡配对后的涡层厚度为单个展向涡核的2倍,三涡合并后近似为单个涡核的3倍。在统计结果中,平均纵向速度U与实验值较好地吻合,U在下游x/h=11.5处达到自相似状态。雷诺应力分布有很大的特殊性,即轴线的横向应力被强烈地抑制,而纵向应力却有不断变大的趋势,这是由于特殊的对称涡量场结构引起的。 本文的第二部分为二维气固两相强剪切平面湍流射流的DNS研究。由于两剪切层的对称性在下游x/h=6.0处即被破坏,之后涡街完全按非对称模式发展。典型的拟序大涡结构包括单个的展向大涡、有相同符号的两涡组成的配对涡(与弱剪切射流中类似)及由符号相反的两涡组成的涡对组合。一阶平均量均与实验及前人的计算结果符合得很好。U达到自相似状态的起始点与涡街的对称性被破坏的起始点是几乎重合的,即在x/h=6.0处。平均速度V发展到自相似状态的过程要缓慢得多。本文中,约至x/h=11.5处,可认为V已有自相似性。射流半宽δ_U与流向坐标x的线性关系及轴线速度平方倒数的线性增长率均与实验中得到的线性规律吻合,区别仅在式中的截距值。由于各种实验及计算工况下的不可测因素,存在这种差别是正常的。雷诺应力分布的峰值量级和曲线性态与实验数据是一致的,但由于缺少展向应力R_(ww)~(1/2)对湍流能量的分流,R_(uu)~(1/2)和R_(vv)~(1/2)均积蓄了较多的能量,峰值较大。 在上述两种流场中,相同尺寸的颗粒在其中的扩散行为是类似的。St=0.01的颗粒惯性小,完全受涡结构的控制,因此可视为示踪粒子,再现了下游流场中各种大涡结构。St=0.1的颗粒行为介于St=0.01和1之间,没有鲜明的自身特征。St=1的颗粒表现出有趣的自组织性,即有规律地散步在单个大涡的涡核外围,同时有少量颗粒几乎按直线排列在配对涡涡核的内部,并且在中心两侧,颗粒的瞬时滑移速度大小相等,方向相反,并且颗粒排列所成的直线滞后于实际的两配对涡边界。作者认为这是由于两涡交界处的涡量作用引起的,滞后角由颗粒的动力响应性质决定。St=10和100的颗粒可视为大尺寸颗粒,由于惯性较大,因此受拟序大涡结构的影响较小。 在两种流场中,S片1的颗粒的网格均方根颗粒数Nrms均为最大,这也证实了该种尺寸 的颗粒在流场中的扩散范围最大,浓度分布最不均匀。并且,不同于弱剪切射流的是,在强 剪切射流中,Nrms(St==0.l)的值比其它尺寸颗粒的Nrms值大很多[=2 xNrms(S卜0.01), =3xNrms(St=IO)】,作者认为这是由于所有的颗粒均已参与到涡结构的发展中。而在弱 剪切射流中,始终有部分颗粒在轴线附近以直线行进的方式通过计算区域,影响了Nrms的 统计结果.因此,作者认为强剪切射流中的Nrms值更有参考价值。颗粒的瞬时滑移速度(Vr) 场清晰地揭示出上述各种颗粒在流场中的特殊行为,不失为一种研究颗粒的瞬时扩散性态和 运动趋势的很好的定呈指标。


知网文化
【相似文献】
中国期刊全文数据库 前20条
1 徐大伟,卢平,钟文琪;带压条件下气固两相输送的试验研究[J];锅炉技术;2004年01期
2 李长志,王莹,董鹏;湍流气固两相流动数值模拟理论研究的最新进展[J];电站系统工程;2002年03期
3 杨任刚,张东东,何枫,许宏庆;气固两相自由射流的瞬态流场研究[J];实验流体力学;2005年01期
4 刘向军,徐旭常;稠密气固两相反应流的拉格朗日法[J];北京科技大学学报;2004年02期
5 高铁瑜;章利特;夏庆锋;徐廷相;;超音速气固两相分离试验研究[J];西安交通大学学报;2007年11期
6 林江;李国能;胡桂林;;横流中气固两相射流的数值模拟[J];工程热物理学报;2010年07期
7 黄士铮,严祯荣;竖直上升转水平圆形弯管气固两相流动数值模拟[J];锅炉技术;2004年05期
8 李洁;任兵;陈伟芳;;稀薄流过渡区气固两相喷流的建模与数值模拟[J];空气动力学学报;2005年04期
9 高雪莲,王国玉,刘淑艳,阎为革,潘书杰;旋转自洁式空气滤清器内部气固两相流场的计算与分析[J];北京理工大学学报;2003年06期
10 周昊,孙国俊,郑立刚,岑可法;带有侧边风的气固两相射流混合特性的实验研究[J];中国电机工程学报;2003年08期
11 赵玉晓,李瑞扬,孙斌,秦裕琨;六角切圆锅炉燃烧器区域气固两相流场的试验研究[J];燃烧科学与技术;2004年04期
12 姚文俊;何正浩;邓鹤鸣;;气固两相体放电对电力系统外绝缘的影响[J];中南民族大学学报(自然科学版);2012年02期
13 李瑞扬,赵玉晓,吕薇,王士军,秦裕琨;六角燃烧器出口区域气固两相流场的研究[J];中国电机工程学报;2002年04期
14 邵瑰玮,李劲,杨咏林,安萍,叶齐政,何政浩;静态模拟气固两相体起晕和击穿的实验研究[J];高电压技术;2005年04期
15 周樟华;周昊;李建中;王子兴;岑可法;;撞击式浓淡气固两相射流的扩散特性[J];动力工程;2006年03期
16 刘仲军;徐通模;惠世恩;;等速及固定式取样装置粉尘收集率的数值计算[J];热力发电;1990年06期
17 赵良举;李斌;高虹;李德胜;袁悦祥;曾丹苓;;气固两相介质音速研究[J];工程热物理学报;2007年03期
18 洪江,沈颐身;水平气固两相高浓度输送流动机理[J];化工冶金;1994年03期
19 张东东,许宏庆,何枫;气固两相射流瞬时速度场和浓度场的PIV研究[J];清华大学学报(自然科学版);2003年11期
20 吴浩;张裕中;;针盘式粉碎机内气固两相流场数值研究[J];郑州大学学报(工学版);2008年02期
中国重要会议论文全文数据库 前3条
1 魏耀东;时铭显;;垂直气固两相流动体系中气固两相之间的协同效应探讨[A];中国颗粒学会2004年年会暨海峡两岸颗粒技术研讨会会议文集[C];2004年
2 许宏庆;何文奇;李良杰;张东东;;应用PIV技术对气固两相射流场进行瞬时测量[A];流动显示’2002第五届全国流动显示学术会议论文集[C];2002年
3 冯毓清;许保华;张思军;余艾冰;Zulli P;;不同粒径颗粒在流化床中偏析及混合行为的模拟[A];中国颗粒学会2002年年会暨海峡两岸颗粒技术研讨会会议论文集[C];2002年
中国博士学位论文全文数据库 前8条
1 程易;气固两相流动数值模拟及其非线性动力学分析[D];清华大学;2000年
2 夏钧;气固两相平面湍流射流的直接数值模拟[D];浙江大学;2002年
3 姚军;气固两相圆柱绕流的直接数值模拟和肋条弯管抗磨机理的数值试验研究[D];浙江大学;2002年
4 罗坤;气固两相自由剪切流动的直接数值模拟和实验研究[D];浙江大学;2005年
5 曾卓雄;密相可压气固两相湍流流动的理论分析及其数值模拟[D];西安交通大学;2001年
6 易富兴;气相/煤粉气固两相圆射流燃烧的直接数值模拟[D];浙江大学;2012年
7 赵勇纲;短管道中气相及气固两相平均流速测量方法研究与应用[D];华北电力大学;2014年
8 查旭东;流化床内高浓度气固多相流动和燃烧的数值试验研究[D];浙江大学;2002年
中国硕士学位论文全文数据库 前10条
1 许怀丽;气固两相体放电及其应用研究[D];华中科技大学;2006年
2 周樟华;浓淡气固两相射流扩散规律研究和炉内防结渣技术[D];浙江大学;2006年
3 吴剑波;静电网格法的气固两相流场测量[D];浙江大学;2012年
4 蔡芬;细粉加料装置料腔内气固两相流场的数值模拟分析[D];华中科技大学;2011年
5 褚开维;流化床中气固两相相互作用行为的数值模拟[D];西安建筑科技大学;2001年
6 楚海建;含气固两相夹杂复合材料宏观本构关系的研究[D];扬州大学;2002年
7 刘朋;流化床内稠密气固两相湍流流动特性的研究[D];大连海事大学;2012年
8 刘虎;PDA接收方式和参数对气固两相流场测量影响的研究[D];哈尔滨工业大学;2006年
9 陈玉忠;煤粉自由射流气固两相燃烧的数值模拟[D];重庆大学;2003年
10 兰中秋;电站锅炉尾部烟道湍流气固两相流场的数值模拟[D];重庆大学;2003年
中国知网广告投放
 快捷付款方式  订购知网充值卡  订购热线  帮助中心
  • 400-819-9993
  • 010-62982499
  • 010-62783978