Application of Mass Lumped Higher Order Finite Elements

doi:10.2172/934516

Title: Application of Mass Lumped Higher Order Finite Elements

Full Record
Similar

Abstract

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.

Authors:: Chen, J.; Strauss, H. R.; Jardin, S. C.; Park, W.; Sugiyama, L. E.; G. Fu; Breslau, J.

Publication Date:: Tue Nov 01 00:00:00 EST 2005

Research Org.:: Princeton Plasma Physics Lab. (PPPL), Princeton, NJ (United States)

Sponsoring Org.:: USDOE Office of Science (SC)

OSTI Identifier:: 934516

Report Number(s):: PPPL-4128
TRN: US0804644

DOE Contract Number:: AC02-76CH03073

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.

Copy to clipboard


                    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.

Copy to clipboard


                    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.

Copy to clipboard


                    
@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 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.},

  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}

}

Copy to clipboard

Technical Report:

View Technical Report (0.65 MB)

DOI: 10.2172/934516

Save / Share:

Export Metadata

Save to My Library

Higher Order Lagrange Finite Elements In M3D

J. Chen ; H.R. Strauss ; S.C. Jardin ; ...

The M3D code has been using linear finite elements to represent multilevel MHD on 2-D 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 »« less
DOI: 10.2172/836490

Full Text Available
Visualizing higher order finite elements. Final report

Thompson, David C ; Pebay, Philippe Pierre

This report contains an algorithm for decomposing higher-order 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 »« less
DOI: 10.2172/876232

Full Text Available
Visualizing Higher Order Finite Elements: FY05 Yearly Report.

Thompson, David ; Pebay, Philippe Pierre

This report contains an algorithm for decomposing higher-order 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 mo-mentim. 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 free-dom, simulations will typically fail to produce a valid approximation to themore »« less
DOI: 10.2172/1143395

Full Text Available
Modeling 3D PCMI using the Extended Finite Element Method with higher order elements

Jiang, W. ; Spencer, Benjamin W.

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.
DOI: 10.2172/1409274

Full Text Available
Visualization of higher order finite elements.

Thompson, David C ; Pebay, Philippe Pierre ; Crawford, Richard H ; ...

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 axis-aligned regions of space, but when theremore »« less
DOI: 10.2172/919127

Full Text Available