NASA CEA化学平衡计算程序实战：从火箭燃烧室到比冲估算

简介：NASA开发的CEA（Chemical Equilibrium with Applications）程序是一款面向火箭发动机研究者、燃烧工程专家及能源化工从业者的热化学计算工具，专用于求解高温高压环境下化学反应达到平衡时的组分浓度、温度、压力及热力学性质，能够为发动机燃烧室设计、燃料选取与排放评估提供定量依据。该zip压缩包共含7个文件，整体大小约3.06MB，其中有3个MATLAB脚本、2个PPT教程、1个程序源码文件和1个MAT格式热力学数据库文件，既覆盖了CEA的基础原理与操作步骤，也提供了可直接运行/调用的计算代码和数据，便于用户从零上手。目前已有804人浏览或学习过这份资料。借助包内的CEA for MATLAB快速上手课件和示例脚本，用户可以逐步掌握输入初始温度、压力、混合比等参数的方法，理解输出结果中各组分浓度、比冲、特征速度等关键指标的含义，并据此完成燃烧产物预测、发动机整机性能评估或燃烧不稳定性分析，为实际工程和科研工作提供可靠参考。 搞火箭发动机、喷气推进、爆轰计算的人，手头几乎都留着一个叫CEA的程序。CEA全称 Chemical Equilibrium with Applications，是NASA Glenn研究中心（早年间叫Lewis研究中心）开发的化学平衡计算程序。做燃烧室设计、算理论比冲、预测燃烧产物，它几乎是行业默认参考工具。哪怕你只是做热力学数据库校核、燃料评估，或者写课程作业，这套程序也能省掉大量手算平衡组分的痛苦。

这篇文章不是NASA编程原理源码精读，而是一份“怎么把CEA用起来”的实战笔记。我会从它到底能算什么、怎么拿到手、怎么输入一个火箭算例，到输出文件怎么读、哪些坑我踩过，一次说清楚。适合刚接触CEA的学生，也适合从网页版换到本地版、想批量跑工况的工程师。

1. CEA到底是何方神圣：一个化学计算器如何撑起航天发动机设计

1.1 从名字拆解：Chemical Equilibrium with Applications

CEA的核心是“化学平衡”。在高温燃烧室里，几十种组分同时存在，它们之间的反应速率极快，经过极短时间后达到热力学平衡态。CEA做的事情，就是给定反应物的种类、配比、温度和压力，通过最小化体系吉布斯自由能，算出平衡时各组分的摩尔分数。

听起来像是一个物理化学教科书问题，但真实现实远比教科书复杂。真实燃烧涉及几十种可能出现的中间产物，比如氢氧燃烧体系里不只H2O、H2、O2，还有H、O、OH、HO2、H2O2等；碳氢燃料体系里还要加上CO、CO2、C（石墨）以及各种碳氢自由基。手解平衡方程基本不可能，所以必须靠程序做数值求解。

NASA开发CEA的初衷，就是给火箭发动机和航空发动机设计人员提供一个通用、可靠、更新及时的热化学计算平台。它最早可以追溯到上世纪60年代的Gordon和McBride那套化学平衡程序，后来的CEA在数值方法、热力学数据库、用户接口上都做了大量升级。目前公开版本被全世界科研机构、高校和商业公司用作参考基准，很多商业软件（比如某些计算流体力学软件里的燃烧模型）在给化学机理做简化时，也会拿CEA的平衡结果做对标。

1.2 它能算什么：从燃烧室温度到比冲的一条龙计算

CEA不是一个“只能算一个温度点”的小工具。输入不同问题类型，它能直接给出工程上常用的一整套性能参数：

最常用的还是rocket模式。在这个模式下，你告诉程序燃烧室压力、混合比、喷管面积比，它不仅能算燃烧室内的平衡组分和温度，还能分别用“平衡流”和“冻结流”两种假设推算燃气在喷管中膨胀后的状态，最后给出理论比冲和特征速度。比冲是火箭发动机最重要的性能指标，CEA能在没做实验之前就给一个可信的理论上限。

另外一个很有用的功能是“分批计算”。你可以在一个输入文件里写多个case，每个case用不同的压力或混合比，这样一条命令就能扫描出性能随参数变化的曲线。我一开始不知道这个功能，傻乎乎地一个参数一个参数手动改文件，效率极低。

2. 获取和运行CEA：别被Fortran源码吓住

2.1 官方渠道与在线Web版

CEA的源代码和数据文件是公开的，NASA Glenn Research Center官网提供了下载入口，里面包含Fortran源程序、热力学数据库、传输物性数据库，还有一份非常详细的用户手册。对于有编译基础的人，用gfortran或Intel Fortran都能编过，不过要注意旧代码在某些现代编译器上会有一些兼容调整，需要按照README里的提示设置编译选项。

