建筑LES模拟脉动入口的适用性研究.PDF

第29卷第1期 空气动力学学报 v01．29。No．1 201 1 1年2月 ACTA Feb．，201 AERODYNAMICASINICA 文章编号：0258—1825(2011)01-0124．05 建筑LES模拟脉动入口的适用性研究 朱伟亮1，杨庆山1，曹曙阳2 (1．北京交通大学土木工程学院，北京100044；2．同济大学土木工程防灾国家重点实验室上海200092) 摘要：采用周期人流边界，在展向等间距网格系统中对平板湍流边界层(驱动单元)进行LES模拟，并将顺压力梯 度引入流体控制方程以维持湍流边界层厚度的稳定。记录驱动单元脉动时程并将其引入虚拟模拟网格系统(验证 单元)以研究LES脉动入口方法的适用性。数值结果表明，采用周期人流边界可成功实现脉动输入，边界层在顺压 力梯度下的自保持性良好；脉动入流特性在驱动单元内保持良好，可满足工程需求；顺流向平均流速剖面以及湍流 强度剖面均满足低湍流度下的风场特性。数值结果可为我国B类地貌下抗风研究采用。 关键词：结构工程；结构抗风；大涡模拟；脉动入口；压力梯度 中图分类号：TU973+．213 文献标识码：A 时程，并利用其研究了各向同性湍流运动。该方法较第 0 引 言 一类更为直接，但生成的速度场难于满足NS方程，人 土木工程结构一般处于大气边界层甚至城市下垫 流面的风场统计特性有可能在计算域内衰减，作者并未 层内，结构周围常伴随着分离、再附着、漩涡脱落等复杂 对该问题进行说明。方法一的优势在于获得的流场时 流动现象。考虑到该特征湍流场及结构表面风压的时 程满足驱动单元中离散后的NS方程及连续方程。然 变特性，大涡模拟成为采用CFD方法对其进行模拟的 而一般情况下，驱动单元与模拟单元的(展向)网格划 良好选择。而符合近地脉动风场特性的脉动入口条件 分并不匹配，前者的流场解能否在后者中良好保持需要 则是这一模拟过程的必要前提条件，其对结构上脉动风 数值计算的验证。 压的模拟以及强风作用下局部瞬时风荷载的预测有着 本文基于第一类方法生成脉动时程，首先对驱动单 极为重要的作用。 元得到流场进行分析，继而对其作为模拟单元入口条件 脉动入口条件的生成一般分为两类。第一类，预先 的有效性(自保持性)进行检验。 在人流驱动单元(以下简称驱动单元)进行风场(湍流 边界层)的LES模拟，生成并记录所模拟风场的速度压 力时程，将其作为入口条件引入绕流模拟单元(以下简 称模拟单元)，如图1所示。基于该思想，Lund…对平板 脉动人流进行了模拟，验证了流场在平板上的自保持 性。NozawaBl等将Lund的方法发展到粗面边界层，但 图1脉动输入方法 并未对流场在模拟单元内的保持性进行探讨。第二类 Turbulentinflowmethod Fig．1 即在绕流模拟单元的人口边界，通过逆傅里叶变换结合 零散度修正，人工生成入口边界节点速度，其将随空间 1脉动入口生成方法 和时间变化。Kondopl基于蒙特卡洛模拟，考虑风速谱1．1入流驱动单元 及其空间相关性，经零散度修正后生成离散点上的风速 对于驱动单元湍流边界层的模拟，为缩减计算时 ·收稿日期：2010一Ol一嚣；修订日期：20