风环境数值模拟中模拟植被的数学模型与应用研究.pdfVIP

第十四届全国结构风工程学术会议论文集 风环境数值模拟中模拟植被的数学模型与应用 杨易 金新阳 杨立国 (中国建筑科学研究院风工程研究中心，北京100013) 摘要在建筑风环境数值模拟中，植被的影响经常不可忽略，但模拟流经植被区表面的湍流流动非常复杂．此 类问题通常涉及大尺度的模拟，同时需要模拟重点关注的近地表面的流动特性。在进行大尺度的模拟时，可以通过 修正的壁面率以适当的粗糙度参数来考虑植被的影响．然后，这种方法往往和近地表的网格离散有矛盾：更主要的， 不能捕捉植被区的湍流结构．在本文中，作者采用在动量和湍流榆运方程中附加源／汇项的方法，结合相糙壁面率修 正，以模拟植被绕流。基于SSTk-to模型，首先从理论上推导出端动能耗散频率tO方程的附加源项表达式，然后采 用数值模拟的方法进行了验证．结合作者提出的模拟平衡边界层的数学模型【1-3】，通过和试验结果【4】的比较，表 明本文提出的模型和方法可以较准确的模拟植被绕流．在此基础上，将本文导出的附加源／汇项模型应用到一高层建 筑群风环境的方案优化设计中，以改善建筑群不利风环境品质． 关键词：CFI)、风环境、植被绕流、SST k-co模型，源／汇项 1．引 言 建筑风环境主要研究近地表建筑与周围环境的空气绕流问题。地表的植被覆盖对风环境的影响 通常不可忽略，例如，植被形成的屏障增加了风的阻力，可以有效降低下游遮蔽区的风速，增加湍 流度，影响下游流场中介质的输运。对于数值模拟来说，植被覆盖以及其对地表粗糙度的影响，给 扰流场的数值模拟增加了复杂性。 通常，在进行大尺度边界层流动模拟的时候，一般只关心整体的流动特征，近地表面的流动可 以忽略，这样可以通过考虑粗糙影响的修正的壁面函数模拟地表环境。从数值模拟的角度，使用壁 面函数，也减小了模拟粘性亚层所需的过量的计算网格消耗。然而，这种模拟方法在遇到需要模拟 近地表植被的影响时候，就显得力不从心。在壁面函数中使用粗糙度参数修正的方法，不能考虑植 被的重要特征参数，例如叶面面积密度(LeafArea Density,LAD)，也无法得到所关心的植被遮蔽区 的湍流结构信息。为研究植被覆盖区复杂的流动特性，植被效应的模拟，不仅可以通过表面粗糙度 参数，也可以通过在流动输运方程中引入附加的源／汇项加以考虑【5—8】。 使用粗糙修正的壁面函数方法返回的是整个计算域的整体流动特征。这种方法在粘性亚层需要 的网格数目较少，计算收敛速度更快。在计算量大的复杂的3维流动数值模拟中尤其方便。使用附 加源／汇项模拟可以提供植被遮蔽区详细的流动信息，它需要对输运方程中额外增加的附加源项模型 仔细考虑。文献[8】中采用SSTk-∞模型模拟植被覆盖的地表风场环境，但其中耗散频率∞的附加源 项S∞并不恰当。 本文将这两种方法相结合，以模拟大尺度流动下，流经复杂植被表面的粗糙度效应。首先基于 计算风工程模拟钝体绕流具有较好适用性的SST k．CO模型，从理论上推导出湍动能耗散频率f．O方程 的附加源项表达式，然后采用数值模拟的方法，研究一植被绕流的算例，以此进行验证。在数值模 拟中，结合作者提出的模拟平衡边界层的数学模型【I．3】，通过和试验结果【4】的比较表明，本文提出 的模型和方法可以较准确的模拟植被绕流效应，具有一定的理论和实用意义。 2．湍流模型与附加源／汇项 植被效应通过在动量和湍流模型输运方程中附加源项进行模拟。对于标准k．￡模型，完整的流 第十四届全国结构风工程学术会议论文集 (1)标准k-e模型 连续性方程： 誓：0 (1) 动量方程： c2， 鲁+等=去B+；七)+专{u(考+鲁))+瓯 湍动能k方程： 警+百Ou,k』Ox／f[,旦ok卦即+瓯