Implicit SUPG solution of Euler equations using edge-based data structures Original Research Article

In this work we present an implicit, edge-based implementation of the semi-discrete SUPG formulation with shock-capturing for the Euler equations in conservative variables. By disassembling the resulting finite element matrices into their edge contributions, sparse matrix coefficients, residuals and matrix-vector products needed in Krylov-update techniques are computed based on edge data structures. The resulting solution method requires less memory and CPU time than element-based implementations.