Skip to main content Accessibility help

Quasi-static free-boundary equilibrium of toroidal plasma with CEDRES++: Computational methods and applications

  • H. Heumann (a1), J. Blum (a1), C. Boulbe (a1), B. Faugeras (a1), G. Selig (a1), J.-M. Ané (a2), S. Brémond (a2), V. Grandgirard (a2), P. Hertout (a2) and E. Nardon (a2)...


We present a comprehensive survey of various computational methods in CEDRES++ (Couplage Equilibre Diffusion Résistive pour l'Etude des Scénarios) for finding equilibria of toroidal plasma. Our focus is on free-boundary plasma equilibria, where either poloidal field coil currents or the temporal evolution of voltages in poloidal field circuit systems are given data. Centered around a piecewise linear finite element representation of the poloidal flux map, our approach allows in large parts the use of established numerical schemes. The coupling of a finite element method and a boundary element method gives consistent numerical solutions for equilibrium problems in unbounded domains. We formulate a new Newton method for the discretized nonlinear problem to tackle the various nonlinearities, including the free plasma boundary. The Newton method guarantees fast convergence and is the main building block for the inverse equilibrium problems that we can handle in CEDRES++ as well. The inverse problems aim at finding either poloidal field coil currents that ensure a desired shape and position of the plasma or at finding the evolution of the voltages in the poloidal field circuit systems that ensure a prescribed evolution of the plasma shape and position. We provide equilibrium simulations for the tokamaks ITER and WEST to illustrate the performance of CEDRES++ and its application areas.


Corresponding author

éEmail address for correspondence:


Hide All
Albanese, R., Blum, J. and De Barbieri, O. 1986 On the solution of the magnetic flux equation in an infinite domain. In: EPS. 8th Europhysics Conf. on Computing in Plasma Physics, European Physical Society, Mulhouse, France, pp. 41–44.
Albanese, R., Blum, J. and De Barbieri, O. 1987 Numerical studies of the Next European Torus via the PROTEUS code. In: 12th Conf. on Numerical Simulation of Plasmas, San Francisco.
Albanese, R. and Villone, F. 1998 The linearized CREATE-L plasma response model for the control of current, position and shape in tokamaks. Nucl. Fusion 38 (5), 723.
Ané, J. M., Grandgirard, V., Albajar, F. and Johner, J. 2000 Design of next step tokamak: consistent analysis of plasma performance flux composition and poloidal field system. In: 18th IAEA Fusion Energy Conf., Sorrento,
Ariola, M. and Pironti, A. 2008 Magnetic Control of Tokamak Plasmas. London: Springer.
Artaud, J. al. 2010 The CRONOS suite of codes for integrated tokamak modelling. Nucl. Fusion 50 (4), 043 001.
Bielak, J. and MacCamy, R. C. 1991 Symmetric finite element and boundary integral coupling methods for fluid-solid interaction. Q. Appl. Math. 49 (1), 107119.
Blum, J. 1989 Numerical Simulation and Optimal Control in Plasma Physics: With Applications to Tokamaks. Paris: Wiley/Gauthier-Villars.
Blum, J., Boulbe, C. and Faugeras, B. 2012 Reconstruction of the equilibrium of the plasma in a tokamak and identification of the current density profile in real time. J. Comput. Phys. 231 (3), 960980.
Blum, J. and Heumann, H. 2014 Optimal control for quasi-static evolution of plasma equilibrium in tokamaks. Technical Report, INRIA, Sophia Antipolis, HAL Id: hal-00988045, version 1.
Blum, J., Le Foll, J. and Thooris, B. 1981 The self-consistent equilibrium and diffusion code SCED. Comput. Phys. Commun. 24, 235254.
, Bucalossi, et al. 2011 Feasibility study of an actively cooled tungsten divertor in tore supra for ITER technology testing. Fusion Eng. Des. 86 (6–8), 684688.
Chen, G. and Zhou, J. 1992 Boundary Element Methods (Computational Mathematics and Applications). London: Academic Press, Ltd.
Ciarlet, P. G. 1978 The Finite Element Method for Elliptic Problems. Amsterdam: North-Holland Publishing Co.
Cooper, W. A., Hirshman, S. P., Merkel, P., Graves, J. P., Kisslinger, J., Wobig, H. F. G., Narushima, Y., Okamura, S. and Watanabe, K. Y. 2009 Three-dimensional anisotropic pressure free boundary equilibria. Comput. Phys. Commun. 180 (9), 15241533.
Costabel, M. 1987 Principles of boundary element methods. In: Finite Elements in Physics (Lausanne, 1986), Amsterdam: North-Holland, pp. 243274.
Costabel, M. and Stephan, E. P. 1990 Coupling of finite and boundary element methods for an elastoplastic interface problem. SIAM J. Numer. Anal. 27 (5), 12121226.
Coster, D. P., Basiuk, V., Pereverzev, G., Kalupin, D., Zagórksi, R., Stankiewicz, R., Huynh, P. and Imbeaux, F. 2010 The european transport solver. IEEE Trans. Plasma Sci. 38 (9), 20852092.
Davis, T. A. 2011 Suitesparse: a suite of sparse matrix software.
Degtyarev, L. M. and Drozdov, V. V. 1985 An inverse variable technique in the MHD-equilibrium problem. Comput. Phys. Rep. 2 (7), 341387.
Degtyarev, L. M. and Drozdov, V. V. 1991 Adaptive mesh computation of magnetic hydrodynamic equilibrium. Int. J. Mod. Phys. C 02 (01), 3038.
Delfour, M. C. and Zolésio, J.-P. 2011 Shapes and Geometries, 2nd edn., Advances in Design and Control, Vol. 22. Philadelphia, PA: Society for Industrial and Applied Mathematics (SIAM).
DeLucia, J., Jardin, S. C. and Todd, A. M. M. 1980 An iterative metric method for solving the inverse tokamak equilibrium problem. J. Comput. Phys. 37 (2), 183204.
Falchetto, G. al. 2014 The European Integrated Tokamak Modelling (ITM) effort: achievements and first physics results. Nucl. Fusion 54 (4), 043 018.
Feistauer, M. and Sobotikova, V. 1990 Finite element approximation of nonlinear elliptic problems with discontinuous coefficients. ESAIM: Math. Modelling Numer. Anal. - Modlisation Mathématique et Analyse Numrique 24 (4), 457500.
Feneberg, W. and Lackner, K. 1973 Multipole tokamak equilibria. Nucl. Fusion 13 (4), 549.
Fitzgerald, M., Appel, L. C. and Hole, M. J. 2013 EFIT tokamak equilibria with toroidal flow and anisotropic pressure using the two-temperature guiding-centre plasma. Nucl. Fusion 53 (11), 113 040.
Freidberg, J. P. 1987 Ideal Magnetohydrodynamics. US: Plenum.
Gatica, G. N. and Hsiao, G. C. 1995 The uncoupling of boundary integral and finite element methods for nonlinear boundary value problems. J. Math. Anal. Appl. 189 (2), 442461.
Glowinski, R. and Marrocco, A. 1974 Analyse numérique du champ magnétique d'un alternateur par éléments finis et sur-relaxation ponctuelle non linéaire. Comput. Methods Appl. Mech. Eng. 3 (1), 5585.
Goedbloed, J. P., Keppens, R. and Poedts, S. 2010 Advanced Magnetohydrodynamics: with Applications to Laboratory and Astrophysical Plasmas. Cambridge: Cambridge University Press.
Goedbloed, J. P. and Lifschitz, A. 1997 Stationary symmetric magnetohydrodynamic flows. Phys. Plasmas (1994–present) 4 (10), 35443564.
Goedbloed, J. P. and Poedts, S. 2004 Principles of Magnetohydrodynamics: with Applications to Laboratory and Astrophysical Plasmas. Cambridge: Cambridge University Press.
Grad, H. 1967 Toroidal containment of a plasma. Phys. Fluids (1958–1988) 10 (1), 137154.
Grad, H. and Rubin, H. 1958 Hydromagnetic equilibria and force-free fields. In: Proc. 2nd UN Conf. on the Peaceful Uses of Atomic Energy, United Nations Publications, Geneva, Vol. 31, p. 190.
Grandgirard, V. 1999 Modélisation de l'équilibre d'un plasma de tokamak. PhD thesis, l'Université de Franche-Comté.
Gruber, R., Iacono, R. and Troyon, F. 1987 Computation of MHD equilibria by a quasi-inverse finite hybrid element approach. J. Comput. Phys. 73 (1), 168182.
Guazzotto, L., Betti, R., Manickam, J. and Kaye, S. 2004 Numerical study of tokamak equilibria with arbitrary flow. Phys. Plasmas (1994–present) 11 (2), 604614.
Helton, F. J. and Wang, T. S. 1978 MHD equilibrium in non-circular tokamaks with field-shaping coil systems. Nucl. Fusion 18 (11), 1523.
Hertout, P., Boulbe, C., Nardon, E., Blum, J., Bremond, S., Bucalossi, J., Faugeras, B., Grandgirard, V. and Moreau, P. 2011 The cedres++ equilibrium code and its application to ITER, JT-60SA and Tore Supra. Fusion Eng. Des. 86, 10451048.
Hinton, F. L. and Hazeltine, R. D. 1976 Theory of plasma transport in toroidal confinement systems. Rev. Mod. Phys. 48, 239308.
Hiptmair, R. 2003 Coupling of finite elements and boundary elements in electromagnetic scattering. SIAM J. Numer. Anal. 41 (3), 919944.
Hirshman, S. P. and Betancourt, O. 1991 Preconditioned descent algorithm for rapid calculations of magnetohydrodynamic equilibria. J. Comput. Phys. 96 (1), 99109.
Hirshman, S. P. and Jardin, S. C. 1979 Two-dimensional transport of tokamak plasmas. Phys. Fluids 22 (4), 731742.
Hofmann, F. and Tonetti, G. 1988 Tokamak equilibrium reconstruction using Faraday rotation measurements. Nucl. Fusion 28 (10), 1871.
Huysmans, G. T. A. and Czarny, O. 2007 MHD stability in X-point geometry: simulation of ELMs. Nucl. Fusion 47 (7), 659.
Huysmans, G. T. A., Goedbloed, J. P. and Kerner, W. 1991 Isoparametric bicubic Hermite elements for solution of the Grad–Shafranov equation. In: Proc. CP90 Conf. on Comp. Phys., Vol. 2, pp. 371–376.
ITM 2013 Integrated tokamak modelling., integrated Tokamak Modelling.
Jardin, S. C. 2010 Computational Methods in Plasma Physics. Boca Raton, Florida: CRC Press/Taylor and Francis.
Jardin, S. C., Pomphrey, N. and DeLucia, J. 1986 Dynamic modeling of transport and positional control of tokamaks. J. Comput. Phys. 66, 481507.
Johnson, J. al. 1979 Numerical determination of axisymmetric toroidal magnetohydrodynamic equilibria. J. Comput. Phys. 32 (2), 212234.
Lackner, K. 1976 Computation of ideal MHD equilibria. Comput. Phys. Commun. 12 (1), 3344.
Lao, L. L., Ferron, J. R., Geoebner, R. J., Howl, W., St. John, H. E., Strait, E. J. and Taylor, T. S. 1990 Equilibrium analysis of current profiles in Tokamaks. Nucl. Fusion 30 (6), 1035.
Lao, L. L., Hirshman, S. P. and Wieland, R. M. 1981 Variational moment solutions to the Grad–Shafranov equation. Phys. Fluids 24 (8), 14311440.
Lao, L. L., John, H. St., Stambaugh, R. D., Kellman, A. G. and Pfeiffer, W. 1985 Reconstruction of current profile parameters and plasma shapes in tokamaks. Nucl. Fusion 25 (11), 1611.
Li, J., Melenk, J. M., Wohlmuth, B. and Zou, J. 2010 Optimal a priori estimates for higher order finite elements for elliptic interface problems. Appl. Numer. Math. 60 (1–2), 1937.
Ling, K. M. and Jardin, S. C. 1985 The Princeton spectral equilibrium code: PSEC. J. Comput. Phys. 58 (3), 300335.
Lüst, R. and Schlüter, A. 1957 Axialsymmetrische magnetohydrodynamische Gleichgewichtskonfigurationen. Z. Naturforsch. A12, 850854.
Luxon, J. L. and Brown, B. B. 1982 Magnetic analysis of non-circular cross-section tokamaks. Nucl. Fusion 22 (6), 813.
Maschke, E. K. and Perrin, H. J. 1984 An analytic solution of the stationary MHD equations for a rotating toroidal plasma. Phys. Lett. A 102 (3), 106108.
Mc Carthy, P. J., Martin, P. and Schneider, W. 1999 The CLISTE interpretive equilibrium code. Technical Report IPP Report 5/85. Max-Planck-Institut fur Plasmaphysik.
Murat, F. and Simon, J. 1976 Sur le contrôle par un domaine géométrique. Technical Report 76015. Laboratoire d'Analyse Numérique, Université de Paris 6.
Nédélec, J.-C. 2001 Acoustic and Electromagnetic Equations Applied Mathematical Sciences, 144). New York: Springer-Verlag.
Nocedal, J. and Wright, S. J. 2006 Numerical Optimization, 2nd edn.Springer Series in Operations Research and Financial Engineering. New York: Springer.
Parail, al. 2013 Self-consistent simulation of plasma scenarios for iter using a combination of 1.5d transport codes and free-boundary equilibrium codes. Nucl. Fusion 53 (11), 113 002.
Park, W., Belova, E. V., Fu, G. Y., Tang, X. Z., Strauss, H. R. and Sugiyama, L. E. 1999 Plasma simulation studies using multilevel physics models. Phys. Plasmas (1994-present) 6 (5), 17961803.
Pechstein, C. and Jüttler, B. 2006 Monotonicity-preserving interproximation of B-H-curves. J. Comput. Appl. Math. 196 (1), 4557.
Pustovitov, V. D. 2010 Anisotropic pressure effects on plasma equilibrium in toroidal systems. Plasma Phys. Control. Fusion 52 (6), 065 001.
Renard, Y. and Pommier, J. 2014 GetFEM++, an open-source finite element library.
Schwab, C. 2004 p- and hp- Finite Element Methods: Theory and Applications in Solid and Fluid Mechanics. Oxford: Clarendon Press.
Shafranov, V. D. 1958 On magnetohydrodynamical equilibrium configurations. Sov. J. Exp. Theor. Phys. 6, 545.
Shewchuk, J. R. 1996 Triangle: Engineering a 2D quality mesh generator and delaunay triangulator. In: Applied Computational Geometry: Towards Geometric Engineering, Vol. 1148 (ed. Lin, M. C. and Manocha, D.), Lecture Notes in Computer Science, Springer-Verlag, from the First ACM Workshop on Applied Computational Geometry, Springer, Berlin, pp. 203222.
Stephan, E. P. 1992 Coupling of finite elements and boundary elements for some nonlinear interface problems. Comput. Methods Appl. Mech. Eng. 101 (1–3), 6172.
Turkington, B., Lifschitz, A., Eydeland, A. and Spruck, J. 1993 Multiconstrained variational problems in magnetohydrodynamics: equilibrium and slow evolution. J. Comput. Phys. 106 (2), 269285.
Wesson, J. and Campbell, D. J. 2004 Tokamaks. The International Series of Monographs in Physics, 3rd edn. Vol. 118, Oxford: Clarendon Press.
Zhao, K., Vouvakis, M. N. and Lee, J.-F. 2006 Solving electromagnetic problems using a novel symmetric fem-bem approach. IEEE Trans. Magn. 42 (4), 583586.
Zwingmann, W., Eriksson, L.-G. and Stubberfield, P. 2001 Equilibrium analysis of tokamak discharges with anisotropic pressure. Plasma Phys. Control. Fusion 43 (11), 1441.
MathJax is a JavaScript display engine for mathematics. For more information see

Quasi-static free-boundary equilibrium of toroidal plasma with CEDRES++: Computational methods and applications

  • H. Heumann (a1), J. Blum (a1), C. Boulbe (a1), B. Faugeras (a1), G. Selig (a1), J.-M. Ané (a2), S. Brémond (a2), V. Grandgirard (a2), P. Hertout (a2) and E. Nardon (a2)...


Full text views

Total number of HTML views: 0
Total number of PDF views: 0 *
Loading metrics...

Abstract views

Total abstract views: 0 *
Loading metrics...

* Views captured on Cambridge Core between <date>. This data will be updated every 24 hours.

Usage data cannot currently be displayed