# Keywords

- A6. Modeling, simulation and control
- A6.1. Methods in mathematical modeling
- A6.1.1. Continuous Modeling (PDE, ODE)
- A6.1.4. Multiscale modeling
- A6.1.5. Multiphysics modeling
- A6.2. Scientific computing, Numerical Analysis & Optimization
- A6.2.1. Numerical analysis of PDE and ODE
- A6.2.6. Optimization
- A6.2.7. High performance computing
- A6.2.8. Computational geometry and meshes
- A6.3. Computation-data interaction
- A6.3.1. Inverse problems
- A6.3.2. Data assimilation
- A6.3.4. Model reduction
- A6.4. Automatic control
- A6.4.1. Deterministic control
- A6.4.4. Stability and Stabilization
- A6.5. Mathematical modeling for physical sciences

- B4. Energy
- B4.2.2. Fusion

# 1 Team members, visitors, external collaborators

## Research Scientists

- Hervé Guillard [Team leader, Inria, Senior Researcher, HDR]
- Florence Marcotte [Inria, Researcher]
- Sebastian Minjeaud [CNRS, Researcher]
- Richard Pasquetti [CNRS, Senior Researcher, HDR]

## Faculty Members

- Didier Auroux [Univ Côte d'Azur, Professor, HDR]
- Jacques Blum [Univ Côte d'Azur, Emeritus, HDR]
- Cédric Boulbe [Univ Côte d'Azur, Associate Professor]
- Boniface Nkonga [Univ Côte d'Azur, Professor]
- Francesca Rapetti [Univ Côte d'Azur, Associate Professor]
- Afeintou Sangam [Univ Côte d'Azur, Associate Professor]

## Post-Doctoral Fellow

- Paul Mannix [Inria, from Sep 2020]

## PhD Students

- Ashish Bhole [Univ Côte d'Azur]
- Ali Aboudou Elarif [Inria]
- Louis Lamerand [Univ Côte d'Azur]

## Technical Staff

- Blaise Faugeras [CNRS, Engineer]

## Administrative Assistant

- Montserrat Argente [Inria]

# 2 Overall objectives

In order to fulfill the increasing demand, alternative energy sources have to be developed. Indeed, the current rate of fossil fuel usage and its serious adverse environmental impacts (pollution, greenhouse gas emissions, ...) lead to an energy crisis accompanied by potentially disastrous global climate changes.

Controlled fusion power is one of the most promising
alternatives to the use of fossil resources, potentially
with a unlimited source of fuel. France with the ITER
(http://

The other route for fusion plasma is inertial confinement. In this latter case, large scale hydrodynamical instabilities prevent a sufficiently large energy deposit and lower the return of the target. Therefore, for both magnetic and inertial confinement technologies, the success of the projects is deeply linked to the theoretical understanding of plasma turbulence and flow instabilities as well as to mathematical and numerical improvements enabling the development of predictive simulation tools.

Castor gathers the activities in numerical simulation
of fusion plasmas with the activities in control and optimisation done in the laboratory
Jean-Alexandre Dieudonné of Université Côte d'Azur.
The main objective of the Castor team is to contribute
to the development of innovative numerical tools to improve
the computer simulations of complex turbulent or unstable flows in plasma physics and to develop
methods allowing the real-time control of these flows or the optimisation of scenarios of plasma discharges in tokamaks.
Castor is a common project between Inria (http://

# 3 Research program

## 3.1 Plasma Physics

The main reseach topics are:

- Modelling and analysis
- Fluid closure in plasma
- Turbulence
- Plasma anisotropy type instabilities
- Free boundary equilibrium (FBE)
- Coupling FBE – Transport
- MHD instabilities

- Numerical methods and simulations
- High order methods
- Curvilinear coordinate systems
- Equilibrium simulation
- Pressure correction scheme
- Anisotropy
- Solving methods and parallelism

- Identification and control
- Inverse problem: Equilibrium reconstruction
- Open loop control

# 4 Application domains

## 4.1 MHD and plasma stability in tokamaks

Participants: Hervé Guillard, Boniface Nkonga, Afeintou Sangam, Ali Elarif, Ashish Bhole.

The magnetic equilibrium in tokamaks results from a balance between the Lorentz force and the pressure gradient. Using Ampère law, a convenient description of this equilibrium is provided by the Grad-Shafranov equation. Of course, the magnetic equilibrium solution of the Grad-Shafranov equation is required to be stable. Actually any loss of MHD stability can lead to the end of the existence of the plasma, the so-called disruptions that can affect negatively the integrity of the machine. The primary goal of MHD (Magneto-Hydro-Dynamics) studies is therefore to determine the stability domain that constraints the operational range of the machine.A secondary goal of MHD studies is to evaluate the consequences of possible disruptions in term of heat loads and stresses on the plasma facing components. In modern machines in the so-called H-mode some mild instabilities leading to a near oscillatory behavior are also known to exist. In particular, the so-called ELMs ( Edge Localized Modes) are of particular importance since they can have large effects on the plasma facing components. The control and understanding of these instabilities is therefore of crucial importance for the design of future machines as ITER. Unfortunately, ELM occur in the edge plasma and their modeling requires to take in account not only the intricate magnetic topology of this region where co-exist both open and closed field lines but also the existence of molecular and atomic processes involving neutrals.

At present, the linear theory of MHD stability is relatively well understood. However, the description of the non-linear behavior is far from being complete. As a consequence and due to the intrinsic difficulty of the subject, only a few numerical codes worldwide have been developed and validated for non linear MHD in tokamaks. One of these codes is the JOREK code developped since 2006 from a collaborative work between CEA-Cadarache (main developer), LABRI Bordeaux, LJAD-UCA and Inria. A comprehensive description of JOREK is given in 17

## 4.2 Long term plasma evolution and optimization of scenarii

Participants: Didier Auroux, Jacques Blum, Cédric Boulbe, Blaise Faugeras, Hervé Guillard.

The magnetic equilibrium evolves in time due to diffusion processes on the slow resistive diffusive time scale and moreover it has to be monitored with active and passive control based on external coils, current drive, heating system, particle or pellets injections. This set of control mechanism has to be modeled and this is the goal of real time codes or global evolution codes.In the same order of ideas, the steering and control of the plasma from the beginning to the end of the discharge require the research of optimal trajectories through the space of operational parameters. This is usually performed in an empirical way in present Tokamaks, but the complexity of the problem requires today the use of optimization techniques for processes governed by MHD and diffusion-type equations.

## 4.3 Turbulence and models for the edge region of tokamaks

Participants: Didier Auroux, Louis Lamerand, Francesca Rapetti.

The edge region of the plasma is characterized by low temperature and density leading to an increase of the collision frequency that makes the edge plasma nearly collisional. This combined with the intricate magnetic topology of this region makes the development of kinetic codes adapted to the edge regions a real long term adventure. Consequently the fluid approach remains a standard one to study edge plasma turbulence. The use of optimal control theory to derive simplified models matching data either experimental or derived from direct numerical simulations is part of the objectives of the team.## 4.4 Understanding magnetogenesis in stellar systems

Participants: Didier Auroux, Paul Mannix, Florence Marcotte.

The considerable diversity of long-lived magnetic fields observed in the Universe raises fundamental questions regarding their origin. Although it is now widely accepted that such fields are sustained by a dynamo instability in the electrically conducting fluid layers of astrophysical bodies, in most cases the very nature of the flow motions powering the dynamo is essentially unknown, and the conditions required for amplifying large-scale magnetic fields in non-convective stellar systems are poorly understood. We claim that optimal control represents a powerful tool to investigate the nonlinear stability of fully 3D, unsteady magnetohydrodynamic flows with respect to the dynamo instability. Nonlinear optimisation can be also used as a physical diagnostic to gain novel understanding of the mechanisms that are most favorable to dynamo action in a natural system.

# 5 New results

## 5.1 An overview paper on the equilibrium code NICE

Participants: Blaise Faugeras.

The code NICE (Newton direct and Inverse Computation for Equilibrium) enables to solve numerically several problems of plasma free-boundary equilibrium computations in a tokamak: plasma free-boundary only reconstruction and magnetic measurements interpolation, full free-boundary equilibrium reconstruction from magnetic measurements and possibly internal measurements (interferometry, classical linear approximation polarimetry or Stokes model polarimety, Motional Stark Effect and pressure), direct and inverse, static and quasi-static free-boundary equilibrium computations.

NICE unifies and upgrades 3 former codes VacTH, EQUINOX and CEDRES++. The strength of NICE is to gather in a single finite element framework different equilibrium computation modes. It makes intensive use of Newton method and Sequential Quadratic Programming method to solve non linear problems.

NICE is used routinely for WEST tokamak operation. It is also adapted to the IMAS (ITER Modelling and Analysis Suite) format which makes it usable on many different fusion tokamak reactors.

In this document we give a general overview of the numerical methods implemented in NICE as well as a number of computation examples 7.

## 5.2 Introduction of C1 Clought-Tocher finite elements in NICE

Participants: Ali Elarif, Blaise Faugeras, Francesca Rapetti.

The numerical simulation of the equilibrium of the plasma in a tokamak as well as its self-consistent coupling with resistive diffusion should benefit from higher regularity of the approximation of the magnetic flux map. In this work, we propose a finite element approach on a triangular mesh of the poloidal section, that couples piece-wise linear finite elements in a region that does not contain the plasma and reduced Hsieh-Clough-Tocher finite elements elsewhere. This approach gives the flexibility to achieve easily and at low cost higher order regularity for the approximation of the flux function in the domain covered by the plasma, while preserving accurate meshing of the geometric details in the rest of the computational domain. The continuity of the numerical solution at the coupling interface is weakly enforced by mortar projection. A new technique for the computation of the geometrical coefficients is also presented (15)

## 5.3 Introduction of a resistive diffusion equation for the diamagnetic function in NICE evolutive equilibrium computations

Participants: Blaise Faugeras.

This work is based on an unpublished paper: A weak formulation of current diffusion in tokamak plasmas by H. Heumann, former member of the team. First numerical tests are encouraging and the method needs now to be further tested and validated.

## 5.4 Equilibrium reconstruction at JET

Participants: Blaise Faugeras.

The brand-new RAPTOR suite of codes is here presented. The suite has been developed for the JET tokamak to combine real-time model-based predictions of the plasma state with the available diagnostic measurements. The suite embeds: the upgraded equilibrium reconstruction EQUINOX code, the FLUXMAP algorithm, which maps the diagnostic measurements from geometric to normalized magnetic flux coordinates; the RABBIT code for the NBI reconstruction and eventually RAPTOR state observer, which combines the output from all these codes with the predictions of 1D control-oriented transport code. The suite is both implemented in MATLAB/Simulink and it is being integrated in the C++ real-time MARTe2 framework. Thanks to its user-friendly interfaces, which are based on the MDSplus I/O and visualization tools, the RAPTOR suite can be used both offline, for a fast reconstruction of the plasma state, and in integrated control algorithms once it will be deployed in the JET real-time data network.

## 5.5 Coupling NICE-METIS

Participants: Cédric Boulbe, Blaise Faugeras, Jean François Artaud(CEA).

Started last year, the coupling between the free boundary equilibrium code and the fast transport solver Metis (CEA) has been improved. In particular new outputs have been implemented in Nice in order to implement a version of the worklow using mixed boundary conditions in the diffusion equation sovled by Metis. This new version of the workflow has been tested on an ITER test case.

In order to use this workflow in the ITER IMAS framework, the static, evolutive and inverse static mode of Nice have also been adapted to IMAS.

## 5.6 High-order Whitney finite elements

Participants: Francesca Rapetti.

In this work, a local complete fully discrete de Rham sequence on polygonal and polyhedric elements has been defined, in collaboration with Daniele Di Pietro (Univ. Montpellier, France) and Jerome Droniou (Monash Univ., Australia) 6. In collaboration with Ana Alonso Rodriguez (Univ. Trento, Italy), the definition of the Lebesgue constant to the case of field interpolation by high-order Whitney forms on simplices has been generalized 4.

## 5.7 Extended MHD Modeling of the Tokamak plasmas

Participants: Ashish Bhole, Boniface Nkonga, Stanislas Pamela (UKAEA-Culham).

JOREK has been a useful numerical tool to simulate Magneto-hydrodynamic (MHD) instabilities with reduced-MHD models in tokamaks. One of the assumptions behind reduced-MHD model is remove fast magneto-sonic waves from the model. In full-MHD model, these fast magneto-sonic waves are retained and hence its implementation in JOREK is thought to offer a more complete physics modeling. Preliminary simulations for full-MHD models have been performed in Haverkort et al. (2016 JCP). Further advancement is made in the context of full-MHD to demonstrate the capabilities of full-MHD at production level in JOREK 13. Additionally, extended physics terms are implemented in full-MHD model. The simulations of full-MHD model are validated and compared with reduced-MHD models and benchmark is performed.The contributions are made from the point of view of modeling (reduced and full-MHD) and numerical aspects to the extensive review paper 17. This review paper includes some representative results obtain using reduced-MHD model with a wide range of physical scenarios. The next task is to be able to simulate such physical scenarios using full-MHD model and compare them with results of reduced-MHD model.

## 5.8 Mesh singularity treatment and corrected Gears Method

Participants: Ashish Bhole, Boniface Nkonga.

For most MHD simulations, JOREK uses polar grids with a singularity at the grid-center. The C1 continuity achieved in Finite Element Method of JOREK is lost at center of polar grid. This singular grid-center is known to be a source of spurious nonphysical waves that can pollute the simulation or even crash it. The strategy is developed to achieve regularity at the polar grid-center that reduces the noise in numerical simulations significantly.In realistic JOREK simulations, time steps maybe increased or reduced depending upon numerical requirements. The second order backward difference method (BDF2 or Gears) for time integration in JOREK was not adapted for changing time steps. The correction to adapt Gears method in JOREK for changing time-steps is now done.

Validation of the Variable Multi-scale Stabilized (VMS) method in JOREK is an ongoing task.

## 5.9 A path conservative finite volume method for a shear shallow water model

Participants: Praveen Chandrashekar(Tata Institute Bangalore), Boniface Nkonga, Ashish Bhole.

The shear shallow water model provides an approximation for shallow water flows by including the effect of vertical shear in the model. This model can be derived from the depth averaging process by including the second order velocity fluctuations, which are neglected in the classical shallow water approximation. The resulting model has a non-conservative structure, which resembles the 10-moment equations from gas dynamics. This structure facilitates the development of path conservative schemes and we construct HLL, 3-wave and 5-wave HLLC-type solvers. An explicit and semi-implicit MUSCL-Hancock type second order scheme is proposed for the time integration. Several test cases including roll waves show the performance of the proposed modeling and numerical strategy 5.

Ongoing work will derive an exact solution when approximated jump conditions are assumed.

## 5.10 Nonlinear stability of MHD flows: minimal dynamo seeds

Participants: Didier Auroux, Paul Mannix, Florence Marcotte.

This project consists in the development of a first adjoint-based optimisation code for fully nonlinear, 3D, unsteady MHD flows (governed by the fully nonlinear Navier-Stokes equations and induction equation, subject to solenoidal constraints on both the velocity and magnetic fields). This code is currently in development, and partially builds on the open-source framework Dedalus (dedalus-project.org) for spectrally solving PDEs. It is aimed at optimizing initial magnetic fields perturbations leading to maximal field amplification by a given target time, in a cylindrical Taylor-Couette flow. This optimisation will be used as an intermediate step toward performing nonlinear stability analysis of the considered flow with respect to the dynamo instability, a study which we anticipate will be transposable to various kinds of MHD instabilities and geometries in the future (e.g. possible long-term application to the ELMs instabilities in fusion devices). The direct and adjoint solvers of the developed code have been validated and various optimization algorithms are currently being tested.

## 5.11 Effect of spatially variable electrical conductivity on dynamo instability threshold

Participants: Florence Marcotte.

Temperature or chemical composition gradients play an important role in the dynamics of astrophysical plasmas. They are also associated with considerable spatial variations of the electrical conductivity, which are likely to impact on the development of magnetohydrodynamic instablities and the generation of long-lived astrophysical magnetic fields. This numerical study shows that finite-amplitude conductivity variations (with an amplitude of the same order as the mean value) can significantly lower the kinematic threshold for dynamo instability, allowing for large-scale magnetic fields to be exponentially amplified in a parameter regime where they would be entirely dissipated away otherwise. This work relies on a staggered-grid, semi-implicit, finite-differences solver for kinematic induction in 3D cylindrical geometry developed by F.Marcotte, and is to be submitted to Europhysics Letters in the weeks to come.

## 5.12 Dynamo instabilities in radiative stellar interiors

Participants: Florence Marcotte.

This work is a collaboration with Ludovic Petitdemange (LERMA, Observatoire de Paris) and Christophe Gissinger (Ecole Normale Supérieure, Paris).

Understanding the conditions of magnetic field generation in radiative (stably-stratified) stars is a long-lasting challenge, which has important implication for the reconstruction of stellar rotation profiles and evolution. So far, the only models that have been considered to parametrize the effect of magnetism in available stellar evolution codes are either fossil magnetic fields or Tayler-Spruit dynamo fields. Yet, direct numerical simulations have consistently failed so far in exhibiting Tayler-Spruit dynamos, or in fact any type of dynamo mechanism that could operate in a stratified, hydrodynamically stable stellar layer. We numerically demonstrate how a dynamo instability can build up and trigger MHD turbulence in a radiative star, and show that the dynamo mechanism actually bears striking similarities with the elusive Tayler-Spruit model. A first article is to be submitted soon; a second article is in preparation to present a parametric study of the dynamo transition in the presence of stable temperature stratification.

## 5.13 Adjoint code for a reduced MHD model

Participants: Hervé Guillard, Ali Elarif, Laurent Hascoet(ECUADOR).

Non linear MHD simulations for tokamaks are now mature enough to be used for control, optimization or data assimilation purposes. The use of these tool can largely benefit from the possibility to compute gradients of the results of the simulations with respect to initial data or parameters. In collaboration with the ECUADOR team, we have begun to test the use of the automatic differentiation tool TAPENADE 16 to obtain the adjoint code of a reduced MHD model.

## 5.14 Sensitivity of Numerical Simulations of airborne droplets

Participants: Christophe Henry(CALISTO), Kerlyns Martinez-Rodriguez(CALISTO), Mireille Bossy(CALISTO), Nicolas Rutard(ONERA), Angelo Murrone(ONERA), Hervé Guillard.

To participate in the answers to the pandemic issues, researchers from Inria and the French Aerospace Lab ONERA have conducted a joint project to assess the relevance of the advices for social distancing. Numerical simulations of airborne droplet dispersion have been realized and analyzed with sensitivity analysis tools. It was found that the variability of the results depends on a number of factors, including: physical variables (e.g. droplet size, ejection velocity), modelling methods used (e.g. turbulence model) and numerical aspects (mesh) 12# 6 Partnerships and cooperations

## 6.1 International initiatives

### 6.1.1 Inria international partners

#### Informal international partners

The team collaborates with TUC (Technical University of Crete, Prof. Argyris Delis) on the modelling of acoustic streaming phenomena and the simulation of fast flood events.

## 6.2 Visits of international scientists

Stanislas Pamela from UKAEA (UK Atomic Energy Authority) visited the team in January. The collaboration was centered on the development of the JOREK code.

## 6.3 European initiatives

### 6.3.1 FP7 & H2020 Projects

CASTOR participates to the following EuroFusion consortium projects :

- EUROfusion WPCD (Working Package Code Development):
- EuroFusion Enabling Research CfP-AWP19-ENR-01, Strengthening the non-linear MHD code JOREK for application to key questions of the fusion roadmap.
- EUROfusion WPSA(Work Package JT-60SA)

## 6.4 National initiatives

### 6.4.1 ANR Sistem

Member of the ANR SISTEM , Oct. 2019 - Sept. 2023 coordinated by the M2P2 Institute of Aix-Marseille Univ. "SImulations with high-order schemes of tranSport and TurbulencE in tokaMak" programme Modeles numeriques 2019, Contact : Francesca Rapetti

Participants: Didier Auroux, Jacques Blum, Cédric Boulbe, Francesca Rapetti, Blaise Faugeras.

### 6.4.2 ANR DYNSEED

"Graines minimales de dynamos célestes" May 2020-2022, PI : Florence Marcotte.

Participants: Didier Auroux, Paul Mannix, Florence Marcotte.

# 7 Dissemination

## 7.1 Promoting scientific activities

### 7.1.1 Scientific events: organisation

- F. Marcotte: co-organizer of the Rencontres Nicoises de Dynamique des Fluides, budget: 4400
(with Giorgio Krstulovic, Yannick Ponty (OCA) ; Hélène Politano (LJAD); Sergey Nazarenko (INPHYNI) ; Jérémie Bec (CEMEF)) — event postponed due to Covid pandemic

### 7.1.2 Journal

#### Member of the editorial boards

- J. Blum is member of the editorial board of the Journal of Scientific Computing by Springer.
- F. Rapetti is member of the editorial board of the Advances in Computational Mathematics (ACOM) journal by Springer
- C. Boulbe is managing editor of the SMAI Journal of Computational Mathematic

#### Reviewer - reviewing activities

- H. Guillard has been reviewer for the Journal of Computational physics, Computers and Fluids and International Journal for Numerical methods in Fluids.
- F. Marcotte has been reviewer for the Journal of Fluid Mechanics, Physical Review Fluids, and Proceedings of the Royal Society of London A.

### 7.1.3 Invited talks

- High order Whitney forms on simplices, Seminar at the Montpellier Univ., 24 January 2020. Francesca Rapetti.
- High order Whitney forms on simplices, web Seminar, Oxford Univ., 12 November 2020. Francesca Rapetti.

### 7.1.4 Research administration

- Boniface Nkonga is member of the executive committee of the ECCOMAS (European Community on Computational Methods in Applied Sciences) association and treasurer of the association.
- H. Guillard represents Inria and UCA in the board of the "Fédération de Recherche, Fusion par Confinement Magnétique" (FR-FCM).
- H. Guillard is coordinator of the topic “Turbulence and transport of edge plasma” within the Fédération FR-FCM.

## 7.2 Teaching - Supervision - Juries

### 7.2.1 Teaching

- Licence : F. Rapetti, Mathématiques 2,30h équivalent TD, L2, Université Côte d'Azur
- Master : F. Rapetti, Méthodes numériques, 70h équivalent TD, M1, Université Côte d'Azur
- Ecole d'ingénieur, C. Boulbe, Analyse Numérique 1, 45.5h équivalent TD, niveau L3, Université Côte d'Azur
- Ecole d'ingénieur, C. Boulbe, Analyse numérique 2, 45.5h équivalent TD, niveau L3, Université Côte d'Azur
- Ecole d'ingénieur, C. Boulbe, Algèbre linéaire et Scilab, 26h équivalent TD, Université Côte d'Azur
- Ecole d'ingénieur, C. Boulbe, Projet 1, 48h équivalent TD, Université Côte d'Azur
- Licence: A. Sangam, Mathématiques Compléments d'Analyse, 54h, Semestre 3 de la Licence, Université Côte d'Azur
- Licence: A. Sangam, Analyse Numérique, 86h, L3 Mathématiques, Université Côte d'Azur
- Master: A. Sangam, Équations aux Dérivées Partielles & Différences Finies, 72h, Semestre 1 du Master M1 Mathématiques, Université Côte d'Azur
- Ecole d'ingénieur/Master: B. Nkonga, Méthode des éléments finis, 24h, M2, Polytech Nice Sophie, Université Côte d'Azur
- Ecole d'ingénieur/Master: B. Nkonga, Eléments finis mixtes, 24h, M2, Polytech Nice Sophia, Université Côte d'Azur
- Ecole d'ingénieur/Master, H. Guillard,Développement durable et enjeux de gouvernance, 16h, M2, Polytech Nice Sophia, Université Côte d'Azur

### 7.2.2 Supervision

- PhD : Ali Elarif, “Approximation par éléments finis C1 de modèles magnétohydrodynamique pour les plasmas de fusion”, December 17, 2020, Hervé Guillard
- PhD in progress : Ashish Bhole, supervisor : Boniface Nkonga. "Hermite-Bezier finite element method for Fusion plasma"
- PhD in progress: Louis Lamerand, PhD student fully payed by the ANR project SISTEM, who is working with Didier Auroux and Francesca Rapetti on “data assimilation” and “model reduction”: it consists in developing a simplified model for a physical phenomenon, here the heat transport in a tokamak, whose parameters are calibrated through either experimental measurements or accurate long-run computations with existing codes. With respect to the existing codes, the one based on the reduced model should provide an accurate and physically meaningful solution in a much shorter time, since october 2020.
- Postdoctoral researcher: Paul Mannix (Inria), started September 1st 2020, fully funded by ANR DYNSEED (PI: Florence Marcotte)

### 7.2.3 Juries

- Blaise Faugeras was part of the PhD defense jury of Francesco Carpanese, EPFL, September 30, 2020
- Boniface Nkonga was part of the PhD defense jury of Ahmed Blida, Université Côte d’Azur, July 2, 2020.

## 7.3 Popularization

F. Marcotte: Seminar in Rendez-vous des Jeunes Mathématiciennes et Informaticiennes, Nice, December 5th 2020, organized by Université Côte d’Azur, Animath, and Femmes & Mathémathiques

# 8 Scientific production

## 8.1 Major publications

- 1 articleReconstruction of the equilibrium of the plasma in a Tokamak and identification of the current density profile in real timeJournal of Computational Physics2312012, 960-980
- 2 article Fully discrete polynomial de Rham sequences of arbitrary degree on polygons and polyhedra Mathematical Models and Methods in Applied Sciences August 2020
- 3 articleExtended full-MHD simulation of non-linear instabilities in tokamak plasmasPhysics of Plasmas2710October 2020, 102510

## 8.2 Publications of the year

### International journals

- 4 article On a generalization of the Lebesgue's constant Journal of Computational Physics 2020
- 5 articleA path conservative finite volume method for a shear shallow water modelJournal of Computational Physics413July 2020, 109457
- 6 article Fully discrete polynomial de Rham sequences of arbitrary degree on polygons and polyhedra Mathematical Models and Methods in Applied Sciences August 2020
- 7 articleAn overview of the numerical methods for tokamak plasma equilibrium computation implemented in the NICE codeFusion Engineering and Design1602020, 112020
- 8 articleStationary shock-like transition fronts in dispersive systemsNonlinearity3310October 2020, 5477-5509
- 9 articleAdvances in the physics studies for the JT-60SA tokamak exploitation and research planPlasma Physics and Controlled Fusion621January 2020, 014009
- 10 articleAn asymptotic preserving scheme on staggered grids for thebarotropic Euler system in low Mach regimesNumerical Methods for Partial Differential Equations3652020, 1098-1128
- 11 articleAn explicit MUSCL scheme on staggered grids with kinetic-like fluxes for the barotropic and full Euler systemCommunications in Computational Physics2732020, 672-724
- 12 article Social Distancing: The Sensitivity of Numerical Simulations ERCIM News 2021 124 2021
- 13 articleExtended full-MHD simulation of non-linear instabilities in tokamak plasmasPhysics of Plasmas2710October 2020, 102510
- 14 articleElliptic problems in curved domains using cubature points based triangular spectral elements and isoparametric mappingsJournal of Mathematical Study533September 2020, 237-246

### Reports & preprints

- 15 report Tokamak free-boundary plasma equilibrium computation using finite elements of class C0 and C1 within a mortar element approach INRIA Sophia Antipolis October 2020

## 8.3 Cited publications

- 16 articleThe Tapenade Automatic Differentiation tool: Principles, Model, and SpecificationACM Transactions On Mathematical Software3932013, URL: http://dx.doi.org/10.1145/2450153.2450158
- 17 misc The JOREK non-linear extended MHD code and applications to large-scale instabilities and their control in magnetically confined fusion plasmas 2020