A parallel implementation of an implicit discontinuous Galerkin finite element scheme for fluid flow problems. (November 2017)

• Record Type: Journal Article
• Authors: Vimmr, Jan; Bublík, Ondřej; Pecka, Aleš
• Is Part Of: Advances in engineering software. Volume 113(2017)
• Page Start: 108
• Page End: 119
• Publication Date: 2017-11
• Subjects: Discontinuous Galerkin finite element method -- Implicit scheme -- Compressible Navier–Stokes equations -- Parallel computing -- Overlapping Schwarz method -- Shock capturing -- Java RMI; 76N10 -- 65N30 -- 65Y05 -- 35Q35; Computer-aided engineering -- Periodicals; Engineering -- Computer programs -- Periodicals; Engineering -- Software -- Periodicals; Periodicals; 620.0028553
• DOI: 10.1016/j.advengsoft.2016.11.007
• Languages: English
• ISSNs: 0965-9978

Abstract

Highlights: Interior penalty discontinuous Galerkin method for Navier–Stokes equations. Overlapping Schwarz method applied to parallel computing. Distributed computing on cluster/grid systems with Java Remote Method Invocation. Simple shock capturing technique based on sub-cell discontinuous basis functions. Benchmarking the developed parallel algorithm on heterogeneous network of computers. 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.