隐晦的胡言乱语

亲爱的所有,

我写了一个基于隐式crancknicolson公式的有限差分编码的一维扩散方程(带时间和空间的抛物型偏微分方程)的代码。该问题是为了解决具有初始和边界条件的降雨引起的土坡稳定问题。我正在将我的结果与已经发表的工作进行比较，但是，在6小时的持续时间内，结果在某种程度上是匹配的。当模拟时间超过6小时时，观察到的结果差异大于20%。我检查了公式，初始条件和边界条件，看起来没问题。我一直保持稳定的单位。

请告诉我哪里错了。曲尼康逊隐式方法(中心差分法)是否存在稳定性和收敛性问题?我甚至保持了时间步长与空间步长的平方比几乎等于1。

请建议我对任何稳定性问题与CN方法或任何方法，可以给出精确的结果或有什么我可以改变我的代码。

提前谢谢你。