développement du code de simulation pour la fusion par confinement magnétique TOKAM3X
-en- Publications scientifiques au M2P2
2023
Thomas Cartier-Michaud, Philippe Ghendrih, Virginie Grandgirard, Eric Serre. Verification and accuracy check of simulations with PoPe and iPoPe. Journal of Computational Physics, 2023, 474, pp.111759. ⟨10.1016/j.jcp.2022.111759⟩. ⟨hal-03871954⟩ Plus de détails...
The theoretical background of the PoPe and iPoPe verification scheme is presented. Verification is performed using the simulation output of production runs. The computing overhead is estimated to be at most 10%. PoPe or iPoPe calculations can be done offline provided the necessary data is stored, for example additional time slices, or online where iPoPe is more effective. The computing overhead is mostly that of storing the necessary data. The numerical error is determined and split into a part proportional to the operators, which are combined to form the equations to be solved, thus modifying their control parameters, completed by a residual error orthogonal to these operators. The accuracy of the numerical solution is determined by this modification of the control parameters. The PoPe and iPoPe methods are illustrated in this paper with simulations of a simple mechanical system with chaotic trajectories evolving into a strange attractor with sensitivity to initial conditions. We show that the accuracy depends on the particular simulation both because the properties of the numerical solution depend on the values of the control parameter, and because the target accuracy will depend on the problem that is addressed. One shows that for a case close to bifurcations between different states, the accuracy is determined by the level of detail of the bifurcation phenomena one aims at describing. A unique verification index, the PoPe index, is proposed to characterise the accuracy, and consequently the verification, of each production run. The PoPe output allows one to step beyond verification and analyse for example the numerical scheme efficiency. For the chosen example at fixed PoPe index, therefore at fixed numerical error, one finds that the higher order integration scheme, comparing order 4 to order 2 Runge-Kutta time stepping, reduces the computation cost by a factor 4.
Thomas Cartier-Michaud, Philippe Ghendrih, Virginie Grandgirard, Eric Serre. Verification and accuracy check of simulations with PoPe and iPoPe. Journal of Computational Physics, 2023, 474, pp.111759. ⟨10.1016/j.jcp.2022.111759⟩. ⟨hal-03871954⟩
Raffaele Tatali, Eric Serre, Patrick Tamain, Davide Galassi, Philippe Ghendrih, et al.. Impact of collisionality on turbulence in the edge of tokamak plasma using 3D global simulations. Nuclear Fusion, 2021, ⟨10.1088/1741-4326/abe98b⟩. ⟨hal-03182318⟩ Plus de détails...
Collisionality is one of the key parameters in determining turbulent transport in the plasma edge, regulating phenomena such as "shoulder formation", separation of scale lengths in the scrape-off layer, turbulence damping and zonal flow dynamics. Understanding its role is therefore of primary importance for future reactors like ITER. Obtaining reliable predictions and a better characterization of plasma flow properties when varying collisionality remains, however, a critical challenge for the simulations. This paper focuses on the impact of varying collisionality in a nonisothermal three-dimensional fluid model of the plasma edge. A high field side limited configuration encompassing open and closed magnetic field lines with parameters typical of a medium-sized tokamak is considered. The present model can consistently account for the variations of collisionality and its impact on both the parallel resistivity η and the ion and electron parallel thermal conductivities χ e,i. Details on mean flow and turbulence properties are given. Changing collisionality leads to significant changes in the flow properties both on the mean and fluctuating quantities. In particular, lowering collisionality decreases the size of coherent structures, the fluctuation levels of turbulence, and steepens the density and temperature equilibrium profiles around the separatrix leading to a global reduction of the turbulent transport. The scrape-off layer (SOL) width is observed to increase with collisionality, eventually resulting in the disappearance of the scale lengths separation between near and far SOL, consistently with previous experimental observations. At low collisionality, where the presence of narrow feature is well-established, a contribution of heat conduction increases up to compete with heat convection.
Raffaele Tatali, Eric Serre, Patrick Tamain, Davide Galassi, Philippe Ghendrih, et al.. Impact of collisionality on turbulence in the edge of tokamak plasma using 3D global simulations. Nuclear Fusion, 2021, ⟨10.1088/1741-4326/abe98b⟩. ⟨hal-03182318⟩
T. Cartier-Michaud, D. Galassi, Ph Ghendrih, P. Tamain, F. Schwander, et al.. A posteriori error estimate in fluid simulations of turbulent edge plasmas for magnetic fusion in tokamak using the data mining iPoPe method. Physics of Plasmas, 2020. ⟨hal-02613800⟩ Plus de détails...
Progressing towards more reliable numerical solutions in the simulation of plasma for magnetic confinement fusion has become a critical issue for the success of the ITER operation. This requires developing rigorous and efficient methods of verification of the numerical simulations in any relevant flow regimes of the operation. The paper introduces a new formulation of the PoPe 1 method, namely the independent Projection on Proper elements method (iPoPe) to quantify the numerical error by performing a data-driven identification of the mathematical model from the simulation outputs. Based on a statistical postprocessing of the outputs database, the method provides a measure of the error by estimating the distance between the (numerical) effective and (analytical) theoretical weights of each operator implemented in the mathematical model. The efficiency of the present method is illustrated on turbulent edge plasma simulations based on a drift-reduced Braginskii fluid model in realistic magnetic geometries. Results show the effective order of the numerical method in these multiscale flow regimes as well as the values of the plasma parameters which can be safely simulated with respect to a given discretization. In this sense, the method goes one step further than the Method of Manufactured Solution (MMS 2-4), recently introduced in fusion, and provides an efficient verification procedure of the numerical simulations in any regimes, including turbulent ones that could be generalized to other scientific domains.
T. Cartier-Michaud, D. Galassi, Ph Ghendrih, P. Tamain, F. Schwander, et al.. A posteriori error estimate in fluid simulations of turbulent edge plasmas for magnetic fusion in tokamak using the data mining iPoPe method. Physics of Plasmas, 2020. ⟨hal-02613800⟩
Maxime Lesur, Chabha Djerroud, Kyungtak Lim, Etienne Gravier, Malik Idouakass, et al.. Validity limits of the passive treatment of impurities in gyrokinetic tokamak simulations. Nuclear Fusion, 2020, 60 (3), pp.036016. ⟨10.1088/1741-4326/ab6e48⟩. ⟨hal-02485083⟩ Plus de détails...
In gyrokinetic simulations of turbulent impurity transport, trace impurity species are often treated as passive species, in the sense that they are not included in Maxwell equations. This is consistent with the assumption that impurities with low enough concentrations are impacted by turbulence generated by electrons and main ions, but do not impact it significantly in return. In this work, we relax this assumption, and investigate the active impacts of impurity on impurity transport as a function of its concentration, in the presence of trapped-particle-driven turbulence. We focus on W 40+ tungsten, which is relevant for modern tokamaks, and adopt a reduced gyrokinetic bounce-averaged model for trapped particles in a simplified tokamak geometry. The impacts depend on the relationship between equilibrium density gradient and temperature gradient. When these gradients are equal, we observe that tungsten can be treated as a passive species for concentrations below . Above this concentration, the impurity significantly impacts both density and heat transport, essentially quenching them for concentrations above 10−3. This quenching occurs as electric potential fluctuations become in phase with impurity density fluctuations.
Maxime Lesur, Chabha Djerroud, Kyungtak Lim, Etienne Gravier, Malik Idouakass, et al.. Validity limits of the passive treatment of impurities in gyrokinetic tokamak simulations. Nuclear Fusion, 2020, 60 (3), pp.036016. ⟨10.1088/1741-4326/ab6e48⟩. ⟨hal-02485083⟩
Philippe Ghendrih, Y. Asahi, E. Caschera, Guilhem Dif-Pradalier, P. Peter Donnel, et al.. Generation and dynamics of SOL corrugated profiles. Journal of Physics: Conference Series, 2018, 1125, pp.012011. ⟨10.1088/1742-6596/1125/1/012011⟩. ⟨hal-02196677⟩ Plus de détails...
The staircase transport regime reported in kinetic simulations of plasma turbulent transport in magnetic confinement is recovered with a simple 2D fluid model allowing for reduced damping of the zonal flows. Some of the complex dynamics of the kinetic zonation regime are recovered but the pattern of the corrugation appears to be sinusoidal with a characteristic scale comparable to that of turbulence modes with largest spectral energy, in contrast to regimes observed in global and flux-driven kinetic simulations. Enhanced zonal flows govern both an overall reduction of the SOL width and a gradual steepening of the gradients with distance to the separatrix.
Philippe Ghendrih, Y. Asahi, E. Caschera, Guilhem Dif-Pradalier, P. Peter Donnel, et al.. Generation and dynamics of SOL corrugated profiles. Journal of Physics: Conference Series, 2018, 1125, pp.012011. ⟨10.1088/1742-6596/1125/1/012011⟩. ⟨hal-02196677⟩
D Rigamonti, L. Giacomelli, G Gorini, M. Nocente, M. Rebai, et al.. Neutron spectroscopy measurements of 14 MeV neutrons at unprecedented energy resolution and implications for deuterium–tritium fusion plasma diagnostics. Measurement Science and Technology, 2018, 29 (4), pp.045502. ⟨10.1088/1361-6501/aaa675⟩. ⟨hal-02177077⟩ Plus de détails...
An accurate calibration of the JET neutron diagnostics with a 14 MeV neutron generator was performed in the first half of 2017 in order to provide a reliable measurement of the fusion power during the next JET deuterium-tritium (DT) campaign. In order to meet the target accuracy, the chosen neutron generator has been fully characterized at the Neutron Metrology Laboratory of the National Physical Laboratory (NPL), Teddington, United Kingdom. The present paper describes the measurements of the neutron energy spectra obtained using a high-resolution single-crystal diamond detector (SCD). The measurements, together with a new neutron source routine 'ad hoc' developed for the MCNP code, allowed the complex features of the neutron energy spectra resulting from the mixed D/T beam ions interacting with the T/D target nuclei to be resolved for the first time. From the spectral analysis a quantitative estimation of the beam ion composition has been made. The unprecedented intrinsic energy resolution (<1% full width at half maximum (FWHM) at 14 MeV) of diamond detectors opens up new prospects for diagnosing DT plasmas, such as, for instance, the possibility to study non-classical slowing down of the beam ions by neutron spectroscopy on ITER.
D Rigamonti, L. Giacomelli, G Gorini, M. Nocente, M. Rebai, et al.. Neutron spectroscopy measurements of 14 MeV neutrons at unprecedented energy resolution and implications for deuterium–tritium fusion plasma diagnostics. Measurement Science and Technology, 2018, 29 (4), pp.045502. ⟨10.1088/1361-6501/aaa675⟩. ⟨hal-02177077⟩
Thomas Cartier-Michaud, Philippe Ghendrih, Guilhem Dif-Pradalier, Xavier Garbet, Virginie Grandgirard, et al.. Verification of turbulent simulations using PoPe: quantifying model precision and numerical error with data mining of simulation output. Journal of Physics: Conference Series, 2018, 1125, pp.012005. ⟨10.1088/1742-6596/1125/1/012005⟩. ⟨hal-02196674⟩ Plus de détails...
Verification of a 1D-1V kinetic code with the PoPe method [1] is presented. Investigation of the impact of reducing the precision of the numerical scheme is analysed by following 3 indicators of the physics solved by the code, namely the plasma response to an external high frequency electric field wave. The response of the distribution function in the vicinity of the particle-wave resonance is found to be most sensitive to the resolution. Consistently, a rapid growth of the error indicator determined with PoPe is observed. However, no critical value of this indicator allowing us to retain the physics in a situation of degraded precision could be observed. The response of the amplitude of the electric potential fluctuations is characterised by a transient growth followed by a plateau. It is found that the loss of this plateau is governed by the resolution in v-space, but due to the generation of a symmetry in the problem rather than to errors in the numerical scheme. The analysis of the transient indicates that the growth rate of the amplitude of the electric potential is very robust down to very low resolution, step in velocity of 2 thermal velocities. However, a transition prior to this resolution, with step 0.5 thermal velocity, can be identified corresponding to a PoPe indicator of order zero, namely for errors of order 100 %.
Thomas Cartier-Michaud, Philippe Ghendrih, Guilhem Dif-Pradalier, Xavier Garbet, Virginie Grandgirard, et al.. Verification of turbulent simulations using PoPe: quantifying model precision and numerical error with data mining of simulation output. Journal of Physics: Conference Series, 2018, 1125, pp.012005. ⟨10.1088/1742-6596/1125/1/012005⟩. ⟨hal-02196674⟩
Eric Serre, Sandrine Hugues, Emilia Crespo del Arco, Anthony Randriamampianina, Patrick Bontoux. Axisymmetric and three-dimensional instabilities in an Ekman boundary layer flow. International Journal of Heat and Fluid Flow, 2001, 22 (1), pp.82-93. ⟨hal-01023080⟩ Plus de détails...
Eric Serre, Sandrine Hugues, Emilia Crespo del Arco, Anthony Randriamampianina, Patrick Bontoux. Axisymmetric and three-dimensional instabilities in an Ekman boundary layer flow. International Journal of Heat and Fluid Flow, 2001, 22 (1), pp.82-93. ⟨hal-01023080⟩
Journal: International Journal of Heat and Fluid Flow
