Title of article :
p = 2 Continuous finite elements on tetrahedra with local mass matrix inversion
Author/Authors :
Brazell، نويسنده , , Michael J. and Helenbrook، نويسنده , , Brian T.، نويسنده ,
Issue Information :
روزنامه با شماره پیاپی سال 2012
Pages :
10
From page :
289
To page :
298
Abstract :
When using continuous finite element methods for unsteady simulations, one must invert a global mass matrix even for an explicit time advancement scheme. In this paper we develop a local mass matrix inversion procedure for p = 2, C0 tetrahedral elements and also show that, unlike on triangles [B.T. Helenbrook, SIAM J. Numer. Anal., 47 (2009) 1304–1318], this procedure does not generalize to higher order. For the p = 2 basis, it is demonstrated that the accuracy of the local inversion is better than any previous mass lumping scheme. This allows the development of fast and low memory iterative algorithms as well as accurate explicit computations to be performed. The efficiency of using the approach as an iterative relaxation scheme for solving the Helmholtz equation is shown by comparing to more standard approaches such as p-multigrid with Jacobi relaxation. Lastly, the local mass matrix inversion is used to march the unsteady heat equation explicitly with third order spatial accuracy.
Keywords :
tetrahedra , hp-finite element method , Polynomial basis , Mass lumping , Mass matrix
Journal title :
Computer Methods in Applied Mechanics and Engineering
Serial Year :
2012
Journal title :
Computer Methods in Applied Mechanics and Engineering
Record number :
1595264
Link To Document :
بازگشت