Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Scrape-off-layer”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2

Core-edge integrated predictive studies of ST40 and NSTX plasmas with the scrape-off layer box model

The ability to model the interplay between the core and edge of tokamak plasmas is crucial to designing both the plasma operating scenario of a fusion pilot plant and the design of the tokamak itself. Scrape-off-layer (SOL) models that are tailored to integrated scenario modeling need to have fast turn-around time and minimal computational burden to enable wide parameter-space coverage for design scoping. The SOL 0-D Box model is a reduced SOL model based on global power and particle balance that captures the essential physics of SOL transport with little computational cost. The usage of the 0-D Box model in core-edge coupled simulations has been demonstrated in both interpretive and predictive modes on a variety of devices. This paper presents a sensitivity study of the 0-D Box model to the input SOL heat-flux width for an ST40 plasma. This study demonstrates that accurate prediction of this width is crucial to predicting global performance parameters of a plasma scenario, such as energy confinement time and flux consumption. We also present an extension of the Box model to 1-D to allow for parallel variation of plasma parameters along the magnetic field lines. The 1-D Box model is then compared with SOLPS-ITER simulations of an NSTX plasma. Advantages and limitations of the Box model are discussed, and future directions are outlined.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Measurements and modeling of type-I and type-II ELMs heat flux to the DIII-D divertor

Type-I and type-II edge-localized-modes (ELMs) heat flux profiles measured at the DIII-D divertor feature a peak in the vicinity of the strike-point and a plateau in the scrape-off-layer (SOL), which extends to the first wall. The plateau is present in attached and detached divertors and it is found to originate with plasma bursts upstream in the SOL. The integrated ELM heat flux is distributed at ~65% in the peak and ~35% in this plateau. The parallel loss model, currently used at ITER to predict power loads to the walls, is benchmarked using these results in the primary and secondary divertors with unprecedented constraints using experimental input data for ELM size, radial velocity, energy, electron temperature and density, heat flux footprints and number of filaments. The model can reproduce the experimental near-SOL peak within ~20%, but cannot match the SOL plateau. Employing a two-component approach for the ELM radial velocity, as guided by intermittent data, the full radial heat flux profile can be well matched. The ELM-averaged radial velocity at the separatrix, which explains profile widening, increases from ~0.2 km s –1 in attached to ~0.8 km s –1 in detached scenarios, as the ELM filaments' path becomes electrically disconnected from the sheath at the target. The results presented here indicate filaments fragmentation as a possible mechanism for ELM transport to the far-SOL and provide evidence on the beneficial role of detachment to mitigate ELM flux in the divertor far-SOL. However, these findings imply that wall regions far from the strike points in future machines should be designed to withstand significant heat flux, even for small-ELM regimes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Analysis and design of fast flow liquid Li divertor for fusion nuclear science facility (FNSF) using coupled plasma boundary and LM MHD/heat transfer codes *

The SOLPS-ITER code is utilized to analyze the boundary plasma associated with a fast-flow lithium (Li) divertor configuration in the fusion nuclear science facility (FNSF) tokamak and identify operational regimes with acceptable divertor and core conditions. Plasma transport from the SOLPS-ITER code has been coupled with a liquid metal (LM) MHD/heat transfer code to model a Li open-surface divertor design and assess its impact on the scrape-off-layer (SOL) and core plasma performance. Simulations with only Neon (Ne) impurity seeding have been conducted to evaluate its impact on meeting FNSF design demands for the divertor and upstream plasma parameters. Simulation results indicate that Ne seeding significantly mitigates divertor heat flux but potentially reduces both upstream electron and main ion density due to fuel dilution. The combined application of Ne seeding and deuterium (D 2 ) puffing is required to satisfy the FNSF design requirements on upstream density ($n^{OMP}_{e,sep}$~1× 10 20 m -3 ) and divertor energy flux ($q^{Odiv}_{\bot,max}$<10 MW m -2 ). D 2 puffing plays a role in counteracting upstream density drops and augmenting energy and momentum losses through atomic and molecular processes. The inlet Li flow velocity is systematically varied across a wide range to identify acceptable flows and corresponding LM surface temperatures. This comprehensive analysis identifies the acceptable Li flow parameters, LM surface temperature, and emitted Li fluxes necessary to meet the major design constraints. The emitted Li fluxes exhibit minimal impact on the main plasma at surface temperatures up to approximately ~525 °C, corresponding emitted Li fluxes of up to φ Li ~2X10 23 atoms s -1 . Uncertainties in the Li emission processes from the surface are also investigated, primarily influencing Li loss in the lower surface temperature range (<525 °C), with simulation results indicating a minor impact on the divertor and upstream plasma. Conversely, evaporation predominantly drives the Li loss processes at higher surface temperature ranges (>525 °C), contaminating both the divertor and upstream plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Exploring the effect of ELM and code-coupling frequencies on plasma and material modeling of dynamic recycling in divertors

