[论文解读] A fully discrete Calderon Calculus for two dimensional time harmonic waves

Sijiang L. Lu|arXiv (Cornell University)|Jan 1, 2014

Numerical methods in inverse problems参考文献 16被引用 11

一句话总结

本文提出了一种用于二维时 harmonic Helmholtz 方程的全离散 Calderón 微积分，采用交错网格、Dirac delta 分布表示源密度，并使用分段常数函数近似偶极子。当连续问题适定性成立时，所得格式达到二阶精度（O(h²)），并通过数值实验得到验证。

ABSTRACT

In this paper, we present a fully discretized Calderon Calculus for the two dimensional Helmholtz equation. This full discretization can be understood as highly non-conforming Petrov-Galerkin methods, based on two staggered grids of mesh size h, Dirac delta distri- butions substituting acoustic charge densities and piecewise constant functions for approxi- mating acoustic dipole densities. The resulting numerical schemes from this calculus are all of order h 2 provided that the continuous equations are well posed. We finish by presenting some numerical experiments illustrating the performance of this discrete calculus.

研究动机与目标

• 开发二维时 harmonic Helmholtz 方程 Calderón 微积分的全离散格式。
• 解决边界积分算子离散化过程中保持 Calderón 微积分结构的挑战。
• 在连续问题适定的前提下，确保离散格式具有高阶精度（O(h²)）。
• 为通过二维边界积分方程求解散射问题提供一个稳健的数值框架。

提出的方法

• 该方法采用两个网格尺寸为 h 的交错笛卡尔网格来离散化问题。
• 声学电荷密度通过放置在其中一个网格上的 Dirac delta 分布表示。
• 声学偶极子密度通过在第二个交错网格上的分段常数函数近似。
• 该方法被表述为一种高度非协调的 Petrov-Galerkin 方法，以确保稳定性和一致性。
• 通过将这些离散密度与适当的积分算子结合，构建离散 Calderón 微积分。
• 所得格式在保持连续 Calderón 微积分结构的同时，实现了完全离散化。

实验结果

研究问题

• RQ1 能否在保持底层数学结构的前提下，为二维 Helmholtz 方程构建一个全离散 Calderón 微积分？
• RQ2 在连续问题适定的前提下，所得离散格式的收敛速率是多少？
• RQ3 交错网格与非协调函数空间如何影响数值解的稳定性和精度？
• RQ4 能否通过 Dirac delta 分布与分段常数函数有效实现离散微积分中的偶极子近似？

主要发现

• 所提出的全离散 Calderón 微积分在连续 Helmholtz 问题适定的前提下，实现了二阶精度，误差被 O(h²) 所控制。
• 使用交错网格与非协调函数空间，实现了边界积分算子的稳定且一致的离散化。
• 该方法在离散设置中保持了 Calderón 微积分的数学结构，确保与连续理论的一致性。
• 数值实验验证了理论收敛速率，并展示了该格式在散射问题中的鲁棒性。