Skip to Main content Skip to Navigation

Deux étapes majeures pour le développement du code XTOR : parallélisation poussée et géométrie à frontière libre.

Abstract : The XTOR-2F code simulates the 3D dynamics of full bi-fluid MHD instabilities in tokamak plasmas.The first part of the thesis was dedicated to the parallelisation of XTOR-2F code. The code has been parallelised significantly despite the numerical profile of the problem solved, i.e. a discretisation with pseudo-spectral representations in all angular directions, the stiffness of the two-fluid stability problem in tokamaks, and the use of a direct LU decomposition to invert the physical pre-conditioner. The execution time of the parallelised version is an order of magnitude smaller than the sequential one for low-resolution cases, with an increasing speedup when the discretisation mesh is refined. Moreover, it allows to perform simulations with higher resolutions, previously forbidden because of memory limitations.The second part of the thesis was dedicated to the development of free boundary condition. The original fixed boundary computational domain of the code was generalised to a free-boundary one, thus approaching closely the geometry of today’s and future large experimental devices. The initial conditions are given by the CHEASE equilibrium code inside the plasma. Outside the plasma, fitting the magnetic potential at the CHEASE computation domain boundary with a set of external poloidal magnetic coils extends the solution. The boundary conditions use Green functions to construct a response matrix matching the normal and tangential components of the outside magnetic field with the inside solution. A thin resistive wall can be added to the computational domain. This new numerical setup generalises the investigation field from internal MHD instabilities towards external instabilities. The code linear behaviour is validated with two families of instabilities, n = 0 axisymmetric modes and n = 1/m = 2 external kinks. In order to validate the nonlinear behaviour, nonlinear resistive MHD simulations of tearing modes at zero beta evolving to a stationary state have been performed.
Complete list of metadata

Cited literature [58 references]  Display  Hide  Download
Contributor : ABES STAR :  Contact
Submitted on : Monday, February 19, 2018 - 4:16:05 PM
Last modification on : Monday, June 14, 2021 - 10:41:28 AM
Long-term archiving on: : Sunday, May 6, 2018 - 7:17:58 AM


Version validated by the jury (STAR)


  • HAL Id : tel-01712594, version 1


Alain Marx. Deux étapes majeures pour le développement du code XTOR : parallélisation poussée et géométrie à frontière libre.. Physique des plasmas [physics.plasm-ph]. Université Paris Saclay (COmUE), 2017. Français. ⟨NNT : 2017SACLX095⟩. ⟨tel-01712594⟩



Record views


Files downloads