Implicit method for the solution of supersonic and hypersonic 3D flow problems with Lower-Upper Symmetric-Gauss-Seidel preconditioner on multiple graphics processing units

被引:16
作者
Bocharov, A. N. [2 ]
Evstigneev, N. M. [1 ]
Petrovskiy, V. P. [2 ]
Ryabkov, O., I [1 ,2 ]
Teplyakov, I. O. [2 ]
机构
[1] Russian Acad Sci, Fed Res Ctr Comp Sci & Control, Pr 60 Letiya Oktyabrya,Bldg 9, Moscow 117312, Russia
[2] Russian Acad Sci, Joint Inst High Temp, Izhorskaya 13 Bldg 2, Moscow 125412, Russia
关键词
LU-SGS preconditioning; Implicit methods; Hypersonic flows; Computional fluid dynamics; Multiple GPUs; Krylov methods; SCHEME; ACCURATE; GRIDS;
D O I
10.1016/j.jcp.2019.109189
中图分类号
TP39 [计算机的应用];
学科分类号
081203 ; 0835 ;
摘要
The paper describes a numerical method for the solution of stationary gas dynamics 3D spatial equations on unstructured grids that is designed for multiple graphics processing unit (GPU) computational architecture. Discretization of governing equations is done using first and second order TVD schemes. The Newton's method with simple pseudo time-step homotopy is used to solve the problem. Each iteration step involves solution of the linear system originated from the linearization of gas dynamics equations. Krylov subspace iterative methods are used to solve the linear system. The main aim of the paper is to describe a preconditioning Lower-Upper Symmetric-Gauss-Seidel (LU-SGS) method and its adaptation on multiple GPU computational systems. It is shown that deliberately reordered matrices with rearranged solution process of arising lower and upper triangular linear systems allow one to obtain close algebraic properties to the original single threaded LUSGS. The method is benchmarked against published results. The analysis of computational efficiency and acceleration is presented for different flows with Mach number ranging from 1.2 up to 25. (C) 2019 Elsevier Inc. All rights reserved.
引用
收藏
页数:18
相关论文
共 33 条
[1]   Toward a GPU-aware comparison of explicit and implicit CFD simulations on structured meshes [J].
Aissa, Mohamed ;
Verstraete, Tom ;
Vuik, Cornelis .
COMPUTERS & MATHEMATICS WITH APPLICATIONS, 2017, 74 (01) :201-217
[2]  
[Anonymous], [No title captured]
[3]  
[Anonymous], 2007, MEC COMPUT
[4]  
[Anonymous], [No title captured]
[5]  
[Anonymous], [No title captured]
[6]  
[Anonymous], [No title captured]
[7]  
Berger AA, 2012, UNDERSTANDING AMERICAN ICONS: AN INTRODUCTION TO SEMIOTICS, P43
[8]   Coupled gas heat and mass transfer problem solution using fully implicit method on graphics processing units [J].
Bocharov, A. N. ;
Bityurin, V. A. ;
Evstigneev, N. M. ;
Fortov, V. E. ;
Golovin, N. N. ;
Petrovskiy, V. P. ;
Ryabkov, O. I. ;
Solomonov, Yu S. ;
Shustov, A. A. ;
Teplyakov, I. O. .
XXXIII INTERNATIONAL CONFERENCE ON EQUATIONS OF STATE FOR MATTER, 2019, 1147
[9]   Computational gas dynamics in a wide range of Mach number on heterogeneous cluster architecture [J].
Bocharov, A. N. ;
Evstigneev, N. M. ;
Ryabkov, O. I. .
XXX INTERNATIONAL CONFERENCE ON INTERACTION OF INTENSE ENERGY FLUXES WITH MATTER (ELBRUS 2015), 2015, 653
[10]   The correct use of the Lax-Friedrichs method [J].
Breuss, M .
ESAIM-MATHEMATICAL MODELLING AND NUMERICAL ANALYSIS-MODELISATION MATHEMATIQUE ET ANALYSE NUMERIQUE, 2004, 38 (03) :519-540