“他们用的这种算子分裂法在处理热传导和烧蚀质量损失这两个模块的耦合时，分裂误差在时间步长上应该是一阶的，而不是二阶。

    要做到他们声称的二阶精度，那就需要在每一步分裂之后额外加一个校正步骤，但是这个步骤他们没有写，公式里也没有体现。

    要么是他们漏写了，要么是他们实际跑的时候根本没有用这个分裂法。”

    韩柏岩凑过来仔细看了看肖宿指出来的内容。

    算子分裂法他知道，但在数值方法的收敛精度上他还没那么敏锐，毕竟这也不是他的专业方向。

    “可之前我们在超算上跑过他们的模型，”韩柏岩说，“按照他们论文里公开的公式重新搭了一遍，跑出来的结果确实是收敛的，最后算出来的壁面温度跟他们公布的数据相差不大。如果他们的方法有问题，那为什么我们复现出来的结果是对的呢？”

    肖宿把论文翻到数值仿真那部分的参数表那块，扫了一眼，然后抬起头说道：“因为他们在另一个地方做了手脚。”

    他的手指点在参数表里一行标着“等离子体边界层化学反应的活化能修正因子”的数据上。

    这个修正因子论文里给的值是一点三五，旁边引了一篇参考文献说这个值是根据量子化学计算得出的。

    但是肖宿顺着那篇参考文献的标题回忆了一下，那是一篇做分子动力学仿真的论文，里面计算的活化能体系跟艾瑞昂论文里用的材料体系根本不是同一种东西，化学反应路径都不同，活化能的数值在物理上根本不可能通用。

    “他们引的参考文献是做碳化硅基材料表面催化的，活化能计算结果是针对硅氧键断裂的。

    而他们自己论文里用的材料是一种酚醛树脂基烧蚀材料，烧蚀过程中的主要化学反应是碳氢键断裂和碳氧键重组。

    硅氧键的断裂活化能和碳氢键差了将近一个数量级，把前者的修正因子拿来用在后者身上，属于张冠李戴了。”