APS Logo

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.

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.<br>[2] Ongoing paper including the last upgrades in ITER and NSTX showcases.<br>[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.<br>[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