Abstract

A high-order Galerkin Least-Squares (GLS) finite element discretization is combined with massively parallel implicit solvers. The Balancing Domain Decomposition by Constraints (BDDC) algorithm is applied to the linear system arising from the two-dimensional, high-order discretization of the advection-diffusion equation and the Euler equation. The Robin-Robin interface condition is extended to the Euler equations using entropy-symmetrized variables. The BDDC method maintains scalability for the high-order discretization of the diffusion-dominated flows, and achieves low iteration count in the advection-dominated regime. The BDDC method based on inexact local solvers with incomplete factorization and p = 1 coarse correction maintains the performance of the exact counterpart for the wide range of the Peclet numbers considered.

Full Text
Published version (Free)

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call