如果你不想折腾编译，NASA还维护了一个Web版CEA界面。在网页上直接选择问题类型、填压力、混合比、燃料氧化剂名称，点运行就能拿到输出文件。网页版对新手极其友好，验证一个小算例、查一组数据都很快。唯一的问题是Web版不太适合批量扫描参数，也不方便嵌入到你自己的脚本里。如果你只是偶尔用一下，网页版足够；如果你要做参数优化、批量算几十上百个工况，还是本地版更顺手。

2.2 本地编译运行的两种常见姿势

本地运行有两种主流姿势。第一种是下载官方Fortran源码自己编译。编译后通常会得到一个可执行文件（比如fce.exe），运行时按提示输入输入文件名，程序会生成一个输出文件。你也可以把输入文件重定向给程序，写成批处理脚本循环跑多个文件。这种方式的优点是“原汁原味”，完全可控；缺点是需要自己处理编译环境和文件IO，刚上手时有点门槛。

第二种姿势是用第三方封装，比如Python包里就有人把CEA的核心数据库和计算逻辑封装成了rocketcea。你不需要会Fortran，只要pip install rocketcea，然后写几行Python就能算一个火箭发动机的性能。本质上它调用的是CEA的底层数据和求解算法，但用起来更现代，也能方便地画图、做优化。这个包不是NASA官方发布的，但在我测试的绝大多数案例里，结果跟官方源码一致，作为快速验证工具很实用。

运行官方程序的交互逻辑很简单：启动后程序提示输入文件名，你把之前编辑好的.inp文件名字敲进去，回车，几秒钟后同一个目录下会出现一个.out文件，所有结果都在里面。如果你重定向输入，比如执行fce.exe < h2o2.inp，程序会把输出打到屏幕上或写到默认输出，具体看版本设定。我个人习惯每个算例一个独立.inp文件，文件名用h2o2_20bar.inp这种带含义的命名，避免批量跑的时候把输出文件覆盖了。

2.3 输入输出文件的生成逻辑

CEA的输入文件是纯文本，看起来有点像Fortran的namelist，但有自己的规则。一个输入文件可以包含多个case，每个case都以problem开头，用end结束。problem后面可以带case=xxxx给这个算例起个名字，输出文件里会用这个名字区分不同case。

这种“块状结构”非常方便批量计算。我经常用一个Python脚本生成一整个.inp文件，比如压力从10bar到100bar，每5bar一个case，自动生成20个块，然后一次运行，最后再写一个解析脚本把.out里的温度、比冲、分子量全部抽出来画曲线。这种做法让CEA变成一个完全自动化的计算引擎，而不是一个只能手点一下的“计算器”。

输出文件则是规整的文本表格。里面有输入条件回显、平衡组分摩尔分数、燃烧室热力学参数、不同面积比下的喷管出口参数等。信息量很大，第一次看可能有点懵，但只要知道几个关键字段的位置，后续解析就很顺畅。

3. 实战演练：用CEA算氢氧火箭发动机燃烧室参数

3.1 一个最小可用的输入文件

用最经典的氢氧组合来演示。假设我们要计算一个燃烧室压力20bar、混合比8（氧化剂质量/燃料质量）的算例，输入文件如下：

problem case=h2o2_20bar rocket p,bar=20 o/f=8 fuel H2 298.15 oxid O2 298.15 end

解释一下每一行。rocket表明这是火箭发动机问题；p,bar=20设定燃烧室压力，单位是bar；o/f=8是混合比；fuel H2 298.15和oxid O2 298.15分别指定燃料和氧化剂，后面的数字是初始温度，单位K。这里用气态氢和气态氧，所以不需要指定相态。如果你要用液氢液氧，需要在物质名后面加(L)，比如fuel H2(L) 20，括号里可以给温度，否则程序按该物质在给定温度下的饱和态处理。

把这个文件保存成h2o2_20bar.inp，然后用编译好的CEA可执行文件运行，选择输入这个文件，等待输出文件生成。

3.2 输出结果逐行解读

输出文件的前半部分会把你的输入条件回显一遍，接着就是计算结果的正文。核心区域大概长这样：

不同版本、不同热力学数据库计算结果会有微小差异，但量级基本如此。注意比冲通常以“秒”为单位，在CEA输出里有时也给“m/s”单位，需要换算时除以标准重力加速度9.80665。

