有限差分法解电磁场方程的共轭梯度法三角阵预处理器

摘要

通过考察电磁场微分方程经非均匀网格离散后的有限差分方程组矩阵,建立了函数偏微分运算与离散向量矩阵相乘运算的对应关系,给出了差分方程组矩阵对应于微分算子的分解式,并据此提出了共轭梯度法的三角阵预处理器。此外,还提出了对不同的边界条件、求解域内部边界、介质分界面和时谐场方程的处理技术以便应用该预处理器。数值计算结果验证了本文算法的正确性,展示了其十分明显的加速收敛效果,表明了本文算法有线性的存储复杂度和几乎线性的计算复杂度,可有较广泛的应用。本文中将算子细节和矩阵细节对应的基本思想对构造其它高效预处理器具有借鉴作用。