STU大气边界层风洞流场特性的模拟
大气边界层模拟风洞研究综述

(结构动力试验论文)大气边界层模拟风洞研究综述Performance of Simulation of Atmospheric BoundaryLayer in Wind Tunnels学生姓名:指导教师:学院:专业班级:大气边界层模拟风洞研究综述姓名(大学学院)摘要:本文介绍了大气边界层风洞的发展过程和模拟方法。
大气边界层的模拟方法主要有主动模拟方法和被动模拟方法,前者包括多风扇风洞技术与振动尖塔技术,后者采用尖劈、粗糙元、挡板、格栅等装置进行模拟。
被动模拟技术较为经济、简便,所以得到了广泛采用。
关键词:风洞;大气边界层;主动模拟;被动模拟.Performance of Simulation of Atmospheric Boundary Layerin Wind TunnelsNAME(University)Abstract:In this paper , the simulation of atmospheric boundary layer are introducted from the history of the development and the methods of the technology. The methods of atmospheric boundary layer simulation contain active simulation and passive simulation. The active simulation mainly include multiple fans wind tunnel technology and vibratile spire technology. The equipments of the passive simulation main include spire, roughness element, apron and gridiron. The passive simulation technology is simple and economical, so it has been widely used.Key words:wind tunnel; atmospheric boundary lay er; active simulation; passive simulation.一、引言1940年,美国塔科马悬索桥由于风致振动而破坏的风毁事故,首次使科学家和工程师们认识到了风的动力作用的巨大威力[1]。
平衡大气边界层自保持及TTU风压数值模拟

