摘要

作为近年来广受关注的一种数值方法,虚拟元方法具有很多优势。但在求解实际问题导出的一些辐射扩散方程时,该方法可能无法保证数值解的非负性及一般多边形网格上的局部守恒性。针对辐射扩散方程,利用非线性两点流逼近方法作为后处理措施,提出了一种基于虚拟元方法的保正守恒格式。该格式通过最低阶虚拟元方法得到数值解的单元顶点值,再利用非线性两点流逼近方法得到数值解的非负单元中心值,同时使格式满足局部守恒性。任意多边形网格上的数值结果表明,该格式具有保正性和解的近似二阶收敛速度,对于处理含强间断或非线性扩散系数的辐射扩散问题均有较强的适应性。