基于数值风洞的低矮建筑钝体绕流特征的研究-湖北土木建筑学会.docVIP

基于数值风洞的低矮建筑钝体绕流特征的研究 王鑫，吉柏锋，谷倩 (武汉理工大学土木建筑工程学院，湖北 武汉 430072)1 引 言 [1]统计，联合国公布的资料显示1937~1980年全球死于台风灾害的人数为49.9万人，居十大灾害死亡人数之首。全球每年发生热带风暴和台风80~100个，造成经济损失60~70亿美元和2万人丧生。中国是世界上少数几个受台风影响最严重的国家之一，历史上就是台风重灾区，尤其是近年来台风灾害造成的直接经济损失呈明显上升趋势。进入90年代以来，平均每年损失达100亿元以上。 科研工作者们从上个世纪60年代开始系统地对建筑物风荷载作用开展大规模研究，形成了一门专业的学科，即结构风工程。结构风工程学是综合气象学、空气动力学和气动弹性力学、建筑结构力学、振动工程学、试验力学等多学科的兴交叉学科，其特点是工程应用背景强，理论研究难度大[2]。结构风工程经过几十年的发展，形成了比较完善的体系，研究方法包括理论分析、现场实测、风洞模拟和计算风工程方法。其中计算风工程方法（Computational Wind Engineering，简称CWE）的核心内容是计算流体动力学，通过在计算机上对建筑物周围风流动所遵循的流体动力学方程进行数值求解，并且可借助计算机图形学技术将模拟结果形象地描述出来，以对建筑物周围风场进行仿真模拟。计算风工程方法不受建筑物形状及周围环境的限制，可获得详尽、直观的流场信息，其优点表现为：成本低、速度快、资料完备等[3]。基于数值模拟方法的优点，本文运用计算流体力学知识，采用计算风工程方法对标准低矮建筑低矮表面及周围风流场进行研究。 2 数值模拟理论 图3.1 Realizable 模型（现场实测数据拟合） 由图3.1，将Realizable 模型所得到的测点平均风压系数与现场实测数据[4]以及TJ-2风洞试验数据[5]拟合程度较高，说明模型建立正确，所选择的Realizable模型比较适用。 4 三维低矮建筑物钝流风场绕流研究 4.1 平均风压分布 风对建筑物绕流使得建筑物周围发生分离、再附、回流等复杂现象，导致建筑物表面各个位置的风速不同，进而各个位置的风压值也不相同，因此作用在建筑物上的风压在空间上表现出不均匀性。对于不同体型的建筑物，在相同风速条件下，平均风压在建筑物上的分布规律也会不相同。 对于建筑物表面风压分布，如图4.1所示。在建筑物下部中心位置压力最大，往四周均匀减小。在屋面迎风面的前缘处的负压梯度变化剧烈，最大的负风压系数都出现在屋坡面的边缘以及附近位置处，在几何尺寸具有明显突变处计算负风压力值较大，因而这些部分承受较高的吸力。 （a）迎风面 （b）背风面 图4.1 风向角90°TTU模型表面风压系数等值线图 对于建筑物周围风压分布，如图4.2所示。H=3m处压力分布建筑物周围的流场里，正面压力为正，背面为负，压力变化阶梯沿建筑物表面均匀展开。建筑物正面中心处压力最大，背面则负压力差不多大，建筑物侧边压力最小。建筑物中心切面上的风压分布图可知，建筑物顶面迎风尖点处压力变化剧烈，从最大正压力突变到最大负压力，风流过建筑物，在建筑物上面压力较小，很不规则。 （a）H=3m处压力分布 （b）建筑物中心切面上的风压分布 图4.2 建筑周围风压分布 4.2 建筑周围湍动能分布 对于建筑物周围风压分布，如图4.3所示。由H=3m处的湍动能图可知，在流体钝体绕流过程中，气流分离严重，变化剧烈，伴随有较高的湍动能，但湍动能分布很不规则，钝体绕流产生的湍动能的分布粗具对称性。由建筑物中心切面湍动能图可以看出，在建筑物迎风侧的边角处，阻塞的流体聚集于此，速度加剧，具有较高的湍动能，表明气流分离处的脉动程度较大。 （a）H=3m处的湍动能图 （b）建筑物中心切面湍动能图 图4.3 建筑周围湍动能分布 4.3 建筑物周围风速分布 对于建筑物周围速度分布，如图4.4所示。从H=3m处的速度云图中，看到在建筑物的附近形成了许多的小旋涡还有环绕现象，流体流经建筑物时，经过撞击减速，气流向建筑物两侧分离，速度加剧，绕建筑物产生环绕效果，在建筑物后面产生回流，形成不规则的漩涡。从建筑物中心切面速度云图中，在建筑物顶面迎风尖角处速度变化剧烈，从无限小到非常大，再到无限小，速度变化阶梯从此点向周围发散，流体越过建筑物产生的风速涡流具有不规则性，回流形状也很随机。建筑物对流体产生阻塞作用，流体在接近建筑物上表面后速度陡增，越过后则平缓减速；在建筑物表面速度则陡降，直到为零。 （a）H=3m处的速度云图