关于建筑物表面风压的数值模拟
本人初学FLUENT,现在需要模拟建筑物表面风压与风洞实验的结果相对照，希望各位大虾能够提供一、二个完整的算例，谢谢。（时间太紧，我们这边无人做数值仿真，太着急了）回复：（cpcx）菜鸟紧急求助：关于建筑物表面风压的数...
严寒地区住宅小区室外风场的数值模拟分析1 引言随着我国严寒地区低能耗住宅建筑的发展，住宅室内通风换气问题已不容忽视。一般情况下，室内自然通风的形成，既有热压通风的因素，也有风压通风的原因，从自然通风改善室内空气品质角度来看，风压通风对室内气候条件的效果比较显著，故应首先考虑如何组织建筑物室外的风压通风来改善室内热环境。
2 室外风场的物理模型和CFD数值模拟2.1 物理模型
哈尔滨市位于严寒地区，冬季持续时间长，且室内空气质量与室外环境相差较大，故节能住宅建筑的通风关键在冬季，本文以哈尔滨地区气象参数中冬季的主导风向和风速为依据，以哈尔滨市泰海小区44号楼及其周围4栋建筑物作为室外风场模拟对象，分析住宅小区室外风场的气流流动情况。模拟建筑物及其周围四栋楼均为高度为22m的建筑物，如图1，图中相应地给出各建筑物在泰海小区中的位置及其建筑物布局。
为建立数学模型，对物理模型作以下假设和简化：
（1）建筑物外气流分布取决于来风风速以及风向，建筑尺寸及形状，以及建筑物开口大小和位置。若开口尺寸小于建筑物立面面积的1/6，三栋建筑可简化为混凝土块。
（2）室外气流为风速梯度分布的低速流范围，据Boussinesq假设，空气一般为粘性不可压缩流体。一次简化为稳态的紊流气流流动，考虑到计算机的硬件设备(RAM256M，CPUPⅣ2.4GHz)有限，仅分析最大风速的稳态紊流情况。
2.2 CFD数值模拟
FLUENT软件设计基于“CFD计算机软件群的概念”，针对每一种流动的物理问题的特点，采用适合于它的数值解法，从而高效率的解决各个领域的复杂流动的计算问题。
FLUENT中提供了下列可供选取的湍流模型：Spalart-Allmaras模型、标准k-ε模型、RNG（重组化群）k-ε模型、可实现k-ε模型、雷诺应力模型（RSM）和大涡模拟模型（LES）。湍流模型选取取决于诸多因素，如流动物理机理、特定类型问题以往的经验、精度级别的要求、现有的计算机资源和模拟所用时间等。对于住宅小区这样具有较大的建筑物尺寸和较高的风速的特定条件，室外流动的Re从50.000到100.000变化，为完全发展流动，因此，采用标准k-ε湍流模型。
参见前人对计算模拟区域的经验设定，室外流动模型模拟区域如下：当所着重模拟的建筑物外表尺寸为1时，模拟区域为上风侧为建筑物长度的3倍，下风侧为建筑物长度的12倍，两侧宽度为建筑物的3倍，高度为建筑物高度的4倍。几何建模和网格划分采用FLUENT的前置处理器－GAMBIT。
建筑物室外风场的来流为哈尔滨地区冬季主导风：风向西南，平均风速为按10米高处风速3.8http://www.shuigong.com/d/file/papers/dianqi/2005-11-16/987ffeb3edd3ca9b72de36e2594bb290.gif计算的沿高度递增的梯度风速。上空面、地面及建筑物表面按光滑壁面设定。方程求解中压力与速度的耦合采用压力耦合的半隐方法（SIMPLE），除压力采用二阶迎风格式进行离散外，其他如动量、紊流脉动动能和紊流脉动动能耗散率均采用一阶迎风格式进行离散。
http://forum.vibunion.com/neteaseivp/resource/paper/doc/20059291127958820080/image002.jpghttp://www.shuigong.com/d/file/papers/dianqi/2005-11-16/fe360e540f3f34eba74426b683fb8383.jpg
图1 哈尔滨泰海小区44号楼及其周围4栋建筑物平面图3 结果分析与讨论3.1 室外风速矢量场分析
为了研究建筑物周围不同朝向不同高度处的室外气流流动情况，分别计算了位于44号楼中的两个算例：（1）平面高度5.94m（以地面为基准的送风高度）；（2）平面高度19.94m。
由图2中的速度矢量分布来看，在西南风向的影响下，建筑物群的西南向建筑物处于迎风侧，而东北向建筑处于背风侧。在建筑物群外侧的西北角和东南角以及建筑物群的入口处，速度梯度达到最大值；并在建筑物群背风侧的西北角和东南角产生背风涡流区。
建筑物群外侧，速度沿南向建筑物的变化规律为：由西向东逐渐增大，在建筑物的拐角处达到最大值；速度沿西向建筑变化规律为：由南向北逐渐增大，在建筑物的拐角处达到最大值。沿西南向建筑物的速度绝对值较大，速度方向变化不大。
在建筑物群外侧，速度沿北向建筑物的变化规律为：40号楼侧，速度由西向东速度先变小后变大，在建筑物拐角处均达到最大值，速度方向发生180°变化；42号楼侧，速度大小始终由西向东增大，且速度大小和方向变化较平缓。速度沿东向建筑物的变化为：由北向南速度大小稍有增加，速度方向基本不变。
http://forum.vibunion.com/neteaseivp/resource/paper/doc/20059291127958820080/image004.jpghttp://www.shuigong.com/d/file/papers/dianqi/2005-11-16/6736ac22f9c85a4eaae897f94aad1e68.jpg
算例1：平面高度5.94m 算例2：平面高度19.94m
图2 44号楼及其周围4栋建筑物室外风场速度矢量图及等动压线图
注：图中网格为■的位置分别是44号楼3单元202、702户南、北向房间位置。在建筑群内，速度大小变化较小，但方向沿围护结构变化很大。因此，在左右两个马蹄形建筑群内形成了两个强度相似，但旋转方向相反的旋涡。
由此可见，在建筑物群外侧拐角等锐缘处，来流的速度大小和方向都发生剧烈变化，且在建筑物群背风侧形成的涡流区内，速度梯度大，风向不稳定。在建筑物群内，易形成强度较小的旋涡区。
3.2 室外风场沿建筑物表面风压分析
建筑物处于大气流场中，由于建筑物形状和空气粘性等因素的影响，使气流速度在建筑物的前后发生变化而引起压强的变化。当风吹响建筑物正面时，因受到建筑物表面的阻挡而在迎风面上产生正压区，气流再向上偏转同时绕过建筑物各侧面及背面，在这些面上产生负压区。因此，当建筑物围护结构存在开口时，由于压差作用，室内就会形成自然通风。建筑物周围的压力分布通常由无因次风压系数描述，及建筑物外表面某点的风压与建筑物同高度出来流风压之比。
在对44号楼三单元不同平面高度的风压系数及风压值的计算结果中，如表1所示。建筑物迎风面的风压系数及风压均随着建筑物高度的增加而增加，且风压从1.6http://forum.vibunion.com/neteaseivp/resource/paper/doc/20059291127958820080/image006.gif变化到17.77http://www.shuigong.com/d/file/papers/dianqi/2005-11-16/03e0f4163f535e538f6708c5fd1d19a2.gif增幅较大；而建筑物背风面则处于很弱的负压作用下，风压系数及风压均相对较小，风压作用很弱。对于西向建筑物，由于其同样具有迎风面与背风面风压差大，风压系数变化明显的特征，因此，风压系数及风压变化规律同上。
由此，对于处于建筑物群迎风侧的建筑物，沿建筑物垂直方向上的风压系数和风压值具有风压差大，风压系数变化明显的特点，因而建筑物高处的通风效果较好。而沿建筑物水平方向上，尽管存在相对不同的风压系数和风压值，但变化规律由对速度场的分析可知，亦存在一定变化规律，即：通风方向均为由建筑物群外侧到建筑物群内侧，且通风效果强。
对于建筑物群背风侧的东向和北向建筑物，结合速度矢量分布和风压分布，采取对个别点的采样分析计算可知：由于建筑物两侧速度绝对值小，方向变化复杂，风压系数和风压沿水平和垂直方向变化均不大，因此，背风面的东北向建筑物具有通风强度较弱，通风方向复杂，规律性不明显等特点。
表1 44号楼不同平面高度的风压系数及风压值 44号楼房间朝向通风器高度Z(m)高度Z处梯度风风速V(Z)(m/s)模拟风速(m/s)动压P(X,Y,Z))(pa)静压Phttp://www.shuigong.com/d/file/papers/dianqi/2005-11-16/f7a8ff962f13bd6a3fe62806e5efaedd.gif(Z)(pa)风压系数Cp(X,Y,Z)3单元102南向2.842.31.651.830.230.45北向2.842.3-0.350.080.23-0.053单元202南向5.943.092.443.990.40.55北向5.943.09-0.470.150.4-0.053单元302南向8.743.62.845.410.610.55北向8.743.6-0.550.20.61-0.053单元402南向11.544.023.186.790.790.55北向11.544.02-0.910.560.79-0.023单元502南向14.344.393.789.590.950.79北向14.344.39-10.670.95-0.023单元602南向17.144.715.0617.181.11.08北向17.144.71-1.070.711.1-0.023单元702南向19.944.955.3218.991.221.08北向19.944.95-1.130.861.22-0.02综上所述，哈尔滨泰海小区44号楼及其周围4栋建筑物，在冬季为西南向主导风的作用下，即：风向投射线与建筑围护结构法线的交角－风向投射角为45°，综合考虑风场和涡流区的关系，认为投射角较恰当，建筑物间距（33.66m≈1.5H）适宜。建筑物群迎风侧的建筑物通风作用明显，通风方向稳定，且应根据以上分析合理地布置建筑物周围环境，改变建筑物周围的气流流场，创造良好的建筑物室内外通风环境。对于建筑物群背风侧的建筑物，也应通过数值模拟计算分析，研究前栋建筑物的阻挡状况以及周围建筑物，寻找特定环境下的通风特点，采取不同的方法和措施，使建筑物室内外获得良好的自然通风环境。
4 总结对于受多种因素和条件影响的住宅小区室外环境，以及应用广泛、功能强大的FLUENT软件，本文仅分析和应用了一小部分内容，随着计算机技术的发展，综合考虑室外太阳辐射、建筑周围绿化等因素将成为生态建筑环境数值研究的一个新方向，而大涡模拟、直接模拟也将会应用的越来越多，使数值模拟技术在实际工程应用中发挥重要作用。
5 参考文献1． Yi Jiang, Qingyan Chen. Effect of fluctuating wind direction oncross natural ventilation in buildings from large eddy simulation.Building and Environment 37(2002)379-386.
2. James W.Axley, Steven J Emmerich, George N.Walton. Modeling thePerformance of a Naturally Ventilated Commercial Building with aMultizone Coupled Thermal/Airflow Simulation Tool. ASHRAE 2002HI-02-21-4
3. Samir S.Ayad. Computational study of nature ventilation. Journal of engineering industrial aerodynamics 82(1999)49-68
4. 孙一坚主编. 简明通风设计手册. 中国建筑工业出版社 2000.2
5. 赵福云,汤广发,周安伟.住宅小区热环境数值分析－风环境的数值模拟.全国暖通空调制冷2002年学术年会论文集（下册）2002.
页:
[1]