0 引　言

二冲程柴油机在船舶上大量应用，具有低油耗、高效率的特点，有着良好的经济性和动力性。随着全球大气污染的加重和人们环保意识的增强，船舶柴油机的尾气污染引起了大众的广泛关注，IMO对船舶柴油机排放提出了更高的要求^{[1 – 2]}。因此，船用二冲程柴油机的动力性和排放性能成为船舶运营商关注的焦点，故对柴油机的燃烧排放优化提出了迫切的要求^[3]。

环境湿度、环境温度以及环境压力都会对柴油机的燃烧过程产生影响，但是进气湿度对柴油机燃烧和排放的影响具有较大的影响，且其具有难度低、成本小、耐久性好、减排潜力可观的优势，是未来船舶柴油机达到法规的主要技术手段。进气湿度对柴油机燃烧的影响主要在热效应、化学反应、稀释效应等方面。

1）水的高比热容能降低燃烧温度，水相对湿度的增加会提高进气的热容值，进而吸收更多的热量降低燃烧温度，从而减少NOx的产生。

2）燃烧过程中，水的存在会导致“水煤气反应”进而降低Soot的排放，其本质是水受热分解为H基和OH基团，OH基和N的反应较迟钝，且促进Soot的氧化。

3）进气湿度的增大会稀释进气氧浓度，可抑制NOx的产生，但同时会提高Soot的生成^{[4 – 6]}。

本文以GT-Power为工具建立船用柴油机的一维模型，分析不同的进气湿度对柴油机NOx和Soot排放的影响。

1 计算模型介绍与验证 1.1 缸内燃烧的基本方程

1）质量守恒方程

$ \frac{{{\rm d}m}}{{{\rm d}\varphi }} = \frac{{{\rm d}{m_s}}}{{{\rm d}\varphi }} + \frac{{{\rm d}{m_e}}}{{{\rm d}\varphi }} + {g_f}\frac{{{\rm d}x}}{{{\rm d}\varphi }}{\text{。}} $

式中：m为喷入气缸燃油质量；m_S为喷入气缸空气质量；m_e为流出气缸尾气质量；gf，kg/cyc为柴油机的循环喷油量； $\dfrac{{{\rm d}x}}{{{\rm d}\varphi }}$ 为燃烧速率，或者称为放热率。

2）能量守恒方程

由热力学第一定律dq=du+pdv有：

$ \frac{{{\rm d}U}}{{{\rm d}\varphi }} = \frac{{{\rm d}{Q_B}}}{{{\rm d}\varphi }} + \frac{{{\rm d}{m_s}}}{{{\rm d}\varphi }}{h_s} + \frac{{{\rm d}{m_e}}}{{{\rm d}\varphi }}{h_e} - \frac{{{\rm d}{Q_w}}}{{{\rm d}\varphi }} - p\frac{{{\rm d}V}}{{{\rm d}\varphi }}{\text{。}} $

式中：φ为曲轴转角；U为系统的内能；Q_B为燃料在燃烧时放出的能量；Q_W为气缸与外界交换的热量；h_s、h_e为进气阀和排气阀处工质的比焓；p为气缸内工质压力；V为气缸工作容积；

3）动量守恒方程

微元体中流体的动能对时间的变化率等于外界作用在该微元体上各种力之和，由此定律可导出动量守恒方程：

