图7为不同C风风率下燃尽率的变化曲线。从图中可以看到，随着C风风率的减小燃尽率不断降低。这是因为随着C风风率的减小，炉内回流区减小，拱上气流下冲深度减小，部分煤粉停留时间变短，而煤粉的着火距离变长，从而使得燃尽率降低。另一方面由于C风风率的减小，使得空气分级程度增加，因而燃尽率降低。
4结论
本文利用数值模拟的方法，研究了某低N O x燃烧新系统W 火焰锅炉的C风风率对燃烧特性及N O排放特性的影响。得到的主要结论有:
4.1随着C风风率的减小，对煤粉气流的托举作用减弱，拱上气流下冲深度减小，炉内燃烧剧烈程度减弱使得温度水平降低。
4.2随着C风风率的减小，空气分级程度增加，主燃烧区的氧含量降低，还原性气氛增强；且炉内温度水平降低，均有利于降低N O排放量。
4.3C风风率对煤粉燃尽率有较大的影响；随着C风风率的降低，炉内回流区减小，部分煤粉停留时间变短，而煤粉的着火距离变长，使得煤粉燃尽率不断降低。
参考文献
[1]任枫.FW型W火焰锅炉高效低NO x燃烧技术研究[D].哈尔滨：哈尔滨工业大学,2010.
[2]赵斯楠，方庆艳，马仑，陈刚.燃烧初期化学当量比对锅炉NOx 生成与排放特性的影响[J].燃烧科学与技术,2017,23(03):236-241.
[3]Ma L,Fang Q,Tan P,et al.Effect of the parated overfire air location on the combustion optimization and NOx reduction of a600MWe FW down-fired utility boiler with a novel combustion system[J].Applied Energy,2016,180:104-115.
[4]马仑，方庆艳，张成，陈刚，吕当振，段学农.600MW W型火焰锅炉拱上二次风低NO x燃烧特性的数值模拟及优化[J].燃烧科学与技术,2016,22(01):64-70.
[5]周安鹂.W火焰锅炉无烟煤掺烧煤泥的试验与数值模拟[D].武汉：华中科技大学,2019.
[6]吕当振，马仑，段学农，方庆艳.600MW亚临界W型火焰锅炉低氮燃烧特性数值模拟[J].热能动力工程,2015,30(04):598-604+ 654-655.
作者简介：周安鹂(1993,4-)，女，籍贯：湖北襄阳，硕士，助教，研究方向：电力生产技术、节能减排技术、电气自动化。
含导热块封闭方腔自然对流格子玻尔兹曼模拟研究Lat t i ce B ol t zm ann s i m ul at i on of nat ural convect i on i n a cl os ed s quare cavi t y w i t h heat conduct i ng bl ock
郑竟波杨伟林陈磐*韩奎廉正阳倪虎
(南京工程学院能源与动力工程学院，江苏南京211167)
封闭方腔自然对流现象广泛存在于太阳能集热器的集热过程、温室型太阳房的设计、电子元件的散热、双层玻璃窗的生产等[1-3]工业领域中，是传热与流动学科中的一个经典问题，具有广阔的实际应用背景。在方腔自然对流这一问题
摘要：采用格子Boltzmann方法，构建了内含导热块方腔自然对流的数学模型，进行了模型验证，并重点分析了瑞利数、导热块大小及位置对腔体流线分布和等温线分布的影响规律。结果表明，随着瑞利数的增加，其流场和温度场都会随之变化。但导热块大小变化时，其温度场变化不大。导热块位置对腔内流场温度场的影响很大，导热块位置偏上或者偏下，其涡产生的位置明显不同，腔内温度分布在导热块区域变为倾斜分布。
关键词：自然对流；瑞利数；导热块；格子Boltzmann模拟
Abstract:The lattice Boltzmann method is ud to construct a mathematical model of natural convecti
on in a square cavity with heat conduction blocks.The developed code is validated against published works by numerical simulation for natural convection.The influence of Rayleigh number,size and location of heat conduction blocks on the streamline and isotherm distribution of the cavity is analyzed.The results show that the flow field and temperature field will change with the increa of Rayleigh number.However,when the size of the heat conducting block changes,the temperature field changes little.The position of the heat conduction block has a great influence on the temperature field of the flow field in the cavity.The position of the heat conduction block is higher or lower,and the position of the vortex is obviously different.The temperature distribution in the cavity becomes inclined in the heat conduction block area.
Key words:Natural Convection;Rayleigh Number;Heat conducting block;Lattice Boltzmann Method
中图分类号:TK12文献标识码:A文章编号:2096-4390(2021)17-0041-05
(转下页)
中，速度场和温度场相互耦合，流动情况比较复杂，可呈现不同
的流动结构，因此，对方腔内部的流动及传热问题的研究一直是各国学者普遍关注的研究热点。
至今为止已有不少研究者致力于研究封闭腔体内的自然对流特性，由于其影响因素繁多，机理复杂，故研究者们都力求一种更为高效的数值计算方法。基于分子动理论的格子玻尔兹曼方法[4,5]由于其清晰的物理背景、高效的并行效率及算法简单等优异特征，近年来受到越来越多的关注，并在自然对流流动换热领域得到了很好的应用。如王立海[6]采用格子Bol t z m ann 方法研究了方腔内自然对流换热问题，并对该换热现象进行了理论分析。董长青[7]以压力分布函数和内能密度分布函数为演化变量，构建了新的双分布热格子Bol t z m ann 模型，对封闭方腔内自然对流进行了数值研究。曹先齐等[8]采用格子Bol t z m ann 方法构建了非稳态方腔自然对流的数学模型，并进行模型验证，重点讨论了瑞利数Ｒa 、倾斜角度θ对方腔内非稳态流动特性和传热特性的影响。文献[9]采用格子Bol t z m ann 方法，分析研究侧壁面为正弦温度分布条件下的方腔自然对流，并研究了瑞利数以及高宽比变化时其对腔体流线分布和等温线分布的影响规律。文献[10]采用格子Bol t z m ann 方法对倾斜二维封闭方腔自然对流换热进行研究，重点讨论了倾斜角对自然对流换热的影响。文献[11]用格子Bol t z m ann 方法对方腔内二维多孔介质自然对流传热问题进行了研究，并详细讨论了孔隙度及孔隙度变化时其对自然对流传热的影响。文献[12]另辟蹊径，应用格子Bol t z m ann 方法(LBM )研究了三角形腔内的水的自然流动问题，分析了不同瑞利数和三角形方腔倾角对流场和温度场的影响规律。文献[13]采用格子Bol t z m ann 方法研究了多孔介质方腔底部局部加热的自然对流传热，主要分析了高温热源位置及尺寸对流场、温度场的影响规律。
然相比于封闭方腔自然对流问题，实际工程应用中的方腔结构更为复杂，方腔中含有导热块更是非常普遍的情况，因此，本文以内含导热体的方腔自然对流为研究对象，采用格子玻尔兹曼方法，以探讨瑞利数、导热块大小及位置等因素对方腔自然对流流动规律的影响，分析自然对流过程中的流场特性和温度分布规律，以期获得对带有导热块的方腔自然对流研究及相似领域研究提供有益参考。
1问题描述及控制方程
方腔内含导热块自然对流的物理模型如图1所示。腔体的水平宽度为L ，垂直高度为H ，腔体高宽比K =H /L=1.0。腔体左侧壁面为高温壁面，温度为T h ，右侧壁面为低温壁面，温度为T c ，上下壁面绝热边界。在腔体内部正中心处有一正方形导热固体块，其宽度为W 。模拟时假设导热块的热扩散系数与腔内介质热扩散系数相同。
图1自然对流示意图
宏观控制方程如下：荦·u=0(1)(2)
纸醉金迷的意思(3)
(4)方程(1-4)中u 、T 、ρ、p 、ｖ、λ、c p 、β、T r e f 、t 分别为宏观速度、温度、密度、压力、运动粘性系数、导热系数、比热、热膨胀系数、参考温度、时间。
控制输运过程的无量纲准则数是普朗特数Pr 和瑞利数R a ，分别为:
(5)式中ａ、ΔT 为热扩散系数、特征温度差。2格子B ol t zm ann 方程模型采用基于Bous s i nes q 假设的耦合双分布函数模型来模拟封闭方腔自然对流[14]，速度场和温度场的耦合通过速度分布函数演化方程中增加一外力项来实现。其速度场和温度场分布函数分别为:
(6)(7)
式(6-7)中r 为空间矢量位置,t 为时间,e i 为格子离散速
度,Δt 为时间步长,τ、τT 为无量纲松弛时间,f i e q 、g i e q
、F i 为平衡态速度分布函数、平衡态温度分布函数和外力项，其表达式分别为:
(8)(9)
满月宝宝祝福语
(10)式中w i 为权系数,ρ密度，格子速度c=Δx/Δt ,Δx 为格子步长。
宏观密度和速度可计算如下:
(11)
自助链
松弛时间系数τ、τT (12)
采用非平衡外推方法[15]进行固体壁面速度边界的处理，其整体精度为2阶，具有较好的数值稳定性。
f i (r w )=f i e q (r w )+[f i (r n )-f i e q
(r n
)](13)式中r w 代表边界上的点,r n 位离边界点最近一个内部点。3模拟结果讨论3.1模型验证
图2为方腔内不含导热块的流线图和等温线图。从图2(a )
可以看出，瑞利数R a 为103
时，流动的特点是方腔中央出现一个近似圆形的涡，而在图2(b )等温线图中，瑞利数R a 为103时其等温线近似垂直。这与文献[6,7,15]的计算结果基本一致。
为
21
()()p t
u u u u F 2
[]p
T T T t c u ref T T F g 3
v g TL Pr Ra钮扣结
i i i i i 1
(,t)(,)[(,)(,)]eq i
f r e t t f r t f r t f r t tF i 1
远古动物(,)(,)[(,)(,)]
橙子英语eq i i i i T g r e t t t g r t g r t g r t 22
i i i 242()[13 4.5 1.5]
eq
i e u e u u f w c c c
22
i i 242()[13 4.5 1.5]
eq i i e u e u u
g Tw c c c
i i i i i 24
2()()1(1[393]2F c c c
e F e F u e u F i i i i
i i i
,/2t
f f T
g u e F,2233/()0.5
/()0.5T p
c dt c dt c
(a )Ra=103(b )Ra=105(c )Ra=106
图3不同瑞利数下流线分布情况
了进一步验证模型的准确性，还将本文获得的方腔中央涡的位
置与文献[6,7,15]中的涡位置进行比较，其涡位置也基本一致，见表1所示。说明本文构建的格子Bol t z m ann 模拟程序是可靠的。
3.2瑞利数R a 数的影响
图3-4展示了方腔内含导热块自然对流在不同R a 数下流动稳定后的流线、等温线变化趋势图，其中固体块宽度W =H /3。对于流线图，由图可以看出，随着瑞利数R a 的变化，其流线变化显著。当R a 数较小(R a=103)时，其流线表现为绕着中心导热块的多层同心圆。当瑞利数增加至105时，在导热块左右附近出现了两个涡。当R a=106时，浮升力增加，竖壁边界层的作用也大大增加，因而涡的卷起作用增强，在导热块左右附近产生的两个涡偏离方腔中心的导热块，分别向左上壁面和右下壁面移
动，且两个涡都分别出现了两个次级涡。
对于等温线图，在R a 数较小(R a=103)时，观察图4(a )等温线图，可知其等温线几乎是垂直的，表明在瑞利数偏低时冷热壁之间的传热主要由热传导引起的，即热传导占主要地位；观察图4(b )-(c )等温线图，随着瑞利数R a 不断增加，等温线在方腔中央近似呈水平分布，而热壁面和冷壁面的等温线保持垂直，且越来越密集，形成一层薄边界层，说明腔体内对流加剧，热传导逐渐减弱，对流占据主导地位。其中，由于固体块位于方腔中央，且模拟时假设固体块的热扩散系数和腔体内介质一致，因此腔体内加入导热块和无导热块时其等温线区别不大。
3.3导热块宽度W 的影响
图5-6展示了不同导热块宽度下的流线、等温线变化趋势图，其中R a 数固定为106。对于流线图的变化，对比图3(c )和
图
(a )流线图(b )等温线图
图2瑞利数R a=103时，方腔内不含导热块时流线图和等温线图
涡的位置
文献[6] 文献[7] 文献[15] 本文 x 0.2872 0.5063 0.5024 0.5011 y
0.5691
电脑截图快捷键是哪个0.5125
0.5087
0.5031
表1自然对流涡的位置(R a=103
)
(a )Ra=103(b )Ra=105(c )Ra=106
图4不同瑞利数下等温线分布情况
5(a )-(c )，可以看出，导热块宽度W 比较大(W =0.6H )时，在左右两侧高低温壁面条件作用下，形成了一个绕导体块的大涡流，且从图5(a )还可发现在导体块左右附近的两个次级涡。当W 减小至H /3时，固体块附近的两个次级涡向左上和右下移动，且分别又产生了两个三级涡，见图3(c )。随着W 的继续减小，流线变得不规则，左右两侧的两个三级涡其中一个逐渐消失。当W 减小至H /15时，在固体块左右分裂出两个小涡。对于等温线图，对比图4(c )和图6，可以看出，随着W 的减小，其等温线没有显著区别。只有在W 特别大时，其水平分布的区域比较大一些，也更平直一些，见图6(a )所示。3.4导热块在方腔内位置的影响图6-7展示了导热块在方腔不同位置下的流线、等温线
变
(a )W=0.6H (b )W=H/5(c )W=H/15
图5不同导热块宽度时流线分布情况
(a )W=0.6H (b )W=H/5(c )W=H/15
图6不同导热块宽度时等温线分
布情况
化趋势图，其中R a 数固定为106。对于流线图的变化，对比图3(c )和图7，可以发现，导热块在方腔不同位置，其流线变化显著不同。当导热块位于方腔上部时，方腔右下部靠近低温壁面处出现一个涡心，包含两个小涡。当导热块位于方腔下部时，涡心则在方腔左上部靠近高温壁面出现。另外，导热块偏上或者偏下，由于导热块的存在，流线在导热块附近都发生了扭曲。
对于等温线图，对比图4(c )和图8，可以发现，当导热块位置偏上或者偏下时，方腔中部等温线水平分布区域减小，导热块区域等温线倾斜分布，热壁面和冷壁面等温线保持垂直，形成的薄边界层仍相同，没有发生变化。
4结论
本文建立了内含导热块的方腔自然对流格子Bol t z m ann 方程模型并进行了数值模拟，重点分析了不同瑞利数R a 和导热块宽度及位置对腔体流线分布和温度分布的影响，具体结论如下:
4.1随着瑞利数的增加，在导热块左右位置处出现两个次级涡，随着瑞利数的继续增加，导热块附近的两个涡分别向左上、右下移动，且两个次级涡都又分裂出两个三级涡。
4.2不同的导热块大小对腔体内流体流型影响较大。当腔内导热块比较大时，形成了绕导体块的大涡流，以及导体块左右附近的两个次级涡，流线比较规则。随着导热块宽度减小，流线变得不规则，腔内漩涡的数量增加。
4.3改变瑞利数及导热块大小其温度场与不含导热块方腔自然对流的温度场相比，区别不大。
4.4导热块位置对腔内流型及温度场影响都比较大。当导热块位于方腔上部时，方腔右下部靠近低温壁面处出现一个涡心。当导热块位于方腔下部时，涡心则在方腔左上部靠近高温
壁面出现。当导热块位置偏上或者偏下时，等温线在
导热块区域是倾斜分布的。
参考文献
[1]Chen S Y,Doolen G D.Lattice Boltzmann method for fluid flows [J].v.fluidMech,2012,30(1):329-364.
[2]Davis G D V.Natural convection of air in a square cavity:A bench mark numerical solution [J].International Journal for Numerical Methods in Fluids,1983,3(3):249-264.
[3]徐丰，崔会敏.侧加热腔内的自然对流[J].力学进展,2014,44:98-136.
[4]Succi S.Applied Lattice Boltzmann Method for Transport Phenomena,Momentum,Heat and Mass Transfer.[J].Canadian Journal of Chemical Engineering,2010,85(6):946-947.
[5]何雅玲，李庆，王勇，等.格子Boltzmann 方法的工程热物理应用[J].科学通报,2009,54(18):2638-2656.[6]王立海，李春光.基于LBM 的自然对流换热问题的数值模拟[J].宁夏师范学院学报,2016,37(3):20-25.[7]童长青，何雅玲，王勇，等.封闭方腔自然对流的格子-Boltzmann 方法动态模拟[J].西安交通大学学报,2007,41(1):32-36.
[8]曹先齐，杲东彦，蔡宁，等.倾斜角对方腔内非稳态自然对流影响的格子Boltzmann 方法模拟[J].南京工程
学院学报(自然科学版),2019,17(2):55-60.
[9]陈磐，杲东彦，王亮，等.正弦温度分布边界条件下腔体内自然对流的格子Boltzmann 方法模拟[J].南京工程学院学报(自然科学版),2019,17(1):64-69.
[10]朱建奇，侯佳煜，杲东彦，等.基于格子Boltzmann 方法的倾斜方腔自然对流模拟[J].南京师范大学学报(工程技术版),2018,18(04):25-32.
[11]严微微.用格子Boltzmann 研究多孔介质内的自然对流换热问题[J].西安石油大学学报(自然科学版),2007,22(2):149-152+182.
[12]Mejri.LBM simulation of natural convection in an inclined triangular cavity filled with water [J].Alexandria Engineering Journal,2016,55(2):1385-1394.
[13]何宗旭，严微微，张凯，等.底部局部加热多孔介质自然对流传热的格子Boltzmann 模拟[J].物理学报,2017,66(20):156-165.[14]郭照立，郑楚光.格子Boltzmann 方法的原理及应用[M].北京：科学出版社,2009.
[15]何雅玲，王勇，李庆.格子Boltzmann 方法的理论及应用[M].北京：科学出版社,2009.
基金项目：江苏省大学生实践创新训练计划项目(项目号201911276023Y )。
作者简介：陈磐，女，博士，讲师，研究方向为传热传质学；郑竟波，女，本科在读，就读于南京工程学院能源与动力工程专业
。
图8导热块在方腔位置不同时等温线分布情况
图7导热块在方腔位置不同时流线分
我心中的桥
布情况