Abstract Integrated modeling of plasma-surface interactions provides a comprehensive and self-consistent description of the system, moving the field closer to developing predictive and design capabilities for plasma facing components. One such workflow, including descriptions for the scrape-off-layer plasma, ion-surface interactions and the sub-surface evolution, was previously used to address steady-state scenarios and has recently been extended to incorporate time-dependence and two-way information flow. The new model can address dynamic recycling in transient scenarios, such as the application presented in this paper: the evolution of W samples pre-damaged by helium and exposed to ELMy H-mode plasmas in the DIII-D DiMES. A first set of simulations explored the effect of ELM frequency. This study was discussed in detail in this conference’s proceedings and is summarized here. The 2nd set of simulations, which is the focus of this paper, explores the effect of code-coupling frequency. These simulations include initial SOLPS solutions converged to the inter-ELM state, ion impact energy ( E in ) and angles ( A in ) calculated by hPIC2, and an improved heat transfer description in Xolotl. The model predicts increases in particle fluxes and decreases in heat fluxes by 10%–20% with the coupling time-step. Compared with the first set of simulations, the less shallow impact angle leads to smaller reflection rates and significant D implantation. The higher fraction of implanted flux (and deeper), in particular during ELMs, increases the accumulated D content in the W near-surface region. Future expansion of the workflow includes coupling to hPIC2 and GITR to ensure accurate descriptions of E in and A in , and W impurity transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Divertor-safe nonlinear burn control based on a SOLPS parameterized core-edge model for ITER

Abstract For ITER operations, the range of desirable burning-plasma regimes with high fusion power output will be restricted by various operational constraints. These constraints include the saturation of ITER’s various heating and fueling actuators such as the neutral beam injectors, the ion and electron cyclotron heating systems, the gas puffing system, and the deuterium–tritium pellet injectors. In addition to these actuator constraints, the H-mode power threshold, divertor detachment, and the heat load on the divertor targets may apply limitations to ITER’s operational space. In this work, Plasma Operation Contour (POPCON) plots that map the aforementioned constraints to the temperature-density space are used to investigate which constraints are most limiting towards accessing regimes with high fusion power output. The presented POPCON plots are based on a control-oriented core-edge model that couples the nonlinear density and energy response models for the core-plasma region with SOLPS4.3 parameterizations for conditions in the edge-plasma regions (scrape-off-layer and divertor). Using this control-oriented core-edge model, a nonlinear burn controller, which aims to regulate the plasma temperature and density in the core-plasma region, is constructed in this work. This controller is augmented with an online optimization scheme that governs the control references such that the plasma can be guided towards regimes with high fusion powers while protecting the divertor targets from dangerously high heat loads. A closed-loop simulation study illustrates the capability of this burn control scheme.

Physics↗

Interpretive modeling of tungsten divertor leakage during experiments with neon gas seeding

Abstract Many existing and future tokamaks with tungsten divertors operate, or will operate, with low- Z impurity seeding, but the direct effect of these seeded impurities on tungsten Scrape-off-Layer (SOL) transport has not been explored in detail. This paper reports on a DIII-D experiment designed to test how tungsten divertor leakage from the Small-Angle Slot V-Shaped, tungsten-coated divertor is impacted by neon seeding at a variety of injection rates and poloidal injection locations. Measurements from the experiment show an inverse relationship between the neon injection rate and the tungsten core penetration factor. Interpretive modeling is performed with a combination of the SOLPS-ITER and DIVIMP codes to assess the underlying tungsten behavior. The modeling results show that the reduction in tungsten divertor leakage is driven by both an increase in the divertor collisionality as well as a reduction in the ion temperature gradient near the divertor target. Collisions between low- Z impurities and tungsten impurities are found to have a significant impact on the tungsten SOL transport, such that ignoring the low- Z impurity collisional effects on the tungsten transport can result in an overestimate of the divertor leakage by an order-of-magnitude. Given the importance of these localized interactions, neon seeding from the closed, slot-like divertor has a clear advantage in being able to reduce tungsten divertor leakage without the high levels of neon core contamination that occur when seeding from other poloidal locations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Evaluation of SPARC divertor conditions in H-mode operation using SOLPS-ITER

