首页 >

适用于非结构混合网格的改进线隐式BLU-SGS迭代算法

张耀冰1,李明2,张健1,周乃春2,龚小权1   

  1. 1. 中国空气动力研究与发展中心
    2. 中国空气动力研究与发展中心计算空气动力研究所
  • 收稿日期:2025-07-17 修回日期:2025-08-26 出版日期:2025-08-28 发布日期:2025-08-28
  • 通讯作者: 张健
  • 基金资助:
    国家自然科学基金

An improved line-implicit BLU-SGS iteration algorithm for unstructured hybrid grids

  • Received:2025-07-17 Revised:2025-08-26 Online:2025-08-28 Published:2025-08-28
  • Contact: Jian ZHANG

摘要: LU-SGS方法作为当前计算流体力学(CFD)中的主流时间推进算法,其通过采用谱半径近似雅可比矩阵来降低计算开销,但同时也牺牲了收敛效率。尤其在基于非结构混合网格的高雷诺数流动模拟中,由于附面层网格在壁面法向方向的加密导致数值刚性显著增强,使得计算的收敛速度与稳定性进一步下降。提出了一种改进的线隐式BLU-SGS算法。该算法基于局部网格几何特征与流场物理特性并行构造隐式线,通过隐式线内控制体的联立求解,有效缓解了附面层法向小网格尺度对时间步长的限制,并提升了远场区域上下游流场信息的交互效率,实现计算收敛速度的显著提升。同时,沿隐式线推导了与离散方程右端项协调一致的雅可比矩阵,在前向与后向扫描过程中充分考虑了非对角线块的隐式贡献,通过保持隐式时间推进的部分非线性特性进一步增强了收敛效率。数值验证采用二维平板、三维M6机翼及CHN-T1飞机等典型算例,测试结果表明,新算法相较传统LU-SGS方法,收敛速度提升可达5至10倍。

关键词: 计算流体力学, 非结构混合网格, LU-SGS, 线隐算法, 收敛效率

Abstract: LU-SGS method is a mainstream time marching algorithm in computational fluid dynamics (CFD). It reduces computational costs by employing low-order Jacobian matrix approximations such as the maximum spectral radius. However, this approach sacrifices convergence efficiency. Especially in high Reynolds number flow simulations based on unstructured hybrid grid, the numerical stiffness is significantly exacerbated due to the wall-normal direction grid refinement in boundary layer regions, leading to further deterioration in convergence speed and stability. To address these issues, an improved line-implicit BLU-SGS algorithm is proposed. Based on local grid geometric characteristics and flow field physics, this algorithm constructs implicit lines in parallel. By simultaneously solving the coupled equations of control volumes within these implicit lines, the method effectively alleviates the time-step restrictions imposed by small grid scales in the wall-normal direction of boundary layers and enhances the efficiency of upstream-downstream information exchange in far-field regions, thereby achieving a significant acceleration of convergence. Additionally, the Jacobian matrix is derived along the implicit lines to ensure consistency with the discretized right-hand side terms of the governing equations. During forward and backward sweeps, the algorithm fully accounts for implicit contributions from off-diagonal blocks to preserve certain nonlinear characteristics of the implicit time marching, further improving convergence efficiency. Numerical validation using benchmark cases such as a 2D flat plate, a 3D M6 wing, and the CHN-T1 aircraft demonstrates that the new algorithm achieves 5 to 10 times faster convergence compared to the traditional LU-SGS method.

Key words: CFD, unstructured hybrid grids, LU-SGS, line-implicit, convergence efficiency

中图分类号: