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

《256定风场中庭院建筑资用风压模拟研究.docx》由会员分享，可在线阅读，更多相关《256定风场中庭院建筑资用风压模拟研究.docx（8页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、定风场中庭院建筑资用风压模拟研究*基金项目：国家外国专家局重点项目（DZ20062100070）；建设部科技攻关项目（05-k3-26）.作者简介：李威，男，硕士研究生，主要从事建筑节能和建筑环境保障技术研究.沈阳建筑大学 李 威，冯国会，梁传志，江向阳，宁经洧摘 要：以计算流体力学理论为基础，采用数值的方法定性的分析其迎风面与背风面的总压力之差，并与宽度和高度相同的底层非全架空庭院建筑周围风环境进行了比较分析得到了底层全架空的庭院建筑和底层非全架空的庭院建筑周围速度场、压力场、迹线图的分布情况，并对结果进行了对比、分析同一位置高度上，底层全架空庭院北楼资用风压的均值随着风速增加先增大后减小;

2、底层非全架空庭院北楼资用风压的均值随着风速增加而增大.三种风速下，底层全架空庭院北楼资用风压的均值大体上随着位置高度增加而增大;底层非全架空庭院北楼资用风压的均值随着位置高度的增加而减小;底层全架空北楼的资用风压基本上大于底层非全架空结构北楼，并且同层资用风压的波动较小.该结论可作为底层架空庭院建筑周围风环境研究的前期指导和风洞验证.关键词：底层架空庭院建筑；风场；资用风压；数值模拟多层或高层民用建筑如办公楼、教学楼、住宅楼等，这些建筑物中大多数没有机械通风式的中央空调系统或在每个房间装有各种窗式、分体式的空调器等1.在这些建筑物中，充分利用自然通风是十分重要的，但如何合理、有效地利用自然通风

3、是一个极其复杂的课题，影响自然通风效果的因素很多2.建筑底层架空空间作为城市公共开放空间已受到人们的普遍欢迎和关注3，其研究具有现实意义.以计算流体力学理论为基础，采用数值的方法定性的分析底层全架空庭院建筑(如图1(a)所示)迎风面与背风面的总压力之差，并与宽度和高度相同的底层非全架空庭院建筑(如图1(b)所示)周围风环境进行了比较分析，分析底层架空结构对通风的影响，得到了底层全架空庭院建筑和底层非全架空庭院建筑周围速度场、压力场、迹线图的分布情况，为研究底层架空庭院建筑周围风环境的前期指导和风洞验证提供理论基础.1 数值模拟方法1.1 物理模型以沈阳建筑大学的架空庭院建筑的数值模拟对象.建筑

4、楼体共5层，底层为架空层.整楼高为21.3m，架空底层高为4.3m.(a) 底层全架空 (b) 底层非全架空图1 底层全架空及非全架空模型图2 庭院平面图把庭院模型的四栋围楼根据方位分别称为东楼、南楼、西楼和北楼,如图2所示.两楼之间由一个大厅连接.大厅的尺寸为18.6m18.6m21.3m.即整个庭院模型由东楼、南楼、西楼、北楼及连接两楼的4个大厅组成.底层全架空庭院模型中，东楼、南楼、西楼和北楼都是底层架空的建筑；而底层非全架空庭院模型中，东楼、南楼、西楼和北楼不都是底层架空的，只有南楼和北楼是底层架空的，东楼和西楼是底层非架空的.这两个庭院模型的外观尺寸基本一致（97.2m97.2m21

5、.3m）.对比底层全架空庭院和底层非全架空庭院北楼的周围风场.1.2 计算区域确定根据文献4-8，确定计算区域：上游区域、下游区域分别为建筑物宽度的3倍和10倍，宽度为建筑物宽度的6倍，高度为建筑物高度的4倍.以L、B和H表示计算区域与建筑（）相平行边的长度，则计算区域可表示为L=，B=和H=4h.1.3 控制方程室外大气底层一般属于低速湍流流动范围，并且风速呈典型的梯度分布.常温空气为黏性、不可压缩流体，符合Boussinesq假设.故选取如下控制方程和湍流模型.湍流流动控制方程采用张量形式的雷诺平均Navier-Stokes方程，描述如下：质量守恒 （1）式中：，为无量纲速度和无量纲长度.

