综述：基于DWSIM实现液态金属反应器中甲烷干重整详细动力学模型的方法学模拟框架

《International Journal of Hydrogen Energy》：A methodological simulation framework for implementing a detailed kinetic model using DWSIM for dry reforming of methane in a liquid metal reactor

　　

甲烷干重整（Dry Reforming of Methane, DRM）是实现合成气（syngas）生产的有前景途径，同时能够资源化利用CO₂和CH₄等温室气体。然而，其高度吸热特性要求高效的热传递，且催化剂因积炭导致的失活仍是主要挑战。液态金属作为一类新兴的反应介质，具有优异的热导率和稳定性。本研究提出了一种将复杂多步朗缪尔-欣谢尔伍德（Langmuir-Hinshelwood, L-H）动力学表达式集成到开源过程模拟器DWSIM中的方法学框架。该框架通过一个理论案例研究进行展示，即适用于假设性液态金属催化环境的DRM与逆水煤气变换（Reverse Water-Gas Shift, RWGS）反应的耦合。液态金属特有的物理现象被纳入并进行分析，以确保与模型简化在机理上的一致性。推导的速率方程使用文献中的动力学参数作为基线，通过自定义Python脚本实现。研究结果展示了该框架在指导实验设计方面的潜力，所有定量预测均受研究中量化的结构性和参数性不确定性约束。然而，计算结果受制于缺乏可用于验证的实验数据，这成为评估反应本身化学动力学的关键制约因素。

1 引言

DRM是合成气生产的最有前景途径之一，能够同时资源化利用CH₄和CO₂两种主要温室气体。该反应为高度吸热反应，通常需要1073–1173 K以上的温度才能实现足够的转化率和工艺性能。因此，高效传热流体（Heat Transfer Fluid, HTF）的选择是DRM反应器设计和运行中的关键问题。此外，快速催化失活因积炭形成而成为工业部署的重大障碍。

传统HTF如空气、蒸汽、有机导热油或熔盐在此语境下存在显著局限：空气和其他气体热导率低，需要高流速才能实现有效传热；有机导热油和有机流体在673 K以上发生化学分解；熔盐虽已成功应用于聚光太阳能热发电系统，但其操作温度范围很少超过1073 K，且面临高黏度和腐蚀问题。这些不足突显了对能够满足DRM极端温度要求的先进HTF的需求。

液态金属作为新一代高性能HTF在核应用和聚光太阳能热发电技术中已崭露头角，近期更应用于高温化学反应器。其可分为碱金属（如Na）、重金属（如Pb-Bi或铅铋共晶LBE）以及低熔点金属（如Sn）。其吸引力在于以下优越性能组合：宽广的操作温度范围（多数候选金属可达1273 K以上），高的热化学稳定性，极高的热导率，高的允许热通量，以及在某些情况下的热能储存能力。与熔盐相比，液态金属黏度显著更低，且具有很高的沸点。

然而，液态金属的应用也伴随着特定局限。最相关的挑战包括高温腐蚀问题、相对较高的成本、安全问题（尤其是碱金属如液态钠遇水或空气时的高反应性），以及工业应用中最大的挑战之一——积炭形成导致的催化剂快速失活。

本研究的目标明确为方法学性质而非预测性质。鉴于目前液态金属反应器中DRM实验数据的缺失，本研究追求五个目标：（i）建立将详细表面机理转化为DWSIM可实施速率表达式的严格可复现框架；（ii）将液态金属特有物理现象纳入动力学框架并提供理论依据；（iii）量化参数从固相催化剂域向液态金属域转移时模拟输出的不确定性；（iv）识别将本框架转化为经验证预测工具所最关键需要的实验测量。模拟结果以不确定性约束的理论估计形式呈现，定量的转化率和选择性预测应被理解为框架能力的说明而非设计规格。

2 DRM研究的趋势与挑战

2.1 热力学与反应路径

DRM为高度吸热反应，其主要缺点在于实现可观转化率所需的高温。主反应为：

CH₄(g) + CO₂(g) ? 2CO(g) + 2H₂(g) ΔH°_298K = +247.3 kJ/mol

通常需要1073 K以上的温度；然而该过程不仅受表面动力学限制，还受热力学约束和副反应影响，主要包括逆水煤气变换反应（RWGS）、甲烷裂解、布多阿尔反应（Boudouard反应）以及气化反应。在大多数DRM实验中，产物中可检测到H₂O。某些液态金属（如Ga、Sn）对氧和氢具有不同的亲和力，可能使RWGS更为显著。

