Crank Nicolson 不同离散规范背后的直觉

计算科学 pde 稳定 抛物线pde
2021-11-27 02:25:41

我正在用 Crank-Nicolson 有限差分法项的常用离散化矩阵。我想谈谈整个网格上的整个误差向量。我可以将该方法的离散迭代写为因为是正常的,所以我有并且我可以看到已知特征值我可以找到的特征值,因此我有其中中整个网格的矢量误差ut=AuAuxx

un+1=Bun=I+τ/2AIτ/2Aun
Bρ(B)=||B||2AB||B||21||e||2=O(h2)eL2规范。

因此,我的第一个问题是我是否也可以为 Crank-Nicolson为此,我需要估计而我不知道该怎么做。||e||=O(h2)||B||=maxi,jbi,j

第二个问题,如果第一个问题是一个真实的陈述,我对两个离散范数的误差进行了估计:L2L具有二阶误差向量的L^2离散误差的直觉是什么L2,在我看来,“平均”误差减少了 2,所以如果我选择网格上的点,我不必有二次收敛?我可以暗示任何关于它将如何在网格上的特定点收敛吗?关于L范数的相同问题。我应该使用哪一个来衡量错误,一个会暗示另一个吗?

1个回答

关于您的第一个问题:我假设“通常的离散化矩阵”是指一维中的 3 点有限差分离散化,或者您使用线性有限元得到的结果。无论哪种情况,实际上都不是 Crank-Nicolson 方案决定了这一点。确实,对于上面提到的两个空间离散化,空间误差,如果你选择那么整体误差也有这种收敛性命令。但这仅在 1d 中是正确的。O(h2)L2Lτ=h

另一方面,如果你在 2d 中,那么你会范数中错误,但在规范。这表明仅查看 Crank-Nicolson 方案的属性是不够的 - 事实上,还需要查看底层空间离散化的属性来证明您感兴趣的估计值。您的方法只是因此,只看的矩阵范数是不够的。O(h2)L2O(h2|logh|)LB

这也回答了你第二个问题的一部分。在任意数量的空间维度中,误差总是一样衰减误差并非如此。换句话说,第一个意味着平均而言,空间误差以二次方衰减,但这不一定是逐点的情况。另一方面,如果中衰减是真的,那么证明误差也一定是这样的就很简单了。L2O(h2)LLO(h2)L2