6、动量方程 （2）式中：P为无量纲压力；Re为雷诺数；为无量纲空气黏度；为无量纲湍流动能9.湍流动能输运方程 （3）式中：为湍流模型参数10， 取1.0 ；为湍流动能耗散率.湍流动能耗散率输运方程 （4）式中：为湍流模型参数，取1.3；系数分别取1.44，1.0，1.92，1.0 .平均剪切应力产生项，由下式计算 （5）式中：系数分别取0.09，1.0 .上述无量纲量的定义为 11 （6）分别代表有量纲的坐标张量、速度张量、压力、湍流动能及湍流动能耗散率， 为动力黏度.1.4 计算方法采用有限体积法对上述偏微分方程组进行离散.其中，对扩散项采用中心差分方法，而对流项采用QUICK差分方法.采用S

7、IMPLE算法和交错网格技术来实现压力和速度的解耦，以避免不合理压力、速度场的出现12.2 边界条件确定2.1 来流面边界条件平均风速沿高度变化规律按风工程中常用的指数率形式给出： （7）为参考高度处的水平方向分速度，；为粗糙指数，取为0.3013.分析近年沈阳地区气象资料，得到具有代表性的3个对应于参考高度处的值（见表1）表1 风速样本点的选择序号风速(ms-1)备注13.22过渡季节平均风速21.95月平均风速最小值33.93月平均风速最大值2.2 上空面、侧面及出流面边界条件上空面、侧面及出流面边界条件由于计算区域尺寸较实际大许多，故认为出流面处流动已充分发展，各变量的法向梯度为零.上空

8、面、侧面与来流面及出流面一样，均无实际的物理边界，需作边界条件的假设.现计算控制区域取得较大，上空面及侧面离建筑物壁面较远，并且来流为水平风向，可以认为速度沿切线方向的梯度为零，沿法向的值为零，而对于其他标量如等则取为沿法向的梯度为零.2.3 对计算区域内地面及建筑的处理建筑物表面做无滑移处理.将地面处用平均粗糙高度设置为4mm，粗糙系数设置为0.5.建筑物表面平均粗糙高度设置为2mm，粗糙系数设置为0.5.2.4 近壁面边界条件标准模型对充分发展的高速流有效.对于近壁面附近的流动，由于流速较低，雷诺数较小，湍流的发展并不充分.对此，在近壁面处采用壁面函数法来处理.虽然其存在着一定的不足，但其

9、精度足以满足建筑工程上的需求.3 架空庭院建筑对资用风压的影响3.1 数据采集在三种速度场中，采集从二层到五层窗前的总压，根据统计规律求出各层窗前的总压均值，求对应高度的迎风面和背风面的总压压差，即资用风压.3.2 资用风压3.2.1 底层全架空庭院北楼图3(a)、(b)、(c)分别显示了在1.95m/s、3.22m/s、3.93m/s的风场中，底层全架空庭院北楼二层到五层的资用风压（迎风面和背风面的总压压差）随着与中轴面距离的不同而变化.(a) (b) (c)图3 定风速底层全架空庭院北楼由图3可知，风速1.95m/s，底层全架空庭院北楼二层同层资用风压波动最为明显，其次依次是三层、四层和五

10、层，其资用风压的最大值出现在离中轴面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.4

11、99310.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.53737

12、0.18122643.564363.80363.18379对表2的统计数据，分别按定风速和定楼层两种方式来分析.(a)定风速三种风速下，底层全架空庭院北楼资用风压的均值大体上随着位置高度增加而增大；风速1.95m/s,北楼资用风压的标准差大体上随着位置高度增加而减小;风速3.22m/s，北楼资用风压的标准差最大值出现在五层高度，其次依次为三层、四层和二层；风速3.93m/s，北楼资用风压的标准差大体上随着位置高度增加而增大.风速3.22m/s，二层高度位置资用风压的均值为3.49931 Pa，三层为3.60693 Pa，四层为3.83646 Pa，五层为4.45336 Pa.二层到三层资用风压

13、的均值增加了0.10762 Pa，增长3.08%；三层到四层，资用风压增加了0.22953Pa，增长6.36%；四层到五层，资用风压增加了0.6169Pa，增长16.08%.二层到五层，资用风压增加了0.95405 Pa，增长27.26%.在风速为1.95m/s和3.93m/s下，二层到五层的资用风压分别增加了0.05593 Pa和1.09604 Pa，即增加34.59%和44.90%.(b)定楼层同一位置高度上，底层全架空庭院北楼资用风压的均值随着风速增加先增大后减小.当风速从1.95m/s增加到3.22m/s，增大值为1.27m/s，风速增大了65%，北楼二层到五层资用风压均值分别增大3.

14、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，受影响而使资用风压减小最多的是三楼，

