基于SST k-ω湍流模型的平衡大气边界层模拟.pdfVIP

第30卷第6期 空气动力学学报 V01．30．No．6 2012年12月 ACTAAERODYNAMICASINICA DeC．。2012 文章编号：0258 07 1825(2012)06—0737 基于SST k-co湍流模型的平衡大气边界层模拟 胡 朋，李永乐，廖海黎 (西南交通大学桥梁工程系，四川成都610031) 摘 要：为实现SST融m湍流模型下空流场大气边界层的平衡，提出通过在湍动能输送方程与比耗散率输送方程 中添加源项的方法使来流边界条件与湍流模型相协调，并通过相应的CFD数值算例分别考察了不考虑源项与考 虑源项时沿流场不同位置处的平均速度、湍动能及比耗散率剖面的一致性。计算结果表明，当考虑源项时沿流场 不同位置处的平均速度剖面、湍动能剖面及比耗散率剖面均有较高程度的一致性，较好地实现了大气边界层的平 衡。研究结论为平衡大气边界层的研究提供了新的思路和研究方法，对结构风工程中CFD数值模拟有一定的理 论及应用价值。 关键词：平衡大气边界层；SST肛m湍流模型；源项；计算流体动力学 中图分类号：0357．5 文献标识码：A 相应的风洞试验作对比，结果表明两建筑物中心线上 U 引 百 行人高度处风速放大系数的计算值与试验值相差达 Wind 计算风工程(ComputingEngineering，简15％左右，其原因主要是由于入口处的风速封面在流 称CWE)是结构风工程中极具前景的一个方向，也是动过程中发生了改变，即不满足水平均匀性条件，从 当前国际风工程研究的热点。CWE的核心内容是计 而使计算风速放大系数的参考风速发生了变化，当根 Fluid 算流体动力学(ComputationalDynamics，简称据风速剖面变化对参考风速进行修正后，风速放大系 CFD)，与传统的风洞试验相比，CFD方法具有费用数的计算值与试验值吻合较好，文中还建议在进行数 低、周期短、便于模拟真实环境、描述流场细节及给出 值模拟之前，应首先对无模型的空流场进行水平均匀 流场定量结果等优点。自20世纪60年代兴起以来，性检验。严格说来，当且仅当计算区域入口处的来流 随着数值求解方法的发展和计算机硬件技术的进步， 边界条件与湍流模型(包括壁面函数)及地表粗糙特 CFD在结构风工程中的应用越来越广泛。 性相协调时，才能得到水平均匀性的平衡大气边界 由于结构风工程研究的对象位于大气边界层内， 层[3]，但由于湍流的复杂性，在数值模拟中精确实现 研究的重点是钝体绕流，其流动可视为低速、不可压 缩粘性流体的湍流运动，因此准确模拟大气边界层内 基于融e湍流模型提出了满足平衡大气边界层的几 的湍流流动是精确预测结构风荷载的前提，这在 个原则，给出了一组来流边界条件，并对湍流模型参 CFD数值模拟中就首先要求其满足水平均匀性【1]的数的取值进行了探讨，但文中所给出的湍动能边界条 条件，即在没有任何障碍物的简单边界层流动下，数 件为一常数，这与实测结果及试验结果不符。张建和 值模拟区域内的平均风剖面和湍流风剖面自人El到 杨庆山口1基于标准融￡湍流模型，将近壁区第一层网 出口保持不变，这种满足水平均匀性条件的边界层称 格的湍流量进行平均处理并施加剪切应力，从而使边 为平衡大气边界层，边界层内流场的水平均匀性程度 界层内沿流场不同位置处的流动特征量吻合程度较 直接影响数值模拟的精度。Blocken等[21采用融￡湍高。方平治等[53通过修正的壁面函数并采用 流模型模拟了两平行建筑物内的风速分布情况并与 是一e湍流模型