The predicted divertor conditions for the SPARC tokamak are calculated using SOLPS-ITER for a range of scrape-off-layer (SOL) heat flux widths λ q , input powers, and particle fueling locations. Under H-mode scenario conditions with an upstream separatrix density of 1 x 10 20 m -3 , the most conservative range of λ q extrapolations ( 0.15 mm) results in extremely high unmitigated particle and energy fluxes to the divertor, both under full field (12.2 T) and power (P SOL = 29 MW) conditions, and 2/3 field with P SOL = 10 MW. Increasing the cross-field SOL diffusivities by 2–10× reduces the magnitude of the mitigation challenge, however strategies such as impurity seeding or strike-point-sweeping will likely still be required. A combination of steady-state and time-dependent SOLPS-ITER simulations are used to map out phase space diagrams of upstream and divertor conditions. The simulations include parallel currents but neglect cross-field drifts. At low upstream density the inner and outer divertor conditions are highly asymmetric, with a large temperature difference and significant heat fluxes driven by parallel currents. The solution has sharp bifurcations with a region of hysteresis, depending on whether the initial state is at a low or high density. This behavior is observed even when the fueling location, cross-field diffusivity, and impurity level is changed, although the density window with asymmetry is reduced with increasing diffusivity. The addition of neon impurity seeding reduces the divertor heat fluxes, but also causes a drop in the upstream electron density with fixed particle throughput. This drop can be counteracted by increased main ion throughput, however too much neon results in a back transition into the asymmetric divertor regimes suggesting a need for control of both main ion and impurity seeding levels to achieve a desired divertor state.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Dependence of divertor turbulence on plasma density and current in TCV

To reliably predict the distribution of heat and particle fluxes at the target plates of tokamaks, a comprehensive understanding of turbulence throughout the entire Scrape-Off-Layer (SOL) is imperative. This study examines divertor turbulence systematically across a broad parameter range on the TCV tokamak, including variations in magnetic field direction, plasma current I p ∈ [140,320] kA, edge safety factor q 95 ∈ [2.6,4.7] and Greenwald fraction f G ∈ [0.18,0.6]. The TCV X-point Gas Puff Imaging (GPI) system is used to measure 2D filament properties in the inner and outer divertor region. The fluctuation levels in the divertor are found to strongly increase with density (to 80% over most of the SOL) while remaining insensitive to I p . The previously identified divertor-localized filaments (DLF), located on the bad curvature side of the outer divertor leg, are found to be a common feature on TCV, while no filaments are observed in the PFR. DLFs are present over most of the parameter space and in both field directions. However, they are absent, or appear only closer to the target, for sufficiently large Λ div ≳ 10 or q 95 ≳ 3.7. Across both I p and f G scans, some clear trends with Λ div are found for divertor filament sizes and velocities, and with target fall-off lengths of density and heat flux profiles at the outer target. This study provides important experimental insights to turbulent transport in the divertor also for comparison with self-consistent, turbulence simulations and extrapolation to future reactor conditions.

cross-field transport↗

Integrated modeling of RF-induced tungsten erosion at ICRH antenna structures in the WEST tokamak *

