Application of Mass Lumped Higher Order Finite Elements
Authors: J. Chen, H.R. Strauss, S.C. Jardin, W. Park, L.E. Sugiyama, G. Fu, and J. Breslau
Date of PPPL Report: November 2005
There are many interesting phenomena in extended-MHD such as anisotropic transport, mhd, 2-fluid effects stellarator and hot particles. Any one of them challenges numerical analysts, and researchers are seeking for higher order methods, such as higher order finite difference, higher order finite elements and hp/spectral elements. It is true that these methods give more accurate solution than their linear counterparts. However, numerically they are prohibitively expensive. Here we give a successful solution of this conflict by applying mass lumped higher order finite elements. This type of elements not only keep second/third order accuracy but also scale closely to linear elements by doing mass lumping. This is especially true for second order lump elements. Full M3D and anisotropic transport models are studied.