Showing 1–17 of 17 results for author: Churchill, M

Flexible Stellarator Physics Facility

Authors: F. I. Parra, S. -G. Baek, M. Churchill, D. R. Demers, B. Dudson, N. M. Ferraro, B. Geiger, S. Gerhardt, K. C. Hammond, S. Hudson, R. Jorge, E. Kolemen, D. M. Kriete, S. T. A. Kumar, M. Landreman, C. Lowe, D. A. Maurer, F. Nespoli, N. Pablant, M. J. Pueschel, A. Punjabi, J. A. Schwartz, C. P. S. Swanson, A. M. Wright

Abstract: We propose to build a Flexible Stellarator Physics Facility to explore promising regions of the vast parameter space of disruption-free stellarator solutions for Fusion Pilot Plants (FPPs). We propose to build a Flexible Stellarator Physics Facility to explore promising regions of the vast parameter space of disruption-free stellarator solutions for Fusion Pilot Plants (FPPs). △ Less

Submitted 4 July, 2024; originally announced July 2024.

Comments: White paper submitted to FESAC subcommittee on Facilities, 8 pages

Unraveling Diffusion in Fusion Plasma: A Case Study of In Situ Processing and Particle Sorting

Authors: Junmin Gu, Paul Lin, Kesheng Wu, Seung-Hoe Ku, C. S. Chang, R. Michael Churchill, Jong Choi, Norbert Podhorszki, Scott Klasky

Abstract: This work starts an in situ processing capability to study a certain diffusion process in magnetic confinement fusion. This diffusion process involves plasma particles that are likely to escape confinement. Such particles carry a significant amount of energy from the burning plasma inside the tokamak to the diverter and damaging the diverter plate. This study requires in situ processing because of… ▽ More This work starts an in situ processing capability to study a certain diffusion process in magnetic confinement fusion. This diffusion process involves plasma particles that are likely to escape confinement. Such particles carry a significant amount of energy from the burning plasma inside the tokamak to the diverter and damaging the diverter plate. This study requires in situ processing because of the fast changing nature of the particle diffusion process. However, the in situ processing approach is challenging because the amount of data to be retained for the diffusion calculations increases over time, unlike in other in situ processing cases where the amount of data to be processed is constant over time. Here we report our preliminary efforts to control the memory usage while ensuring the necessary analysis tasks are completed in a timely manner. Compared with an earlier naive attempt to directly computing the same diffusion displacements in the simulation code, this in situ version reduces the memory usage from particle information by nearly 60% and computation time by about 20%. △ Less

Submitted 2 November, 2023; originally announced November 2023.

2022 Review of Data-Driven Plasma Science

Authors: Rushil Anirudh, Rick Archibald, M. Salman Asif, Markus M. Becker, Sadruddin Benkadda, Peer-Timo Bremer, Rick H. S. Budé, C. S. Chang, Lei Chen, R. M. Churchill, Jonathan Citrin, Jim A Gaffney, Ana Gainaru, Walter Gekelman, Tom Gibbs, Satoshi Hamaguchi, Christian Hill, Kelli Humbird, Sören Jalas, Satoru Kawaguchi, Gon-Ho Kim, Manuel Kirchen, Scott Klasky, John L. Kline, Karl Krushelnick , et al. (38 additional authors not shown)

Abstract: Data science and technology offer transformative tools and methods to science. This review article highlights latest development and progress in the interdisciplinary field of data-driven plasma science (DDPS). A large amount of data and machine learning algorithms go hand in hand. Most plasma data, whether experimental, observational or computational, are generated or collected by machines today.… ▽ More Data science and technology offer transformative tools and methods to science. This review article highlights latest development and progress in the interdisciplinary field of data-driven plasma science (DDPS). A large amount of data and machine learning algorithms go hand in hand. Most plasma data, whether experimental, observational or computational, are generated or collected by machines today. It is now becoming impractical for humans to analyze all the data manually. Therefore, it is imperative to train machines to analyze and interpret (eventually) such data as intelligently as humans but far more efficiently in quantity. Despite the recent impressive progress in applications of data science to plasma science and technology, the emerging field of DDPS is still in its infancy. Fueled by some of the most challenging problems such as fusion energy, plasma processing of materials, and fundamental understanding of the universe through observable plasma phenomena, it is expected that DDPS continues to benefit significantly from the interdisciplinary marriage between plasma science and data science into the foreseeable future. △ Less