平衡大气边界层自保持及TTU风压数值模拟杨磊;李磊;杨向龙【摘要】构造合理且满足水平均匀性的大气边界层是计算风工程的一个重要难题.基于雷诺平均方法,采用标准k-ε湍流模型,结合Richards&Hoxey入口条件,对数值风洞的顶面边界条件施加剪切应力修正,进行二维空流场的平衡大气边界层模拟,给出了一类能够较好满足大气边界层自保持要求的流场边界条件.并将此类边界条件同时运用到标准k-ε湍流模型和SST(shear stress transport) k-ω湍流模型中,针对美国德克萨斯技术大学模型场地试验进行三维流场全尺寸数值模拟,验证了所提边界条件对上述两类双方程湍流模型在低矮建筑钝体绕流的适用性.【期刊名称】《深圳大学学报(理工版)》【年(卷),期】2016(033)005【总页数】6页(P470-475)【关键词】计算风工程;平衡大气边界层;水平均匀性;标准k-ε湍流模型;顶部边界条件;德克萨斯技术大学模型【作者】杨磊;李磊;杨向龙【作者单位】深圳大学土木工程学院,广东深圳518060;深圳大学土木工程学院,广东深圳518060;深圳大学土木工程学院,广东深圳518060【正文语种】中文【中图分类】TU973.32计算风工程(computational wind engineering, CWE)是利用计算流体动力学(computational fluid dynamics, CFD)方法研究不同类型建筑结构风荷载和风环境的学科.其中,准确模拟大气边界层内的湍流流动是利用数值方法研究结构风特性的重要前提.计算风工程应首先满足大气边界层(atmosphere boundary layer, ABL)的水平均匀性(horizontal homogeneity)要求.水平均匀性是指在没有障碍物的简单边界层流动下,数值风洞内风剖面的物理参数与测量地点无关,即计算域入口处定义的流动物理量(速度、湍动能和湍流耗散率等),从入口到出口保持一致.满足水平均匀性要求的边界层称为平衡大气边界层.大气边界层的水平均匀性对CFD模拟结果的准确性有显著影响.在数值风洞内平衡大气边界层,影响自保持水平均匀性的主要因素包括边界条件、湍流模型、壁面处理和网格划分等方面.Richards等[1-3]提出入流边界、湍流模型和地面粗糙度的综合作用会产生一个内边界层(internal boundary layer, IBL),使得风剖面在到达建筑物之前发生较大的改变;同时研究了来流边界条件的合理给定问题,给出了一组基于标准k-ε湍流模型的入口边界条件(Richards & Hoxey (RH)入口条件).Juretic'等[4]通过研究雷诺应力结合壁面函数的修正,给出了适用于标准k-ω湍流模型的来流边界条件.曾锴[5]基于shear stress transport (SST) k-ω湍流模型,考察了不同入口边界条件的设置对钝体绕流计算结构的影响,提出了入口湍流风剖面自保持的表达方法.杨伟等[6-8]从模型方程本身出发,基于风洞实验数据,推导出一类近似满足k-ε模型自保持边界要求的入口湍动能表达式,并定义了模型常数.顾明等[9]研究了壁面Y+值对平衡大气边界层数值模拟结果的影响,结果表明相比于近壁面函数,通过近壁面网格来预测湍流的产生能提高模拟的准确性.孙玉航等[10]针对平衡大气边界层,基于标准k-ε湍流模型,通过推导耗散率表达式,并在输运方程中添加源项的方法,得到了一类较为适用的边界条件.对于计算域顶部,常用的边界条件有对称边界条件[11]、速度边界条件[12]、压力出口边界条件[13]和剪切应力边界条件[1,11].顶部边界条件对计算结果会产生较大影响.从理论上来说,剪切应力边界条件更符合物理实际,并且在计算域顶部施加剪切应力边界条件确实可以得到很好的结果.然而,正如Hargreaves等[11]指出的,这种边界条件在通用的CFD模拟软件(如Fluent和CFX等)中难以实现,因此,被很多研究者忽略.Richards等[1]提出了通过在靠近顶部边界的一层网格中增加动量和湍流耗散率源项的方法来施加剪切应力的思想,但他们并没有给出相应的源项公式.目前已有的数值模拟研究中,尚未构造出全部主要特征量均能够很好满足水平均匀性的数值风洞.本研究基于标准k-ε湍流模型,采用RH入口条件,提出并验证了一种便于实现并能较好满足大气边界层自保持要求的顶部边界条件.采用此类边界条件对美国德克萨斯技术大学(Texas Tech University, TTU)建筑周围风场进行了数值模拟,将模拟结果与实测数据进行分析对比,验证了该方法的可行性.基于ANSYS—Fluent 15.0平台,利用标准k-ε湍流模型,研究顶部边界条件对平衡大气边界层自保持特性的影响.针对无建筑物扰动的空流场特征,采用二维数值模拟以突出重点和提高效率.计算区域尺寸为1 000 m×300 m,计算域内未放置任何障碍物.采用结构网格对计算域进行离散化处理,近地面沿高度方向最小网格大小为0.01 m,竖向增长因子为1.05,网格单元总数约为40 000个,如图1.计算域入口统一采用velocity-inlet边界条件,速度和湍动能以及湍流耗散率的数值利用RH入口条件确定.来流风速剖面u(z)采用对数率表达为湍动能k为湍流耗散率ε为其中,u*为摩擦速度,根据参考高度处的参考风速确定;κ为Von Karman常数,取κ= 0.41; z为离地面的高度; z0为粗糙度长度,取z0=0.01 m;Cμ为k-ε湍流模型常数,取0.09.计算域地面采用固壁边界条件,引入粗糙壁面修正.粗糙度常数Cs、粗糙度高度Ks以及粗糙度长度z0满足式(4),CsKs≈9.793z0其中, Cs=0.5, Ks=20, z0=0.2.出流面采用完全发展出流边界条件为其中, x表示顺流方向; u、 v和w分别为3个方向的速度分量.对于计算域顶部边界条件,本研究提出了一种容易在通用CFD模拟软件中实现的施加剪切应力的方法.对紧靠顶部边界的两层网格,直接施加满足式(1)至式(3)的速度、湍动能和湍流耗散率的条件.通过给定顶部边界附近的物理梯度的方法,达到施加剪切应力的目的.为比较分析顶部边界条件对大气边界层自保持特性的影响,采用4种不同顶部边界条件计算流域顶部,分别为对称边界条件、指定速度和湍动能的速度入口边界条件、静压边界条件以及在计算域顶部施加剪切应力边界条件.流动视为定常不可压缩流,采用速度-压力耦合方式(semi-implict method for pressure-linked equations consistent,SIMPLEC),对流项采用二阶迎风格式进行离散化.当所有相关变量的迭代残差达到平稳且监测的某些典型物理量不再显著变化时,判定计算收敛.为比较各工况下平衡大气边界层效果,选取入流面、入流面下游300和600 m及出流面4个不同观测位置相应的速度、湍动能和湍流耗散率等关键物理量进行比较.图2为4种不同计算工况下,4个不同观测位置平均速度随高度变化曲线.从图2可见,工况1和工况3的速度曲线在顶部出现了一定的分离,相比工况1的轻微分离,工况3的分离现象更加明显,表明在RH入口条件下,顶部采用压力出口边界条件对速度剖面的水平均匀性有着不利的影响.而工况2和工况4的条件下,各观测位置的速度曲线都能很好的吻合,表现出很好的水平自保持效果.图3为4种不同计算工况下,4个不同观测位置湍动能随高度变化曲线.由图3可见,工况1至工况3的下游截面湍动能在壁面附近均大于入口处,这种现象在很大程度上是由于数值离散格式所引起的.随着高度的增加,湍动能逐渐减小,顶部位置的湍动能均小于入口处.特别是工况1和工况3,顶部位置的湍动能出现了明显的衰减,这表明顶部设置为对称边界条件和压力出口边界条件,对湍动能剖面的水平均匀性有不利影响,图3(d)所示工况4中,边界层湍动能剖面在大部分区域保持良好,特别是在顶部,基本与入口处完全吻合,而在近地面小范围区域存在一定的误差.这种误差可能与壁面粗糙度的设置有关.入口边界条件所设定的来流中包含的湍流在一定程度上反映出地面粗糙性质,而该性质与计算域地面的固壁边界条件所包含的粗糙性质存在一定的差异,因此,流场在向下游运动的过程中会反映出该差异性.相比工况1至工况3,工况4在顶部施加剪切应力边界条件后,不但能有效减小计算域顶部的湍动能差异,且能有效降低近地面小范围区域内的湍动能误差.综上分析表明,入口采用RH入口条件,顶部施加剪切应力边界条件,能够获得较好的平衡大气边界层,实现速度、湍动能和湍流耗散率的自保持性.TTU实际尺寸场地试验由美国德克萨斯理工大学的风工程研究现场实验室于1991年完成[14-15].该实验室建造了一个永久的全尺寸金属建筑用以进行风压的测量.该试验建筑原始尺寸为13.7 m×9.2 m×3.9 m,建筑周围地势开阔平坦,建筑物表面光滑.作为具有代表性的现场实测数据,该试验已成为国际公认的权威的评估数值模拟方法的基准试验,如图4.其中,1#~11#为测试点.文献[16-17]研究结果表明,SST k-ω湍流模型由于能够有效计算湍流切应力在逆压梯度边界层的输运,更适用于钝体分离流模拟.因此,本研究分别采用标准k-ε湍流模型和SST k-ω湍流模型,边界条件采用第2节工况4的设定,针对TTU试验的90°风攻角进行了三维全尺寸数值模拟,并将数值模拟结果与试验结果进行对比.计算域大小为200 m×100 m×50 m.建筑物放置在入口下游80 m处.设建筑物的长、宽、高分别为L、 B和H,模型距上游为9L,距下游为13L,计算域宽度为7.5B,在竖直方向上为12.5H,模型阻塞率为1.1%.计算域采用六面体网格进行划分,网格总数约150万个.图5为两种湍流模型的数值模拟结果与TTU模型场地实测数据的对比.可见,在90°风攻角的典型计算工况下,采用标准k-ε和SST k-ω两种湍流模型得到的计算结果相差不大,与试验数据总体吻合,均落在TTU建筑现场实测数据的包络线范围内.在模型的迎风面,标准k-ε湍流模型略低于试验结果的平均值,而SST k-ω湍流模型略高,最大偏差为7.4%.对于4#测点之后的钝体绕流分离区,两种湍流模型得到的数值模拟结果均略低于TTU试验结果的平均值.最大偏差出现在4#测点,即建筑物迎风面与顶面相交的尖角分离处,此测点处标准k-ε和SST k-ω两种湍流模型的计算结果基本位于试验结果的下限,与试验平均值的误差分别为9.3%和7.9%.在数值模拟通常较为困难的背风面,计算结果整体与试验数据吻合良好.基于ANSYS—Fluent 15.0中的标准k-ε湍流模型,利用RH入口条件,通过本研究方法对顶面施加剪切应力边界条件,能够获得较好的平衡大气边界层,实现整个流域内速度、湍动能和湍流耗散率的自保持特性.在平衡大气边界层的条件下,采用标准k-ε湍流模型和SST k-ω湍流模型均能够较好地对TTU模型场地试验进行数值模拟,分辨钝体流动特征.[1] Richards P J, Norris S E. Appropriate boundary conditions forcomputational wind engineering models using the k-ε turbulence model[J]. Journal of Wind Engineering and Industrial Aerodynamics, 1993,s46/s47(93): 145-153.[2] Richards P J, Younis B A. Comments on ‘Prediction of wind-generated pressure distribution around buildings’ by Mathews E H[J]. Journal of Wind Engineering and Industrial Aerodynamics, 1990, 34(1): 107-110. [3] Richards P J, Quinn A D, Parker S. A 6 m cube in an atmospheric boundary layer flow. Part 2, computational studies[J]. Journal of Wind and Structures, 2002, 5(2/3/4): 177-192.[4] Juretic' F, Kozmar H. Computational modeling of the neutrally stratified atmospheric boundary layer flow using the standard k-ε turbulence model[J]. Journal of Wind Engineering and Industrial Aerodynamics, 2013, 115:112-120.[5] 曾锴. 计算风工程入口湍流条件改进与分离涡模拟[D]. 上海: 同济大学, 2007. Zeng Kai. The improvement of inlet turbulence boundary condition in CWE and detached eddy simulation[D]. Shanghai: Tongji University, 2007.(in Chinese)[6] 杨伟, 顾明, 陈素琴. 计算风工程中的k-ε模型的一类边界条件[J]. 空气动力学学报, 2005, 23(1): 97-102.Yang Wei, Gu Ming, Chen Suqin. A set of turbulenceboundary condition of k-ε model for CWE[J]. Acta Aerodynamica Sinica, 2005, 23(1): 97-102.(in Chinese)[7] Yang Wei, Quan Yong, Jin Xinyang, et al. Influences of equilibrium atmosphere boundary and turbulence parameter on wind loads of low-risebuildings[J]. Journal of Wind Engineering and Industrial Aerodynamics, 2008, 96(10/11): 2080-2092.[8] 杨伟, 金新阳, 顾明, 等. 风工程数值模拟中平衡大气边界层的研究与应用[J]. 土木工程学报, 2007, 40(2): 1-5.Yang Wei, Jin Xinyang, Gu Ming, et al. A study on the self-sustaining equilibrium atmosphere boundary layer in computational wind engineering and its application[J]. China Civil Engineering Journal,2007,40(2): 1-5.(in Chinese)[9] 顾明, 李孙伟, 周印. 平衡大气边界层的数值模拟和风洞实验[J]. 同济大学学报自然科学版, 2009, 37(3): 298-302.Gu Ming, Li Sunwei, Zhou Yin. Numerical and experimental studies on self-sustaining equilibrium atmosphere boundary layer[J]. Journal of Tongji University Natural Science, 2009, 37(3): 298-302.(in Chinese)[10] 孙玉航, 王国砚, 王梦虹. 一类平衡大气边界层边界条件研究[J]. 力学季刊, 2015, 36(2): 328-335.Sun Yuhang, Wang Guoyan, Wang Menghong. A boundary condition of equilibrium atmospheric boundary layers[J]. Chinese Quarterly of Mechanics, 2015, 36(2): 328-335.(in Chinese)[11] Hargreaves D M, Wright N G. On the use of the k-ε model in commercial CFD software to model the neutral atmospheric boundary layer[J]. Journal of Wind Engineering and Industrial Aerodynamics, 2007,95(5) : 355-369.[12] Blocken B, Stathopoulos T, Carmeliet J. CFD simulation of the atmospheric boundary layer-wall function problems[J]. AtmosphericEnvironment, 2007,41(2): 238-252.[13] Kose D A, Fauconnier D, Dick E. ILES of flow over low-rise buildings: Influence of inflow conditions on the quality of the mean pressure distribution prediction[J]. Journal of Wind Engineering and Industrial Aerodynamics, 2011,99(10): 1056-1068.[14] Levitan M L, Mehta K C, Vann W P. Field measurements of pressures on the texas tech building[J]. Journal of Wind Engineering and Industrial Aerodynamics, 1991, 38(10/11): 227-234.[15] Levitan M L, Mehta K C. Texas tech field experiments for wind loads part 1: building and pressure measuring system[J]. Journal of Wind Engineering and Industrial Aerodynamics, 1992, 43(1/2/3): 1565-1576. [16] Wilcox D C. Turbulence modeling for CFD[M]3rd ed. Los Angeles, USA: DCW Industries, Inc., 1994.[17] Taghinia J, Rahman M, Siikonen T. Simulation of indoor airflow with RAST and SST-SAS models: a comparative study[J]. Building Simulation, 2015, 8(3): 297-306.。
大气边界层的风洞模拟

