airpak-建筑群流场模拟中边界条件的设定
- 格式:doc
- 大小:222.50 KB
- 文档页数:14
边界条件定义边界条件概述边界条件包括流动变量和热变量在边界处的值。
它是FLUENT分析得很关键的一部分,设定边界条件必须小心谨慎。
边界条件的分类:进出口边界条件:压力、速度、质量进口、进风口、进气扇、压力出口、压力远场边界条件、质量出口、通风口、排气扇;壁面、repeating, and pole boundaries:壁面,对称,周期,轴;内部单元区域:流体、固体(多孔是一种流动区域类型) ;内部表面边界:风扇、散热器、多孔跳跃、壁面、内部。
(内部表面边界条件定义在单元表面,这意味着它们没有有限厚度,并提供了流场性质的每一步的变化。
这些边界条件用来补充描述排气扇、细孔薄膜以及散热器的物理模型。
内部表面区域的内部类型不需要你输入任何东西。
)下面一节将详细介绍上面所叙述边界条件,并详细介绍了它们的设定方法以及设定的具体合适条件。
周期性边界条件在本章中介绍,模拟完全发展的周期性流动将在周期性流动和热传导一章中介绍。
使用边界条件面板边界条件(Figure 1)对于特定边界允许你改变边界条件区域类型,并且打开其他的面板以设定每一区域的边界条件参数菜单:Define/Boundary Conditions...Figure 1: 边界条件面板改变边界区域类型设定任何边界条件之前,必须检查所有边界区域的区域类型,如有必要就作适当的修改。
比方说:如果你的网格是压力入口,但是你想要使用速度入口,你就要把压力入口改为速度入口之后再设定。
改变类型的步骤如下::1.在区域下拉列表中选定所要修改的区域2.在类型列表中选择正确的区域类型3.当问题提示菜单出现时,点击确认确认改变之后,区域类型将会改变,名字也将自动改变(如果初始名字时缺省的请参阅边界条件区域名字一节),设定区域边界条件的面板也将自动打开。
!注意:这个方法不能用于改变周期性类型,因为该边界类型已经存在了附加限制。
创建边界条件一节解释了如何创建和分开周期性区域。
2.6.边界条件2.6.1概述所有CFD问题都需要有边界条件,对于瞬态问题还需要有初始条件。
流场的解法不同,对边界条件和初始条件的处理方式也不一样。
所谓边界条件,是指在求解域的边界上所求解的变量或其一阶导数随地点及时间变化的规律。
只有给定了合理边界条件的问题,才可能计算得出流场的解。
因此,边界条件是使CFD 问题有定解的必要条件,任何一个CFD问题都不可能没有边界条件。
在CFD模拟时,常用的基本边界条件包括:.流动进口边界,.流动出口边界,.给定压力边界,壁面边界,.对称边界和.周期性(循环)边界2.6.2边界条件2.6.2.1流动进口、出口边界条件FLUENT提供了10种类型的流动进、出口条件,它们分别是:1.速度进口:给出进口速度及需要计算的所有标量值,适用于不可压缩流动2.压力进口:给出进口的总压和其它需要计算的标量进口值3.质量流进口:主要用于可压缩流动,给出进口的质量流量。
对于不可压缩流动,没有必要给出该边界条件,因为密度是常数,我们可以用速度进口条件。
4.压力出口:给定流动出口的静压。
对于有回流的出口,该边界条件比outflow 边界条件更容易收敛。
5.压力远场:该边界条件只对可压缩流动适合。
6.outflow:该边界条件用以模拟在求解问题之前,无法知道出口速度或者压力;出口流动符合完全发展条件,出口处,除了压力之外,其它参量梯度为零。
该边界条件不适合可压缩流动。
7.inlet vent:进口风扇条件需要给定一个损失系数,流动方向和环境总压和总温。
8.intake fan:进口风扇条件需要给定压降,流动方向和环境总压和总温。
9.out let vent:排出风扇给定损失系数和环境静压和静温。
10.exhaust fan.:排除风扇给定压降,环境静压。
2.6.2.2压力进口边界条件压力进口边界条件通常用于给出流体进口的压力和流动的其它标量参数,对计算可压和不可压问题都适合。
压力进口边界条件通常用于不知道进口流率或流动速度时候的流动,这类流动在工程中常见,如浮力驱动的流动问题。
计算流体力学中的边界条件处理在计算流体力学中,边界条件处理是一个至关重要的步骤。
边界条件是指在数值计算中,对于流场的边界处所设定的条件,用于模拟真实流动情况,并保证数值计算的准确性和可靠性。
本文将对计算流体力学中的边界条件处理进行综述,包括常见的边界条件类型和其在不同应用中的处理方式。
一、边界条件类型1. 进口边界条件进口边界条件是指流场的进口边界,即外部流体进入计算区域的边界。
在进口边界处需要设定流体的入口流速、温度、浓度等参数。
常用的进口边界条件有恒定流速、恒定温度和恒定浓度等。
进口边界条件的处理方式通常采用指定数值来模拟实际流动情况。
2. 出口边界条件出口边界条件是指流场的出口边界,即计算区域的外部流体离开的边界。
出口边界条件需要设定出口处的压力、速度等参数。
常见的出口边界条件有静压出口、出流出口等。
出口边界条件的处理方式主要是通过迭代计算来确定达到稳定状态的数值解。
3. 壁面边界条件壁面边界条件是指流场与实际物体接触的部分,需要考虑流体在壁面上的速度、温度等的变化。
通常情况下,流体在壁面上的速度是零,即无滑移边界条件;温度则可根据壁面材料的传热性质进行设定。
壁面边界条件的处理方式通常采用无滑移条件和指定壁面温度条件。
4. 对称边界条件对称边界条件是指流场的某个边界面对称分布的情况。
在对称边界处,流动的物理量具有对称分布的特点,例如速度分量、压力等。
对称边界条件的处理方式是将对称面上的物理量进行相等的设定,以模拟对称分布情况。
二、边界条件处理方式1. 插值法插值法是一种常用的边界条件处理方式。
通过在已知的边界节点上求解物理量的值,然后通过插值方法计算出其他边界节点上物理量的近似值。
插值法能够通过边界条件的已知值预测其他未知值,从而实现对流场的模拟和计算。
2. 外推法外推法是一种基于已知的数值求解方法,通过已知节点上的物理量值来预测边界处未知节点上的物理量。
外推法的基本思想是根据已知节点处的物理量值,利用数值计算方法来迭代求解其他未知边界节点上的值。
fluent压力出口边界条件设置FLUENT是一种基于CFD(计算流体动力学)数值计算的软件,可以用来模拟各种物理场景,如流体流动、热传导、化学反应等。
在模拟流体流动时,FLUENT可以通过设置边界条件来模拟各种不同的情况,其中较为重要的一个条件就是压力出口边界条件。
本文将从步骤、作用以及注意事项等方面来详细介绍如何设置fluent压力出口边界条件。
1.概述压力出口边界条件是FLUENT中十分重要的一个边界条件,它是用来规定在流体通过边界时的压力变化。
设置正确的压力出口边界条件可以使得模拟结果更加准确,这对于各种流体流动问题的研究都具有重要意义。
2.步骤(1)首先打开FLUENT软件,选择要进行计算的模型(通常是一个几何体),打开模型,在界面的BCs(边界条件)标签中找到压力出口。
(2)在压力出口边界条件中,需要设置出口的压力值。
根据实际问题可以选择不同的类型,如静态压力、总压力或者平均压力等。
(3)在压力出口的另外一个设置中,需要设置出流方向,通常FLUENT 会根据模型的几何形状自动识别出出流的方向,可根据实际情况进行调整。
(4)设置完毕后,需要点击Calculate按钮,并选择解算参数,包括数量、误差、最大迭代次数等参数,通常根据实际需要进行设置即可。
(5)最后,点击Run按钮,开始计算。
计算时间根据模型规模不同,可能需要几分钟到几个小时不等。
3.注意事项在设置压力出口边界条件时,需要注意以下几点:(1)在实际计算中,需要确定出口的实际压力值,可以通过实地测量、经验公式或者其他软件计算来确定。
(2)需要根据实际问题选择合适的边界条件类型以及合适的出流方向来进行设置。
(3)在计算过程中,需要注意迭代次数、计算结果的收敛情况以及计算结果的稳定性等问题。
(4)需要注意FLUENT软件的版本和计算参数对于计算结果的影响。
以上就是关于FLUENT压力出口边界条件设置的详细介绍,可以仔细阅读并按照步骤进行设置。
flow simulation边界条件在流体模拟中,边界条件是指在流场中确定流体流动特性的条件。
边界条件分为两大类:入口边界条件和出口边界条件。
入口边界条件用于确定流体进入流场的初始状态。
常见的入口边界条件有:1. 速度边界条件(Velocity boundary condition):指定流体进入流场的速度分布。
2. 压力边界条件(Pressure boundary condition):指定流体进入流场的压力分布。
3. 温度边界条件(Temperature boundary condition):指定流体进入流场的温度分布。
4. 流量边界条件(Mass flow rate boundary condition):指定流体进入流场的质量流量。
除了入口边界条件外,出口边界条件还用于确定流体离开流场时的特性。
常见的出口边界条件有:1. 压力出口边界条件(Pressure outlet boundary condition):指定出口处的压力。
2. 流量出口边界条件(Mass flow outlet boundary condition):指定出口处的流体质量流量。
3. 常压出口边界条件(Atmospheric pressure outlet boundary condition):将出口处视为常压条件。
其他常见的边界条件还包括:1. 对称边界条件(Symmetry boundary condition):假设流场对称,使流过对称边界时速度、压力等参数满足对称要求。
2. 壁面边界条件(Wall boundary condition):指定流场与固体墙壁接触时的流动状态,通常要考虑壁面摩擦对流体的影响。
3. 切向旋转平均边界条件(Tangential average boundary condition):常用于旋转机械中,用于对流场的切向部分进行平均处理。
4. 界面边界条件(Interface boundary condition):用于模拟不同相流体之间的界面,常用于多相流动模拟中。
fluent中的边界条件在Fluent中,边界条件是用来定义问题的边界和限制条件,以便进行数值模拟和求解。
边界条件对于模拟结果的准确性和可靠性至关重要。
下面我将从多个角度来回答关于Fluent中边界条件的问题。
1. 类型,Fluent提供了多种类型的边界条件,以适应不同的模拟需求。
常见的边界条件包括,速度入口边界条件、压力出口边界条件、壁面边界条件、对称边界条件等。
每种边界条件都有特定的物理意义和数学表达方式。
2. 物理意义,边界条件反映了流体在模拟过程中与模拟区域边界的相互作用。
例如,速度入口边界条件用于指定流体从哪个方向进入模拟区域,压力出口边界条件用于指定流体从模拟区域中的哪个位置流出。
壁面边界条件用于模拟流体与实际物体表面的相互作用。
3. 数学表达,每种边界条件在Fluent中都有相应的数学表达方式。
例如,速度入口边界条件可以通过指定流体的速度分量来定义,压力出口边界条件可以通过指定出口处的压力值来定义。
壁面边界条件可以通过指定表面的摩擦系数或温度来定义。
4. 设置方法,在Fluent中,设置边界条件可以通过图形界面或者命令行界面来完成。
在图形界面中,用户可以通过选择相应的边界条件类型,并输入相应的参数值来设置边界条件。
在命令行界面中,用户可以使用相应的命令来设置边界条件。
5. 边界条件的影响,边界条件的设置对模拟结果有着重要的影响。
合理选择和设置边界条件可以保证模拟结果的准确性和可靠性。
不恰当的边界条件设置可能导致模拟结果的偏差或不收敛。
总结起来,Fluent中的边界条件是用来定义问题边界和限制条件的重要参数。
合理选择和设置边界条件对于模拟结果的准确性和可靠性至关重要。
在设置边界条件时,需要考虑物理意义、数学表达和设置方法等因素,并根据具体模拟需求进行选择和调整。
Airpak对教室热环境的数值模拟研究作者:郑万冬苟湘来源:《消费导刊·理论版》2008年第23期[摘要]教室为重要学习场所,其环境质量的优劣直接影响学习效率和授课质量。
本文利用Airpak软件对教室空调系统的温度分布、气流组织、空气龄、热舒适性等情况进行模拟研究,得出结果为教室热环境舒适性的改善提供参考。
[关键词]Airpak 数值模拟热舒适性一、引言有学者应用Airpak软件对办公室、客车等空间的气流组织、热舒适性等情况进行预测,并实践证实了Airpak软件模拟效果良好[1-4]。
二、模型的建立和边界条件的确立(一)物理模型选用的具有典型特征的阶梯教室(长×宽×高=11m×7.5m×4m),如图1所示,有一面外墙,外墙上有6个窗户,两扇内门,教室座位144个,一个电脑桌,一个投影仪,照明灯具28只,6个送风口,2个回风口。
室内人员、电脑、投影仪、灯为室内负荷,均按实际运行计算。
(二)数学模型及网格划分为了简化问题,作如下假设[5]:(1)室内气体低速流动,可视为不可压缩流体,且符合Bossinesq假设。
(2)流动为稳态流动。
(3)室内空气为辐射透明介质。
(三)边界条件1.入口边界条件:送风温度为18℃,送风方式为顶部垂直送风,风速为3m/s。
2.出口边界条件:回风方式为顶部回风,压力为室内环境压力。
3.壁面边界条件:教室有一面外墙,窗户6个,门2扇,均按实际情况计算负荷。
其余三面墙、地面、房顶设为绝热边界。
三、数值模拟结果及分析(一)温度图2和图3分别是地面1m和2m高度的截面温度分布,可以看出,正对送风口下方的温度较低,阶梯教室前面的温度较低,这与阶梯教室后面较高、人员较集中,以及前面教室空间较大有关。
但总的看来,在同学座位附近温度较均匀,基本维持在22℃~26℃。
(二)速度图4和图5分别给出了教室送风口和回风口处截面速度分布,在正对送风口下方的座位,送风速度超过1m/s,但在其它大部分座位,速度在0.35m/s以下,即在正对送方口下方的座位,会有一定的不舒适吹风感,其它座位因风速小而吹风感不明显。
fluent 翼型边界条件设置
在 Fluent 中,设置翼型的边界条件涉及以下几个步骤:
1. 导入翼型的几何模型:在 Fluent 的“File”菜单中选择“Import”->“Geometry”,然后选择翼型的几何模型文件进行导入。
2. 定义边界条件:在 Fluent 的“Define”菜单中选择“Boundary Conditions”。
选择翼型表面上的边界,例如翼型的上表面和下表面,然后给它们分配适当的边界条件。
常见的翼型边界条件有:
- 装壁面:对于翼型的表面,可以选择“Wall”作为边界条件,并指定壁面的摩擦系数。
这样可以模拟气体在壁面处的粘性效应。
- 远场:对于翼型周围的远场区域,可以选择“Far Field”作为边界条件,并指定远场参数,例如空气的压力和速度。
- 入流:对于翼型前端(进气流动方向的一侧),可以选择“Inlet”作为边界条件,并指定入流的参数,例如入流速度、入流温度和入流湍流参数。
- 出流:对于翼型后端(气流流动方向的一侧),可以选择“Outlet”作为边界条件。
可以根据实际情况指定出流的压力、速度和湍流参数。
3. 设置求解器参数:在 Fluent 的“Solver”菜单中选择“Solver Settings”。
根据需要调整求解器的参数,例如迭代次数、收敛准则等。
4. 运行求解器:在 Fluent 的主界面上点击“Calculate”按钮,开始运行求解器进行计算。
通过以上步骤,可以在 Fluent 中设置翼型的边界条件,并进行流场计算和分析。
具体的设置方法还可以根据实际情况进行调整和优化。
Airpak 2.0 使用指南办公室通风•简介•前言•问题描述•第一步:打开并定义一个新任务•第二步:创建模型•第三步:加入辐射和模型•第四步:计算网格•第五步:检查风状况•第五步:将模型保存到任务文件中•第六步:计算一个解决方案•第七步:检查结果•汇总•参考文献简介这篇教程演示了如何使用Airpak在计算机上创建有两个人工作的办公室的模型教程中你将学会:•打开一个新任务•包括计算机模拟相对湿度影响的分布•创建blocks, openings, vents, partitions, 和 walls•创建辐射影响的模型•改变计算机计算的迭代次数•计算一个结果•检查模型表面轮廓和物体表面的矢量,以及横断面图•空气入口的粒子流轨迹•通过计算机计算检查房间里的舒适度的预测平均值(PMV)和不满百分比(PPD)必装条件这个教程假设你对airpak没有经验,但是你很熟悉接口,如果您不是,请检查第一张的用户指南问题描述这个办公室分为两个部分,每部分包含一个人,坐在电脑前工作,如图 1.1。
这个办公室还包括六个荧光灯、暖气片、空气入口、回风口,窗户。
表面温度和空气流速剖面,目的是为了确定整个舒适的房间。
Table 1.1:几何参数,导热性,窗户和扩散器的流动边界条件尺寸温度速度13.5 C 0.85 m/s进气口扩散0.2 m 0.3 m30.9 C ---窗户3.65 m 1.16 mTable 1.2:热源的尺寸和性能热源Size Power1500 W墙挂式取暖1.2 m 0.1 m 0.2 m75 W人0.4 m 0.35 m 1.1 m108 WComputer 10.4 m 0.4 m 0.4 m173 WComputer 20.4 m 0.4 m 0.4 m34 W灯0.2 m 1.2 m 0.15 mFigure 1.1:通风问题所研究的办公室的几何参数Step 1: Open and Define a New Job1.启动airpak,如简介1.5中说明的那样,打开任务面板将自动开启2.为你的任务取一个名字(a)In the Select the job to open text box, type /office at the end of the path.在选择任务打开文字框内,输入/office在路径的最后(b)Click Accept. 点击接受Airpak will create a default room with the dimensions 10 m 3 m 10 m, anddisplay the room in the graphics window.Airpak将创建一个默认的尺寸为10x3x10的房间,并放置在图形窗口内。
Airpak 2.0 使用指南办公室通风∙简介∙前言∙问题描述∙第一步:打开并定义一个新任务∙第二步:创建模型∙第三步:加入辐射和模型∙第四步:计算网格∙第五步:检查风状况∙第五步:将模型保存到任务文件中∙第六步:计算一个解决方案∙第七步:检查结果∙汇总∙参考文献简介这篇教程演示了如何使用Airpak在计算机上创建有两个人工作的办公室的模型教程中你将学会:∙打开一个新任务∙包括计算机模拟相对湿度影响的分布∙创建blocks, openings, vents, partitions, 和 walls∙创建辐射影响的模型∙改变计算机计算的迭代次数∙计算一个结果∙检查模型表面轮廓和物体表面的矢量,以及横断面图∙空气入口的粒子流轨迹∙通过计算机计算检查房间里的舒适度的预测平均值(PMV)和不满百分比(PPD)必装条件这个教程假设你对airpak没有经验,但是你很熟悉接口,如果您不是,请检查第一张的用户指南问题描述这个办公室分为两个部分,每部分包含一个人,坐在电脑前工作,如图 1.1。
这个办公室还包括六个荧光灯、暖气片、空气入口、回风口,窗户。
表面温度和空气流速剖面,目的是为了确定整个舒适的房间。
Table 1.1:几何参数,导热性,窗户和扩散器的流动边界条件尺寸温度速度13.5 C 0.85 m/s进气口扩散0.2 m 0.3 m30.9 C ---窗户3.65 m 1.16 mTable 1.2:热源的尺寸和性能热源Size Power1500 W墙挂式取暖1.2 m 0.1 m 0.2 m75 W人0.4 m 0.35 m 1.1 m108 WComputer 10.4 m 0.4 m 0.4 m173 WComputer 20.4 m 0.4 m 0.4 m34 W灯0.2 m 1.2 m 0.15 mFigure 1.1:通风问题所研究的办公室的几何参数Step 1: Open and Define a New Job1.启动airpak,如简介1.5中说明的那样,打开任务面板将自动开启2.为你的任务取一个名字(a)In the Select the job to open text box, type /office at the end of the path.在选择任务打开文字框内,输入/office在路径的最后(b)Click Accept. 点击接受Airpak will create a default room with the dimensions 10 m 3 m 10 m, anddisplay the room in the graphics window.Airpak将创建一个默认的尺寸为10x3x10的房间,并放置在图形窗口内。
建筑群流场模拟中边界条件的设定摘要:CFD(Computational Fluid Dynamic)模拟技术近几年来随着计算机性能的大幅提高和模拟技术的进步,越来越得到广泛重视和应用。
边界条件的设定是直接关系到模拟成败的关键。
本文对边界设定中几个常见的问题分别进行了分析。
计算区域根据实际情况选取,水平方向参考取10倍最高建筑的高度,高度方向取3-4倍建筑物高度。
计算区域的顶面和侧面对于对称边界,速度边界和壁面这3种边界可以任意选择。
出流面可设为完全发展出流边界或速度边界。
对于大气层结中性情况,给出了风速廓线的计算公式,并介绍了一种为节约计算能力而分阶段划分网格的方法。
关键词:边界条件建筑外流场数值模拟1 引言风是构成室外环境的重要因素。
建筑群造成的流场变化会对其周围环境会产生一定的影响,特别是高层建筑更容易引起建筑周围行人高度局部风速的增大,产生不良的风环境,造成人的不适,降低环境质量。
如果在建筑规划设计阶段进行模拟研究,问题就可以得到预防和解决。
另外流场模拟对于结构风荷载以及自然通风等领域的研究也有价值。
风洞模型试验和实测是以往建筑流场研究的主要手段,但都费时费力。
CFD(Computational Fluid Dynamic)模拟技术近几年来随着计算机性能的大幅提高和模拟技术的进步,越来越得到重视和广泛应用,并取得了不少的成果。
作为有力的工具,数值模拟也已被用于建筑内外的流场研究中。
在模拟的过程中,除了湍流模型的选择,网格的划分等因素,边界条件也是直接关系到模拟成败的关键环节。
目前,数值模拟技术还处于发展阶段,只有不断与风洞试验和实测的结果进行比较,总结经验,对各方面进行改进,才能使之更加成熟完善。
本文结合模拟的例子,对边界设定中几个常见的问题分别进行了分析讨论。
2 模拟例子的情况简介本文通过结合一个模拟的例子来说明一些有关边界条件方面的问题。
例子是模拟上海某繁华商业街上一个供行人休息的小广场周围的行人高度(1.5至2米)的流场,结合日照等因素综合研究建筑对人及人的行为的影响。
小广场周围的大楼高度都在100米左右,形体及其造成的流场情况复杂。
建筑分布图如图1所示:图1 广场周围建筑平面图城市中下垫面结构复杂。
因在同一高度上增温冷却情况不均一,往往引起热力湍流,在高温区有上升气流,在低温区有下沉气流。
热力环流只在盛行风速极小或静风时才表现出来[文献1]。
因此,在本模拟中可忽略辐射及浮升力的影响。
由于建筑的存在而改变了流动结构,流动不能充分发展。
标准k-e模型对于某些情况过大估计了k值,但对于速度场,标准k-e模型与风洞试验结果较为接近[文献2,3],此外其波动小,有较高的精度和效率,在低速湍流工程上使用广泛。
本文选用标准k-e湍流模型,近壁面则采用标准壁面函数法计算。
3 梯度风的边界模拟任一时刻风速都是随高度变化的。
在摩擦层中,对于风随高度的变化规律,主要考虑摩擦力的影响。
风吹过地面时,受到地面上各种粗糙元(草、庄稼、森林、房屋建筑等等)产生的摩擦阻力作用而使风的能量减少即风速减小,这一层受地球表现摩擦阻力影响的大气层称为大气边界层,它的厚度约为1千米,建筑物模拟的计算区域一般不会超出这个高度。
对于垂直方向风速分布的描述常用以下两类公式:3.1 指数律3.1.1 假定大气边界层风速剖面符合简单的幂指数分布规律,指数a在梯度高度δ内保持不变,而δ本身只是a的函数,即:(1)式中Vo为在Zo高度的平均风速。
参考高度Zo一般采用10 米。
Vo即取10米高度的平均风速。
当Z大于某一高度ZH 时,风速趋于均匀V=V(ZH),ZH就是地转风出现的高度。
不同的地面条件,幂指数a不同,一般来说,a在0.14-0.40内取值。
3.1.2 指数a的选取对于指数a的值各研究者提供的值差异很大。
例如,波锐辛科等曾研究过在近地面100米高度以下,指数a与粗糙长度Zo(计算方法见下文公式4,5)具有线性关系:a=0.12Zo+0.18 (2)M.Hussain [文献4] 指出指数也可以用a=Ae-bv式得出,式中A,b是由当地实际条件取得的经验值。
结果与当地实测结果相符。
我国“建筑结构载荷规范”GBJ 9-87中不同类型地表面下的a值与梯度风高度(即大气速度边界层厚度)的关系如表1所示:地面类型适用区域指数a梯度风高度A近海地区,湖岸,沙漠地区0.12300mB田野,丘陵及中小城市,大城市郊区0.16350mC有密集建筑的大城市区0.2400m3.2 对数律3.2.1 对数律是由理论推导得出的。
具体推导见[文献5] :(3)式中表示在地转风高度以下,摩擦层内Z高度出的平均风速;V*表示摩擦速度;k表示卡曼常数,一般取0.4;Z表示高度;Z表示粗糙度长度(在中性稳定平衡下,垂直风速廓线外延到风速等于零的高度)。
上式的适用范围可由0.02 V*/f确定,强风时在200米的范围内能满意地表达风速分布。
粗糙度长度Z随地点的不同,相差很大,从零点零几到几米不等。
介绍两种常用的经验公式:1)(4)式中V1,V2分别为Z1,Z2高度的平均风速。
只要测出两高度及两个相应的风速,就能算出粗糙度长度Z0。
2)(5)式中H表示障碍物高度(米);a表示气流接触到的障碍物侧面积(米2);A表示障碍物覆盖面积。
常见的地形Z值如下表2所示:海面海滨矮草地,干旷草原田野休耕地森林0.00 10.0050.01-0.040.050.02-0.032.松树林果树林稀疏建成市郊密集建成市郊大都市中心0.9-1.0 1.00.2-0.40.8-1.22.0-3.03.2.2 对数律存在的问题从公式3可以看出,对数律在近地面小于粗糙长度Z时求得的风速为负值,这显然是不合理的。
对于研究建筑风荷载,这点问题尚可以忽略,但对于行人高度的流场情况研究,这是无法忍受的。
纽陶就Ekman边界层内风速对数廓线形式提出以下公式:(6)式中表示在Z高度的平均风速;VG 表示地转风速,VG=,式中为空气密度;f为科氏力参数,可根据所在地的纬度查得;为该地的水平压力梯度,从气象台可查得;c表示阻力系数,c阻力系数与无量纲的罗斯贝数(Rossy Number Ro)有关。
Csanady给出其关系式为:经观测和模拟试验证实,上述公式符合实况。
与公式3相比,公式6避免了风速出现负值的情况。
笔者认为在有条件取得地转风的大小的情况下,可以用上式代替对数律使用。
3.3 对数律与指数律的比较和共同存在的问题目前,对于近地面风速的描述到底是用对数律[文献6]还是指数律[文献7]尚无定论,两种公式都有实测支持。
还有待微气象方面的学者进行进一步的观测结果来验证。
目前大部分使用CFD工具的学者倾向于使用指数律。
值得注意的是对于不同的大气层结情况,风速廓线是不同的。
只有在中性时,才近似符合对数律或指数律。
因此对于夜晚等大气层结非中性情况,使用指数律或对数律描述风速会有很大的偏差。
平均风速与对数廓线的偏差如表3所示:情况1情况2情况3高度502001515不稳定层结-1%-4%-4%≈-22%稳定层结1%4%2%≈12%注:情况1,2,3,分别为10米高度时均风速约25m/s,12m/s,5m/s时的情况。
表3数据表明情况3偏差达10%以上。
需根据实际情况进行修正。
但是迄今为止,尚无可用于任意层结条件下的边界层风场的公式,虽有些工作计入了层结,但对层结处理过于简单且主要用于理论研究。
3.4 本例中的设置本例根据上海市气象资料冬季主导风向为北风和西北风,平均风速为3m/s。
对于舒适性研究风速不应只用平均风速,这样的结果偏保险。
最常见的最大风速等于平均风速的1.1倍[文献8]。
采用指数律,考虑到地处市中心和前方的绿地,指数a值折衷取0.35,风速廓线为:,风向西北。
4 计算区域边界条件的设定4.1 顶面和侧面的边界设定在无法对边界进行实测的情况下,可参考风洞模型试验的做法,将建筑群完全浸没在风里。
整个计算区域就如同一个大风洞,计算区域的顶面和侧面相当于风洞的壁面。
笔者在查阅一些文献后发现,顶面和侧面有3种边界可以选择:对称边界,速度边界和壁面。
这三种边界都有人使用,有将顶面和侧面都设成三者之一的,也有分别使用其中两种的。
其中速度边界设置为k,ε和速度三项的梯度为零,法向速度为零。
笔者分别将边界设成以上三种情况进行了模拟,三种情况1.5米高度截面速度矢量图如图2,图3,图4所示:?图2 使用对称面边界图3 使用速度界面边界图4 使用壁面边界从结果可以看出,除了靠外的一些部分流场有所不同,在内部主要关心的区域计算结果是基本相同的。
以上三种边界对内部流场的影响不大,即对最后结果影响不大。
虽然三者形式不同,共同的特点是通过对称平面没有扩散流量。
在建筑与计算区域各面之间有足够距离的情况下,三者是等效的,任取一种都可以。
4.2 入口面的边界设定对入口边界的设定没有异议,笔者查阅的文献中都设为速度边界。
但是入口湍流条件的设定意见并不统一。
4.2.1 描述大气湍流的最简单的参数是湍流强度。
[文献5]给出如下的湍流强度剖面:式中V*表示摩擦速度;V(z)表示Z高度的平均速度;一般假设β不随高度变化,对于不同粗糙长度的β值Z0.0050.070.31 2.5β 6.56 5.25 4.8544.2.2 [文献9]根据上海和广州的实测提出如下的湍流强度剖面。
4.2.3 [文献10] 给出如下的湍流条件:式中V*表示摩擦速度;Z表示高度;常数Cμ=0.09;卡曼常数k=0.4。
笔者将几种条件带入,对计算结果进行比较,发现入口湍流条件,对结果影响不大。
4.3 出流面的边界设定对于出流面的设定,可设成完全发展出流边界,即各项的梯度都为0,这就要求在建筑群的下游有足够的长度;或设成与入口相同的速度边界,当然是基于风速在流过建筑后不变的前提,水平剪切使边界层内流动减速,因而边界层增厚,但在大气边界层内,在梯度高度之下,水平压力梯度只是部分地被科里奥利力平衡,它将再赋能量给流体,从而抵消边界层的增长,因此流动的水平均匀性就能够保持[文献5 ]。
笔者试验结果是两者对结果的影响不大,使用两种边界的速度等值线图如图5,图6所示。
可任选其一。
图5 使用充分发展边界的剖面图图6 使用速度边界的剖面图4.4 地面的边界设定路面和建筑的墙面可以设成壁面。
然而除了建筑本身,周围栽种的植物尤其是枝叶茂盛的树木对近地面流场风速的阻碍作用很大,种植枝叶茂盛的行道树的人行道的风速比街道中心风速减低20-30%;公园浓荫下的风速为街道中心风速的50%[文献1]。
如不进行修正,风速结果可能偏大,对于行人高度的舒适性研究这点事不能忽略的。
由于植物形态的复杂性和多样性,对其进行建模是很困难的。
[文献11] 对树木对流场影响的模型进行研究,使用图7 所示薄片组成的树木模型。