2.2 催化系统与动力学建模

在液态金属催化反应器中，金属同时作为溶剂和催化剂。其关键特征是活性位点位于液态金属表面，反应通过吸附物种进行。液态金属表面能有效解离分子（尤其是CH₄），碳在金属中的溶解度可能显著，形成潜在的碳储库。

本研究采用的动力学模型基于吸附-反应-脱附机理，涉及CH₄、CO₂及中间物种在活性催化位点上的反应。该模型借鉴了多项参考文献，包括Boudart和Djéga-Mariadassou（1984）以及Rostrup-Nielsen（2002）针对Rh/γ-Al₂O₃催化剂开发的DRM一维多相动力学模型，以及Snock等人提出的积炭沉积和气化反应模型。

3 方法学

3.1 模型假设与热力学分析

建立可靠的表面反应机理本质复杂；然而文献中已提出多种DRM反应的基元机理。在构建动力学模型时，必须分析与积炭形成相关的反应路径。本研究中，由于使用液态金属时积炭倾向显著降低（因碳可扩散至体相、再分散或更易被氧化），积炭形成被忽略。

对于传热分析，液态金属被假设为热导体。系统在高温（1073–1273 K）和常压下运行。相关活性位点位于液态金属表面，该表面被视为能够吸附和活化CH₄和CO₂的催化相。该催化表面由等效位点阵列表示。

在推导表面反应速率定律时做出以下假设：液态金属表面包含等效位点阵列；吸附和反应在这些位点上均匀发生。一维模型在径向温度和浓度梯度可忽略时通常足够。本研究中径向传热阻力和径向传质限制均被忽略，假设稳态、层流、常压流动。这些假设显著简化了系统，其有效性在4.3节中进行了解析评估。

3.2 朗缪尔-欣谢尔伍德（L-H）动力学模型开发

本研究开发了一个涵盖第2.1节中所有全局反应的详细L-H动力学模型，包括23步基元反应，涉及表面吸附物种和气相物种。机理顺序遵循：物种的吸附和脱附、表面扩散、吸附中间体的化学反应。

基元步骤包括：CH₄的解离吸附（DRM的速率控制步骤，Rate-Determining Step, RDS）；CO₂的分子吸附（拟平衡步骤）；吸附氧与含碳自由基的反应；吸附氢与其他相邻位点自由基的反应；以及RWGS反应的并行发生。此外还包括导致积炭的气化反应。

液态金属表面吸附CH₄和CO₂，但最大覆盖率低于固体催化剂。表面扩散快速，表面物种可部分再溶解到液态金属中，有助于缓解碳积聚。在相对较低温度范围（873–1043 K），CH₄吸附被视为RDS；而在更高温度（1043–1223 K），吸附CH₄的热裂解成为RDS。本研究采用后者作为动力学适应模型。RWGS反应的RDS为CO₂活化。

完整的23步机理通过表格详细呈现，包括吸附、表面反应和脱附步骤。为动力学建模，采用满足整体热力学平衡约束的L-H形式速率表达式，RDS在所研究反应条件下被视为实际上不可逆，但包含热力学驱动力以确保在平衡附近的行为正确。

3.3 液态金属特有的L-H框架修正

虽然第3.2节的 kinetic 模型遵循为固体多相催化剂开发的经典L-H形式，捕捉了DRM的基本表面化学，但它并未固有地考虑液态金属反应介质的独特物理性质。本节将液态金属特有现象纳入L-H框架。

3.3.1 表面更新动力学与有效位点空位

在经典固体催化剂中，活性表面位点空间固定。在液态金属中，对流流动和热波动持续更新气-液界面。该表面更新由更新速率s（s^-1）表征。无量纲表面达姆科勒数（Damkohler number, Da_s）比较更新速率与吸附物自然脱附速率：Da_s = s/k_des。对于液态金属上的弱吸附物种，脱附速率常数k_des在1273 K反应温度下约为10⁸–10¹² s^-1，因此Da_s << 1，表明脱附远快于表面更新。界面组成主要由吸附热力学而非流体动力学周转控制。

3.3.2 碳向金属体相的溶解

