水煤浆在双锥逆喷燃烧器内燃烧过程的数值分析

0 引 言

水煤浆含有30%～40%水分,其燃烧过程比煤粉燃烧需要更长时间来完成水分蒸发。在这段时间内需要保持较高的温度和氧化剂浓度[1]。水煤浆流化床锅炉通过在炉膛内流化大量石英砂为水煤浆燃烧提供温度场[2]。其点炉过程需要消耗大量的动力和燃油,耗时、耗资巨大,且锅炉负荷可调范围较小。直接配置在锅炉前端的水煤浆旋流燃烧器,在空间较大的炉膛内无法有效组织燃烧,普遍存在点火困难和燃尽率低的问题。而设置了卫燃带的水煤浆预燃室虽然解决了点火问题,但是燃烧组织已不在控制范围内,导致预燃室内部结渣严重[3]。

煤炭科学研究总院自主研发的煤粉双锥燃烧器是一种先进的煤炭燃烧利用技术体系[4]。经工程实践验证其具有热效率高、与天然气燃烧可比拟的低污染物排放、宽负荷调节、启停时间短等优势。这些特点应用于水煤浆燃烧,符合其技术发展方向,对污泥、污水的减量化、资源化处理具有积极意义[5]。

从燃烧组织分析,煤粉双锥燃烧所采用的强化着火燃烧组织形式大幅降低点火功率[6],同时燃烧器内强旋流能够提高燃料在高温区的停留时间。基于以上设计原理有效地解决了水煤浆燃烧停留时间和温度之间的矛盾,但是目前相关研究较欠缺。本文以14 MW双锥燃烧器为研究对象进行数值模拟,采用组分输运模型对燃用煤粉和水煤浆分别进行三维数值模拟,通过比较2种燃料的流场特点、温度、组分分布特点,讨论了水煤浆双锥燃烧组织的特性。

（3）预习复习模块：提供每周的课程教学方案。该模块的课程教学方案以“星期”为周期，提供给学生上周、本周、下周的需要预习和复习的章节内容、目的、重难点。学生可以自行安排适合的时间查看，经调查问卷统计，该模块学生访问量最多，访问时间在周末和课程即将开始前。

二甲双胍是治疗2型糖尿病的基石，它的优势不仅体现在降糖方面，还可以降低胆固醇和甘油三酯，升高高密度脂蛋白胆固醇，抗血小板聚集，改善纤溶酶活性，对心脑血管具有保护作用。二甲双胍这已被使用了60年的“年过半百”的老药，其安全性值得肯定，患者完全可以按照医嘱放心地长期使用！

1 网格划分及模型假设

双锥逆喷燃烧器的设计运行程序为：通过点火燃油将燃烧器内的温度提升至300～500℃。一次风携带燃料经一次风管到达回流帽,并被强制回流进入燃烧室前锥;二次风通过旋流叶片以旋流形式进入燃烧室前锥,并在前锥中部与逆流的一次风和燃料相遇。在很短的时间内,燃料被加热到点火温度。着火后到的燃料随二次风边旋转边燃烧,并沿着二次风的流向继续向前运动,经过后锥的加速后喷出燃烧器进入炉膛。

底下的网友们一边倒地喷那个司机。大家纷纷留言，“这场雪太大了，环卫工人凌晨就起来工作，多不容易，竟然还好意思破坏环境”、“那是中央大街，如果有垃圾没及时处理环卫工人是要被罚款的”、“真想人肉那个司机，什么素质啊”

14 MW双锥燃烧器的几何拓扑结构为5层圆柱体和圆锥的嵌套(A～E),如图1所示。由外至内通过创建5层O-block(O-1～5)描述拓扑结构,可解决圆弧形状Block边界处网格扭曲,并且在近壁面处生成较理想的边界网格。燃烧器网格划分及边界条件如图2所示。考虑到连续相入口为带压射流,离散相为颗粒射流,故将一次风和二次风的入口设置为质量入口边界条件和反射离散相壁面;出口为压力出口边界条件和逃逸离散相壁面。该模型计算域采用结构化网格,网格节点数为55万,2×2×2行列式检查最小值为0.693,网格角度最小值为45°。

  

