Predicting core transport in ITER baseline discharges with neon injections

This study utilizes integrated modeling to identify a restricted compatibility window for neon-seeded ITER baseline discharges, determining that core transport predictions align with divertor protection targets only when the effective charge is approximately 1.6–1.75 and auxiliary heating is maintained between 75% and 100% of nominal levels.

The Gist

Imagine you are trying to bake the perfect, massive cake in a giant oven (the ITER fusion reactor). Your goal is to make the cake so hot and dense that it generates its own heat forever (a "burning plasma"), creating enough energy to power a city.

However, you have two major problems:

1. The Oven Door is Melting: The heat escaping the oven is so intense it could melt the walls and the exhaust pipe (the divertor).
2. The Recipe is Delicate: If you add too much "flour" (impurities) to cool the oven down, the cake stops rising and loses its power. If you add too little, the oven door melts.

This paper is like a master chef's new recipe guide. It tries to figure out the exact balance of ingredients and oven settings needed to keep the cake rising while keeping the oven door safe.

The Ingredients: Neon and "Effective Charge"

To stop the oven door from melting, the scientists plan to spray a special gas called Neon into the edge of the oven. Think of Neon as a "heat shield" or a "radiator." It glows and radiates heat away before it hits the walls, protecting the oven.

But here's the catch: Neon is heavy. If you put too much in, it dilutes the fuel (the hydrogen isotopes), making the cake weak and unable to generate enough power. If you put too little, the oven door burns.

The scientists measured this balance using a number called $Z_{eff}$ (Effective Charge).

• Low $Z_{eff}$: Not enough neon. The oven is too hot; the door might melt.
• High $Z_{eff}$: Too much neon. The cake is too "watery" and won't cook properly.

The Simulation: A Digital Test Kitchen

Since we can't actually build a real fusion reactor in a garage to test this, the team used a super-computer simulation called OMFIT/STEP.

Think of this as a virtual reality cooking show. They fed the computer:

• The shape of the oven.
• The amount of fuel.
• The amount of neon spray.
• The heat coming from the oven's own "self-cooking" (fusion).

They then ran the simulation to see what happens to the heat flow. They asked: "If we spray this much neon, how much heat is actually hitting the oven door?"

The Big Discovery: The "Goldilocks Zone"

The computer ran thousands of scenarios, changing the amount of neon and the amount of extra heat (auxiliary heating) we put in. They found a very narrow "sweet spot" where everything works:

1. The Neon Amount ($Z_{eff} \approx 1.6 - 1.75$): You need just the right amount of neon.

   • If you have 1.75 (a bit more neon), the oven door is safe, but the cake is slightly weaker.
   • If you have 1.6 (a bit less neon), the cake is stronger, but the door is hotter.
   • The Magic: At these specific levels, the heat hitting the door is about 100 Megawatts. This is the exact limit the engineers say the oven door can handle without melting.

2. The Heating Adjustment:

   • If you use the "standard" amount of extra heat, you need that slightly higher neon level (1.75).
   • If you turn down the extra heat to about 75% of normal, you can get away with slightly less neon (1.6).
   • Analogy: It's like turning down the stove burner so the pot doesn't boil over, allowing you to use less water (neon) to cool it down.

The "Spin" Factor

The scientists also wondered: What if the cake spins? In fusion reactors, the plasma spins. Does this spin help or hurt?

• They tested spinning the plasma faster and slower.
• Result: It didn't matter much! Whether the plasma spun fast or slow, the heat hitting the door only changed by about 20%. The "Neon vs. Heat" balance was the real boss, not the spin.

The "Ghost" Particles

They also checked if invisible gas particles (neutrals) floating inside the oven would cause extra radiation.

• Result: No. The oven is so hot and the gas so thin inside that these particles are practically ghosts. They don't affect the recipe.

The Conclusion: A New Roadmap

This paper gives the ITER team a clear instruction manual for their first years of operation:

• Don't guess: Don't just spray neon randomly.
• Aim for the window: Keep the impurity level ($Z_{eff}$) between 1.6 and 1.75.
• Adjust the heat: If the oven gets too hot, turn down the auxiliary heaters to about 75% instead of just adding more neon.

By following this "Goldilocks" recipe, ITER can hopefully achieve its goal of creating a self-sustaining, super-hot plasma that generates massive energy without melting its own exhaust pipes. It's the difference between a successful, world-changing cake and a burnt, ruined mess.

Technical Summary

1. Problem Statement

Achieving the ITER mission goal of a sustained burning plasma with a fusion gain $Q \ge 10$ requires a delicate balance between core performance (high temperature and confinement) and divertor protection (limiting heat flux to engineering limits).

• The Conflict: To protect the divertor, ITER will likely use impurity seeding (specifically neon) to radiate power in the Scrape-Off Layer (SOL) and divertor. However, increasing impurity concentration ($Z_{eff}$) dilutes the fuel, increases collisionality, and alters turbulence, potentially degrading core confinement and fusion reactivity.
• The Gap: Previous studies often treated the core and edge separately. Divertor models (SOLPS-ITER) prescribed core profiles to satisfy heat flux limits, while core transport models (TGYRO) did not account for the specific impurity constraints required by the divertor. There was a lack of self-consistent predictions linking impurity seeding levels directly to core power balance and fusion performance.