大气边界层的风洞模拟
王兆印
【期刊名称】《实验力学》
【年(卷),期】1998(13)3
【摘要】本文通过比较大气边界层与不同比尺的模型实验,研究了风洞模型实验
中的比尺效应。
结果表明,大气边界层粗糙高度的比尺应当和边界层厚度比尺相等。
边界层厚度的比尺在200~400之间时可以得到最好的速度分布的相似性,表面压强分布主要依赖于边界层比尺而与建筑物几何比尺关系不大。
如果正确选择了粗糙高度比尺,扩散系数比尺和边界层厚度比尺相同。
【总页数】11页(P283-293)
【关键词】大气边界层;风洞模型;比尺效应;粗糙高度
【作者】王兆印
【作者单位】国际泥沙研究培训中心
【正文语种】中文
【中图分类】P433
【相关文献】
1.大气边界层湍流特性的风洞模拟实验研究 [J], 贺林;王军;Demetri Bouris
2.大气边界层流场的模拟与大气边界层风洞 [J], 黄东群;马健
3.TJ-2风洞大气边界层风场数值模拟 [J], 雷文治;吉柏锋
4.下击暴流出流段非稳态风场的大气边界层风洞模拟 [J], 钟永力;晏致涛;李妍;杨
小刚;蒋森
5.海上采油平台锚泊状态和拖航状态海面上大气边界层模拟以及海面表层中海流模拟的风洞试验 [J], 施宗城;张大春
因版权原因,仅展示原文概要,查看原文内容请购买。
大气边界层风场特征模拟与分析研究