图1 O-block嵌套与燃烧器的拓扑结构Fig.1 O-block and the topology of the burner

  

图2 网格划分及边界条件Fig.2 Meshing and boundary conditions

为建立稳定可靠的计算模型,采用如下假设：①假定水煤浆以充分雾化的形式进入计算域,水煤浆液滴的运动过程不发生结团、黏壁等现象。②使用两步总包反应原理描述燃烧反应过程;假定所述燃烧反应是快速化学反应,即混合时间尺度远大于反应时间尺度;燃烧反应所处流域为完全湍流状态,即高雷诺数。③模拟是稳态计算过程,燃料量、空气量及燃烧器所有其他运行参数均不随时间变化。不包含点火、熄火等动力学细节描述,不涉及混合与反应时间尺度相近的其他现象。

2 数学模型及计算步骤

2.1 基本模型

双锥逆喷燃烧器中煤粉和水煤浆在z=0平面的速度分布如图4所示。可知,前锥扩散角降低了二次风的速度,并提供了足够的空间以实现动量的充分传递。一、二次风速度趋于一致后,经后锥收缩角的加速,形成高速射流从预燃室喷出,双锥燃烧器这一特性对保持火焰刚性具有重要意义。一次风和二次风在前锥中部形成低速区,为煤粉或水煤浆的加热和着火提供了时间,这对稳燃和点火十分关键。由于水煤浆中水分加速对系统动量的消耗,水煤浆的出口速度水煤浆比煤粉的低5%。但是2者的速度场形态类似,且在z=0平面的面积分平均速度均为47 m/s。说明,在相同的设计速度下若可以实现水煤浆的充分雾化,则煤粉的流场适用于水煤浆。

根据煤粉燃烧的模拟与工业运行参数的对比研究,本文采用Realizable k-ζ模型对湍流模型进行封闭[10]。该模型对圆形射流的扩散现象和强压力梯度边界层的模拟效果较好。对于Boltzmann数＜10的煤颗粒燃烧,需要考虑辐射传热。Fluent定义光学厚度为特定长度路径气体对特定波长辐射的减弱能力,其值越大表示越强的辐射吸收能力。燃烧器内富含水蒸气和CO2,均对辐射热具有强烈的吸收作用[11],可以视作光学厚介质。据此,本文选定P1模型计算燃烧器内的辐射传热[12]。该模型被广泛应用于颗粒燃烧的计算研究。

煤粉和水煤浆均经历水分蒸发、挥发分析出、挥发分燃烧、焦炭异相反应等。采用基于质量分数的输运模型进行描述(式(1));连续相反应使用涡耗散模型[13]计算反应速率(式(2)、(3)),反应速率受大涡混合时间尺度控制;离散相颗粒反应包括挥发分析出反应和焦炭的异相表面反应(式(4)),反应速率动力学方程见式(5),扩散作用控制方程见式(6)[14]。热解反应方程设置为两步总包反应：第1步反应是挥发分析出并氧化为CO、CO2和H2O;第2步反应为CO氧化为CO2。总包反应参数为0.5,即第1步反应CO和CO2的生成比例为1。

 

式中,Yi为i组分的质量分数,%;ρ为气相混合物的密度,kg/m3;v为速度矢量,m/s;Ji为i组分的扩散通量,kg/(m2·s);Ri为i组分的化学反应生成速率,mol/(L·s);Si为i组分自定义源相生成速率和离散相生成速率,mol/(L·s)。

 

其中,Yp为任意产物质量分数,%;Mω,f为反应物的摩尔质量,kg/kmol;下标f表示反应物;下标j取值1～N,表示任意组分;minf为所有对反应物计算中的最小值;Yf为反应物质量分数,%;A和B为经验常数,分别为4.0和0.5;ε和κ为封闭湍流方程的参数;v′i,r为反应 r 中反应物 i的化学计量数;v″r为反应r中产物的化学计量数;Mω,i为i组分的的摩尔质量,kg/kmol。Ri,r值为上述2个方程计算结果中的较小值,mol/(L·s)。

