航空学报 > 2026, Vol. 47 Issue (6): 132588-132588   doi: 10.7527/S1000-6893.2025.32588

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

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

  1. 1.中国空气动力研究与发展中心 计算空气动力研究所,绵阳 621000
    2.空天飞行空气动力科学与技术全国重点实验室,绵阳 621000
  • 收稿日期:2025-07-17 修回日期:2025-08-06 接受日期:2025-08-15 出版日期:2025-08-29 发布日期:2025-08-28
  • 通讯作者: 张健 E-mail:zhangjian@cardc.cn
  • 基金资助:
    国家自然科学基金(12402347);国家数值风洞(NNW)工程

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

Yaobing ZHANG1,2, Ming LI1, Jian ZHANG1(), Naichun ZHOU1,2, Xiaoquan GONG1   

  1. 1.Computational Aerodynamics Institute,China Aerodynamics Research and Development Center,Mianyang 621000,China
    2.State Key Laboratory of Aerodynamics,Mianyang 621000,China
  • Received:2025-07-17 Revised:2025-08-06 Accepted:2025-08-15 Online:2025-08-29 Published:2025-08-28
  • Contact: Jian ZHANG E-mail:zhangjian@cardc.cn
  • Supported by:
    National Natural Science Foundation of China(12402347);National Numerical Windtunnel (NNW) Project

摘要:

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 con-vergence. 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-10 times faster convergence compared to the traditional LU-SGS method.

Key words: computational fluid dynamics, unstructured hybrid grids, LU-SGS, line-implicit algorithm, convergence efficiency

中图分类号: