Bézier surfaces and finite elements for MHD simulations

O. Czarny, G. Huysmans

Research output: Contribution to journalArticleAcademicpeer-review

113 Citations (Scopus)
2 Downloads (Pure)


A finite element method based on bicubic Bézier surfaces is applied to the simulation of MHD instabilities relevant to magnetically confined fusion. The major advantage of the new technique is that it allows a natural way to implement mesh refinement strategy, which is not supported by a pure Hermite formulation. Compared to a Lagrangian formulation the number of degrees of freedom is significantly reduced. The use of an isoparametric representation of the space coordinates allows an accurate alignment of the finite elements to the magnetic field line geometry in a tokamak plasma. The Bézier finite elements have been implemented in a MHD code using the non-linear reduced MHD model in toroidal geometry. As an illustration, results for Soloviev equilibrium and time-dependent current-hole computations are presented and discussed.
Original languageEnglish
Pages (from-to)7423-7445
JournalJournal of Computational Physics
Issue number16
Publication statusPublished - 1 Aug 2008
Externally publishedYes


Dive into the research topics of 'Bézier surfaces and finite elements for MHD simulations'. Together they form a unique fingerprint.

Cite this