孔板后流场的数值模拟研究

1、孔板后流场的数值模拟研究郭婷婷，邹晓辉，王坤玉 (东北电力大学能源与机械工程学院，吉林省 吉林市 132012)摘 要：本文采用基于雷诺时均N-S方程标准k-紊流模型结合SIMPLE算法对夹角150°正六边形孔、180°圆孔和150°圆孔孔板的挡风特性进行了数值研究；对三种板型在Re=2.34×105时的速度场进行了比较，结果表明：在开孔率相近的情况下，各种板型板后附近都会有回流出现，夹角150°、正六边形孔孔板板后的回流区沿主流方向会逐渐消失，速度也将达到一个比较稳定的数值；结果表明：夹角150°正六边形孔孔板的挡风效果最好。关键词

2、：紊流模型；数值模拟；开孔率；回流Numerical simulation on the shelter effect of porous wind fencesGuo Tingting, Zou Xiaohui, Wang Kunyu (Department of Energy and Mechanical Engineering, Northeast Dianli University , Jilin 132012, Jilin Province, China; )ABSTRACT: Numerically simulation was performed to investigate t

3、he shelter characteristic of the shelter effect of porous wind fences which angle ()of 150°and 180°,hole shape of round and hexagon, using control volume method and standard k- model, integrated with SIMPLE algorithm. The velocity field of three fences at the same Re=2.34×105 have bee

4、n compared. The results suggest that as the close porosity, there was re-circulating flow region at downstream behind the fences, along with stream-flow, the re-circulating flow region of the angle ()of 150°and hexagon hole fences would disappear gradually and the velocity trended to a fix valu

5、e. The conclusion was that the angle =150°and hexagon hole fences had better shelter characteristic than others.Key word: turbulence model; numerically simulated; porosity; re-circulating flow引言近二十多年来，日本、美国、澳大利亚、英国、新西兰等国都对孔板后流场进行了研究，其中日本将孔板用于港口露天煤堆场的防尘抑尘工程，其它方面的研究还包括孔板对抑制土壤起尘、飞雪和飞沙也有研究1-3，Bofah

6、和Alhinai4 研究孔板对流沙的防治，Lee5-8等人对POSCO露天煤堆安装的孔板、流体流过孔板后有三角菱形物体的流场进行了数值模拟、风洞监测孔板后沙粒的运动轨迹和实验测量孔板后速度和湍动能场进行了研究，Kim9等人研究了在开孔率相同的情况时孔开孔率相同的情况时孔径对尾迹区流动特性的影响，Park10作者简介：郭婷婷（1973），女，博士，副教授，研究方向：叶轮机械内部流动的空气动力学和计算流体力学。等人研究了孔板对底部有间隙和不同开孔率情况下板后三角菱形物体表面压力的影响。 本文采用基于雷诺时均N-S方程、标准k-紊流模型结合SIMPLE算法对不同工况下的孔板进行数值研究，并与实验结果

7、进行比较。许多的学者在不同的研究方面运用了雷诺时均N-S方程、标准k-紊流模型来研究问题11-14。1物理模型和数值计算方法及边界条件1.1几何模型本文的几何模型如图1所示，其上设有正阻风面（101），侧阻风面（102）；侧阻风面（102）与正阻风面（101）形成角（如图2所示），角为150°；侧阻风面（102）上有通风孔（103），正阻风面（101）上有通风孔（104）。本文几何模型的尺寸为长1400mm，宽40 mm，高875 mm，正阻风面的宽度为80 mm，侧阻风面的宽度为150 mm，正阻风面底面开2排内切圆直径为25 mm的六边形大孔，侧阻风面上开8排内切圆直径8 mm的

8、六边形小孔，开孔率18.79，开孔率关系式如下：×100式中：Ah开孔孔总面积； Ab板的总面积。图1 几何模型示意图图2 夹角Z101102102103104UXY103 1.2.1数值计算方法及边界条件本文采用标准 k-紊流模型，通用控制方程的离散采用控制容积法，对流项差分格式采用二阶迎风，流体压力-速度耦合基于SIMPLE算法。根据几何模型的特点，边界条件包括主气流速度入口边界条件、压力出口边界条件、壁面为无滑移壁面。主气流速度入口边界， ， 压力出口边界给定沿流动方向各参数的一阶导数等于零 Z/H孔103X/H图3网格局部图 孔104 1.2.2网格的划分根据孔板几何模型的特

