2DESR: a two-dimensional Fourier-space gyrokinetic eigenvalue code for the ion-temperature-gradient modes in tokamaks

The paper introduces 2DESR, a new two-dimensional gyrokinetic eigenvalue solver designed to study ion-temperature-gradient (ITG) modes in tokamaks, which has been validated against existing initial-value codes and demonstrates the coexistence of two ITG mode branches.

The Gist

The Cosmic Soup Stirrer: Understanding the 2DESR Code

Imagine you are trying to study a giant, swirling pot of soup. This soup isn't just any soup; it’s a plasma—a super-hot, electrified "gas" that scientists use to try and create clean energy through nuclear fusion (the same process that powers the sun).

The problem is that this "soup" is incredibly turbulent. It has tiny, invisible whirlpools called ITG modes (Ion-Temperature-Gradient modes). If these whirlpools get too big and wild, they "leak" the heat out of the pot, making it impossible to keep the fusion reaction going. To build a working fusion reactor, we need to predict exactly how these whirlpools behave.

Here is how this paper explains a new tool designed to do just that.

---

1. The Problem: The "Too Much Information" Trap

To understand these whirlpools, scientists use computer simulations. There are two main ways they usually do this:

• The Movie Method (Initial-Value Codes): This is like filming a long movie of the soup. You press "play" and watch how the bubbles move over time. It’s very accurate, but it takes a massive amount of computer power and time to watch the whole movie just to see one tiny detail.
• The Snapshot Method (1D Eigenvalue Codes): This is like taking a quick photo to find the "vibe" of the whirlpool. It’s much faster, but it’s like looking at a flat photograph of a 3D object—you lose the depth and the true shape of the swirl.

The Gap: We needed something that has the speed of a snapshot but the depth of a movie.

---

2. The Solution: Enter "2DESR"

The researchers created a new piece of software called 2DESR.

Think of 2DESR as a High-Tech 3D Scanner. Instead of filming a long movie (slow) or looking at a flat photo (inaccurate), 2DESR uses advanced math to "scan" the soup. It looks at the two most important dimensions: how far out the swirl goes (radial) and how it twists around the center (poloidal).

By focusing on these two dimensions in a special mathematical "space," the code can find the "DNA" of the whirlpools—their frequency and how fast they grow—without having to simulate every single second of the soup's life.

---

3. The Discovery: The "Double Trouble" Whirlpools

When the scientists ran their new 2DESR scanner on a simulated plasma, they found something very interesting.

Previously, other tools only saw one main type of whirlpool. But 2DESR revealed that there are actually two different branches of these ITG whirlpools coexisting at the same time!

It’s like discovering that a storm isn't just one big wind, but actually two different types of gusts overlapping each other. One branch might be small and localized, while the other is wider and more aggressive. If we only look for one, we might miss the one that actually causes the heat to leak!

---

4. Why Does This Matter?

If we want to build a "star in a bottle" (a fusion reactor) to provide endless clean energy, we have to be master chefs of this plasma soup.

2DESR is like a new, high-speed diagnostic tool for the kitchen. It allows scientists to:

1. See the full shape of the turbulence (not just a flat version).
2. Work much faster than the "movie" method.
3. Identify hidden dangers (like that second branch of whirlpools) before they ruin the recipe.

In short: This paper introduces a faster, smarter way to map the invisible storms inside fusion reactors, helping us get one step closer to mastering the power of the stars.

Technical Summary

Technical Summary: 2DESR – A Two-Dimensional Fourier-Space Gyrokinetic Eigenvalue Code

1. Problem Statement

In tokamak plasma physics, Ion-Temperature-Gradient (ITG) modes are a primary driver of anomalous transport, which degrades plasma confinement. To study these instabilities, researchers typically use two types of numerical tools:

• Initial-value codes: Powerful for nonlinear turbulence but computationally expensive for linear stability analysis as they require time evolution.
• Eigenvalue codes: Efficient for computing growth rates ($\gamma$) and real frequencies ($\omega_r$) in the linear regime.

Existing eigenvalue solvers often face a trade-off: 1D solvers (using ballooning theory) are efficient but cannot capture the global 2D radial structures of modes, which are essential for understanding zonal flow evolution and radial correlation lengths. Conversely, 2D solvers in ballooning space often rely on the "well-circulating particle" assumption, which simplifies the physics by neglecting certain kinetic effects. There is a need for a 2D solver that operates in poloidal Fourier space to retain full kinetic effects (including trapped ions) while capturing global mode structures.

2. Methodology

The authors developed 2DESR (2D eigenvalue solver in real space), a code designed to solve the 2D gyrokinetic eigenvalue problem in the poloidal Fourier space.

• Mathematical Framework: The code is derived from the Vlasov–Poisson system. It considers electrostatic perturbations of the ion distribution function with an adiabatic approximation for electrons.
• Coordinate System: To balance computational efficiency with physical accuracy, the authors employ a specific set of coordinates:
  • $(z, m)$ space: Where $z = nq(r) - m$ (radial/rational surface coordinate) and $m$ is the poloidal harmonic number. This allows for radial localization of modes.
  • $(v_\parallel, \mu)$ space: Parallel velocity and magnetic moment. Unlike previous 1D codes that used energy coordinates, this choice allows for a uniform treatment of both passing and trapped ions within the poloidal Fourier decomposition.
• Numerical Implementation:
  • The 2D eigenvalue problem is transformed into an algebraic eigenvalue equation $[M(\omega)] [\delta\hat{\phi}] = 0$.
  • Discretization: The Vlasov equation is discretized using the finite difference method, and the velocity space integral is computed using Gauss–Laguerre quadrature.
  • Solver: The code uses a sparse linear solver (PARDISO) to handle the large matrices and Newton’s method to find the eigenvalues $\omega$. It is parallelized along the $\mu$ direction using MPI.

3. Key Contributions

• Full Kinetic Treatment: Unlike many ballooning-space 2D codes, 2DESR retains full kinetic effects for ions, including the complex dynamics of trapped particles.
• 2D Fourier-Space Approach: By solving in poloidal Fourier space rather than ballooning space, the code avoids the high-$n$ (toroidal mode number) asymptotic approximations required by ballooning theory, making it applicable to both low-$n$ and high-$n$ regimes.
• Global Structure Capture: The code provides the full 2D radial envelope and the coupling between neighboring poloidal harmonics, which is critical for studying the interaction between ITG modes and zonal flows.

4. Results

The code was benchmarked using the "Cyclone" test case against established initial-value codes (GENE and NLT):

• Eigenvalue Accuracy: 2DESR successfully identified two distinct branches of ITG modes. It showed excellent agreement with GENE for the first branch and with NLT for the most unstable branch.
• Mode Structure: The 2D eigenmode structures (electrostatic potential) produced by 2DESR match the "tilted ballooning mode" profiles seen in NLT.
• Resolution of Discrepancies: The study explains a previously observed mismatch in real frequencies between the GENE and GKW codes, attributing it to the presence of a strong sub-dominant mode that 2DESR is capable of resolving.
• Computational Efficiency: For a standard $n=20$ simulation, the code solved a matrix of dimension $\approx 57,000$ in approximately 140 seconds using 64 cores, demonstrating its readiness for rapid experimental analysis.

5. Significance

The development of 2DESR provides the fusion community with a robust tool for the linear analysis of ITG instabilities. Its ability to capture the radial localization and poloidal harmonic coupling of modes makes it indispensable for studying the mechanisms of anomalous transport and the evolution of zonal flows. Furthermore, its applicability to low-$n$ modes fills a gap left by traditional ballooning-based eigenvalue solvers.