超临界燃烧仿真实战：从物性模型到收敛判断

最近在跟进一个编号090的超临界燃烧仿真项目前前后后折腾了快两个月总算把从物性预处理到收敛判断的一整套流程跑通了。这个方向的燃烧仿真跟常规的喷雾燃烧完全是两个玩法很多在常压常温下好用的模型一放到超临界工况下就露馅。这篇博文就把这段时间积攒的经验捋一遍——超临界燃烧仿真到底在求解什么物理、为什么常规模型会失效、物性处理和求解策略怎么搭才靠谱。准备做高压燃烧室、超燃冲压发动机或者火箭推力室相关仿真验证的朋友这篇能帮你省掉大量试错时间。1. 超临界燃烧仿真到底在解决什么问题1.1 什么是超临界态为什么燃烧偏偏要往这个区域推先说基本概念。纯物质有一个临界点温度和压力分别超过这个点之后物质就进入超临界态。以甲烷为例临界温度约190.6K临界压力约4.6MPa煤油的临界压力通常在2.0~2.4MPa临界温度在640~720K左右。工程上更常遇到的是跨临界状态——喷射进入燃烧室时流体压力已经高于临界压力但温度还可能在临界温度以下随后在加热过程中穿越临界温度线进入超临界区。很多人第一次接触超临界射流时最大的困惑是它到底是气体还是液体答案很微妙超临界流体既不完全是气态也不完全是液态。超过临界压力之后液相和气相之间的分界面彻底消失加热液相燃料不会像常压下那样先沸腾产生气泡而是密度连续下降变成一种稠密流体。形象一点说常压下的沸腾就像教室里有人喊了一声下课大家陆陆续续往外走有明显的波次而超临界态下的连续变密度更像晚高峰地铁站散场所有人同时往出口挪动没有抵抗力弱的气液分界面可依赖。燃烧室为什么要在这种状态下工作核心原因是循环效率和推力性能。透平进口温度、燃烧室压力越高布雷顿循环的热效率越可观。航空发动机的燃烧室压力已经从早期的0.5MPa级别抬升到目前的4~6MPa新型预冷组合动力和超燃冲压方案的燃烧室压力往往更高。燃料以超临界或近临界状态喷射有助于强化掺混、缩短火焰长度让燃烧更稳定更紧凑。代价是数值模拟的难度陡增因为经典的喷雾燃烧框架在这里几乎失效。1.2 常规燃烧仿真框架在超临界工况下的失灵点传统的气液两相燃烧仿真核心是一套离散相液滴连续气相的描述方式。液滴有表面张力有蒸发率公式有液滴直径分布用拉格朗日方法追踪颗粒轨道质量和能量源项在欧拉网格中做耦合。这套框架在航空发动机常压或者中等压力工况下用了很多年成熟可靠。可一旦燃烧室压力超过燃料的临界压力这套框架的两个根基就悬空了。第一表面张力和相变机制消失了。临界态下表面张力趋近于零液滴不再是一个有明确界面的圆球而是变成一团密度的局部集中区。没有界面就没有经典的蒸发模型可挂靠液滴尺寸分布这种输入参数也变得没有物理意义。第二密度变化极端剧烈。从喷射出口的稠密流体到混合后的高温燃气密度可以相差一到两个数量级而且是连续变化不存在一个清晰的相界面把液和气分开。流体压缩性、非理想气体行为、真实流体的比热变化全面接管了物理过程理想气体状态方程在超临界工况下做出来的结果几乎就是错的。所以超临界燃烧仿真这条技术路线本质上是把喷射、混合、点火、燃烧整个过程当作单相连续介质中的密度层化流动来对待。不用喷雾模型不做离散液滴而是基于真实流体物性直接求解跨临界/超临界射流在高温高压环境中的扩散、混合和化学反应。这看起来是回避了两相流的麻烦实际上是把麻烦转移到了物性模型和湍流-化学耦合这两个方向上。2. 物性模型是超临界燃烧仿真的命根子2.1 状态方程与热物性计算的选型逻辑在超临界条件下理想气体状态方程的误差大到不能接受。原因很简单高压下分子间距缩小分子间的引力和分子体积效应占据了不可忽略的比例真实气体行为偏离理想模型非常明显。工程上最常用的替代方案是立方型状态方程——Peng-RobinsonPR和Soave-Redlich-KwongSRK。选哪个更多取决于流体类型和个人习惯。碳氢燃料体系做超临界燃烧仿真PR方程用得更多它在临界点附近的体积-压力行为还原度比SRK好一些尤其对煤油这类多组分混合物PR配合合适的混合规则能比较合理地复现临界温度、临界密度附近的物性变化趋势。SRK在烃类相平衡计算中也有大量应用但临界区附近的密度预测偏差通常比PR略大。光有状态方程还不够。燃烧仿真需要密度、比热、焓、熵、声速、输运系数黏度、热导率、扩散系数随压力和温度连续变化的物性场。物性模型必须自洽也就是说由状态方程算出的p-V-T关系要和由热力学关系式积分得到的焓、熵、比热相互一致否则能量方程里的源项和状态方程的耦合会出现系统性的虚假振荡。工程上建物性数据库的常规做法有三种。第一种是全部调用商用软件或者开源库的实时物性计算接口比如Refprop、CoolProp精度高但计算开销极大大规模网格的CFD计算几乎跑不动。第二种是用查表法事先把物性数据离散到压力-温度网格上仿真时插值获取速度极快但对网格分辨率敏感临界点附近如果离散太粗会出现插值跳变导致的数值振荡。第三种是折中方案用扩展对应态原理ECS或者基于对应态原理的简化关联式把关键物性拟合成连续函数比如密度用PR方程解析计算比热和焓用工况表格插值输运系数用对应态关联式估算。我在项目里用的是第二种加第三种混合——状态方程解析求解输运系数查高精度表格这样速度能接受精度也保得住。一个不得不强调的细节是临界点附近物性的剧烈变化这是超临界燃烧仿真中大量数值问题的根源。压力在临界压力附近时流体的等温压缩率、比热、热扩散系数都会出现尖峰或者极值。有人做过统计二氧化碳在临界点附近的热扩散系数会跌到接近零而比热可能出现数倍甚至更高的峰值。这就像开车过连续减速带——物理上本来就是急剧变化数值格式如果不够稳直接就发散给你看。2.2 气液界面处理的物理视角转变传统喷雾燃烧中非常关心的雾化质量Sauter平均直径破碎模型标定在超临界燃烧仿真中基本不再适用。这不是说雾化不重要而是说雾化的形态完全改变了。低于临界压力喷射时液体射流经历完整的破碎序列完整液核→袋状破碎→丝状破碎→液滴分布。这个过程可以用喷嘴几何、韦伯数、Ohnesorge数等来关联经验模型。但一旦压力超过临界压力喷射出口处的工质只要高于临界温度射流就不会再有经典的破碎过程而是以类气态射流的方式向环境中扩散——没有液滴没有表面波破碎混合过程主要由射流动量和高密度比带来的大尺度涡结构控制。跨临界的情况更有意思。射流初始温度低于临界温度、但环境压力高于临界压力时射流核心区可能存在一个稠密液核但这个液核的边界不是稳定的相界面而是随着流动不断从两相区过渡到单相区传统意义上的液核长度也很难定义。数值上处理这种状态工程上更常用的办法是放弃两相界面追踪直接在整个计算域内用统一的真实流体物性求解N-S方程让密度场自己去刻画哪里稠密、哪里稀疏。这样的处理方式干净利落省掉了界面追踪、相变模型、表面张力源项这一大堆复杂环节代价是要求CFD求解器有非常稳的密度基数值格式能够在密度从800kg/m3一路降到5kg/m3的计算域内不出现负密度、不产生压力锯齿。用我同事的话说我们不模拟雾化我们模拟的是密度场如何自我演化成云雾状结构。这句话是对超临界燃烧仿真物理图像最浓缩的概括。3. 湍流与化学反应的耦合光有湍流模型不够3.1 大涡模拟在超临界燃烧仿真中的主导地位超临界射流从喷嘴喷出后与环境气流之间的密度差极大近场剪切层会形成大尺度的拟序结构这些结构直接决定燃料和氧化剂的掺混效率进而控制火焰的稳定位置和释放形式。雷诺平均RANS对这种强间歇性、强各向异性的混合层预测能力有限时间平均抹掉了很多对燃烧起关键作用的瞬态混合事件。大涡模拟LES成了超临界燃烧仿真事实上的主力工具。LES直接解析大尺度涡结构仅对小尺度湍流做模化对射流近场的卷吸、夹带、大尺度振荡的还原度远高于RANS。代价是网格量暴涨近场剪切层要求网格尺度进入泰勒微尺度量级再加上高压工况密度变化剧烈时间步长被声速限制压得很小算一个中等复杂度算例往往需要几千核时起步。做LES时网格设计有一个常见误区——只关注壁面附近的y却忽略了射流近场的网格各向异性。超临界射流的密度梯度集中在剪切层内如果这个区域的网格长宽比过大伪数值扩散会把混合锋面抹平进而让点火位置判断偏移。我的建议是射流出口附近的网格尽量保持各向同性加密范围覆盖到射流核心区两侧至少各3倍射流直径的范围剪切层内的网格尺度控制在喷嘴直径的1/20以下。3.2 湍流-化学相互作用模型的选择边界湍流-化学相互作用TCI是燃烧仿真的核心难点超临界工况下这个问题被进一步放大。涡耗散模型家族EDM、EDC在常压湍流扩散火焰中使用广泛核心假设是化学反应速率受湍流混合速率控制化学动力学时间尺度远小于湍流混合时间。但超临界燃烧室中高压环境让化学反应速率加快同时流体物性变化剧烈湍流混合与化学反应的速率比值不再有清晰的量级差EDM类模型的假设前提被打破直接套用的结果往往是火焰位置过于靠前或者燃烧效率被严重高估。工程上更稳妥的选择集中在两类模型。一类是火焰面模型把湍流火焰看成一系列层流小火焰面的统计集合通过混合分数和标量耗散率来建表。稳态火焰面Steady Flamelet在超临界燃烧中应用较广计算量适中对扩散火焰形态较好的算例精度不错。但火焰面模型隐含假设化学反应时间尺度远大于湍流小尺度时间超临界工况下这个条件也在边界上用得时候需要评估。另一类是输运PDF模型用拉格朗日蒙特卡洛粒子描述组分和温度的概率密度函数演化直接构造化学反应源项不依赖混合时间假设理论适用范围最广。代价是每步都要跟踪大量随机粒子计算成本极其昂贵工程大算例很少敢全程用更多是作为基准结果验证简化模型的正确性。如果你做的是工程验证而不是前沿研究我的务实方案是用火焰面模型做主算但用输运PDF或者精细机理在冷态流场抽几个代表性位置做对标。这样既保住了主算例的计算可行性又能对火焰面模型在这个工况下的适用范围心里有底。火焰稳定位置、上游预热区温度分布这两个指标是对标时最值得盯的它们对TCI模型的选择最敏感。4. 实操设置从网格、边界条件到求解策略4.1 算例设计与物性预处理超临界燃烧仿真不能上来就拎一个复杂的全尺寸燃烧室开算那是自找麻烦也不利于排查问题。稳妥的做法是先设计一个物理过程简化的代表性算例。以我做高压煤油超临界射流燃烧的例子来说计算域选取一个带环形外流的轴对称射流燃烧室燃烧室压力设定为5MPa空气入口温度1100K煤油从中心喷嘴以喷前温度650K注入。这个工况下煤油处于跨临界-超临界过渡区间压力已经超过临界压力煤油临界压力约2.2MPa温度正接近临界温度煤油临界温度约700K很多典型的超临界效应都能体现出来。物性预处理在整个项目流程中工作量占比最高也是最容易被低估的环节。煤油是多组分混合物不能简单套纯物质的物性表。我的做法是先根据已知的组分构成正构烷烃、环烷烃和芳香烃的典型质量分数分布用PR状态方程和经典混合规则van der Waals混合规则加二元交互系数搭建一个替代流体模型把它在0.5MPa到10MPa、300K到2500K的温度区间内离散成物性查表节点。温度靠近临界区时节点加密到5K一个间隔远离临界区放宽到20K这样插值误差可控。这个替代流体模型的焓值、比热必须单独校核。方法很简单把PR方程算出的物性数据和NIST数据库中相近组成的煤油替代模型的公开数据做对比偏差阈值控制在5%以内。燃烧室高压工况下比热不仅随温度变化剧烈随压力也在显著改变如果物性表没有包含压力维度后期结果一定出问题。网格策略上射流近场是绝对的关键区。射流出口直径初定0.5mm从出口往下游15倍直径范围内做逐级加密近场网格尺度控制在喷嘴直径的1/30到1/50壁面边界层第一层网格高度按目标y在1左右设置。超临界流动中壁面热流对火焰驻定有显著影响近壁网格不能太稀。边界条件设置有一个容易遗漏的点——超临界射流进口的湍流参数设定。喷入燃烧室的稠密流体虽然初始密度高但流速也高进口湍流强度不能按常压低速射流那套给。我在算例中湍流强度设为5%湍流长度尺度按喷嘴直径的10%来初始化这样近场剪切层能更快进入充分发展状态避免初始过渡段过长。4.2 数值格式、时间推进与收敛判据超临界燃烧仿真的求解器选择密度基求解器是主流。压力基求解器的压力修正算法在高度可压缩、密度大范围变化的流场中收敛性和稳定性都吃力。密度基求解器配合AUSM类或者Roe类迎风格式对强激波、强密度间断的捕捉能力要好很多。超临界射流近场的密度梯度非常陡钝格式的数值耗散会把火焰结构抹平高阶格式如三阶MUSCL或五阶WENO虽然计算量大但值得在近场区域使用。时间推进方面如果做稳态或者准稳态分析可以用隐式时间推进配合局部时间步长加速收敛如果目标是捕捉瞬态火焰振荡细节就必须走显式或者低耗散隐式格式做LES时间推进时间步长受CFL条件控制同时声速在高压高温燃气区域较高这会进一步压缩时间步长。算例中主燃烧区声速大约在700~900m/s量级网格尺度0.1mm时纯显式推进的时间步长只能给到约1e-7秒量级要跑几十毫秒的物理时间就意味着几十万步的迭代量级这就是超临界燃烧仿真真正昂贵的原因。收敛判据不能只看残差曲线下降几个数量级。超临界燃烧算例中残差可能在某些区域出现周期性波动这是射流近场涡脱落的正常表现不等于发散。我建议同时监控三类信号来判断收敛一是流场关键点温度序列是否进入统计定常二是出口截面的组分流量是否稳定三是燃烧区内的当量比分布是否不再显著漂移。三者都稳定才算真正收敛。火焰初始化是另一个常见的收敛拦路虎。冷态流场直接点燃高压超临界射流往往要在很大计算量下才能等火焰发展起来而且中间可能因为局部过高温升导致数值震荡。我的做法是分两步走先算冷态混合流场确认速度场、密度场合理然后在射流下游某个稳定混合区给一个高温点火斑温度抬升到1800K左右持续数毫秒再用人工增大的反应速率常数把火焰诱导驻定最后逐步降温恢复真实反应速率。这个过程听上去简单但每一步都要反复确认流场没有发散趋势再往下走。5. 常见问题与排查技巧实录5.1 典型问题速查表跑了这个算例前前后后踩了不少坑把典型问题汇总成一张速查表方便大家遇到类似问题时快速定位方向。问题现象可能原因排查方向发散早期负密度或负压力物性表插值越界或密度基格式对高密度比间断捕捉不足检查物性表覆盖区间加密临界区节点近场改用高阶格式温度场出现棋盘状锯齿压力-速度失耦或物性跳变导致的数值振荡检查网格正交性近场加密改用AUSM类格式降低振荡火焰驻定位置偏上游湍流-化学模型过度强化混合评估火焰面模型标量耗散率的截断上限或改用输运PDF验证燃烧全区温度偏低反应机理活化能过高或初始点火斑能量不足校核机理在高压下的点火延迟增大点火斑面积和温度延长点火时间冷态OK点燃后残差飙升点火过程过陡局部温升速率超出格式承受极限采用人工增大反应速率分步诱导法减缓点火释放出口组分流量不稳射流近场大尺度涡结构未统计收敛延长统计平均时间确认监测点序列进入定常分布5.2 几个重要的调参经验第一个经验是冷态流场必须先走到统计定常再点燃火焰。很多人图省事初始化之后直接开反应结果火焰上游区域的速度脉动尚未充分发展燃烧释放的反馈效应和冷态流场耦合在一起就会出现一种假性振荡——看起来燃烧在进行实际上热释放和湍流结构互相打架最后要么火焰熄灭要么出现非物理的周期性脉动。第二个经验是物性曲线的平滑处理远比想象重要。临界点附近比热和密度随温度的变化曲率很大如果物性表里的数据点分布不均插值时会产生一阶导数不连续数值格式就会感受到虚假的物性冲击波。我在实际项目中对临界点附近采用单调三次Hermite插值替代线性插值这种看起来不起眼的细节对抑制压力锯齿有奇效。第三个经验是关于网格无关性验证的方向。超临界燃烧的网格敏感性往往表现为火焰位置对近场网格密度的依赖。如果两套网格比如基础网格和加密1.5倍的网格算出来的火焰位置相差超过一个射流直径问题基本出在近场剪切层分辨率不够而不是化学模型问题。加网格要优先加射流剪切层而不是均匀加密整个计算域。第四个经验是实验锚点的用法。超临界燃烧的实验数据通常只有特定位置的点测或者光学诊断数据但哪怕数据再稀疏也要在后处理阶段把数值结果的对应位置切片提取出来做对比。射流近场的混合分数分布或者火焰中下游的径向温度分布这些是判断整个计算是否物理的锚点。我见过太多算例全局指标很好看跟流场内部结构的实验数据一对中间掺混特性完全对不上。6. 结尾经验补充最后再分享一点个人体会吧。超临界燃烧仿真这个方向最核心的门槛其实不是求解器也不是湍流模型而是物理图像的转变——你能不能接受没有液滴、没有相界面、所有事情都发生在连续密度场里这套世界观。一旦接受了这个设定很多决策就顺理成章了物性表必须精细处理状态方程必须用真实流体模型湍流-化学耦合不能套常压结论网格和格式要向密度梯度区倾斜。实际项目的推进路径建议从冷态超临界射流混合算起不要一上来就点火等摸清了物性和密度场的脾气再加反应这是性价比最高的一条路线。这个主题后续还可以向多组分燃料真实化学反应机理、辐射传热耦合、旋流燃烧室全环模拟等方向扩展但核心的物理认知和数值功底在把090号算例完整跑通的过程中就都会打好底子了。