2. Methodology

The authors developed a fully integrated, predictive modeling framework using the OMFIT STEP workflow to close the loop between core transport and divertor power exhaust.

• Workflow Integration:
  • Equilibrium: Used CORSICA for the 15 MA H-mode inductive baseline equilibrium.
  • Transport Solver: TGYRO (coupled with TGLF for turbulent transport and NEO for neoclassical transport) was used to predict stationary temperature and density profiles self-consistently.
  • Divertor Constraint: Boundary conditions (pedestal values, neutral fueling, and impurity profiles near the separatrix) were constrained by existing SOLPS-ITER solutions for neon-seeded, no-ELM scenarios.
• Simulation Strategy:
  • $Z_{eff}$ Scan: A controlled scan of effective charge ($1.5 \le Z_{eff} \le 2.5$) was performed by varying neon and helium concentrations while keeping the D-T fuel ratio fixed.
  • Auxiliary Power Scan: At a fixed $Z_{eff}$, auxiliary heating power ($P_{aux}$) was varied from 60% to 120% of the nominal 51 MW (33 MW NBI + 18 MW EC).
  • Rotation Sensitivity: Toroidal rotation amplitudes ($0$ to $30$ krad/s) were tested to assess the impact of intrinsic rotation shear on transport.
  • Charge Exchange (CX) Check: The AURORA module was used to verify if charge-exchange radiation significantly impacts the core power balance.

3. Key Contributions

• First Self-Consistent Map: This study provides the first quantitative link between impurity seeding requirements (derived from divertor protection) and core transport predictions for the ITER baseline.
• Identification of a Compatibility Window: The work defines a narrow operational window where core transport predictions align with divertor heat flux limits without compromising fusion gain.
• Validation of SOLPS Constraints: It demonstrates that SOLPS-ITER predictions of $\sim 100$ MW crossing the separatrix are only achievable if the core $Z_{eff}$ and auxiliary heating fall within specific ranges.

4. Key Results

A. Sensitivity to Effective Charge ($Z_{eff}$)

• Power Flow Variation: The power crossing the separatrix ($P_{sep}$) is highly sensitive to impurity content. As $Z_{eff}$ increases from 1.5 to 2.5, $P_{sep}$ drops by more than a factor of 1.7 (from $\sim 120$ MW to $\sim 70$ MW).
• The "Sweet Spot": A value of $Z_{eff} \approx 1.75$ yields $P_{sep} \approx 100$ MW, which matches the target condition for the SOLPS-ITER neon-seeded divertor solution.
• Trade-off: While higher $Z_{eff}$ aids divertor cooling, it dilutes the plasma and reduces alpha heating, threatening the $Q \ge 10$ goal.

B. Sensitivity to Auxiliary Heating ($P_{aux}$)

• Alternative Pathway: For a fixed $Z_{eff} = 1.6$ (consistent with SOLPS edge profiles), reducing auxiliary heating to $\sim 75\%$ of nominal ($f_{Paux} \approx 0.75$) also results in $P_{sep} \approx 100$ MW.
• Implication: Operators can satisfy divertor limits either by optimizing impurity seeding ($Z_{eff} \approx 1.75$) or by slightly reducing external heating power, though the latter reduces fusion gain.

C. Sensitivity to Rotation and Charge Exchange

• Rotation: Variations in toroidal rotation amplitude (representing low-torque scenarios) modify $P_{sep}$ by $\lesssim 20\%$. While rotation shear stabilizes turbulence and slightly increases core temperatures, its impact is secondary compared to $Z_{eff}$ and $P_{aux}$.
• Charge Exchange: Modeling confirmed that charge-exchange radiation inside the separatrix is dynamically negligible ($< 1\%$ of total core radiation) due to extremely low neutral densities in the ITER core.

D. The Compatibility Window

The study identifies a restricted operating region where core transport and divertor protection are simultaneously satisfied:

• $Z_{eff} \approx 1.6 - 1.75$
• $0.75 \lesssim f_{Paux} \le 1.0$ (where $f_{Paux}$ is the fraction of nominal auxiliary heating).

5. Significance and Outlook

• Actionable Guidance: The results provide concrete targets for early ITER operation regarding impurity control (neon seeding rates) and heating schedules to ensure the machine does not overheat the divertor or lose burning plasma capability.
• Whole-Device Optimization: This framework moves beyond isolated component modeling, demonstrating that "whole-device" consistency places non-trivial constraints on viable operating points.
• Future Work: The authors note that while the current study focuses on neon, future iterations must include tungsten (the actual wall material) transport and accumulation, as well as more detailed coupling of pedestal physics and momentum transport to refine the compatibility window.

In summary, this paper establishes that achieving ITER's dual goals of high fusion gain and safe divertor operation requires a precise, self-consistent balance of impurity concentration and heating power, with a narrow margin for error that must be managed through advanced predictive modeling.