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 modelling 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", stated PPPL physicist Ammar Hakim, lead author of aPhysics 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, co-authored 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", Ammar Hakim said.
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 Programme 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 BritishJournal of Plasma Physicshas published and highlighted his report as a featured article.
Noah 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.
Noah 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 stated. "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."
What makes this finding unique is that previous gyrokinetic codes have simulated SOL blobs but assumed that the field lines were rigid, Noah 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", Noah 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", Ammar Hakim stated. "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.
A video of first computer simulations of kinetic plasma turbulence near the edge of fusion devices.