燃烧室温度是热设计的关键。3350K已经很接近很多材料熔点以上了，所以实际发动机必须考虑再生冷却或烧蚀冷却。如果燃料混合比偏离化学当量（氢氧当量比约8），燃烧温度会下降，比冲也会变化。你可以在输入文件里改o/f=10再跑一下，马上就能看到温度降几百度，比冲也可能降低。

3.3 把结果用于估算喷管出口速度与比冲

CEA的rocket模式还能计算喷管流动。要在输入文件里加喷管参数，最常用的是面积比ae/at。比如：

problem case=h2o2_20bar_ae50 rocket p,bar=20 o/f=8 ae/at=50 fuel H2 298.15 oxid O2 298.15 end

ae/at=50表示喷管出口面积与喉部面积之比为50。程序会从燃烧室平衡态出发，分别按平衡流和冻结流两种假设，等熵膨胀到出口，给出出口温度、压力、马赫数，以及比冲。这时候你会看到平衡流假设下的比冲比冻结流高一点，因为燃气在膨胀过程中继续释放化学能，给了额外的推力。

在实际工程估算里，真实比冲往往介于两者之间。CEA的好处是你一次就能拿到上下界，设计初期的性能包络就有了。我一般都会同时跑一个平衡流和一个冻结流工况，然后告诉结构工程师用哪个范围去定贮箱增压需求。

4. 输入卡片语法和热力学数据库：细节决定成败

4.1 基本语法结构

CEA输入语法其实很宽容，大小写不敏感，逗号可以加空格，但有几个规则必须遵守：

• 每一个case以problem开头，以end结束。
• problem后面可以跟case=名字，建议每次都加，方便区分输出。
• 一个输入文件里可以连续放多个case，程序会依次计算。
• 除了计算类型和参数设置外，不会的语句默认忽略吗？不，程序会报错。所以不要写多余的注释行，除非你知道注释语法。老版本CEA里行首的感叹号!一般被当作注释，但我不建议依赖它，命名的case已经足够说明问题。

4.2 燃料、氧化剂和惰性组分怎么填

物质名称必须与CEA热力学数据库中的名称完全一致。常见物质如下：

• 单质或简单化合物：H2、O2、N2、CO、CO2、H2O、CH4、C2H5OH等。
• 相态：气体一般不加后缀，液体在名称后加(L)，如H2(L)、O2(L)，固体加(S)，如C(gr)代表石墨。
• 如果燃料是混合物，比如某种配方燃料，可以用reagent把额外组分加进去，并用质量分数或摩尔分数指定比例。

我做过一个用煤油和液氧的算例，燃料不是单一化合物，就需要在fuel行后面用reagent添加组成。语法大概是：

fuel RP-1 reagent C 0.86 H 0.14

这里的C和H指的是碳和氢的元素质量占比。具体格式可以参考用户手册里的“reagent”章节。这种自定义反应物能力很强大，但也是最容易出错的地方——我之前把质量分数和摩尔分数搞混，算出来的温度差了800多K。

4.3 温度和压力背后的物理陷阱

温度参数在CEA里有两个位置：一个是反应物初始温度，一个是计算得到的燃烧温度，前者是输入，后者是输出。不要试图让CEA“直接设定燃烧温度”来算平衡组分，程序会通过能量方程自己求绝热火焰温度，除非你指定t模式做等温平衡计算。

压力单位是经典坑。CEA支持p,bar=、p,psi=、p,atm=，有些版本还支持MPa。不同case里可以混用，但同一个case内部一定不要混合使用。我见过有人前面写p,bar=20，后面在另一个参数里写p,atm=20，导致结果完全偏离。建议统一用bar，因为工程里常用，数值也容易检查。

混合比o/f的定义是氧化剂总质量除以燃料总质量。对氢氧来说，化学当量比是8左右；但如果燃料不是H2，这个数字没有直观意义，要换算成当量比phi。CEA还支持用equivalence ratio直接输入，某些版本里是phi=1.0，这样无论什么燃料，都能直接反映富燃还是贫燃。对于不熟悉混合物化学的人，用当量比会安全很多。

5. 我踩过的坑和排查思路

5.1 单位混用导致的“不可能”结果

有一次我算一个甲烷-氧气算例，输入了p,psi=300，本意是300psi，但后面又写了一个ae/at=20，没注意之前的压力单位。程序跑出来燃烧室温度只有2400K，怎么看都不对。后来我把压力换算成bar一算，发现300psi大约是20.68bar，跟我想要的100bar差距巨大。问题就出在我上一个case用的p,bar=100，复制下来改单位时只改了数字，忘了改单位后缀。

从那以后，我在每个输入文件名里都会标注单位，比如ch4_o2_100bar.inp，跑完先看输出文件里回显的压力是否正确，再往下看温度。这个习惯帮我避免了很多无效计算。

