技术文章和通讯

降低非道路车辆排放的先进催化剂建模与模拟

蒂姆·沃特林(Tim Watling)、约翰逊·马修(Johnson Matthey)著


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

在Johnson Matthey,我们使用MATLAB®和仿真软金宝app件®在确定原型之前确定最有希望的设计。我们在Simulink中开发了一个完整的后处理系统模型,其中包含了对单个催化剂组件的优化MATLAB模型。金宝app

MATLAB和Simulink中的仿真使我们能够了解金宝app催化剂内发生的复杂相互作用,进行灵敏度分析,以确定哪些参数对输出影响最大,并根据结果进行设计权衡。通过模拟各种驱动循环的Simulink模型,我们可以快速且廉价地评估多种设计选项。我们还使用该模型系统地检查配置和参数范围,以找到最优设计。因此,我们需要的后处理系统原型要少得多。

非道路催化剂设计挑战

催化剂广泛应用于发动机,包括发电机、客车、采矿、农业和建筑设备。虽然催化剂设计的基本原则在不同的应用中保持一致,但优化非道路车辆的催化剂设计带来了一些独特的挑战。

非道路车辆的催化剂比乘用车的催化剂生产的数量要少得多,这意味着制造更少的原型,以最大限度地降低开发成本。此外,乘用车催化剂(图1)可以针对特定车辆进行设计和优化,这意味着设计的许多方面,包括催化剂与发动机的距离,都是提前知道的。这不是非公路发动机的情况。

图1所示。汽车排放控制催化剂与外壳切割显示内部结构。

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

催化剂组分的MATLAB建模

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

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

我们为每种成分创建了MATLAB模型。这些模型捕捉了相互关联的物理过程和动力学的复杂组合。物理过程包括气体流动以及催化剂内的传热和传质。动力学描述了化学反应发生的速率,并显示了速率如何随时间变化ng与温度和气体成分有关。

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

图3。不同气体混合物的一氧化碳氧化随温度变化的示意图。这些点代表测量数据;这些线,模拟数据。

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

每个MATLAB组件在Simulink库中实现为一个S功能块(图4)。除催化剂组件模型外,该库还包括排气管的热损失模型、双层皮管的金宝app热损失模型和进料块。进料块根据监管机构使用的驱动循环(包括非道路瞬态循环(NRTC))向Simulink后处理系统模型提供气体流量、温度和其他输入。我们通过在实际柴油发动机执行驱动循环时捕获发动机排气数据来获取供油块的数据。

图4.柴油发动机催化金宝app剂组件的Simulink库。

模拟完整的后处理系统

我们从催化剂库块快速组装了完整的后处理系统模型(图5)。这只需要几分钟,比构建真实系统的时间要短得多。我们可以配置任何块通过设置其长度,直径,初始温度,初始烟灰装载(一个过滤器模型),贵金属装载,和其他参数。

图5。金宝app后处理系统的Simulink模型如图2所示。

我们在Simulink中运行仿真,以评估任何金宝app给定驱动周期的各种系统配置和参数的有效性。我们可以检查链上任何一点的中间输出。例如,我们可以绘制第一阶段的模拟一氧化碳和总碳氢化合物(THC)排放量,并将结果与测量数据进行比较,以验证模型的该阶段(图6)。

图6。图中比较了DOC的催化剂产量与CO(上)和THC(下)的模型预测。排放表现为非道路瞬变循环的累积排放。

在某些情况下,我们客户的设计要求是灵活的,例如,他们可以通过将催化剂移近发动机或改变发动机的标定来提高催化剂的进口温度。为了评估设计备选方案,我们对驱动循环进行了多次模拟,改变每次模拟的入口温度,并绘制结果(图7)。然后,客户可以就催化剂的放置位置做出明智的决定。通常,我们使用MATLAB脚本自动化多个仿真运行,该脚本以编程方式调整每次运行的Simulink模型中的关键参数,启动仿真,并捕获结果进行分析。金宝app

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

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

图8。由DOC + DPF + SCR组成的完整后处理系统的排放结果:累积NOX(上)和NH3滑移(下)。

当我们构建原型时,我们将其测量输出与仿真输出进行比较,以验证模型。然后,我们可以使用模型和模拟来调整原型,然后再投入生产。

为什么我们选择MATLAB而不是定制的过程建模软件包

在使用MATLAB和Simulink对催化剂建模金宝app之前,Johnson Matthey的工程师尝试使用商业软件包来开发定制过程模型。我们用这个包开发的模型不够灵活,无法处理我们经常遇到的场景。例如,求解器通常足以满足稳态条件和恒定温度,但不能处理我们处理的输入的瞬态性质,包括进入催化剂的气体混合物的温度范围和变化。有了这个包,我们无法更改源代码以使其更精确,也无法克服模拟无法收敛等问题。

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

发布于2012年- 92054v00