基于自适应区域分解的电磁场有限元求解方法研究

为了高效利用现有并行计算资源、提高电磁场数值分析效率,提出一种自适应区域分解有限元法用以电磁场求解。该方法是在加性Schwarz区域分解法中引入基于度量张量的各向异性网格自适应方法,将h型自适应有限元方法与区域分解方法的优点结合起来,在求解过程中生成对场分布响应更好的网格。为避免子区域独立网格自适应可能导致的悬点,在每步求解结束后对整个定义域的网格进行统一调整。计算过程中的多次子区域划分通过METIS实现,以基于当前网格对计算量在多处理器间进行合理分配。为进一步提高求解效率,在代数方程求解中将区域分解法用作Krylov子空间法的预处理算子。最后通过两个典型数值算例对所提方法在精度和效率方面的性能做了验证。结果表明:所提方法在保证求解精度的同时,所需自由度数目有显著降低(两个算例分别降低50%和36%),计算效率得到有效提高,验证了所提方法的有效性。