[论文解读] Conformal modulus on domains with strong singularities and cusps

Harri Hakula, Antti Rasila+1 more|arXiv (Cornell University)|Jan 27, 2015

Advanced Mathematical Modeling in Engineering参考文献 3被引用 9

一句话总结

本文提出了一种hp-FEM方法，用于数值计算具有强奇点和尖点的环形域与四边形域的共形模，将问题简化为求解拉普拉斯方程的狄利克雷问题和狄利克雷-Neumann边界值问题。该方法实现了近乎最优的收敛率，并表明辅助误差估计与互惠误差估计能可靠地跟踪真实误差，从而实现了对具有已知精确解的几何复杂域上FEM软件的精确基准测试。

ABSTRACT

We study the problem of computing the conformal modulus of rings and quadrilaterals with strong singularities and cusps on their boundary. We reduce this problem to the numerical solution of the associated Dirichlet and Dirichlet-Neumann type boundary values problems for the Laplace equation. Several experimental results, with error estimates, are reported. In particular, we consider domains with dendrite like boundaries, in such cases where an analytic formula for the conformal modulus can be derived. Our numerical method makes use of an $hp$-FEM algorithm, written for this very complicated geometry with strong singularities.

研究动机与目标

• 开发一种鲁棒的数值方法，用于计算因强几何奇点和尖点而使标准FEM方法因需要极高网格加密而失效的域中的共形模。
• 引入一族由(r, m, p)参数化的环形域，可在保持解析已知共形模值的同时提供可扩展的计算挑战。
• 在缺乏已知解析解的情况下，评估两种误差估计器——辅助空间误差与互惠误差——对真实误差的跟踪可靠性。
• 展示hp-FEM方法在具有角奇点和尖点的域上实现近乎最优的收敛率，验证Babuška等人理论预期的正确性。
• 为FEM软件提供一个基准测试套件，结合几何复杂性与精确解，以供验证。

提出的方法

• 将共形模的计算简化为在具有尖点和奇异边界的域上求解拉普拉斯方程的狄利克雷问题和狄利克雷-Neumann问题。
• 采用hp-FEM算法，利用高阶多项式逼近（p-型加密）和局部网格加密（h-型加密）以处理几何奇点。
• 借助Babuška及其合作者的理论，即使在角奇点存在的情况下，也能实现近乎最优的收敛率。
• 在缺乏精确解时，使用先前工作中提出的互惠误差估计作为真实误差的可靠替代指标。
• 从hp-空间理论推导辅助空间误差估计，并与真实误差和互惠误差进行比较。
• 应用莫比乌斯变换将无界或复杂域映射到单位圆盘以进行数值计算，同时保持共形结构。

实验结果

研究问题

• RQ1 hp-FEM方法是否能在具有尖点和强奇点的域上实现共形模计算的近乎最优收敛率？
• RQ2 在缺乏精确解的情况下，辅助空间误差估计与互惠误差估计是否能可靠地跟踪真实误差？
• RQ3 具有递归树状边界的环形域族是否可作为具有已知精确模值的可扩展基准？
• RQ4 当引入如尖点或小角度等几何特征时，hp-FEM方法的收敛行为如何变化？
• RQ5 在实际FEM应用中，互惠误差估计与辅助误差估计在多大程度上可作为真实误差的实用替代？

主要发现

• hp-FEM方法在具有尖点和奇点的域上实现了近乎最优的收敛率，与Babuška及其合作者的理论框架一致。
• 在所有测试案例中，互惠误差估计与辅助空间误差估计的表现与真实误差完全一致，验证了其作为可靠替代指标的有效性。
• 对于在s = t = √2 − 1 ≈ 0.41421处具有尖点的四边形，使用hp-FEM方法（p = 21，误差编号 = 10）计算得到的模为2.7823418086，与Bergweiler和Eremenko的结果在6位有效数字内一致。
• 在存在小几何特征的情况下，收敛性从指数收敛转变为代数收敛，表明由于奇点的存在导致谱收敛性的丧失。
• 由(r, m, p)参数化的新型环形域族提供了具有已知精确模值的可扩展基准测试套件，适用于FEM软件的验证。
• 该方法成功计算了具有类似尖点奇点的四边形的模，与Bergweiler和Eremenko的先前数值结果保持一致。