In order to meet the three objectives of the project over its four-years span, namely achieve the scale up of our three solvers to tackle simulations in realistic ITER configurations, increase the reliability of low-fidelity models, and promote an interdisciplinary expertise and strong, the project has been structured in 6 tasks involving different partners.
Chart of the SISTEM project organisation showing how the main scientific tasks contribute
Task 1: Development of a new high-order finite element solver [G. Giorgiani, P1]
Contributed partners: 1 & 2
Objectives: Due to the fast decrease of the discretization error (minimization of the number of degree of freedom), and the great flexibility and accuracy in the description of any magnetic and tokamak geometries, this task aims to develop a high-order finite-element solver based on Hybrid-Discontinuous Galerkin method to succeed the scale-up of fluid simulations to ITER. This task in the core of the project.
Deliverables:
1: A 2D high-order non-isothermal fluid solver including neutral dynamics to simulate turbulent transport in ITER size edge plasma (WP1.1, WP1.2)
2: A 2D high-order non-isothermal fluid solver coupled to realistic source of energy in the core (WP1.3)
3: A 3D high-order isothermal fluid solver (WP1.4)
4: Publications in international journals: J. Comp. Phys., Phys. Of Plasma, ...
Task 2: Development of a high-order HDG version of CEDRES++/FEEQ.M [F. Rapetti (P3)]
Contributed partners: 1 & 3
Objectives: Upgrade the magnetic equilibrium solver with a high-order DG method will yield to a better description of the poloidal magnetic flux and this will increase the accuracy of the geometric coefficients calculation which will be of key importance in the future coupling between the fluid solver (that better simulates the interaction plasma-wall) and the equilibrium one (that better describes the plasma profile). That will be done in close connection with task1 for a future coupling between codes.
Deliverables:
1:A high-order release of CEDRES++/FEEQ.M based on HDG and preparing future coupling with fluid codes of task 1.
2:Publications in international journals: J. Comp. Phys., Phys. Of Plasma, ...
Task 3: Numerical and mathematical issues related to codes developments [M. MehrenbergerP1].
Contributed partners: 1, 2 & 3
Objectives: The development of new solvers able to simulate ITER size edge plasmas requires the optimization of different parts of the algorithm to achieve targeted performances. To ensure the reliability of such developments a strong mathematical background is required.
Deliverables:
1:A semi-implicit time integration scheme adapted to HDG formulation for codes of task 1.
2: Report on the different linear solvers evaluated during the task
3: A mesh mixing triangles and quadrangles to optimize the domain discretization
4: Publications in international journals
Task 4: Data assimilation to reduced models optimization [D. Auroux, P3]
Objectives: Explore the efficiency of the assimilation of data obtained from experimental measurements and/or 3D simulationstofixfree parameters of the transverse turbulent transport models in the set of averaged equations in 2D (an example of such a model is given in Eq. 1). The success of this task will require a close collaboration between all partners to strength the synergy between plasma physics and mathematics.
Contributed partners: 1, 2 & 3
Deliverables:
1:Reports on data assimilation techniques
2: Publications
Task 5: Simulations with non-steady magnetic configurations and related heat exhaust physics [Ph. Ghendrih, P2].
Objectives:They concern plasma simulations and the investigation of heat exhaust physics in original configurations allowed by the new capabilities of the codes developed in tasks 1 & 2. Objectives will be addressed in relation with the evolution of these codes during the project. They are more particularly related to heat fluxes estimate in large size tokamak like ITER and to the interplay between heat fluxes and magnetic equilibrium evolution.
Milestone 1: Non isothermal limiter/X-point transition 2D simulations. Using non steady CEDRES++/FEEQ.M equilibrium, 2D plasma simulations will be carried out. No coupling between the two codes. Analyse of the time evolution of the density.
Milestone 2: Non isothermal exploration of strike point sweeping and analysis on heat deposition
Milestone 3: First ITER size tokamak 3D simulation and turbulent transport properties analysis.