固体催化剂与液态金属操作之间最显著的机理区别是：吸附碳具有第三种去除途径——溶解到液态金属体相中。这与气化（主要是步骤10）和表面积累（通过甲烷分解）竞争。

碳溶解速率遵循对表面覆盖率θ_C的依赖关系，并受溶解热力学驱动力控制：r_diss = k_dissθ_C(1?C_bulk/C_sat)，其中k_diss为溶解速率常数，(1?C_bulk/C_sat)捕获趋近饱和的程度。在1273 K操作温度下，碳在液态金属中的溶解度非常有限，例如Sn的饱和浓度C_sat ≈ 20–30 mol/m³。碳形成、气化和溶解之间的平衡使稳态表面碳覆盖率θ_C,SS保持在10^?3量级或更低，为在L-H速率表达式的分母中忽略θ_C提供了定量机理依据。

碳在体相液态金属中累积的净速率决定了再生前反应器的操作寿命：r_C,net^bulk = k_crackp_CH4θ_v ? k_dissθ_C,SS(1?C_bulk/C_sat)。达到饱和所需时间为t_sat = ρ_LMV_LMC_sat/(r_C,net^bulkM_C)。

3.3.3 气-液金属界面的传质

在催化剂表面为液态金属的气泡柱反应器中，整体反应速率可能受气态反应物向气-液界面传输的限制。适当的理论框架是双膜模型，其中物种i从体相气体到液体表面的通量与其体相分压和界面值之差成正比。当气相阻力可忽略时，传质速率由液相传质系数k_L（m/s）和比界面积a（m²/m³）表征。

传质与表面动力学的相对重要性由Hatta数（Ha）量化：Ha = √(k_rxnD_LM)/k_L，其中k_rxn为伪一级速率常数，D_LM为限制性反应物（CH₄）在液态金属中的扩散系数。当Ha << 1时，系统为动力学控制，总有效因子η_GL → 1；当Ha >> 1时，反应主要在液膜中进行，速率受传质限制。平面薄膜不可逆一级反应的有效因子为η_GL = tanh(Ha)/Ha。

由于缺乏k_L和D_LM的实验数据，无法进行Ha的可靠估算。基线模拟中假设动力学控制极限η_GL = 1。

3.3.4 液态金属表面的吸附平衡

液态金属上反应物的较弱吸附（低覆盖率行为的基础）目前缺乏直接实验验证。因此，对从Ni/Al₂O₃ kinetic 模型（Benguerba等人，2017）的参数进行了系统转移。所有物种的吸附焓使用线性缩放因子S = 0.60 ± 0.10进行缩放，该因子源自支持液态金属催化剂的密度泛函理论（Density Functional Theory, DFT）研究。此缩放反映了与固体过渡金属的局域d轨道相比，吸附物与液态金属离域电子态之间的较弱相互作用。

3.4 L-H kinetic 模型的推导

使用完整模型需要计算资源和非常详细的表面覆盖率信息，这超出了DWSIM的便捷处理能力。因此，在推导反应速率定律时进行了简化，对最丰富的物种和处于平衡的物种做出合理假设。

基于准平衡近似（除RDS步骤4和11外的所有步骤均处于平衡），速率表达式推导如下：对于DRM，以CH₄解离为RDS，在弱吸附假设下（∑K_ip_i << 1），速率表达式为r_DRM = k₄K₁p_CH4(1?K_DRM^?1p_CO²p_H2²/(p_CO2p_CH4))/(K₁p_CH4+K₂p_CO2+K_?21p_CO+K₃p_H2+K_?23p_H2O+1)²。

对于RWGS反应，以CO₂与吸附氢的活化（步骤11）为RDS，速率表达式为r_RWGS = k₁₁K₂K₃p_CO2p_H2^1/2(1?K_RWGS^?1p_COp_H2O/(p_CO2p_H2^1/2))/(相同分母)²。

碳形成通过甲烷裂解方程（式(3)）计算，采用幂律动力学：r_crack = k_crack·p_CH4。布多阿尔反应（式(4)）为放热反应，在高温下平衡强烈向左移动，有利于CO而非碳。在操作温度下，布多阿尔反应在热力学上非自发，对本研究中的碳形成贡献可忽略。

3.5 速率控制步骤（RDS）指定的验证

3.5.1 速率控制度（Degree of Rate Control, DRC）分析

