×

A coupled FEM-BEM approach for the solution of the free-boundary axi-symmetric plasma equilibrium problem. (English) Zbl 07493156

Summary: In this paper we present a coupled Finite Element Method - Boundary Element Method (FEM-BEM) approach for the solution of the free-boundary axi-symmetric plasma equilibrium problem. The proposed method, obtained from an improvement of the Hagenow-Lackner coupling method, allows to efficiently model the equilibrium problem in unbounded domains by discretizing only the plasma region; the external conductors can be modelled either as 2D or 3D models, according to the problem of interest. The paper explores different iterative methods for the solution of the nonlinear Grad-Shafranov equation, such as Picard, Newton-Raphson and Newton-Krylov, in order to provide a robust and reliable tool, able to handle large-scale problems (e.g. high resolution equilibria). This method has been implemented in the FRIDA code (FRee-boundary Integro-Differential Axisimmetric – https://github.om/matteobonotto/FRIDA), together with a suitable Adaptive Integration Technique (AIT) for the computation of the source term. FRIDA has been successfully tested and validated against experimental data from RFX-mod device, and numerical equilibria of an ITER-like device.

MSC:

82-XX Statistical mechanics, structure of matter
65-XX Numerical analysis
PDFBibTeX XMLCite
Full Text: DOI

References:

[1] J. Wesson and D. Campbell, Tokamaks. International Series of Monogr, OUP Oxford, 2011.
[2] J. P. Freidberg, Ideal MHD. Cambridge University Press, 2014.
[3] H. Grad and H. Rubin, “Hydromagnetic equilibria and force-free fields,” in Proc. Second UN Intern. Conf. on Peaceful Uses of Atomic Energy 31, p. 190-197, 1959.
[4] V. D. Shafranov, “Plasma Equilibrium in a Magnetic Field,” Reviews of Plasma Physics, vol. 2, p. 103, Jan. 1966.
[5] B. Faugeras and H. Heumann, “FEM-BEM coupling methods for tokamak plasma axisym-metric free-boundary equilibrium computations in unbounded domains,” Journal of Compu-tational Physics, vol. 343, pp. 201-216, 2017. · Zbl 1380.78003
[6] J. Blum, Numerical simulation and optimal control in plasma physics. United States: John Wiley and Sons Inc, 1989. · Zbl 0717.76009
[7] K. Lackner, “Computation of ideal MHD equilibria,” Computer Physics Communications, vol. 12, no. 1, pp. 33-44, 1976.
[8] F. Hofmann, “FBT -a free-boundary tokamak equilibrium code for highly elongated and shaped plasmas,” Computer Physics Communications, vol. 48, no. 2, pp. 207-221, 1988.
[9] P. Barabaschi, “The maxfea code,” Proceedings Plasma Control Technical Meeting, 1993.
[10] Y. M. Jeon, “Development of a free-boundary tokamak equilibrium solver for advanced study of tokamak equilibria,” Journal of the Korean Physical Society, vol. 67, p. 843-853, Sep 2015.
[11] R. Albanese, J. Blum, and O. D. Barbieri, “On the solution of the magnetic flux equation in an infinite domain,” EPS. 8th Europhysics Conference on Computing in Plasma Physics, 1986.
[12] R. Albanese, J. Blum, and O. D. Barbieri, “Numerical studies of the next european torus via the PROTEUS code,” 12th Conf. on Numerical Simulation of Plasmas, 1987.
[13] R. Albanese and F. Villone, “The linearized CREATE-l plasma response model for the control of current, position and shape in tokamaks,” Nuclear Fusion, vol. 38, pp. 723-738, may 1998.
[14] R. Albanese, R. Ambrosino, and M. Mattei, “CREATE-NL+: A robust control-oriented free boundary dynamic plasma equilibrium solver,” Fusion Engineering and Design, vol. 96-97, pp. 664-667, 2015. Proceedings of the 28th Symposium On Fusion Technology (SOFT-28).
[15] H. Heumann, J. Blum, C. Boulbe, B. Faugeras, G. Selig, J. Anè, S. Bremond, V. Grand-girard, P. Hertout, and E. Nardon, “Quasi-static free-boundary equilibrium of toroidal plasma with cedres++: Computational methods and applications,” Journal of Plasma Physics, vol. FirstView, pp. 1-35, 06 2015.
[16] J. L. Johnson, H. Dalhed, J. Greene, R. Grimm, Y. Hsieh, S. Jardin, J. Manickam, M. Ok-abayashi, R. Storer, A. Todd, et al., “Numerical determination of axisymmetric toroidal mag-netohydrodynamic equilibria,” Journal of Computational Physics, vol. 32, no. 2, pp. 212-234, 1979. · Zbl 0411.76098
[17] F. Villone, L. Barbato, S. Mastrostefano, and S. Ventre, “Coupling of nonlinear axisymmetric plasma evolution with three-dimensional volumetric conductors,” Plasma Physics and Con-trolled Fusion, vol. 55, p. 095008, Jul 2013.
[18] P. Piovesan, D. Bonfiglio, F. Auriemma, F. Bonomo, L. Carraro, R. Cavazzana, G. De Masi, A. Fassina, P. Franz, M. Gobbin, L. Marrelli, P. Martin, E. Martines, B. Momo, L. Piron, M. Valisa, M. Veranda, N. Vianello, B. Zaniol, M. Agostini, M. Baruzzo, T. Bolzonella, A. Canton, S. Cappello, L. Chacón, G. Ciaccio, D. F. Escande, P. Innocente, R. Lorenzini, R. Paccagnella, M. E. Puiatti, P. Scarin, A. Soppelsa, G. Spizzo, M. Spolaore, D. Terranova, P. Zanca, L. Zanotto, and M. Zuin, “RFX-mod: A multi-configuration fusion facility for three-dimensional physics studies,” Physics of Plasmas, vol. 20, no. 5, p. 056112, 2013.
[19] J. Luxon and B. Brown, “Magnetic analysis of non-circular cross-section tokamaks,” Nuclear Fusion, vol. 22, pp. 813-821, Jun 1982.
[20] C. Johnson and J. C. Nedelec, “On the coupling of boundary integral and finite element methods,” Mathematics of Computation, vol. 35, no. 152, pp. 1063-1079, 1980. · Zbl 0451.65083
[21] J. BIELAK and R. C. MacCAMY, “An exterior interface problem in two-dimensional elasto-dynamics,” Quarterly of Applied Mathematics, vol. 41, no. 1, pp. 143-159, 1983. · Zbl 0519.73021
[22] C. Geuzaine and J.-F. Remacle, “Gmsh: A 3-D finite element mesh generator with built-in pre-and post-processing facilities,” International Journal for Numerical Methods in Engineering, vol. 79, pp. 1309-1331, 09 2009. · Zbl 1176.74181
[23] D. Pepper and J. Heinrich, The Finite Element Method: Basic Concepts and Applications with MATLAB, MAPLE, and COMSOL, Third Edition. Computational and Physical Processes in Mechanics and Thermal Sciences, CRC Press, 2017. · Zbl 1360.65004
[24] D. A. Dunavant, “High degree efficient symmetrical Gaussian quadrature rules for the trian-gle,” International Journal for Numerical Methods in Engineering, vol. 21, no. 6, pp. 1129-1148, 1985. · Zbl 0589.65021
[25] A. Cenedese, P. Bettini, and M. Bonotto, “Model-based approach for magnetic reconstruc-tion in axisymmetric nuclear fusion machines,” IEEE Transactions on Plasma Science, vol. 46, pp. 636-644, March 2018.
[26] C. T. Kelley, Iterative Methods for Linear and Nonlinear Equations. No. 16, SIAM, 1995. · Zbl 0832.65046
[27] S. Banach, “Sur les opérations dans les ensembles abstraits et leur application aux équations intégrales,” Fundamenta Mathematicae, vol. 3, no. 1, pp. 133-181, 1922. · JFM 48.0201.01
[28] M. Crisfield, “Accelerating and damping the modified newton-raphson method,” Computers & Structures, vol. 18, no. 3, pp. 395-407, 1984. · Zbl 0526.65039
[29] W. Xu and T. F. Coleman, “Solving nonlinear equations with the Newton-Krylov method based on automatic differentiation,” Optimization Methods and Software, vol. 29, no. 1, pp. 88-101, 2014. · Zbl 1282.65052
[30] S. Jardin, Computational Methods in Plasma Physics. USA: CRC Press, Inc., 1st ed., 2010. · Zbl 1198.76002
[31] Y. Saad, Iterative Methods for Sparse Linear Systems. USA: Society for Industrial and Applied Mathematics, 2nd ed., 2003. · Zbl 1031.65046
[32] D. Knoll and D. Keyes, “Jacobian-free Newton-Krylov methods: a survey of approaches and applications,” Journal of Computational Physics, vol. 193, no. 2, pp. 357-397, 2004. · Zbl 1036.65045
[33] D. Knoll, V. Mousseau, L. Chacon, and J. Reisner, “Jacobian-free Newton-Krylov methods for the accurate time integration of stiff wave systems,” Journal of Scientific Computing, vol. 25, pp. 213-230, 11 2005. · Zbl 1203.65071
[34] L. Urankar, “Vector potential and magnetic field of current-carrying finite arc segment in analytical form, part iii: Exact computation for rectangular cross section,” IEEE Transactions on Magnetics, vol. 18, pp. 1860-1867, November 1982.
[35] D. Abate, G. Marchiori, and F. Villone, “Modelling and experimental validation of RFX-mod tokamak shaped discharges,” Fusion Engineering and Design, vol. 146, pp. 135-138, 2019. SI:SOFT-30.
[36] O. Kudlacek, P. Zanca, C. Finotti, G. Marchiori, R. Cavazzana, and L. Marrelli, “Real time measurement of plasma macroscopic parameters on RFX-mod using a limited set of sen-sors,” Physics of Plasmas, vol. 22, no. 10, p. 102503, 2015.
[37] J. H. Kane, D. E. Keyes, and K. G. Prasad, “Iterative solution techniques in boundary element analysis,” International Journal for Numerical Methods in Engineering, vol. 31, no. 8, pp. 1511-1536, 1991. · Zbl 0825.73901
[38] N. Jovanovic, D. Keyes, and G. Prasad, “Drop tolerance ilu preconditioners for iterative solution techniques in boundary element analysis,” in Boundary Elements XV (C. Brebbia and J. Rencis, eds.), pp. 501-516, Publ by Computational Mechanics Publ, 1993. Proceedings of the International Conference on Boundary Element Methods (BEM XV); Conference date: 10-08-1993 Through 13-08-1993.
[39] L. Chacón, “An optimal, parallel, fully implicit Newton-Krylov solver for three-dimensional viscoresistive magnetohydrodynamics,” Physics of Plasmas, vol. 15, no. 5, p. 056103, 2008.
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.