我正在使用 Galerkin 方法(准确地说是不连续)在空间中离散化标量线性波动方程和显式二阶中心有限差分方案以在时间上离散化,从而产生一种半离散系统:
米ü¨+ ķü= 0 ,ü( 0 ) =ü0,ü˙( 0 ) =ü1,
在哪里米, ķ是 Galerkin 方法的质量和刚度矩阵。现在研究稳定性的想法是观察如果双线性形式一种从中ķ起源有一个L2-离散空间上特征函数的正交基Vh:a(wi,v)=λi(wi,v),∀v∈Vh,i=1…N,dimVh=N,
然后在使用时间积分方案并表示后uk=∑Nj=1ukjwj:
(un+1−2un+un−1,wi)+Δt2a(un,wi)=0
变为,使用以下事实(wj,wi)=δi,j:
un+1i=(2−λiΔt2)uni−un−1i,∀i=1…N
现在的重点基本上是求解这个二阶差分方程并要求它不会振荡或爆炸,这会导致这种 CFL 条件Δt≤Ch(绑定后λi与网格尺寸h)。
因此,作为二阶差分方程,可能存在判别为正、负(两个复根)或零的情况。这是进行此类稳定性分析的正确途径吗?这些要求是拥有 CFL 条件所需要的全部吗?非常感谢您的澄清。