Campbell的速率控制度X_RC,i定义为X_RC,i = (k_i/r_overall)·(?r_overall/?k_i)，其中偏导数在保持所有平衡常数K_i和其他速率常数k_j≠i恒定时取。对于DRM路径，求导得到X_RC,4 = 1，确认了速率定律的数学结构与RDS假设的一致性。同样，对于RWGS路径，X_RC,11 = 1。

3.5.2 准平衡假设的表面覆盖率验证

对于非RDS步骤的准平衡假设要求反应中间体的覆盖率足够小。使用BEP缩放后的吸附常数在1273 K和典型操作分压下，CH_x中间体的覆盖率比值为：θ_CH3/θ_v ~ 10^?1，θ_CH2/θ_v ~ 10^?2，θ_CH/θ_v ~ 10^?3。若将三者纳入分母，总和增加约10^?1，对Ω的增加不到10%。对于RWGS路径，关键中间体HCOO*的覆盖率θ_HCOO/θ_v ≈ 10^?6，确认了步骤11远离其自身平衡。

3.5.3 RWGS速率控制步骤的机理验证

步骤11而非直接CO₂活化（步骤8）被确定为RWGS的RDS，基于液态金属的电子结构。DFT机理研究表明，RWGS反应在过渡金属上的优选路径由表面的氧亲和力控制。在强氧结合金属（如Ni、Rh、Pt）上，直接CO₂解离有利；在弱氧结合金属（如Cu、Ag、Pd）上，CO₂解离在能量上不利，通过甲酸盐（HCOO*）的氢辅助路径成为主导路径。液态Sn具有极低的氧亲和力，不形成稳定体相氧化物，因此直接CO₂解离（步骤8）在Sn上动力学不利，而CO₂氢化为甲酸盐（步骤11）构成CO₂活化的最低能垒路径。

3.6 参数可转移性：域转移评估

在本研究中使用的 kinetic 参数最初针对Ni/Al₂Al₃固定床催化剂确定，现已适应于液态金属表面。这构成了域转移：在固体催化剂域（局域活性位点、强吸附、多孔载体）中校准的参数被应用于根本不同的域（液态金属表面、可移动位点、弱吸附、体相溶解）。通过三个组分评估此转移引入的系统偏差：

机理可转移性：23步机理对大多数步骤转移时仅有轻微修改，因为DRM的表面化学是保守的。最不确定的转移涉及RWGS的RDS（步骤11），该步骤依赖甲酸盐中间体的形成。在固体表面上，HCOO通过相邻固定位点的双齿配位稳定化；在可移动的液态金属表面上，这种明确定义的配位几何可能无法实现，可能削弱甲酸盐结合并使机理向直接CO₂活化路径（步骤8）转移。

速率常数可转移性：BEP原理将表面反应的活化能与反应焓联系起来。当所有金属-吸附物键能均匀从Ni缩放到Sn时，解离活化步骤显著变得更吸热。然而，此估计忽略了液态金属特有的若干现象：氢可能部分溶解到Sn体相中，提供超出纯表面缩放关系的额外稳定化；液态界面的动态非晶性质可能改变过渡态稳定化。鉴于这些影响的定量影响尚未解决，本研究采用70–90 kJ/mol的较高活化能作为工作假设。

吸附常数可转移性：液态金属上的平衡吸附常数通过线性因子S = 0.60 ± 0.10缩放Ni参考值的吸附焓获得，该因子来自液态金属模型表面的DFT研究。前置因子保持与Ni相同，隐含类似吸附熵的假设。

基于域转移，模拟输出预测的可靠性可分为：最可靠的是H₂/CO < 1的定性预测；中等可靠的是温度趋势方向；最不可靠的是绝对转化率值。

4 模拟

DWSIM是提供稳态质量和能量平衡的开源化学工程模拟器。本研究使用DWSIM 9.0.2版本进行 kinetic 模拟，采用用户自定义的活塞流反应器（Plug-Flow Reactor, PFR），通过嵌入式Python脚本在每个迭代步骤计算转化率。

4.1 DWSIM反应分析

假设的设计条件和参数包括：反应器温度固定为1273 K（等温）；流量为0.02 mol/s；CH₄/CO₂流量比为1。为简化模型，假设CH₄和CO₂为纯气体，不含杂质；RWGS反应是重整过程中唯一发生的副反应；由于催化剂即为液态金属，将整个反应器体积视为催化剂。