大气边界层风场特征模拟与分析研究在大气科学领域中,大气边界层是指地球表面与上层大气之间的过渡层,对于气象预测、空气质量评估等领域具有重要意义。
而其中的风场特征是其研究的重点之一。
本文将探讨大气边界层风场特征的模拟与分析方法。
一、背景介绍大气边界层的风场特征是受多种因素综合影响的,在不同的地理环境中表现出差异性。
由于大气边界层风场的复杂性,为了更好地理解和模拟其特征,研究者们提出了不同的方法和技术。
二、数值模拟方法数值模拟方法是研究大气边界层风场特征的常用手段之一。
该方法通过数学模型和计算机算法来模拟大气边界层风场的演变过程。
常见的数值模拟方法包括气象模式、雷达回波模拟等。
1. 气象模式气象模式是通过一系列参数和方程来描述大气运动、能量传递和物质运输的数值模型。
通过运行气象模式,可以模拟大气边界层的风场特征。
常用的气象模式包括欧洲中期天气预报中心的欧洲中期天气预报模式(ECMWF)、美国国家环境预报中心的全球预报系统(GFS)等。
2. 雷达回波模拟雷达回波模拟是利用雷达回波数据和数学模型,通过计算和推算得到大气边界层的风场信息。
通过分析雷达回波的特征,可以获取大气边界层风场的分布和运动情况。
三、实测资料分析方法除了数值模拟方法外,实测资料的分析也是研究大气边界层风场特征的重要手段之一。
通过各种地面、航空、卫星观测站点所获取的实测数据,可以对大气边界层的风场特征进行分析。
1. 地面观测站点地面观测站点是通过建立气象观测站网络,采集并记录大气各种要素的实测资料。
通过对地面观测站点资料的分析,可以得到不同地理环境中大气边界层风场的特征。
2. 航空观测资料航空观测资料是通过飞机或无人机等航空平台所采集的数据。
通过对航空观测资料的分析,可以获取大气边界层风场在垂直方向上的变化情况,进而揭示其垂直结构特征。
3. 卫星观测资料卫星观测资料是通过卫星对地球表面进行遥感探测所获取的数据。
卫星观测资料具有广覆盖区域、高时空分辨率的特点,通过对卫星观测资料的分析,可以更全面地认识大气边界层风场特征。
大气边界层风洞流场特性的模拟

