Fluid Mechanics and Flight Mechanics

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

  • Yaobing ZHANG ,
  • Ming LI ,
  • Jian ZHANG ,
  • Naichun ZHOU ,
  • Xiaoquan GONG
Expand
  • 1.Computational Aerodynamics Institute,China Aerodynamics Research and Development Center,Mianyang 621000,China
    2.State Key Laboratory of Aerodynamics,Mianyang 621000,China
E-mail: zhangjian@cardc.cn

Received date: 2025-07-17

  Revised date: 2025-08-06

  Accepted date: 2025-08-15

  Online published: 2025-08-28

Supported by

National Natural Science Foundation of China(12402347);National Numerical Windtunnel (NNW) Project

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.

Cite this article

Yaobing ZHANG , Ming LI , Jian ZHANG , Naichun ZHOU , Xiaoquan GONG . An improved line-implicit BLU-SGS iteration algorithm for unstructured hybrid grids[J]. ACTA AERONAUTICAET ASTRONAUTICA SINICA, 2026 , 47(6) : 132588 -132588 . DOI: 10.7527/S1000-6893.2025.32588

References

[1] SLOTNICK J, KHODADOUST A, ALONSO J, et al. CFD Vision 2030 Study: A path to revolutionary computational aerosciences: NASA/CR-2014-218178[R]. Washington, D.C.: NASA. 2014.
[2] WANG G, JIANG Y W, YE Z Y. An improved LU-SGS implicit scheme for high Reynolds number flow computations on hybrid unstructured mesh[J]. Chinese Journal of Aeronautics201225(1): 33-41.
[3] JAMESON A, YOON S. Lower-upper implicit schemes with multiple grids for the Euler equations[J]. AIAA Journal198725(7): 929-935.
[4] MAVRIPLIS D J. On convergence acceleration techniques for unstructured meshes[C]∥29th AIAA, Fluid Dynamics Conference. Reston: AIAA, 1998.
[5] 洪俊武, 李伟, 岳皓, 等. 超大规模结构网格CFD数值模拟初步研究[J]. 航空学报202445(20): 129866.
  HONG J W, LI W, YUE H, et al. A pilot study on CFD numerical simulation of ultra-large-scale structured grid[J]. Acta Aeronautica et Astronautica Sinica202445(20): 129866 (in Chinese).
[6] 张健, 李瑞田, 邓亮, 等. 面向多核CPU/众核GPU架构的非结构CFD共享内存并行计算技术[J]. 航空学报202445(7): 128888.
  ZHANG J, LI R T, DENG L, et al. Shared-memory parallelization technology of unstructured CFD solver for multi-core CPU/many-core GPU architecture[J]. Acta Aeronautica et Astronautica Sinica202445(7): 128888 (in Chinese).
[7] MAVRIPLIS D J. Multigrid strategies for viscous flow solvers on anisotropic unstructured meshes[J]. Journal of Computational Physics1998145(1): 141-165.
[8] 曹文博, 刘溢浪, 张伟伟. 基于降阶模型和梯度优化的流场加速收敛方法[J]. 航空学报202344(6): 127090.
  CAO W B, LIU Y L, ZHANG W W. Accelerated convergence method for fluid dynamics solvers based on reduced-order model and gradient optimization[J]. Acta Aeronautica et Astronautica Sinica202344(6): 127090 (in Chinese).
[9] 胡姝瑶, 蒋崇文, 李椿萱. 可压缩流动模拟的多重网格-扰动域推进方法[J]. 航空学报202344(11): 127649.
  HU S Y, JIANG C W, LEE C H. Disturbance region update method with multigrid for compressible flows[J]. Acta Aeronautica et Astronautica Sinica202344(11): 127649 (in Chinese).
[10] 陈龙, 夏健, 田书玲. 适用于非定常流模拟的分布式并行GMRES方法[J]. 空气动力学学报201937(5): 834-843.
  CHEN L, XIA J, TIAN S L. A GMRES method on distributed parallel computers for unsteady flow simulation[J]. Acta Aerodynamica Sinica201937(5): 834-843 (in Chinese).
[11] ZHANG J, DENG L, LI R T, et al. Achieving high performance and portable parallel GMRES algorithm for compressible flow simulations on unstructured grids[J]. The Journal of Supercomputing202379(17): 20116-20140.
[12] PANDYA M J, DISKIN B, THOMAS J L, et al. Improved convergence and robustness of USM3D solutions on mixed-element grids[J]. AIAA Journal201654(9): 2589-2610.
[13] PANDYA M J, DISKIN B, THOMAS J L, et al. Assessment of USM3D hierarchical adaptive nonlinear method preconditioners for three-dimensional cases[J]. AIAA Journal201755(10): 3409-3424.
[14] PANDYA M J, JESPERSEN D C, DISKIN B, et al. Efficiency of mixed-element USM3D for benchmark three-dimensional flows[J]. AIAA Journal202159(8): 2997-3011.
[15] KALKOTE N, ASSAM A, ESWARAN V. Acceleration of later convergence in a density-based solver using adaptive time stepping[J]. AIAA Journal201957(1): 352-364.
[16] 阎超. 航空CFD四十年的成就与困境[J]. 航空学报202243(10): 526490.
  YAN C. Achievements and predicaments of CFD in aeronautics in past forty years[J]. Acta Aeronautica et Astronautica Sinica202243(10): 526490 (in Chinese).