15、其次是二楼、四楼和五楼.这里出现的随着来流风速增大而资用风压却减小的现象与南楼的阻挡效应有关.3.2.2 底层非全架空庭院北楼图4(a)、(b)、(c) 显示了不同风速下，底层非全架空庭院北楼资用风压（迎风面和背风面的总压差）随着与中轴面距离的不同而变化. (a) (b) (c)图4 定风速底层非全架空庭院北楼由图4可知，三种风速下，随着与建筑中轴线的距离的增加，底层非全架空庭院北楼各层的资用风压均逐渐增大，但在靠近两端时，二层资用风压突然下降，三层至五层则突然上升.同一水平位置上，底层非全架空庭院北楼的资用风压随着高度的增加而减小.表3 非全架空庭院北楼统计数据位置风速MeanSDSizeM

16、edianMaxMin二层1.950.347340.06453460.316950.465570.28134三层1.950.221340.051460.203370.332170.16606四层1.950.147790.05622460.136710.290080.08561五层1.950.138620.07286460.129350.297680.05426二层3.221.020340.15011460.947241.301810.84737三层3.220.623610.14072460.586540.922550.45506四层3.220.375390.15748460.349870.77

17、7810.20171五层3.220.317540.19717460.285810.769770.09791二层3.931.585840.20451461.50141.963821.32002三层3.930.933170.22465460.889051.388930.67254四层3.930.542650.24223460.500891.158520.28367五层3.930.438560.27751460.386091.104660.13939对表3的统计数据，分别按定风速和定楼层两种方式来分析.三种风速下，底层非全架空庭院北楼资用风压的均值随着位置高度的增加而减小；底层非全架空庭院北楼资用风

18、压的标准差的最大值都出现在五层高度位置.(a)定风速三种风速下，底层非全架空庭院北楼资用风压的均值随着位置高度的增加而减小.底层非全架空庭院北楼资用风压均值的最大值是二层，其次依次是三层、四层和五层.风速3.22m/s下，北楼二层到五层资用风压的均值分别为1.02034Pa、0.62361Pa、0.37539Pa和0.31754Pa.从二层到五层资用风压均值的逐层减小值分别为0.39673Pa、0.24822Pa和0.05785Pa，分别降低了38.88%、39.80%和15.41%.风速1.95m/s下，北楼二层到五层资用风压均值的逐层减小值分别为0.126Pa、0.07355Pa和0.00

19、917Pa，分别降低了36.28%、33.23%和6.20%.风速3.93m/s下，北楼二层到五层资用风压均值的逐层减小值分别为0.65267Pa、0.39052Pa和0.10409Pa，分别降低了41.16%、41.85%和19.18%.风速3.22m/s下，底层非全架空庭院北楼资用风压的标准差的最大值为五层, 其次依次是四层、二层和三层；风速1.95m/s下，底层非全架空庭院北楼资用风压的标准差的最大值为五层, 其次依次是二层、四层和三层；风速3.93m/s下，底层非全架空庭院北楼资用风压的标准差的最大值为五层,其次依次是四层、三层和二层.即五层资用风压的同层波动最明显.(b)定楼层同一位

20、置高度上，底层非全架空庭院北楼资用风压的均值随着风速增加而增大；底层非全架空庭院北楼各层的资用风压的标准差也都随着风速增加而增大.当风速从1.95m/s增加到3.22m/s，风速增长了65%，增大值为1.27m/s，底层非全架空庭院北楼二层到五层资用风压均值分别增大0.673Pa、0.40227Pa、0.2276Pa和0.17892Pa，即分别增长了193.76%、181.74%、154.00%和129.07%.当风速从3.22m/s增加到3.93m/s，风速增长了22%，增大值为0.71m/s，北楼二层到五层资用风压均值分别增大0.5655Pa、0.30956Pa、0.16726Pa和0.1

21、2102Pa，即分别增长了55.42%、49.64%、44.56%和38.11%.即当风速从1.95m/s增加到3.22m/s，最大的受益楼层是二层，其次是三层、四层和五层；当风速从3.22m/s增加到3.93m/s，最大的受益楼层是二层，其次是三层、四层和五层.若把1.95m/s到3.93m/s当成一个大的过程，则其资用风压的均值分别增加1.2385Pa、0.71183Pa、0.39486Pa和0.29994Pa，即分别增长了356.57%、321.60%、267.18%和216.38%，则可知最大的受益楼层是二层，其次是三层、四层和五层.底层非全架空庭院北楼各层的资用风压的标准差也都随着风