5.2 相态缺失导致物性偏差

CEA热力学数据库对气体、液体、固体的热力学性质分别存储。如果你要算液体燃料，却没有加(L)后缀，程序会默认按气态算。气态甲烷和液态甲烷的总焓差别很大，燃烧温度自然不同。

我刚开始用CEA时算液氢液氧火箭，直接把燃料写成fuel H2，氧化剂写成oxid O2，算出来的温度和比冲跟资料对比偏低了。后来发现是没有处理相态，加上(L)并指定温度后，结果立刻合理了。这里顺便说一下，指定液态物质的温度时，最好用正常沸点附近的值，比如液氧在90K，液氢在20K。CEA会根据这些温度调用对应的焓值。

5.3 反应物名称拼写与数据库匹配问题

CEA不是搜索引擎，它对别名没有容忍度。写methyl alcohol会报错，正确名称是CH3OH；写ethanol也可能不认，要写C2H5OH。如果你不确定名称，去用户手册后面的“Thermochemical Data”表格里慢慢找。

有一次我想算某种高能燃料，手册上没有现成名称，只能用元素组成和生成焓自定义。CEA允许你通过h,cal/mol=或h,kcal/mol=指定物质的生成焓，然后用reagent添加。这时候务必搞清楚单位，我因为漏掉了kcal/mol和cal/mol之间的1000倍换算，结果燃烧温度高出天际，检查了很久才发现。

5.4 收敛失败处理技巧

当燃烧室压力极高（比如几百bar）或者混合物组分接近临界状态时，CEA可能告诉你“无法收敛”或者输出奇怪的负数。这时候别急着怀疑程序。

我的排查顺序是：先降低压力测试，看是不是压力超出数据库适用范围；再检查混合比是不是写反了；最后看反应物名称和相态是否正确。都排除之后，可以尝试给problem行后面加一个n=1或调整迭代参数，让求解器更稳健。老版本CEA还有一些求解控制变量，但不同版本命名不同，建议直接看用户手册。

最笨但有效的方法是把复杂混合物简化：先用单一燃料算一个基准case，然后逐步加入reagent。这样一旦不收敛，能很快定位到是新加入的哪种组分引起的。

6. 从火箭到气候：CEA在排放研究中的扩展应用

6.1 燃烧产物的完整物种清单

CEA虽是为航天发动机设计的，但在航空发动机和气涡轮机的排放估算里也经常被引用。因为它的平衡计算不会漏掉高温下那些重要的小分子和自由基，比如NO、OH、H、O等。燃烧室温度越高，热力型NOx生成越明显，CEA算出的NO摩尔分数可以作为定性参考。

我曾用CEA扫描航空煤油替代燃料（正十二烷）在不同当量比下的CO2和H2O产量，得到一个很漂亮的理论曲线。虽然真实发动机燃烧室不是完全平衡状态，但理论极值能告诉我们这个燃料体系在最理想情况下能生成多少CO2，这对排放清单的上限估算很有价值。

6.2 把CEA结果接入排放清单和气候模型

NASA这些年做的气候数据研究中，很重要的一块是航空排放对大气的影响。要估算航班排放总量，除了需要燃油流量和航班轨迹，还需要“单位燃油燃烧产生的各物种质量”。CEA给的就是这个系数。

你可以用CEA算航空煤油完全燃烧的理论产物，然后结合真实发动机的燃烧效率进行修正，得到CO2、H2O、NOx的排放因子。把这些排放因子乘上全球航班耗油数据库，就能估算出整个航线网络向大气中排放的水汽和碳氧化物量，进而输入气候模型评估高空卷云等影响。CEA在这里不是主角，但它提供了一个快速、一致的热化学计算基础，比我们自己去翻庞大的化学机理数据库方便得多。

我自己在实际项目里的做法是：用CEA生成一系列当量比下的产物组分表，然后用一个Python小工具把这些表转化为插值函数，再嵌入到排放计算模块。这样整个链路至少从“理论燃烧”这一步就和NASA官方程序保持一致，汇报的时候也更有说服力。可以说，CEA不只是一个“老古董”火箭工具，它在现代碳中和、排放评估里依然有不可替代的复用价值。

最后分享一个小经验：别一上来就求新求全，先从最简单的氢氧火箭case跑通，输出文件能看懂、结果量级合理之后，再去碰液体燃料、自定义组分、批量扫描这些东西。CEA的文档写得非常详实，但几百页手册没人愿意从头啃。你只需要一个最小算例作为锚点，后面所有扩展都能在它的基础上一点点加。这种“先跑通，再复杂”的思路，在做任何工程计算工具的时候都适用。

本文还有配套的精品资源，点击获取