Simulating the GEM challenge problem using a SUPG finite element method for the two fluid plasma model.

In this paper we present the results of simulations of the GEM challenge problem performed using the SUPG finite element method for the two fluid plasma model. Simulations performed with and without the use of discontinuity capturing terms in the fluid dynamics equations are presented, as well as simulations in which the term is switched based on a shock sensor function to obtain higher order accuracy away from discontinuities. Results of simulations performed with the discontinuity terms included in the electrodynamics equations are also presented. The use of discontinuity capturing terms in the fluid dynamics equations, at minimum, appears necessary to obtain a full duration simulation of this problem, however, inclusion of these terms was found to substantially increase the mesh resolution necessary to obtain a solution, and to induce nonphysical energy production. The rate of this nonphysical energy production was found to decrease with mesh refinement.