[论文解读] Shape derivatives of boundary integral operators in electromagnetic scattering

Martin Costabel, Frédérique Le Louër|arXiv (Cornell University)|Feb 8, 2010

Numerical methods in inverse problems参考文献 29被引用 10

一句话总结

本文通过边界积分方程建立了电磁散射解的形状可微性，证明解在障碍物边界之外无限次形状可微，而在边界上导数的正则性会降低。它将一阶形状导数表征为一个新散射问题的解，从而实现了利用积分方程方法进行电磁学中的基于梯度的形状优化。

ABSTRACT

We develop the shape derivative analysis of solutions to the problem of scattering of time-harmonic electromagnetic waves by a bounded penetrable obstacle. Since boundary integral equations are a classical tool to solve electromagnetic scattering problems, we study the shape differentiability properties of the standard electromagnetic boundary integral operators. Using Helmholtz decomposition, we can base their analysis on the study of scalar integral operators in standard Sobolev spaces, but we then have to study the Gâteaux differentiability of surface differential operators. We prove that the electromagnetic boundary integral operators are infinitely differentiable without loss of regularity and that the solutions of the scattering problem are infinitely shape differentiable away from the boundary of the obstacle, whereas their derivatives lose regularity on the boundary. We also give a characterization of the first shape derivative as a solution of a new electromagnetic scattering problem.

研究动机与目标

• 分析穿透性障碍物在时 harmonic 电磁散射中解的形状可微性。
• 在能量空间 $\boldsymbol{\mathsf{T}}\boldsymbol{\mathsf{H}}^{-\frac{1}{2}}(\operatorname{div}_{\Gamma},\Gamma)$ 上，建立作用于切向量场的电磁边界积分算子的 Gâteaux 可微性。
• 将一阶形状导数表征为具有修改边界条件的新电磁散射问题的解。
• 通过提供显式导数表示，实现基于梯度的形状优化在雷达和天线设计等应用中的应用。

提出的方法

• 使用 Hodge 分解将空间 $\boldsymbol{\mathsf{T}}\boldsymbol{\mathsf{H}}^{-\frac{1}{2}}(\operatorname{div}_{\Gamma},\Gamma)$ 分解为梯度与旋度分量，将分析简化为标量积分算子。
• 分析标准 Sobolev 空间中伪齐次核的 Gâteaux 可微性，将声学结果推广至电磁设置。
• 研究在 $\mathscr{C}^1$ 变形下表面微分算子（如 $\nabla_\Gamma$, $\operatorname{curl}_\Gamma$）的行为。
• 利用换位子恒等式与曲面向量微积分，推导电场与磁场一阶形状导数的显式表达式。
• 通过修改传输条件，构造一个新的边值问题，其解等于一阶形状导数。
• 利用积分表示与逆算子，将可微性通过解映射传播。

实验结果

研究问题

• RQ1 电磁边界积分算子是否对障碍物边界的 $\mathscr{C}^k$ 变形具有无限形状可微性？
• RQ2 介电散射问题的解在内部是否保持无限可微性，还是在边界附近正则性会退化？
• RQ3 电磁场的一阶形状导数能否被表征为具有修改边界条件的新散射问题的解？
• RQ4 在关于域变形的 Gâteaux 微分下，如 $\nabla_\Gamma$ 和 $\operatorname{curl}_\Gamma$ 等表面微分算子如何表现？
• RQ5 何种最优边界变形正则性可保持积分算子的 Gâteaux 可微性？

主要发现

• 电磁边界积分算子在区域内部无限次形状可微，且不损失正则性。
• 散射问题的解在内部无限次形状可微，但其导数在边界 $\Gamma$ 上正则性会降低。
• 电磁场的一阶形状导数被表征为一个具有涉及形变矢量场表面导数的修改传输条件的新散射问题的解。
• 电场的形状导数满足一个包含源项 $\nabla_\Gamma\xi$ 和 $\boldsymbol{\mathsf{n}}\times\partial_\mathbf{n}(\boldsymbol{\mathsf{E}}^i - \boldsymbol{\mathsf{E}}^s)$ 的 Maxwell 方程组。
• 对于 $\mathscr{C}^1$ 变形，边界积分算子的 Gâteaux 导数是良定义的，即使在 Lipschitz 域上也成立，但高阶导数需要完整的 $\mathscr{C}^2$ 正则性。
• Helmholtz 分解可将向量值问题约化为标量算子，从而在 Sobolev 空间中系统地分析可微性。