一种基于局部间断Galerkin方法的IC互连线电容提取策略

求解椭圆方程的局部间断Galerkin(LDG)方法具有精度高、并行效率高的优点,且能适用于各种网格。文章提出采用LDG方法来求解IC版图中电势分布函数满足的Laplace方程,从而给出了一个提取互连线电容的新方法。该问题的求解区域需要在矩形区域内部去掉数量不等的导体区域,在这种特殊的计算区域上,通过数值测试验证了LDG方法能达到理论的收敛阶。随着芯片制造工艺的发展,导体尺寸和间距也越来越小,给数值模拟带来新的问题。文章采用倍增网格剖分方法,大幅减小了计算单元数。对包含不同数量和形状导体的七个电路版图,用新方法提取互连线电容,得到的结果与商业工具给出的结果非常接近,表明了新方法的有效性。

对流扩散方程的隐式全离散局部间断Galerkin方法

研究了对流扩散方程的隐式全离散局部间断Galerkin方法的稳定性和误差分析.将三阶隐式Runge-Kutta时间离散和具有广义交替数值流通量的LDG方法相结合得到全离散LDG格式,通过广义交替数值流通量,建立数值解和辅助解内积之间的关系,证明了全离散LDG格式的无条件稳定,同时引入广义Gauss-Radau投影,通过投影的逼近性质和一些基本不等式建立了数值方法的最优误差估计,最后通过数值实验验证该方法理论分析的正确性.

Caputo型时间分数阶变系数扩散方程的局部间断Galerkin方法

提出一种带有Caputo导数的时间分数阶变系数扩散方程的数值解法.方程的解在初始时刻附近通常具有弱正则性,采用非一致网格上的L1公式离散时间分数阶导数,并使用局部间断Galerkin(local discontinuous Galerkin,LDG)方法离散空间导数,给出方程的全离散格式.基于离散的分数阶Gronwall不等式,证明了格式的数值稳定性和收敛性,且所得结果关于α是鲁棒的,即当α→1^(-)时不会发生爆破.最后,通过数值算例验证理论分析的结果.

非定常线性双曲方程的差分——间断Galerkin方法的精度估计

本文对非定常线性双曲方程的差分间断———Ｇａｌｅｒｋｉｎ 格式进行了误差分析，作了收敛阶估计，得到了Ｌ∞（０ ，Ｔ；Ｈｒ＋ １（Ω）） 模的次最优估计。

非定常线性双曲方程的差分──间断Galerkin方法及其稳定性

本文针对非定常线性双曲方程建立了对时间进行差分离散的间断──Galerkin格式，并对该格式进行了稳定性分析。

abcd Boussinesq系统的高阶中心间断Galerkin-有限元方法

Boussinesq类型的方程是一类重要的非线性色散浅水波方程,广泛应用于研究小振幅浅水波的传播现象。我们考虑一个含有四个参数的Boussinesq系统:abcd Boussinesq系统,并针对该系统设计了一个高阶中心间断Galerkin-有限元方法。该方法首先将abcd Boussinesq系统改写为守恒律方程与椭圆型方程的一个耦合系统,然后采用中心局部间断Galerkin方法求解守恒律方程,使用有限元方法求解椭圆型方程,最后我们通过一系列的数值算例验证所提出方法的高阶精度和有效性。

基于雅可比矩阵精确计算的GMRES隐式方法在间断Galerkin有限元中的应用