本 试 验 在汕头 大 学风 洞 实 验 室 的 S D T X一1大气
边 界层 风 洞进 行 ,T X一1 一座 具 有 串置 双 试 验段 SD 是 的全 钢 结 构 的 闭 口 回流 低 速 风 洞 。其 主 试 验 段 宽 3 高 2 长 2 m。 风 速 连 续 可 调 , 大 风 速 可 达 m、 m、 0 最
大气 边界 层 的风 洞 模 拟 是 研究 风对 建 筑 结 构 的 作 用 以及 建筑 物 周 围 的风 环 境 等 风 工程 问题 的基 础 工作 。在 进行 建 筑结 构 的风 洞试 验 时 , 了必 须 满足 除 建筑 结构 的几 何 特性相 似 外 , 还要 准 确地模 拟 出建筑 结构 所 在 区域 的 大 气 边 界 层 特 性 。 大气 边 界层 的特 性 主要包 括平 均 风 速 剖 面 、 流 度 剖 面 、 功 率 谱 以 湍 风 及湍 流 积分尺 度 等 。尤 其 是 对 于 研 究 建 筑 结构 风 致
大 气 边 界 层 风 洞 流 场 特 性 的模 拟
石碧青, 洪海波, 谢壮宁, 倪振华
( 头 大 学 土木 工程 系 , 东 汕 头 汕 广 556 ) 103
摘
要 : 洞 中 大 气 边 界 层 模 拟 的 准 确 度 是 保 证 风 洞 试 验 结 果 准 确 度 的重 要 因 素 。本 文 采 用 了 一 种 新 型 曲边 梯 形 风
12 模 拟装 置 .
采用 被 动模 拟方 法模 拟 流场 , 般 情况 是采 用 在 一 风洞 试验 段 人 口加尖 塔 、 加挡 板 和沿 风 洞 地板 分布 粗
在4m×3m风洞中模拟大气边界层

在4m×3m风洞中模拟大气边界层
李会知;樊友景;庞涛;王婕
【期刊名称】《郑州大学学报(工学版)》
【年(卷),期】2003(024)004
【摘要】通过适当地延长风洞的试验段长度,并利用尖塔、挡板、粗糙元等边界层发生装置,在4 m×3 m航空风洞中模拟建立了大气边界层流场,测量了模拟流场的平均风速剖面、湍流强度剖面、风谱等.结果表明:在4 m×3 m风洞中建立的模拟大气边界层流场适宜于建筑物风荷载、风环境模型试验.
【总页数】3页(P90-92)
【作者】李会知;樊友景;庞涛;王婕
【作者单位】郑州大学土木工程学院,河南,郑州,450002;郑州大学土木工程学院,河南,郑州,450002;河南省建筑设计研究院,河南,郑州,450005;河南省建筑设计研究院,河南,郑州,450005
【正文语种】中文
【中图分类】V211.73
【相关文献】
1.风洞中模拟大气边界层的实验要求 [J], 金文;王元
2.风洞中模拟大气边界层的相似参数分析 [J], 金文
3.数值模拟大气边界层风洞中湍流模型的比较 [J], 梁村民;陈治全
4.风洞短试验段中基于被动技术的大气边界层模拟 [J], 李永乐;卢伟;李明水;廖海
黎
5.在1.4米航空风洞中模拟大气边界层 [J], 李会知;刘忠玉;郑冰;吴义章
因版权原因,仅展示原文概要,查看原文内容请购买。
数值模拟大气边界层风洞中湍流模型的比较

