A Fast, Matrix-free Implicit Method for Compressible Flows on Unstructured Grids

Abstract

A fast, matrix-free implicit method has been developed to solve the three-dimensional compressible Euler and Navier–Stokes equations on unstructured meshes. An approximate system of linear equations arising from the Newton linearization is solved by the GMRES (generalized minimum residual) algorithm with a LU-SGS (lower–upper symmetric Gauss–Seidel) preconditioner. A remarkable feature of the present GMRES+LU-SGS method is that the storage of the Jacobian matrix can be completely eliminated by approximating the Jacobian with numerical fluxes, resulting in a matrix-free implicit method. The method developed has been used to compute the compressible flows around 3D complex aerodynamic configurations for a wide range of flow conditions, from subsonic to supersonic. The numerical results obtained indicate that the use of the GMRES+LU-SGS method leads to a significant increase in performance over the best current implicit methods, GMRES+ILU and LU-SGS, while maintaining memory requirements similar to its explicit counterpart. An overall speedup factor from eight to more than one order of magnitude for all test cases in comparison with the explicit method is demonstrated.