Submitted 31 May, 2022; originally announced May 2022.

Comments: 112 pages (including 700+ references), 44 figures, submitted to IEEE Transactions on Plasma Science as a part of the IEEE Golden Anniversary Special Issue

Report number: Los Alamos Report number LA-UR-22-24834

Journal ref: IEEE Transactions on Plasma Science 51, 1750 - 1838 (2023)

Electromagnetic total-f algorithm for gyrokinetic particle-in-cell simulations of boundary plasma in XGC

Authors: Robert Hager, Seung-Hoe Ku, A. Y. Sharma, C. S. Chang, R. M. Churchill

Abstract: The simplified delta-f mixed-variable/pull-back electromagnetic simulation algorithm implemented in XGC for core plasma simulations by M. Cole et al. [Phys. Plasmas 28, 034501 (2021)] has been generalized to a total-f electromagnetic algorithm that can include, for the first time, the boundary plasma in diverted magnetic geometry with neutral particle recycling, turbulence and neoclassical physics… ▽ More The simplified delta-f mixed-variable/pull-back electromagnetic simulation algorithm implemented in XGC for core plasma simulations by M. Cole et al. [Phys. Plasmas 28, 034501 (2021)] has been generalized to a total-f electromagnetic algorithm that can include, for the first time, the boundary plasma in diverted magnetic geometry with neutral particle recycling, turbulence and neoclassical physics. The delta-f mixed-variable/pull-back electromagnetic is based on the pioneering work by Kleiber and Mischenko et al. [R. Kleiber et al., Phys. Plasmas 23, 032501 (2016); A. Mishchenko et al., Comput. Phys. Commun. 238, 194 (2019)]. An electromagnetic demonstration simulation is performed in a DIII-D-like, H-mode boundary plasma, including a corresponding comparative electrostatic simulation, which confirms that the electromagnetic simulation is necessary for a higher fidelity understanding of the electron particle and heat transport even at the low-beta pedestal foot in the vicinity of the magnetic separatrix. △ Less

Submitted 26 April, 2022; v1 submitted 12 February, 2022; originally announced February 2022.

Machine learning accelerated particle-in-cell plasma simulations

Authors: R. Kube, R. M. Churchill, B. Sturdevant

Abstract: Particle-In-Cell (PIC) methods are frequently used for kinetic, high-fidelity simulations of plasmas. Implicit formulations of PIC algorithms feature strong conservation properties, up to numerical round-off errors, and are not subject to time-step limitations which make them an attractive candidate to use in simulations fusion plasmas. Currently they remain prohibitively expensive for high-fideli… ▽ More Particle-In-Cell (PIC) methods are frequently used for kinetic, high-fidelity simulations of plasmas. Implicit formulations of PIC algorithms feature strong conservation properties, up to numerical round-off errors, and are not subject to time-step limitations which make them an attractive candidate to use in simulations fusion plasmas. Currently they remain prohibitively expensive for high-fidelity simulation of macroscopic plasmas. We investigate how amortized solvers can be incorporated with PIC methods for simulations of plasmas. Incorporated into the amortized solver, a neural network predicts a vector space that entails an approximate solution of the PIC system. The network uses only fluid moments and the electric field as input and its output is used to augment the vector space of an iterative linear solver. We find that this approach reduces the average number of required solver iterations by about 25% when simulating electron plasma oscillations. This novel approach may allow to accelerate implicit PIC simulations while retaining all conservation laws and may also be appropriate for multi-scale systems. △ Less

Submitted 24 October, 2021; originally announced October 2021.

Near real-time streaming analysis of big fusion data

Authors: Ralph Kube, R. Michael Churchill, CS Chang, Jong Choi, Jason Wang, Scott Klasky, Laurie Stephey, Minjun Choi, Eli Dart

