耦合势有限体积法模拟海洋移动导线源三维频率域电磁响应

摘要

本文采用耦合势有限体积法建立了一套快速计算海洋移动导线源电磁响应的三维数值模拟技术.基于电场矢势标势分解将Maxwell方程转换为关于矢势和标势的亥姆霍兹方程;采用洛伦兹规范整理得到对称形式的离散系数矩阵.对导线源采用直接离散方法,基于ILUT(Incomplete LU Preconditioner with Threshold)预处理的BICGSTAB(Biconjugate Gradient Stabilized)迭代算法实现离散线性方程的求解.为了减少移动源总的计算时间,提出采用解的初值重置技术加快正演模拟的速度,即利用发射源移动前的电磁场响应作为发射源移动后的离散方程的初始解.数值计算结果表明:对导线源采用直接离散方法,当导线源的长度大于源周围网格长度的4倍以上时,源附近离散误差显著减小;相比常规的直接迭代算法,采用解的初值重置技术能够减少迭代次数、提高正演计算速度.