[论文解读] Creating domain mappings

Kendall Atkinson, Olaf Hansen|arXiv (Cornell University)|Jun 16, 2011

Advanced Numerical Methods in Computational Mathematics参考文献 11被引用 11

一句话总结

本文提出了一种数值方法，用于在给定边界对应关系的情况下，从单位球到 ℝᵈ 中有界且单连通的区域构造一个光滑且可逆的映射。该方法使用调和延拓作为初始猜测，并通过优化改进以确保单射性和非奇异雅可比行列式，在二维和三维星形区域中取得了成功结果。

ABSTRACT

Consider being given a mapping ϕfrom the unit sphere S^{d-1}, d>2, to the smooth boundary of a simply-connected region Ωin R^d. We consider the problem of constructing an extension Φfrom the unit ball B_d to Ω. The mapping is required to be 1-1 and continuously differentiable with a nonsingular Jacobian matrix. We discuss ways of obtaining initial guesses for such a mapping Φand of then improving it by an iteration method.

研究动机与目标

• 在给定边界对应关系 φ 的前提下，构造从单位球 B_d 到 ℝᵈ 中有界且单连通区域 Ω 的连续可微双射映射 Φ。
• 确保 Φ 的雅可比行列式在 B_d 内处处非零，从而保证可逆性和光滑性。
• 开发一种迭代优化框架，以在保持边界条件和单射性的同时改进 Φ 的初始猜测。
• 通过误差度量 E₁,K 和 E₂,K,L 对方法进行数值验证，以评估单射性和像的包含性。
• 将该方法从二维推广到三维，使用球谐函数和 L² 空间中的正交投影，在高维中实现逼近。

提出的方法

• 通过极坐标或球坐标系中的傅里叶级数展开，利用边界数据 φ 的调和延拓生成初始映射 Φ₀。
• 将初始映射投影到调和多项式空间 Πₙ 的有限维子空间中，以确保光滑性和可微性。
• 应用优化过程以改进初始猜测，最小化一个惩罚非单射性和边界误差的泛函。
• 使用误差度量 E₁,K(Φₙ) 和 E₂,K,L(Φₙ) 监控单射性和像的包含性，其中 K 和 L 控制采样分辨率。
• 在 L² 空间中使用正交投影，从 φ 的非光滑径向延拓中计算初始近似 Φₙ⁽⁰⁾。
• 在三维情况下，通过球坐标系推广径向延拓，并将 ρ(θ,ϕ) 定义为 S² 上的光滑正函数。

实验结果

研究问题

• RQ1 当仅给出边界对应关系 φ: S^{d−1} → ∂Ω 时，如何构造一个从 B_d 到 Ω 的光滑可逆映射 Φ: B_d → Ω？
• RQ2 何种初始近似策略可确保可微性并满足边界条件，同时避免原点处的奇点？
• RQ3 优化技术能否改进初始映射，使其满足单射性和非奇异雅可比行列式条件？
• RQ4 如何利用数值度量 E₁,K 和 E₂,K,L 验证所得映射的单射性和像的包含性？
• RQ5 该方法在具有复杂边界形状的三维星形区域中，其成功程度如何？

主要发现

• 对于由 ρ(θ,ϕ) = 2 + cos²θ 定义的星形区域，初始映射 Φ₆⁽⁰⁾ 已经是局部最优解，其中 Λ(Φ₆) = 3.0576，E₁,₄₀(Φ₆) = 0.7486，E₂,₄₀,₁₀(Φ₆) = 0.6626。
• 对于更复杂的边界 ρ(θ,ϕ) = 2 + cosθ + ½ sinθ sinϕ，优化将 Λ(Φ₆) 从 394.37 降低至 43.88，表明在光滑性和单射性方面有显著改进。
• 度量 E₁,₄₀ 和 E₂,₄₀,₁₀ 均保持远离零（例如 E₁ ≈ 0.20），表明未违反单射性或像的包含性。
• 两个示例的边界误差均为零，证实边界条件 φ = Φ|_{S^{d−1}} 被精确满足。
• 初始猜测 Φₙ⁽⁰⁾ 虽在原点不可微，但通过投影进入 Πₙ 后，获得了光滑且可用的优化起点。
• 该方法成功生成了雅可比行列式不为零的映射，并在所有测试的二维和三维示例中保持了单射性。