近地脉动风场数值模拟可行性研究.pdfVIP

第十四届全国结构风工程学术会议论文集 近地脉动风场数值模拟的可行性研究 朱伟亮 杨庆山 (北京交通大学土木工程学院，北京100044) 摘要符合近地脉动风场特性的LES脉动入口条件是准确模拟结构表面时变风压的必要条件．本文基于LES 模型对平板湍流边界层进行模拟，采用周期边界条件实现脉动输入，并将顺压力梯度引入流体控制方程以维持湍流 边界层的稳定．数值结果表明，采用周期边界条件可成功实现脉动输入，边界层在顺压力梯度下的自保持性良好， 可满足工程需求．顺流向平均流速剖面以及湍流强度剖面均满足低湍流废下的风场特性，数值结果可为我国B类地 貌下抗风研究采用． 关键词结构抗风，大涡模拟、湍流边界层、压力梯度、入流边界 1引 言 土木工程结构一般处于大气边界层甚至城市下垫层内，所处风场是低速、不可压缩粘性流体的 湍流运动，结构周围常伴随着分离、再附着、漩涡脱落等复杂流动现象。考虑到该特征湍流场及结 构表面风压的时变特性，大涡模拟(LargeEddy Fluid Dynamics)方法对其进行模拟的良好选择。而这一过程尚有两个关键问题亟待解决，即：建立对结 构周围特征湍流具有普适性的亚格子模型以及开发具有近地脉动风特性的入流边界条件。目前有关 钝体绕流的LES模拟多采用均匀流入口条件，即忽略了来流的湍流特性。该条件对模拟结构上脉动 风压以及强风作用下局部瞬时风荷载的预测尚显不足，因此有必要对近地湍流风场进行数值模拟， 使其满足平均风剖面、湍流强度剖面以及风场的频谱特性等方面的要求，并将其与结构表面风压模 拟相结合，以得到更符合实际情况的平均风压、脉动风压及风压瞬时值，为主体结构及围护结构的 风致响应提供分析依据。 脉动入口边界是大涡模拟中的一个重要分支，存在其特有的研究难点。一方面，由于LES的数 学模型直接求解大尺度涡的瞬时运动，故不能比拟RANS(雷诺平均)模型的入流边界处理方法， 即通过施加统计量的方式给定如湍动能、耗散率、雷诺应力等参数。另一方面，如果直接将现场或 风洞的风速时程作为已知量施加于入流边界，由于其一般不满足离散后流体控制方程，且CFD中壁 面边界条件处理方法基本不能满足结构风工程对大气边界层的模拟要求，故时程的统计特性有可能 在求解域内迅速衰减或变异，故实验数据也难于直接作为入口条件施加于数值模型。 目前，有关脉动输入的研究尚较少。Kondo”基于蒙特卡洛模拟，考虑风速的能谱及其空间相关 性，生成离散点上脉动风速时程，该风速时程经过处理使其满足零散度要求后被施加于入流面的离 散点上。基于该方法并结合LES，Kondo研究了各向同性湍流运动，并与实验结果作了对比。但该 方法生成的速度场难于满足NS方程，入流面的风场统计特性有可能在计算域内衰减，故尚存在改 入流面内生成二维的脉动漩涡场，以此作为顺流向速度的扰动量。该方法被商用流体求解程序Fluent 引入，用于生成LES的脉动入流，其特点是考虑了入流面上各点的空间相关性，随机脉动生成较快， 节省计算资源。但该方法生成的脉动速度场不满足零散度(连续方程)要求，也不能同时满足NS 方程，故其在计算域内的湍流特性难以得到保持。 952 第十四届全国结构风工程学术会议论文集 如前所述，LES的脉动入口难于通过统计以及实验的方式施加，且人工合成的方式也存在衰减 或变异问题。因此，本文提出通过数值方法直接求解近地湍流边界层，将计算域内某平面的脉动数 值解施加于入流面，以此实现脉动输入。其次，对顺压力梯度下边界层的自保持性进行探讨，并将 该压力梯度引入数值模型，进而求解得到满足近地脉动风场特性的自平衡湍流边界层。 2脉动入口生成方法 近地脉动风场指靠近地球表面、受地面摩擦阻力影响的大气湍流边界层，就整个地球大气而言， 该层只不过是紧贴地球表面的很薄的一层大气。大气边界层与普通物体附面边界层在结构上类似， 均存在对数率层，而该层亦即土木工程结构抗风研究所关心的边界层范围。因此，可将求解物体附 面边界层的思想引入大气边界层的模拟中，并将控制方程无量纲化，从而在普遍意义上实现大气湍 流边界层，即近地脉动风场的模拟。 地面粗糙度是摩擦阻力形成的直接原因。在大气边界层