Implicit solution of harmonic balance equation system using the LU-SGS method and one-step Jacobi/Gauss-Seidel iteration

This paper presents efficient alternative numerical methods for an implicit solution of the harmonic balance equation system for analysing temporal periodic unsteady flows. The proposed method employs approximate factorisation to decouple the common residual term and the time spectral source term of a harmonic balance equation system when it is discretised implicitly. With this approximate factorisation, the complexity of implicit solution of the discrete system is greatly reduced. The common residual term can be dealt with using a lower-upper symmetric-Gauss-Seidel (LU-SGS) method and the time spectral source term is integrated using a Jacobi iteration (JI) or one step Gauss-Seidel (GS) iteration, leading to the LU-SGS/JI method or LU-SGS/GS method. The NASA stage 35 compressor and the 1.5 stage Aachen turbine were used to demonstrate the effectiveness of the proposed methods in stabilising solution and its advantages in comparison with the existing lower-upper symmetric-Gauss-Seidel/block Jacobi (LU-SGS/BJ) method. The LU-SGS/GS method and the LU-SGS/JI method are more robust than the LU-SGS/BJ method in stabilising solution. The LU-SGS/GS method also has faster and tighter convergence and lower memory consumption in comparison with the LU-SGS/BJ method.