二维交错网格上新型间断伽罗金有限元方法-new discontinuous galerkin finite element method on two-dimensional staggered grid.docxVIP

Makridakis[17]，Machiels[18]，Estep和Freund[19]，Schwab[10]等人也进行了一些后续的研究和探索。1989年，DG方法首次被使用在对粘弹性流的数值计算上[20]。随后Baaijiens，Boga- erds和Verbeeten[24]探讨了该方法在粘弹性流体分析使用中的优缺点。Baranger和Sandri[22]， Wardi[23]，Mamoum[25]，Bahhar[26]等人对该领域的DG方法进行了分析研究。DG方法在线性双曲问题上的优点促使人们探索用其解决非线性双曲问题:diut + γ(fi(u))x = 0i=1这里为了简化，假设u是标量函数，给上式两边同乘以试验函数v，并在K上积分，采用分部积分原理，得到下式:ddi(ut, v) ? γ(fi(u), ?x v)K + γ(fi(u)(nK )i, v)?K = 0i=1近似解uh是满足如下弱形式的解:?K ∈ Th :i=1i=1i((uh)t, v)K ? )，d(fi(uh), ?x v)K + (h?, v)?K = 0, ?v ∈ P k (K)i=1这里h?是)，dfi(u)(nK )i在K边界上的近似。1982年，Chavent和Salzano[27]首次采用如上的D-G空间离散建立起解决非线性守恒律方程的框架。这之后所剩下的主要问题就是如何对 上面的方程在时间上进行离散。如何找到一种时间离散并最终导出一种稳定，高效并且 具有高阶精度的的DG方法，这是DG解决双曲守恒律方程的主要困难。对于全局的时间离散，由于非线性部分fi的出现，使得逐元地计算数值解变得非常困 难，因此不得不采用隐式的时间离散，并且在每一个新的时间步要求解一个新的非线性方 程组。这就使得数值计算双曲方程时效率比较低。Reed和Hill[1]，Bar-Yoseph[28]，Elata[29] 等人在这方面进行了一些尝试和探索。而对于局部的时间离散，通过Richter[31]，Lowrie[30]，Roe[32]等人的研究，发现通过 建立时空元，使得局部逐元地计算有可能实现。1990年，Hulbert和Hughes[33]提出了一种 空间时间有限元方法，其中在时间上使用了DG方法，在空间上使用了连续的空间近似。 为发展出完整的时空DG方法求解双曲方程提供了一种思路。对于非线性情形下的DG方法的分析比较困难，目前为止著述不多，第一个结果是在1994年，由Jiang和Shu[34]证明的:对一维的满足严格凸(或凹)非线性的守恒律方程时，如果DG方法所得的近似解收敛，它会收敛到脑解，对任意阶的逼近多项式都成立。另外两个结果是对求解标量非线性守恒律方程的空间时间DG方法进行分析得到的。为了避免隐式时间离散所带来的不便，1982年，Chavent和Salzano[27]对一维标量守恒 律方程建立了显式的DG方法，在空间上使用DG 方法的离散，在时间上采用向前的Euler方 法，但是该方法在 ?t 为常数时，是无条件不稳定的，而在 ?t 与√?x同阶时，该方法是稳?x?x定的，对于双曲方程而言，这是一个很严格的条件。为了改进上述方法的稳定性，在吸取了van Leer[35]的思想的基础上，Chavent和Cockbu- rn[36]在原方法中引入了slope限制器，对上述方法进行了修正，从而得到了新的方法，并?x≤时该方法还是TVB的。虽然数值结果2证明了该方法具有TVD性质，而且在满足f l ?t1表明该方法的近似解收敛到脑解，但是该方法在时间上只有一阶精度。并且slope 限制器 会对由线性不稳定性引起的光滑区域中的伪振荡现象进行调整，平衡，因此对这些区域 的逼近程度产生了不利的影响。3随着Runge Kutta Discontinuous Galerkin(RKDG)[37]方法的提出，上面的问题得到了 解决，该方法仍然使用分片的线性的DG方法对空间离散，使用显式的TVD的2阶Runge Kutta形式，修正slope限制器来保证该方法在极值处的正常精度。该显式方法被证明在CFL 数小于 1 时是线性稳定的，并且在时间和空间上都达到了2阶精度，是TVB的。在1989年，Cockburn和Shu[3]对标量的双曲守恒律方程，建立了具有高阶精度的RKDG 方法:采用了k阶分片多项式的DG空间离散，TVD且具有k +1阶精度的Runge-Kutta时间离 散和一般化的slope限制器。数值结果表明，对于k = 1, 2，在光滑区域以及在大幅震荡出现 的区域，都达到了k+1阶精度，并且都收敛到了脑解。同年，RKDG方法在Cockburn，Lin和Sh- u[4]的努力下推广到了1维方程组。并且在1990年，RKDG方法被推广到了标量多维情