(7)安全要点。① 爆破作业做好巷道管路、电缆防护，做好设备安全，并做好人员警戒，严格执行“一炮三检”制度，做好防尘、降尘工作。② 工作台搭设要牢固可靠，爆破施工情况下工作台安全可靠，保证施工安全。

于是，在传统的数据库教学中，往往形成了教师讲授例题，学生按照教材上的例题代码原封不动地上机实践的模式，以此验证例题结果的正确性。当然，用心理解的学生也可以从例题中学到知识，再次巩固，并非不可以，可是这样便削弱了实践的真正价值。

计算过程分为冷态计算和热态计算。冷态计算中,关闭了能量方程和辐射方程,离散相和连续相为单向耦合。以一阶差分格式开始计算,待计算稳定后改为二阶差分格式。设置燃烧器出口速度为监控参数;以入口、出口质量差低于0.1%、速度残差低于1%、连续性方程残差低于10-3、其他控制方程残差低于10-3作为判断冷态计算收敛的指标。获得冷态收敛解后,激活能量方程和辐射方程,并激活离散相和连续相的双向耦合,计算直到获得收敛解。热态计算的收敛判定依据在冷态判定依据基础上增加出口温度变化监视,判定值为0.1%。

 

式中,kc为质量传递速率,mol/(m3·s);Yi,s和Yi,∞为i组分自定义源相的生成速率和离散相的生成速率,mol/(L·s)。

对于水煤浆的模拟计算,本文将包含水分的煤粉颗粒视为水煤浆颗粒。为了便于比较煤粉和水煤浆燃烧组织特性的区别,煤粉的计算工况设置为10 MW(表1)。水煤浆的计算工况采用同样的煤种和煤消耗量,煤质数据见表2。

中亚和东亚2017年的核电装机容量为101.8 GWe。该地区核电容量在2030年、2040年和2050年的高值情景预测值分别为173 GWe、258 GWe和310 GWe，低值情景预测值分别为116 GWe、131 GWe和145 GWe。

 

表1 计算工况Table 1 Operating condition of calculation

  

模拟介质 一次风风量/(Nm3·h-1)二次风风量/(Nm3·h-1)进料量/(t·h-1)煤粉 580 14 000 1.7水煤浆 580 14 000 2.4

 

表2 煤质分析Table 2 Proximate and ultimate analysis of coal

  

工业分析/%元素分析/%Qnet,ar/MadVdAdFCd Cad Had Oad Nad (MJ·kg-1)2.90 28.689.6061.72 74.794.03 7.93 0.69 24.63

水煤浆浓度设定值为60%。水煤浆水分的25%为煤工业分析水分(受颗粒表面热解规律支配),其他水分定义为射流液态水,受液滴汽化规律支配方程见式(7)、(8)[15]。

目前，临床上尚无特效药物用于治疗面部激素依赖性皮炎。常规西药或外用综合治疗该病可达到百分之八十以上，但治疗结束后复发风险较高，临床治愈率也较低。激素依赖性皮炎是随着现代医学的发展而发生的皮肤病。糖皮质激素具有抗炎，免疫抑制和抗过敏等作用，对皮肤病有一定治疗作用，但若选择不当和药物滥用，可引起激素依赖性皮炎的发生，出现脸部灼热和瘙痒，甚至影响美观，对患者社交、生活和工作均产生不良影响，需要给予有效治疗[3]。

 

 

式中,D0为扩散系数;C1为密度,kg/m3;Tp为水煤浆颗粒的温度,K;T∞为混合物温度,K;dp为水煤浆颗粒的直径,m;k为动力学反应速率,mol/(L·s);CZ为指前因子;E为活化能,J/mol;Ap为水煤浆颗粒的表面积,m2;mp为水煤浆的质量,kg;pox为氧化剂在气氛中的分压,Pa;R为反应速率。

2.2 材料特性及边界条件

如前假设,燃烧器点火不在讨论的范围,所以计算过程将煤的热解温度设定为空气温度,使煤粉颗粒进入计算域后,随即发生热解。混合物导热系数、黏性系数按照气体质量组成计算。对于气相中含有颗粒的模型,由于介质吸收和散射导致的辐射强度在经过每单位长度介质后的改变不容忽视,所以计算过程定义混合物的吸收系数模型为灰气体加权模型[16],散射系数设定为0.5。

