[论文解读] Theoretical Considerations on the Computation of Generalized Time-Periodic Waves

Dirk Pauly, Tuomo Rossi|arXiv (Cornell University)|May 20, 2011

Numerical methods in inverse problems参考文献 42被引用 7

一句话总结

本文提出了一种新颖的最小二乘控制方法，用于通过在微分形式中构建的一阶系统，计算广义波动方程的时间周期解。通过利用一种 $ L^2 $-型能量范数，基于一阶麦克斯韦型系统，该方法实现了无需预条件的共轭梯度最小化，实现了 $ O(n) $ 的计算复杂度，并在数值测试中展现出二阶误差减少的鲁棒收敛性。

ABSTRACT

We present both, theory and an algorithm for solving time-harmonic wave problems in a general setting. The time-harmonic solutions will be achieved by computing time-periodic solutions of the original wave equations. Thus, an exact controllability technique is proposed to solve the time-dependent wave equations. We discuss a first order Maxwell type system, which will be formulated in the framework of alternating differential forms. This enables us to investigate different kinds of classical wave problems in one fell swoop, such as acoustic, electro-magnetic or elastic wave problems. After a sufficient theory is established, we formulate our exact controllability problem and suggest a least-squares optimization procedure for its solution, which itself is solved in a natural way by a conjugate gradient algorithm operating in the canonical Hilbert space. Therefore, it might be one of the biggest advances of this approach that the proposed conjugate gradient algorithm does not need any preconditioning.

研究动机与目标

• 开发一种鲁棒且高效的方法，用于计算非均匀介质中广义波动方程的时间周期解。
• 将基于控制的波动方程方法推广至使用微分形式的统一框架，涵盖声学、电磁学和弹性波问题。
• 通过在 $ L^2$-型希尔伯特空间中重新表述问题，消除共轭梯度最小化中对预条件的依赖。
• 通过良好的条件性最小二乘公式，实现现代多核与分布式架构上的高效并行实现。
• 通过具有精确解的三维散射问题对方法进行数值验证，重点关注收敛性与准确性。

提出的方法

• 使用交替微分形式将波动方程表述为一阶系统，实现对声学、电磁学和弹性波问题的统一处理。
• 应用精确能控性理论，将时间周期解重新表述为一个二次泛函的最小化问题，该泛函衡量一个周期后初始状态与最终状态之间的差异。
• 在 $ L^2$-型希尔伯特空间中采用最小二乘公式，确保最小化问题条件良好，避免对预条件的需求。
• 直接在 $ L^2$-空间中实现共轭梯度算法，该算法天然无预条件，适用于并行计算。
• 使用带良好中心化单纯形网格和基于外心的对偶网格的离散外微分（DEC）进行系统离散化，以确保离散霍奇星算子为对角形式。
• 在微分形式框架中应用吸收边界条件（ABC）、广义狄利克雷-到-诺伊曼（DtN）算子和完美匹配层（PML），以实现域截断。

实验结果

研究问题

• RQ1 在微分形式中对波动方程的一阶表述是否能导致一个条件良好的最小二乘问题，从而避免预条件？
• RQ2 所提出的共轭梯度方法在实际中是否实现 $ O(n) $ 的计算复杂度，且与网格尺寸无关？
• RQ3 该方法在具有已知精确解的三维时间周期波动散射问题中的收敛性和准确性表现如何？
• RQ4 该框架是否可通过非结构化单纯形网格扩展至处理变系数介质和复杂几何形状？
• RQ5 残差与误差在共轭梯度算法中的收敛行为在不同网格细化下如何表现？

主要发现

• 共轭梯度算法在140次迭代后将残差降低了 $ 10^{-5} $ 倍，且收敛性与网格单元尺寸无关。
• 模拟电场的相对误差随网格细化而减小，表明其相对于平均边长具有二阶收敛性。
• 通过 $ L^2$-型能量范数实现了条件良好的最小化问题，消除了共轭梯度求解器中对预条件的需求。
• 在具有精确解的三维电磁散射问题上的数值结果表明，六次网格细化下均表现出一致的收敛性与准确性。
• 采用基于外心的对偶网格的离散外微分（DEC）实现成功处理了非结构化三维网格与异步时间推进。
• 该框架通过统一的微分形式公式，自然地兼容各种波问题——声学、电磁学、弹性波。