An accurate SUPG-stabilized continuous Galerkin discretization for anisotropic heat flux in magnetic confinement fusion

Golo A. Wimmer,Ben S. Southworth,Koki Sagiyama,Xian-Zhu Tang

2024-12-17

Abstract:We present a novel spatial discretization for the anisotropic heat conduction equation, aimed at improved accuracy at the high levels of anisotropy seen in a magnetized plasma, for example, for magnetic confinement fusion. The new discretization is based on a mixed formulation, introducing a form of the directional derivative along the magnetic field as an auxiliary variable and discretizing both the temperature and auxiliary fields in a continuous Galerkin (CG) space. Both the temperature and auxiliary variable equations are stabilized using the streamline upwind Petrov-Galerkin (SUPG) method, ensuring a better representation of the directional derivatives and therefore an overall more accurate solution. This approach can be seen as the CG-based version of our previous work (Wimmer, Southworth, Gregory, Tang, 2024), where we considered a mixed discontinuous Galerkin (DG) spatial discretization including DG-upwind stabilization. We prove consistency of the novel discretization, and demonstrate its improved accuracy over existing CG-based methods in test cases relevant to magnetic confinement fusion. This includes a long-run tokamak equilibrium sustainment scenario, demonstrating a 35% and 32% spurious heat loss for existing primal and mixed CG-based formulations versus 4% for our novel SUPG-stabilized discretization.

Numerical Analysis

What problem does this paper attempt to address?

### What problem does this paper attempt to solve?

This paper aims to solve the problem of anisotropic heat flux of plasma in magnetic confinement fusion. Specifically, the author proposes a new spatial discretization method to improve the numerical simulation accuracy under extremely high anisotropic conditions (such as magnetized plasma in magnetic confinement fusion).

#### Main problem description

1. **Extremely high anisotropy**: In the magnetic confinement fusion environment, plasma exhibits extremely anisotropic heat flux, that is, the thermal conductivity along the magnetic field line direction is several orders of magnitude (up to 10 orders of magnitude) higher than that perpendicular to the magnetic field line direction. If the numerical simulation cannot accurately represent this parallel heat flux, it will lead to excessive transverse diffusion, resulting in too short energy confinement time.

2. **Limitations of existing methods**: The existing continuous Galerkin (CG) discretization method has the problem of decreasing accuracy when dealing with high anisotropy, especially when the computational grid cannot be aligned with the magnetic flux surface. This may lead to significant spurious heat loss and affect the accuracy of the simulation results.

#### Solution

To solve the above problems, the author proposes a continuous Galerkin discretization method based on SUPG (Streamline Upwind Petrov - Galerkin) stabilization. This method introduces the directional derivative of temperature along the magnetic field line direction as an auxiliary variable, and performs CG discretization on it and the temperature field, and at the same time uses the SUPG method to stabilize the two directional derivatives. Main features include:

- **Hybrid formula**: Introduces the directional derivative of temperature along the magnetic field line direction as an auxiliary variable.
- **SUPG stabilization**: Applies the SUPG method to the directional derivative equation to ensure better representation of the directional derivative, thereby improving the overall solution accuracy.
- **Consistency proof**: Proves the consistency of the new discretization method and verifies its superior performance through test cases.

#### Experimental results

In a series of numerical tests, the new method shows significant advantages:
- In the long - time tokamak equilibrium maintenance scenario, the existing original and hybrid CG discretization methods produce 35% and 32% spurious heat loss respectively, while the newly proposed SUPG - stabilized discretization method only produces 4% spurious heat loss.
- The new method shows higher accuracy and stability in the MHD model test with actual parameter scales.

In conclusion, this paper solves the accuracy problem of anisotropic heat flux numerical simulation in magnetic confinement fusion by proposing a new SUPG - stabilized continuous Galerkin discretization method, providing a powerful tool for more accurate simulation.