Abstract

A new finite element (FE) scheme, based on the decomposition of a second order differential equation into a set of first order symplectic (Hamiltonian) equations, is presented and tested on a one-dimensional, driven SturmdashLiouville problem. Error analysis shows improved cubic convergence in the energy norm for piecewise linear “tent” elements, as compared to quadratic convergence for the standard and symplectic hybrid (i.e. ‘tent’ and piecewise constant) FE methods. The convergence deteriorates in the presence of a regular singular point, but can be recovered by appropriate mesh node packing. Optimal mesh packing exponents are derived to ensure cubic (respectively quadratic for the hybrid FE method) convergence with minimal numerical error. The symplectic hybrid FE scheme is shown to be approximately 30–40 times more accurate than the standard FE scheme, for an exact test problem based on determining the nonideal magnetohydrodynamic stability of a fusion plasma. A further suppression of the error by one order of magnitude is achieved for the symplectic tent element method.

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call

Disclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.