一种船用柴油机基于现象学过程的放热率快速预测方法与流程

1.本发明涉及柴油机燃烧性能分析领域，具体涉及一种船用柴油机基于现象学过程的放热率快速预测方法。


背景技术：

2.船用柴油机是世界上各类大型远洋船舶的主要动力装置之一，其性能的高低对航运业的效益有着重大影响。柴油机的放热率是表征燃烧与工作性能的核心参数之一，寻找能够快速预测放热率的方法，对于船用柴油机提高瞬态性能预测能力与运行可靠性有着重要作用。
3.目前运用较多的柴油机放热率预测方法主要包括零维模型预测与cfd模型预测。其中，零维模型搭建简单，预测速度与响应快，但是由于燃烧部分多采用韦伯或者双韦伯等半经验公式进行替代，因此无法反应柴油机内部物理化学过程，所预测的放热率缺少真实性。cfd模型结构详尽，与真实柴油机工作过程最为接近，但是由于计算时考虑的因素过多，计算负荷大，运行效率低，难以对放热率实现快速的预测。


技术实现要素：

4.本发明的目的是为了实现基于真实的喷雾发展与燃烧机理，对放热率进行快速计算，预测工况性能，为船用柴油机瞬态性能预测提供更好的数值分析方法的一种船用柴油机基于现象学过程的放热率快速预测方法。
5.一种船用柴油机基于现象学过程的放热率快速预测方法，包括以下步骤：
6.步骤1:根据柴油机几何参数与工况参数确定边界条件与初始条件，将整个柴油机工作过程从压缩开始到燃烧结束划分为许多小的计算步长，并采用能量守恒方程式计算从压缩到喷油开始每个曲轴转角步长的缸内温度变化
[0007][0008]
式中为一个步长内的缸内传热量，为工质所作体积功，m为工质质量，c
v1
为工质比热容。
[0009]
依照上式，根据初始条件的温度t1求解出第一个步长终点的温度t2，将t2作为下一步长的初始温度迭代计算，得到压缩至喷油开始时刻各步长的温度。已知各步长温度，由理想气体状态方程pv＝rt可以得到相应各步长的压力p。
[0010]
步骤2:喷油开始后，采用下式计算各个时间步长内的喷雾贯穿距s、空气卷吸率ma：
[0011]
[0012][0013]
式中，δp
inj
为喷油压力与缸内压力之差，ρg为缸内工质密度，t为从喷油开始到该计算步长为止的时间，d
noz
为喷孔直径，t为缸内温度，mf为喷油质量流量。
[0014]
步骤3:计算步骤2中喷雾的发展带来的缸内湍动能k的变化：
[0015][0016]
式中m代表喷雾区的总质量，uj为喷孔出口速度，∈为湍能耗散率，采用如下公式计算：
[0017][0018]
式中，lj为缸内湍能尺度，用下式计算：
[0019][0020]
式中，ρf为燃油密度。
[0021]
步骤4:计算每一步长内喷雾的蒸发率，只有蒸发的喷雾才能参与最终燃烧：
[0022][0023][0024]
式中，ra为单位质量燃油的吸热速率，p为缸内压力，rv为燃油蒸发速率，m
liq
为当前步长内液态燃油的质量，e为单位质量燃油蒸发所需的热量。
[0025]
步骤5:采用如下公式计算滞燃期τi：
[0026][0027]
式中，s
p
为当前工况的活塞平均速度，r为气体常数，ea为燃油的活化能，可通过下式计算：
[0028]
ea＝618840/(cn 25)
[0029]
式中，cn为燃油的辛烷值。
[0030]
步骤6:计算步骤4中蒸发燃油的燃烧速率，计算过程如下采用层流-湍流特征时间模型：
[0031][0032]
式中，x代表燃油的质量分数，τc是燃烧的特征时间。特征时间τc由层流特征时间τ
l
与湍流特征时间τ
t
以及延迟系数f组成：
[0033]
τc＝τ
l
fτ
t
[0034]
式中层流特征时间τ
l
采用下式计算：
[0035][0036]
其中a为标定系数，[fuel]为燃油摩尔浓度分数，[o2]为氧气摩尔浓度分数。
[0037]
湍流特征时间τ
t
则是根据步骤3中所计算的湍动能k与湍能耗散率∈计算：
[0038]
τ
t
＝0.1k/∈
[0039]
延迟系数f的计算公式为：
[0040][0041]
式中，r是所有燃烧产物的质量与所有反应物质量的比值，表示为：
[0042][0043]
步骤7:根据步骤6中所计算的燃烧速率，可计算出缸内整体的放热率dq/dt：
[0044][0045]
式中，lhv是燃料的低热值，n
noz
是柴油机喷油器上喷孔的数目。
[0046]
与现有技术相比，本发明的有益效果是：
[0047]
本发明基于现象学建模机理，考虑从喷油时刻开始，喷雾的贯穿、空气卷吸、蒸发以及燃烧，以尽可能简单的公式还原柴油机的真实燃烧过程，实现对放热率的快速预测。根据对喷油条件与柴油机几何参数的调整，该模型预测方法不局限于特定型号柴油机与特定运行工况。除放热率外，该方法还可以对滞燃期、缸压以及空气卷吸率等其他燃烧特征参数进行预测。
附图说明
[0048]
图1是本发明预测方法流程图；
[0049]
图2是本发明针对某型号柴油机缸压与放热率计算结果与实验结果对比图；
[0050]
图3是本发明某型号柴油机参数。
具体实施方式
[0051]
下面结合附图与具体实施方式对本发明作进一步详细描述。
[0052]
针对现有的柴油机放热率预测方法难以反应真实过程、计算速度过慢、瞬态预测性能差等问题，本发明以现象学建模理论为基础，考虑从喷油时刻开始，喷雾的贯穿、空气卷吸、蒸发以及燃烧的过程，提出了基于现象学过程的放热率快速预测方法。该方法计算速度快，精度高，适用范围广，便于对船用柴油机瞬态工作性能的预测与评估。
[0053]
本发明的目的通过以下技术方案实现：
[0054]
步骤1：根据柴油机几何参数与工况参数确定边界条件与初始条件，将整个柴油机工作过程从压缩开始到燃烧结束划分为许多小的计算步长，并采用如下能量守恒方程式计算从压缩到喷油开始每个曲轴转角步长的缸内温度变化
[0055][0056]
式中为一个步长内的缸内传热量，为工质所作体积功，m为工质质量，c
v1
为工质比热容。依照上式，根据初始条件的温度t1求解出第一个步长终点的温度t2，将t2作为下一步长的初始温度迭代计算，得到压缩至喷油开始时刻各步长的温度。已知各步长温度，由理想气体状态方程pv＝rt可以得到相应各步长的压力p。
[0057]
步骤2：喷油开始后，采用下式计算各个时间步长内的喷雾贯穿距s、空气卷吸率ma：
[0058][0059][0060]
上式中，δp
inj
为喷油压力与缸内压力之差，ρg为缸内工质密度，t为从喷油开始到该计算步长为止的时间，d
noz
为喷孔直径，t为缸内温度，mf为喷油质量流量。
[0061]
步骤3：计算步骤2中喷雾的发展带来的缸内湍动能k的变化：
[0062][0063]
式中m代表喷雾区的总质量，uj为喷孔出口速度，∈为湍能耗散率，采用如下公式计算：
[0064][0065]
式中，lj为缸内湍能尺度，用下式计算：
[0066][0067]
式中，ρf为燃油密度。
[0068]
步骤4：计算每一步长内喷雾的蒸发率，只有蒸发的喷雾才能参与最终燃烧：
[0069][0070][0071]
式中，ra为单位质量燃油的吸热速率，p为缸内压力，rv为燃油蒸发速率，m
liq
为当前步长内液态燃油的质量，e为单位质量燃油蒸发所需的热量。
[0072]
步骤5：采用如下公式计算滞燃期τi：
[0073][0074]
式中，s
p
为当前工况的活塞平均速度，r为气体常数，ea为燃油的活化能，可通过下
式计算：
[0075]
ea＝618840/(cn 25)
[0076]
式中，cn为燃油的辛烷值。当缸内曲轴转角超过滞燃期公式计算值时，步骤4中蒸发的燃油开始燃烧。
[0077]
步骤6：计算步骤4中蒸发燃油的燃烧速率，计算过程如下采用层流-湍流特征时间模型：
[0078][0079]
式中，x代表燃油的质量分数，τc是燃烧的特征时间。特征时间τc由层流特征时间τ
l
与湍流特征时间τ
t
以及延迟系数f组成：
[0080]
τc＝τ
l
f
τ
t
[0081]
式中层流特征时间τ
l
采用下式计算：
[0082][0083]
其中a为标定系数，[fuel]为燃油摩尔浓度分数，[o2]为氧气摩尔浓度分数。
[0084]
湍流特征时间τ
t
则是根据步骤3中所计算的湍动能k与湍能耗散率∈计算：
[0085]
τ
t
＝0.1k/∈
[0086]
延迟系数f的计算公式为：
[0087][0088]
式中，r是所有燃烧产物的质量与所有反应物质量的比值，表示为：
[0089][0090]
步骤7，根据步骤6中所计算的燃烧速率，可计算出缸内整体的放热率dq/dt：
[0091][0092]
式中，lhv是燃料的低热值，n
noz
是柴油机喷油器上喷孔的数目。
[0093]
以上仅为对本发明的优选具体实施方式进行的描述，并不限定本发明的构思和范围。在不脱离本发明的创新构思和权利要求所保护范围的前提下，本领域的技术人员基于本发明的技术方案做出的任何修改、变型和改进等，均应落入到本发明的保护范围。本发明请求保护的技术内容,已经全部记载在权利要求书中。