解决问题的灵感，来自一个意想不到的方向：
彩色电视机的信号处理算法。
在航天防护项目期间，温卿参与过彩色电视解码电路的优化。
当时的难题是如何从复合视频信号中，干净地分离出色度信号和亮度信号——
两者频带重叠，分离不好会导致色彩串扰和图像模糊。
工程团队采用了一种“自适应滤波”算法：
不是用固定的滤波器，而是根据信号特性动态调整滤波参数。
在信号平稳时用锐利滤波器提高清晰度，在信号突变时用平滑滤波器抑制噪声。
温卿突然想到：
冲击波不也是信号的“突变”吗？
在冲击波附近，物理量剧烈变化，需要平滑处理保证稳定；
在平滑区域，需要锐利处理保持精度。
为什么不把自适应滤波的思想，用到数值算法中？
她立刻开始推导。
传统的godunov格式是“一视同仁”的——整个计算区域用同一种通量计算方法。
她要设计的是“自适应godunov格式”：
在识别出冲击波区域时，自动切换到更耗散但更稳定的算法；
在平滑区域，切换到更精确但更敏感的算法。
关键是如何“识别”冲击波。
温卿想到了一个简单的判据：
相邻格子间物理量的梯度。
梯度大的地方，可能是冲击波；
梯度小的地方，是平滑区域。
但梯度阈值怎么定？
设得太低，会把正常变化误判为冲击波，导致过度平滑；
设得太高，会漏掉弱冲击波。
温卿从记忆碎片中找到一个概念：
“特征值分析”。在双曲型方程中，信息沿着特征线传播。
冲击波对应的是特征线汇聚的区域。
通过分析方程的特征值，可以更精确地识别冲击波的位置和强度。
这个方法的计算量不大，但需要深入理解方程的数学结构。
温卿花了一周时间，重新推导了内爆控制方程组的特征系统。
当最后一套特征值和特征向量出现在草稿纸上时，她长舒一口气。
有了理论准备，温卿开始撰写改进方案。
她刻意采用了非常技术性的语言，聚焦在数值算法层面，避免触及更敏感的物理模型问题。
方案的标题很朴实：
《关于内爆仿真中冲击波捕捉算法的改进建议》。
开篇先肯定现有算法的优点：
“godunov格式在保证计算稳定性方面表现优异，为我国的核武器设计提供了可靠的数值工具。”
然后委婉指出不足：
“但在处理非对称内爆等复杂情况时，数值耗散可能掩盖真实的物理不对称性，导致设计余量估计偏乐观。”
接着提出改进思路：
“借鉴信号处理中的自适应滤波思想，建议开发一种‘特征值引导的自适应godunov格式’。基本思想是……”
她详细描述了算法框架：
第一步，在每个时间步，计算每个格子边界处的特征值。
第二步，根据特征值判断该位置是否可能发生冲击波。判断准则基于特征线的汇聚程度。
第三步，对可能发生冲击波的边界，采用低阶耗散大的通量计算方法；
对平滑区域的边界，采用高阶耗散小的计算方法。
第四步，在冲击波区域和平滑区域的过渡带，采用加权平均，避免算法切换引起的数值振荡。
为了证明可行性，温卿用简化的一维模型做了小规模验算。
结果令人鼓舞：
她把方案和验算结果整理成一份十页的技术报告，按照规定程序，提交给理论组算法小组。
算法小组的负责人是赵研究员，一位五十多岁的计算数学专家，以严谨甚至有些古板着称。
收到温卿的报告，赵研究员的第一反应是皱眉：
“自适应？这个概念在流体力学里有人提过，但在我们这种极端条件下的应用……”
他组织了小组讨论会。
除了温卿，还有三位资深的算法工程师。
会议开始，温卿先做简要汇报。
她特意放慢语速，在黑板上一步步推导，确保每个人都跟上思路。
“这里的关键是特征值判据。”
她指着公式。
“我们不是简单地看梯度大小，而是看特征线的汇聚。这是更本质的物理判断。”
一位工程师提出问题：
“特征值计算本身有计算成本。每个时间步、每个格子边界都要算，总开销会增加多少？
温卿调出预计算数据：
“根据一维模型测试，特征值计算的开销，约占单时间步总计算量的3。
考虑到我们通过自适应减少了平滑区域的高阶计算，整体开销增加控制在5以内。”
“稳定性呢？”
另一位工程师问。
“算法在不同区域用不同格式，切换时会不会引起数值振荡？”
“这是过渡带加权平均要解决的问题。”
温卿展示过渡函数的数学形式。
“我们设计了一个光滑的过渡函数，保证格式变化连续可导。在一维测试中，没有观察到数值振荡。”
赵研究员一直没说话，仔细看着报告。
等所有人都问完了，他才开口：
“温卿同志，你的思路很新颖。但我有个根本问题：
我们怎么知道新算法真的‘更准’？数值仿真没有绝对的真值，只能和实验比。但我们的实验数据……”
他没说完，但意思很清楚：
核试验数据极其有限，而且很多是综合效应，很难单独验证冲击波捕捉的精度。
温卿早有准备：
“赵老师，我们可以用‘方法验证’的思路。虽然不能直接和实验比，但我们可以比较不同算法在‘基准问题’上的表现。”
她调出几个经典基准问题：
sod激波管问题、双马赫反射问题、noh问题……
这些都是计算流体力学中公认的测试案例，有解析解或高精度数值解。
“我用新算法算了这几个问题，和现有算法对比。”
温卿展示对比图。
“在sod问题中，新算法捕捉接触间断的精度提高40；在双马赫反射中，波系结构更清晰；在noh问题中……”
一幅幅对比图显示，新算法在几乎所有基准问题上都表现更好。
赵研究员的眼神开始变化。
从最初的怀疑，到认真，再到感兴趣。
“这些基准问题……和我们的实际内爆还有差距。”
他说，但语气已经松动。