第90章 刚性方程问题
  坐在陈景明左侧的王海峰副教授，此时满眼血丝，头髮乱得像个鸡窝，显然已经熬了好几个通宵。
  他苦笑著摊开手，声音沙哑：“陈院士，不是我们不努力，这实在是————算不动啊。
  “”
  他指著屏幕上的数据解释道：“大家请看，这是我们团队將李教授的34步化学反应模型，和周教授的高精度边界层模型耦合进求解器后的运行情况。”
  【写到这里我希望读者记一下我们域名101??????.??????】
  “刚开始计算的前几微秒，一切正常。但一旦激波扫过，化学反应开始剧烈发生，刚性问题就立刻爆发了。”
  王海峰切换了一张ppt，上面显示著一个令人绝望的数字。
  “为了保证雅可比矩阵不奇异，我们的时间步长被自適应算法自动压缩到了10^—13
  秒。”
  “这是什么概念？这意味著，我们要模擬飞行器飞行1秒钟的物理过程，即使用上咱们超算中心全部的算力，也需要算上整整100年！”
  会议室里响起了一阵倒吸凉气的声音。
  “能不能尝试用隱式时间积分？比如全隱格式？”陈院士问道，“隱式格式不是可以突破cfl条件限制，放大步长吗？”
  “试过了，陈老。”王海峰无奈地摇摇头，“这就涉及到了第二个死结。一旦放大步长，非线性叠代的矩阵条件数直接飆升到了10^16以上，而这已经是双精度浮点数的极限了，计算机根本分不清这是数据还是误差。结果就是矩阵变得极度病態，线性方程组求解器直接发散，程序崩溃，直接nan了。”
  “那————能不能用算子分裂法？”陈院士团队中的一人说道，“把流动和化学反应分开算？”
  “绝对不行！”