Abstract: While experiments on fusion plasmas produce high-dimensional data time series with ever increasing magnitude and velocity, data analysis has been lagging behind this development. For example, many data analysis tasks are often performed in a manual, ad-hoc manner some time after an experiment. In this article we introduce the DELTA framework that facilitates near real-time streaming analysis of bi… ▽ More While experiments on fusion plasmas produce high-dimensional data time series with ever increasing magnitude and velocity, data analysis has been lagging behind this development. For example, many data analysis tasks are often performed in a manual, ad-hoc manner some time after an experiment. In this article we introduce the DELTA framework that facilitates near real-time streaming analysis of big and fast fusion data. By streaming measurement data from fusion experiments to a high-performance compute center, DELTA allows to perform demanding data analysis tasks in between plasma pulses. This article describe the modular and expandable software architecture of DELTA and presents performance benchmarks of its individual components as well as of entire workflows. Our focus is on the streaming analysis of ECEi data measured at KSTAR on NERSCs supercomputers and we routinely achieve data transfer rates of about 500 Megabyte per second. We show that a demanding turbulence analysis workload can be distributed among multiple GPUs and executes in under 5 minutes. We further discuss how DELTA uses modern database systems and container orchestration services to provide web-based real-time data visualization. For the case of ECEi data we demonstrate how data visualizations can be augmented with outputs from machine learning models. By providing session leaders and physics operators results of higher order data analysis using live visualization they may monitor the evolution of a long-pulse discharge in near real-time and may make more informed decision on how to configure the machine for the next shot. △ Less

Submitted 19 August, 2021; originally announced August 2021.

Constructing a new predictive scaling formula for ITER's divertor heat-load width informed by a simulation-anchored machine learning

Authors: C. S. Chang, S. Ku, R. Hager, R. M. Churchill, J. Hughes, F. Köchl, A. Loarte, V. Parail, R. Pitts

