预条件共轭梯度法求解三维地电场有限元方程的网格分析

摘要

为了优化预条件共轭梯度法求解三维地电场有限元方程的效率,通过系统的网格分析,提出在常规六面体网格基础上二次剖分得到四面体网格的方案。模型分析结果表明,对于均匀网格,不完全Cholesky共轭梯度法(ICCG)和超松弛预条件共轭梯度法(SORCG)均可成功求解;对于非均匀网格,六面体剖分会导致ICCG的预条件因子不符合条件而求解失败,但采用新的四面体剖分ICCG不仅成功求解而且相对六面体网格可以节省约50%的内存需求。