22、速增加而增大.当风速从1.95m/s增加到3.22m/s，北楼各层的资用风压的标准差分别增大0.08558Pa、0.08972Pa、0.10126Pa和0.12431Pa，分别增长了132.62%、175.92%、180.11%和170.61%；当风速从3.22m/s增加到3.93m/s，北楼各层的资用风压的标准差分别增大0.0544Pa、0.08393Pa、0.08475Pa和0.08034Pa，分别增长了36.24%、59.64%、53.82%和40.75%.风速1.95m/s增大到3.93m/s，北楼各层的资用风压的标准差分别增大0.13998Pa、0.17365Pa、0.18601Pa

23、和0.20465Pa，分别增长了216.92%、340.49%、330.86%和280.88%.即风速的增长对二层的同层资用风压的波动的影响最小，对三层的同层资用风压的波动的影响最大.3.3 底层全架空庭院北楼与底层非全架空庭院北楼的资用风压对比表4 底层全架空与非全架空庭院北楼对比位置风速MeanSD全架空非全架空差值全架空非全架空差值二层1.950.16170.34734-0.185640.209910.064530.14538三层1.950.208050.22134-0.013290.110540.0510.05954四层1.950.229920.147790.082130.064320

24、.056220.0081五层1.950.217630.138620.079010.048340.07286-0.02452二层3.223.499311.020342.478970.084480.15011-0.06563三层3.223.606930.623612.983320.132860.14072-0.00786四层3.223.836460.375393.461070.123850.15748-0.03363五层3.224.453360.317544.135820.144110.19717-0.05306二层3.932.441331.585840.855490.051620.20451-0

25、.15289三层3.932.432840.933171.499670.089740.22465-0.13491四层3.932.832370.542652.289720.108620.24223-0.13361五层3.933.537370.438563.098810.181220.27751-0.09629由表4的统计数据，风速1.95m/s，底层全架空结构北楼除二层、三层外其余各层的资用风压的均值都大于底层非全架空对应层的资用风压的均值；底层全架空结构北楼除五层外其余各层的资用风压的标准差都大于底层非全架空对应层的资用风压的标准差.风速3.22m/s、3.93m/s，底层全架空结构北楼各层的资

26、用风压的均值都大于底层非全架空对应层的资用风压的均值；底层全架空结构北楼的各层资用风压的标准差都小于底层非全架空对应层的资用风压的标准差.总得来说，三种风速下，底层全架空北楼的资用风压基本上大于底层非全架空结构北楼，并且同层资用风压的波动较小.4 结 论（1）同一位置高度上，底层全架空庭院北楼资用风压的均值随着风速增加先增大后减小.（2）同一位置高度上，底层非全架空庭院北楼资用风压的均值随着风速增加而增大.（3）三种风速下，底层全架空庭院北楼资用风压的均值大体上随着位置高度增加而增大.（4）三种风速下，底层非全架空庭院北楼资用风压的均值随着位置高度的增加而减小.（5）三种风速下，底层全架空北楼

27、的资用风压基本上大于底层非全架空结构北楼，并且同层资用风压的波动较小.参考文献: 1 贾庆贤，赵荣义.吹风对舒适性影响的主观调查与客观评价J.暖通空调，2000，30(3) :15-17. 2 龚光彩，李红祥，李玉国.自然通风的应用与研究J.建筑热能通风空调，2003，22(4):4-6. 3 李艳卓，任斌.谈建筑底层架空式开放空间设计J.辽宁建材，2005，5:60-61. 4 杨洁等.设有空中花园的高层住宅建筑自然通风的研究J.暖通空调. 2004，34(5):1-5. 5 Fluent Inc.，FLUENT Users GuideM. NH: Fluent Inc.，2003. 6 F

28、luent Inc.，FLUENT Users Defined Function ManualM. NH: Fluent Inc.，2003. 7 Fluent Inc.，GAMBIT Modeling GuideM. NH: Fluent Inc.，2003. 8 冯国会等.非对称建筑外形对周围风环境影响的数值模拟J.沈阳建筑大学学报:自然科学版.2006，22(6 ):945-950. 9 Versteeg H K， Malalasekera W. An Introduction to Computational Fluid Dynamics : the Finite Volume MethodM. Essex : Longman Scientific & Technical，1995.10 Launder B E， Spalding D B. Mathematical Models of TurbulenceM. New York : Academic Press， 1972.11 高学平.高等流体力学M.天津:天津大学出版社，2005.12 陶文铨.数值传热学(第2版) M.西安:西安交通大学出版社，2001.13 黄本才.结构抗风分析原理及应用M.上海:同济大学出版社，2001.