对于燃烧器壁面,假定外界平均温度为300 K,炉内平均温度为1 400 K,结合已有工业炉运行数据,考虑壁面沾污情况,计算得到双锥的壁面热通量为60 W/m2。壁面双锥壁面均设置为离散相反弹边界,内部发射率为0.8。

2.3 计算过程

如前所述,采用离散相模型来考虑颗粒相和连续相的相互作用,并通过随机轨道模型追踪颗粒运动。本文计算过程设置每次离散型注入计算域后连续相计算40步,并追踪48 000条颗粒路径粒径以描述离散相的运动,计算过程如图3所示。

  

图3 模拟计算过程Fig.3 Simulation Calculation process

4)计算环境影响负荷。经加权后的各种环境影响潜值具有可比性，反映了相对重要性，因此可以将其综合为一个简单的指标，即环境影响负荷(EIL),它反映了所研究产品系统在其整个生命周期中对环境系统的压力大小。用公式表示为

将模型应用于7 MW试验平台,结果见表3。燃烧器出口温度误差为3%,说明前述假定和简化可确保数学模型的合理性和计算的可靠性。

由于水煤浆中水分汽化为水蒸气,使得燃烧器内的轴向阻力变大,因此在同样的二次风入口旋流条件下,水煤浆连续相获得了更大的旋流强度,这可以从离散相迹线积分出的颗粒停留时间看出。水煤浆的颗粒平均停留时间为0.33 s,而煤粉的颗粒停留时间为0.28 s。平均停留时间提高了18%,这可有效提高水煤浆在双锥燃烧器内燃烧的稳定性。

 

表3 工业运行数据与模型预测值比较Table 3 Measurement data and prediction data

  

项目 空气过量系数 负荷/(kg·h-1) 出口温度/℃模拟值 1.00 1 000 1 260测试值 0.73 960 1 300

3 计算结果及分析

3.1 速度场及迹线分析

为便于区分,煤粉和水煤浆指代所研究的2种燃料,煤粉颗粒指代煤粉的单颗粒,水煤浆颗粒指代含水的煤粉的单颗粒。在仅描述颗粒行为时,对煤粉颗粒和水煤浆颗粒不作区分,统一用颗粒表述。

  

图4 速度云图Fig.4 Contours of velocity magnitude

燃料和空气分别视作离散相和连续相。整个计算域简化为由颗粒离散相和气体连续相组成的弥散多相流体系[7]。其中,连续相的数学描述采用欧拉方法,并求解体积平均N-S方程,包括连续性方程、动量方程和能量方程[8]。能量方程中包括热流、外力、扩散黏性耗散等形式的能量平衡。分散相采用拉格朗日方法描述[9]。通过在连续相内对大量质点的追踪和路径积分运算得到颗粒的轨迹。分散相和离散相之间进行动量、质量和能量的交换,即双向耦合求解。这一过程通过交替求解连续相控制方程和离散相控制方程实现。直到连续相和离散相各参数不在随迭代的进行发生变化后,输出结果。

通过双锥逆喷燃烧器中煤粉和水煤浆的连续相迹线(图5)可以看出,二次风的旋流产生了旋转速度场使得一次风也形成了旋流,旋流在燃烧器内基本没有衰减,并且在回流帽后方产生了强湍流区。结合离散相迹线分析可知,在双旋流条件下,水煤浆或煤粉中较大颗粒获得的切向速度更大,离心力也更大,从而与O2接触的几率变大。大颗粒进入到二次风旋流通道,在燃烧器内的行程增多,其燃尽效果将得到提升。较小的颗粒在回流帽后方的湍流区富集,及时为火炬核心提供燃料,维持其高温。

会上，研发三部还与当地经销商签订“中国农资”系列新产品的销售代理协议。参会的一位经销商表示：“这几年我们一直与中农合作，当地农民对中农产品非常认可。我们与中农控股这样的大企业合作感到很放心。跟着中农走，小康路上加速走！”

  

