A parallel implementation of an implicit discontinuous Galerkin finite element scheme for fluid flow problems

Abstract The discontinuous Galerkin (DG) method is frequently used in computational fluid dynamics for its stability and high order of accuracy. A disadvantage of the DG method is its high computational demands. The aim of this paper is to weaken this drawback by means of parallelization of the DG algorithm. The computation is performed on a network of computers with distributed memory using the Java Remote Method Invocation, which is included in the Java programming language. The partition of the boundary value problem into n subproblems, which is then solved by n computers separately, is based on the overlapping Schwarz method. On basis of physical nature of the problem, the present paper proposes minimal size of the overlap that allows for only one Schwarz iteration thereby increasing efficiency of parallelization. The scalability and efficiency of the presented parallelization approach is demonstrated on several test problems. In order to stabilize the DG method in presence of shocks, a recently developed technique by Huerta et al. (Int. J. Numer. Meth. Fluids 69(10), 2012, 1614–1632), which introduces discontinuities in basis functions in regions with a shock, is adopted here. A modification of this approach, which lowers the computational and implementational demands, is presented here.

[1]  J. Douglas,et al.  Interior Penalty Procedures for Elliptic and Parabolic Galerkin Methods , 1976 .

[2]  Chi-Wang Shu,et al.  TVB Runge-Kutta local projection discontinuous galerkin finite element method for conservation laws. II: General framework , 1989 .

[3]  Pierluigi Amodio,et al.  Numerical solution of differential algebraic equations and computation of consistent initial/boundary conditions , 1997 .

[4]  Douglas N. Arnold,et al.  Unified Analysis of Discontinuous Galerkin Methods for Elliptic Problems , 2001, SIAM J. Numer. Anal..

[5]  Chi-Wang Shu,et al.  The Runge-Kutta local projection discontinuous Galerkin finite element method for conservation laws. IV. The multidimensional case , 1990 .

[6]  Lilia Krivodonova,et al.  An efficient local time-stepping scheme for solution of nonlinear conservation laws , 2010, J. Comput. Phys..

[7]  D. Arnold.  An Interior Penalty Finite Element Method with Discontinuous Elements , 1982 .

[8]  W. H. Reed,et al.  Triangular mesh methods for the neutron transport equation , 1973 .

[9]  Chi-Wang Shu,et al.  Runge–Kutta Discontinuous Galerkin Methods for Convection-Dominated Problems , 2001, J. Sci. Comput..

[10]  Claus-Dieter Munz,et al.  Efficient Time Integration for Discontinuous Galerkin Methods for the Unsteady 3D Navier-Stokes Equations , 2012 .

[11]  Jan Vimmr,et al.  Comparison of discontinuous Galerkin time integration schemes for the solution of flow problems with deformable domains , 2015, Appl. Math. Comput..

[12]  O. Widlund,et al.  Some Domain Decomposition Algorithms for Elliptic Problems , 2011 .

[13]  Y. Saad,et al.  GMRES: a generalized minimal residual algorithm for solving nonsymmetric linear systems , 1986 .

[14]  J. Tinsley Oden,et al.  A discontinuous hp finite element method for the Euler and Navier–Stokes equations , 1999 .

[15]  Chi-Wang Shu,et al.  The Local Discontinuous Galerkin Method for Time-Dependent Convection-Diffusion Systems , 1998 .

[16]  S. Rebay,et al.  A High-Order Accurate Discontinuous Finite Element Method for the Numerical Solution of the Compressible Navier-Stokes Equations , 1997 .

[17]  J. Peraire,et al.  Sub-Cell Shock Capturing for Discontinuous Galerkin Methods , 2006 .

[18]  Sanjeev Sharma,et al.  AN OVERVIEW OF NON PORTABLE DISTRIBUTED TECHNIQUES. , 2013 .

[19]  J. Remacle,et al.  Shock detection and limiting with discontinuous Galerkin methods for hyperbolic conservation laws , 2004 .

[20]  Andreas Meister,et al.  A positivity preserving and well-balanced DG scheme using finite volume subcells in almost dry regions , 2016, Appl. Math. Comput..

[21]  O. Widlund,et al.  Additive Schwarz Methods for Elliptic Finite Element Problems in Three Dimensions , 2018 .

[22]  Jaime Peraire,et al.  Discontinuous Galerkin Solution of the Navier-Stokes Equations on Deformable Domains , 2007 .

[23]  Y. Maday,et al.  Optimized Schwarz methods without overlap for highly heterogeneous media , 2007 .

[24]  E. Toro.  Riemann Solvers and Numerical Methods for Fluid Dynamics , 1997 .

[25]  Antonio Huerta,et al.  A simple shock‐capturing technique for high‐order discontinuous Galerkin methods , 2012 .

[26]  Kalim Qureshi,et al.  A Performance Evaluation of RPC, Java RMI, MPI and PVM , 1970 .

[27]  Frédéric Magoulès,et al.  A Schur Complement Method for Compressible Navier-Stokes Equations , 2013, Domain Decomposition Methods in Science and Engineering XX.

[28]  Vít Dolejší,et al.  On the discontinuous Galerkin method for the numerical solution of the Navier–Stokes equations , 2004 .

[29]  J. Mark Bull,et al.  Benchmarking Java against C and Fortran for scientific applications , 2001, JGI '01.

[30]  Ewald Krämer,et al.  A parallel, high-order discontinuous Galerkin code for laminar and turbulent flows , 2008 .