封闭方腔自然对流换热

封闭方腔自然对流换热描述该物理模型的无量纲方程组为：连续性方程：()()0d U d U dx dyρρ+= 动量方程：2222U V P U U U V X Y X X Y ∂∂∂∂∂+=-++∂∂∂∂∂ 2222PrU V P U U RaU V X Y X X Y ∂∂∂∂∂+=-+++Θ∂∂∂∂∂ 能量方程：22221Pr U V X Y X Y ⎛⎫∂Θ∂Θ∂Θ∂Θ+=+ ⎪∂∂∂∂⎝⎭其中，无量纲几何参数,x y XY l l == ；无量纲速度ul U v =，vlV ν= ；无量纲压力()02/p gy p v l ρρ+=，无量纲温度0h c T T T T -Θ=-；普朗特数Pr pc v a lμ==；瑞利数()3h c g T T l Ra va β-=，空气的体胀系数1pT ρβρ∂⎛⎫=- ⎪∂⎝⎭，λ 为空气的导热系数。

、具体模拟计算参数：55353,500,360,0.0033331.74510,Pr 0.712, 2.36101.11/, 1.9310h c L m T K T K v a kg m βρμ---=====⨯==⨯==⨯对方腔划分网格，采取的是60⨯60网格，，壁面处加密。

在FLUENT 软件中，使用分离求解器求解控制方程组。

材料的物性设置密度使用Boussinesq 假设。

本例主要分别计算了数为34456110,110,510,110,110⨯⨯⨯⨯⨯的情况。

压力插值方案选择Body Force Weighted 格式；压力-速度耦合方程用SIMPLE 算法；动量、能量方程选择二阶迎风格式。

有公式：()3h c g T T l Ra vaβ-=可得对应的g 入下表所示本模拟与文献中的Nu 比较模拟图对比：Ra=3110⨯ T VRa=4110⨯ T VRa=4510⨯ T VRa=5110⨯ T V。

