fluent炉膛仿真教学教程文档

.炉膛仿真过程及其其中的问题一、（Gambit）几何建模部分1.大体尺寸在本次设计中，（实际标高-5=图中的标高）锅炉的尺寸为：锅炉高度为26890mm，精品文档放心下载宽度为7570mm，深度为7570mm。燃烧器的高度为2.105m，最底层的燃烧器低端距冷灰斗距离为2.1775m。谢谢阅读采用四角切圆（顺时针切圆，假想切圆直径0.8m）的均等配风燃烧方式。其中一次风谢谢阅读2层，二次风3层。由低到高燃烧器风口布置依次为二、一、二、一、二。燃烧器宽度为谢谢阅读.0.4m，一次风口高度0.2405m，二次风口高度0.352/0.315m，风口间距为谢谢阅读0.21/0.12/0.155m。2.简化处理将水冷壁简化成一个恒温平面；将燃烧器简化成一个平面，各次风口为平面中的一个矩形区域，作为速度入口；感谢阅读忽略屏式过热器，将折焰角上方与水平烟道相连结的平面作为出口（outflow）。感谢阅读.3.几何建模过程及网格划分为了方便锅炉的网格划分，我们将整个计算域划分为5个区域：冷灰斗下端至燃烧区感谢阅读域下端、燃烧区域、燃烧区上端至折焰角下端、折焰角区域、折焰角上端至炉膛出口。感谢阅读3.1点线面的生成几何建模的方法通常可以是自下而上的，即先生成体的各个点（通过坐标确定位置）；精品文档放心下载将生成的点依次连接成线；将线围成体的各个面；最后将面组合成一个实体。谢谢阅读当然建模时也可以通过设置实体（面）的长宽高（长宽）直接生成。精品文档放心下载3.2实体分割块的划分方法如下：先产生一个面，并将该面平移至该实体要切割的位置，splitvolume选卡中，splitwith谢谢阅读选择face（real），然后选中要切割的实体（对应splitvolume中的volume）以及用来精品文档放心下载切割这个体的面（对应face栏）（注意：在切割时需要选中Connected，保证切割产生的感谢阅读两个体之间的面是公共面，而不是两个重合的面。因为公共面可以通过物质和能量，而重精品文档放心下载合的面不加定义时是wall），最后点击APPLY确定。精品文档放心下载.根据这种方法，我们可以在Z方向将燃烧区分为很多层，方便以后设置一、二次风入感谢阅读口的边界条件。同时，在xy平面内燃烧区被分为8份，如图所示：感谢阅读3.3网格划分网格划分的最后记结果如图所示：.这种网格的特点是：四个角的地方网格比较密，而中间网格比较稀疏。同时网格线的方感谢阅读向与流动合速度方向重合度比较高。这样的网格划分可以很好的抑制伪扩散的发生。谢谢阅读这种网格的划分步骤如下：在将区域分块的基础上对实体按照线、面、体的顺序进行依次划分。精品文档放心下载Edge的划分：为了形成这种对称的网格，我们需要对edge进行划分，如图所示：谢谢阅读.其中，softlink采用maintain形式，Spacing选用Intervalcount（划分数目）。在精品文档放心下载本设计中，将每条线均分为30份，即ratio为1，intervalcount为30，其他保持不变。谢谢阅读Face的划分：由于前面对每条边进行了划分，所以对面的网格划分就只需要设置网格精品文档放心下载的形式和类型如图所示：.其中，Elements采用Quad形式，Type采用Map形式（映射成结构化网格）。此时精品文档放心下载不需要对Spacing进行设置了。Volume的划分：对volume的划分，我们采用Cooper（制桶）方式。采用这种划分谢谢阅读方式时，有一点需要注意，就是上下两个Face的网格划分要完全一样，也就是说组成Face感谢阅读的Edge的划分也要一样。如图所示：.其中，Element采用Hex/Wedge形式。Sources表示需要选择制“桶”的上下两个谢谢阅读面。Intervalcount表示两个面之间划分的数目。本设计中，根据风口和墙面的高度进行划精品文档放心下载分，每个网格高度在0.1m左右。最后依照上面的方法和步骤对燃烧区的每一层进行这样的网格划分。感谢阅读对于除了燃烧器区的其他区域的网格划分，要求就比较低一些了。对我们依旧采用感谢阅读COOPER的方式对体进行划分。不过其他地方的Sources是沿y轴方向的两个面（燃烧器感谢阅读区域的sources是沿z方向的）。.最后的网格为：冷灰斗30*30*30；燃烧器30*30*8*（3+2+3+1+3+1+3+2+3）；燃谢谢阅读烧器上端至折焰角：50*50*66；折焰角：50*50*16；折焰角上方：40*50*40。最后网格精品文档放心下载数目大概在480000个，其中燃烧器区域网格为151200个。感谢阅读3.4交接面处的处理在划分计算域的时候会涉及到interface的设置。在燃烧区的上下两个端面，我们需要感谢阅读分别将这个面与其相重合的那个面设置成一对interface。因为燃烧器区与相邻的两个实体感谢阅读并不是通过分割而来，是3个独立的实体，为了能让物质和能量通过该重合的面，需要通谢谢阅读过设置interface来实现，如图所示：.由于燃烧器区域上端的xy平面被划分为了8块，所以需要将这8个面一起设置为感谢阅读interface11，然后将与燃烧器区上端重合的面设置为interface12。对于interface21和精品文档放心下载22的设置和上述一样。4.边界条件设置在gambit中需要预先设置边界条件。将折焰角上方与水平烟道相连接的那个面设置为outflow边界条件。感谢阅读.Interface的设置上面已经说过了，下面我们进行一、二次风入口的设置。根据燃烧器谢谢阅读的结构确定各次风口在模型中的位置，然后将边界条件的Type设置成velocity_inlet（速精品文档放心下载度入口）。名字格式为ofa/pa/sa+两位数字，数字前一位表示在xy平面所处的象限，后一感谢阅读位表示自高向低同类型风口的层数。如图所示：二、Fluent仿真过程0.网格导入、Interface设置以及网格检查.在完成Gambit中的工作后，需要将生成的.msh文件导入到Fluent中。谢谢阅读0.1网格导入、检查以及解法器设置在General中点击Check完成网格检查（网格检查中不能出现网格体积为负数的情况，感谢阅读否则会出错，需要重新进行稽核建模）。点击ReportQuality进行网格质量检查。精品文档放心下载在解法器中选择Pressure-Based、Absolute、Steady的情况。勾选Gravity，建立重感谢阅读力场（z=-9.81m/s2），设置如图所示：0.2Interface设置点击MeshInterfaces中的Create。在InterfaceZone1中选择interface11，在精品文档放心下载InterfaceZone2中点选interface12，MeshInterface名称为interface1，点击Create精品文档放心下载设置完成。按照同样的方法设置interface2，如图所示：感谢阅读.1.燃料及边界条件参数确定1.1燃料特性及风煤计算1.1.1燃料计算工况符号单位3#T-017#T-018#T-01适用标准全水分Mt%GB/T211-2007空气干燥基水分Mad%2.531.822.17收到基灰分Aar%27.8733.2236.39GB/T212-2008干燥无灰基挥发分Vdaf%41.0144.0039.92收到基碳Car%51.3647.5045.49收到基氢Har%3.493.433.09DL/T568-1995收到基氮Nar%0.820.760.72.收到基氧Oar%6.676.175.86全硫St,ar%0.690.220.25GB/T214-2007收到基高位发热量Qgr,v,arMJ/kg20.4819.1818.19GB/T213-2008收到基低位发热量Qnet,v,arMJ/kg19.5518.2717.36根据表格，我们将元素分析数据转换成干燥无灰基的挥发分的元素组成。由于干燥无灰感谢阅读基无水、无灰，故剩下的成分不受水分和灰分的影响，是表示碳、氢、氧、氮、硫成分百分谢谢阅读函数最稳定额基准，所以通常选择转换为干燥无灰基来计算。感谢阅读各种煤不同分析基之间的换算公式为X=0×其中，X0，X分别为某成分原基准与新基准的质量百分数，%；K为换算系数。精品文档放心下载收到基与干燥无灰基之间的换算系数为K=

