基于风雷软件的壁面距离场的快速迭代并行求解方法

A parallel fast iteration method for wall distance implemented in PHengLEI

  • 摘要: 计算流体力学(computational fluid dynamics, CFD)中,网格的壁面距离场对雷诺平均Navier-Stokes(Reynolds-averaged Navier-Stokes, RANS)模型计算、重叠网格组装、富集函数方法等起到基础性作用。非结构网格上的距离场计算主要分为搜索方法、偏微分方程方法以及迭代方法等。快速迭代方法(fast iteration method, FIM)是一种适用于非结构网格的迭代算法,具有效率高、稳定性好的特点。本文基于国产开源CFD软件风雷平台(PHengLEI),探索并实现了这一算法。为适配风雷软件的数据结构并提高计算效率,本文将FIM方法改造为格心格式,设计了一种基于格心插值的几何对偶控制体算法,充分复用CFD并行数据结构;同时,为加速三维计算,本文将FIM的三维局部优化算法改为解析形式,提升了计算效率。在测试中,本文发现FIM方法与风雷软件原有ADT(alternating digital tree)方法相比,距离场求解误差很小,对基于SA(Spalart-Allmaras)模型模拟的CFD计算结果影响可忽略不计;同时,当网格规模超过千万且壁面网格数量超过数十万时,FIM方法相比ADT方法加速可达4倍。总体而言,FIM方法在保证计算精度的同时,在并行扩展能力与计算效率方面相比风雷软件原有方法具有明显的优势,具备良好的应用前景。本文所用代码与算例已上传至开源平台。

     

    Abstract: In CFD simulations, the wall-distance field plays a fundamental role in RANS-model computations, overset-grid assembly, enrichment-function methods, and related procedures. On unstructured meshes, distance-field computation typically falls into three categories: search-based methods, partial-differential-equation (PDE) methods, and iterative methods. Within the class of iterative methods, the Fast Iterative Method (FIM) is well suited to unstructured meshes and offers both high efficiency and good robustness. This work therefore presents the implementation of FIM on the domestic open-source CFD platform PHengLei. To accommodate PHengLei’s data structures and improve computational performance, we reformulate FIM into a cell-centered scheme, designing a parallel FIM algorithm that relies solely on the degrees of freedom located at cell centers by introducing a dual control volume based on cell-centered interpolation. In addition, to accelerate three-dimensional computations, we derive a closed-form analytical expression for the local optimization step in 3D FIM, further enhancing efficiency. Tests show that, compared with PHengLei’s original ADT (Alternating Digital Tree)-based method, the present FIM approach yields sufficiently small wall-distance errors and introduces negligible differences in CFD solutions based on the Spalart–Allmaras (SA) model. When the mesh contains more than ten million cells and the wall-cell count exceeds several hundred thousand, FIM achieves up to a fourfold speedup over the ADT method. Overall, the FIM method demonstrates acceptable accuracy while exhibiting substantially better parallel scalability and computational efficiency than PHengLei’s existing approach, showing promising potential for practical application. The source code and test cases used in this work have been made available on an open-source platform.

     

/

返回文章
返回