Application of Mass Lumped Higher Order Finite Elements
Abstract
There are many interesting phenomena in extendedMHD such as anisotropic transport, mhd, 2fluid 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.
 Authors:
 Publication Date:
 Research Org.:
 Princeton Plasma Physics Lab. (PPPL), Princeton, NJ (United States)
 Sponsoring Org.:
 USDOE Office of Science (SC)
 OSTI Identifier:
 934516
 Report Number(s):
 PPPL4128
TRN: US0804644
 DOE Contract Number:
 AC0276CH03073
 Resource Type:
 Technical Report
 Country of Publication:
 United States
 Language:
 English
 Subject:
 70 PLASMA PHYSICS AND FUSION TECHNOLOGY; ACCURACY; STELLARATORS; TRANSPORT; Higher order finite elements; mass lumping; highly anisotropic transport
Citation Formats
Chen, J., Strauss, H. R., Jardin, S. C., Park, W., Sugiyama, L. E., G. Fu, and Breslau, J. Application of Mass Lumped Higher Order Finite Elements. United States: N. p., 2005.
Web. doi:10.2172/934516.
Chen, J., Strauss, H. R., Jardin, S. C., Park, W., Sugiyama, L. E., G. Fu, & Breslau, J. Application of Mass Lumped Higher Order Finite Elements. United States. doi:10.2172/934516.
Chen, J., Strauss, H. R., Jardin, S. C., Park, W., Sugiyama, L. E., G. Fu, and Breslau, J. Tue .
"Application of Mass Lumped Higher Order Finite Elements". United States.
doi:10.2172/934516. https://www.osti.gov/servlets/purl/934516.
@article{osti_934516,
title = {Application of Mass Lumped Higher Order Finite Elements},
author = {Chen, J. and Strauss, H. R. and Jardin, S. C. and Park, W. and Sugiyama, L. E. and G. Fu and Breslau, J.},
abstractNote = {There are many interesting phenomena in extendedMHD such as anisotropic transport, mhd, 2fluid 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.},
doi = {10.2172/934516},
journal = {},
number = ,
volume = ,
place = {United States},
year = {Tue Nov 01 00:00:00 EST 2005},
month = {Tue Nov 01 00:00:00 EST 2005}
}

The M3D code has been using linear finite elements to represent multilevel MHD on 2D poloidal planes. Triangular higher order elements, up to third order, are constructed here in order to provide M3D the capability to solve highly anisotropic transport problems. It is found that higher order elements are essential to resolve the thin transition layer characteristic of the anisotropic transport equation, particularly when the strong anisotropic direction is not aligned with one of the Cartesian coordinates. The transition layer is measured by the profile width, which is zero for infinite anisotropy. It is shown that only higher order schemesmore »

Visualizing higher order finite elements. Final report
This report contains an algorithm for decomposing higherorder finite elements into regions appropriate for isosurfacing and proves the conditions under which the algorithm will terminate. Finite elements are used to create piecewise polynomial approximants to the solution of partial differential equations for which no analytical solution exists. These polynomials represent fields such as pressure, stress, and momentum. In the past, these polynomials have been linear in each parametric coordinate. Each polynomial coefficient must be uniquely determined by a simulation, and these coefficients are called degrees of freedom. When there are not enough degrees of freedom, simulations will typically fail tomore » 
Visualization of higher order finite elements.
Finite element meshes are used to approximate the solution to some differential equation when no exact solution exists. A finite element mesh consists of many small (but finite, not infinitesimal or differential) regions of space that partition the problem domain, {Omega}. Each region, or element, or cell has an associated polynomial map, {Phi}, that converts the coordinates of any point, x = ( x y z ), in the element into another value, f(x), that is an approximate solution to the differential equation, as in Figure 1(a). This representation works quite well for axisaligned regions of space, but when theremore » 
Visualizing Higher Order Finite Elements: FY05 Yearly Report.
This report contains an algorithm for decomposing higherorder finite elementsinto regions appropriate for isosurfacing and proves the conditions under which thealgorithm will terminate. Finite elements are used to create piecewise polynomialapproximants to the solution of partial differential equations for which no analyticalsolution exists. These polynomials represent fields such as pressure, stress, and momentim. In the past, these polynomials have been linear in each parametric coordinate.Each polynomial coefficient must be uniquely determined by a simulation, and thesecoefficients are called degrees of freedom. When there are not enough degrees of freedom, simulations will typically fail to produce a valid approximation to themore » 
Modeling 3D PCMI using the Extended Finite Element Method with higher order elements
This report documents the recent development to enable XFEM to work with higher order elements. It also demonstrates the application of higher order (quadratic) elements to both 2D and 3D models of PCMI problems, where discrete fractures in the fuel are represented using XFEM. The modeling results demonstrate the ability of the higher order XFEM to accurately capture the effects of a crack on the response in the vicinity of the intersecting surfaces of cracked fuel and cladding, as well as represent smooth responses in the regions away from the crack.