Fully implicit multiple graphics processing units’ schemes for hypersonic flows with lower upper symmetric Gauss-Seidel preconditioner on unstructured non-orthogonal grids

Abstract

The paper discusses the simulation of the hypersonic flow problems on unstructured grids using fully implicit scheme and the analysis of the orthogonal corrections on the quality of heat fluxes on the surfaces of the bodies subject to hypersonic flow. The governing equations are the Navier-Stokes equations for viscous calorically perfect gas. Second order finite volume discretization is applied to the problem on arbitrary unstructured grids with implicit temporal treatment. The problem is solved using the Newton-Raphson method. Such methods require the solution of large linear systems with iterative solvers and require a preconditioning operator to converge. In this paper a well known Lower-Upper Symmetric Gauss Seidel (LU-SGS) preconditioner is applied. The new method is based on the tricky reordering for the factored Jacobi matrix that allows one to execute block triangular matrix solvers on GPUs in parallel without the loss of algebraic properties of the original non-factored operator. Good convergence properties are demonstrated for large scale problems of external aerodynamics with implicit Courant number around 1000 - 10000 for flows with Mach number around 13-25. The influence of the non-orthogonal corretions on the fluxes for highly scewed anisotropic grids is demonstrated at Mach numbers around 20. It is shown that the non-orthogonality corrections are of the second infinitesimal order compared to the effects of the shock fix influence on the heat fluxes on the body.