Abstract: Understanding and predicting divertor heat-load width $λ_q$ is a critically important problem for an easier and more robust operation of ITER with high fusion gain. Previous predictive simulation data for $λ_q$ using the extreme-scale edge gyrokinetic code XGC1 in the electrostatic limit under attached divertor plasma conditions in three major US tokamaks [C.S. Chang et al., Nucl. Fusion 57, 11602… ▽ More Understanding and predicting divertor heat-load width $λ_q$ is a critically important problem for an easier and more robust operation of ITER with high fusion gain. Previous predictive simulation data for $λ_q$ using the extreme-scale edge gyrokinetic code XGC1 in the electrostatic limit under attached divertor plasma conditions in three major US tokamaks [C.S. Chang et al., Nucl. Fusion 57, 116023 (2017)] reproduced the Eich and Goldston attached-divertor formula results [formula #14 in T. Eich et al., Nucl. Fusion 53, 093031 (2013); R.J. Goldston, Nucl. Fusion 52, 013009 (2012)], and furthermore predicted over six times wider $λ_q$ than the maximal Eich and Goldston formula predictions on a full-power (Q = 10) scenario ITER plasma. After adding data from further predictive simulations on a highest current JET and highest-current Alcator C-Mod, a machine learning program is used to identify a new scaling formula for $λ_q$ as a simple modification to the Eich formula #14, which reproduces the Eich scaling formula for the present tokamaks and which embraces the wide $λ_q^X{GC}$ for the full-current Q = 10 ITER plasma. The new formula is then successfully tested on three more ITER plasmas: two corresponding to long burning scenarios with Q = 5 and one at low plasma current to be explored in the initial phases of ITER operation. The new physics that gives rise to the wider $λq_^{XGC} is identified to be the weakly-collisional, trapped-electron-mode turbulence across the magnetic separatrix, which is known to be an efficient transporter of the electron heat and mass. Electromagnetic turbulence and high-collisionality effects on the new formula are the next study topics for XGC1. △ Less

Submitted 5 January, 2021; v1 submitted 19 December, 2020; originally announced December 2020.

Comments: 23 pages, 16 figures, submitted to Physics of Plasmas

Training neural networks under physical constraints using a stochastic augmented Lagrangian approach

Authors: Alp Dener, Marco Andres Miller, Randy Michael Churchill, Todd Munson, Choong-Seock Chang

Abstract: We investigate the physics-constrained training of an encoder-decoder neural network for approximating the Fokker-Planck-Landau collision operator in the 5-dimensional kinetic fusion simulation in XGC. To train this network, we propose a stochastic augmented Lagrangian approach that utilizes pyTorch's native stochastic gradient descent method to solve the inner unconstrained minimization subproble… ▽ More We investigate the physics-constrained training of an encoder-decoder neural network for approximating the Fokker-Planck-Landau collision operator in the 5-dimensional kinetic fusion simulation in XGC. To train this network, we propose a stochastic augmented Lagrangian approach that utilizes pyTorch's native stochastic gradient descent method to solve the inner unconstrained minimization subproblem, paired with a heuristic update for the penalty factor and Lagrange multipliers in the outer augmented Lagrangian loop. Our training results for a single ion species case, with self-collisions and collision against electrons, show that the proposed stochastic augmented Lagrangian approach can achieve higher model prediction accuracy than training with a fixed penalty method for our application problem, with the accuracy high enough for practical applications in kinetic simulations. △ Less

Submitted 15 September, 2020; originally announced September 2020.

Encoder-decoder neural network for solving the nonlinear Fokker-Planck-Landau collision operator in XGC

Authors: M. A. Miller, R. M. Churchill, A. Dener, C. S. Chang, T. Munson, R. Hager

Abstract: An encoder-decoder neural network has been used to examine the possibility for acceleration of a partial integro-differential equation, the Fokker-Planck-Landau collision operator. This is part of the governing equation in the massively parallel particle-in-cell code, XGC, which is used to study turbulence in fusion energy devices. The neural network emphasizes physics-inspired learning, where it… ▽ More An encoder-decoder neural network has been used to examine the possibility for acceleration of a partial integro-differential equation, the Fokker-Planck-Landau collision operator. This is part of the governing equation in the massively parallel particle-in-cell code, XGC, which is used to study turbulence in fusion energy devices. The neural network emphasizes physics-inspired learning, where it is taught to respect physical conservation constraints of the collision operator by including them in the training loss, along with the L2 loss. In particular, network architectures used for the computer vision task of semantic segmentation have been used for training. A penalization method is used to enforce the "soft" constraints of the system and integrate error in the conservation properties into the loss function. During training, quantities representing the density, momentum, and energy for all species of the system is calculated at each configuration vertex, mirroring the procedure in XGC. This simple training has produced a median relative loss, across configuration space, on the order of 10E-04, which is low enough if the error is of random nature, but not if it is of drift nature in timesteps. The run time for the Picard iterative solver of the operator scales as order n squared, where n is the number of plasma species. As the XGC1 code begins to attack problems including a larger number of species, the collision operator will become expensive computationally, making the neural network solver even more important, since the training only scales as n. A wide enough range of collisionality is considered in the training data to ensure the full domain of collision physics is captured. An advanced technique to decrease the losses further will be discussed, which will be subject of a subsequent report. Eventual work will include expansion of the network to include multiple plasma species. △ Less

Submitted 17 December, 2020; v1 submitted 14 September, 2020; originally announced September 2020.

Comparison of edge turbulence characteristics between DIII-D and C-Mod simulations with XGC1

Authors: I. Keramidas Charidakos, J. R. Myra, S. Ku, R. M. Churchill, R. Hager, C. S. Chang, S. Parker

Abstract: The physical processes taking place at the edge region are crucial for the operation of tokamaks as they govern the interaction of hot plasma with the vessel walls. Numerical modeling of the edge with state-of-the-art codes attempts to elucidate interactions between neoclassical drifts, turbulence, poloidal and parallel flows that control the physical set-up of the SOL region. Here, we present pos… ▽ More The physical processes taking place at the edge region are crucial for the operation of tokamaks as they govern the interaction of hot plasma with the vessel walls. Numerical modeling of the edge with state-of-the-art codes attempts to elucidate interactions between neoclassical drifts, turbulence, poloidal and parallel flows that control the physical set-up of the SOL region. Here, we present post-processing analysis of simulations from the gyrokinetic code XGC1, comparing edge turbulence characteristics from a simulation of DIII-D against one of C-Mod. We find that the equilibrium $E \times B$ flux across the separatrix has a similar poloidal pattern in both discharges which can be explained by magnetic drifts and trapped ion excursions. However, collisionality is noted to play a major role in that it prevents local charge accumulations from having global effects in C-Mod. In both cases, turbulent electron heat flux is higher than the ion one. This seems to be a universal characteristic of the tokamak edge. We identify turbulent frequencies and growth rates of the dominant mode in both simulations. In C-Mod, these numbers point to the presence of a drift wave. In DIII-D, linear simulations with Gene reveal a trapped electron mode. Furthermore, we present the amplitude and size distributions of the blobs from both simulations. Amplitude distributions are in qualitative agreement with experimental observations while size distributions are consistent with the fact that most blobs are not connecting to the divertor plates and suggest that they are generated by the shearing of the turbulent modes. △ Less

Submitted 25 March, 2020; originally announced March 2020.

Comments: 14 pages, 11 figures

Deep convolutional neural networks for multi-scale time-series classification and application to disruption prediction in fusion devices

Authors: R. M. Churchill, the DIII-D team

Abstract: The multi-scale, mutli-physics nature of fusion plasmas makes predicting plasma events challenging. Recent advances in deep convolutional neural network architectures (CNN) utilizing dilated convolutions enable accurate predictions on sequences which have long-range, multi-scale characteristics, such as the time-series generated by diagnostic instruments observing fusion plasmas. Here we apply thi… ▽ More The multi-scale, mutli-physics nature of fusion plasmas makes predicting plasma events challenging. Recent advances in deep convolutional neural network architectures (CNN) utilizing dilated convolutions enable accurate predictions on sequences which have long-range, multi-scale characteristics, such as the time-series generated by diagnostic instruments observing fusion plasmas. Here we apply this neural network architecture to the popular problem of disruption prediction in fusion tokamaks, utilizing raw data from a single diagnostic, the Electron Cyclotron Emission imaging (ECEi) diagnostic from the DIII-D tokamak. ECEi measures a fundamental plasma quantity (electron temperature) with high temporal resolution over the entire plasma discharge, making it sensitive to a number of potential pre-disruptions markers with different temporal and spatial scales. Promising, initial disruption prediction results are obtained training a deep CNN with large receptive field (~30k), achieving an $F_1$-score of ~91% on individual time-slices using only the ECEi data. △ Less

Submitted 21 November, 2019; v1 submitted 31 October, 2019; originally announced November 2019.

Comments: Accepted at the Second Workshop on Machine Learning and the Physical Sciences (NeurIPS 2019), Vancouver, Canada

Analysis of equilibrium and turbulent fluxes across the separatrix in a gyrokinetic simulation

Authors: I. Keramidas Charidakos, J. R. Myra, S. Parker, S. Ku, R. M. Churchill, R. Hager, C. S. Chang

Abstract: The SOL width is a parameter of paramount importance in modern tokamaks as it controls the power density deposited at the divertor plates, critical for plasma-facing material survivability. An understanding of the parameters controlling it has consequently long been sought (Connor et al. 1999 NF 39 2). Prior to Chang et al.(2017 NF 57 11), studies of the tokamak edge have been mostly confined to r… ▽ More The SOL width is a parameter of paramount importance in modern tokamaks as it controls the power density deposited at the divertor plates, critical for plasma-facing material survivability. An understanding of the parameters controlling it has consequently long been sought (Connor et al. 1999 NF 39 2). Prior to Chang et al.(2017 NF 57 11), studies of the tokamak edge have been mostly confined to reduced fluid models and simplified geometries, leaving out important pieces of physics. Here, we analyze the results of a DIII-D simulation performed with the full-f gyrokinetic code XGC1 which includes both turbulence and neoclassical effects in realistic divertor geometry. More specifically, we calculate the particle and heat ExB fluxes along the separatrix, discriminating between equilibrium and turbulent contributions. We find that the density SOL width is impacted almost exclusively by the turbulent electron flux. In this simulation, the level of edge turbulence is regulated by a mechanism we are only beginning to understand: $\nabla B$-drifts and ion X-point losses at the top and bottom of the machine, along with ion banana orbits at the low field side (LFS), result in a complex poloidal potential structure at the separatrix which is the cause of the ExB drift pattern that we observe. Turbulence is being suppressed by the shear flows that this potential generates. At the same time, turbulence, along with increased edge collisionality and electron inertia, can influence the shape of the potential structure by making the electrons non-adiabatic. Moreover, being the only means through which the electrons can lose confinement, it needs to be in a balance with the original direct ion orbit losses to maintain charge neutrality. △ Less

Submitted 26 April, 2018; originally announced April 2018.

Comments: 14 pages, Submitted for publication to POP

Pedestal and edge turbulence characteristics from an XGC1 gyrokinetic simulation

Authors: R. M. Churchill, C. S. Chang, S. Ku

Abstract: Understanding the multi-scale neoclassical and turbulence physics in the edge region (pedestal + scrape-off layer) is required in order to reliably predict performance in future fusion devices. We explore turbulent characteristics in the edge region from a multiscale neoclassical and turbulent XGC1 gyrokinetic simulation in a DIII-D like tokamak geometry, here excluding neutrals and collisions. Fo… ▽ More Understanding the multi-scale neoclassical and turbulence physics in the edge region (pedestal + scrape-off layer) is required in order to reliably predict performance in future fusion devices. We explore turbulent characteristics in the edge region from a multiscale neoclassical and turbulent XGC1 gyrokinetic simulation in a DIII-D like tokamak geometry, here excluding neutrals and collisions. For an H-mode type plasma with steep pedestal, it is found that the electron density fluctuations increase towards the separatrix, and stay high well into the SOL, reaching a maximum value of $δn_e / \bar{n}_e \sim 0.18$. Blobs are observed, born around the magnetic separatrix surface and propagate radially outward with velocities generally less than 1 km/s. Strong poloidal motion of the blobs is also present, near 20 km/s, consistent with $E \times B$ rotation. The electron density fluctuations show a negative skewness in the closed field line pedestal regions, consistent with the presence of "holes", followed by a transition to strong positive skewness across the separatrix and into the SOL. These simulations indicate that not only neoclassical phenomena, but also turbulence, including the blob-generation mechanism, can remain important in the steep H-mode pedestal and SOL. Qualitative comparisons will be made to experimental observations. △ Less

Submitted 8 March, 2017; originally announced March 2017.

Electrostatic gyrokinetic simulation of global tokamak boundary plasma and the generation of nonlinear intermittent turbulence

Authors: S. Ku, R. M. Churchill, C. S. Chang, R. Hager, E. S. Yoon, M. Adams, E. D'Azevedo, P. H. Worley

Abstract: Boundary plasma physics plays an important role in tokamak confinement, but is difficult to simulate in a gyrokinetic code due to the scale-inseparable nonlocal multi-physics in magnetic separatrix and open magnetic field geometry. Neutral particles are also an important part of the boundary plasma physics. In the present paper, noble electrostatic gyrokinetic techniques to simulate the flux-drive… ▽ More Boundary plasma physics plays an important role in tokamak confinement, but is difficult to simulate in a gyrokinetic code due to the scale-inseparable nonlocal multi-physics in magnetic separatrix and open magnetic field geometry. Neutral particles are also an important part of the boundary plasma physics. In the present paper, noble electrostatic gyrokinetic techniques to simulate the flux-driven, low-beta electrostatic boundary plasma is reported. Gyrokinetic ions and drift-kinetic electrons are utilized without scale-separation between the neoclassical and turbulence dynamics. It is found that the nonlinear intermittent turbulence is a natural gyrokinetic phenomenon in the boundary plasma in the vicinity of the magnetic separatrix surface and in the scrape-off layer. △ Less

Submitted 24 January, 2017; v1 submitted 20 January, 2017; originally announced January 2017.

Comments: 14 pages, 5 figures, submitted to Nuclear Fusion

Total Fluid Pressure Imbalance in the Scrape-Off Layer of Tokamak Plasmas

Authors: R. M. Churchill, J. M. Canik, C. S. Chang, R. Hager, A. W. Leonard, R. Maingi, R. Nazikian, D. P. Stotler

Abstract: Simulations using the fully kinetic neoclassical code XGCa were undertaken to explore the impact of kinetic effects on scrape-off layer (SOL) physics in DIII-D H-mode plasmas. XGCa is a total-f, gyrokinetic code which self-consistently calculates the axisymmetric electrostatic potential and plasma dynamics, and includes modules for Monte Carlo neutral transport. Previously presented XGCa results… ▽ More Simulations using the fully kinetic neoclassical code XGCa were undertaken to explore the impact of kinetic effects on scrape-off layer (SOL) physics in DIII-D H-mode plasmas. XGCa is a total-f, gyrokinetic code which self-consistently calculates the axisymmetric electrostatic potential and plasma dynamics, and includes modules for Monte Carlo neutral transport. Previously presented XGCa results showed several noteworthy features, including large variations of ion density and pressure along field lines in the SOL, experimentally relevant levels of SOL parallel ion flow (Mach number~0.5), skewed ion distributions near the sheath entrance leading to subsonic flow there, and elevated sheath potentials [R.M. Churchill, Nucl. Mater. & Energy, submitted]. In this paper, we explore in detail the question of pressure balance in the SOL, as it was observed in the simulation that there was a large deviation from a simple total pressure balance (the sum of ion and electron static pressure plus ion inertia). It will be shown that both the contributions from the ion viscosity (driven by ion temperature anisotropy) and neutral source terms can be substantial, and should be retained in the parallel momentum equation in the SOL, but still falls short of accounting for the observed fluid pressure imbalance in the XGCa simulation results. △ Less

Submitted 19 January, 2017; originally announced January 2017.

A fast low-to-high mode bifurcation dynamics in a tokamak edge plasma gyrokinetic simulation

Authors: C. S. Chang, S. Ku, G. R. Tynan, R. Hager, R. M. Churchill, I. Cziegler, M. Greenwald, A. E. Hubbard, J. W. Hughes

Abstract: Transport barrier formation and its relation to sheared flows in fluids and plasmas are of fundamental interest in various natural and laboratory observations and of critical importance in achieving an economical energy production in a magnetic fusion device. Here we report the first observation of an edge transport barrier formation event in a gyrokinetic simulation carried out in a realistic tok… ▽ More Transport barrier formation and its relation to sheared flows in fluids and plasmas are of fundamental interest in various natural and laboratory observations and of critical importance in achieving an economical energy production in a magnetic fusion device. Here we report the first observation of an edge transport barrier formation event in a gyrokinetic simulation carried out in a realistic tokamak edge geometry. The results show that turbulent Reynolds stress driven sheared ExB flows act in concert with neoclassical orbit loss to quench turbulent transport and form a transport barrier just inside the last closed magnetic flux surface. △ Less

Submitted 8 May, 2017; v1 submitted 19 January, 2017; originally announced January 2017.

Comments: Accepted, Physical Review Letter

Journal ref: Phys. Rev. Lett. 118, 175001 (2017)

Towards Real-Time Detection and Tracking of Spatio-Temporal Features: Blob-Filaments in Fusion Plasma

Authors: Lingfei Wu, Kesheng Wu, Alex Sim, Michael Churchill, Jong Y. Choi, Andreas Stathopoulos, Cs Chang, Scott Klasky

Abstract: A novel algorithm and implementation of real-time identification and tracking of blob-filaments in fusion reactor data is presented. Similar spatio-temporal features are important in many other applications, for example, ignition kernels in combustion and tumor cells in a medical image. This work presents an approach for extracting these features by dividing the overall task into three steps: loca… ▽ More A novel algorithm and implementation of real-time identification and tracking of blob-filaments in fusion reactor data is presented. Similar spatio-temporal features are important in many other applications, for example, ignition kernels in combustion and tumor cells in a medical image. This work presents an approach for extracting these features by dividing the overall task into three steps: local identification of feature cells, grouping feature cells into extended feature, and tracking movement of feature through overlapping in space. Through our extensive work in parallelization, we demonstrate that this approach can effectively make use of a large number of compute nodes to detect and track blob-filaments in real time in fusion plasma. On a set of 30GB fusion simulation data, we observed linear speedup on 1024 processes and completed blob detection in less than three milliseconds using Edison, a Cray XC30 system at NERSC. △ Less

Submitted 2 July, 2016; v1 submitted 13 May, 2015; originally announced May 2015.

Comments: 14 pages, 40 figures