我正在使用一个代码来解决二维非结构化网格上的扩散/反应方程。由于某些过程的僵硬，我从接近 1e-13 的时间步长开始，并以毫秒范围内的最终时间结束，时间步长接近 1e-8。

问题出现在这个问题上，我想称之为几何僵硬。在最初的 ~250 ns 模拟中，我需要一个非常精细的网格来解析强制函数的几何特征。然而，在这段时间之后，系统的强迫大约为零，我不再需要空间分辨率。

精细网格需要大约 2 周的时间来运行前 250 ns。粗网格需要 2-4 周才能将整个模拟运行到 2.5 毫秒（加速一个数量级）。如果不认真重写代码（目前也超出了我的技能范围），就不可能运行自适应网格，但是有一些实用程序可供我重新启动模拟。我正在考虑做的是将细网格运行到不再​​需要为止，获取重新启动转储文件，将其插入到粗网格上，然后重新启动以使其运行到我们想要的最后时间。

根据我的插值方法，我可以很容易地估计插值中的误差，但是我不确定如何估计在我再次开始时间步进时传播的误差。我特别担心温度等变量的平滑，因为我的许多加工速率都以指数方式依赖它。我的问题是，如何在不运行多个案例的情况下估计一些固定错误如何通过我的僵硬系统传播（目前在计算上不可行）。