建模和模拟高级催化剂,减少非公路车辆排放量

蒂姆瓦特林,约翰逊马瑟


世界各地的监管机构正在对拖拉机、挖掘机和其他非道路柴油动力机械实施越来越严格的排放要求。为了满足这些要求,制造商使用了复杂的后处理系统,其中包括专为非道路车辆设计的催化剂。使用物理原型设计这些系统既昂贵又耗时,不仅因为所有部件都需要制造,而且还因为每个催化剂必须运行很长一段时间,或至少人工老化以代表这段时间,然后才能进行评估。汽车或机器必须在使用寿命结束时达到排放目标,根据法律的定义,使用寿命为3000至10000小时,取决于发动机功率和使用情况。虽然基于实验室的加速老化方法可以将所需时间减少到200小时,但这仍然是相当长的一段时间。

在约翰逊Matthey,我们使用Matlab®和模拟万博1manbetx®在交付原型之前确定最有希望的设计。我们在Simulink中开发了一个完整的后处理系统模型,该模型包含了单个催化剂组分的优化MATLAB模型。万博1manbetx

Matlab和Simulink中的模拟使我们能够了解万博1manbetxCatalyst内发生的复杂交互,执行灵敏度分析,以查看哪些参数对输出具有最大的影响,并根据结果进行设计权衡。通过模拟各种驱动循环的Simulin万博1manbetxk模型,我们可以快速且廉价地评估多种设计选项。我们还使用该模型来系统地检查配置和参数范围以找到最佳设计。结果,我们需要较少的后处理系统原型。

非公路催化剂设计挑战

催化剂用于广泛应用的发动机,包括发电机,乘用车和采矿,农业和建筑设备。虽然催化剂设计的基本原理在应用中保持一致,但优化非公路车辆的催化剂设计引入了一些独特的挑战。

用于非道路车辆的催化剂比用于乘用车的催化剂的产量要小得多,这意味着要制造更少的原型车,从而将开发成本降至最低。此外,乘用车催化剂(图1)可以针对某一特定车辆进行设计和优化,这意味着设计的许多方面,包括催化剂与发动机的距离,都是提前知道的。对于非道路发动机来说,情况就不一样了。

图1.具有外壳切割的汽车排放控制催化剂,以显示内部结构。

一个明显的解决方案是模拟。然而,模拟非道路车辆的催化剂有其自身的困难。与固定发动机(包括备用发电机发动机)不同,非道路车辆上的发动机具有广泛的工作条件。例如,拖拉机拉着犁,可以在田里耕作,也可以在路上行驶。模拟必须考虑流量的变化、温度的变化和其他瞬态变化,以在催化剂运行的所有条件下保持准确性。

MATLAB中的催化剂组分

为了满足排放法规,柴油发动机的完整后处理系统必须取出一氧化碳,未反应的烃,氮氧化物(没有X)和颗粒物。因此,一个完整的Johnson Matthey后处理系统包括柴油氧化催化剂(DOC)、柴油颗粒过滤器(DPF)、氨选择性催化还原(NH)3.SCR)催化剂和氨滑催化剂(ASC)(图2)。

图2。由DOC、DPF、NH3 SCR催化剂和ASC组成的后处理系统示意图。

我们为每个组件创建了MATLAB模型。模型捕获相互关联的物理过程和动力学的复杂组合。物理过程包括气流,以及催化剂内的热量和质量传递。动力学描述了化学反应发生的速率,并展示了速率如何根据温度和气体组成而变化。

为了开发催化剂模型,我们从描述系统的物理学的方程开始,包括气体和固体(催化剂)相的能量和质量平衡,以及描述这些相之间的热量和质量传递的方程。然后,我们在实验室运行实验,使我们能够准确测量催化剂的输出,同时精确地控制输入和催化剂参数。例如,我们测量一氧化碳转化为各种气体混合物的温度函数(图3)。

图3.表明一氧化碳氧化如何随温度变化各种气体混合物的曲线。点代表测量数据;线条,模拟数据。

