FLIPEC, an Ideal MHD free-boundary equilibrium solver for plasma with flows in axisymmetric geometries.

POSTER

Abstract

The ability of calculating ideal MHD equilibria is needed not only during the design of a magnetic fusion device, but also in its normal operation. In the case of tokamak axisymmetric configurations with significant plasma flow, the plasma edge is maintained fixed as the code searches for the equilibrium in the majority of codes. However, plasma flows may induce changes such as displacements of the position of the X-point, the magnetic axis or the shape of the plasma boundary, to name a few, that cannot be quantified properly in a fixed-boundary setup. In this contribution we present a new equilibrium code and its most significant and recent upgrades. In this context, FLIPECode was written to obtain free-boundary axisymmetric ideal MHD equilibria for arbitrary plasma flows. Particularly, the free-plasma-boundary scheme has been previously applied to 3D, non-axisymmetric configurations in the absence of flow within the SIESTA equilibrium code.

The upgrades open the code to many other possibilities. FLIPEC was originally written in toroidal coordinates, (r,θ,ξ), providing the code with a circular domain section, what made the code incompatible with many tokamaks, since the coils and solenoid leads to a singularity in the magnetic field when included in the domain. Currently, general coordinates have been implemented and the computational mesh can be adapted to any plasma shape. Additionally, the onset of a vertical instability has been successfully addressed by the use of a control stability scheme which utilizes two stabilization coils.

Examples of application of the code to ITER and NSTX configurations will be used to showcase the capabilities of the code. Additional developments currently in progress, such as the possibility of keeping fixed some geometrical locations (magnetic axis, x-point) and the targeting to significant parameters (integrated current) will be discussed as well.

*Thanks to: ITER Scientist Fellow Network; Spanish National Research Project No. PID2019-110734RB-I00; Excellence of University Professors (EPUC3M14); Uranus, a cluster funded jointly by EU FEDER funds and by the Spanish Government via national projects UNC313-4E-23612, ENE2009-12213-C03-03, ENE2012-33219, ENE2012-31753 and ENE2015-68265; EUROFUSION.

Publication: [1] G. F-Torija Daza, J.M. Reynolds-Barredo, R. Sanchez, A. Loarte, and G. Huijsmans. Free-plasma boundary solver for axisymmetric ideal mhd equilibria with flow. Nuclear Fusion, 62(12):126044, oct 2022.
[2] Ongoing paper including the last upgrades in ITER and NSTX showcases.
[3] J.M. Reynolds-Barredo, H. Peraza-Rodríguez, R. Sanchez, and V. Tribaldos. A novel efficient solver for ampere's equation in general toroidal topologies based on singular value decomposition techniques. Journal of Computational Physics, 406:109214, 2020.
[4] H. Peraza-Rodriguez, J. M. Reynolds-Barredo, R. Sanchez, J. Geiger, V. Tribaldos, S. P. Hirshman, and M. Cianciosa. Extension of the siesta mhd equilibrium code to free-plasma-boundary problems. Physics of Plasmas, 24(8):082516, 2017.

Presenters

  • Gonzalo Fernandez-Torija Daza

    • UC3M (Universidad Carlos III de Madrid)

Authors

  • Gonzalo Fernandez-Torija Daza

    • UC3M (Universidad Carlos III de Madrid)
  • Jose M Reynolds Barredo

    • Universidad Carlos III de Madrid
    • UC3M (Universidad Carlos III de Madrid)
  • Raul Sanchez

    • University Carlos III De Madrid
    • Univ Carlos III De Madrid
    • Universidad Carlos III de Madrid
  • Victor Tribaldos Macía

    • UC3M (Universidad Carlos III de Madrid)
  • Alberto Loarte

    • ITER Organization