A multi-block viscous flow solver based on GPU parallel methodology

Abstract A multi-block viscous flow solver for steady and unsteady turbulent flows based on GPU parallel methodology under the finite volume frame is presented in this paper. Both numerical accuracy and computational efficiency are concerned. Numerical flux scheme for all speeds SLAU is adopted because of its wide adaptability and strong robustness; high order reconstruction schemes like MLP and WENO are chosen to evaluate the inviscid terms while a set of fully conservative 4th-order central differencing schemes are utilized to deal with the viscous terms. Second-order temporal accuracy is forfeited for unsteady simulations by coupling DP-LUR with dual time-stepping strategy. Furthermore, heterogeneous multiple CPU + GPU coprocessing system is well established with CUDA and MPI methodology. Design details about GPU implementation are analyzed and discussed. Impressive speedup factor is achieved on our GPU platform compared with CPU indicating the bright feature of these algorithms. Numerical results of several complex configurations have demonstrated the validity and reliability for aerospace engineering applications.

[1]  A. Jameson,et al.  Lower-upper Symmetric-Gauss-Seidel method for the Euler and Navier-Stokes equations , 1988 .

[2]  Yong-Tao Zhang,et al.  Resolution of high order WENO schemes for complicated flow structures , 2003 .

[3]  Meng-Sing Liou,et al.  A sequel to AUSM, Part II: AUSM+-up for all speeds , 2006, J. Comput. Phys..

[4]  P. Spalart,et al.  A New Version of Detached-eddy Simulation, Resistant to Ambiguous Grid Densities , 2006 .

[5]  Naga K. Govindaraju,et al.  A Survey of General‐Purpose Computation on Graphics Hardware , 2007 .

[6]  Graham Pullan,et al.  An Accelerated 3D Navier–Stokes Solver for Flows in Turbomachines , 2009 .

[7]  Carl Ollivier-Gooch,et al.  Effect of discretization order on preconditioning and convergence of a high-order unstructured Newton-GMRES solver for the Euler equations , 2008, J. Comput. Phys..

[8]  J. Dacles-Mariani,et al.  Numerical/experimental study of a wingtip vortex in the near field , 1995 .

[9]  O. J. Boelens,et al.  Analysis of Transonic Flow on a Slender Delta Wing Using CFD , 2006 .

[10]  F. Nicoud,et al.  Large-Eddy Simulation of the Shock/Turbulence Interaction , 1999 .

[11]  Kyu Hong Kim,et al.  Accurate, efficient and monotonic numerical methods for multi-dimensional compressible flows Part II: Multi-dimensional limiting process , 2005 .

[12]  E. Shima,et al.  On New Simple Low-Dissipation Scheme of AUSM-Family for All Speeds , 2009 .

[13]  Alistair Revell,et al.  DESider A European Effort on Hybrid RANS-LES Modelling: Results of the European-Union Funded Project, 2004 - 2007 , 2009 .

[14]  Shia-Hui Peng,et al.  A Comparison of Turbulence Models in Prediction of Flow Around the DLR-F6 Aircraft Configuration , 2004 .

[15]  Nobuyuki Satofuka,et al.  Calculation by artificial compressibility method and virtual flux method on GPU , 2011 .

[16]  Graham Pullan,et al.  Acceleration of a two-dimensional Euler flow solver using commodity graphics hardware , 2007 .

[17]  Gecheng Zha,et al.  Implicit WENO Scheme and High Order Viscous Formulas for Compressible Flows , 2007 .

[18]  Graham V. Candler,et al.  Data-parallel lower-upper relaxation method for the navier-stokes equations , 1996 .

[19]  Z. Wang.  High-order methods for the Euler and Navier–Stokes equations on unstructured grids , 2007 .

[20]  Russell R. Boyce,et al.  Numerical Investigation of Wall Temperature and Entropy Layer Effects on Double Wedge Shock / Boundary Layer Interactions , 2006 .

[21]  Yoonsik Kim,et al.  Drag Prediction of DLR-F6 Using the Turbulent Navier-Stokes Calculations with Multigrid , 2004 .

[22]  Florian R. Menter,et al.  The Scale-Adaptive Simulation Method for Unsteady Turbulent Flow Predictions. Part 1: Theory and Model Description , 2010 .

[23]  W. K. Anderson,et al.  Comparison of Finite Volume Flux Vector Splittings for the Euler Equations , 1985 .

[24]  Mark J. Harris.  Fast fluid dynamics simulation on the GPU , 2005, SIGGRAPH Courses.

[25]  P. Spalart.  A One-Equation Turbulence Model for Aerodynamic Flows , 1992 .

[26]  Hubert Nguyen,et al.  GPU Gems 3 , 2007 .

[27]  F. Menter.  Two-equation eddy-viscosity turbulence models for engineering applications , 1994 .