间断Galerkin方案中的CFL条件


9

我已经实施了ADER-间断Galerkin方案,用于解决以下类型的守恒律线性系统 tU+AxU+ByU=0并观察到CFL条件非常严格。在参考书目中,时间步长的上限Δthd(2N+1)λmax 可以在哪里找到 h 是单元格的大小 d 是维数, N 是多项式的最大次数。

有什么办法可以避免这个问题?我一直在使用WENO-ADER有限音量方案,而CFL限制则更加宽松。例如,对于5阶方案,使用DG时必须将CFL设置为低于0.04,而CFL = 0.4仍可以在WENO-ADER FV方案中使用。

例如,为什么在计算航空气动(线性欧拉方程)或类似应用(气体动力学,浅水,磁流体动力学)中使用DG方案而不是ADER-FV?尽管时间步长低得多,但是该方案的总体计算成本是否与ADER-FV相似?

欢迎对此提出意见和建议。

Answers:


6

DG方案的限制性CFL通常来自高阶精度和紧凑型模板的结合(例如,请参阅此参考资料)。CFL依赖于以L2解的范数,它取决于多项式的导数和迹线。这些数量中的每一个的界限(使用Bernstein或Markov兄弟不等式和离散迹线不等式)给出的常数反比于h 并按顺序二次 N,导致整体CFL为 O(h/N2)

仅供参考-我已经看过您提到的CFL,但我不记得它的证明。我想知道他们如何避免二次依赖N 在他们的范围内。

有限差分和WENO方案(以及周期网格上基于B样条的有限元方法)具有较宽松的CFL条件,因为相似边界中的常数在 N。这又是因为模板尺寸倾向于随着顺序增加N,从而减少了其中一些问题。

DG方法更昂贵,但它们可以轻松处理非结构化网格,并且可以高效地实施。对于非结构化网格,存在WENO的高阶版本(或类似的重构),尽管它们会带来额外的数学或实现复杂性。


非常感谢您的详细回答杰西(Jesse),它为我提供了关于此问题的更广阔的视野。在DG-ADER的数值试验中,我注意到,当使用结构化四边形网格(具有任意四边形形状,例如正方形,梯形或平行四边形...)时,数值解是非振动性的,并且收敛于精确解。但是,当移动到非结构化网格时,即使对于准结构化网格,也会出现振荡,这种振荡是通过将结构化网格的节点随机移动一小段距离而产生的。这是预期的行为吗?
Adr

1
@Adrian-一旦远离均匀网格,就会出现振荡。一旦使用了通用网格,也将不再清楚网格尺寸的确切含义。h。它可以是像元直径,最短边的长度,面积的平方根(以2d为单位),或任何其他定义“网格大小”的方式。
Wolfgang Bangerth '17

感谢您对沃尔夫冈的澄清。到目前为止,我正在设定h作为最短边的长度。但是无论如何,即使将CFL的数量从公式给出的规定的CFL减少一个数量级或更多,仍然令人生厌。
2015年
By using our site, you acknowledge that you have read and understand our Cookie Policy and Privacy Policy.
Licensed under cc by-sa 3.0 with attribution required.