[讨论]计算发散

NASA

<P>计算发散现象常见的, 在设计算法的过程中,需要注意算法的收敛性,稳定性,相容性</P>
<P>一般在很少的时间内结果较好,而在长时间模拟过程时,出现发散的现象,可能是稳定性不好,试着用较小步长看看.如果仍然不行,那需要考虑算法的问题了,它的相容性怎么样?是否是原方程的适当的离散形式?它是否收敛到原模型的解析解等等?</P>

<P>根据别人的论述开个头，看看大家都是什么观点？</P>

FSI

算法的收敛性、稳定性、相容性应该是在计算之前就必须验证的。当然对于复杂的问题来说，很难对算法进行分析。<br>
<br>
步长大小都是相对的，小步长带来的累积误差也不可小看，特别是在长时间模拟的情况下，记得看过一篇文章，对保结构算法和runge-kutta法作了对比
( 显式runge-kutta算法是经常使用的了，其稳定性要求时间步长相对较小)
，结果显示在长时间计算后，Runge-Kutta的累积误差导致计算结果的完全不同。<br>
<br>
<br>