图5 连续相和离散相的迹线Fig.5 Path line of continuous phase and discrete phase

（4）施工方案动态反馈。对施工过程中出现的问题，由队长向监控人提出申请，监控人根据实际情况制定新的技术方案并由专家组汇商，技术领导同意后组织实行。

3.2 温度场分析

双锥逆喷燃烧器中煤粉和水煤浆在z=0平面计算获得的温度场对比如图6所示。在二次风旋流卷吸和后锥强制回流的共同作用下,部分高温烟气回流至前锥,在前锥形成高温段,实现了燃料的加热和点火。由于燃料和一次风的温度较低,在回流帽出口处形成低温区,低温区左侧是挥发分和煤粉颗粒着火形成的高温区。该区域内煤粉颗粒温度经过回流烟气加热达到着火温度,并与旋流二次风充分混合接触形成高温、高燃料浓度的稳定点火区。回流帽右侧是火炬高温区,煤粉或水煤浆和空气在此范围内速度趋于一致,形成具有一定刚性的火炬。

[3]王恒,王世均,赵立合.一种高浓度水煤浆燃烧器的开发与应用[J].钢铁,1994(1)：62-66.WANG Heng,WANG Shijun,ZHAO Lihe.A highly concentrated coal-water slurry burner[J].Iron and Steel,1994(1)：62-66.

水煤浆比相同含碳量的煤粉增加了40%的水分,这部分水蒸发消耗的热量使水煤浆着火比煤粉延迟且温度降低[17]。相同入口速度下,水煤浆低温区的径向尺寸比煤粉增大300 mm,轴向脱火距离增大110 mm。煤粉燃烧和水煤浆燃烧在z=0平面的面积分平均温度分别为1 151 K和989 K;煤粉和水煤浆的出口面积分平均温度分别为1 457 K和1 258 K。此外,煤粉火炬对后锥的充满度较水煤浆火炬要好,这是单位能量密度大的表现[18]。

3.3 水分分布

双锥逆喷燃烧器中煤粉和水煤浆燃烧工况下水蒸汽体积分数如图7所示。在前文所述低温区内,水蒸汽浓度最大。在双锥逆喷燃烧器中煤粉的水分高浓度区呈现为三角形;而水煤浆的则为纺锤形,且偏向前锥。该现象主要是2者水分含量上的区别在动量传递上的体现。结合温度的对比分析可知,低温区大小与水分含量一致,水分浓度梯度最大的地方也是温度梯度最大的位置。逆喷的煤粉或水煤浆中的颗粒经过加热,水分汽化为蒸汽。由于水蒸汽和煤粉密度上的差异,两者沿轴向前进的距离不同,因而颗粒可以与O2充分接触。双锥燃烧器速度场的这一特点同样适用于煤粉和水煤浆,这也是影响点火和稳燃的关键因素之一。

  

图6 温度云图Fig.6 Contours of static temperature

  

图7 水分分布云图Fig.7 Contours of H2O

3.4 基于组分分布的火炬特性分析

CO2和O2体积分数可以反映水煤浆或煤粉燃烧过程中火炬特性。由双锥逆喷燃烧器中煤粉和水煤浆燃烧工况下O2体积分数如图8所示。可以看出,水煤浆的火炬锋发散较小,摇曳闪烁较小,火炬更加坚挺。

  

图8 组分云图Fig.8 Contours of species

水煤浆火炬燃烧最剧烈的位置较煤粉提前了438 mm。但其衰减速度也更快,在燃烧器出口位置火炬截面衰减了近50%。结合图6可知,水煤浆火炬核心温度低于煤粉火炬约150℃。由CO体积分数分布可知(图8),水煤浆火炬中CO高浓度区明显较煤粉大。这说明由于水煤浆水分含量高,氧含量相对降低,点火后第1步反应进行较充分,第2步反应明显较煤粉燃烧的情况滞后,所以生成了更多的CO等可燃气体,气体的燃烧速率快于焦炭颗粒。这也说明了水煤浆火炬核心区较煤粉火炬提前且衰减速度更快的原因。通过以上分析可以发现,双锥燃烧器对水煤浆的稳燃效果同样明显。与此同时,火炬核心温度降低和低空气密度的特性对燃烧过程的低NOx排放具有积极意义。

