基于格子Boltzmann方法的固定方柱绕流研究

!第"#卷第$期郑州大学学报!理学版"%&’("#)&($!*#+,年-月./0123451&67389/!):;/<=8/>?/"
<2@/*#+,
收稿日期!*#+AB
+*B +C 基金项目!国家自然科学基金项目!*+EAE*"A "/
作者简介!靳遵龙!+-A$%"#男#河南郑州人#教授#主要从事换热设备强化传热及流体流动研究#>B F :8’&5’R 83G556/2?6/=3$通信作者&王永庆!+-A-%"#男#河南郑州人#副教授#主要从事流动与传热数值分析及过程传热强化与节能技术研究#>
B F :8’&]:34a HG556/2?6/=3/基于格子$
AF 4U H :==方法的固定方柱绕流研究靳遵龙!!王元凯!!王永庆
!郑州大学化工与能源学院!河南郑州C"###+"
摘要!研究了格子D
&’;5F :33方法的基本原理和边界条件处理方法#并利用该方法对不同雷诺数下固定单方柱绕流流场进行分析#探究了方柱产生卡门涡街的临界雷诺数和范围#验证了格子D
&’;5F :33方法边界处理和数值模拟的正确性和便捷性/对雷诺数为*##时并列双方柱不同分布间距的流场进行了模拟#结果表明#当柱间距为*倍方柱边长时#流体绕流方柱的涡流彼此影响最为明显/
关键词!方柱绕流$卡门涡街$格子D
&’;5F :33方法$数值模拟中图分类号!U
$"A(+文献标志码!M 文章编号!+EA+B E,C+!*#+,"#$B #++EB #"&’(!+#(+$A#"Q R /8I I 3/+EA+B E,C+(*#+A$E+
)*引言
实际生活中很多情形都可以归结为钝体绕流问题#如河流绕过桥墩的流动以及建筑物的风载等/流体流经钝体时#由于过流界面减小#会产生一系列复杂的物理现象#特别是边界层的分离进而在钝体后产生周期脱落的涡/目前钝体绕流研究主要以圆柱绕流为主#对于一般形状柱体的研究还更多的处于实验阶段/
对于复杂流体绕流#传统的K L O 方法是从宏观):982X B <;&[2I 方程出发#利用各种数值方法进行模拟/文献’+(分别使用有限差分法和离散涡方法计算了低雷诺数和高雷诺数下的方柱绕流$文献’*(通过有限元法求解了*个圆柱左右并排时和前后排列时的绕流/格子D &’;5F :33方法作为一种介观模拟方法#将宏观运动和微观运动的统计平均联系起来#在流体力学问题的数值模拟方面具有广阔的前景/格子D &’;5F :33方法不但可以方便地模拟复杂几何形状边界的流体流动#还可以模拟系统的时间演化#并且只需要简单的算法就可以在计算机上实现并行计算#具有更高的计算效率/本文利用格子D &’;5F :33方法对不同雷诺数下固定单方柱绕流流场进行分析#并对并列双方柱不同分布间距的流场进行了模拟#验证了格子D &’;5F :33方法边界处理和数值模拟的正确性和便捷性
/
图+!O *g -模型离散速度"./0+!O 8I =X 2;292’&=8;82I &W O *g -F
&?2’+*格子$AF 4U H :==方法
格子D &’;5F :33方法中的基本变量是格点的分布函数)-#)-为沿-方向的粒子分布函数/本文采用O *g -模型#-个离散速度分布如图+所示/其演化方程!不含外力项"为
)-!;/+-!,#,/!,"7)-!;#,"’7!,3’)-!;#,"7)2H
-!;
#,"(#式中&!,为时间步长$3为无量纲松弛时间#3’3#K !,$离散速度+-表示粒子的运动方向#可以表示为
+#’!###"#++#$#+*#C ’!e+##"#!##e+"#+"#E #+A #,’!e+
#+"#!+#e+"#其中&+为当地声速/局部平衡态分布函数)2H
-!;
#,"仅与当地的密度和动量有关#由熵增加原理推出#可以表
!第$期靳遵龙#等$基于格子D
&’;5F :33方法的固定方柱绕流研究示为
)2H
-
’%0-’+/+-*H !R M /!+-*H "**!*R M
*
/H **!R M (/;!H $"5!!由平衡态分布函数通过K 1:@F :3>3I [&4展开#可以导出模型的宏观密度和宏观速度分别为
%’!-
)-#H ’+
%!-
)-+-5!!格子D &’;5F :33方法在物理空间上将系统粒子运动分为迁移和碰撞两个相对独立的过程#使其具备很好的并行特性以及较强的复杂边界处理能力/
-*格子$AF 4U H :==方法的初始条件和边界格式处理
流体的运动总是在一定的初始条件和边界条件下进行的#因此必须设置相应的初始条件和边界条件/格子D &’;5F :33方法中的基本变量是分布函数/对于稳态或者准稳态流动#初始条件对最终计算结果影响不大#可以直接将初始分布函数设为其平衡态分布函数/实际问题中#边界条件往往基于宏观物理量#如何根据宏观量合理地确定分布函数是格子D &’;5F :33方法中的重要问题/常用的边界主要分为启发格式边界)动力学格式边界)外推格式边界以及复杂格式边界$根据实际问题又可分为速度边界和压力边界等/-,+*非平衡态外推格式边界
借鉴了传统计算流体力学方法!如有限差分法"中的边界处理方法#文献’$(在此基础上提出了非平衡
态外推格式#将边界节点;2上的分布函数分解为平衡态和非平衡态两部分#即)-#;2’)2H -#;2/)32
-#;2
#-’*#"#E 5平衡态函数中用临近流体点;)的密度%)代替壁面密度%%#而非平衡态部分则使用临近流体点;)的非平衡态
部分进行近似#即)32-#;2’)32-#;)’)-#;)7)32
-#;)
#-’*#"#E 5理论研究结果表明#非平衡态外推格式具有二阶精度/-,-*非平衡态反弹格式边界
实际应用中#边界处的速度分量或者压力分布往往是已知的#如管道流动等/文献’C (在+--A 年提出的非平衡态反弹格式则可以很好地应用于上述两种边界条件/对于速度边界条件#有
%%’++7H %’)#/)*/)C /!)$/)E /)A "(#)+’)$/*$%%H %
#)"’)A 7+*!)*7)C "/+E %%H %/+*%%0%
#),
’)E /+*!)*7)
C "/+E %%H %7+*%%0%
5!!通过非平衡态反弹格式可以直接计算垂直于速度边界上的$个未知平衡态分布函数#相似的方法也可以用于压力边界和其他边界条件/
3*固定方柱绕流的数值模拟
与分析
图-!单方柱绕流模型示意图"./0-!<=12F :;8=?8:4X :F &W =1:332’W ’&]
@:I ;:I 834’2I H6:X 2=a ’83?2X
3,+*单方柱绕流模拟
柱体绕流是一类经典的非定常流体力学问题#数值模拟是这类问题最有力的研究工具/大量理论研究结果表明#柱体绕流中对流场起决定性作用的参数是雷诺数!!$"5!$不同#黏性不可压缩流体会呈现不同的流体状态/本文采用单松弛!c D \S "模型#利用Koo 语言对单方柱绕流问题进行计算机模拟#研究不同!$对方柱绕流状态的影响/图*为单方柱绕流模型示意图/可以看出#模型长宽比为$w
+#入口采用非平衡态反弹格式速度边界#入口速A
++
郑州大学学报!理学版"第"#卷
度为均匀来流E _
#(+$出口采用非平衡态反弹格式压力边界!设出口密度为+(#"#方柱表面和管道上下壁面则采用非平衡态外推格式/此处#雷诺数!
$’E X
#其中方柱的特征尺寸X 取"个格子单位/为了避免管道上下壁面对尾流的影响#这里取管道直径!即垂直方向"为+##个格子单位/
图$为不同雷诺数下流场涡量云图/可以看出#当!$较小时#流场是定常的#方柱后面出现一对上下对称的尾涡/随着!$的增加#当!$_$,时#流场逐渐转变为非定常#尾涡也逐渐脱落#尾流区中形成周期性摆动和交错的漩涡#即出现卡门涡街/随着!$的进一步增加#当!$_*$,时#尾涡有紊乱的趋势/进一步的模拟显示#当!$大于*$,时#尾流区由层流状态过渡到紊流状态
/
图3!不同雷诺数下流场涡量云图
"./03!%&X ;8=8;a =&3;&6X I :;?8W W 2X 23;Y 2a 3&’?I 36F Z2X I
图C 显示了!
$_*##时阻力系数和升力系数随时间的变化曲线#当前平均阻力系数和最大升力系数分别为+("+和#($+#与文献’"(中实验研究的数据相比具有一致性#可见利用格子D &’;5F :33方法模拟是可行的
/
图7!!
$_*##时阻力系数和升力系数随时间的变化曲线"./07!P 8F 218I ;&X a &W ?X :4=&2W W 8=823;:3?’8W ;=&2W W 8=823;:;!$_*##
!$_$"及!$_$A 时的流场流线分析如图"所示#可以看出#上下对称的尾涡逐渐扩大#上涡有分离的趋势/因此#定常流失稳的临界!$在$A 左右#这与文献’E (中的结论是相符合的/
与理论研究及文献中圆柱绕流的结果对比后发现#方柱绕流整体上符合一般的钝体绕流趋势#并且由于方柱的特殊性#临界!$又与圆柱绕流不同/由以上分析可知#方柱绕流中卡门涡街现象在!$为C#‘*$,时较为明显/本文将采用!$_*##进行下一步的模拟研究/3,-*双方柱绕流模拟
实际生活中更多的是流体绕过多个柱体流动的情形#柱体的个数)形状)排列方式等均会影响流场的结构/本文研究双方柱并列模型下不同分布间距对流场的影响#控制的变量为柱间距G 与方柱特征尺寸X 的比值
,
++
!第$期靳遵龙#等$基于格子D
&’;5F :33方法的固定方柱绕流
研究图8!不同雷诺数下流场流线图
"./08!<;X 2:F ’832@:;;2X 3I :;?8W W 2X 23;Y 2a 3&’?I 36F Z2X
I
图<!并列双方柱绕流模型示意图"./0<!<=12F :;8=?8:4X :F &W =1:332’W ’&]
@:I ;@:X :’’2’I H6:X 2=a ’83?2X I
!GK X "#并列双方柱绕流模型示意图如图E 所示/
采用相同的边界条件#V
)W 方向格子数分别为$E#和+*##在!$_*##的条件下分别对并列方柱G Q X _+)G Q X _*)G Q X _$)G Q X _C 四种条件进行模拟#其涡量云图如图A 所示/
图,显示了并列方柱在G Q
X _*和G Q X _$时的流场流线图/从图A 和图,中可以看出#当分布间距较小!G Q
X _+"时#两柱形成的漩涡彼此影响较小#还保留着部分单柱绕流的流场特征$随着分布间距的增加#流场逐渐变得复杂#漩涡分布呈现不规则状态#当
GK X _*时#流体绕流方柱的涡流彼此影响最为明显$分布间距继续增加#当G Q X 大于$时#两个涡的相互影响越来越小#方柱下游形成两个反向同步脱落的涡街#各自涡形逐渐接近单方柱卡门涡街状态
/
图V !不同分布间距下流场涡量云图
"./0V !%&X ;8=8;a =&3;&6X I :;?8W W 2X 23;?8I ;X 8Z6;8&3?8I ;:3=2I
7*小结
通过对单方柱和并列双方柱的模拟研究#验证了格子D
&’;5F :33方法边界处理和数值模拟的正确性和便捷性/通过对不同!$下单方柱流场的模拟分析#探究了方柱产生卡门涡街的临界雷诺数!!$_$A "#并得到了产生明显卡门涡街的雷诺数范围!!$_C#‘*$,"/此外#雷诺数!
$_*##下并列双方柱不同分布间距-
++
郑州大学学报!理学版"第"#
卷
图W!不同分布间距下流场流线图
"./0W!<;X2:F’832@:;;2X3I:;?8W W2X23;?8I;X8Z6;8&3?8I;:3=2I
!G Q X_+)G Q X_*)G Q X_$)G Q X_C"的流场模拟结果表明#当柱间距为*倍方柱边长时#流体绕流方柱的涡流彼此影响最为明显/
参考文献!
’+(!U S M.N T MM/)6F2X8=:’I8F6’:;8&3&W W’&]:X&63?X2=;:346’:X=a’83?2X I’.(/.&6X3:’&W]83?2348322X834:3?83?6I;X8:’:2X&?a B 3:F8=I#+--##$$!+"&+A+V+,#/
’*(!T>)>\d N)N.Y#<M c P M Y ML#<N g7>N Y MKcY#2;:’/)6F2X8=:’I8F6’:;8&3&W W’&]83;2X W2X23=2Z2;]223;]&=8X=6’:X=a’83?2X I 83;:3?2F:3?I8?2B Za B I8?2:X X:342F23;I’.(/.&6X3:’&W W’68?I:3?I;X6=;6X2I#*##+#+"!*"&$*A V$"#/
’$(!\7U0c#0d>)\K\#<d N DK/)&3B2H68’8ZX86F2^;X:@&’:;8&3F2;1&?W&X92’&=8;a:3?@X2I I6X2Z&63?:X a=&3?8;8&3I83;12’:;B ;8=2D&’;5F:33F2;1&?’.(/K1832I2@1a I8=I D#*##*#++!C"&$EE V$AC/
’C(!0U7g<#d>ie/U3@X2I I6X2:3?92’&=8;aZ&63?:X a=&3?8;8&3I W&X;12’:;;8=2D&’;5F:33D\S F&?2’’.(/J1a I8=I&W W’68?I# +--A#-!E"&+"-+V+"-,/
’"(!<M d MMS#D N<fM<\#T7Y M c N O d M Y S/P1X22B?8F23I8&3:’I;6?a&W W’&]@:I;:I H6:X2=a’83?2X:;’&]Y2a3&’?I36F Z2X I ’.(/N3;2X3:;8&3:’R&6X3:’&W12:;:3?W’68?W’&]#*##$#*C!+"&"C V EE/
’E(!周云龙#邓冬#刁成东#等/方柱绕流和圆柱绕流旋涡脱落的数值模拟’K(Q Q中国工程热物理学会多相流学术会议/青岛# *##,&+V E/
24?E M A=4D5"F AJ%:@424:4.B2^?:>5!MF.=E5>@$:@5EA=
4D5N:44.B5$AF4U H:==C54D AE
.N)063’&34#fM)\e6:3[:8#fM)\e&34H834
!P+Q..*.)D Q$:-+J*A#R-#$$(-#R J#F A#$(R3#?Q$#R@Q.H E#-0$(1-,3#?Q$#R@Q.H C"###+#D Q-#J"
#;@4>:B4&P12Z:I8=@X83=8@’2:3?;12Z&63?:X a=&3?8;8&3I&W;12’:;;8=2D&’;5F:33F2;1&?]2X2I;6?82?/ f8;1;18I F2;1&?#;12W’&]W82’?&W W’&]:X&63?:I834’2I H6:X2=a’83?2X:;?8W W2X23;Y2a3&’?I36F Z2X I]:I I;6?82?#:3?;12=X8;8=:’Y2a3&’?I36F Z2X:3?X:342&W;12K:X F239&X;2^I;X22;@X&?6=2?Za;12I H6:X2 =a’83?2X]2X22^@’&X2?#I&;1:;;12=&X X2=;32I I:3?=&3923823=2&W;12Z&63?:X a@X&=2I I834:3?36F2X8=:’I8F6’:;8&3&W;12’:;;8=2D&’;5F:33F2;1&?]2X292X8W82?/P12W’&]W82’?&W@:X:’’2’I H6:X2=a’83?2X I]8;1 8W W2X23;?8I;X8Z6;8&3?8I;:3=2I]:I I8F6’:;2?]8;1Y2a3&’?I36F Z2X:;*##/P12X2I6’;I I1&]2?;1:;;12 83;2X@’:a:F&34;129&X;2^:X&63?;12I H6:X2=a’83?2X]:I;12F&I;&Z98&6I]123;12=a’83?2X I@:=834]:I *;8F2I;12’234;1&W;12I H6:X2=a’83?2X/
L5M J A>E@&W’&]@:I;I H6:X2=a’83?2X$K:X F239&X;2^I;X22;$’:;;8=2D&’;5F:33F2;1&?$36F2X8=:’I8F6B ’:;8&3
!责任编辑&孔!薇"
#*+。