9、点,使用一套网格计算比较困难,本文将计算区域分为若干个子区域,由于气流流过孔板后扰动加剧,并存在回流现象,所以本文对孔板附近区域进行了网格局部加密,在孔板孔103部分采用楔形网格,在开孔104部分采用三角形网格,在板101未开孔部分采用三角形网格,板102未开孔部分采用结构化网格,板的网格划分如图3所示混合网格数目为263940个。数值计算值实验值u/UY/H图4 数值计算与实验值的比较2计算结果的分析和讨论2.1数值模拟与实验值的比较文献9中对夹角180°，圆孔，开孔率为20，Re=8360的孔板进行了实验研究，本文对相同条件下的孔板做了数值模拟。图4给出了X/H=1时实验和数值模

10、拟的速度分布，主流方向为X轴方向，板高方向为Y轴方向，u是流向速度，U是主流入口平均速度, H为孔板高，从图中可以得到数值计算的结果和实验结果吻合较好，发展趋势基本相同，平均误差为25.7%，结果表明：数值计算的结果是可靠的。2.2三种板型数值计算结果的比较表1给出了三种板的板型及板型参数：表1：孔板类型和板型参数孔板类型夹角孔型开孔率Re文献9中的孔板：板180圆孔202.34×105国内常用孔板：板150圆孔17.04%2.34×105新型孔板：板150正六边型孔18.79%2.34×105图5为三种板型在X/H=1时X轴方向上的速度分布，从图中可以得到板和板

11、在X方向上的速度分布曲线基本一致。在Y/H在00.2之间时，板后u/U范围在00.5之间，板后u/U范围在-1.2 -0.9之间，板的u/U范围在-0.85-0.7之间，速度值都大于板的速度值；在Y/H范围为0.40.9之间时，板的u/U范围在-0.20.2, 板和板的u/U范围在-0.70.5，板、板和板后都会有回流出现，相比之下板后的回流不太明显，结果表明：在X/H=1、X轴方向上板的挡风效果最好。图6是三种板型在X/H=1时Z轴方向上的速度分布，从图中可以得到Y/H的范围在01时，板的v/U范围在-0.00010.00062之间，板的v/U范围在-0.00090.0011之间，板的v/U

12、范围在-0.000670.00062之间，这几种板型都出现了反向速度，但板的v/U正值最小，负值最大，结果表明：在X/H=1、Z轴方向上板的挡风效果最好。图7是三种板型在X/H=1.5时X轴方向上的速度分布，从图中可以得到Y/H在00.2之间时，板的u/U范围在0.250.9之间，板的u/U范围在-1.4-1.2之间，板的u/U范围在-1-0.8之间，当Y/H在0.20.9之间时板的u/U数值达到一个比较稳定的状态，不再变化，为0.25左右，即速度达到一个比较稳定的状态，板和板的u/U范围在-0.850.5之间，有回流现象存在，板在这个区域不再有回流，速度也较小，结果表明：在X/H=1.5、X

13、轴方向上板的挡风效果最好。图8是X/H=1.5时Z方向上的速度的分布图，从图中可以得到，板的v/U范围在-0.000440.00028之间，板的v/U的范围在-0.000280.00022之间，板的v/U范围在-0.000150.00025之间,结果表明：在X/H=1.5、Z轴方向上板的挡风效果最好。比较图5和图7，从图中可以得到，随着孔板后流场的发展，板的Y/H在0.20.9之间时，u/U将趋近一个定值，结果表明：流场趋于稳定；在图5中有回流现象，但在图7中不再回流现象，板和板的u/U的变化趋势不太明显，结果表明：板的回流现象逐渐消失，板和板的回流现象一直存在。比较图6和图8，我们可以得到，

