256定风场中庭院建筑资用风压模拟研究

定风场中庭院建筑资用风压模拟研究沈阳建筑大学李威，冯国会，梁传志，江向阳，宁经洧摘要：以计算流体力学理论为基础，采用数值的方法定性的分析其迎风面与背风面的总压力之差，并与宽度和高度相同的底层非全架空庭院建筑周围风环境进行了比较分析．得到了底层全架空的庭院建筑和底层非全架空的庭院建筑周围速度场、压力场、迹线图的分布情况，并对结果进行了对比、分析．同一位置高度上，底层全架空庭院北楼资用风压的均值随着风速增加先增大后减小;底层非全架空庭院北楼资用风压的均值随着风速增加而增大.三种风速下，底层全架空庭院北楼资用风压的均值大体上随着位置高度增加而增大;底层非全架空庭院北楼资用风压的均值随着位置高度的增加而减小;底层全架空北楼的资用风压基本上大于底层非全架空结构北楼，并且同层资用风压的波动较小.该结论可作为底层架空庭院建筑周围风环境研究的前期指导和风洞验证.关键词：底层架空庭院建筑；风场；资用风压；数值模拟多层或高层民用建筑如办公楼、教学楼、住宅楼等，这些建筑物中大多数没有机械通风式的中央空调系统或在每个房间装有各种窗式、分体式的空调器等[1].在这些建筑物中，充分利用自然通风是十分重要的，但如何合理、有效地利用自然通风是一个极其复杂的课题，影响自然通风效果的因素很多[2].建筑底层架空空间作为城市公共开放空间已受到人们的普遍欢迎和关注[3]，其研究具有现实意义.以计算流体力学理论为基础，采用数值的方法定性的分析底层全架空庭院建筑(如图1(a)所示)迎风面与背风面的总压力之差，并与宽度和高度相同的底层非全架空庭院建筑(如图1(b)所示)周围风环境进行了比较分析，分析底层架空结构对通风的影响，得到了底层全架空庭院建筑和底层非全架空庭院建筑周围速度场、压力场、迹线图的分布情况，为研究底层架空庭院建筑周围风环境的前期指导和风洞验证提供理论基础.1数值模拟方法1.1物理模型以沈阳建筑大学的架空庭院建筑的数值模拟对象.建筑楼体共5层，底层为架空层.整楼高为21.3m，架空底层高为4.3m.(a)底层全架空(b)底层非全架空图1底层全架空及非全架空模型基金项目：国家外国专家局重点项目（DZ20062100070）；建设部科技攻关项目（05-k3-26）.作者简介：李威，男，硕士研究生，主要从事建筑节能和建筑环境保障技术研究.图2庭院平面图把庭院模型的四栋围楼根据方位分别称为东楼、南楼、西楼和北楼,如图2所示.两楼之间由一个大厅连接.大厅的尺寸为18.6m×18.6m×21.3m.即整个庭院模型由东楼、南楼、西楼、北楼及连接两楼的4个大厅组成.底层全架空庭院模型中，东楼、南楼、西楼和北楼都是底层架空的建筑；而底层非全架空庭院模型中，东楼、南楼、西楼和北楼不都是底层架空的，只有南楼和北楼是底层架空的，东楼和西楼是底层非架空的.这两个庭院模型的外观尺寸基本一致（97.2m×97.2m×21.3m）.对比底层全架空庭院和底层非全架空庭院北楼的周围风场.1.2计算区域确定根据文献[4-8]，确定计算区域：上游区域、下游区域分别为建筑物宽度的3倍和10倍，宽度为建筑物宽度的6倍，高度为建筑物高度的4倍.以L、B和H表示计算区域与建筑（hbl）相平行边的长度，则计算区域可表示为L=lll5.25.2，B=lll103和H=4h.1.3控制方程室外大气底层一般属于低速湍流流动范围，并且风速呈典型的梯度分布.常温空气为黏性、不可压缩流体，符合Boussinesq假设.故选取如下控制方程和湍流模型.湍流流动控制方程采用张量形式的雷诺平均Navier-Stokes方程，描述如下：质量守恒0iiXU（1）式中：iU，iX为无量纲速度和无量纲长度.动量方程]32)Re1[()(kXUXXPXUUiitiiiji（2）式中：P为无量纲压力；Re为雷诺数；t为无量纲空气黏度；k为无量纲湍流动能[9].湍流动能输运方程PXkXXkUiktiii])Re1[()(（3）式中：k为湍流模型参数[10]，取1.0；为湍流动能耗散率.湍流动能耗散率输运方程kfcPkfcXXXUitiii22211])Re1[()(（4）式中：为湍流模型参数，取1.3；系数2211,,,fcfc分别取1.44，1.0，1.92，1.0.平均剪切应力产生项][ijjijitXUXUXU，t由下式计算2kfct（5）式中：系数fc,分别取0.09，1.0.上述无量纲量的定义为[11]0020202000Re,,,,,LUUUkkUpPUuULxXiiii（6）,,,,kpuxii分别代表有量纲的坐标张量、速度张量、压力、湍流动能及湍流动能耗散率，为动力黏度.1.4计算方法采用有限体积法对上述偏微分方程组进行离散.其中，对扩散项采用中心差分方法，而对流项采用QUICK差分方法.采用SIMPLE算法和交错网格技术来实现压力和速度的解耦，以避免不合理压力、速度场的出现[12].2边界条件确定2.1来流面边界条件平均风速沿高度变化规律按风工程中常用的指数率形式给出：)(bbZZUU（7）bU为参考高度bZ处的水平方向分速度，mZb10；为粗糙指数，取为0.30[13].分析近年沈阳地区气象资料，得到具有代表性的3个对应于参考高度bZ处的值bU（见表1）表1风速样本点的选择序号风速(m·s-1)备注13.22过渡季节平均风速21.95月平均风速最小值33.93月平均风速最大值2.2上空面、侧面及出流面边界条件上空面、侧面及出流面边界条件由于计算区域尺寸较实际大许多，故认为出流面处流动已充分发展，各变量的法向梯度为零.上空面、侧面与来流面及出流面一样，均无实际的物理边界，需作边界条件的假设.现计算控制区域取得较大，上空面及侧面离建筑物壁面较远，并且来流为水平风向，可以认为速度沿切线方向的梯度为零，沿法向的值为零，而对于其他标量如,k等则取为沿法向的梯度为零.2.3对计算区域内地面及建筑的处理建筑物表面做无滑移处理.将地面处用平均粗糙高度设置为4mm，粗糙系数设置为0.5.建筑物表面平均粗糙高度设置为2mm，粗糙系数设置为0.5.2.4近壁面边界条件标准k模型对充分发展的高速流有效.对于近壁面附近的流动，由于流速较低，雷诺数较小，湍流的发展并不充分.对此，在近壁面处采用壁面函数法来处理.虽然其存在着一定的不足，但其精度足以满足建筑工程上的需求.3架空庭院建筑对资用风压的影响3.1数据采集在三种速度场中，采集从二层到五层窗前的总压，根据统计规律求出各层窗前的总压均值，求对应高度的迎风面和背风面的总压压差，即资用风压.3.2资用风压3.2.1底层全架空庭院北楼图3(a)、(b)、(c)分别显示了在1.95m/s、3.22m/s、3.93m/s的风场中，底层全架空庭院北楼二层到五层的资用风压（迎风面和背风面的总压压差）随着与中轴面距离的不同而变化.051015202530-0.20.00.20.40.6v=1.95m/s资用压差/Pa距离/m二层三层四层五层0510152025303.43.63.84.04.24.44.6v=3.22m/s资用压差/Pa距离/m二层三层四层五层0510152025302.22.42.62.83.03.23.43.63.8v=3.93m/s资用压差/Pa距离/m二层三层四层五层(a)(b)(c)图3定风速底层全架空庭院北楼由图3可知，风速1.95m/s，底层全架空庭院北楼二层同层资用风压波动最为明显，其次依次是三层、四层和五层，其资用风压的最大值出现在离中轴面16m左右的位置.风速3.22m/s、3.93m/s，基本上随着高度的增加资用风压增大，并且北楼两侧资用风压略大于中间位置.表2全架空庭院北楼统计数据位置风速MeanSDSizeMedianMaxMin二层1.950.16170.20991640.133790.49718-0.18299三层1.950.208050.11054640.191580.381090.06486四层1.950.229920.06432640.230570.322830.12211五层1.950.217630.04834640.232380.279310.1228二层3.223.499310.08448643.505193.796873.37897三层3.223.606930.13286643.598923.944713.43092四层3.223.836460.12385643.837254.082023.66427五层3.224.453360.14411644.46484.685624.24823二层3.932.441330.05162642.43542.670052.3699三层3.932.432840.08974642.432742.758532.30499四层3.932.832370.10862642.832423.049032.646五层3.933.537370.18122643.564363.80363.18379对表2的统计数据，分别按定风速和定楼层两种方式来分析.(a)定风速三种风速下，底层全架空庭院北楼资用风压的均值大体上随着位置高度增加而增大；风速1.95m/s,北楼资用风压的标准差大体上随着位置高度增加而减小;风速3.22m/s，北楼资用风压的标准差最大值出现在五层高度，其次依次为三层、四层和二层；风速3.93m/s，北楼资用风压的标准差大体上随着位置高度增加而增大.风速3.22m/s，二层高度位置资用风压的均值为3.49931Pa，三层为3.60693Pa，四层为3.83646Pa，五层为4.45336Pa.二层到三层资用风压的均值增加了0.10762Pa，增长3.08%；三层到四层，资用风压增加了0.22953Pa，增长6.36%；四层到五层，资用风压增加了0.6169Pa，增长16.08%.二层到五层，资用风压增加了0.95405Pa，增长27.26%.在风速为1.95m/s和3.93m/s下，二层到五层的资用风压分别增加了0.05593Pa和1.09604Pa，即增加34.59%和44.90%.(b)定楼层同一位置高度上，底层全架空庭院北楼资用风压的均值随着风速增加先增大后减小.当风速从1.95m/s增加到3.22m/s，增大值为1.27m/s，风速增大了65%，北楼二层到五层资用风压均值分别增大3.33761Pa、3.39888Pa、3.60654Pa和4.23573Pa，即分别增长了2064.08%、1633.68%、1568.61%和1946.30%.当风速从3.22m/s增加到3.93m/s，风速增长了22%，增大值为0.71m/s，北楼二层到五层资用风压均值分别减小1.05798Pa、1.17409Pa、1.00409Pa和0.91599Pa，即分别降低了30.23%、32.55%、26.17%和20.57%.即当风速从1.95m/s增加到3.22m/s，最大的受益楼层是二楼，其次是五楼、三楼和四楼；当风速从3.22m/s增加到3.93m/s，受影响而使资用风压减小最多的是三楼，其次是二楼、四楼和五楼.这里出现的随着来流风速增大而资用风压却减小的现象与南楼的阻挡效应有关.3.2.2底层非全架空庭院北楼图4(a)、(b)、(c)显示了不同风速下，底层非全架空庭院北楼资用风压（迎风面和背风面的总压差）随着与中轴面距离的不同而变化.0510152025300.10.20.30.40.5v=1.95m/s资用压差/Pa距离/m二层三层四层五层0510152025300.00.20.40.60.81.01.21.4v=3.22m/s资用压差/Pa距离/m二层三层四层五层0510152025300.00.51.01.52.02.5v=3.93m/s资用压差/Pa距离/m二层三层四层五层(a)(b)(c)图4定风速底层非全架空