山东师大新校区综合教学楼风环境模拟研究.docVIP

PAGE PAGE 1 山东师大新校区综合教学楼风环境模拟研究 　　摘要：基于fluent计算流体力学平台，采用standardk-ε模型simple算法对山东师大新校区综合教学楼周围的风环境进行了数值模拟，结果表明：当计算域尺度从15h×10h×6h调整到15w×10d×6h（其中w为建筑物长度，d为建筑物宽度，h为建筑物高度）时，计算精度明显提高。另外当入流风方向为西北风，来流风速为1.6m时，距建筑物上风向200m位置受建筑物影响西北风变为北风，在教学楼区域形成多湍窝及高风区的复杂风场，在距建筑物下风向500～600m位置，建筑物对风的影响基本消除。 　　关键词：小区风场；数值模拟；计算流体力学；建筑绕流 　　中图分类号：P461.1文献标识码：A文章编号2013　　1引言 　　随着社会、经济的高速发展，城市的规模日益扩大，城市中高大建筑的数量和高度也与日俱增，这些建筑显著改变了城市的风环境。一方面高大密集的建筑群，降低了城市的通风、自净能力，加剧了在低风速条件下城市的空气污染和热岛效应；而另一方面在风速较大时，高大建筑周围会产生局地强风，影响到行人的舒适与安全，引出行人风环境（PedestrianLevelWindEnvironment）问题＼[1，2＼]。对于建筑物周围风场情况的获取通常来自三种方法：现场监测、物理风洞模拟和利用计算机进行数值模拟＼[3＼]。现场监测费时费力，物理风洞花费较高而且我国目前物理风洞实验室较少，而随着计算机技术的不断发展，通过数值模拟了解建筑物周围风场成为趋势。 　　山东师大新校区综合教学楼是新校区的主要教学楼，教学楼风环境跟同学们的日常学习和生活息息相关。建筑主体分为A、B、C、D、E5个区，每个区域通过3层以上的走廊连接，相互连通。教学区为5层建筑，楼高约20m。教学楼处于校园西北部，四周宽阔，周围无明显遮挡物，因此运用数值模拟的方法研究建筑物局地风环境有良好的适用性（图1）。 　　图1教学区楼体示意（东北视角）2模拟方案 　　建筑物周围风场的数值模拟是将建筑物致于流动的风场中，以流动风的质量守恒方程、动量守恒方程、湍流输运方程作为控制方程，采用离散化的数值方法获得风场中离散风压的解，可以确定出风场内建筑物表面的风压力值，同时获得流畅中其他相关的物理量状态值。本文以fluent软件包为平台，对综合教学楼周围风流场进行数值模拟。 　　2.1物理建模 　　教学楼本身体型不规则，四周较开阔，无明显遮挡物。在CAD中按照建筑物原尺寸建立三维物理模型，避免尺寸效应可能对结果产生的影响，并将CAD中的物理模型导入到fluent6.2的前处理模块gambit中。 　　2.2确定计算域 　　位于大气边界层中的建筑物，当风对其绕流时，处于一个完全开口的流动风场中，但风对建筑物作用的影响有一定的范围，在数值模拟时可以给定有限的三维计算区域，并确定计算区域各边界的边界条件，以模拟实际的流动风场。计算流域的确定需要考虑到计算量，计算流域越大，网格越多，计算耗时越长，对计算机的配置要求越高。然而计算区域过小则会导致模拟失真。本文方案一入口距建筑物4～5h（h为建筑物高度）的长度，出口距建筑物9～10h，侧面及顶面距建筑物5h的方案确定计算区域，总范围为450m×280m×120m；方案二入口距建筑物5w（w为建筑物长度），出口距建筑物10w；侧面距建筑物5d（d为建筑物宽度），顶面距建筑物5h确定为计算域，总范围为2400m×880m×120m。 　　2.3网格划分 　　为了避免计算时计算机出现内存不足的错误，将整个计算区域划分为4个区域，分别进行网格划分。网格的类型和数量对计算结果的准确性有很大的影响，非结构网格具有良好的灵活性和适应性，易于网格的自适应，因为建筑物的体型不规则，因此采用非结构网格对其周围流域进行网格划分。方案一得到网格总数1658929，EquiAngleSkew小于0.83，AspectRadio为3.4∶1小于5∶1，网格质量良好以交付fluent计算；方案二得到网格总数2285771，EquiAngleSkew小于0.82，AspectRadio为3.9∶1小于5∶1，网格质量良好可以交付fluent计算。 　　2.4边界条件的选取 　　计算流域的入口边界条件采用velocity-inlet（速度入口边界条件），入口速度1.6m/s，入流风向为西北风；出口采用outflow（完全发展出流边界条件），计算域侧面采用symmetry（对称边界层，定义靠近壁面速度为0），建筑物表面及地面故采用无滑移壁面条件，即壁面处流体速度为0。 　　2.5湍流模型的选择 　　k-ε湍流模型是应用最为广泛的湍流模型，standardk-ε湍流模型采用布西内斯克