100 ar−ar利用上述公式将煤的收到基转换为干燥无灰基的元素组成，如下表所示，这些数据也是感谢阅读Fluent中计算PDF元素分数所需要的。PDF中元素分析元素分数元素 C H O N元素组成 0.820971 0.055766 0.105757 0.017506精品文档放心下载PDF中工业分析分析组分组分 V FC A M组分分数 0.2212 0.3329 0.3639 0.082精品文档放心下载.1.1.2风煤计算锅炉实际燃煤量t/h26.015设计值一次风流速m/s24.3设计值一次风份额%28.28设计值一次风温K303温风份额%64.09设计值二次风温K600一次风口面积m24*0.1924二次风口面积m24*0.4075注：由于不知道乏气送粉的位置，将乏气份额归并到二次风中，即二次风份额为71.72%。精品文档放心下载根据克拉伯龙方程PV=nRT和表格数据可知：标况下一次风速 =. ⁄；再根据一、二次风的份额和面积可以得到：二次风速 = ∗二、一次风份额之比∗一、二次风口面积之比=. ⁄；感谢阅读根据克拉伯龙方程可知：实际二次风速 =. ⁄；根据煤量，可以知道每个一次风口煤的质量流量：coal= ÷ ÷ ÷ =. ⁄根据切圆直径和炉膛尺寸可以知道风煤的入口方向：夹角θ=40.55°；cosθ=0.76；sinθ=0.65.谢谢阅读1.2边界条件设置.现以一次风pa11为例，介绍对流场数据的设置。在Fluent中BoundaryConditions菜单下找到pa11项目，如图所示。然后这个风口感谢阅读进行设置。点击Edit进入设置页面，如图所示。.VelocitySpecificationMethod选项中选择MagnitudeandDirection（速度大小和方向）；在VelocityMagnitude中填入24.3m/s；精品文档放心下载CoordinateSystem中选择Cartesian（X,Y,Z）笛卡尔直角坐标系，然后在下面依次填入流体流动的方向（X轴为-sin40.55°，Y轴为-cos40.55°）。在第一象限的风口方向为（-sin40.55°，-cos40.55°）；在第二象限的为（cos40.55°，-sin40.55°）；第三象限的为（sin40.55°，cos40.55°）；第四象限的为（-cos40.55°，sin40.55°）。谢谢阅读在设置湍流参数时，我们选用IntensityandHydraulicDiameter（湍流强度和水力直径）方式。湍流强度I我们设置成10%，为强湍流状态，水力直径D的设置根据公式D=4⁄谢谢阅读设置成0.3m（二次风口的水力直径为0.374m/0.352m）。精品文档放心下载Species选卡中将MeanMixtureFraction设置为0（氧化剂入口），MixtureFractionVariance设置为0，如图所示。感谢阅读.在DPM选卡中将DPBCType设置成reflect（反射）精品文档放心下载其他设置保持不变。依照上面的方法，可以完成对PA和SA流场参数设置。谢谢阅读2.与流动和燃烧相关的模型设置在Fluent中我们打开Models选项。在中意菜单中，我们可以设置包括流动、传热、精品文档放心下载燃烧等方面的模型。2.1.气相流动模型本文采用标准k-e双方程湍流流动模型，同时采用标准壁面函数处理近壁面的流动问题，精品文档放心下载其中的流动参数保持默认参数。其设置如下：.2.2气相湍流燃烧模型模拟气相湍流燃烧过程的关键在于如何模化湍流燃烧反应率。针对扩撒火焰的模型有精品文档放心下载k-ε-g模型，混合数-概率密度函数模型。为了减少计算量，采用但混合数PDF模型。对于煤粉燃烧，我们在Models->Species谢谢阅读中选用Non-PremixedCombustion（非预混燃烧）模型，然后再弹出的菜单中进行相关感谢阅读.的参数设置，如图所示。在PDFTableCreation栏目中点击Chemistry选卡。在StateRelation中选择谢谢阅读Equilibrium（化学平衡法），EnergyTreatment中选择Non-Adiabatic（非绝热）形式，精品文档放心下载StreamOptions选择EmpiricalFuelStream（经验燃料流）。谢谢阅读点击CoalCalculator会弹出对燃料特性进行设置的对话框，根据燃料特性表的中数据，谢谢阅读可以设置完燃料工业分析和元素分析的参数。其中物料的名称为coal-particle，HCV为感谢阅读1.819e+07j/kg，其他的保持默认数值，点击Apply和OK确认，我们可以看到Model感谢阅读Settings中的数据发生了相应的变化。如图所示。谢谢阅读.点击Boundary选卡将燃料温度设置为303K，氧化剂的温度设置为600K。精品文档放心下载点击Table选卡中的CalculatePDFTable进行燃烧的计算。计算完成后，我们可以点感谢阅读击DisplayPDFTable查看关于煤粉燃烧的数据，如图所示。谢谢阅读.计算结果显示燃烧形成的成分有20种，点击Control选卡可以查看成分名称，如图所精品文档放心下载示。设置完成后点击Apply和OK。同时我们可以看出Models中的Species变成了感谢阅读Non-PremixedCombustion模型。谢谢阅读.2.3煤粉燃烧模型煤粉燃烧可分为煤粉预热、挥发分析出和燃烧过程、焦炭燃烧等过程。谢谢阅读在本设计中，挥发分析出模型采用单速率析出模型，焦炭燃烧模型选用扩散-动力控制感谢阅读模型。在Fluent中首先需要对煤粉颗粒的喷射进行相关的设置。感谢阅读在Models菜单中点击DiscretePhase，弹出对话框，如图所示。精品文档放心下载在Interaction中勾选InteractionwithContinuousPhase（对连续相的影响，但是谢谢阅读在仿真时候需要先建立无颗粒相的流场，即在仿真开始时不勾该选项）。谢谢阅读在Trackking选项中，Max.NumberofSteps设置为3000，勾选SpecifyLength精品文档放心下载Scale，其中，LengthScale设置为0.01m；感谢阅读.在PhysicalModels选项和Numerics选项中的参数和选项均为默认设置。谢谢阅读为了让煤粉能够喷射入炉膛，我们需要进行对颗粒相的设置，点击 DiscretePhase谢谢阅读Model中的Injections弹出对话框，如图所示。谢谢阅读点击Create对煤粉颗粒进行参数设置。以某一次风口进入的煤粉为例，如图所示。感谢阅读.在InjectionType中选择surface，然后在ReleaseFromSurface中选择相应的煤粉精品文档放心下载喷射地点。在ParticleType中勾选Combusting（燃烧）。感谢阅读在设置DiameterDistribution时，可以选择uniform（均匀）形式，当然最好选择感谢阅读rosin-rammler形式。在PointProperties中设置好速度方向（与空气速度矢量一致）、颗粒直径（1e-06m）、感谢阅读颗粒温度（303K）、质量流量（0.903kg/s）等参数；精品文档放心下载在TurbulentDispersion中的StochasticTracking选择DiscreteRandomWalk精品文档放心下载Model（随机轨道模型），如图所示。.依照上述方法对其他的煤粉喷射源进行参数和模型的设置选择。谢谢阅读煤粉喷射设置完后可以发现Models->Materials选卡中多了CombustingParticle谢谢阅读（coal-particle）一项。接下来我们对这个成分进行相关设置，如图所示。感谢阅读.在弹出的对话框中的Properties中设置DevolatilizationModel（挥发分析出模型）精品文档放心下载为single-rate（单速率）；同时设置CombustionModel（焦炭燃烧模型）为谢谢阅读kinetic/diffusion-limited（扩散-动力控制）；谢谢阅读对于该颗粒相的其他物性参数保持默认参数。2.3辐射传热模型炉内的能量主要通过辐射的形式进行传递。Fluent软件提供了多种辐射模型，在本设感谢阅读计中，我们选用P1模型，如图所示。.在Model中点选P1即可。P1模型考虑了辐射散射作用，更适用于光学厚度较厚以及几何结构复杂的燃烧设备。精品文档放心下载但是该模型也有缺点，包括对来自内部热源的辐射热通量有过高估计的趋势等。精品文档放心下载在选用P1模型后，Models中Energy选卡自动开启。感谢阅读2.4边界条件的设置在BoundaryConditions中设置Supersonic/InitialGaugePressure为0pascal；一感谢阅读次风入口温度为303K；二次风入口温度为600K；DPBCType为reflect。精品文档放心下载墙壁为固定无滑移壁面；热力条件为很稳壁面600K，内部发射率为1，壁厚及生热率感谢阅读均为0；壁面材料为Al（铝）；DPM边界条件数类型BCType为reflect。其他保持默认设精品文档放心下载置。出口边界条件设置为outflow，保持默认设置。3求解过程（Solution）.3.1SolutionMethods点选Solution->SolutionMethods。谢谢阅读压力-速度耦合采用SIMPLE格式；空间离散方法中，梯度采用格林-高斯单元法（Green-GaussCellBased）；感谢阅读压力采用Standard方法；其他均采用二阶迎风格式（SecondOrderUpwind）。精品文档放心下载3.2SolutionInitialization谢谢阅读初始化方法采用StandardInitialization方法；感谢阅读初始温度为600K；其他保持默认参数和设置，点击Initialize完成初始化，如图所示：感谢阅读.3.3求解在求解时，需要先建立连续相的流场，然后再加入颗粒相耦合修正流场，即在开始时精品文档放心下载Models->DiscretePhase中不勾选InteractionwithContinuousPhase，迭代100步先谢谢阅读建立连续相流场；迭代完成后，勾选InteractionwithContinuousPhase，其他保持默认设置，然后进谢谢阅读行迭代，如图所示：.三、仿真结果中出现的问题在迭代步数为2670步基础上的仿真结果1.炉膛内切圆直径与假想切圆直径相差不大，说明有可能迭代步数不够，如图所示：精品文档放心下载.2.在燃烧器上端至出口的区域内，温度变化不大，所散点图所示：精品文档放心下载3.现场测试数据显示乙侧温度都高于甲侧温度，根据经验表明是顺时针切圆；但在仿真精品文档放心下载中采用顺时针切圆，结果却是乙侧温度低于甲侧温度，标高24.2m处温度情况如图所示：感谢阅读标高27.5m处温度情况如图所示：.4.各标高所在平面温度等值图4.1下二次风4.2下一次风4.3中二次风.4.4上一次风4.5上二次风4.6标高24.2处温度.X=-1.685处Y坐标（m）温度（K）-3.785 600-3.63361449.03-3.48221511.38-3.33081522.37-3.17941530.99-3.028 1538.02-2.87661543.94.-2.72521549-2.57381553.35-2.42241557.08-2.271 1560.24-2.11961562.91-1.96821565.17-1.81681567.08-1.66541568.69-1.514 1570.07-1.36261571.25-1.21121572.24-1.05981573.07-0.90841573.77-0.757 1574.35-0.60561574.83-0.45421575.22-0.30281575.52-0.15141575.75-7.35393e-18 1575.910.1514 15760.3028 1576.020.4542 150561575.780.757 1575.470.90841575.041.05981574.441.21121573.631.36261572.611.514 1571.331.66541569.761.81681567.851.96821565.542.11961562.782.271 1559.532.42241555.752.57381551.372.72521546.332.87661540.573.028 1533.933.17941526.223.33081517.183.48221506.053.63361449.113.785 600.4.6标高27.5处温度X=-1.467m处Y坐