为改善高阶间断Galerkin有限元方法(DG)时间推进效率,在三维非结构网格下针对该方法建立了并行广义最小残差(Generalized Minimal Residual,GMRES)隐式时间迭代方法,GMRES方法基于科学计算工具包PETSc中的Krylov子空间求解器实现。为进一步提高GMRES的计算效率,发展了方程组右端项残值雅可比精确计算方法,针对无黏通量Roe格式和黏性通量BR2(Bassi Rebay 2)黏性计算方法,分别解析给出其对守恒变量多项式自由度的雅可比矩阵。基于建立的方法首先采用NACA0012翼型研究了GMRES的重启次数及收敛参数对方法收敛性影响,然后采用无黏及黏性算例对比研究了基于雅可比矩阵不同计算方法的GMRES计算效率,同时对比研究了雅可比矩阵完全近似求解下GMRES和LU-SGS(Lower Upper-Symmetric Gauss-Seidel)的计算效率。结果表明,建立的基于右端项残值雅可比矩阵精确求解的GMRES方法能够大幅提高不同精度DG方法的CFL(CourantFriedrichs-Lewy)数,相比前面提到的其它方法具有更高的计算效率,其收敛速度实现量级以上的提高。

Navier-Stokes方程间断Galerkin有限元方法研究

通过引入全局提升算子和局部提升算子,发展了求解Navier-Stokes方程的间断Galerkin(discontinuousGalerkin,DG)有限元方法的一般框架,并在此框架下给出了几种典型黏性离散格式的具体表达形式.对局部提升算子的求解给出了详细的计算步骤.同时还给出了一种简单有效的计算方法来对物面边界进行高阶近似.为了能够对NS方程进行精度测试,采用对原始系统添加源项的方法构造精确解.二维Euler和NS系统的精度测试表明该方法达到了DG方法的理论精度.二维圆柱无黏绕流的计算结果表明关于物面边界的高阶近似方法能够保持DG方法原有的精度.卡门涡街数值模拟则进一步验证了该方法的正确性并且显示出DG方法较高的计算精度和分辨率.

求解RANS方程的高阶间断Galerkin方法研究

基于二维结构网格,对间断Galerkin方法(DGM)求解雷诺平均Navier-Stokes(RANS)方程进行了研究。根据混合有限元方法的思想,引入辅助变量,对粘性项中的高阶导数项进行降阶处理,再应用DGM对辅助变量方程和降阶后的RANS方程进行联立求解。对平板层流流动进行了数值模拟,结果表明,随着逼近多项式次数的增加,速度型和摩擦系数更加接近Blasius解。此外,还引入了Bald-win-Lomax湍流模型,对NACA0012翼型跨音速粘性流场进行了数值模拟,与实验结果进行了对比,进一步验证了算法的可靠性。

一维溃坝问题的间断Galerkin方法

将间断Galerkin方法应用于一维溃坝问题中,采用了两种不同数值流函数,对其进行比较。使用TVDM,TVBM型的限制器,同时利用一种改进的限制器来消除振荡。与一种高分辨率TVD性质的差分方法做了比较。结果表明间断Galerkin方法得到了更加逼真的图像。

间断的Galerkin方法在地震波场数值模拟中的应用概述

通过数值求解描述地震波传播的微分方程来模拟波的传播角度而言,在日趋繁荣的地震波场数值模拟方法的探索中,间断的Galerkin方法(Discontinuous Galerkin,DG)是继有限差分法、伪谱法、有限元法以及谱元法之后出现的另一种有效的数值模拟方法,该方法从2005年以来在地震波场模拟中得到了迅速发展,给这一学科注入了新的活力.但是从国内的发展来看,在地震波场数值模拟方面还没有得到足够的关注.尤其是现有的综述性文献都忽略了这一方法得到的有益的成果,本文将重点介绍间断的Galerkin方法的发展及其在地震波场数值模拟的最新应用研究.

对流占优的扩散问题的局部间断Galerkin方法

针对具有周期性边界条件对流占优的扩散问题中的二阶导数,引入辅助变量,构造了局部间断Galerkin(LDG)方法,并给出了方法的稳定性结果和误差估计式.局部间断Galerkin方法是Runge-Kutta间断Galerkin方法的推广,具有高阶精度,能够灵活处理复杂区域,易于处理复杂边界的边值问题,能够有效去除近似解在间断、大梯度处产生的虚假振荡.数值实验说明,当有限元空间取为一次多项式空间时,LDG方法具有二阶收敛,误差满足理论估计式.该方法可以推广到更高阶的方程,如Korteweg-de Vries方程、重调和方程等.