数值模拟大气边界层风洞中湍流模型的比较梁村民;陈治全【摘要】风洞模拟试验是风工程研究的一种重要手段,在对大气边界层风洞进行数值模拟的时候,湍流模型的选取是影响模拟计算结果的至关重要的因素之一,本文针对大气边界层的风洞模型为考察对象对湍流模型进行了考察,考虑了各种湍流模型对计算结果的影响,并与实验结果进行了对比,旨在得出大气边界层风洞模拟中比较适合的湍流模型.研究发现在模拟大气边界层的风洞时,一般宜采用RNG k-ε模型,在计算资源和时间允许的情况下,建议使用雷诺应力模型.【期刊名称】《四川建材》【年(卷),期】2012(038)002【总页数】2页(P273-274)【关键词】空气动力学粗糙度;大气边界层厚度;湍流强度【作者】梁村民;陈治全【作者单位】中国建筑西南设计研究院有限公司,四川成都610042【正文语种】中文【中图分类】TU120 前言数值风工程就是通过计算机数值模拟的方法来研究风工程问题,是近年来兴起的发展很快的一门学科,采用数值模拟方法来预测和评估建筑物的风压分布和周边的风环境问题已越来越被学术界和工程界所接受。
风洞模拟试验是风工程研究的一种重要手段,在边界层风洞中正确复现大气边界层流动特性,是试验结果可信的必要条件,也是风工程研究的重要基础工作。
本文就是针对这种情况用数值模拟的方法(CFD)对风洞中的大气边界层进行了研究,模拟了在风洞中粗糙元和大气边界层之间的关系。
进行数值模拟时的主要手段是CFD技术,但在进行模拟时有诸多因素会直接影响计算结果的收敛性和精确性,如湍流模型、边界条件、网格划分等等。
其中湍流模型的合理选择就是其中重要的因素之一,也引起研究人员的广泛关注。
本文就是通过选取钝体绕流模拟过程中常见的湍流模型,包括RNG κ-ε模型、realizabl eκ-ε模型、SSTκ-ω模型、标准κ-ω模型、RSM模型,来研究风洞中第二试验段得到的风速廓线,通过与试验中得到的风速廓线的比较,选取合适的湍流模型对大气边界层风洞进行模拟。
第五章:大气边界层风特性及其风洞模拟技术