[17] 黄宇, 阎超, 袁武. 适用于混合网格的改进雅可比迭代法及其应用[J]. 北京航空航天大学学报201642(3): 551-561.
  HUANG Y, YAN C, YUAN W. Improved Jacobi iterative method for hybrid grid and its application[J]. Journal of Beijing University of Aeronautics and Astronautics201642(3): 551-561 (in Chinese).
[18] LANGER S. Application of a line implicit method to fully coupled system of equations for turbulent flow problems[J]. International Journal of Computational Fluid Dynamics201327(3): 131-150.
[19] OTERO E, ELIASSON P. Acceleration on stretched meshes with line-implicit LU-SGS in parallel implementation[J]. International Journal of Computational Fluid Dynamics201529(2): 133-149.
[20] REMPKE A. Parallel line identification for line-implicit-solvers[J]. BIT Numerical Mathematics202363(4): 44.
[21] CHEN R F, WANG Z J. Fast, block lower-upper symmetric Gauss-seidel scheme for arbitrary grids[J]. AIAA Journal200038(12): 2238-2245.
[22] LIU J, CHEN J Q, YANG Y F, et al. Improvements for matrix-type implicit temporal scheme regarding entropy fix and local time step[J]. Chinese Journal of Aeronautics202437(11): 138-146.
[23] NIVED M R, REDDY P P K, ESWARAN V. Time-accurate solution of unsteady flows in an implicit solver using block LUSGS method[J]. International Journal of Computational Fluid Dynamics202337(3): 218-233.
[24] 陈坚强, 吴晓军, 张健, 等. FlowStar: 国家数值风洞(NNW)工程非结构通用CFD软件[J]. 航空学报202142(9): 625739.
  CHEN J Q, WU X J, ZHANG J, et al. FlowStar: General unstructured-grid CFD software for National Numerical Windtunnel (NNW) project[J]. Acta Aeronautica et Astronautica Sinica202142(9): 625739 (in Chinese).
[25] SHAROV D, NAKAHASHI K, SHAROV D, et al. Reordering of 3-D hybrid unstructured grids for vectorized LU-SGS Navier-Stokes computations[C]∥13th Computational Fluid Dynamics Conference. Reston: AIAA, 1997: 2102.
[26] 冯选燕, 燕振国, 朱华君, 等. 非精确Newton方法中线性迭代收敛判据研究[J]. 空气动力学学报202341(12): 28-36.
  FENG X Y, YAN Z G, ZHU H J, et al. Study on the convergence criterion of linear iteration in inexact Newton methods[J]. Acta Aerodynamica Sinica202341(12): 28-36 (in Chinese).
[27] 德)吉瑞·布拉泽克. 计算流体力学原理和应用[M]. 张耀冰, 陈江涛, 陈喜兰, 等, 译. 3版. 北京: 国防工业出版社, 2023: 92-95.
  JIRI. Computational fluid dynamics: Principles and applications[M]. ZHANG Y B, CHEN J T, CHEN X L, et al, trans. Beijing: National Defense Industry Press, 2023: 92-95 (in Chinese).
[28] NASA Langley Research Center. VERIF/2DZP: 2D zero pressure gradient flat plate verification case—Intro page[EB/OL]. (2021-03-24)[2025-06-10] .
[29] SAXENA S, NAIR M. Implementation and testing of Spalart-Allmaras model in a multi-block code[C]∥40th AIAA Aerospace Sciences Meeting & Exhibit. Reston: AIAA, 2002.
[30] SCHMITT V, CHARPIN F. Pressure distributions on the onera-m6 wing at transonic mach numbers, experimental data base for computer program assessment: AGARD AR138, 1979.
[31] WANG Q Z, STEVEN J M, KHALED S. Implementation of advanced two equation turbulence models in the USM3D unstructured flow solver: NASA CR-210202[R]. Washington, D.C.: NASA, 2000.
[32] 余永刚, 周铸, 黄江涛, 等. 单通道客机气动标模CHN-T1设计[J]. 空气动力学学报201836(3): 505-513.
  YU Y G, ZHOU Z, HUANG J T, et al. Aerodynamic design of a standard model CHN-T1for single-aisle passenger aircraft[J]. Acta Aerodynamica Sinica201836(3): 505-513 (in Chinese).
[33] 张耀冰, 唐静, 陈江涛, 等. 基于非结构混合网格的CHN-T1标模气动特性预测[J]. 空气动力学学报201937(2): 262-271.
  ZHANG Y B, TANG J, CHEN J T, et al. Aerodynamic characteristics prediction of CHN-T1 standard model with unstructured grid[J]. Acta Aerodynamica Sinica201937(2): 262-271 (in Chinese).
Outlines

/