14、图8中板v/U的负值最小值小于图6中负值最小值，图8中板v/U的正值最大值小于图6中正值最大值，板和板v/U的负值和v/U的正值成减小的趋势。3结论(1)夹角180°、圆孔孔板，Re=8360时的数值模拟得出的结果与实验值吻合的较好，结果表明：数值计算的结果是可靠的。(2)夹角150°、正六边形孔孔板风速的变化沿主流方向逐渐趋近于一个稳定值，其它两种板型则会维持一个比较稳定状态，局部风速没有明显的变化。(3)当风通过夹角150°、正六边形孔孔板，会有回流现象出现，沿主流方向，回流逐渐消失。图6 三种板型在Z轴方向上v/U的速度分布v/UY/H150度、圆孔150度

15、、正六边型孔180度、圆孔150度、圆孔150度、正六边型孔180度、圆孔u/UY/Hv/UY/H150度、圆孔150度、正六边型孔180度、圆孔图7 三种板型X/H1.5时的u/U分布图8 X/H1.5时Z方向上v/U的速度分布u/UY/H150度、圆孔150度、正六边型孔180度、圆孔图5 三种板型在X/H=1时X方向上的u/U分布(4)在开孔率相近的情况下，夹角150°、正六边形孔孔板的挡风效果好于其它两种板型的挡风效果。参考文献：1 D.A. Gillette. On the production of soil wind erosion aerosols having th

16、e potential for long range transportJ. Res Atmos, 1974, 8(3/4):735774.2 R.J. Kind. Mechanics of aeolian transport of snow and sand J. Wind Eng. Ind. Aerodyn, 1990, 36:855866.3 J.A. Peterka, R.L. Petersen. On the relaxation of saltation length as a modeling criteria for particulate transport by windJ

17、. Wind Eng. Ind. Aerodyn,.1990, 36:867876.4 K.K. Bofah, K.G. Alhinai. Field tests of porous fences in the regime of sand laden windJ. Wind Eng.Ind. Aerodyn,1986,23:309319. 5 S. J. Lee , C. W. Park. The shelter effect of porous wind fences on coal piles in POSCO open storage yard . J .Wind Eng. Ind.

18、Aerodyn, 2000,84:101-118.6 S. J. Leea, K. C. Parkb, C. W. Park. Wind tunnel observations about the shelter effect of porous fences on the sand particle movements. Atmospheric Environment, 2002,36: 1453-1463.7 S. J. Leea, H. C. Limb. A numerical study on flow around a triangular prism located behind

19、a porous fence. Fluid Dynamics Research, 2001, 28:209-221. 8 S. J. Leea, H. B. Kim. Laboratory measurements of velocity and turbulence field behind porous fences. J .Wind Eng. Ind. Aerodyn , 1999,80:311-326.9 H. B. Kim, S. J. Lee. Hole diameter effect on flow characteristics of wake behind porous fe

20、nces having the same porosity. Fluid Dynamics Research, 2001, 28:449-464.10 C. W. Park, S. J. Lee. The effects of a bottom gap and non-uniform porosity in a wind fence on the surface pressure of a triangular prism located behind the fence. Journal of Wind Engineering and Industrial Aerodynamics, 200

21、1, 89:1137-1154. 11 李彦鹏,张强,顾璠,徐益谦.二次风反切对大容量切向燃烧锅炉烟气偏差影响的数值研究.中国电机工程学报J, 2001,21(9):33-37. LI Yan-peng,ZHANG Qiang,GU Fan,XU Yi-qian. Numerical study on the effect of reverse swirl of secondary air on flue gas imbalance in large tangentially fired boilerJ. Proceedings of the CSEE, 2001, 21(9):33-37(i

22、n Chinese). 12 李军,邓清华,丰镇平.蜂窝汽封和迷宫式汽封流动性能比较的数值研究.中国电机工程学报J, 2005, 25(16):108-112. LI Jun, DENG Qing-hua, FENG Zhen-ping. Comparison of the flow characteristics for the honeycomb and labyrinth seal using numerical simulation J. Proceedings of the CSEE, 2005,25(16):108-112(in Chinese). 13 赵鸣志,浦兴国,周洁,袁镇福.高温两相流风洞收缩段和试验段的数值模拟研究.中国电机工程学报J, 2006, 26(9