《Extreme Materials》：Volatility sets high-entropy carbide service life above 3000
K, and both wall temperature and propellant chemistry set the required hafnium fraction
编辑推荐：
块体相稳定性并不决定高熵碳化物在火箭应用中的服役寿命，表面氧化尺度的氧化与退移才是关键。研究人员将正规溶液稳定性模型、平衡热化学与退移动力学结合的筛选方法应用于(Hf,Zr,Ti,Ta,Nb)C体系。预测无组分式自旋分解（compositional spino
块体相稳定性并不决定高熵碳化物在火箭应用中的服役寿命，表面氧化尺度的氧化与退移才是关键。研究人员将正规溶液稳定性模型、平衡热化学与退移动力学结合的筛选方法应用于(Hf,Zr,Ti,Ta,Nb)C体系。预测无组分式自旋分解（compositional spinodal），故选择落于氧化物尺度；喉部所有候选氧化物均熔融，因此挥发性而非熔点控制退移。损失介于解离汽化（下限）与蒸汽氢氧化物通道（上限）之间，后者对水蒸气为二级反应，使氢发动机喉部退移提升两个数量级。壁温对退移呈指数控制，55次飞行寿命要求壁温低于约1870 K。固定壁温下主导杠杆为水蒸气分压；氢与煤油发动机的蒸汽载荷互换使退移差22倍，而其105 K气体温差仅差1.4倍。铪富集在解离通道上价值8.8倍，而在氢氧化物通道上仅1.06倍。以采用的100 μm/h判据计，七候选中五者通过前缘、无者通过氢发动机喉部；富铪锆碳化物在两项均居首。
研究背景与立项依据
超高温陶瓷（ultra high temperature ceramics, UHTC）包括早期过渡金属硼化物、碳化物与氮化物，具有已知化合物最高熔点，是极端热环境首要候选。最严苛条件超过这些陶瓷所能形成一切氧化物的熔点，此时难熔性不再区分候选材料，熔融尺度的挥发性取而代之。两类代表应用界定本研究：高超声速前缘暴露于小曲率半径解离空气，驻点热流与前缘半径反比；火箭喷管喉部处于声速面高温高压燃气，对流热流全程最大。两者均需材料在极端热机载荷下保持结构与尺寸稳定。高熵碳化物（high entropy carbides, HEC）将UHTC成分空间大幅拓展，Sarker等引入熵形成能力描述符，证实十二种五金属碳化物中九种近2473 K烧结为单一岩盐相，阴离子亚晶格碳空位额外贡献构型熵。后续研究报道硬度超混合法则预测。既有设计多聚焦块体单相稳定性，但该性质不决定服役寿命；陶瓷可3000 K热力学稳定仍因表面降解失效，氧化尺度演化决定抗氧化性。ZrC在1570至1770 K数分钟即实测氧化，而过渡金属二硼化物形成保护性硼硅酸盐玻璃，约1370 K以上保护劣化，动力学由抛物转准线性再转线性。碳化物无硅不能成硼硅酸盐玻璃。(HfZrTiTaNb)C及其二硼化物直接研究已确认氧化尺度化学而非块体相稳定性主宰性能与寿命。另一局限是计算研究多始于密度泛函理论（density functional theory, DFT）精修块体稳定性，昂贵仿真前置。本研究反转顺序，沿Olson过程-结构-性质-性能框架，以动力学作为热力学稳定性实际边界；相可热力学稳定但因氧化尺度挥发过快而不适用。本研究发展低计算强度筛选：闭式稳定性判据缩搜空间，平衡热化学定氧化物尺度及挥发物种，双区退移模型转挥发性为退移率。体系为(Hf,Zr,Ti,Ta,Nb)C及低组元子集。HfO?与ZrO?抗氧化优于Ta?O?与Nb?O?已定性确立，本研究将其转为氢发动机喉部定量退移率，量化铪含量效应，并由JANAF热化学函数推Nb氧化物分压。方法可迁至其他氧化受限陶瓷。
主要技术方法
研究人员采用三阶段筛选。其一，组态熵稳定性屏幕：各碳化物作金属亚晶格混合单一岩盐相，混合摩尔Gibbs能写正规溶液G=0.5x?Lx+RTΣ?x?lnx?，十对金属交互参数由DFT（Quantum ESPRESSO 7.3.1, PBE, ONCV SG15）算等原子二元混合焓L??=4ΔH???(0.5)，0 K下混合Gibbs能等于混合焓，由Gibbs自由能Hessian本征值判单相稳定，最小本征值零温即等原子自旋分解温。其二，氧化物尺度与挥发压：弃用Reaktoro+NASA CEA（缺Hf/Zr/Ta/Nb凝聚相与挥发物种），改用第三定律热化学，依各氧化物实际汽化反应（非一致升华）由表列Gibbs自由能函数与298 K标准反应焓算平衡分压；HfO?经HfO(g)+O(g)，ZrO?经ZrO(g)+ZrO?(g)，Ta?O?不一致汽化为TaO?(g)+O，Nb在喉部还原态作NbO?、前缘氧化态作Nb?O?；水蒸气引入MO?(c)+nH?O(g)=MO?(OH)?(g)二级氢氧化通道。尺度成分取块体比极限与稳态富集两极限。其三，退移模型分双区：前缘自由分子Knudsen数未达阈值作Hertz-Knudsen-Langmuir极限；喉部边界层传质控制，J?=h?,?p?/(RT)，h?由Bartz湍流关联+Chilton-Colburn类比Sh=0.026Re?·?Sc1?3，喉部直径作特征长，二元扩散Chapman-Enskog+Neufeld积分；RS-25（氧氢分级燃烧富蒸汽还原喉）与F-1（氧煤富碳喉）两发动机参数取自公开规格，壁温由气侧热流hg(Taw?Tw)传导平衡定。蒸发系数α取0.01至0.1工作估计并与α=1上限对照。
研究结果
3.1 稳定性屏幕：正规溶液模型下等原子五元碳化物等原子自旋分解温2240 K，低于2473 K合成温，与实验单相一致；全局共切线法混溶隙2253 K仅高13 K。最低能失稳为Ti-Zr分相（正混合焓17.3 kJ/mol），Ti-Hf次之12.0 kJ/mol；偏离等原子约2.8 at.%向Ti+Zr富集即令自旋分解温超合成温。去Ti后HfZrTaNb四元自旋分解温降至381 K，三元HfZrNb为55 K，二元HfZr无混溶隙，五元二硼化物<50 K无隙。所有候选在合成与服役期满足单相稳定判据，故排序由氧化定。
3.2 两环境氧化物尺度与挥发性：喉部燃气3248 K超所有候选氧化物（含HfO?/ZrO?）熔点，退移由氧化尺度汽化定而非熔点；熔点只决熔融尺度抗气剪机械移除，在热化学模型外。挥发性序明确：富Ta、Nb氧化物蒸气压最高退移最大；水蒸气经硼硅酸陶瓷转挥发氢氧化物、低氧分压促SiC活性氧化成SiO(g)，故硼硅系喉部差。增Hf、Zr降退移最多。表面保护指数SPI=?log??(v/v_crit)，v_crit=100 μm/h。前缘2200 K干空气下SiO?仍保护，ZrB?+SiC居中，纯二硼化物因B?O?(g)挥发最差；碳化物同序，富HfZr最优。七候选中五者过前缘判据、无者过喉部判据。
3.3 实测氧化物挥发性推绝对退移：RS-25喉3248 K下(Hf,Zr)C尺度挥发压132 Pa、退移265 μm/h；换Ta得933 Pa/1880 μm/h，换Nb得1580 Pa/3190 μm/h，含Ta+Nb四元尺度1820 Pa/3670 μm/h；F-1同序约半值。纯ZrO?尺度438 μm/h，耦合等原子HfO?-ZrO? 265 μm/h，纯HfO? 50 μm/h。前缘2200 K干空气(Hf,Zr)C尺度5.2×10?? Pa、α=1时0.22 μm/h。
3.4 与实测烧蚀一致性核对：大气氧乙炔炬（同热化学机制、低温、异成分）退移数百μm/h量级，HfC+20wt%SiC约100 μm/h、Cf/ZrB?-SiC-ZrC约430 μm/h；固体火箭发动机因50至100 bar对流增强+Al?O?颗粒机械冲蚀高三个数量级（碳碳约1至2×10? μm/h，HfC+SiC约5.7×10? μm/h）。模型265 μm/h居清洁燃气热化学贡献位，作边界而非验证。
3.5 设计图与敏感性：以100 μm/h判据，RS-25喉需Hf组分0.91，F-1喉需0.50（等原子居阈）。壁温指数控制：RS-25与F-1喉气温差仅105 K（3%），仅升温自F-1至RS-25在不变传质下Hf需由0.50升至0.87，余至0.91来自RS-25传质系数高18%（室压/喉径/燃气定）。故煤油机标定成分不等同氢机标定。水蒸气压平方律：氢/煤油气载互换退移差22倍，远大于105 K温差1.4倍。反应焓不确定传递后HfZr退移不确定因子1.3（238至297 μm/h），Nb系2.4，但成分序不变。Hf富集在解离通道8.8倍、氢氧化通道1.06倍。
3.6 非理想尺度混合影响：HfO?-ZrO?近理想固溶，退移对活度模型不敏；五氧化物Ta?O?/NbO?活度系数扫两数量级，低至0.01时Ta/Nb系退移194/207 μm/h仍高于二元265 μm/h以下不成立，且稀释亦稀HfZr基质，地板仍由HfZr基质定，降地板需增Hf而非改活度模型。
3.7 极端温复现实验设计规则：3873 K下W无凝聚氧化物，还原炬氧势<5×10?? bar时W最不挥发，火箭喉氧势高数级则W成高挥发气态氧化物。四金属挥发压序Hf≈Zr?Ta
讨论部分总结
研究人员指出服役环境定成分可行性而非仅偏好：前缘靠干富氧空气硼硅酸盐玻璃，喉部水蒸气+低氧分压令其失效；喉部优材靠HfO?-ZrO?难熔混合层，该材前缘亦佳但不享硼硅自愈。各合金元素角色分化：Ta、Nb供HEC构型熵却成最挥发氧化物大幅增退移；Hf反向。抗氧性随去Ta/Nb/Ti/B/Si、增HfO?-ZrO?单调升，但同元素供韧性/烧结性/熵稳定，故抗氧与最大熵需不同成分——解为结构而非成分：高熵碳化物基体供力学与工艺，表覆富Hf涂层成HfO?-ZrO?层；或沿壁富燃注射降热流兼降局部水蒸气压（退移∝水汽压平方）。方法学上DFT仅算十对二元混合焓，氧化压取实测热化学，退移取动理论+传质关联，贵仿真留予不可测者。局限：Nb喉部压由JANAF第三定律推，与前人外推差6%；前缘Nb?O?自1988 K外推至2200 K；尺度成分由稳态通量平衡解而非假定块体，去主自由参；单斜HfO?标准熵文献59.33/56.15/57.7 J·mol?1·K?1三值，取最低增HfO?压约20%、Hf/Zr挥发比自9降至7，Hf需量变化≤0.01；未表低氧化态通道令预测偏保守低估；熔尺度剪切排液μm厚时与热化学率可比，故预测为总退移下限。构型熵作理想处理，Hew递归熵法显五元<1%修正、2000 K上可略。两发动机为连续设计面两点，Fig.4横轴跨液机喉温，质量转移系数差异使等值线微移，温度为主效。三区：<3060 K全HfC-ZrC连结满足100 μm/h；3060至3340 K Hf需量随温陡升，优化强发动机特定；>3340 K纯HfO?亦超限，须主动冷却/涂层/梯度架构。RS-25喉距上限于92 K，氢机类碳化物表面近成分单独可达极限，氧化物挥发性比块体相稳定性施更严约束。
结论部分翻译
单相稳定性不决定高熵碳化物火箭服役寿命，寿命由氧化层化学与挥发性控。本热力学-动力学筛选得三主论：一、所有候选碳化物合成期保单相，氧化分析选出的富HfZr成分全服役温程稳定，选材由氧化而非块体相稳定性控；二、优选成分依环境变，富HfZr碳化物两环境均最优，而环境定是否有成分可行——七候选前五过前缘判据、无过喉部判据，喉部退移由壁温与水蒸气压而非底层碳化物成分定，故梯度架构+富燃壁膜联用；三、分析由采用服役判据而非本征材料限给出定量热化学退移率与对应成分需量，DFT限于二元交互焓余取实测评估热化学，富Hf碳化物在表征燃气中直接退移测量仍待实验验证。因主导挥发通道由壁面气体定，推进剂化学成材料设计变量，本框架供探此依赖之基。
（论文发表于《Extreme Materials》）