封闭腔内自然对流的数值研究金涛【摘要】在本文的研究当中,笔者利用SMPLE算法,采用二阶迎风格式的对流扩散项,,建立了封闭腔内的自然对流物理模型,并通过该物理模型计算与研究了在封闭腔内部的自然对流换热.本文的研究最终得出,在一定的Ra下,长度不同的阻流件的平均Nu数,有水平阻流件的封闭方腔与无阻流件时相比,相同条件下自然对流的换热系数随阻流件长度的增加先略减少,然后增加.同时在阻流件的长度发生变化的前提之下,封闭腔内部的环流也会随着增加,一般会增加2个或者3个,另一方面,不同壁面上的阻流件的布置方式也会对换热产生不同的影响.【期刊名称】《赤峰学院学报（自然科学版）》【年(卷),期】2017(033)005【总页数】2页(P6-7)【关键词】不可压Navier-Stokes;Boussinesq方程组;自然对流;数值算法【作者】金涛【作者单位】中国矿业大学银川学院,宁夏银川 750021【正文语种】中文【中图分类】O35在实际的工程应用当中，如建筑的暖通空调与制冷装置的设计当中，当封闭腔内存在一定量的空气时，内壁与外壁之间会存在着一定的温差，从而导致空气会产生自然对流的现象.许多专家与学者对此问题进行了研究，研究的主要方式都是通过建立必要的数学分析方程组，并对其进行求解，从而得出封闭腔内部气体的运动规律，这些研究也被广泛的运用在最近几年的工程实际当中[1].为了使得封闭腔内的换热得以有效的减弱，往往会在封闭腔内部设置竖直隔板，从而使得封闭腔可以被完全或者部分隔断.在实际的工程运用当中，人们发现对存在水平等温阻流件的封闭腔进行研究具有非常重要的实际应用价值，基于此种情况本文的研究是针对存在水平等温阻流件的情况下进行的，通过研究之后得到了水平阻流件对封闭腔内部自然对流传热的影响.2.1 物理模型的建立本文当中所研究的物理模型如下图1当中所示，在该模型当中封闭腔的长度为L，腔壁为上下绝热，两侧等温，其中热壁的温度为Th，冷壁的温度为Tc，同时在腔壁上还存在着水平阻流件，阻流件的长度为h，厚度为Ø，因此就可以将该问题的物理模型简化在二维当中，在封闭腔内部的气体的密度变化与不可压气体的理想气体方程之间存在着高度的一致性，所有的参数均取做常数，Ra=gβΔTL3/(αv),在该方程当中，α表示的是热扩散系数，g表示的则是重力加速度，β为流体容积膨胀系数，L为特征尺度,v为运动粘度,温度的差值则表示为ΔT=Th-Tc,在本文的研究当中Ra=1.0× 105.2.2 控制方程的建立在封闭腔内部发生自然对流换热的速度以及温度的边界条件为：绝热壁面的速度与温度的梯度都需要保持为0，同时隔热板与左壁面的温度保持在Th，右壁面的温度则保持在Tc，水平阻流件的长度则分别为L1与L2.则可以采用以下数学方程对其进行描述[2].2.3 计算结果采用有限容积法建立相关方程，方程的求解采用SMP L E方法，采用二阶迎风格式的对流扩散项，将100×100的均不网格设定为计算网格.2.3.1 阻流件长度变化而导致的流场和温度场的变化在阻流件的长度发生变化的过程当中，封闭腔内部的流场也会产生较为明显的变化，尤其是当封闭腔没设置的阻流件的长度较长时，开始的环状流场会被不断的压缩，最终形成两个流场，在这个过程当中等温线则会受到压缩，在下图当中，给出了当隔板的相对长度为0.1、0.4与0.6时的温度场与流场的变化.2.3.2 N u数的变化在上图五当中，笔者对存在水平隔板的情况下，封闭腔内部自然对流随着阻流件长度的变化N u数所发生的变化，从上图当中我们可以发现，在阻流件长度发生变化的前提下，封闭腔内部的自然对流的强度也会发生一定的变化，随着长度的增加，强度首先是从高到低变化，随后则不断增高，同时上升的趋势也逐渐平缓，在阻流件的长度超过0.5时，上升的幅度不断增加.导致这一现象的主要原因在于长度较低的阻流管会对环形流场造成一定的破坏，但是随着阻流件长度的增加，封闭腔内部的环流被压缩成两个，则使得封闭腔内部的气体流动速度得以增强，随着长度的进一步加大，环流被进一步压缩，从而使得对流明显增强[3].（1）随着阻流件长度的增加，封闭腔内部的环流数量会不断增加，先是增加一个，最后增加两个.（2）随着阻流件长度的增加，在封闭腔内部所发生的对流换热的强度在刚开始时会发生一定的下降，随后开始稳步上升.（3）不同的阻流件布置方式会对封闭腔内部的流场以及对流换热的强度都会发生一定的影响[4].〔1〕Oztop H F,Abu-Nada E.Numerical study of natural convection in partially heated rectangular enclosures filled w ithnanofluids[J].International Journal of Heat& Fluid Flow,2008,29(5):1326-1336.〔2〕Davis G D V.Natural convection of air in a square cavity,a benchmark numerical solution.Int J Numer Methods Fluids 3:249-264[J].International Journal for Numerical Methods in Fluids,1983,3(3):249-264.〔〕International Journal for Numerical Methods in Fluids, 1994,18(7):695–719.〔3〕王宇飞,徐旭,王文龙,范利武,俞自涛.封闭腔内Al_2O_3-EG纳米流体自然对流传热特性的数值研究[J].能源工程,2014（01）:1-6.〔4〕阳祥,陶文铨.高瑞利数下封闭腔内自然对流的数值模拟[J].西安交通大学学报,2014（05）:27-31.。

封闭腔内水自然对流换热数值模拟
自然对流换热是一种重要的热传递方式，它在许多工程和科学
领域都有着广泛的应用。

在封闭腔内，水的自然对流换热特性对于
工业设备的设计和运行具有重要意义。

为了更好地理解和优化这一
过程，数值模拟成为了一种重要的研究手段。

通过数值模拟，我们可以利用计算机模拟封闭腔内水的自然对
流换热过程，从而研究其传热特性。

在模拟过程中，我们需要考虑
腔体的几何形状、水的流动状态、温度分布等因素，以及流体的物
性参数。

通过数值方法，我们可以计算出不同条件下水的温度分布、传热速率等关键参数，从而为工程实践提供重要的参考。

在实际工程中，封闭腔内水自然对流换热数值模拟的研究成果
可以为工程设计和优化提供重要依据。

通过模拟分析，我们可以评
估不同工况下的换热性能，指导设备的优化设计和运行参数的选择。

同时，数值模拟还可以帮助我们理解自然对流换热的机理，为工程
实践提供科学依据。

总之，封闭腔内水自然对流换热数值模拟是一种重要的研究手段，它为工程设计和优化提供了有力的支持。

通过模拟分析，我们
可以更好地理解和控制自然对流换热过程，为工程实践提供科学依据。

随着计算机技术的不断发展，数值模拟将在工程领域发挥越来越重要的作用。

图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 separated 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-),女,籍贯:湖北襄阳,硕士,助教,研究方向:电力生产技术、节能减排技术、电气自动化。