Discontinuous Galerkin Method for MHD
Abstract
While the discontinuous Galerkin (DG) method has received considerable attention recently, there is considerable uncertainty about its applicability to MHD calculations especially for long integrations times, t >> L/v_A, due to the discontinuous approximation space. We investigate the DG method applied to a simple equilibrium in two dimensions with no flows and constant total pressure. Our hp implementation of the DG method uses unstructured hierarchal triangular elements in which each element can individually select a polynomial order from one to ten, providing for both grid refinement (h), and order refinement (p). We observe the behavior of numerically induced nablaot B as a function of time and the dependence on grid refinement and polynomial order. We also explore several ``divergence cleaning'' techniques and their dependence on polynomial order.
- Publication:
-
APS Division of Plasma Physics Meeting Abstracts
- Pub Date:
- October 2003
- Bibcode:
- 2003APS..DPPFP1110W