初始尝试在DWSIM PFR中直接包含碳形成失败，因为该反应器模型设计用于单相（气体）反应，无法容纳固体碳作为产物相。这一技术限制 necessitated 解耦方法。碳形成通过独立的 kinetic 模型评估，发现其比主DRM速率低两个数量级，因此未纳入反应器模拟。

4.2 结果与讨论

以下结果为名义参数值下动力学框架生成的理论估计，不应被解释为经验证的预测。

碳形成在所研究条件下可忽略（<0.05%的CH₄消耗），不影响合成气组成。模型参数源自Benguerba等人（2017）和Abanades与Martinez（2023）的实验数据，使用BEP原理进行参数适配。

在假设η_GL = 1（动力学控制）和等温操作于1215 K的条件下，模型给出甲烷转化率在70–85%范围，二氧化碳转化率为85–100%（ consistently 高于甲烷转化率，因其同时通过DRM和RWGS路径消耗）。CO选择性估计在所有条件下高于80%，H₂/CO比率在所有模拟条件下均低于1，表明RWGS副反应的显著贡献。

模拟理论结果为液态金属DRM反应器设计提供了若干关键洞察：最佳温度为1150–1250 K，催化剂负载量为200–300 kg/m³，CO₂/CH₄比为1.0–1.2时可最大化转化率同时最小化碳风险。

4.3 模型假设的严格评估

对简化系统的五项主要假设进行了严格评估，每项假设单独放宽以量化对预测转化率的影响。

4.3.1 等温反应器操作

基线反应器因液态金属的高热导率和气泡柱中的强烈混合而保持均匀温度。数学松弛通过耦合能量平衡与质量平衡实现：dT/dz = [p_catA_c[(?ΔH_DRM)r_DRM + (?ΔH_RWGS)r_RWGS] + UaA_c(T_w?T)] / ∑_jF_jC_p,j。

对于绝热操作（Ua = 0），温度在反应器前几厘米内下降超过200 K，使反应淬灭。对于Ua = 5000 W/(m³·K)，温度降约为24 K，甲烷转化率比等温基线低4.7个百分点。该假设引入的预测偏差较小，予以保留。

4.3.2 无传质限制（η_GL = 1）

基线反应为动力学控制；气相反应物以体相分压到达液态金属表面。通过引入气-液有效因子η_GL进行数学松弛。Hatta数和η_GL的界限估计为：k_rxn为10^?1–10¹ s^?1，D_LM为10^?9–10^?8 m²/s，k_L为10^?4–10^?3 m/s，导致Ha为0.03–3.2，η_GL为0.28–1.0。

对于η_GL = 0.75，转化率降低14.7个百分点；对于η_GL = 0.28，转化率降低40.9个百分点。这是模型中最大的单一假设不确定性。

4.3.3 恒定物性和无压降

在基线模型中，气相压力在反应器全段为1 bar，忽略压降。考虑1 m高气泡柱的静压，底部入口压力约为1.64 bar，线性降至顶部出口的1 bar。更高的压力增加反应物分压从而提高反应速率，但也增加产物分压使平衡向反应物移动。使用温度相关关联式（NASA多项式、Chung模型、修正Eucken公式）和耦合非等温能量平衡，性质变化导致转化率变化小于0.5个百分点，该简化完全合理。

4.3.4 低表面覆盖率（Ω ≈ 1）— 弱吸附极限

基线假设L-H分母中所有吸附项之和可忽略（∑K_ip_i << 1）。使用BEP缩放后的吸附常数在1273 K和典型进出口条件下计算，入口处分母Ω ≈ 1.24，使DRM速率降低至Ω^?2 ≈ 0.65；出口处Ω ≈ 1.06。耦合模拟给出的甲烷估计转化率为60.2%，比Ω ≈ 1基线低8.5个百分点。

4.3.5 可忽略的积炭形成

基线模型假设碳表面覆盖率θ ≈ 0，通过快速碳溶解论证。为验证，纳入表面碳平衡自洽求解：dθ_C/dt = r_crack ? r_diss ? r_gas = 0。对于新鲜液态金属至约50%饱和容量，碳覆盖率增加Ω不到2%，速率降低几个百分点（< 5%的模拟不确定性）。

