Publications

Mar, 2026

Conference Paper

An implicit-explicit time splitting strategy for the far SOL plasma fluid model with DG-FEM discretization

We consider a far scrape-off layer (SOL) plasma fluid model of ions that is governed by a Braginskiitype model: a one-dimensional, nonlinear system of advection-diffusion equations coupled with a diffusion equation for neutral particles. Our motivation for studying this system arises from the coupling between the edge plasma and radio-frequency (RF) heating, where solving…

Mar, 2026

Conference Paper

MAPS: the MFEM Anisotropic Plasma Solver

Simulating magnetically confined fusion plasmas presents a uniquely challenging problem due to the nonlinear anisotropic heat conduction. We introduce the MAPS (MFEM Anisotropic Plasma Solver) tool, which uses a high-order finite element method to compute transport solutions on unstructured meshes. We show results for a set of three 2-D verification tests, two of which dem…

Mar, 2026

Conference Paper

Analysis of RF Sheath-Driven Tungsten Erosion at RF Antenna in the WEST Tokamak

This study applies the newly developed STRIPE (Simulated Transport of RF Impurity Production and Emission) framework to analyze tungsten (W) erosion at RF antenna structures in the WEST tokamak. STRIPE integrates SolEdge3x for edge plasma backgrounds, COMSOL for 3D RF sheath potentials, RustBCA for sputtering yields, and GITR for impurity transport and ion energy–angle dis…

Feb, 2026

Journal

Towards self-consistent integrated modeling of the tokamak pedestal, scrape-off layer, and divertor using SOLPS-ITER and EPED

A predictive, self-consistent integrated model of the tokamak edge plasma and neutral system has been developed and validated against experimental data from DIII-D, spanning from the divertor target and scrape-off layer to the top of the pedestal. The model uses the SOLPS-ITER fluid boundary plasma and neutral code to constrain the EPED prediction of the pedestal height an…

Jan, 2026

ORNL Report

Description of FY25 Theory and Simulation Performance Target: Development of an integrated modeling framework for fusion reactor design and assessment

Jan, 2026

ORNL Report

Q1 Report for FY25 Theory and Simulation Performance Target: Development of an integrated modeling framework for fusion reactor design and assessment

Jan, 2026

ORNL Report

Q2 Report for FY25 Theory and Simulation Performance Target: Development of an integrated modeling framework for fusion reactor design and assessment

Jan, 2026

ORNL Report

Q3 Report for FY25 Theory and Simulation Performance Target: Development of an integrated modeling framework for fusion reactor design and assessment

Jan, 2026

Journal

Gyrotropy-induced symmetry breaking in RF-specific plasma wall interaction: Mechanisms and mitigation strategies

Ion cyclotron range of frequencies (ICRF) heating is essential for creating plasma in next-generation fusion devices. ICRF antennas often produce localized wall erosion, reducing reliability, material survivability, and overall plasma performance. We show that the wave physics' gyrotropy is among the causes of toroidal and poloidal asymmetries in these erosion spots, and t…

Jan, 2026

ORNL Report

Final Report for FY25 Theory and Simulation Performance Target: Development of an integrated modeling framework for fusion reactor design and assessment

Nov, 2025

Journal

Challenges and approaches to interpretive modeling of boundary plasma and neutral transport in a closed, pumped divertor

An experimental discharge from the DIII-D tokamak is modeled using the SOLPS-ITER code suite and compared against measurements in the pumped and relatively closed upper divertor. Uncertainties of boundary plasma simulations are identified by attempting to match code inputs to experimental conditions, including iteratively solving transport coefficients to match upstream ex…

Sep, 2025

Journal

Divertor characterization and access to dissipative divertor conditions in negative triangularity discharges in DIII-D

Reduced radial transport, short midplane-to-target parallel connection lengths, and a strong effect of cross-field drifts were responsible for the high densities required for detachment in strong negative triangularity (NT) configurations in DIII-D. Dissipative divertor conditions were achieved in NT discharges at different triangularity, injected power, plasma current, an…