研究目的
为高效求解积分形式的时谐麦克斯韦方程组,开发并应用一种稀疏化预条件子,将标量情形的方法推广至矢量情形以降低计算成本。
研究成果
稀疏化预条件子能有效将求解时谐麦克斯韦方程组稠密积分系统的计算成本降至与稀疏偏微分方程系统相当的水平,其迭代次数基本不受问题规模和频率影响。未来工作可拓展至非均匀磁导率情形,并与扫描预条件子结合以进一步提升效率。
研究不足
该方法假设介质是平滑的,以便在均匀网格上进行精确的Nystr?m离散化并应用FFT;对于具有急剧过渡的介质,其性能可能不佳。计算成本虽有所改善,但在三维情况下仍涉及O(N^2)的设置时间,并且需要依赖特定软件(MATLAB)和硬件才能实现最佳性能。
1:实验设计与方法选择:
该方法采用均匀笛卡尔网格和Nystr?m离散化对时谐麦克斯韦方程的积分方程进行离散化,通过最小化非局部相互作用应用稀疏化预条件子将稠密系统转化为稀疏系统,并使用GMRES求解器进行迭代求解。理论模型包括亥姆霍兹方程和积分公式。
2:样本选择与数据来源:
假设介质为各向同性且具有非均匀介电常数的光滑介质,在紧凑区域内的均匀网格上进行离散化。示例包括二维和三维中的高斯透镜、光滑障碍物以及随机扰动。
3:实验设备与材料清单:
计算资源包括配备英特尔至强E7-4830-v3 CPU的服务器、用于实现的MATLAB软件以及用于高效卷积评估的FFT。
4:实验步骤与操作流程:
步骤包括离散化域、通过SVD计算模板α和β以消除非局部相互作用、形成稀疏系统、应用嵌套剖分法求解以及使用带预条件子的GMRES。参数如网格尺寸h设置为每个波长六个点。
5:数据分析方法:
通过迭代次数(Niter)、设置时间(T_setup)、应用时间(T_apply)和求解时间(T_solve)分析不同波数下的性能,并与理论复杂度进行缩放比较。
独家科研数据包,助您复现前沿成果,加速创新突破
获取完整内容