核心价值发现：对于新鲜液态金属，碳覆盖率的 kinetic 反馈使DRM速率降低约2%；即使在90%饱和时，该影响仍可管理（<20%速率损失）。

综合耦合误差（表7）：对于新鲜金属，总速率损失为1.9%；50%饱和时为3.7%；90%饱和时为16.4%。这证明了在整个有用操作窗口内（直至再生成为必要）基线模型中解耦碳处理的合理性。

4.4 关键参数的敏感性分析

4.4.1 解析敏感性系数

DRM总速率对参数p_i的归一化敏感性系数定义为S_pi = (p_i/r_DRM)·(?r_DRM/?p_i)。解析结果为：S_k4 = 1.0；S_K1 = 1 ? 2K₁p_CH4/Ω；S_Ki = ?2K_ip_i/Ω（i = CO, CO₂, H₂O）；S_K3 = ?K₃p_H2^1/2/Ω。

在1273 K和典型中反应器条件下（Ω ≈ 1.146），各参数敏感性为：k₄为+1.00；K₁为+0.86；K₂为?0.077；K₃为?0.018；K_?21为?0.0015；K_?23为?0.0023。

4.4.2 有效因子η_GL的影响

η_GL是模型中最大的单个不确定性来源。由于DRM和RWGS总速率均与η_GL成正比：r_DRM^eff = η_GLr_DRM^int，r_RWGS^eff = η_GLr_RWGS^int。

甲烷转化率在转化率低于约50%时与η_GL近似线性。关于H₂/CO比率，若两个反应的η_GL相同，该比率保持不变；若传质阻力因扩散系数或溶解度不等而不同，则H₂/CO会相应变化。H₂/CO < 1的定性预测对于η_GL在合理范围内的任何值都是稳健的。

4.4.3 全局敏感性分析：框架与未来实施

提出的解析敏感性系数为在名义参数处的局部导数。为超越局部线性并在完整操作空间上量化参数交互，建议将全局方差敏感性分析（Sobol方法）和筛选分析（Morris方法）作为模型开发路径的自然下一步。

4.5 液态金属反应器的挑战与未来工作建议

甲烷和二氧化碳在液态金属反应器中反应的发展程度处于必须从实验室概念验证向相关规模的首个工业示范演进的阶段。面临的重要技术不确定性包括：气体停留时间与局部温度梯度及合成气形成 progress 的相关性；现有 kinetic 数据仅从传统催化系统（如Ni或Rh基催化剂）获得，不对应于液态金属反应介质的条件；模拟工具缺乏 robust 实验验证；液态金属对不锈钢的腐蚀性。

为提升DRM技术成熟度，未来工作应聚焦于：在更宽条件下的实验验证；多组分液态金属合金的扩展；与经济分析模块的集成。

按预期对输出不确定性影响排序的实验测量优先级为：第一优先级（关键）为确定液态Sn上DRM的活化能E_a,DRM和前指因子A_DRM；第二优先级（关键）为测量液相传质系数k_L和CH₄在液态Sn中的扩散系数D_LM；第三优先级（高）为确定RWGS反应的机理；第四优先级（高）为直接测量CH₄和CO₂在液态Sn表面的吸附平衡；第五优先级（中）为确定液态Sn在操作温度下的碳溶解度极限C_sat。

5 结论

本研究成功建立并展示了将多步L-H动力学集成到开源过程模拟中的方法学框架，明确定位为理论基础和探索工具而非经验证的预测模型。该框架基于23步表面机理，严格识别RDS，并通过竞争性吸附纳入，可通过自定义Python脚本在DWSIM中有效实现。

液态金属特有修正被引入，模型从固体催化剂域适应到液态金属表面。消融分析量化了每项建模假设的影响，证实了基线模拟代表理论上限。碳形成的解耦处理通过三级耦合分析检验，对于整个有用操作窗口直至接近碳饱和，预测转化率的总误差小于5%。

所有定量转化率和选择性结果均为名义参数值下假设动力学控制（η_GL = 1）的说明性估计。单一最具影响力的结构性不确定性来源是气-液传质有效因子η_GL；H₂/CO比率对参数变化稳健，始终预测低于1。当前案例研究的主要局限在于缺乏液态金属特有实验 kinetic 数据。优先级排序的实验计划识别了将框架转化为经验证设计工具所需的关键测量。