4 结 论

1)水煤浆的燃烧特性处于逆喷双锥燃烧组织控制的有效范围,证明该燃烧组织原理是解决水煤浆燃烧过程中,温度与效率矛盾的有效方法之一。

2)水煤浆和煤粉的速度场类似。水煤浆中水分汽化增大了燃烧器的阻力,使离散相颗粒在燃烧器内的停留时间增加。由于水分含量高,水煤浆着火位置较煤粉延后110 mm。燃用水煤浆的燃烧器内平均温度和出口温度分别比燃用煤粉低162 K和199 K,水煤浆火炬的核心温度区比煤粉提前438 mm,且火炬刚性强于煤粉火炬,但其衰减速度快于煤粉火炬。这是由于水分的增加,氧含量相对降低,使第1步反应生成的可燃气体富集所致。

参考文献(References)：

[1]岑可法,姚强,曹欣玉,等.煤浆燃烧、流动、传热和气化的理论与应用技术[J].燃料化学学报,1997,50(1)：113.CEN Kefa,YAO Qiang,CAO Xinyu,et al.Theory and application of combustion,flow,heat transfer,gasification of coal slurry[J].Journal of Fuel Chemistry and Technology,1997,50(1)：113.

[2]姜秀民,马玉峰,崔志刚,等.水煤浆流化悬浮高效洁净燃烧技术研究与应用[J].化学工程,2006,34(1)：62-65.JIANG Xiumin,MA Yufeng,CUI Zhigang,et al.Study and application of fluidization-suspension high efficiency cleaning combustion technology of coal-water slurry[J].Chemical Engineering,2006,34(1)：62-65.

其四，学校应该营造特色的教师文化。“学校应该有向心力，也就是组织文化。现在教师凝聚力还是比较差一点，学校有学生学习的文化，但是教师文化还是比较欠缺的。例如有的学校重视教科研，全校教师都在参与教科研。有的学校比较重视课外活动，教师们则都在参与活动。例如咱们学校教师逐渐通过做科研获得一些荣誉，变得更加自信、高兴，所以也更乐于参与教科研。”[HDW—TY]

[4]王金华,王乃继.散煤资源清洁利用工程示范——现代煤粉工业锅炉[J].煤炭工程,2016,48(9)：1-5,10.WANG Jinhua,WANG Naiji.Demonstration project of high efficiency pulverized coal fired industrial boiler for bulk coal clean utilization[J].Coal Engineering,2016,48(9)：1-5,10.

[5]胡勤海,熊云龙,朱妙军,等.城市污泥掺制水煤浆燃烧动力学特性[J].环境科学学报,2008,28(6)：1149-1154.HU Qinhai,XIONG Yunlong,ZHU Miaojun,et al.Study on combustion kinetics of coal water slurry prepared with sewage sludge[J].Acta Scientiae Circumstantiae,2008,28(6)：1149-1154.

[6]王乃继,尚庆雨,张鑫,等.德国煤粉工业锅炉实践及我国研发、推广的现状与发展[J].工业锅炉,2016(2)：1-10.WANG Naiji,SHANG Qingyu,ZHANG Xin,et al.Practice of pulverized-coal fired industrial boilers in Germany and statue of its research,promotion in China[J].Industrial Boiler,2016(2)：1-10.

[7]吴玉新,张建胜,王明敏,等.简化PDF模型对Texaco气化炉的三维数值模拟[J].化工学报,2007,58(9)：2369-2374.WU Yuxin,ZHANG Jiansheng,WANG Mingmin,et al.3D numerical simulation of Texaco gasifier using assumed PDF model[J].Journal of Chemical Industry and Engineering,2007,58(9)：2369-2374.

人类乳头状瘤病毒（HPV）感染：HPV是通过性传播感染的，有100多个亚型。对许多人来说，人乳头状瘤病毒不会对人造成伤害，多数不经治疗就可能消失。只有很小比例的HPV感染病人会患上口腔癌或口咽癌，与口腔癌和口咽癌相关的主要是16型。