五:大气边界层风特性及其风洞模拟技术
5.2 大气边界层风特性
5.2.1 平均风剖面(廓线) 幂函数律
¾ 幂函数律是历史上最早被用来
描述水平均一地貌上的平均风 剖面的函数形式:
U ( z ) = U r ( z zr ) = U G ( z H G )
α α
z α — 幂函数指数,常被称为地表粗糙度指数,
U (Z)
结构抗风试验
五:大气边界层风特性及其风洞模拟技术
5.2 大气边界层风特性
5.2.2 阵风因子
¾ 阵风风:
Gv = U g U
¾ 阵风因子是一个反映风速阵性程度的参数,影响大小的因
素有许多:
z 平均风速的时距越短,平均风速越大,阵风因子越小 z 天气过程变化越剧烈,阵风因子越大 z 雷暴大风的阵风因子最大,热带风暴次之,寒潮大风(冷空气)最小 z 地表越粗糙阵风因子越大 z 离地越近阵风因子越大 z 平均风速越小阵风因子越大
0.844 0.958
结构抗风试验
五:大气边界层风特性及其风洞模拟技术
5.2 大气边界层风特性
5.2.1 平均风剖面(续) 幂函数律-对数律比较
600 550 500 450 400 350 300 250 200 150 100 50 0
A 类:HG=300m, UG=50m/s:
幂函数律 α= 0.12 对数律 z0=0.01
@ w(t) —与平均风方向垂直的竖向脉动分量,简称为竖向脉动风速
结构抗风试验
五:大气边界层风特性及其风洞模拟技术
5.2 大气边界层风特性
¾ 平均风速计算时距
z 需要足够长 z 据国际上气象领域的实践经验,在大多数国家(包括中国)平均
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
的,但湍流度分布有待改善。 采用曲边梯形尖塔,结合粗糙元进行边界层模拟,提高了边界层中上部的湍流度,使湍流度分
布更加合理。应该指出的是:湍流积分尺度也是衡量一个重要….,,模拟有相当的难度。
参考文献
[1] Nicholas Isyumov , Wind tunnel studies of buildings and structures / Aerospace Division of the American Society of Civil Engineers, Reston, Virginia, 1999, pp65-81, ISBN 0-7844-0319-8
自 Counihan(1969)和 Standen(1972)提出两种模拟大气边界层的装置以来,大气边界层的风 洞模拟技术和方法有了较大的发展。Irwin(1981)对 Standen 尖塔进行了改进,并进一步给出了尖 塔、粗糙元组合模拟梯度风高度的经验公式 1)。相比主动控制技术模拟大气边界层,被动模拟方法 所需的装置简单、造价低,并且比较容易模拟出各种地貌的平均风剖面。因此,风洞中更多地采用 被动模拟方法进行大气边界层的模拟。在国内外同行模拟的大气边界层特性中,平均风剖面的结果 是令人满意的,而湍流度分布、湍流积分尺度等湍流结构则有较明显的差别 2)3)。被动模拟方法中常 用的尖塔形状是三角形,但是三角形尖塔所模拟的湍流度随高度衰减太快,导致风洞大气边界层的 中上部分湍流度过低。本文采用了一种曲边梯形的尖塔,结合粗糙元等被动装置进行风洞大气边界 层的模拟,给出了 A、B、C 和 D 类四种地貌的平均风速剖面、湍流度分布和顺风向脉动风谱和湍 流积分尺度。从模拟的结果可以看出曲边梯形尖塔相对于三角形尖塔的更有利于提高边界层中部以 上高度湍流度,使得湍流度剖面更为合理。
2. 模拟装置
汕头大学风洞的主试验段尺寸为长 20m、宽 3m 和高 2m,风速从 0.5m/s~45m/s 连续可调。大 气边界层的模拟主要使用尖塔旋涡发生器和一定数量的粗糙元来进行。不同地貌所使用的尖塔尺寸 不尽相同,根据多年的实践摸索,得到了用于不同地貌的尖塔,见图 1,尺寸单位为 cm。尖塔迎风 面的基本外形为曲边梯形,高为 200cm。曲边可以采用多项式表示为
[2] Cesar Farell, Arun K.S. Iyengar, Experiments on the wind tunnel simulation of atmospheric boundary layers, J. Wind Eng. Ind. Aerodyn., 1999, 79, pp11-35.
图 3 四类地貌平均风剖面及湍流度分布
(a) 32.6cm 高度
(b) 120.1cm 高度
图 4 C 类地貌不同高度顺风向脉动风谱
4. 讨论
对于在风洞中采取被动模拟方法,Irwin 给出了尖塔和粗糙元设计的经验公式 5),但通过对四种 不同地貌的大气边界层的模拟,发现并不能简单地套用经验公式模拟得到满意的结果,还需要对模
拟装置进一步改进才能得到预定目标。在这些模拟装置中,尖塔的作用最为显著,它的基本外形尺 寸确定了风剖面的大致形状。确定合适的尖塔外形后,通过改变粗糙元的分布可以模拟出所需的风 剖面。三角形尖塔是采用较多的一种方式,它可以方便地模拟出不同地貌的平均风速剖面。通过改 变底部宽度可以提高中下高度的湍流度,但是湍流度沿高度方向衰减过快,使得中上高度的湍流度 则偏小,达不到要求,它在湍流度剖面的模拟就显得不足。挡板对增加湍流度有明显的作用,通过 在中上高度添加挡板将有助于提高上部的湍流度,但是随着挡板数量的增加和挡板宽度的增加,它 们对风剖面的影响也增大,给流场调试带来较大的困难。
粗糙元 数量(排)
8
8
22
31
排列间距 纵向:50cm 横向:60cm 纵向:50cm 横向:60cm 纵向:40cm 横向:50cm 纵向:40cm 横向:50cm
流场模拟时采用图 1 中所示的各类型尖塔和不同的粗糙元,具体布置方案见表 2。粗糙元从转 盘的前沿开始向上游排列,纵横交错,如图 2 所示。按照 1/350 缩尺比模拟出四类地貌的平均风剖 面及湍流度剖面见图 3,对于湍流度剖面,同时给出了一些主要国家的风荷载规范的建议值用于比 较。限于篇幅,文中只给出了 C 类地貌下较低高度 32.6cm 处和梯度风高度以上 120.1cm 处的顺风向 脉动风谱,见图 4。
从图 3 可以看出在所模拟的结果中,四种地貌的平均风剖面都是较理想。在图 4 中,所模拟的 顺风向脉动风谱和理论顺风向脉动风谱进行了比较,可以看出不管是边界层较低高度的顺风向脉动 风谱还是边界层较高高度的顺风向功率谱都和 Karman 谱比较接近。从图 3 也可以看出,我国规范 的湍流度规定相对于世界上一些主要国家规范的湍流度要小得多。西方主要国家的风荷载规范对湍 流度有明文说明,而我国荷载规范并没有这样的说明,只是在风振分析中间接使用到,据此可以导 出我国规范中的湍流度表达式 4)。在汕头大学风洞中模拟的四种地貌的湍流度剖面中,C、D 类地貌 的湍流度剖面基本上处于各主要国家规范中的湍流度剖面之间。B 类地貌的湍流度剖面基本上处于 各主要国家规范中的湍流度剖面的下限,A 类地貌中等高度范围的湍流度剖面要比各主要国家规范 中的湍流度剖面小。
[3] 黄鹏,全涌,顾明,TJ-2 风洞大气边界层被动模拟方法的研究,同济大学学报,1999,27(2),pp136-140。 [4] 洪小健,顾明,顺风向等效风荷载及响应——主要国家建筑风荷载规范比较。建筑结构,2004,34(7),pp39-43。 [5] H.P.A.H. Irwin, The design of spires for wind simulation, J. Wind Eng. Ind. Aerodyn., 1981, 7, pp361–366. [6] Hee J. Ham, Bogusz Bienkiewicz, Wind tunnel simulation of TTU flow and building roof pressure, Wind Eng. Ind.
3. 模拟结果
图 2 风洞中模拟装置的布置图
表 2 四类地貌流场模拟的布置方案
地貌 类型 数量
A A类
尖塔
位置
B B类 C C类
风洞入口下游 0.8m 处 3 距离转盘 14.2m
尖塔之间间隔为 75cm
D D类
尺寸 迎风面 8cm×5cm 顺风向长 6.5cm 迎风面 8cm×5cm 顺风向长 6.5cm 迎风面 10cm×10cm 顺风向长 10cm 迎风面 10cm×10cm 顺风向长 10cm
Aerodyn., 1998, 77&78, pp119-133. [7] 庞加斌,林志兴,陆烨,关于风洞中用尖劈和粗糙元模拟大气边界层的讨论,流体力学实验与测量,2004,
18(2),pp32-37。 [8] 陈凯,毕卫涛,魏庆鼎,振动尖塔对风洞模拟大气湍流边界层的作用,空气动力学报,2003,21(2),pp211-217。 [9] 黄鹏,施宗城,陈伟,顾明,大气边界层风场模拟及测试技术的研究,同济大学学报,2001,29(1),pp40-44。 [10] 黄东群,马健,董力耘,徐有恒,大气边界层流场的模拟与大气边界层风洞,上海力学,1999,20(1),pp10-15。
H.J. Ham6)和庞加斌 7)分别采用了不同的非三角形尖劈,较好地模拟出大缩尺比下的流场特性, 满足了低矮建筑物和桥梁的抗风试验研究。针对风洞大气边界层上部湍流度较小的问题,陈凯 8)采 用了一种具有弹性底座的振动尖塔,改善了风洞大气边界层上部的湍流度,提高了湍流积分尺度, 文中模拟得到的湍流度剖面和风剖面都还有待进一步改善。本文使用的曲边梯形尖塔,其一定宽度 的上底边长保证了在边界层较高的高度能够保持较大的湍流度,尖塔宽度随着高度降低而缓慢增加, 不至于使得风洞底层的堵塞度过大而导致难于协调小地貌指数和高湍流度之间的矛盾。最终模拟得 到的边界层上部具有较高湍流度,整个湍流度分布较为合理。图 5 比较了汕头大学风洞中分别使用 三角形尖塔和曲边梯形尖塔模拟的 C 类地貌的湍流度剖面。表 3 给出了国内几个风洞模拟的 C 类地 貌不同高度湍流度的比较,文献[9]采用三角形尖塔和粗糙元进行边界层模拟,文献[10]在三角形尖 塔上还添加了 3 条格栅,使得边界层上部的湍流度要比文献[8]要高些,本文采用曲边梯形尖塔进行 边界层模拟,得到的湍流度剖面在较高的位置仍能有较大的湍流度。曲边梯形尖塔相对于三角形尖 塔的优势在于提高了边界层中部以上高度湍流度,使得湍流度剖面更为合理。曲边梯形尖塔结合挡 板,模拟出的 D 类地貌在风洞 60cm 高度甚至可以达到 26%的湍流度,为高湍流度下的群体干扰效 应研究提供了流场条件。
关键词 风洞、大气边界层、流场模拟
1. 引 言
大气边界层的风洞模拟是研究风对建筑结构的作用以及建筑物周围的风环境等风工程问题的基 础工作。在进行建筑结构的风洞试验时,除了必须满足建筑结构的几何特性相似外,还要准确地模 拟出建筑结构所在区域的大气边界层特性。大气边界层的特性主要包括平均风速剖面、湍流度剖面、 风功率谱以及湍流积分尺度等。尤其对于研究建筑结构风致响应以及气动弹性效应的风洞试验,准 确模拟湍流结构有着更重要的意义。
大气边界层风洞流场特性的模拟
石碧青 洪海波 谢壮宁 倪振华
(汕头大学土木系,广东汕头,515063)
摘 要 风场模拟的准确程度是保证风洞试验结果准确程度的重要因素。文中采用了一种曲边梯形的尖塔,结 合粗糙元等被动装置进行风洞大气边界层的模拟,模拟出试验对象所在区域的大气边界层特性,包括 A、B、C 和 D 类四种地貌下的平均风速剖面、湍流度剖面和顺风向脉动风谱,为进一步的流场特性研究和风洞试验奠定了基础。 通过对曲边梯形尖塔和三角形尖塔模拟的风洞大气边界层的湍流度剖面的比较,可以发现本文所采用的曲边梯形尖 塔更有利于提高风洞大气边界层中部以上高度的湍流度,使得整个湍流度剖面更为合理。