搞了几年动力总成仿真说实话内燃机这块儿是越做越有意思也越做越觉得水深。尤其到了排放法规和热效率双重压力的今天燃烧系统优化几乎成了每个整机厂和设计公司必须啃的硬骨头。单纯靠台架试错时间和费用都扛不住所以现在的主流打法就是燃烧仿真配合台架验证。今天我就以内燃机燃烧系统优化为切入点聊聊我在实际项目里是怎么用CFD仿真工具把燃烧系统一步步调优的会涉及到模型怎么选、网格怎么处理、参数怎么标定、以及各种让人头大的收敛和失真问题全是实操里趟过坑之后沉淀下来的东西希望能让正在做或者准备做这块的朋友少走点弯路。1. 燃烧系统优化的整体思路与仿真方案选型1.1 先把仿真目标想清楚再谈软件和模型很多人拿到一个内燃机项目就急着开网格、点求解器其实这是最容易翻车的开始。燃烧系统优化的第一件事不是在软件里倒腾而是把优化目标量化出来。比如我们这次的项目是针对一台四缸增压直喷汽油机做燃烧系统优化目标很明确在保证输出的前提下把中低负荷2 bar BMEP到8 bar BMEP区间的BSFC降低3%-5%同时把NOx排放控制在法规限值以内。为什么要先锁定这些因为燃烧仿真的计算资源和时间成本都很高一个三维瞬态整机循环动辄几百万网格、几十个小时的机时如果没有明确目标很容易陷入“算了一堆工况最后不知道该拿哪一版去改缸盖”的困境。目标定义好之后还要梳理边界条件。边界条件不齐算出来的结果连趋势参考价值都不大。比如增压器的特性曲线、进排气道流量系数、喷油器的喷雾特性数据粒径分布、喷雾锥角、贯穿距、缸内壁面温度边界这些都是燃烧仿真的命根子。其中喷油器的喷雾数据建议一定要有试验标定的输入而不是随便拿一个经验分布去蒙因为直喷汽油机的燃烧对喷雾形态极其敏感油滴打在缸壁上直接导致机油稀释和PN排放恶化。我的经验是优化的起点和终点都要落到“参数矩阵”上工况点矩阵转速、负荷、边界条件矩阵进气温度、压力、EGR率、喷油时刻、点火时刻、喷油压力、评价指标矩阵缸压峰值及相位、放热率曲线、压力升高率、未燃碳氢、NOx、碳烟、燃油消耗率。这个矩阵建好了后面的仿真和试验验证都是在往这个矩阵里填数才不会东一榔头西一棒子。1.2 主流的燃烧仿真工具与选型建议燃烧仿真工具圈子里常见的无非是CONVERGE、ANSYS Forte原FIRE、STAR-CCM、FLUENT以及学术圈常用的OpenFOAM搭配cantera做化学动力学。这几个工具各有脾气选型的时候不能只盯着名气要看你手头的团队能力、预算和具体需求。如果做汽油机尤其是直喷点燃式的燃烧优化我自己用得最顺的是CONVERGE。它的优势在于自动网格生成和自适应网格加密AMR你不需要手动去画一堆质量很高的体网格软件会基于流场梯度自动局部加密。我在项目中处理进气门落座、喷油嘴附近液滴破碎、火花塞点火核形成这些关键区域时AMR能精准地把网格加密在需要的空间位置省掉了大量手工修网格的时间。而且CONVERGE内置了一整套经过标定的G方程燃烧模型和SAGE详细化学反应求解器对汽油、柴油、甲醇、天然气等燃料都有成熟的化学机理适配项目启动速度快。如果是柴油机Forte和CONVERGE都做得很好关键是喷雾模型和燃烧模型的配合如果是进气道喷射的PFI汽油机FLUENT或者STAR-CCM配合涡耗散概念模型EDC也不是不能做但精度和调参成本会高一些。另外OpenFOAM是所有工具里最灵活的但也是最折腾人的需要自己写求解器和改代码如果团队里没有专门的CFD开发人员我不建议在项目周期紧张的时候用它做整机燃烧系统优化。算得慢不说一个湍流模型参数设置不对后处理阶段要返工整个开发节奏就垮了。表不同燃烧仿真工具在动力总成领域的适用场景对比工具核心强项相对弱项适合场景CONVERGE自动网格/AMR、喷雾与燃烧模型集成度好授权费用较高直喷汽油机、柴油机三维燃烧优化ANSYS Forte与STAR-CCM生态协同好、化学反应求解快网格前处理工作量偏大柴油机喷雾燃烧、排放预测STAR-CCM网格工具全面、后处理强大燃烧模型库略单一冷却系统、进排气道流动与燃烧耦合FLUENT通用性好、资料多瞬态点火燃烧标定成本高PFI汽油机、基础流动计算OpenFOAM免费、灵活、可深度定制工程化程度低、学习曲线陡高校研究、特定机理开发我个人的项目习惯是梯次配置快速评估用零维/一维整机模型GT-Power或者AMESim把进排气相位、气门升程曲线、EGR率、涡流比这些“大趋势参数”先在循环层面扫一轮筛选出几个有希望的配置组合再用三维CFD精细仿真选定的组合去重点看清楚缸内流动、喷雾撞壁、火焰传播和排放生成的具体细节。零维/一维和三维结合能最大化利用有限的仿真资源。2. 燃烧系统核心细节拆解与建模要点2.1 几何模型处理与计算网格策略从STL到能算的网格燃烧系统的三维仿真几何是不可绕开的第一步。很多新人拿到的CAD模型是带完整水套、螺栓孔、密封面的“全套整车件”恨不得几千万个面都扔进去算这完全没必要。做缸内燃烧模拟只需要进气道的下游段、排气道的上游段、燃烧室活塞顶面缸盖底面缸套壁面这几部分。我建议拿到CAD后先做“瘦身”用SpaceClaim或者ANSA把水套、凸台、螺栓孔、倒角等对缸内流动和喷雾没有直接影响的结构全部拍平保持流道的光顺。然后对气门密封面、火花塞电极、喷油器喷嘴尖端这几个关键区域单独做局部加密区域控制因为这些地方的几何细节会显著影响气流和喷雾。网格策略上我强烈建议采用“基网格加密区控制瞬态AMR”三层思路。基网格尺寸一般设置在2-4 mm就够用因为缸内的大尺度涡流是米量级的基网格太细只会白白消耗算力。加密区主要设置在三处一是进气门座附近气门开启时的高速射流区网格加密到0.5 mm左右二是火花塞间隙和点火区域点火核的形成需要高分辨率加密到0.25 mm三是喷嘴出口到喷雾贯穿的路径上加密到0.5 mm同时配合AMR追踪液滴前锋。AMR的判据一般是速度梯度和温度梯度梯度阈值设得太低会导致网格量暴增设得太高又捕捉不到火焰前锋。我常用的是温度梯度阈值在2-5 K/cm之间速度梯度阈值在1-2 m/s/cm之间起步在此基础上根据收敛情况微调。有没有必要做全周期720°CA瞬态仿真看你要解决什么问题。如果关心的是冷机起动、怠速稳定性或者进气道内残留废气对燃烧的影响那从进气门开启前就开始算整个循环计算代价很大。针对我们的优化目标中低负荷的BSFC和排放主要关心的是压缩上止点附近的流动和燃烧过程所以我通常从进气门关闭时刻比如-140°CA ATDC开始计算到排气门开启前结束大概300多度曲轴转角能在准确性和计算成本之间取得平衡。2.2 燃烧模型与湍流模型选择算得准的核心燃烧系统仿真最核心的部分就是燃烧模型没有之一。燃烧模型选错了后面所有参数都是白调。在三维内燃机CFD中常用的燃烧模型有G方程模型火焰面密度、EDC模型涡耗散概念、SAGE详细化学反应求解器、以及PDF输运方程模型。对点燃式汽油机G方程详细机理的组合普遍兼顾了计算速度和点火精度。G方程的核心思路是把火焰前锋看作一个传播的标量场通过湍流火焰速度模型来计算火焰前锋的传播不需要解析每一根反应链所以算得快。实际项目中我钟情于CONVERGE里的SAGE求解器搭配正庚烷/异辛烷/甲苯等组分的简化机理比如用92种组分、256步反应的PRF机理来描述汽油替代燃料。为了让燃烧预测更加贴近实际汽油我还会在机理里掺入少量乙醇组分模拟E10汽油。SAGE求解器不需要指定着火位置它会根据局部的温度、压力、组分浓度和湍流混合时间自动判断着火条件这就特别适合直喷汽油机的分层燃烧和火花点火后的火焰传播问题。它对网格和时间步长的要求较高但在AMR配合下完全能接受。湍流模型我一般用RANS框架下的RNG k-epsilon变种而非标准k-epsilon因为缸内强压缩、强旋转、剪切流场中RNG k-epsilon对旋流和绕流的预测明显更准确而且计算量比LES小了一个量级。LES在学术上看起来很美好能解析出丰富的涡结构但放在工程优化项目里它对于边界条件和网格分辨率极其敏感动不动就要千万级网格一个工况算一周根本不适用于需要快速迭代的燃烧系统优化。除非是研究爆震机理或者近零排放的超稀薄燃烧否则工程应用上别轻易踩进LES这个坑。2.3 喷雾与壁面相互作用直喷机绕不开的鬼门关直喷汽油机的喷雾过程是整个燃烧仿真里最难伺候的部分。喷雾模型通常包括破碎、蒸发、碰撞聚合、湍流扩散和壁面相互作用这几个环节。液滴破碎模型我常用KH-RT模型因为它在高喷射压力35 MPa以上下对初次破碎和二次破碎的模拟都比较靠谱。在35 MPa喷射压力下喷孔出口处的液滴速度可以超过200 m/s破碎时间和贯穿距离都直接影响油束与缸壁、活塞顶的撞壁程度。如果喷油策略是多次喷射比如预喷主喷后处理喷KH-RT模型也能较好捕捉每段油束的发展。壁面相互作用模型直喷机务必要用心调。因为喷雾撞击缸壁和活塞顶是常态油膜liquid film的形成直接关联到机油稀释、颗粒物排放和积碳。目前主流模型是ORourke的壁面相互作用模型它能模拟液滴反弹、铺展、贴壁和飞溅。我在项目中调整壁面相互作用的参数时有个经验先用喷雾特性试验数据在常压或背压环境下用光学发动机或者定容弹拍到的喷雾贯穿距和锥角照片来标定喷雾模型把喷雾的宏观形态调到和试验吻合再耦合进整机燃烧仿真。千万别跳过这一步直接整机算不然喷雾穿距差30%燃烧相位全错缸压曲线一定会给你颜色看。3. 实操流程与关键环节落地3.1 仿真边界条件的确定从一维到三维的数据传递燃烧仿真的边界条件不是拍脑袋出来的尤其是对四缸增压机进排气系统之间存在强烈的干涉效应各缸的进气量、残余废气系数和进气道温度逐缸都不一样。我建议先用GT-Power建立整机的一维热力学模型并通过台架数据倒拖工况的缸压、全负荷的性能点校核一维模型的标定。然后在一维模型里算出各缸在指定工况点的进气流量、进气温度、进气道压力、EGR率、残余废气系数再把进气道的压力脉动曲线作为三维仿真入口的瞬态边界条件导入。这里有个细节容易被忽略单向耦合时一维模型给三维模型的是“平均意义”上的边界但实际发动机在进排气过程中压力和温度波动非常大尤其在气门重叠期的扫气阶段。因此三维模型最好采用瞬态压力边界把一维模型计算出来的曲轴转角分辨率的压力曲线直接加载到进出口边界面上这样进气涡流、滚流和排气背压的脉动效应才能真实被再现。3.2 燃烧参数标定与优化迭代DOE设计及后处理指标当基准工况模型调试到与试验缸压曲线对得上缸压峰值误差在5%以内峰值相位误差在2°CA以内就可以开始正式的优化迭代。优化迭代绝不是无脑穷举我强烈建议引入试验设计DOE方法。拿我们的项目举例在8 bar BMEP、2000 rpm的典型工况下我们选取了三个变量作为优化因子点火提前角在MBT附近±5°CA、喷油时刻分主喷时刻和预喷比例、EGR率10%-20%。用中心复合设计生成了大约20个仿真工况。每个工况大约需要8-12小时的机时20核并行、120万网格所以总共也就200多小时两周内可以跑完。这比起盲目地单点试算效率提升了数倍。计算完成后后处理只看两个最重要的指标曲线缸压曲线和放热率曲线。但看曲线时不能只看峰值要重点看CA50累计放热50%对应的曲轴转角有没有落在合适的位置一般在8°CA ATDC到10°CA ATDC之间比较多看10%-90%燃烧持续期是否合适太短容易爆震太长热效率偏低。排放方面需要监测NOx生成峰值和碳烟生成的时空分布NOx主要出现在高温、富氧区域的火焰前锋后面所以调整EGR率和点火时刻就是在找碳烟与NOx此消彼长的平衡点。3.3 后处理分析驱动设计改进火要是怎么看的三维仿真最大的优势不是算出几条缸压曲线而是能可视化地看到燃烧室里的“火”是怎么样传的。我在后处理时习惯生成一组剖面云图在火花塞中心位置做一个横切面和纵切面看不同曲轴转角下的温度场和OH自由基分布。OH自由基层是高温火焰锋面的标志物观察它是怎么从火花塞间隙出发、如何绕过气门坑边缘、如何被滚流拉伸扭曲就能判断燃烧室的滚流比设计是否合理。如果火焰前锋到达气缸远端的时刻太晚说明滚流强度不足或者火花塞位置偏置这时候可能需要优化进气道造型、气门导流罩、活塞顶面形状比如加深挤气坑来增强滚流。墙面的热负荷评估也不能忽略。燃烧仿真得到的壁面热流密度分布图可以直接指导冷却水套的设计和铸造壁厚的优化。温度分布特别高的局部区域比如排气侧气门鼻梁区容易导致热疲劳开裂和早燃。我习惯把热流密度都输出到几何对应的表面网格上然后在STAR-CCM或者ParaView里叠加显示一目了然。4. 常见问题排查与踩坑实录4.1 发散与不收敛先压时间步再查网格质量燃烧仿真最常见的故障就是计算发散。现象一般是某个变量温度、压力或速度在迭代几个时间步后溢出NaN或者剧增。我排查的顺序很固定。先看时间步长燃烧化学反应特征时间极短尤其是在着火时刻附近一个时间步推进太多就会导致反应项过冲。点火期间的固定时间步长我通常控制在0.01°CA甚至更小如果用了自适应时间步就把最大时间步长限制在0.05°CA。如果时间步已经压得很小还发散第二大嫌疑是网格质量特别是AMR加密后有一些尺寸过小的网格或者在气门座、活塞环岸附近出现负体积这种网格问题可以通过检查最大扭曲度和最小正交性来筛除。第三大嫌疑是边界条件的突变比如喷油压力曲线如果瞬间从0跳到35 MPa动量源项太猛也会发散此时需要给喷射压力加一个极短时间的斜坡过渡。4.2 喷雾与燃烧预测不准试验标定不可省喷雾贯穿距预测过短或过长是最常见的模型标定问题。贯穿距预测过短一般是因为KH-RT破碎模型的破碎时间常数B1设得过小液滴破碎过快、粒径过小导致动量快速衰减。相反如果贯穿距预测过长则B1设得过大液滴碎不起来一直以很大的惯性往前冲容易撞壁。我建议B1的标定范围在7到30之间变化尤其针对汽油直喷的多孔喷油器并在定容弹喷雾试验数据的基础上通过试算2到3组参数快速拟合。另外环境密度对贯穿距影响很大喷雾仿真一定要在正确的背压下进行在低背压下标定好的参数直接拿去高背压工况大概率是要失真的。燃烧相位对不上试验那是另一个大坑。比如预测的着火时刻比试验晚了好几度CA可能的原因有点火能量模型设置过低、点火持续期过短或者是湍流火焰速度公式中某个常数比如火焰拉伸系数不合适火焰传播速度被低估了。如果是CA50偏早则可能是EGR率在边界条件中给低了缸内残余废气的稀释效果没体现出来着火延迟缩短。我遇到的最诡异的一次是后处理时发现CA50对点火时刻的变化完全不敏感排查了几个小时才发现是火花塞附近网格太粗点火核区域内根本解析不出高温火核火焰传播从一开始就建立不起来。后来把火花塞间隙附近的AMR等级强制提高两级问题立刻消失。这个案例说明模型是否失真的要结合网格分辨率和物理模型一起排查别只盯着一个维度。4.3 计算精度与效率的平衡教你几招节省机时做工程项目的都要讲究效率。燃烧仿真动辄数天的机时对项目排期是巨大压力。我常用的省机时技巧有几个。第一化学反应机理精简。详细的汽油机机理可能有上千步反应在三维CFD中每个网格每一步耦合求解化学反应代价极大。我一般先用机理简化工具比如ANSYS的Response Surface分析、或者CONVERGE自带的机理简化模块把机理精简到100步以内同时用激波管和层流火焰速度数据验证一下简化机理在常温常压和高温高压下的误差。只要着火延迟误差在10%以内就足以支撑燃烧趋势预测。第二分阶段网格策略。在压缩冲程初期缸内流场相对均匀AMR等级可以放低到接近上止点时再把AMR等级迅速提高这样才能既保证精度又控制整体网格规模。很多时候我直接在喷油开始前2°CA设置一个事件触发器用一个更高精度的AMR网格方案去覆盖压缩上止点附近的所有时刻。第三多核并行与动态负载平衡。在集群上运行时经常遇到各核负载不均特别是AMR局部加密后有的核网格量是其他核的几倍导致大量核空闲等待。打开动态负载均衡选项能有效改善并行效率尤其是气体喷射和燃烧这两段最容易出现负载倾斜的阶段。实测下来开和不开动态负载均衡在192核的集群上总耗时差别能达到30%以上这是很多人容易忽略的隐性成本。4.4 常见问题速查表故障现象可能原因排查与解决建议计算在点火时刻附近发散时间步长过大、点火核心区域网格过粗点火阶段固定时间步长降至0.01°CA局部AMR等级2缸压峰值远高于试验壁面温度设置过低、湍流强度过高核对壁温边界缸盖/活塞/缸套分别为450-550 K调低入口湍动能缸压峰值远低于试验燃烧模型火焰速度被低估、EGR边界值偏大标定层流火焰速度系数核对EGR率和残余废气系数NOx预测值系统性偏高NOx生成机理过于详细或温度边界偏高检查热力NOx模型参数、峰值燃烧温度是否超过2500 K喷雾贯穿距不足B1破碎常数偏小、背压设置偏高在定容弹数据下标定B1核对背压条件碳烟排放预测异常碳烟前驱体模型参数不合适、局部过浓区域捕捉不足检查局部当量比2的区域考虑增加近喷嘴区网格分辨率5. 优化效果验证与工程经验延伸5.1 从仿真到台架的闭环验证燃烧仿真做得再漂亮最后也要拿到台架上去验收。我在仿真优化完成后会在台架上重点验证两个内容一是缸压曲线对比看实际台架的缸压峰值、CA50和放热率曲线能不能复现仿真的趋势二是排放对比测NOx、THC、CO和PN看排放是否和仿真预测同向变化。在进入台架验证前我还会把仿真优化出的几组边界条件参数点火时刻、喷油时刻、EGR率、喷油压力直接生成一版ECU标定建议MAP给标定工程师这样可以大大压缩台架标定时间。在我最近一次做增压直喷汽油机燃烧优化的项目中仿真给出的方案是将主喷时刻从-35°CA ATDC推迟到-28°CA ATDC以减少燃油撞壁同时把EGR率从12%提到16%泵气损失降低、爆震倾向下降点火时刻相比MBT延迟了2°CA避免早燃。台架验证结果显示BSFC确实下降了4.2%和仿真预测的下降3.8%只差0.4个百分点NOx排放在法规限制内缸压峰值偏差也只有4%CA50偏差在1.5°CA以内。这个量级的吻合度已经达到了用仿真代替盲目试错的工程目标。5.2 燃烧仿真的边界与未来延伸在项目里摸爬滚打几年后我对燃烧仿真的有更清醒的认知。它并不是万能的在极端工况如高EGR超稀释燃烧、零排放氨燃料、氢内燃机下现有的经验模型往往需要重新标定不确定性仍然很大。但这不意味着仿真无用反而意味着要把仿真模型当成一个持续演进的工程工具去维护。我习惯在每个项目结束后把这次标定好的模型参数、机理文件、边界文件归档入库等下次遇到相似工况时直接调用做到“一次标定、多次复用”。这种积累会让你的燃烧仿真能力越用越顺手每一次项目都不是从零开始。做燃烧仿真这些年我最大的心得体会是软件学起来快真正的门槛在“对燃烧机理和流动本质的理解”上。模具设计、台架试验、标定数据、仿真模型这几样东西只有揉在一起看才能做出真正经得起验证的优化方案。技术文章到这里也接近尾声了最后再分享一个小建议——仿真优化迭代前一定要留出至少30%的时间预算给模型标定和问题排查没有标定的仿真就是一把没有准星的枪打不准目标的。希望大家都能在仿真和试验的相互验证中找到自己最得心应手的那条路。