1、1环境科学导刊2023,42(4)CN 53-1205/XISSN 1673-9655收稿日期:2022-11-15基金项目:湖南省教育厅项目(21C0768)。作者简介:曾志昊(1999-),男,湖南人文科技学院,能机院19级学生,研究方向:数值计算。通信作者:贺广兴(1980-),男,湖南人文科技学院,博士,副教授,研究方向:城市环境。热浮力对不同高宽比街道峡谷内流场影响研究曾志昊,贺广兴(湖南人文科技学院,能源与机电工程学院,湖南 娄底 417700)摘 要:采用二维街道峡谷数值模型,研究墙体加热、不同高宽比、风速条件下对街道峡谷内部流场的影响。结果表明:街道峡谷的高宽比、来流风速及墙体
2、的受热方式对峡谷内部流场影响显著;受热墙体对峡谷内部流场影响与峡谷高宽比及来流风速大小密切相关;背风面受热可加强内部对流;迎风面受热在低风速时有利于置换峡谷内部空气;底部受热可以在局部范围内加强内部对流;三种边界面的受热方式在来流风速大的时候对流场影响有限;高宽比越大,峡谷内部流场越复杂。关键词:街道峡谷;高宽比;数值模拟;流场分布 中图分类号:X169 文献标志码:A 文章编号:1673-9655(2023)04-0001-050 引言随着我国城市的发展,城市建筑日益增多,城市道路与两侧的建筑形成了不同高宽比(建筑高度/道路宽度)的街道峡谷,峡谷内部的空气流场及汽车排放的尾气的分布情况受到众
3、多研究者的关注。一直以来,风洞实验和数值计算模拟成为研究街道峡谷内部空气流场的主要手段1-2。众多研究中,对峡谷内部流场影响因素的研究集中在街道高宽比3、风速4、建筑布置及结构5、屋顶形状6、人行道上植被7、高架桥8-10等因素。现有文献资料中有少数研究峡谷壁面受热对其内部流场的影响11-12,但是大多集中在某一因素,没有对不同因素的综合影响程度进行深入研究。针对此种情况,本文拟从峡谷的高宽比、来流风速和不同受热墙面三个方面综合分析其对峡谷内部空气流场的影响。1 模型建立1.1 计算区域、网格划分及边界条件1.1.1 模型计算区域研究采用经过验证的标准k-模型13,采用Non-Equilibr
4、ium Wall Functions近壁面边界条件。为了使得结果具有更好的对比性,计算模型与验证模型采用相似的物理模型(图1和图3),区别在于模型的物理尺寸。入口边界高度为H,入口长度为a,街道峡谷宽为b,高度为c,街道峡谷的高宽比R=c/b,依次取具有代表性的R=1、2、3作为研究对象,具体尺寸如表1所示。为了流体在达到峡谷前充分发展,入口边界高度H5b,入口长度a10b。当模型宽度远大于峡谷尺寸时,中间流场受宽度方向上边界的流动影响很小,为节省计算资源,采用二维模拟。通过计算,当风速为0.03 m/s、0.3 m/s、3 m/s,流场Re都大于1100010,满足Re独立性。由于一天内不同
5、时刻受到太阳辐射角度影响,当太阳成角度入射时(图2),接受辐射的面温度会比其他面温度高,因此造成不同面的温度差异。根据文献资料14,模拟温差为20,模拟环境温度取300 K。daHb图1 物理模型图2 太阳不同角度辐射表1 各模型尺寸Ra/mb/mc/md/mH/mR=1905515060R=29010R=39015环境科学导刊第 42 卷第 4 期2023 年 8 月21.1.2 边界条件风速入口边界条件为指数函数15:(1)式中:z距离地面高度;u(z)在z高度位置的水平风速;zref 整个计算区域高度;uref 参考风速。出口边界条件设置为0压力梯度,上边界为对称边界条件,其他边界为无滑
6、移边界。2 模型验证CFD模拟结果的可靠性利用Allegrini的实验进行验证12。街道峡谷模型中的建筑高度和街道宽度分别为0.2 m,入口高度为0.855 m,入口长度为7.4 m,实验槽宽度为1.8 m。由于宽度远大于峡谷宽度,可以简化为二维模型。模型示意图及网格划分情况如图3所示。由图4可知,垂直(y)方向风速和湍动能的模拟结果与实测数据吻合理想,其中风速的吻合度更高,湍动能的吻合有一定偏差,模拟湍动能比实测数据稍小,在两端点处偏大,整体模拟结果可行度在可接受范围内,结果证明了模拟结果可靠。(c)实验测量线(b)网格划分(a)实验模型图3 模型示意及网格划分图3 结果及讨论3.1 R=1
7、时,不同风速条件下街道峡谷内部流场分布3.1.1 v=0.03 m/s由图5可知,在街道峡谷没有受热面时,由于流体的剪力作用在峡谷内部形成一个顺时针漩涡,峡谷内部的流速较低;背风面受热时,在浮力作用下,内部流速变大,漩涡中心和底部两个角位置流速偏小;底部受热时,峡谷内部整体速度变大,这是因为底部热空气在浮力作用下有上浮趋势,峡谷内部空气密度变小,在相同的风力作用下,更有利于流体的流动,从而对峡谷内部流动起到促进作用;迎风面受热下,峡谷内部流场发生了很大变化,背风面的速度较低,但迎风面的浮力克服了流体剪切力,导致空气从迎风面向上流动,流出峡谷。3.1.2 v=0.3 m/s由图6可知,当来流风速
8、为0.3 m/s时,无受热面时,峡谷内部形成一个主要的顺时针漩涡;当背风面受热时,由于空气浮力的促进作用,峡谷内部的空气流速变大,中心和底部两个角位置风速依旧较低;底部受热,内部流场较背风面受热变化不大,此情况与风速为0.03 m/s时区别较大,这是由于当来流风速增大时,空气浮力相对作用变小导致。当迎风面受热时,峡谷内部流场发生了显著变化,形成了两个漩涡,迎风面处形成了逆时针漩涡,背风面上部区域形成顺时针漩涡。此种情况可解释为在此风速条件下,迎风面的浮力克服了向下的流体惯性力,但是还不足以对背风面的流体形成影响,因此在背风面处形成了顺时针漩涡。3.1.3 v=3 m/s由图7可知,来流风速达到
9、3 m/s时,峡谷内部出现了相似流场,在高风速条件下,流体的惯性力远大于由温差引起的浮力,导致边界温差浮力对流场的影响可以忽略不计。3.2 R=2时,不同风速条件下街道峡谷内部流场分布3.2.1 v=0.03 m/s由图8可知,在没有受热面条件下,峡谷上部区域形成一个顺时针漩涡,下部区域在流体剪切力作用下形成一个逆时针漩涡。当背风面受热时,峡谷内部空气流速整体增加,下部区域的逆时针漩涡消失,这是由于背风面附近的空气热浮力作用下,加速了空气向上流动的惯性力。当底对照线横坐标(x)(a)y方向风速0.250.200.150.100.050.00-0.05-0.10-0.15-0.20-0.257.
10、407.457.507.557.60v/vf模拟结果实验数据图4 测量线上模拟结果与实验结果对比(b)湍动能对照线横坐标(x)0.0080.0070.0060.0050.0040.0030.0020.0010.000-0.001k/vf27.407.457.507.557.60模拟结果实验结果热浮力对不同高宽比街道峡谷内流场影响研究曾志昊30.0270.0220.0170.0120.0070.0020.0450.0350.0250.0150.0050.0550.0450.0350.0250.0150.005velocityvelocityvelocityvelocity(a)无受热面(b)背风
11、面受热(c)底部受热(d)迎风面受热图5 峡谷内部流场分布图6 峡谷内部流场分布(b)背风面受热(d)迎风面受热0.220.170.120.070.020.270.220.170.120.070.020.270.220.170.120.070.020.270.220.170.120.070.02velocityvelocity(a)无受热面velocityvelocity图7 峡谷内部流场分布2.21.71.20.70.22.21.71.20.70.22.21.71.20.70.22.21.71.20.70.2velocityvelocityvelocityvelocity(a)无受热面(b)
12、背风面受热(c)底部受热(d)迎风面受热0.0270.0220.0170.0120.0070.002(c)底部受热图8 峡谷内部流场分布0.0270.0220.0170.0120.0070.002(a)无受热面0.0550.0450.0350.0250.0150.005velocityvelocity(b)背风面受热(c)底部受热0.0270.0220.0170.0120.0070.002velocity(d)迎风面受热velocity0.220.170.120.070.02velocity0.270.220.170.120.070.02velocity图9 峡谷内部流场分布0.220.170
13、.120.070.02velocity0.270.220.170.120.070.02velocity(a)无受热面(b)背风面受热(c)底部受热(d)迎风面受热图10 峡谷内部流场分布velocityvelocityvelocityvelocity2.21.71.20.70.22.21.71.20.70.22.21.71.20.70.22.21.71.20.70.2(a)无受热面(b)背风面受热(c)底部受热(d)迎风面受热图11 峡谷内部流场分布velocityvelocityvelocity0.0270.0220.0170.0120.0070.0020.0270.0220.0170.01
14、20.0070.0020.0270.0220.0170.0120.0070.002(a)无受热面(b)背风面受热(c)底部受热(d)迎风面受热0.0750.0650.0550.0450.0350.0250.0150.005velocity0.0270.0220.0170.0120.0070.002环境科学导刊第 42 卷第 4 期2023 年 8 月40.270.220.170.120.070.02图12 峡谷内部流场分布velocityvelocityvelocityvelocity2.21.71.20.70.22.21.71.20.70.22.21.71.20.70.22.21.71.20
15、.70.2图13 峡谷内部流场分布velocityvelocityvelocityvelocity(a)无受热面(a)无受热面(b)背风面受热(b)背风面受热(c)底部受热(c)底部受热(d)迎风面受热(d)迎风面受热部受热时,内部流场出现了左右分边的情况,迎风面边出现了一个顺时针漩涡,背风面底部出现一个逆时针漩涡,此种情况可能是由于流场受底部热空气上浮的影响。迎风面受热时,出现与R=1类似流场,原理同3.1.1节所述。3.2.2 v=0.3 m/s由图5可知,当风速为0.3 m/s时,无受热面、背风面受热和底部受热三种情况峡谷内部的流场类似,内部风速在无受热面条件下最低,背风面受热峡谷上部区
16、域出现较大面积的风速提升,此情况主要受由于背风面附件空气浮力影响,底部受热时,峡谷内部风速较无受热时有一定增大。迎风面受热时,由于浮力克服了空气的惯性力,造成了出现逆时针漩涡,且漩涡外部空气流出了峡谷。3.2.3 v=3 m/s由图10可知,来流风速达到3 m/s时,峡谷内部流场在有无受热面都比较相似,如前3.1.3所述,在高风速下,由于流体惯性力远大于热浮力,受热面附件浮力对流场影响可以忽略不计。3.3 R=3时,不同风速条件下街道峡谷内部流场分布3.3.1 v=0.03 m/s由图11可知,在没有受热面条件下,街道峡谷内部流场与R=2类似,上部区域形成一个顺时针漩涡,下部区域在流体剪切力作
17、用下形成一个逆时针漩涡,不同的是下部漩涡所占空间增大。当背风面受热时,下部区域的逆时针漩涡消失,峡谷内部形成一个大的顺时针漩涡,原理同3.2.1,背风面附近的空气热浮力作用下,加速了空气向上流动的速度。当底部受热时,内部流场出现了峡谷内部两边出现漩涡,迎风面边出现了一个顺时针漩涡,背风面底部出现一个逆时针漩涡,此种情况是由于在风速较小的时候,浮力大于惯性力,流体惯性力不足以抵消热浮力作用,造成在背风面下部形成一个次要的逆时针漩涡。迎风面受热时,在热浮力作用下,峡谷内部空气流出峡谷。3.3.2 v=0.3 m/s由图12可知,在没有受热面条件下,上部区域形成一个顺时针漩涡,下部区域在流体剪切力作
18、用下形成一个逆时针漩涡,与v=0.03 m/s不同的是上部漩涡所占空间增大,下部漩涡变小,这主要是由来流风速增大造成。当背风面受热时,下部区域的逆时针漩涡消失,峡谷内部形成一个大的顺时针漩涡,背风面附近的空气热浮力作用下,加速了空气向上流动的速度。底部受热时,内部流场依然出现两个漩涡,与v=0.03 m/s相比较,迎风侧顺时针漩涡所占空间增大,背风面底部逆时针漩涡变小,此种情况是由于在风速增大时,热浮力小于流体惯性力造成。迎风面受热时,峡谷内部出现上部在流体剪切力作用下形成顺时针漩涡,随着深度的增加,热浮力作用逐渐显现,在峡谷底部形成逆时针漩涡。由此可见,随着墙体高度的增加,由温差引起的热浮力
19、也增大。3.3.3 v=3 m/s由图13可知,来流风速达到3 m/s时,峡谷内部流场在有无受热面都比较相似,如前3.1.3所述,在高风速下,由于流体惯性力远大于热浮力,受热面附件浮力对流场影响可以忽略不计。4 结论(1)街道峡谷内部流场受其高宽比、来流风0.220.170.120.070.020.270.220.170.120.070.020.220.170.120.070.02热浮力对不同高宽比街道峡谷内流场影响研究曾志昊5速、受热墙面的影响显著;(2)峡谷高宽比越大,来流风速越小,内部流场越复杂,反之亦然;(3)背风面墙体受热时可增强峡谷内部对流强度;底部受热使得内部流场变得复杂,且增强
20、了局部的对流强度;(4)迎风面墙体受热时,当来流风速低时,可以提高峡谷内部空气的置换效果;当风速超过一定阈值时,受热墙面对流场影响不显著;(5)峡谷两侧墙体越高,由温差引起的热浮力对流场的影响越大。参考文献:1 Chew L,Glicksman L,Norford L.Buoyant flows in street canyons:Comparison of RANS and LES at reduced and full scalesJ.Build-ing and Environment,2018(146):77-87.2 Cui P.,Li Z,Wang T.Wind-tunnel mea
21、surements for thermal effects on the air flow and pollutant dispersion through different scale ur-ban areasJ.Building and Environment,2016(97):137-151.3 Hang J,Chen X,Chen G,et al.The influence of aspect ratios and wall heating conditions on flow and passive pollutant ex-posure in 2D typical street
22、canyonsJ.Building and Environ-ment,2020(168):106-126.4 李晓莉,战乃岩,吴思洁,等.长春市不同风向结构对街道峡谷污染物扩散的影响J.吉林建筑大学学报,2022,39(2):29-35.5 Ramponi,R,Blocken B.,Janssen W.CFD simulation of outdoor ventilation of generic urban configurations with different urban den-sities and equal and unequal street widthsJ.Building
23、and Environ-ment,2015(92):52-166.6 谢海英,陈康民.屋顶形状对街道峡谷内污染物扩散的影响J.环境科学研究,2007(3):27-327 周姝雯,唐荣莉,张育新,等.街道峡谷绿化带设置对空气流场及污染分布的影响模拟研究J.生态学报,2018,38(17):6348-6357.8 任蘇琪,黄远东,崔鹏义.架空建筑街谷内流动与污染物扩散的数值模拟研究J.上海理工大学学报,2021,43(4):368-377.9 黄远东,侯人玮,崔鹏义.道路路障对街道峡谷内污染物扩散的影响J.环境工程学报,2018,12(4):1135-1147.10 Ding,S,Huang,Y,
24、Cui,P et al.Impact of viaduct on flow reversion and pollutant dispersion in 2D urban street canyon with different roof shapes-Numerical simulation and wind tunnel experimentJ.Science of the Total Environment,2018(671):976-991.11 Uenara K,Murkami S,Oikawa S,et al.Wind tunnel experiments on how therma
25、l stratification affects flow in and above urban street canyonsJ.Atmospheric Environment,2000,34(10):1553-1562.12 Fan Y.,Hunt J,Li Y.Buoyancy and turbulence-driven atmospheric circulation over urban areasJ.Journal of Environmental Sci-ences,2017(59):63-71.13 Wang Le,Pan Q,Zheng X,et al.Effects of lo
26、w boundary walls under dynamic inflow on flow field and pollutant dispersion in an idealized street canyonJ.Atmospheric Pollution Re-search,2016(4):1-12.14 Yang H,Cho K,Charlie L,et al.Numerical investigations of Re-independence and influence of wall heating on flow characteristics and ventilation i
27、n full-scale 2D street canyonsJ.Building and Environment,2020(55):2365-2397.15 Kastner-klein,P,Plate,E.J.Wind-tunnel study of concen-tration fields in street canyonsJ.Atmospheric Environ-ment,1999(33):3973-3979.Study on the Impact of Thermal Buoyancy on the Flow Field in Street Canyons with Differen
28、t Aspect RatiosZENG Zhi-Hao,HE Guang-Xing (Hunan University of Humanities,Science and Technology,School of energy and Mechanical and Electrical Engineering,Loudi Hunan 417000,China)Abstract:Numerical calculation was adopted to study the flow field in street canyon under different aspect ratios and w
29、ind speeds.The results showed that aspect ratios,wind speed and wall heating had significant effect on the flow field inside the canyon.Leeward-wall heating could accelerate internal convection inside the canyon.Windward-wall heating was conducive to displacement ventilation.Ground heating could enh
30、ance internal convection locally.Three heating ways had little impact on the flow field under the high incoming wind speed.The greater the aspect ratios,the more complex the flow field in the canyon.The research results implied a good practical reference value for urban planning.Key words:street canyon;aspect ratios;numerical simulation;flow field