[论文解读] Numerical conformal mapping via a boundary integral equation with the adjoint generalized Neumann kernel

Mohamed M. S. Nasser, Ali H. M. Murid+1 more|arXiv (Cornell University)|Aug 19, 2013

Numerical methods in engineering参考文献 30被引用 6

一句话总结

本文提出了一种统一的数值方法，用于从有界多重连通区域到五类经典标准缝合区域的共形映射、其导数及逆映射的计算，采用一种新颖的伴随广义诺伊曼核边界积分方程。该方法将映射问题重新表述为伴随黎曼-希尔伯特问题，推导出边界对应函数 θ′ 的唯一可解积分方程，并通过积分和同一方程计算 θ 和映射参数，从而实现稳定、精确的计算，无需奇异算子或多组方程。

ABSTRACT

This paper presents a new uniquely solvable boundary integral equation for computing the conformal mapping, its derivative and its inverse from bounded multiply connected regions onto the five classical canonical slit regions. The integral equation is derived by reformulating the conformal mapping as an adjoint Riemann-Hilbert problem. From the adjoint Riemann-Hilbert problem, we derive a boundary integral equation with the adjoint generalized Neumann kernel for the derivative of the boundary correspondence function $θ'$. Only the right-hand side of the integral equation is different from a canonical region to another. The function $θ'$ is integrated to obtain the boundary correspondence function $θ$. The integration constants as well as the parameters of the canonical region are computed using the same uniquely solvable integral equation. A numerical example is presented to illustrate the accuracy of the proposed method.

研究动机与目标

• 开发一种统一的数值方法，用于从有界多重连通区域到五类经典标准缝合区域的共形映射计算。
• 克服现有方法中需为不同标准区域分别制定公式或使用多组积分方程的局限性。
• 实现映射函数、其导数及其逆映射的同时计算，提升稳定性和精度。
• 消除因在 θ′(t) 接近零时计算 θ′(t)/|θ′(t)| 而引起的数值不稳定性。
• 提供一个右端项仅随标准区域变化的单一、唯一可解积分方程，简化实现与分析。

提出的方法

• 利用边界对应函数 θ′(t) 的导数性质，将共形映射重新表述为伴随黎曼-希尔伯特问题。
• 为 θ′(t) 推导出带有伴随广义诺伊曼核的边界积分方程，确保唯一可解性。
• 通过积分 θ′(t) 获得边界对应函数 θ(t)，积分常数通过同一积分方程计算。
• 映射参数 R(t) 也使用同一积分方程计算，确保一致性。
• 从 θ(t)、θ′(t) 和 R(t) 重构映射函数 ω(z) 及其边界值，以及 ω′(z) 的边界值。
• 使用柯西积分公式计算 ω(z) 的内部值，逆映射值则通过柯西型积分或牛顿迭代法计算。

实验结果

研究问题

• RQ1 能否为所有五类经典标准缝合区域，统一地推导出一个唯一可解的边界积分方程，以实现共形映射的统一计算？
• RQ2 如何在不依赖 θ′(t)/|θ′(t)| 等数值敏感项的前提下，稳定计算边界对应函数导数 θ′(t)？
• RQ3 伴随广义诺伊曼核在确保不同标准区域下解的存在唯一性并简化公式推导中起到何种作用？
• RQ4 如何利用单一积分方程框架，同时计算映射函数、其导数及其逆映射？
• RQ5 该方法能否避免先前方法中出现的奇异算子和多组方程，从而提升数值稳定性和效率？

主要发现

• 所提方法具有高精度，所有测试标准区域的数值误差均在 10^(-13) 量级，表明其收敛性稳健。
• 该方法仅使用一个积分方程即可同时计算映射函数、其导数及其逆映射，而以往方法需分别使用三个独立方程。
• 积分方程的右端项仅依赖于标准区域，简化了实现，实现了对五类区域的统一处理。
• 该方法避免了在 θ′(t) 接近零时计算 θ′(t)/|θ′(t)| 所引发的数值不稳定性，这是早期方法的常见问题。
• ω(z) 和 ω′(z) 的边界值可直接从 θ(t)、θ′(t) 和 R(t) 计算得出，无需额外的奇异算子或迭代校正。
• 该方法适用于单连通区域和多重连通区域，其适用范围超出了原始标准区域的设定。