$ \frac{\partial }{{\partial t}} (\rho {\mu _i} + \frac{\partial }{{\partial {x_j}}}\left( {\rho {\mu _i}{\mu _j}} \right) = - \frac{{\partial p}}{{\partial {x_i}}} + \frac{{\partial {\tau _{ij}}}}{{\partial {x_j}}} + \rho {g_i} + {F_i}{\text{。}} $

式中：p为静压；τ_ij为应力张量；g_i和F_i分别为i方向上的重力体积力和外部体积力（如离散相相互作用产生的升力）。

1.2 模型的建立

以MAN B&W 6S35ME-B9型船用二冲程柴油机为研究对象，其技术参数如表1所示。在GT-POWER中建立一个一维的仿真计算模型，燃烧模型选择DIJET燃烧模型进行模拟仿真。

1.3 模型的校核

以6S35ME-B9在100%，75%，50%负荷下对实验值和模拟值进行校核，缸压曲线基本吻合。以75%负荷为例，可得出模拟值和实测值差值都小于5%，在工程允许的误差之内。具体对比如表2和图1所示。

由表2和图1可知，所建立的模型在缸内压力的预测和具体运行工况吻合很好，可用来作为燃烧的具体研究。

2 不同环境湿度对单一因素的影响

柴油机环境湿度会影响柴油机的进气湿度，根据实验目的设置不同的环境湿度。数值模拟的基本条件有：发动机恒定转速129 r/min，单缸循环喷油量为10.285 g，扫气压力2.9 bar，以25 ℃，101 kPa的大气条件下，以相对湿度0%，25%，50%，75%，100%的条件下对柴油机进行模拟仿真，以期充分研究环境湿度对柴油机燃烧和排放的影响规律。

2.1 环境湿度对缸内燃烧温度的影响

图2为不同环境湿度下缸内最高燃烧温度的变化趋势，图3为不同进气湿度下最高燃烧温度局部放大。可以看出，在压缩过程中，缸内压力受进气加湿率的影响较小，但随着燃烧的开始，不同进气湿度出现了温度的分离，随着相对湿度的提高，缸内最高温度和平均温度都有降低的趋势，这是因为在燃烧的过程中，掺水蒸气的空气进入燃烧区与燃油蒸气融合，水分子占去了氧气的比例，减缓了反应的速率；低温的掺水的空气与高温的燃油混合发生热传递会减缓高温燃油的混合传递，从而导致燃烧放热率下降，相位后移；水的高比热容能降低燃烧温度，水相对湿度的增加会提高进气的热容值，进而吸收更多的热量降低燃烧温度，从而减少NOx的产生^[7]。另外，水蒸气会吸收热量分解为O基和OH基参与反应，这也会导致燃烧温度的降低。

2.2 环境相对湿度对缸内燃烧压力的影响

图4为不同相对湿度对压缩压力的影响，图5为不同相对湿度对缸内最高燃烧压力的影响。由图4可以看出在压缩过程中，缸内压力随着相对湿度的增加而增加，这是因为水的密度要小于空气的密度，空气中相对湿度的增加会导致气体的密度降低，进而在压缩终点会导致更多的空气进入到缸内，从而提高了压缩空气的压力^[8]。由图5可以看出，燃烧过程中随着进气相对湿度的提高，最高爆发压力会降低，这主要是由于在燃烧过程中水分子参与了燃烧过程，水分子间存在着缔合分子，水在反应中需要吸收热量克服水之间的缔合作用。另外，水的比热较大，会大量吸收反应中昌盛的热量，从而降低反应的速率，进而降低其压力。

2.3 不同湿度对NOx的影响

影响柴油机NOx生成率和排放率的根本因素是燃烧温度、氧原子浓度和工质在高温环境下的持续时间（即形成NOx的反应时间）。柴油机燃料中的含氮量不足0.02%，排气中氮氧化物主要是由空气中所含氮气在高温下生成的，反应机理为：

$ {{\rm N}_{2}}+{{\rm O}_{2}}\underset{\square }{\mathop{\Leftrightarrow }}\,{\rm NO}+{\rm N}{\text{，}} $

(1)

$ {{\rm O}_{2}}+{\rm N}\underset{\square }{\mathop{\Leftrightarrow }}\,{\rm NO}+{\rm O}{\text{，}} $

(2)

$ {\rm N}+{\rm{OH}}\underset{\square }{\mathop{\Leftrightarrow }}\,{\rm NO}+H{\text{。}} $

(3)

图6为不同进气湿度产生的NOx曲线。可知，随着进气湿度的提高，NOx的生成量呈减少趋势。上述反应中OH原子是水在高温分解时产生的，OH原子的存在会促进NO原子的生成，所以OH原子的浓度是NO原子生成的一个因素。另外，由上述温度可知，相对湿度的增加会导致温度的降低，进而导致NOx产生量减少。由实验可知，温度对NOx的产生具有主导作用^[9]。

2.4 不同湿度对Soot产生的影响

Soot源自未燃燃料。在高温富油区域，由气相燃料核变为固相，随后周围的碳氢化合物或者其余可用分子便在核心上凝结或者被吸收。Particulate则是Soot和其余一些液相或者固相物质的结合物。Particulate分为可溶组分和不可溶组分，而Soot为Particulate的不可溶组分，比例一般要高于50%。

虽然目前还没有一个得到公认的Soot生成反应方程式，但实验已经证明Soot的生成主要受到压力、温度和混合气当量比等因素的影响。Soot产生量的常用计算方法为，以HiroyasuSoot模型和NSC(Nagle and Strickland-Constable)Soot氧化模型为基础，通过Arrhenius公式来模拟Soot的生成量^[10]。

1）HiroyasuSoot生成模型

$ \frac{{{\rm{d}}{{{M}}_{{{Sf}}}}}}{{{\rm{d}}t}}{\rm{ = }}{{{K}}_{{f}}}{{{M}}_{{{fv}}}}{\text{，}} $

$ {{{K}}_{{f}}}{\rm{ = }}{{{A}}_{{f}}}{{{P}}^{0.5}}\sqrt {{\rm{ - }}{{{E}}_{{f}}}/{{RT}}} {\text{。}} $

式中：M_Sf为Soot生成总量；M_fv为燃油蒸汽质量；A_f为前置系数；P为缸内压力；E_f为Soot形成的活化能；R为气体常数；T为气缸内温度。

2）NSCSoot氧化模型

假设Soot颗粒为球形，大小一致，并且Soot颗粒表面由易氧化物质和不易氧化物质组成。

$ \frac{{{\rm{d}}{{{M}}_{{{SO}}}}}}{{{\rm{d}}t}}{\rm{ = }}\frac{{6{{{M}}_{{C}}}}}{{\rho {{{D}}_{{S}}}}}\left( {\frac{{{{{K}}_{{A}}}{{\rm{P}}_{{{\rm{O}}_2}}}}}{{1{\rm{ + }}{{{K}}_2}{{\rm{P}}_{{{\rm{O}}_2}}}}}{\rm{ + }}{{{K}}_{{B}}}{{\rm{P}}_{{{\rm{O}}_2}}}\left( {1 - x} \right)} \right){{{M}}_{{S}}}{\text{，}} $

$ {{{K}}_{{A}}}{\rm{ = }}20\sqrt {{\rm{ - 30000}}/R{{T}}} {\text{，}} $

$ {{{K}}_{{B}}}{\rm{ = 4}}{\rm{.46}} \times {\rm{1}}{{\rm{0}}^{{\rm{ - }}3}}\sqrt {{\rm{ - 15200}}/R{{T}}}{\text{，}} $

$ {{{K}}_{{Z}}}{\rm{ = }}21.3\sqrt {4100/R{{T}}}{\text{。}} $

式中：M_SO为Soot氧化质量；Mc为碳原子的摩尔质量；ρ为Soot密度；D_S为Soot的直径； $P_{O_2} $ 为气缸内氧气压力；x为Soot颗粒表面易氧化物质所占的比例。

$ \frac{{{\rm{d}}{{{M}}_{{S}}}}}{{{\rm{d}}t}}{\rm{ = }}\frac{{{\rm{d}}{{{M}}_{{{Sf}}}}}}{{{\rm{d}}t}}{\rm{ - }}\frac{{{\rm{d}}{{{M}}_{{{SO}}}}}}{{{\rm{d}}t}}{\text{。}} $

由图7可知，随着相对湿度的增加，Soot排量呈现出上升的趋势，且在燃烧过程中，其相对增量增加迅速，其主要原因是进气相对湿度提高导致滞燃期延长，Soot的生成时刻也相对延迟。理论上讲，滞燃期延长可以使油气混合的时间延长，使得局部缺氧、富油区减少，从而减少Soot的形成。另外，水分子在高温下分解的OH基可加速Soot的氧化。但是Soot的后期氧化会随着加湿率的增加而减小，其主要是因为水蒸气占据了原本氧的比例，使得燃油的燃烧后期氧的含量不足，导致Soot的产量增加，最终呈现的是随着相对湿度的增加，Soot的最终生成量增加。

3 结　语

1）随着进气湿度的提高，柴油机缸内压力和燃烧温度均出现下降，缸内最高爆发压力降低的幅度约有12.95 bar，缸内最高燃烧温度降低约210 K。

2）随着进气湿度的提高，可以减少局部富氧区的形成，降低燃烧反应温度，抑制NOx的形成。

3）随着进气湿度的增加，Soot的排放呈上升的趋势，主要因为燃烧后期氧气的不足导致氧化能力的降低。

4）在进一步研究相对湿度对柴油机的影响时应考虑动力性和排放的综合效果，继续研究经济性最优的方法。