四边形网格上变系数抛物型方程有限体积元法的L2模误差估计

在h2拟平行四边形条件下,给出变系数抛物型方程一个半离散和两个全离散有限体积元格式的L2模最佳收敛阶误差估计.     

四边形网上双曲型方程的有限体积元法

有限体积元法已引起国内外学者和专家的广泛关注,该方法与有限元方法有着相同的收敛阶,具有计算简单,保持物理量守恒性等优点.讨论一类双曲型方程在四边形网上的有限体积元法,在四边形网格单元满足h2拟平行四边形条件下,给出了双曲型方程半离散有限体积元格式下最优的H1模和L2模误差估计以及两个全离散有限体积元格式下的误差估计.     

一类Poisson-Nernst-Planck方程的边平均有限元计算

针对一类三维Poisson-Nernst-Planck方程, 给出一种边平均有限元离散形式. 在适当的网格条件下, 该离散形式得到的总刚度矩阵为M-矩阵, 从而保证了数值解的非负性. 数值实验结果表明, 边平均有限元方法相比于标准有限元的CPU时间更短, 且误差较小.     

一类Poisson-Nernst-Planck方程的两网格有限元离散方法

应用两网格有限元方法离散求解一类Poisson-Nernst-Planck(PNP)方程. 通过两网格离散, 将耦合PNP系统解耦成较小规模的线性对称系统, 可有效降低计算复杂度. 理论结果表明, 线性对称化的两网格算法具有与传统有限元方法相同的误差阶; 数值结果表明, 相比于传统有限元方法, 该方法计算效率更高.     

标准Criss-Cross剖分下线性有限元方程的快速AMG算法

首先对标准Criss-Cross剖分下的线性有限元空间进行能量正交分解,通过对正交子空间的双尺度分析,获得了一种合适的限制算子,进而构造相应的AMG算法.数值实验结果表明,该方法对求解椭圆方程是非常有效和健壮的,且与通常的代数多重网格法相比较,该AMG算法具有明显的优势.     

双曲型方程的基于外心对偶剖分的有限体积元法

讨论双曲型方程的基于外心对偶剖分的有限体积元法.设原始三角形剖分的任意三角形单元的重心Q和外心C的距离满足|QC|=O(h2),在此条件下,给出了双曲型方程半离散有限体积元格式最优的H1模和L2模误差估计以及两个全离散格式下的误差估计.     

一种三维Poisson-Nernst-Planck方程的虚单元计算

利用虚单元方法在多面体网格上求解一种三维稳态Poisson-Nernst-Planck(PNP)方程,并给出PNP方程的虚单元离散形式,推导电势方程及离子浓度方程的刚度矩阵与荷载向量的矩阵表达式.数值实验结果表明,在3种不同的多面体网格下实现了PNP方程的虚单元计算,数值解在L²和H¹范数下均达到最优阶.     

二阶非线性椭圆问题的自适应有限元方法

有限元方法的计算效率在很大程度上取决于离散网格的好坏,自适应有限元方法能够根据单元上解的误差而自动生成合适网格。对于一类非线性椭圆问题,研究了一种基于重构型后验误差估计子的自适应有限元算法。数值实验表明,这一算法对于求解这类非线性椭圆问题是十分有效的。     

双曲方程基于BB型对偶剖分的有限体积元法

基于三角形剖分和BB型对偶剖分,构造双曲方程半离散及两种全离散的有限体积元法,其中双曲方程的两种全离散格式分别用Grank-Nicolson和向后Euler格式逼近,得到并证明了双曲方程半离散有限体积元格式下最优的H1模和L2模误差估计及两种全离散格式下的误差估计.