This paper introduces STRIPE (Simulated Transport of RF Impurity Production and Emission), an advanced modeling framework developed to analyze material erosion and the global transport of eroded impurities originating from radio-frequency (RF) antenna structures in magnetic confinement fusion devices. STRIPE integrates multiple physics modules: SolEdge3x for scrape-off-layer plasma profiles, COMSOL for 3D RF rectified sheath potentials, RustBCA for erosion yields and surface interactions, and global impurity transport for 3D ion energy-angle distributions and impurity transport. The framework is applied to an ion cyclotron RF-heated L-mode discharge (#57877) in the WEST tokamak, where it predicts a thirty-fold increase in gross tungsten erosion at antenna limiters during the transition from ohmic to ICRH operation. Additionally, under ICRH conditions, a tenfold enhancement in erosion is observed when comparing RF sheath effects to purely thermal sheath conditions. High-charge-state oxygen ions ($\mathrm{O}$ 6+ and above) are identified as the dominant contributors to tungsten sputtering. To validate the model, a synthetic diagnostic tool based on inverse photon efficiency (S/XB coefficients) from the ColRadPy collisional-radiative model enables direct comparison with spectroscopic measurements. Model predictions using a plasma composition of 1% oxygen and 99% deuterium show good agreement with observed W − I (400.9 nm) emission for discharge #57877, supporting the accuracy of the STRIPE framework. This study focuses specifically on gross erosion calculations to demonstrate STRIPE’s capabilities. Future extensions of this work will incorporate net erosion, re-deposition, self-sputtering effects, and whole-device modeling of sputtered tungsten impurity transport. STRIPE is also being applied to other RF-heated linear and toroidal devices, offering valuable insights for antenna design, impurity control, and performance optimization in next-generation fusion reactors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Overview of ST40 results and future: expanding the physics basis of high-field spherical tokamaks

The goal of the ST40 programme is to explore the physics of high-field spherical tokamaks (STs), to validate empirical and theoretical models and, hence, to build confidence in predictions required to support the design of future generations of STs. ST40 is a compact high-field ST that has achieved the following parameters: R 0 = 0.4–0.55 m, I p = 0.20–0.85 MA, B t (R= 0.4 m) = 0.7–2.1 T, κ ⩽ 1.9, and A = 1.6–1.9. Highlights of recent experimental results include (i) H-mode and confinement studies at B t ⩽ 2.1 T, (ii) observation of bifurcation of the scrape-off-layer power fall-off width, λ q , into a ‘wide’ branch that follows existing H-mode scalings and a ‘narrow’ branch that exhibits λ q values that are up to 10 times lower than the predictions of established scalings, (iii) development of high-performance scenarios with plasma current, I p , up to 0.85 MA, (iv) development of highly non-inductive scenarios with high β p , and (v) the first ST40 experiments utilising the newly commissioned impurity powder dropper. The work on all these topics has been supported by a number of advancements in ST40 hardware and software, from plasma control to data analysis and interpretation. At the end of 2025, ST40 embarked on a major upgrade to further expand its capabilities by introducing, among other improvements, all-metal plasma-facing components, 1 MW of electron cyclotron heating, a pellet injector, and a pair of lithium evaporators for wall conditioning.

confinement↗

Overview of Wendelstein 7-X high-performance operation

The Wendelstein 7-X (W7-X) stellarator has completed two consecutive experimental campaigns OP 2.2 (Sep.-Dec. 2024) and OP 2.3 (Feb.-May 2025) under a new operational strategy enabling more than one year of uninterrupted device availability. This approach, supported by exceptionally high subsystem reliability, allowed sustained high-efficiency plasma operations with up to 80–100 discharges per day across a broad range of magnetic configurations. Several key technical upgrades-most notably the first operation of a 1.5 MW class steady-state gyrotron, a new steady-state pellet injector, and advanced real-time feedback control systems significantly enhanced heating, fueling, and plasma control capabilities. Together, these improvements enabled major advances in long-pulse performance, high-β operation, and confinement optimization. Long-pulse discharges achieved 1.8 GJ of injected energy under fully detached divertor conditions, while reduced-field scenarios facilitated record volume-averaged β values approaching 3%. High-performance plasmas with centrally peaked density profiles, created via neutral beam injection (NBI) or sustained pellet fueling, demonstrated strongly reduced turbulent transport and stellarator-record fusion triple products. Complementary studies of power exhaust and divertor heat loads revealed the role of scrape-off-layer drift physics in shaping strike-line patterns under attached conditions. Together, the results from OP 2.2 and OP 2.3 significantly expand the operational space of W7-X and strengthen its role as a leading platform for steady-state stellarator research and reactor-relevant plasma scenarios.

Wendelstein 7-X↗

Physics basis for the reference flat-top plasma scenario in the ST–E1 fusion power plant

As part of the U.S. Department of Energy’s Milestone-Based Fusion Energy Development Program, Tokamak Energy has completed the pre-concept design of the ST–E1 fusion power plant. ST–E1 is envisaged to operate in two phases: a pilot plant phase, targeting sustained net power production of 300 - 500 MWe for a duration >1 hr, followed by a commercial power plant phase targeting steady-state operations and a normalised overnight capital cost of ⩽12 000 $\$$/kWe. The design process adopted was highly iterative, integrating all major plant systems and progressing in a phased fidelity approach. At the pre-conceptual stage, the emphasis has been on exploring the design space, identifying the main system-level trade-offs, and making the key decisions that define the overall plant concept, rather than optimising a single operating point. This paper, part of a focused collection detailing the ST–E1 pre-concept design, addresses the development of a series of reference flat-top plasma operating points for the pilot plant phase. A modelling workflow was established to develop and assess candidate plasma design points and explore key dependencies. The workflow includes integrated core plasma modelling, magnetohydrodynamic (MHD) stability assessment, equilibrium generation, scrape-off-layer and exhaust modelling, heating & current drive design and optimisation, and turbulent transport modelling. Using this framework, the impact of several key parameters on the flat-top operating space was investigated, including the density limit, core radiation fraction and divertor power loading, level of external heating and curent drive power and assumed pedestal characteristics. The MHD stability, controllability and micro-stability characteristics of these plasmas were also analysed. These investigations informed the definition of a set of fully non-inductive, flat-top reference operating points that satisfy the high-level ST–E1 mission, including a low and high density case, a case that is stable to resistive wall modes and a case with reduced divertor power loading.

ST–E1↗

Boundary, SOL, and Divertor Physics Studies on TCV (Final Report for Research Funded by DoE Award DE-SC0020327)

Research supported under this grant was organized around two diagnostic systems installed on the TCV tokamak located at EPFL in Lausanne, Switzerland: Multi-Spectral Imaging (MSI) and Gas-Puff Imaging (GPI). MSI monitors 2-d time-history images of up to 10 spectral lines, in the visible, using identical views for all wavelengths. The primary application of MSI is to infer 2-d time histories of electron temperature and density in the Scrape-Off-Layer (SOL) of diverted plasmas, utilizing ratios of emission from multiple neutral helium lines. The atomic models are validated against 1-d Thomson scattering measurements. SOL density and temperature maps are analyzed to infer divertor dynamics, including particle and power flows, for a variety of divertor configurations, and specifically as flux expansion along the long-legged outboard divertor is systematically varied. GPI is used to measure details of edge fluctuations in the SOL, particularly as a function of plasma triangularity (δ), including both positive and negative δ. A threshold is found, whereby for δ<-0.3, edge turbulence in the far SOL is completely suppressed. Dynamics of the filamentary turbulence was also analyzed using advance machine learning techniques. Research in support of four PhD theses was carried out within the scope of this Award, including 2 at MIT and 2 at EPFL.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Global 3D Simulations of the Tokamak Edge Region

The goal of this project is to better understand magnetized plasma dynamics at the edge region of the tokamak, including both closed-flux-surface region and the scrape-off-layer (SOL), through analytical analysis and advanced numerical simulations. Under the auspices of the U.S. DOE, the PI, along with students and collaborators have made a significant progress. These results were disseminated by 13 peer-reviewed article and many presentations and posters in domestic and international conferences.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Center for Integrated Simulation of Fusion Relevant RF Actuators

This project was part of the “Center for Integrated Simulation of Fusion Relevant RF Actuators” SciDAC-4 project, led by Dr. Paul Bonoli (MIT). Rather than use an acronym (CISFRRFA), the project will be referred to in this document as the “RF-SciDAC4”. The larger SciDAC-4 project goals were to: 1. Develop an integrated simulation of the antenna + sheath + scrape-off-layer + core plasma system which fully utilizes leadership class computing resources to move towards a quantitative predictive capability for the response to RF power. 2. Work closely with the SciDAC-4 Whole Device Modeling (WDM) community to make both our new code development efforts, as well as the established hierarchy of RF tools, available within their environment, and to utilize WDM technologies to implement the couplings below. 3. Validate this predictive capability on appropriately diagnosed experiments including dedicated RF test stands, linear devices, and existing tokamaks. 4. Use these tools to inform design of robust, impurity-mitigating RF heating and current drive sources for future fusion devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Fusion-Enabled Pluto Orbiter and Lander

The Pluto orbiter mission proposed here is credible and exciting. The benefits to this and all outer-planet and interstellar-probe missions are difficult to overstate. The enabling technology, Direct Fusion Drive, is a unique fusion engine concept based on the Princeton Field-Reversed Configuration (PFRC) fusion reactor under development at the Princeton Plasma Physics Laboratory. The truly game-changing levels of thrust and power in a modestly sized package could integrate with our current launch infrastructure while radically expanding the science capability of these missions. During this Phase I effort, we made great strides in modeling the engine efficiency, thrust, and specific impulse and analyzing feasible trajectories. Based on 2D fluid modeling of the fusion reactors outer stratum, its scrape-off-layer (SOL), we estimate achieving 2.5 to 5 N of thrust for each megawatt of fusion power, reaching a specific impulse, Isp, of about 10,000 s. Supporting this model are particle-in-cell calculations of energy transfer from the fusion products to the SOL electrons. Subsequently, this energy is transferred to the ions as they expand through the magnetic nozzle and beyond. Our point solution for the Pluto mission now delivers 1000 kg of payload to Pluto orbit in 3.75 years using 7.5 N constant thrust. This could potentially be achieved with a single 1 MW engine. The departure spiral from Earth orbit and insertion spiral to Pluto orbit require only a small portion of the total delta-V. Departing from low Earth orbit reduces mission cost while increasing available mission mass. The payload includes a lander, which utilizes a standard green propellant engine for the landing sequence. The lander has about 4 square meters of solar panels mounted on a gimbal that allows it to track the orbiter, which beams 30 to 50 kW of power using a 1080 nm laser. Optical communication provides dramatically high data rates back to Earth. Our mass modeling investigations revealed that if current high-temperature superconductors are utilized at liquid nitrogen temperatures, they drive the mass of the engine, partly because of the shielding required to maintain their critical temperature. Second generation materials are thinner but the superconductor is a very thin layer deposited on a substrate with additional layers of metallic classing. Tremendous research is being performed on a variety of these superconducting materials, and new irradiation data is now available. This raises the possibility of operating nearfuture high-temperature superconductors at a moderately low temperature to dramatically reduce the amount of shielding required. At the same time, a first generation space engine may require low-temperature superconductors, which are higher TRL and have been designed for space coils before (AMS-02 experiment for the ISS). We performed detailed analysis of the startup system and thermal conversion system components. The ideal working fluid was determined to be a blend of Helium and Xenon. No significant problems were identified with these subsystems. For the RF system, we conceived of a new, more efficient design using state-of-the-art switch amplifiers, which have the potential for 100% efficiency. This report presents details of our engine and trajectory analyses, mass modeling efforts, and updated vehicle designs.

temperature↗

Compression of tokamak boundary plasma simulation data using a maximum volume algorithm for matrix skeleton decomposition

This report demonstrates satisfactory data compression of SOLPS-ITER simulation output ranging from 2D fields, 1D profiles, and 0D scalar variables with a novel matrix decomposition approach. The singular value decomposition (SVD) scales poorly for large matrix sizes and is unsuited to the application on high dimensional data common to fusion plasma physics simulation. In this work, we employ the columns-submatrix-rows (CUR) matrix factorization technique in order to compute a low-rank approximation up to two orders of magnitude faster than the SVD, but within a nominal L2-norm relative error of ε = 10 –2 . In addition, the CUR approach maintains the original format of the data, in its extracted columns and rows, allowing for interpretable data storage at the original resolution of the simulation. We utilize an iterative algorithm to compute the CUR decomposition of simulation output by maximizing the volume, or linearly independent information content, of a low-rank submatrix contained within the data. Experiments over $\textit{n} × \textit{n}$ randomized test matrices with embedded rank-deficient features show that this maximum volume implementation of CUR matrix approximation has reduced asymptotic computational complexity on the order of n compared to the SVD, which scales approximately as $n^3$. These results show that the CUR technique can be used to effectively select time step snapshots (columns) of over 140 SOLPS-ITER output variables and the associated discretized coordinate timeseries (rows) allowing for reconstruction of the complete simulation dynamics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

1D model of SOL with self consistent calculation of non-coronalradiation and applications to shallow field line angles

A 1D model of the scrape off layer (SOL) is created with a physics based calculation of the impurity radiation and the inclusion of the perpendicular heat flux q ⊥ compared to previous 1D models. The calculation of the impurity radiation includes the addition of a model using parallel and perpendicular motion to calculate the impurity confinement time τ and the inclusion of charge exchange between neutral hydrogen and impurity ions. Reasonable agreement is found between the model and the 2D fluid code SOLPS-ITER. Additionally, an improvement is shown over the constant τ used in previous works The model is used to examine the effects of shallow field line angles at ITER levels of parallel heat flux (q || ). A large increase in radiation results in an increase in the upstream parallel heat flux which can still be in detachment at shallow field line angles for both q ⊥ and the calculated τ for traditional impurities (argon and nitrogen) as well as boron. The physics behind this increase in radiation at shallow field line angles is shown to be a decrease in τ and a change in the temperature profile with the inclusion of q ⊥ .

ITER↗