基于局部间断Galerkin方法的p型有限元

将基于三变量能量原理的局部间断Galerkin方法 (localdiscontinuousGalerkin ,LDG)应用于 p型单元的构造 .该方法采用间断的单元试解 ,不需要满足普通有限元所必须的协调条件 ,就能使构造高阶的插值函数变得更加灵活和容易 .在此基础上 ,对应力和应变场采用Legendre正交多项式进行插值 ,避免了柔度矩阵的求逆过程 .数值算例表明这种方法构造出的 p型单元不仅升阶过程简单 ,而且具有较高的精度 .

二维热传导方程的局部间断Galerkin有限元方法

讨论了局部间断Galerkin有限元方法求解二维热传导方程。通过引入辅助变量将含有二阶导数的热传导方程重新写为一阶偏微分方程组,在空间上用间断有限元离散得到一组常微分方程组,在时间上用显式方法离散,最终给出数值算例验证了该方法的收敛精度。

间断探测器在间断Galerkin方法中的应用

在二维非结构网格上,对高阶间断Galerkin方法求解定常二维欧拉方程进行数值研究。为了很好地抑制解在流场间断附近处产生的伪振荡,采用了间断探测器和限制器相结合的方法,并分别对Krivodonova间断探测器和基于激波物理特征的激波探测器进行了比较和研究。数值结果表明:如果流场中只存在激波时,两者探测效果相似,且后者更适合应用于求解激波问题。如果流场更加复杂,即存在多种间断时,前者依然可以比较准确地用来探测间断单元,能够很好抑制流场间断附近处产生的伪振荡,而后者失效。

非线性对流扩散问题的hp-局部间断Galerkin有限元方法

本文讨论了非线性对流扩散方程的hp-局部间断Galerkin有限元方法,得到了关于网格尺寸h是最优以及关于逼近次数p是次优的hp-误差估计.本文给出了数值算例,数值结果表明,LDG方法比SIPG或NIPG方法更精确且LDG方法对于罚值不敏感,而SIPG及NIPG方法对于罚值敏感.

半线性分数阶扩散方程的时空有限元方法:间断Galerkin方法(英文)

本文研究半线性分数阶扩散问题的Galerkin时空有限元方法,该方法在空间连续,而在时间上间断.将有限元与有限差分方法相结合,充分利用拉格朗日插值多项式在Radau点处的特性,给出弱解的存在唯一性证明,且不需对时空网格施加任何限制.通过引入椭圆投影算子,详细导出了最优阶L∞(L2)模误差估计.

显式Runge-Kutta局部间断Galerkin方法的稳定性分析

针对二阶显式TVD Runge-Kutta局部间断Galerkin方法求解热传导方程的稳定性问题,在方程的解是充分光滑的情况下,通过有限元分析的方法,在理论上严格的证明了对于任意非均匀正则网格和k次分段多项式间断有限元空间,当Courant-Friedrichs-Lewy(CFL)条件取为τ≤λμ-2h2时,算法是L2稳定的,其中τ,h分别是时间步长和空间步长,μ,λ是与h,τ无关的常数。

一种推广的对流扩散方程的局部化间断Galerkin方法

研究求解一种产生于径向渗流问题的推广的对流扩散方程的局部化间断Galerkin方法,对一般非线性情形证明了方法的L2稳定性;对线性情形证明了,当方法取有限元空间为k次多项式空间时,数值解逼近的L∞(0,T;L2)模的误差阶为k.

双曲型守恒律方程的自适应间断Galerkin方法

工程实际中的许多间断问题,例如空气动力学中的激波问题,其数学模型大都是非线性双曲守恒律方程.本文在Runge-Kutta间断Galerkin(RKDG)框架下,结合h型自适应方法处理了一维非线性守恒律方程初值问题和初边值问题.此方法不仅能准确描述间断的出现和位置,而且还能在间断附近适当加密网格,提高计算效率.最后,数值算例验证了算法的有效性.