[8]刘向军,石磊,徐旭常.稠密气固两相流欧拉-拉格朗日法的研究现状[J].计算力学学报,2007,24(2)：166-172.LIU Xiangjun,SHI Lei,XU Xuchang.Research of Euler-lagrange method for dense gas-solid two-phase flow[J].Chinese Journal of Computational Mechanics,2007,24(2)：166-172.

除常规的日常检查之外，学校还重点突击难点问题，比如用“五专”治理食品添加剂的使用，做到专人管理、专柜储存、专人记账、专人领用、专人采购。“我们也在开学季和毕业季等重点时段组织力量全面开展专项服务抽查，规范管理，严格控制无证经营。”该校相关负责人说，今年5月，学校食品安全工作站下发了关于清理校内经营餐饮摊点和出租门面的通知，联合珞南街食药监所突击检查三次，进行摸底排查，对商贸中心等出租门面的餐饮进行集中管理、分类办证，接照生产规范、操作要求进行硬件改造，对教学楼、学生宿舍二次供水系统进行抽查检验，为校园食品安全保驾护航。

[9]张雄,陆明万.任意拉格朗日—欧拉描述法研究进展[J].计算力学学报,1997,14(1)：91-102.ZHANG Xiong,LU Mingwan.Research progress of arbitrary lagrange-Euler description method[J].Chinese Journal of Computational Mechanics,1997,14(1)：91-102.

[10]王月明,曹欣玉,顾战,等.采用反吹射流稳定和控制煤粉和煤浆火炬着火的研究[J].浙江大学学报(工学版),1987(6)：128-138.WANG Yueming,CAO Xinyu,GU Zhan,et al.Study on stabilizing and controlling ignition of pulverized coal and coal slurry torch by reverse jet[J].Journal of Zhejiang University(Engineering Science),1987(6)：128-138.

[11]艾青.热辐射与高速流耦合换热的数值研究[D].哈尔滨：哈尔滨工业大学,2009.

[12]张振奎,钟华贵.多点喷射模型燃烧室性能CFD分析[J].航空发动机,2010,36(4)：17-21.ZHANG Zhenkui,ZHONG Huagui.CFD analysis of performance for the model multipoint Injection combustor[J].Aeroengine,2010,36(4)：17-21.

[13]MAGNUSSEN B F,HJERTAGER B H.On mathematical modeling of turbulent combustion with special emphasis on soot formation and combustion[J].Symposium on Combustion,1977,16(1)：719-729.

[14]BAUM M M,STREET P J.Predicting the combustion behaviour of coal particles[J].Combustion Science&Technology,1971,3(5)：231-243.

[15]MILLER R S,HARSTAD K,BELLAN J.Evaluation of equilibrium and non-equilibrium evaporationmodelsformanydroplet gas-liquid flow simulations[J].International Journal of Multiphase Flow,1997,24(6)：1025-1055.

[16]倪建军,梁钦锋,周志杰,等.气化炉辐射废锅内多相流场和温度场的数值模拟[J].中国电机工程学报,2010,30(5)：59-65.NI Jianjun,LIANG Qinfeng,ZHOU Zhijie,et al.Numerical simulation of multiphase flow field and temperature field in a gasification radiant syngas cooler[J].Proceedings of the CSEE,2010,30(5)：59-65.

[17]周德悟,蔡根才.差热分析测定水煤浆的燃烧特性[J].燃料化学学报,1989(4)：89-94.ZHOU Dewu,CAI Gencai.Determination of combustion characteristics of coal water slurry by differential thermal analysis[J].Journal of Fuel Chemistry and Technology,1989(4)：89-94.

[18]朱瑞,黄定国,吴玉敏,等.新型黑液水煤浆的燃烧特性及动力学分析[J].煤炭转化,2007,30(3)：49-52.ZHU Rui,HUANG Dingguo,WU Yumin,et al.Combustion characteristics and dynamic analysis of a new black liquor coal slurry[J].Coal Conversion,2007,30(3)：49-52.