基于FLUENT的割嘴外部流场火焰燃烧的数值模拟

基于FLUENT的割嘴外部流场火焰燃烧的数值模拟 Numerical Simulation of Flame Combustion at Cutting Nozzle External Flow Field Based on FLUENT first-author 郝　明（1986-），男，讲师，主要从事新型石化设备及火焰切割技术的研发工作。E-mail：haoming198616@163.com。 郝　明（1986-），男，讲师，主要从事新型石化设备及火焰切割技术的研发工作。E-mail：haoming198616@163.com。 haoming198616@163.com 齐齐哈尔大学 机电工程学院，黑龙江 齐齐哈尔 161000 齐齐哈尔大学 机电工程学院，黑龙江 齐齐哈尔 161000 School of Mechanical and Electrical Engineering, Qiqihaer University, Qiqihaer 161000, China School of Mechanical and Electrical Engineering, Qiqihaer University, Qiqihaer 161000, China 利用FLUENT软件对割嘴外部流场的喷雾燃烧过程进行了数值仿真模拟。湍流模型采用标准 k － ε 模型，燃烧模型采用涡耗散（ED）模型，燃油液滴通过离散相模型来描述，尝试利用ED模型探索割嘴外部流场火焰燃烧的可能性，通过模拟得到了合理的温度场、速度场和浓度场，以及火焰燃烧的各种细节。同时，通过对燃料与氧气不同混合比进行对比，发现60%燃料+40%氧气和55%燃料+45%氧气所能达到的最高温度几乎相同，约为2 760 k，进而为提高火焰的燃烧效率提供了理论依据。 In this paper, spray combustion process was simulated numerically at cutting nozzle external flow field by the FLUENT software. Turbulence model was standard model, combustion model was eddy dissipation model, fuel droplet was described by discrete phase model, this article attempts to use the ED model to explore the possibility of flame burning in the flow field outside the cutting nozzle, through simulation, temperature field, velocity field, concentration field and various details of flame combustion reasonable were obtained. At the same time, by comparing the different mixing ratios of fuel and oxygen, it is found that the maximum temperature that 60% fuel and 40% oxygen and 55% fuel and 45% oxygen can reach is almost the same, which is about 2 760 k, and then provide theoretical basis for improve the flame combustion efficiency. 割嘴 火焰燃烧 FLUENT 涡耗散 cutting nozzle flame combustion FLUENT ED 前言 1 火焰燃烧是一种复杂而又迅速的化学反应，是燃料和氧化剂在空间激烈地发生放热化学反应的过程。对于液体燃料，要先将其雾化才能燃烧，而火焰的温度和长度主要取决于燃料与氧气的比例、喷嘴出口的喷射速度以及液滴的初始直径，从而进一步影响切割质量。随着燃烧学与计算流体力学的互相促进与发展，数值模拟在燃烧领域的应用愈加广泛和深入 ［ 1 1 - 5 5 1-5 ］ ，可对液体燃料的湍流雾化燃烧过程进行详细的数学描述 ［ 6 6 ］ ，并进行热态的数值模拟。常用的燃烧模型有涡耗散（ED）模型、涡耗散概念模型、有限速率/涡耗散模型、非预混燃烧模型和组分PDF输运模型，很多学者对此进行了研究 ［ 7 7 - 11 11 7-11 ］ 。Angelo Minotti ［ 7 7 ］ 等采用两个LES模型计算出燃烧室的最高温度并对比了燃烧效率，预测了温度场和不同组分的浓度场；CHENG D等 ［ 8 8 ］ 在研究高速氧气与燃料混合过程中，针对颗粒动量和热能传递建立了一个单向耦合的数学模型，来模拟其对颗粒大小和粒径的影响；Cuoci A等 ［ 10 10 ］ 在研究氧燃料燃烧时，详细描述了燃烧化学的机理，并采用涡耗散模型（EDC）来描述化学-湍流的相互作用，表明模拟结果的可靠性。斛晓飞 ［ 11 11 ］ 等采用FLUENT软件模拟氧-混合气体燃料超音速火焰喷涂（HVOF）过程，分析HVOF的燃烧状态和气体流场，加长喷管对气体流场产生的约束力更大，燃烧室内燃烧充分，同时在喷枪出口处产生了明显的激波。 割嘴结构是影响液体燃料雾化的关键，也是制约液体燃料割嘴发展的主要原因 ［ 12 12 ］ 。本文通过ED模型模拟某氧-乙醇汽油等压式割嘴外部的火焰燃烧状况，研究不同燃料和氧气混合比例对火焰温度变化的影响，从而为进一步优化割嘴尺寸提供理论依据。 计算模型 1 几何建模与网格划分 2 氧-乙醇汽油等压式割嘴立体图如 图1 图1 所示，结构图如 图2 图2 所示。切割氧孔道的入口直径为4 mm，出口直径为1.5 mm；割嘴内芯管的第二层圆台上开有6个预热氧孔道，预热氧孔道的长度取10 mm；燃料孔道与预热氧孔道的夹角为45°，燃料入口直径1.5 mm；该割嘴的外套筒向内缩进到10 mm，内芯管取7.8 mm，以此提高氧气的气动力从而提高乙醇汽油的雾化质量，确保雾滴的燃烧效率和稳定性。 液体燃料与氧气在割嘴内部预混后，经割嘴出口喷出燃料液滴与氧气的混合气，并在一定区域内燃烧。若设定足够大的计算域来实现温度的传递，则计算代价十分高昂，因此给定适当的区域能满足火焰主要部分温度的传递即可；又考虑到割嘴出口由若干分布均匀的条形槽组成，采用GAMBIT建模时考虑计算域的对称性，画出一半即可。割嘴的外套筒域长0.8 m，宽0.2 m，然后生成四边形非结构化网格90 400个，计算域局部如 图3 图3 所示。 模型的建立 2 湍流模型 3 液体燃料液滴与氧气的混合气在割嘴出口会产生湍流现象，用标准 k - ε 湍流方程来描述。标准 k - ε 模型是个半经验公式，主要是基于湍流动能 k 和扩散率 ε ，其公式分别为： <math display="block" id="M1"><mfrac><mrow><mo>∂</mo></mrow><mrow><mo>∂</mo><mi>t</mi></mrow></mfrac><mfenced separators="|"><mrow><mi>ρ</mi><mi>k</mi></mrow></mfenced><mo>+</mo><mfrac><mrow><mo>∂</mo></mrow><mrow><mo>∂</mo><msub><mrow><mi>x</mi></mrow><mrow><mi>i</mi></mrow></msub></mrow></mfrac><mfenced separators="|"><mrow><mi>ρ</mi><mi>k</mi><msub><mrow><mi>μ</mi></mrow><mrow><mi>i</mi></mrow></msub></mrow></mfenced><mo>=</mo><mfrac><mrow><mo>∂</mo></mrow><mrow><mo>∂</mo><msub><mrow><mi>x</mi></mrow><mrow><mi>j</mi></mrow></msub></mrow></mfrac><mfenced open="[" close="]" separators="|"><mrow><mfenced separators="|"><mrow><mi>μ</mi><mo>+</mo><mfrac><mrow><msub><mrow><mi>μ</mi></mrow><mrow><mi>t</mi></mrow></msub></mrow><mrow><msub><mrow><mi>σ</mi></mrow><mrow><mi>k</mi></mrow></msub></mrow></mfrac></mrow></mfenced><mfrac><mrow><mo>∂</mo><mi>k</mi></mrow><mrow><mo>∂</mo><msub><mrow><mi>x</mi></mrow><mrow><mi>j</mi></mrow></msub></mrow></mfrac></mrow></mfenced><mo>+</mo><msub><mrow><mi>G</mi></mrow><mrow><mi>k</mi></mrow></msub><mo>+</mo><msub><mrow><mi>G</mi></mrow><mrow><mi>b</mi></mrow></msub><mo>-</mo><mi>ρ</mi><mi>ε</mi><mo>-</mo><msub><mrow><mi>Y</mi></mrow><mrow><mi>M</mi></mrow></msub><mo>+</mo><msub><mrow><mi>S</mi></mrow><mrow><mi>k</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27427861&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27427860&type= 97.53600311 10.58333302 <math display="block" id="M2"><mfrac><mrow><mo>∂</mo></mrow><mrow><mo>∂</mo><mi>t</mi></mrow></mfrac><mfenced separators="|"><mrow><mi>ρ</mi><mi>ε</mi></mrow></mfenced><mo>+</mo><mfrac><mrow><mo>∂</mo></mrow><mrow><mo>∂</mo><msub><mrow><mi>x</mi></mrow><mrow><mi>i</mi></mrow></msub></mrow></mfrac><mfenced separators="|"><mrow><mi>ρ</mi><mi>ε</mi><msub><mrow><mi>μ</mi></mrow><mrow><mi>i</mi></mrow></msub></mrow></mfenced><mo>=</mo><mfrac><mrow><mo>∂</mo></mrow><mrow><mo>∂</mo><msub><mrow><mi>x</mi></mrow><mrow><mi>j</mi></mrow></msub></mrow></mfrac><mfenced open="[" close="]" separators="|"><mrow><mfenced separators="|"><mrow><mi>μ</mi><mo>+</mo><mfrac><mrow><msub><mrow><mi>μ</mi></mrow><mrow><mi>t</mi></mrow></msub></mrow><mrow><msub><mrow><mi>σ</mi></mrow><mrow><mi>ε</mi></mrow></msub></mrow></mfrac></mrow></mfenced><mfrac><mrow><mo>∂</mo><mi>ε</mi></mrow><mrow><mo>∂</mo><msub><mrow><mi>x</mi></mrow><mrow><mi>j</mi></mrow></msub></mrow></mfrac></mrow></mfenced><mo>+</mo><msub><mrow><mi>C</mi></mrow><mrow><mn mathvariant="normal">1</mn><mi>ε</mi></mrow></msub><mfrac><mrow><mi>ε</mi></mrow><mrow><mi>k</mi></mrow></mfrac><mfenced separators="|"><mrow><msub><mrow><mi>G</mi></mrow><mrow><mi>k</mi></mrow></msub><mo>+</mo><msub><mrow><mi>C</mi></mrow><mrow><mn mathvariant="normal">3</mn><mi>ε</mi></mrow></msub><msub><mrow><mi>G</mi></mrow><mrow><mi>b</mi></mrow></msub></mrow></mfenced><mo>-</mo><msub><mrow><mi>C</mi></mrow><mrow><mn mathvariant="normal">2</mn><mi>ε</mi></mrow></msub><mi>ρ</mi><mfrac><mrow><msup><mrow><mi>ε</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow><mrow><mi>k</mi></mrow></mfrac><mo>+</mo><msub><mrow><mi>S</mi></mrow><mrow><mi>ε</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27427857&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27427863&type= 113.79199982 10.58333302 式中　 ρ 为流体密度； μ 为流体动力粘度； G k 为平均速度梯度引起的湍动能； G b 为由浮力产生的湍动能； μ t 为湍流速度； Y M 为在可压缩湍流中由于过度扩散产生的波动； S k 和 S ε 是自定义的源项，文中均不考虑，取值为0； C 1 ε 、 C 2 ε 、 C 3 ε 为经验常数，取值分别为1.44、1.92、0.99；湍动能 k 和耗散率 ε 的普朗特数 σ k 和 σ ε 分别为1.0和1.3 ［ 13 13 ］ 。 燃烧模型 3 喷嘴外的火焰燃烧过程处于湍流燃烧状态，湍流流动与液体燃料燃烧的化学反应有着密切的联系并互相影响，文中采用湍流-化学反应相互作用模型，即涡耗散模型（ED）来描述火焰的湍流燃烧过程。在Magnussen和Hjertager工作的基础上 ［ 14 14 ］ 提出的ED模型认为，反应速率取决于湍流脉动衰变速率，并能自动选择成分来控制速率，反应 <math id="M3"><mi>r</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27427906&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27427901&type= 1.52400005 2.62466669 中物质 <math id="M4"><mi>i</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27427923&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27427913&type= 1.10066664 2.62466669 的产生速率 <math id="M5"><msub><mrow><mi>R</mi></mrow><mrow><mi>i</mi><mo>,</mo><mi>r</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27427867&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27427878&type= 4.48733330 2.62466669 由 式（3） 式（3） 、 式（4） 式（4） 中较小的一个给出： <math display="block" id="M6"><msub><mrow><mi>R</mi></mrow><mrow><mi>i</mi><mo>,</mo><mi>r</mi></mrow></msub><mo>=</mo><mi>v</mi><msub><mrow><mi mathvariant="normal">'</mi></mrow><mrow><mi>i</mi><mo>,</mo><mi>r</mi></mrow></msub><msub><mrow><mi>M</mi></mrow><mrow><mi>w</mi><mo>,</mo><mi>i</mi></mrow></msub><mi>A</mi><mi>ρ</mi><mfrac><mrow><mi>ε</mi></mrow><mrow><mi>k</mi></mrow></mfrac><munder><mrow><mi mathvariant="normal">m</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi></mrow><mrow><mi>R</mi></mrow></munder><mo stretchy="false">(</mo><mfrac><mrow><msub><mrow><mi>Y</mi></mrow><mrow><mi>R</mi></mrow></msub></mrow><mrow><mi>v</mi><msub><mrow><mi mathvariant="normal">'</mi></mrow><mrow><mi>R</mi><mo>,</mo><mi>i</mi></mrow></msub><msub><mrow><mi>M</mi></mrow><mrow><mi>w</mi><mo>,</mo><mi>R</mi></mrow></msub></mrow></mfrac><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27427870&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27427880&type= 49.02199936 8.63599968 <math display="block" id="M7"><msub><mrow><mi>R</mi></mrow><mrow><mi>i</mi><mo>,</mo><mi>r</mi></mrow></msub><mo>=</mo><mi>v</mi><msub><mrow><mi mathvariant="normal">'</mi></mrow><mrow><mi>i</mi><mo>,</mo><mi>r</mi></mrow></msub><msub><mrow><mi>M</mi></mrow><mrow><mi>w</mi><mo>,</mo><mi>i</mi></mrow></msub><mi>A</mi><mi>B</mi><mi>ρ</mi><mfrac><mrow><mi>ε</mi></mrow><mrow><mi>k</mi></mrow></mfrac><mfrac><mrow><msub><mo>∑</mo><mrow><mi>p</mi></mrow></msub><mrow><msub><mrow><mi>Y</mi></mrow><mrow><mi>p</mi></mrow></msub></mrow></mrow><mrow><msubsup><mo>∑</mo><mrow><mi>j</mi></mrow><mrow><mi>N</mi></mrow></msubsup><mrow><mi>v</mi><msub><mrow><mo>"</mo></mrow><mrow><mi>j</mi><mo>,</mo><mi>r</mi></mrow></msub><msub><mrow><mi>M</mi></mrow><mrow><mi>w</mi><mo>,</mo><mi>j</mi></mrow></msub></mrow></mrow></mfrac></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27427884&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27427882&type= 46.90533447 11.51466751 式中　 Y p 为产物物质的质量分数； Y R 为指定反应物的质量分数； <math id="M8"><mfrac><mrow><mi>k</mi></mrow><mrow><mi>ε</mi></mrow></mfrac></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27427893&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27427889&type= 2.45533323 7.28133297 为大涡混合时间尺度； <math id="M9"><mi>v</mi><msub><mrow><mi mathvariant="normal">'</mi></mrow><mrow><mi>i</mi><mo>,</mo><mi>r</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27427896&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27427894&type= 4.48733330 2.62466669 和 <math id="M10"><mi>v</mi><msub><mrow><mo>"</mo></mrow><mrow><mi>j</mi><mo>,</mo><mi>r</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27427898&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27427897&type= 5.33400011 2.62466669 分别为组分 <math id="M11"><mi>i</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27427923&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27427913&type= 1.10066664 2.62466669 在化学反应 <math id="M12"><mi>r</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27427906&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27427901&type= 1.52400005 2.62466669 中反应物和生成物的化学计量系数； <math id="M13"><mi>N</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27427908&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27427918&type= 2.53999996 2.62466669 为生产物总个数； <math id="M14"><msub><mrow><mi>M</mi></mrow><mrow><mi>w</mi><mo>,</mo><mi>i</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27427911&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27427909&type= 5.84200001 2.62466669 为组分 <math id="M15"><mi>i</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27427923&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27427913&type= 1.10066664 2.62466669 的相对分子质量； <math id="M16"><mi>A</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27427915&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27427914&type= 2.28600001 2.62466669 和 <math id="M17"><mi>B</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27427930&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27427927&type= 2.28600001 2.62466669 均为常数，分别为4和0.5。 燃油颗粒的离散相模型 3 文中用离散相模型来实现液滴与气体的耦合作用。设定割嘴出口为液滴的喷射源，则液滴的作用力在笛卡尔坐标系下的平衡方程为： <math display="block" id="M18"><mfrac><mrow><mi>d</mi><msub><mrow><mi>u</mi></mrow><mrow><mi>p</mi></mrow></msub></mrow><mrow><mi>d</mi><mi>t</mi></mrow></mfrac><mo>=</mo><msub><mrow><mi>F</mi></mrow><mrow><mi>D</mi></mrow></msub><mo stretchy="false">(</mo><mi>u</mi><mo>-</mo><msub><mrow><mi>u</mi></mrow><mrow><mi>p</mi></mrow></msub><mo stretchy="false">)</mo><mo>+</mo><mfrac><mrow><msub><mrow><mi>g</mi></mrow><mrow><mi>x</mi></mrow></msub><mo stretchy="false">(</mo><msub><mrow><mi>ρ</mi></mrow><mrow><mi>p</mi></mrow></msub><mo>-</mo><mi>ρ</mi><mo stretchy="false">)</mo></mrow><mrow><msub><mrow><mi>ρ</mi></mrow><mrow><mi>p</mi></mrow></msub></mrow></mfrac><mo>+</mo><msub><mrow><mi>F</mi></mrow><mrow><mi>x</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27427936&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27427934&type= 51.05400085 8.63599968 <math display="block" id="M19"><msub><mrow><mi>F</mi></mrow><mrow><mi>D</mi></mrow></msub><mo>=</mo><mfrac><mrow><mn mathvariant="normal">18</mn><mi>μ</mi></mrow><mrow><msub><mrow><mi>ρ</mi></mrow><mrow><mi>p</mi></mrow></msub><msubsup><mrow><mi>d</mi></mrow><mrow><mi>p</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup></mrow></mfrac><mfrac><mrow><msub><mrow><mi>C</mi></mrow><mrow><mi>D</mi></mrow></msub><mi mathvariant="normal">R</mi><mi mathvariant="normal">e</mi></mrow><mrow><mn mathvariant="normal">24</mn></mrow></mfrac></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27427939&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27427938&type= 25.99266624 8.89000034 式中　 <math id="M20"><msub><mrow><mi>F</mi></mrow><mrow><mi>D</mi></mrow></msub><mo stretchy="false">(</mo><mi>u</mi><mo>-</mo><msub><mrow><mi>u</mi></mrow><mrow><mi>p</mi></mrow></msub><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27427944&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27427951&type= 15.49400043 3.38666677 为颗粒的单位质量曳力； <math id="M21"><mi>u</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27427953&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27427946&type= 1.86266661 2.62466669 为流体相对速度； <math id="M22"><msub><mrow><mi>u</mi></mrow><mrow><mi>p</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27427958&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27427964&type= 3.04800010 2.62466669 为颗粒相对速度； <math id="M23"><mi>μ</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27427960&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27427966&type= 1.86266661 3.38666677 为流体动力粘度； <math id="M24"><mi>ρ</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27427970&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27427968&type= 1.77800000 3.38666677 为流体密度； <math id="M25"><msub><mrow><mi>ρ</mi></mrow><mrow><mi>p</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27427974&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27427962&type= 2.96333337 3.38666677 为颗粒密度； <math id="M26"><msub><mrow><mi>d</mi></mrow><mrow><mi>p</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27427993&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27427976&type= 3.04800010 2.70933342 为颗粒直径； <math id="M27"><mi mathvariant="normal">R</mi><mi mathvariant="normal">e</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27427996&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27427994&type= 4.14866638 2.62466669 为相对雷诺数。 边界条件与计算方法 2 液体燃料为汽油，割嘴出口处直径为4 mm，采用速度入口给予液滴初始速度为30 m/s，液滴初始直径为800 μm，液体流量为1.5×10 <math id="M28"><msup><mrow/><mrow><mo>-</mo><mn mathvariant="normal">4</mn></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27427981&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=27427998&type= 2.28600001 2.96333337 kg/s，射流成分为60%汽油＋40%氧气；考虑到割嘴在工作时很可能会受到空气流动的影响，所以给予空气0.5 m/s的初速度；出口采用常压压力出口，初始温度为300 K。文中先计算连续相，再将燃料液滴离散于连续相中并进行耦合求解，压力差分格式采用标准离散差分格式，动量方程、湍动能方程和湍动能耗散率均采用一阶迎风差分格式，通过SIMPLE算法耦合求解速度与压力方程。 计算结果与分析 1 割嘴外流场计算域的温度场分析 2 计算域的总温度场分布如 图4 图4 所示。可以看出高温区并不在火焰的焰心，而是随着火焰的扩散集中在焰心四周；由于汽油液滴刚从割嘴出口射出，液滴与氧气的作用并不完全，雾化效果并不明显，所以汽油液滴的燃烧并不完全，温度相对较低；随着燃烧反应的进行，混合气与空气进一步混合，温度也随之增加，最高温度可达2 760 K；由于火焰在空气中燃烧时受到气流的影响，温度有小幅度的摆动。 改变燃料与氧气的混合比例分别为60%汽油+40%氧气、55%汽油+45%氧气、50%汽油+50%氧气，选取距离割嘴出口0.3 m处的截面，温度场对比如 图5 图5 所示。由 图5 图5 可知，由于当氧气含量接近40%时，汽油液滴与氧气已经达到完全燃烧的效果，所以氧气含量为40%和45%时火焰能达到的最高温度（2 760 K）基本相同，而氧气含量进一步提高到50%时，火焰最高温度会低于上述情况，为2 350 K。因此在其他条件不变的情况下，可以通过提高一定的氧气含量来减少燃料含量，从而提高燃烧效率。 割嘴外流场计算域的速度场分布 2 计算域的速度场分布如 图6 图6 所示，0.3 m处速度场对比如 图7 图7 所示，速度矢量图如 图8 图8 所示。由图可知，速度场分布均匀，在燃料液滴行进过程中没有出现回流现象；随着液滴速度的减小，燃烧反应也逐渐减弱。燃料液滴在燃烧尾端处存在较强的回流现象，并且与空气相互作用，致使火焰有向外扩散的趋势。 割嘴外流场计算域各组分的质量分数 2 计算域C 8 H 18 、O 2 和CO 2 的质量分数分布分别如 图9 图9 ~ 图11 图11 所示，其质量分数对比分别如 图12 图12 ~ 图14 图14 所示。可以看出，在C 8 H 18 浓度高处，O 2 浓度较低，这是由于O 2 在其中起到助燃剂的作用，会加速燃烧反应的进行；CO 2 的质量分数与温度场的分布十分相似，这是由于燃烧反应产生CO 2 ，燃烧反应越完全CO 2 所占百分比越大，在一定程度上反映出燃料与氧化剂之间的作用强度，进一步反映出温度场的变化。 结论 1 （1）探索使用涡耗散（ED）模型研究割嘴外部流场火焰燃烧的温度场、速度场和浓度场，得到了合理的温度场、速度场和浓度场；由于文中未考虑辐射对燃烧反应的影响，因此根据计算结果可知，涡耗散模型在一定程度上可以用来模拟割嘴喷射火焰在空气中的湍流燃烧过程。 （2）当燃料与氧气按不同比例混合时，60%燃料+40%氧气和55%燃料+45%氧气所能达到的最高温度几乎相同，约为2 760 K；由于受到空气湍流作用的影响，温度场有向外扩散的趋势，且高温区出现在火焰焰心的四周，在割嘴出口附近火焰温度相对较低。 （3）利用涡耗散（ED）模型在模拟割嘴外部流场喷雾湍流燃烧的过程中，由于该模型认为化学反应作用不明显，而是湍流控制着反应速率，但是大多数火焰的燃烧不能简单地用湍流控制反应速率，且文中未考虑辐射的影响，因而仅能反映出温度场、速度场和浓度场的变化趋势，在下一步工作中，考虑辐射对喷雾湍流燃烧的影响，并对其他几种模型（如组分PDF输运方程、有限速率/涡耗散模型等）进行深入研究。 参考文献 Katsuki Masashi ， Chung Jin-Do ， Kim Jang-Woo ， et al . Development of rapid mixing fuel nozzle for premixed combustion ［J］. Journal of Mechanical Science and Technology ， 2009 （ 23 ）： 614 - 623 . Development of rapid mixing fuel nozzle for premixed combustion T Mahumud ， S K Sangha . Prediction of a Turbulent Non-Premixed Natural Gas Flame in a Semi-Industrial Scale Furance using a Radiative Flamelet Combustion Model ［J］. Flow Turbulence Combustion ， 2010 （ 84 ）： 1 - 23 . Prediction of a Turbulent Non-Premixed Natural Gas Flame in a Semi-Industrial Scale Furance using a Radiative Flamelet Combustion Model A.KLIPFEL ， M.FOUNTI . Numerical Simulation and Experimental Validation of the Turbulent Combustion and Perlite Expansion Processes in an Industrial Perlite Expansion Furnace ［J］. Flow， Turbulence and Combustion ， 1998 （ 60 ）： 283 - 300 . Numerical Simulation and Experimental Validation of the Turbulent Combustion and Perlite Expansion Processes in an Industrial Perlite Expansion Furnace Victor P.Zhukov . The impact of methane oxidation kinetics on a rocket nozzle flow ［J］. Acta Astronautica ， 2019 （ 161 ）： 524 - 530 . The impact of methane oxidation kinetics on a rocket nozzle flow Miroslaw Kowalski . NUMERICAL MODELLING OF COMBUSTION PROCESS WITH THE USE OF ANSYS FLUENT CODE ［J］. Journal of KONES Powertrain and Transport ， 2018 （ 25 ）： 175 - 186 . NUMERICAL MODELLING OF COMBUSTION PROCESS WITH THE USE OF ANSYS FLUENT CODE 周力行 . 湍流两相流动与燃烧的数值模拟 ［M］. 北京 ： 清华大学出版社 ， 1991 ： 8 . 湍流两相流动与燃烧的数值模拟 ZHOU L X . Numerical simulation of turbulent two-phase flow and combustion ［M］. Beijing ： Tsinghua University Press ， 1991 ： 8 . Numerical simulation of turbulent two-phase flow and combustion Angelo Minotti ， Enrico Sciubba . LES of a Meso Combustion Chamber with a Detailed Chemistry Model： Comparison between the Flamelet and EDC model ［J］. Energies ， 2010 （ 3 ）： 1943 - 1959 . LES of a Meso Combustion Chamber with a Detailed Chemistry Model： Comparison between the Flamelet and EDC model D CHENG ， Q XU ， G TRAPAGA ， et al . The Effect of Particle Size and Morphology on the In-Flight Behavior of Particles during High-Velocity Oxyfuel Thermal Spraying ［J］. Metallurgical and Materials Transactions B ， 2001 ， 32 B（ 3 ）： 525 - 535 . The Effect of Particle Size and Morphology on the In-Flight Behavior of Particles during High-Velocity Oxyfuel Thermal Spraying T B GRADINGER ， K BOULOUCHOS . Modeling Turbulent Droplet Motion in Vaporizing Sprays ［J］. Flow， Turbulence and Combustion ， 2000 （ 65 ）： 1 - 29 . Modeling Turbulent Droplet Motion in Vaporizing Sprays A Cuoci ， A Frassoldati ， T Faravelli ， et al . CFD simulation of aturbulent oxy-fuel flame ［J］. Processes and Tec-hnologies for a Sustainable Energy ， 2010 （ 6 ）： 27 - 30 . CFD simulation of aturbulent oxy-fuel flame 斛晓飞 ， 李志勇 ， 李叶 ， 等 . C 3 H 8 +H 2 混合气体超音速火焰喷涂燃烧过程数值模拟 ［J］. 电焊机 ， 2019 ， 49 （ 10 ）： 44 - 49 . C3H8+H2混合气体超音速火焰喷涂燃烧过程数值模拟 HU X F ， LI Z Y ， LI Y ， et al . Numerical simulation of combustion process of C 3 H 8 + H 2 mixed gas supersonic flame spraying ［J］. Electric Welding Machine： 2019 ， 49 （ 10 ）： 44 - 49 . Numerical simulation of combustion process of C3H8 + H2 mixed gas supersonic flame spraying 郝明 ， 胡传顺 ， 蒋应田 . 基于FLUENT的氧乙醇汽油割嘴内部流场分析 ［J］. 热加工工艺 ， 2012 ， 41 （ 17 ）： 201 - 203，213 . 基于FLUENT的氧乙醇汽油割嘴内部流场分析 HAO M ， HU C S ， JIANG Y T . Analysis of internal flow field of oxyethanol gasoline nozzle based on fluent ［J］. Hot working process ， 2012 ， 41 （ 17 ）： 201 - 203，213 . Analysis of internal flow field of oxyethanol gasoline nozzle based on fluent 王福军 . 计算流体动力学分析—CFD软件原理与应用 ［M］. 北京 ： 清华大学出版社 ， 2004 . 计算流体动力学分析—CFD软件原理与应用 WANG Fujun . Computational fluid dynamics analysis-principle and application of CFD software ［M］. Beijing ： Tsinghua University Press ， 2004 . Computational fluid dynamics analysis-principle and application of CFD software HJERTAGER L K ， HJERTAGER B H ， SOLBERG T . CFD modeling of fast chemical reactions in turbulent liquid flows ［J］. Comput. and Chem. Eng. ， 2002 ， 26 （ 4-5 ）： 507 - 515 . CFD modeling of fast chemical reactions in turbulent liquid flows 10.7512/j.issn.1001-2303.2022.01.14 TG482+.3 A 1001-2303（2022）01-0109-06 郝明，王兴国，李腾.基于FLUENT的割嘴外部流场火焰燃烧的数值模拟［J］.电焊机，2022，52（01）：109-114. 郝明，王兴国，李腾.基于FLUENT的割嘴外部流场火焰燃烧的数值模拟［J］.电焊机，2022，52（01）：109-114. HAO Ming, WANG Xingguo, LI Teng.Numerical simulation of flame combustion at cutting nozzle external flow field[J].Electric Welding Machine, 2022, 52(01): 109-114. HAO Ming, WANG Xingguo, LI Teng.Numerical simulation of flame combustion at cutting nozzle external flow field[J].Electric Welding Machine, 2022, 52(01): 109-114. 2021-05-06 2021-07-14 2022-01-20 2022-05-06 电焊机 1 52