Discovery about the edge of fusion plasma could help realize fusion power

First-of-a-kind finding

Among the first-of-a-kind findings has been the discovery that accounting for the turbulent fluctuations in the magnetic fields that confine the plasma that fuels fusion reactions can significantly reduce the turbulent particle flux near the plasma edge. Computer simulations show that the net particle flux can go down by as much as 30 percent, despite the fact that the average magnitude of turbulent particle density fluctuation goes up by 60 percent -- indicating that even though the turbulent density fluctuations are more virulent, they are moving particles out of the device less effectively.

Researchers have developed a specialized code called "Gkeyll" -- pronounced just like "Jekyll" in Robert Louis Stevenson's "The Strange Case of Dr. Jekyll and Mr. Hyde" -- that makes these simulations feasible. The mathematical code, a form of modeling called "gyrokinetics," simulates the orbiting of plasma particles around the magnetic field lines at the edge of a fusion plasma.

"Our recent paper summarizes the Gkeyll group's efforts in the area of gyrokinetic simulation," said PPPL physicist Ammar Hakim, lead author of a Physics of Plasmas paper that provides an overview of the group's achievements, based on an invited talk he gave at the American Physical Society's Division of Plasma Physics (APS-DPP) conference last Fall. The research, coauthored by scientists from six institutions, adapts a state-of-the-art algorithm to the gyrokinetic system to develop the "key numerical breakthroughs needed to provide accurate simulations," Hakim said.

Worldwide effort

Such breakthroughs are part of the worldwide effort to grasp the science behind the production of fusion reactions on Earth. Fusion reactions combine light elements in the form of plasma -- the hot, charged state of matter composed of free electrons and atomic nuclei that makes up 99 percent of the visible universe -- to generate massive amounts of energy that could provide a virtually inexhaustible supply of power to generate electricity for humanity.

Noah Mandell, a graduate student in the Princeton University Program in Plasma Physics, built on the team's work to develop the first gyrokinetic code able to handle magnetic fluctuations in what is called the plasma scrape-off layer (SOL) at the edge of tokamak plasmas. The British Journal of Plasma Physics has published and highlighted his report as a featured article.

Mandell explores how blob-like plasma turbulence bends magnetic field lines, leading to the dynamics of "dancing field lines." He finds that field lines usually move smoothly but when dancing can abruptly reconfigure into reconnection events that cause them to converge and violently snap apart.

Mandell's findings are best described as "proof-of-concept" with regard to the magnetic fluctuations, he said. "We know there are more physical effects that need to be added to the code for detailed comparisons with experiments, but already the simulations are showing interesting properties near the plasma edge," he said. "The ability to handle bending of the magnetic field lines will also be essential for future simulations of edge localized modes (ELMs), which we would like to do better to understand the bursts of heat they cause that must be controlled to prevent tokamak damage."

Very challenging

What makes this finding unique is that previous gyrokinetic codes have simulated SOL blobs but assumed that the field lines were rigid, Mandell noted. Extending a gyrokinetic code to calculate the movement of magnetic fields lines is computationally very challenging, requiring special algorithms to ensure that two large terms balance each other to an accuracy of better than 1 part in a million.

Moreover, while codes that model turbulence in the core of the tokamak can include magnetic fluctuations, such codes cannot simulate the SOL region. "The SOL requires specialized codes like Gkeyll that can handle much larger plasma fluctuations and interactions with the walls of the reactor," Mandell said.

Future steps for the Gkeyll group will include investigating the precise physical mechanism that affects the dynamics of the plasma edge, an effect likely connected to the bending field lines. "This work provides stepping stones that I think are very important," Hakim said. "Without the algorithms that we made, these findings would be very difficult to apply to ITER and other machines."

Support for this work comes from the DOE Office of Science, the DOE Computational Science Graduate Fellowship, and the Partnership for Multiscale Gyrokinetic Turbulence Scientific Discovery through Advanced Computing (SciDAC) project. Computations were performed on the Eddy cluster at PPPL.

PPPL, on Princeton University's Forrestal Campus in Plainsboro, N.J., is devoted to creating new knowledge about the physics of plasmas -- ultra-hot, charged gases -- and to developing practical solutions for the creation of fusion energy.