Energy conservation and numerical stability for the reduced MHD models of the non-linear JOREK code. In this paper we present a rigorous derivation of the reduced MHD models with and without parallel velocity that are implemented in the non-linear MHD code JOREK. The model we obtain contains some terms that have been neglected in the implementation but might be relevant in the non-linear phase. These are necessary to guarantee exact conservation with respect to the full MHD energy. For the second part of this work, we have replaced the linearized time stepping of JOREK by a non-linear solver based on the Inexact Newton method including adaptive time stepping. We demonstrate that this approach is more robust especially with respect to numerical errors in the saturation phase of an instability and allows to use larger time steps in the non-linear phase.
Keywords for this software
References in zbMATH (referenced in 5 articles , 1 standard article )
Showing results 1 to 5 of 5.
- Haverkort, J.W.; de Blank, H.J.; Huysmans, G.T.A.; Pratt, J.; Koren, B.: Implementation of the full viscoresistive magnetohydrodynamic equations in a nonlinear finite element code (2016)
- Minjeaud, Sebastian; Pasquetti, Richard: Fourier-spectral element approximation of the ion-electron Braginskii system with application to tokamak edge plasma in divertor configuration (2016)
- Franck, Emmanuel; Hölzl, Matthias; Lessig, Alexander; Sonnendrücker, Eric: Energy conservation and numerical stability for the reduced MHD models of the non-linear JOREK code (2015)
- Bonnement, A.; Minjeaud, S.; Pasquetti, R.: Towards a high order Fourier-SEM solver of fluid models in tokamaks (2014)
- Després, Bruno; Sart, Rémy: Reduced resistive MHD in tokamaks with general density (2012)