建筑群风场fluent计算分析大作业Word下载.doc
《建筑群风场fluent计算分析大作业Word下载.doc》由会员分享，可在线阅读，更多相关《建筑群风场fluent计算分析大作业Word下载.doc（6页珍藏版）》请在冰点文库上搜索。
在大气边界层中，气流质点运动杂乱无章，气流流动表现为湍流状态。
湍流是由大小不同尺度的涡体叠合而成，对时间和空间都是非线性的随机运动的，因此使用湍流模型解决此问题。
在湍流模型中，基于Reynolds时均的RealizableK-ε模型能在整体上很好地反映出建筑物表面风压的变化趋势，模拟结果与试验值相差较小且计算效率高，所以选用该模型。
二、计算模型设置
1、边界条件选取
计算流域入流处采用FLUENT中的速度进口边界条件（velocity-inlet）。
边界条件用于定义在流动进口处的流动速度及相关其它标量型流动变量。
该边界条件适用于不可压缩流动，对于可压缩流动问题时会使得入口处的总温度或总压有一定的波动，导致非物理结果，所以可压缩流问题不适合采用速度进口边界条件。
本文为不可压缩流，可采用以velocity-inlet边界。
需对流动速度v、k和ε定义。
设建筑物所在地形为B类地形，其风场为B类风场，10m高度处、10min平均的基本风压为w0=0.35kPa，相应的标准高度处平均风速为u0＝23.7m/s；
则人口处的湍流强度I、湍流动能K和湍流耗散率ε的具体表达式如下
；
其中z，分别是流域中任意高度和对应的平均风速，z由模型底部开始算起。
为湍流积分尺度，采用经验公式。
平均风速剖面、湍动能，c和耗散率值采用Fluent提供的UDF编程与Fluent作接口实现。
出口采用完全发展出流边界条件（outflow）。
Outflow边界条件用于出流边界上的压力或速度都未知的情况，适用于出口处流动是完全发展的情况。
计算流域顶部和两侧采用对称边界条件（symmetry），适用于流动及传热场是对称的情况，即具有镜像对称特征。
在对称轴或者是对称面上既无质量交换，又无热量等物理量的交换，在垂直对称面的方向上物理量梯度为零。
可用来描述粘性流动中的自由滑移壁面。
建筑物表面和地面是固定不动，不发生移动的，故采用无滑移的壁面条件（wall），wall是用于限定fluid和solid区域的一种边界条件。
对于粘性流体，采用粘附条件，即认为壁面处流体速度与壁面该处的速度相同，无滑移壁面的速度为零，壁面处流体速度为零。
2、参考压力点选取
对于不可压缩流动，由于边界条件中不包括压力边界条件，必须设置一个参考压力位置。
在计算流域中选一个不受建筑模型影响，并且远离流域边界的位置作为参考压力位置，FLUENT在每次迭代后，都将整个压力场减去这个参考压力位置的压力值，从而使得所有的压力均按照参考压力位置的值来度量，压力绝对值常比相对值高几个数量级，若采用压力绝对值进行计算，则会导致计算存在较大的误差。
因此我们选择点（-50，100,100）为参考点。
3、求解参数的设置
选择合适的离散格式能加速收敛功能，对解的精度和求解的速度起着决定性作用。
一阶精度离散格式包含过多的数值耗散，对于定量分析时应该避免使用。
采用高精度二阶迎风格式对动量、湍动能和湍流耗散率等进行数值离散，采用高精度的离散方法进行离散能得到较高的精度，但容易导致数值计算的失稳和发散。
所以在计算中必须选择较小的松弛系数。
本文需要将风对建筑物的作用进行定量的分析，所以采用二阶迎风格式配合较小的松弛系数。
考虑风对建筑物的作用属于定常状态，由于SIMPLEC可以加快收敛速度，因而采用SIMPLEC算法，配合较小的欠松弛系数使用。
本文进行数值模拟时，经多次尝试后，压力欠松弛系数取0.2，其他项取0.5。
三、计算结果分析
1、平均风压分布
建筑群压力分布云图如下图所示
压力分布云图
（1）
压力分布云图
（2）
从图中可以看到建筑物表面最大正压力为670Pa，主要出现在建筑物的迎风面；
而在两矩形建筑物后侧的圆柱型建筑物两侧出现最大负压，为-2190Pa，超过了迎风面的最大压力。
取z=50平面的压力分布，如下图所示：
z=50平面的压力分布
可以看到距离圆柱形建筑表面负压大小和其与矩形建筑的距离有关，距离矩形建筑较近一侧负压更大。
由此可以看出，随着建筑物间距的减小，风载对建筑物的影响变大。
2、建筑物周围风速矢量和流线分布
仍以50m高度平面为例分析，其速度矢量图如下所示
z=50平面的速度矢量图
其对应的速度迹线图如下图所示
z=50平面的速度迹线图
可以看到在圆柱体建筑物后方无明显的漩涡现象，而在两矩形建筑物后方都有很明显的两个漩涡，且受后方圆柱体建筑影响导致漩涡并不对称，也同样与建筑物之间的距离有关。
在矩形建筑与圆柱形建筑之间气流速度大大增大，在相邻的区域存在着间隙流效应；
而由于矩形建筑与空气流动方向垂直，产生屏蔽效应，直接影响了其迎风面的平均压力值。
而这些漩涡与建筑物的屏蔽效应也直接影响了建筑物风载荷的分布。
这也是之前圆柱形建筑物出现较高负压的原因。
四、结论
通过仿真我们发现，该建筑群模型在建筑物缝隙处产生的负压要大于建筑物迎风面压力，因此在考虑迎风面建筑结构加强问题的同时,还应在建筑房屋时,考虑到建筑物间距的影响,不能使间距过小,否则建筑物将受风载影响过大而发生较大变形,进而影响到寿命。
附录：
UDF自编程序profile.c
/*profile.c*/
/*UDFforspecifyingsteady-statevelocityprofileboundarycondition*/
#include"
udf.h"
#definezo10.0
#defineuo23.7
#definezb5.0
#definezg450.0
#definetur0.31
#definecmu0.09
DEFINE_PROFILE（inlet_x_velocity,thread,index）
{
realx[ND_ND];
realz;
face_tf;
begin_f_loop（f,thread）
{
F_CENTROID（x,f,thread）;
z=x[2];
F_PROFILE（f,thread,index）=uo*pow（z/zo,0.16）;
}
end_f_loop（f,thread）
}
/*PROFILEFORKINETICENERGY*/
DEFINE_PROFILE（k_profile,thread,index）
realx[ND_ND];
realz,lr,um;
face_tf;
begin_f_loop（f,thread）
{
F_CENTROID（x,f,thread）;
z=x[2];
if（z<
zb）
lr=tur*pow（1.0,1.）;
else
lr=0.1*pow（z/zg,-0.25）;
um=uo*pow（z/zo,0.16）;
F_PROFILE（f,thread,index）=1.5*pow（um*lr,2.）;
/*PROFILEFORDISSIPATIONRATE*/
DEFINE_PROFILE（E_profile,thread,index）
realx[ND_ND];
realz,lr,um,l,k;
begin_f_loop（f,thread）
{
F_CENTROID（x,f,thread）;
z=x[2];
if（z<
lr=tur*pow（1.0,1.）;
else
lr=0.1*pow（z/zg,-0.25）;
um=uo*pow（z/zo,0.16）;
k=1.5*pow（um*lr,2.）;
l=100*pow（z/30,0.5）;
F_PROFILE（f,thread,index）=pow（cmu,3/4）*pow（k,3/2）/l;
}
end_f_loop（f,thread）
升级会员 