[论文解读] Towards Unstructured Mesh Generation Using the Inverse Poisson Problem

Guy Bunin|ArXiv.org|Feb 17, 2008

Computational Geometry and Mesh Generation参考文献 24被引用 4

一句话总结

本文提出了一种新颖的方法，通过求解逆泊松（IP）问题，在平面区域中生成非结构化四边形网格。该方法将网格生成建模为全局优化问题，其中不规则顶点对应泊松方程中的点源，解的结果生成具有用户定义分辨率和方向性的共形网格场。主要贡献在于提出了一套完整的网格生成方案，能够生成高质量、边界对齐的网格，具有平滑的过渡和良好形状的单元。

ABSTRACT

A novel approach to unstructured quadrilateral mesh generation for planar domains is presented. Away from irregular vertices, the resulting meshes have the properties of nearly conformal grids. The technique is based on a theoretical relation between the present problem, and the inverse Poisson (IP) problem with point sources. An IP algorithm is described, which constructs a point-source distribution, whose sources correspond to the irregular vertices of the mesh. Both the background theory and the IP algorithm address the global nature of the mesh generation problem. The IP algorithm is incorporated in a complete mesh generation scheme, which also includes an algorithm for creating the final mesh. Example results are presented and discussed.

研究动机与目标

• 开发一种理论基础坚实的全局方法，用于非结构化四边形网格生成，确保高质量的单元形状和平滑的分辨率过渡。
• 通过将网格生成与逆泊松（IP）问题关联，解决其全局特性，使网格属性从势场中自然涌现。
• 构建一个网格生成流水线，通过IP解实现不规则顶点的放置，随后执行最终的网格构造算法。
• 生成边界对齐的网格，其中单元质量和方向性由边界或关键点处的用户定义输入控制。

提出的方法

• 该方法将网格生成建模为求解逆泊松问题：给定目标势场 φ，重构源分布 ρ，使得 ∇²φ = ρ。
• 网格中的不规则顶点被建模为泊松方程中的点源（δ-函数），其位置和阶数（k-值）由IP解确定。
• 势场 φ 被解释为局部网格分辨率的对数，其梯度通过交叉场和测地线曲线定义网格方向性。
• 采用有理函数插值技术求解IP问题，近似用户指定点（包括边界）处所需的分辨率和方向性。
• 最终的网格构造算法从连续场生成有效的四边形网格，确保连通性与单元形状质量。
• 通过从k-值导出的角度偏移，计算奇异点（不规则顶点）处的星形测地线，确保奇异点周围网格边的正确对齐。

实验结果

研究问题

• RQ1 如何利用偏微分方程框架形式化非结构化四边形网格生成的全局特性？
• RQ2 逆泊松问题能否有效用于确定网格中不规则顶点的最优位置和类型？
• RQ3 在逆泊松框架中，通过点源分布在多大程度上可近似实现具有用户定义分辨率和方向性的共形网格场？
• RQ4 如何可靠地将所得连续场转换为具有边界对齐性的有效、高质量四边形网格？

主要发现

• 该方法成功生成了边界对齐的四边形网格，不规则顶点仅在必要位置放置，其余区域高度规则。
• 对于图6中的示例，粗网格（1898个单元）的最小β值为0.576，平均β值为0.945，表明单元形状良好。
• 当分辨率加倍（7630个单元）时，最小β值提升至0.627，平均值提升至0.975，表明分辨率提高后单元质量进一步改善。
• 在图8中，粗网格的最小β值为0.506，平均为0.933；更细的网格分别达到0.562和0.958，证实了质量的一致性提升。
• 该方法在不同域中均保持高质量网格，β值始终高于0.5，平均值接近0.95–0.98，表明单元形状接近正方形。
• 该算法对分辨率缩放具有鲁棒性，因为通过向φ中加常数来加倍分辨率不会改变源分布，从而保持了网格结构。