使用空间中的 Galerkin 方法和时间上的有限差分来找到 CFL 条件的技术?

计算科学 有限元 有限差分 稳定 cfl 加勒金
2021-12-26 10:17:58

我正在使用 Galerkin 方法(准确地说是不连续)在空间中离散化标量线性波动方程和显式二阶中心有限差分方案以在时间上离散化,从而产生一种半离散系统:

MU¨+KU=0,U(0)=U0,U˙(0)=U1,

在哪里M,K是 Galerkin 方法的质量和刚度矩阵。现在研究稳定性的想法是观察如果双线性形式a从中K起源有一个L2-离散空间上特征函数的正交基Vha(wi,v)=λi(wi,v),vVh,i=1N,dimVh=N,

然后在使用时间积分方案并表示后uk=j=1Nujkwj

(un+12un+un1,wi)+Δt2a(un,wi)=0

变为,使用以下事实(wj,wi)=δi,j:

uin+1=(2λiΔt2)uinuin1,i=1N

现在的重点基本上是求解这个二阶差分方程并要求它不会振荡或爆炸,这会导致这种 CFL 条件ΔtCh(绑定后λi与网格尺寸h)。

因此,作为二阶差分方程,可能存在判别为正、负(两个复根)或零的情况。这是进行此类稳定性分析的正确途径吗?这些要求是拥有 CFL 条件所需要的全部吗?非常感谢您的澄清。

1个回答

是的。这就是稳定条件的全部内容。

取材料特性——剪切模量(μ), 体积模量 (κ) 和密度 (ρ) - 考虑到,全局临界时间步被评估为每个元素的临界时间步的最小值 (Δte)

Δte=CFLhe/cκ

其中 CFL 是 Courant-Friedrichs-Lewy 数,he是元素的特征长度,并且cκ是体积波速,计算为

cκ=κ+4μ/3ρ.

有关更多详细信息和数值示例,请参阅我关于弹性动力学和波传播的论文。论文 1论文 2