为了优化模型精度,我们使用来自全局优化工具箱的遗传算法求解器将速率方程的参数拟合到测量数据。在构建催化剂组件的原型之后,我们验证模型的输出,针对从实际组件取出的测量,并根据需要调整模型。

每个MATLAB组件都在SIMULINK库中实现为S函数块(图4)。万博1manbetx除了催化剂组分模型之外,文库还包括用于排气管的热损失模型,双表皮管的热损耗模型和进料块。基于监管机构使用的驱动循环,进料块为Simulink后处理系统模型提供气流,温度和其他输入,包括非公路瞬态周期(NRTC)。万博1manbetx当执行驱动周期时,我们通过从真实柴油发动机捕获引擎排气数据来获得进料块的数据。

图4.柴油发动机催化万博1manbetx剂组件的Simulink库。

模拟完整的后处理系统

我们从催化剂文库块迅速组装完整的后处理系统的型号(图5)。这需要几分钟,远远少于建立真实系统的时间。我们可以通过设置其长度,直径,初始温度,初始烟灰加载(用于滤波器型号),贵金属载荷和其他参数来配置任何块。

图5.图2中显示的后万博1manbetx处理系统的Simulink模型。

我们在Simulink中运行模拟,以评估各种万博1manbetx系统配置和参数的任何给定驱动周期的效力。我们可以在链中的任何点检查中间输出。例如,我们可以绘制第一阶段的模拟一氧化碳和总烃(THC)排放,并将结果与​​测量数据进行比较以验证模型的阶段(图6)。

图6.绘图对CO(顶部)和THC(底部)的模型预测进行比较测量催化剂输出的图。减排作为非公路瞬态周期的累积排放。

在某些情况下,我们的客户的设计要求是灵活的 - 例如,它们可以通过将其更靠近发动机或通过改变发动机的校准来提高催化剂的入口温度。为了评估设计替代方案,我们运行多个模拟驱动周期,改变每个模拟的入口温度,并绘制结果(图7)。然后,客户可以对放置催化剂的位置做出明智的决定。通常,我们自动化多种模拟使用MATLAB脚本运行,以编程方式调整每个运行的Simulink模型中的关键参数,启动模拟,并捕获结果进行分析。万博1manbetx

图7。图显示改变催化剂入口温度对氨SCR催化剂NOX排放的影响。

中间结果对于验证是有用的,但是我们最感兴趣的是排气管的输出,对于图5所示的模型,它是ASC的输出。通过模拟,测量了CO、THC和NO的累积量X排放,以及NH3.滑动,在尾部,评估后处理系统的整体有效性(图8)。

图8.由Doc + DPF + SCR组成的完整后处理系统的排放结果:累积NOx(顶部)和NH3滑动(底部)。

当我们建立一个原型,我们比较它的测量输出与仿真输出来验证模型。然后,我们可以使用模型和仿真来微调原型,在它投入生产之前。

为什么我们选择Matlab通过自定义过程建模包

在使用MATLAB和SIMULINK到型号催化剂万博1manbetx之前,Johnson Matthey的工程师试图使用商业软件包来开发自定义过程模型。我们使用此包开发的模型不足以处理我们定期遇到的场景。例如,溶解器通常足以用于稳态条件和恒定温度,但不能处理我们处理的输入的瞬态性质,包括将气体混合物的温度范围和进入催化剂的气体混合物的变化。使用此包,我们无法更改源代码以使其更准确或克服诸如模拟无法收敛的问题。

相比之下,使用MATLAB,我们可以编写自己的方程和算法,使我们能够完全控制整个模型。我们确切地知道模型是如何工作的,并且可以很容易地确定模型输出和来自真实催化剂的测量数据之间的任何差异的来源。在Simulink系统级模型中集成组件和运行基于时间的仿真的能力节省了时间和成本。万博1manbetx在MATLAB和Simulink中开发我们自己的系统的另一个优势是,我们可以获取Johnson Matthey工程师积累的万博1manbetx组织知识和专业知识,而不是依赖于另一家公司的一刀无垠的解决方案。

发布2012年 - 92054V00