自己写的求解器出现：浮点数例外 (核心已转储)的问题

[小泽同学]

请教一下各位大佬，就是我们课题组在做open foam仿真飞秒激光加工的工作，目前已经有好多毕业的师兄建立了模型和编写求解器，我是在此基础上进一步完善和编写，最近就是在修改变步长的时候遇到了算例算不动的问题，之前是仿真的周期和实际加工的周期不完全一致，现在想通过修改变时间步长的方式将两者统一，但是在我一通修改下，变时间步长的功能实现了，但是算例算一个周期后就出现了下图的错误，
然后我试图寻找原因，发现停止运行的时间步的库朗数130多，然后我减小了十倍的步长，还是报错，但是库朗数变成了0.003，所以我不知道原因出现在那里了，然后我看报错位置是计算求解双温方程的晶格温度时出错的，但我不知道怎么修改或者是什么原因，求助一下大佬们，非常感谢大佬们，第一次发帖，多多指教！
如果有类似研究方向的大佬也请多多交流，谢谢各位！

[李东岳]

自己写的求解器别人很难debug问题在哪。看看有没有类似的研究方向的大佬吧。祝你好运！

[小泽同学]

@李东岳 感谢李老师的回复，也确实是这样，只能自己先摸索摸索了，目前看起来应该是变时间步长的问题，还想问一下李老师，就是一个仿真周期只需要前面一小段时间步长小一点，后续大部分时间不用这个小步长，因为会算的很慢，这种问题是只有变步长这一个解决措施吗？感谢李老师的回复。

[李东岳]

@小泽同学 在 自己写的求解器出现：浮点数例外 (核心已转储)的问题 中说：

就是一个仿真周期只需要前面一小段时间步长小一点，后续大部分时间不用这个小步长，因为会算的很慢，这种问题是只有变步长这一个解决措施吗？

是的老铁

[小泽同学]

@李东岳 好的，李老师，那我自己再修改修改定义步长的代码，感谢李老师，有什么参考的资料吗qaq