no. 5
Energy conservation and numerical stability for the reduced MHD models of the non-linear JOREK code
Franck, Emmanuel1 ; Hölzl, Matthias2 ; Lessig, Alexander2 ; Sonnendrücker, Eric2, 3 
1 Inria Nancy grand Est, TONUS Team, 67000 Strasbourg, France.
2 Max-Planck-Institut für Plasmaphysik, Boltzmannstr. 2, 85748 Garching, Germany.
3 Technische Universität München, Boltzmannstr. 3, 85748 Garching, Germany.
ESAIM: Mathematical Modelling and Numerical Analysis , Tome 49 (2015) no. 5, pp. 1331-1365

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.

Reçu le :
DOI : 10.1051/m2an/2015014
Classification : 65H10, 35Q35, 35Q60, 76W05
Keywords: MHD, instabilities, nonlinear solvers, reduction, toroidal

Franck, Emmanuel 1 ; Hölzl, Matthias 2 ; Lessig, Alexander 2 ; Sonnendrücker, Eric 2, 3

1 Inria Nancy grand Est, TONUS Team, 67000 Strasbourg, France.
2 Max-Planck-Institut für Plasmaphysik, Boltzmannstr. 2, 85748 Garching, Germany.
3 Technische Universität München, Boltzmannstr. 3, 85748 Garching, Germany. 
@article{M2AN_2015__49_5_1331_0,
     author = {Franck, Emmanuel and H\"olzl, Matthias and Lessig, Alexander and Sonnendr\"ucker, Eric},
     title = {Energy conservation and numerical stability for the reduced {MHD} models of the non-linear {JOREK} code},
     journal = {ESAIM: Mathematical Modelling and Numerical Analysis },
     pages = {1331--1365},
     year = {2015},
     publisher = {EDP Sciences},
     volume = {49},
     number = {5},
     doi = {10.1051/m2an/2015014},
     zbl = {1329.76382},
     language = {en},
     url = {https://www.numdam.org/articles/10.1051/m2an/2015014/}
}
TY  - JOUR
AU  - Franck, Emmanuel
AU  - Hölzl, Matthias
AU  - Lessig, Alexander
AU  - Sonnendrücker, Eric
TI  - Energy conservation and numerical stability for the reduced MHD models of the non-linear JOREK code
JO  - ESAIM: Mathematical Modelling and Numerical Analysis 
PY  - 2015
SP  - 1331
EP  - 1365
VL  - 49
IS  - 5
PB  - EDP Sciences
UR  - https://www.numdam.org/articles/10.1051/m2an/2015014/
DO  - 10.1051/m2an/2015014
LA  - en
ID  - M2AN_2015__49_5_1331_0
ER  - 
%0 Journal Article
%A Franck, Emmanuel
%A Hölzl, Matthias
%A Lessig, Alexander
%A Sonnendrücker, Eric
%T Energy conservation and numerical stability for the reduced MHD models of the non-linear JOREK code
%J ESAIM: Mathematical Modelling and Numerical Analysis 
%D 2015
%P 1331-1365
%V 49
%N 5
%I EDP Sciences
%U https://www.numdam.org/articles/10.1051/m2an/2015014/
%R 10.1051/m2an/2015014
%G en
%F M2AN_2015__49_5_1331_0
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. ESAIM: Mathematical Modelling and Numerical Analysis , Tome 49 (2015) no. 5, pp. 1331-1365. doi: 10.1051/m2an/2015014

M. Bécoulet, F. Orain, G.T.A. Huijsmans, S. Pamela, P. Cahyna, M. Hoelzl, X. Garbet, E. Franck, E. Sonnendrücker, G. Dif-Pradalier, C. Passeron, G. Latu, J. Morales, E. Nardon, A. Fil, B. Nkonga, A. Ratnani and V. Grandgirard, Mechanism of Edge Localized Mode mitigation by Resonant Magnetic Perturbations. Phys. Rev. Lett. 113 (2014) 115001. | DOI

P. Cahyna, M. Becoulet, G.T.A. Huijsmans, F. Orain, J. Morales, A. Kirk, A.J. Thornton, S. Pamela, R. Panek and M. Hoelzl, Modelling of spatial structure of divertor footprints caused by edge-localized modes mitigated by magnetic perturbations. Nucl. Fus. (submitted).

L. Chacón, Scalable parallel implicit solver for 3D magnetohydrodynamics. J. Phys.: Conf. Ser. 125 (2008) 012041.

L. Chacón, An optimal, parallel, fully implicit Newton-Krylov solver for three-dimensional viscoresistive magnetohydrodynamics. Phys. Plasmas 15 (2008) 056103. | DOI

L. Chacón and D.A. Knoll, A 2D high β hall MHD implicit nonlinear solver. J. Comput. Phys. 188 (2003) 573–592. | Zbl | DOI

L. Chacón, D.A. Knoll and J.M. Finn, An implicit, nonlinear reduced resistive MHD solver. J. Comput. Phys. 178 (2002) 15–36. | Zbl | DOI

R.S. Dembo, S.C. Eisenstat and T. Steihaug, Inexact Newton methods. SIAM J. Numer. Anal. 19 (1982) 400-408. | MR | Zbl | DOI

E. Deriaz, B. Després, G. Faccanoni, K.P. Gostaf, L.M. Imbert, G. Sadaka and R. Sart, Magnetic equations with FreeFem ++: the Grad−Shafranov equation and the Current Hole. ESAIM: Proc. 32 (2011) 149–162. | Zbl | MR

B. Després and R. Sart, Reduced resistive MHD in Tokamaks with general density. ESAIM: M2AN 46 (2012) 1081–1106. | MR | Zbl | Numdam | DOI

B. Després and R. Sart, Derivation of hierarchies of reduced MHD models in Tokamaks, submitted to Archive of Rational Mechanics and Analysis. Preprint: http://hal.archives-ouvertes.fr/docs/00/79/64/25/PDF/mhdtout12.pdf

S.C. Eisenstat and H.F. Walker, Globally convergent Inexact Newton methods. SIAM J. Sci. Stat. Comput. 6 (1985) 793–832. | MR | Zbl

A. Fil, E. Nardon, M. Bécoulet, G. Dif-Pradalier, V. Grandgirard, R. Guirlet, M. Hoelzl, G.T.A. Huijsmans, G. Latu, M. Lehnen, P. Monier-Garbet, F. Orain, C. Passeron, B. Pǵourié, C. Reux, F. Saint-Laurent and P. Tamain, Modeling of disruption mitigation by massive gas injection. 41st EPS Conference on Plasma Physics. Berlin, Germany (2014), P1.045.

R. Freund, G.H. Golub and N. Nachtigal, Iterative solution of linear systems. Acta Numerica (1992) 57–100. | MR | Zbl

M. Hoelzl, S. Günter, R.P. Wenninger, W.-C. Mueller, G.T.A. Huysmans, K. Lackner and I. Krebs, ASDEX Upgrade Team. Reduced-MHD Simulations of Toroidally and Poloidally Localized ELMs. Phys. Plasmas. 19 (2012) 082505. | DOI

G.T.A. Huysmans and O. Czarny, MHD stability in X-point geometry: simulation of ELMs. Nucl. Fusion 47 (2007) 659–666. | DOI

M. Hoelzl, G.T.A. Huijsmans, P. Merkel, C. Atanasiu, K. Lackner, E. Nardon, K. Aleynikova, F. Liu, E. Strumberger, R. Mcadams, I. Chapman and A. Fil, Non-Linear Simulations of MHD Instabilities in Tokamaks Including Eddy Current Effects and Perspectives for the Extension to Halo Currents. J. Phys.: Conf. Ser. 561 (2014) 012011.

G.T.A. Huijsmans, F. Liu, A. Loarte, S. Futatani, F. Koechl, M. Hoelzl, A. Garofalo,W. Salomon, P.B. Snyder, E. Nardon, F. Orain and M. Bécoulet, Non-linear MHD Simulations for ITER. 25th Fusion Energy Conference (FEC 2014). Saint Petersburg. Russia (2014) TH/6-1Ra.

G.T.A. Huysmans and O. Czarny, Bézier surfaces and finite elements for MHD simulations. J. Comput. Phys. Arch. 227 (2008) 7423–7445. | MR | Zbl | DOI

I. Krebs, M. Hoelzl, K. Lackner, S. Günter, Nonlinear excitation of low-n harmonics in reduced MHD simulations of edge-localized modes. Phys. Plasmas 20 (2013) 082506. | DOI

S.E. Kruger, C.C. Hegna and J.D. Callen, Generalized reduced magnetohydrodynamic equations. Phys. Plasmas 5 (1998). | MR | DOI

F. Liu, G.T.A. Huijsmans, A. Loarte, A.M. Garofalo, W.M. Solomon and M. Hoelzl, Nonlinear MHD simulations of QH-mode plasmas in DIII-D. 41st EPS Conference on Plasma Physics. Berlin, Germany (2014) O5.135.

S.K. Malapaka, B. Després and R. Sart, Unconditionally stable numerical simulations of a new generalized reduced resistive magnetohydrodynamics model. Int. J. Numer. Methods in Fluids 74 (2014) 231–249. | MR | Zbl | DOI

M. Martin, Modélisations fluides pour les plasmas de fusion: approximation par éléments finis C1 de Bell. Ph.D. thesis, University of Nice (2013).

F. Murphy, G.H. Golub and A.J. Wathen, A note on preconditioning for indefinite linear systems. Report (1999). | MR | Zbl

F. Orain, M. Bécoulet, J. Morales, G. Dif-Pradalier, X. Garbet, E. Nardon, C. Passeron, G. Latu, A. Fil, G.T.A. Huijsmans, M. Hoelzl, S. Pamela and P. Cahyna, Non-linear MHD modeling of multi-ELM cycles and mitigation by RMPs. Plasma Phys. Control. Fusion 57 (2014) 014020. | DOI

S.J.P. Pamela, G.T.A. Huysmans, A. Kirk, I.T. Chapman, M. Becoulet, F. Orain and M. Hoelzl and the MAST Team, Influence of Diamagnetic Effects on Resistive MHD Simulations of RMPs in MAST (In preparation).

B. Phillip, L. Chacón and M. Pernice, Implicit adaptive mesh refinement for 2D reduced resistive Magnetohydrodynamics. J. Comput. Phys. 227 (2008). | MR | Zbl

D.D. Schnack, D.C. Barnes, D.P. Brennan, C.C. Hegna, E. Held, C.C. Kim, S.E. Kruger, A.Y. Pankin and C.R. Sovinec, Computational modeling of fully ionized magnetized plasmas using the fluid approxmation. Phys. Plasmas 13 (2006) 058103. | MR | DOI

H. R. Strauss, Reduced MHD in nearly potential magnetic fields. J. Plasma Phys. 57 (1997) 83–87. | DOI

P.B. Snyder, H.R. Wilson, J.R. Ferron, L.L. Lao, A.W. Leonard, T.H. Osborne, A.D. Turnbull, D. Mossessian, M. Murakami and X.Q. Xu, Edge localized modes and the pedestal: A model based on coupled peeling ballooning modes. Phys. Plasmas 9 (2002) 2037–2043. | DOI

H. Zohm, Edge localized modes (ELMs). Plasma Phys. Control. Fusion 38 (1996) 105. | DOI

Cité par Sources :