SEARCH · Engineering Papers

Results for “accelerator modelling and simulations”

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.

SourceResearch connectionFrom yearThrough year

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

At least 379 records · Page 21

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Computational Design of Alloys for Energy Technologies

Advanced materials that maintain their mechanical performance under elevated temperatures, corrosive environments, and a range of static and evolving stresses are needed to improve the efficiency and reduce the environmental impact of a wide spectrum of energy technologies. For instance, cost-efficient alloys that can withstand high temperatures (e.g., 700 °C) have a critical role in improving the efficiency and economics of power generation to support decarbonization of the energy sector; such is true of both the nuclear and fossil energy sectors. Considering both the threats of the energy crisis, namely soaring costs of greenhouse gas emission-producing energy and climate change, it is essential to increase the pace of material discovery and enable rapid paths for material qualification to advance clean energy technologies. Conventionally, alloy development has followed a slow Edisonian process that uses repeated cycles of making, characterizing, and modifying to arrive at optimum composition and processing conditions to achieve the desired component performance. This optimization is followed by the necessary stepwise materials qualification. Furthermore, the increasing adoption of sound data management and physics-informed machine learning represents the next step in the acceleration of materials design and development. In the integrated computational materials engineering (ICME) approach, computational modeling and simulation data from different length and time scales can be combined with complex microstructural details from multimodal experimental characterization and selective property testing to close the design loop for rapid alloy development.

Computational Design Of Materials↗

DOE OSTI2024

A real-time energy and cost efficient vehicle route assignment neural recommender system

Here, this paper presents a neural network recommender system algorithm for assigning vehicles to routes based on energy and cost criteria. In this work, we applied this new approach to efficiently identify the most cost-effective medium and heavy duty truck (MDHDT) powertrain technology, from a total cost of ownership (TCO) perspective, for given trips. We employ a machine learning based approach to efficiently estimate the energy consumption of various candidate vehicles over given routes, defined as sequences of links (road segments), with little information known about internal dynamics, i.e. using high level macroscopic route information. A complete recommendation logic is then developed to allow for real-time optimum assignment for each route, subject to the operational constraints of the fleet. We show how this framework can be used to (1) efficiently provide a single trip recommendation with a top-k vehicles star ranking system, and (2) engage in more general assignment problems where n vehicles need to be deployed over m (m ≤ n) trips. This new assignment system has been deployed and integrated into the POLARIS. Transportation System Simulation Tool for use in research conducted by the Department of Energy's Systems and Modeling for Accelerated Research in Transportation (SMART) Mobility Consortium (SMART, 2024).

Energy consumption↗

DOE OSTI2021

Machine learning accelerated turbulence modeling of transient flashing jets

We report modeling the sudden depressurization of superheated liquids through nozzles is a challenge because the pressure drop causes rapid flash boiling of the liquid. The resulting jet usually demonstrates a wide range of structures, including ligaments and droplets, due to both mechanical and thermodynamic effects. As the simulation comprises increasingly numerous phenomena, the computational cost begins to increase. One way to moderate the additional cost is to use machine learning surrogacy for specific elements of the calculation. This study presents a machine learning-assisted computational fluid dynamics approach for simulating the atomization of flashing liquids accounting for distinct stages, from primary atomization to secondary breakup to small droplets using the Σ - Y model coupled with the homogeneous relaxation model. Notably, the models for thermodynamic non-equilibrium (HRM) and Σ - Y are coupled, for the first time, with a deep neural network that simulates the turbulence quantities, which are then used in the prediction of superheated liquid jet atomization. The data-driven component of this method is used for turbulence modeling, avoiding the solution of the two-equation turbulence model typically used for Reynolds-averaged Navier-Stokes simulations for these problems. Both the accuracy and speed of the hybrid approach are evaluated, demonstrating adequate accuracy and at least 25% faster computational fluid dynamics simulations than the traditional approach. This acceleration suggests that perhaps additional components of the calculation could be replaced for even further benefit. Published under an exclusive license by AIP Publishing.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2022

Modeling ablator grain structure impacts in ICF implosions

High-density carbon is a leading ablator material for inertial confinement fusion (ICF). This and some other ablator materials have grain structure which is believed to introduce very small-scale (~nm) density inhomogeneity. In principle, such inhomogeneity can affect key ICF metrics like fuel compression and yield, by, for example, acting as a seed for instabilities and inducing mix between ablator and fuel. However, assessments of such effects are uncertain due to the difficulty of modeling this small-scale structure in ICF simulations, typically requiring reduced-resolution modeling that scales these features. We present a grain model and show both the impact of de-resolving grains and the complex mixing dynamics such structures can induce. We find that different methods for de-resolving grains can yield both different total deposition of kinetic energy perturbations and different fuel–ablator mixing. We then show a simple-to-implement approach for approximately conserving the deposition of perturbed kinetic energy and demonstrate that, for the present grain model and test cases, this approach yields a reasonably matched time history of mix width between less and more resolved grain models. The simulations here also demonstrate the complex interaction history between grain-induced mixing and instability around the fuel–ablator interface, showing, for example, that the grain-induced perturbations typically trigger instability of conduction-driven density gradients in the DT fuel, enhancing mix penetration early in the acceleration of the shell. Furthermore simulating both microscale and nanoscale grains, we find initial evidence for larger mixing in the microscale case of the present model, despite smaller deposited kinetic energy perturbation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2026

Strong-field QED limitations on TeV-class plasma wakefield accelerators

We demonstrate that quantum and classical radiation effects can become non-negligible for TeV-class beams propagating through plasma channels typical of staged plasma accelerators. Although the quantum nonlinearity parameter χ e remains small under currently envisioned experimental conditions, the cumulative influence of radiation over long acceleration distances can lead to significant modifications to the beam’s energy spread, emittance, and polarization. Our analytic models, validated by particle-in-cell simulations, highlight that for standard Gaussian beams, the orbit-induced energy spread dominates over quantum stochastic effects but can be mitigated by tailoring the beam profile, for example, through ring-shaped transverse distributions. In regimes where the radiation reaction approaches the accelerating force, the emittance may be cooled, forming distinctive ring-shaped phase-space structures. Finally, we analyze the influence of the radiation effect on spin transport inside the wakefield. These findings underscore the importance of considering both classical and quantum radiation dynamics in the design and optimization of future high-intensity plasma accelerators.

accelerator↗

DOE OSTI2025

GPR_calculator: An on-the-fly surrogate model to accelerate massive nudged elastic band calculations

We present GPR_calculator, a package based on Python and C++ programming languages to build an on-the-fly surrogate model using Gaussian Process Regression (GPR) to approximate computationally expensive electronic structure calculations. The key idea is to dynamically train a GPR model during the simulation that can accurately predict energies and forces with uncertainty quantification. When the uncertainty is high, the costly electronic structure calculation is performed to obtain the ground truth data, which is then used to update the GPR model. To illustrate the effectiveness of GPR_calculator, we demonstrate its application in Nudged Elastic Band (NEB) simulations of surface diffusion and reactions, achieving 3-10 times acceleration compared to pure ab initio calculations. The source code is available at https://github.com/MaterSim/GPR_calculator.

Gaussian process regression↗

DOE OSTI2025

Energy Spectrum and Mass Composition of Ultra-high-energy Cosmic Rays Originating from Relativistic Jets of Nearby Radio Galaxies

Relativistic jets of radio galaxies (RGs) are possible sources of ultra-high-energy cosmic rays (UHECRs). Recent studies combining relativistic hydrodynamic simulations with Monte Carlo particle transport have demonstrated that UHECRs can be accelerated to energies beyond 10 20 eV through shocks, turbulence, and relativistic shear in jet-induced flows of Fanaroff–Riley type RGs. The resulting time-asymptotic UHECR spectrum is well modeled by a double power law with an “extended” exponential cutoff, primarily shaped by relativistic shear acceleration. In this study, we adopt this novel source spectrum and simulate the propagation of UHECRs from nearby RGs using the CRPropa code. We focus on Virgo A (Vir A), Centaurus A (Cen A), Fornax A (For A), and Cygnus A (Cyg A), expected to be the most prominent UHECR sources among RGs. We then analyze the energy spectrum and mass composition of UHECRs arriving at Earth. We find that, due to the extended high-energy tail in the source spectrum, UHECRs from Vir A, which has a higher Lorentz factor, exhibit a higher flux at the highest energies and a lighter mass composition at Earth compared to those from Cen A and For A with lower Lorentz factors. Despite Cyg A having an even higher Lorentz factor, the large distance limits its contribution. With a small number of nearby prominent RGs, our findings suggest that if RGs are the major sources of UHECRs, the energy spectrum and mass composition of observed UHECRs would exhibit hemispheric differences between the Northern and Southern skies at the highest energies.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Validation of material models for puncture of 7075-T651 aluminum plate

Plate puncture simulations are challenging computational tasks that require advanced material models including high strain rate and thermal-mechanical effects on both deformation and failure, plus finite element techniques capable of representing large deformations and material failure. The focus of this work is on the material issues, which require large sets of experiments, flexible material models and challenging calibration procedures. In this study, we consider the puncture of 12.7 mm thick, 7075-T651 aluminum alloy plates by a cylindrical punch with a hemispherical nose and diameter of 12.7 mm. The plasticity and ductile failure models were isotropic with calibration data obtained from uniaxial tension tests at different temperatures and strain rates plus quasi-static notched tension tests and shear-dominated tests described here. Sixteen puncture experiments were conducted to identify the threshold penetration energy, mode of puncture and punch acceleration during impact, The punch was mounted on a 139 kg mass and dropped on the plates with different impact speeds. Since the mass was the same in all tests, the quantity of interest was the impact speed. The axis and velocity of the punch were perpendicular to the plate surface. The mean threshold punch speed was 3.05 m/s, and the mode of failure was plugging by thermal-mechanical shear banding accompanied by scabbing fragments. Application of the material models in simulations of the tests yielded accurate estimates of the threshold puncture speed and of the mode of failure. Time histories of the punch acceleration compared well between simulation and test. Remarkably, the success of the simulations occurred in spite of even the smallest element used being larger than the width of the shear bands.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Interpreting experimental measurements of helium bubbles using stochastic cluster dynamics models of heterogeneous nucleation and growth in irradiated ferritic alloys

Among a number of other advantageous features, ferritic/martensitic steels are being considered as fusion reactor structural materials due to their low intrinsic swelling under irradiation. However, under high-energy neutron irradiation, He produced through (n, α) reactions stabilizes vacancy clusters, which then act as seeds for bubble formation and growth, negating the intrinsic swelling resistance of these alloys. Standard models of irradiation damage accumulation and microstructural evolution consider homogeneous nucleation as the basis for bubble formation and growth. However, this generally does not explain the large bubble densities and sizes observed experimentally under a number of different conditions. Here, we propose a more realistic physical model of bubble nucleation, formation, and growth designed to capture recent experimental measurements of He-bubble formation and evolution during co-implantation of He and Fe ions in model ferritic alloys. We find that experimental results are explained only when the following three features are considered simultaneously: (i) heterogeneous nucleation of He-vacancy bubbles at defect sinks (e.g., dislocations, grain boundaries, and second-phase precipitates), (ii) nucleation and growth of bubbles via the ‘trap mutation’ mechanism (i.e., spontaneous production of Frenkel pairs due to absorption of He atoms), and (iii) transition from a growth-limited, He-stabilized bubble regime to a ‘runaway’ void-type growth scenario. The model is implemented into a stochastic cluster dynamics framework that takes advantage of cluster size grouping methods to accelerate the simulations, allowing us to reach 10 dpa of simulated irradiated dose, and to capture cluster sizes in excess of 20 nm. Further, a careful extrapolation exercise conducted assuming classical nucleation theory leads to excellent agreement with the experimental measurements at 50 dpa of irradiation.

36 MATERIALS SCIENCE↗

DOE OSTI2022

A multiphysics model of the versatile test reactor based on the MOOSE framework

The traditional modeling approach for sodium fast reactor cores relies on separate physics models, where the fuel performance, thermal–hydraulics, and neutronics calculations required to predict the core physics characteristics for nominal conditions are decoupled by relying on user-imposed boundary conditions. Here, this paper aims at evaluating the impact of multiphysics simulations for predicting the core characteristics of the Versatile Test Reactor, which is being designed as a 300-MWt sodium-cooled fast reactor. The purpose of the Versatile Test Reactor is to accelerate the testing of advanced nuclear materials in the United States. The proposed multiphysics model relies on the Griffin reactor physics code, the SAM thermal–hydraulic system code, the BISON fuel performance code, as well as generic Multiphysics Object-Oriented Simulation Environment capabilities implemented in the open-source tensor mechanics module. For k eff calculations, the introduction of a tight coupling between the neutronics, thermo-mechanical and thermal–hydraulics models induces a change of around 543 pcm in the eigenvalue, compared to the traditional standalone neutronics calculation where approximate temperature profiles are used. The multiphysics model is then employed for quantifying the impact of the thermal conductivity uncertainties on some of the key figures of merit, such as the fuel centerline temperature, assembly powers, and keff for nominal core conditions. As anticipated, uncertainties on fuel thermal conductivity mostly impact the fuel centerline temperature, and to a lesser extend the k eff .

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2026

From Cell to System: Accelerated hpc Simulations of BESS Aging under Frequency Regulation and Arbitrage use cases

Lithium-ion battery energy storage systems (BESS) packs have emerged as a leading solution for grid-scale energy storage, enhancing resiliency and balancing load fluctuations. Yet, experimental characterization of large-format LIB packs-particularly to assess performance and degradation over hundreds of cycles - demands substantial hardware investment and multi-year testing campaigns. In this work, we couple a hierarchical, physics-based modeling framework agnostic to electrode chemistries with high-performance computing to accelerate systems level evaluation by upto two orders of magnitude. Building on the open-source liionpack platform, we implement cell, module, and pack-scale electrochemical models enriched with mechanistic aging mechanisms and deploy them on an HPC cluster to simulate 150−200kWh systems over 500 - 1,000 cycles with in days. We subject these virtual B ESS to both constant-current cycling and realistic grid service profiles spanning frequency regulation, ramp-rate support, and energy arbitrage-and quantify the resulting degradation patterns. Our results reveal that localized cell aging can induce substantial nonuniformity at module and pack levels, with service-specific cycling protocols driving distinct aging modes. This rapid, multiscale modeling approach provides a powerful design-space exploration tool for optimizing electrical architecture, control strategies, and operational schedules to prolong pack lifetime and lower total cost of ownership.

Ayalasomayajula, Surya [ORNL] (ORCID:0009000860788↗

DOE OSTI2024

Accelerated Irradiation and Qualification of Ceramic Nuclear Fuels

Accelerated neutron irradiation testing is an component of accelerated qualification of new nuclear fuels for light water reactor (LWRs), microreactors, and other special purpose reactors. The qualification and licensing of nuclear fuel is a lengthy process that can take 20-25 years to bring a new fuel into service. Accelerated fuel qualification combines both experimental and modeling work to expedite the total qualification time to 5-10 years timeframe. The experimental aspect of this is accelerated irradiation aims to reduce the total time needed for neutron irradiation to achieve targeted burnup, which can take years using conventional irradiation profiles. The data that results from this irradiation testing can then be entered into BISON models to develop robust and reliable performance simulations to ensure safe operation under normal and off normal conditions. This milestone focused on the fabrication of test articles for accelerated irradiation testing at the Advance Test Reactor (ATR).

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2025

Performance of an optical TPC Geant4 simulation with opticks GPU-accelerated photon propagation

We investigate the performance of Opticks, a NVIDIA OptiX API 7.5 GPU-accelerated photon propagation tool compared with a single-threaded Geant4 simulation. We compare the simulations using an improved model of the NEXT-CRAB-0 gaseous time projection chamber. Performance results suggest that Opticks improves simulation speeds by between 58.47 ± 0.02 and 181.39 ± 0.28 times relative to a CPU-only Geant4 simulation and these results vary between different types of GPU and CPU. A detailed comparison shows that the number of detected photons, along with their times and wavelengths, are in good agreement between Opticks and Geant4.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2021

Simulating low-energy neutrino interactions with MARLEY

Monte Carlo event generators are a critical tool for the interpretation of data obtained by neutrino experiments. Several modern event generators are available which are well-suited to the GeV energy scale used in studies of accelerator neutrinos. However, theoretical modeling differences make their immediate application to lower energies difficult. In this paper, I present a new event generator, MARLEY, which is designed to better address the simulation needs of the low-energy (tens of MeV and below) neutrino community. The code is written in C++14 with an optional interface to the popular ROOT data analysis framework. The current release of MARLEY (version 1.2.0) emphasizes simulations of the reaction 40 Ar ( ν e , e − ) 40 K ⁎ but is extensible to other channels with suitable user input. This paper provides detailed documentation of MARLEY's implementation and usage, including guidance on how generated events may be analyzed and how MARLEY may be interfaced with external codes such as Geant4. Further information about MARLEY is available on the official website at http://www.marleygen.org. Program title:MARLEY 1.2.0 CPC Library link to program files:https://doi.org/10.17632/4v7zxnc8j3.1 Developer's respository link:http://github.com/MARLEY-MC/marley Code Ocean capsule:https://codeocean.com/capsule/9868179 Licensing provisions: GNU General Public License 3.0 Programming language: C++14 External routines/libraries used: GNU Scientific Library [1,2] (required), ROOT [3,4] (optional) Nature of problem: Simulation of neutrino-nucleus scattering events at energies of tens-of-MeV and below Solution method: Initial two-to-two scattering kinematics are sampled using the allowed approximation differential cross section and tables of precomputed nuclear matrix elements. Subsequent de-excitations of the remnant nucleus are simulated using a Monte Carlo implementation of the Hauser-Feshbach statistical model and tabulated γ-ray decay schemes for discrete nuclear levels. Additional comments including restrictions and unusual features: Input data are provided with the code that are suitable for producing simulations of the charged-current reaction 40 Ar ( ν e , e − ) 40 K ⁎ , coherent elastic neutrino-nucleus scattering on spin-zero target nuclei, and neutrino-electron elastic scattering on any atomic target. Preparation of new reaction input files (whose format is documented in Appendix B) would enable other reaction channels and nuclear targets to be handled by the existing code framework. Although there is no maximum neutrino energy enforced by the code itself, realistic neutrino-nucleus scattering events may be generated up to roughly 50 MeV. Above this energy, the effects of forbidden nuclear transitions, which are neglected in the current treatment of the cross sections (see section 2.1), become increasingly important.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2021

Size-dependent radiation damage mechanisms in nanowires and nanoporous structures

Nanostructures with a high density of interfaces, such as in nanoporous materials and nanowires, resist radiation damage by promoting the annihilation and migration of defects. This study details the size effect and origins of the radiation damage mechanisms in nanowires and nanoporous structures in model face-centered (gold) and body-centered (niobium) cubic nanostructures using accelerated multi-cascade atomistic simulations and in-situ ion irradiation experiments. Our results reveal three different size-dependent mechanisms of damage accumulation in irradiated nanowires and nanoporous structures: sputtering for very small nanowires and ligaments, the formation and accumulation of point defects and dislocation loops in larger nanowires, and a face-centered-cubic to hexagonal-close-packed phase transformation for a narrow range of wire diameters in the case of gold nanowires. Smaller nanowires and ligaments have a net effect of lowering the radiation damage as compared to larger wires that can be traced back to the fact that smaller nanowires transition from a rapid accumulation of defects to a saturation and annihilation mechanism at a lower dose than larger nanowires. These irradiation damage mechanisms are accompanied with radiation-induced surface roughening resulting from defect-surface interactions. Comparisons between nanowires and nanoporous structures show that the various mechanisms seen in nanowires provide adequate bounds for the defect accumulation mechanisms in nanoporous structures with the difference attributed to the role of nodes connecting ligaments in nanoporous structures. Altogether, our results shed light on the compounded, size-dependent mechanisms leading to the radiation resistance of nanowires and nanoporous structures.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Prediction of performance and turbulence in ITER burning plasmas via nonlinear gyrokinetic profile prediction

Burning plasma performance, transport, and the effect of hydrogen isotope (H, D, D-T fuel mix) on confinement has been predicted for ITER baseline scenario (IBS) conditions using nonlinear gyrokinetic profile predictions. Accelerated by surrogate modeling (Rodriguez-Fernandez et al 2022 Nucl. Fusion 62 076036), high fidelity, nonlinear gyrokinetic simulations performed with the CGYRO code (Candy et al 2016 J. Comput. Phys. 324 73), were used to predict profiles of T i , T e , and n e while including the effects of alpha heating, auxiliary power (NBI + ECH), collisional energy exchange, and radiation losses inside of $r/a$ = 0.9. Predicted profiles and resulting energy confinement are found to produce fusion power and gain that are approximately consistent with mission goals ($P_\textrm{fusion} = 500$ MW at Q = 10) for the baseline scenario and exhibit energy confinement that is within 1σ of the H-mode energy confinement scaling. The power of the surrogate modeling technique is demonstrated through the prediction of alternative ITER scenarios with reduced computational cost. These scenarios include conditions with maximized fusion gain and an investigation of potential resonant magnetic perturbation (RMP) effects on performance with a minimal number of gyrokinetic profile iterations required (3–6). These predictions highlight the stiff ITG nature of the core turbulence predicted in the ITER baseline and demonstrate that $Q \gt$ 17 conditions may be accessible by reducing auxiliary input power while operating in IBS conditions. Prediction of full kinetic profiles allowed for the projection of hydrogen isotope effects around ITER baseline conditions. The gyrokinetic fuel ion species was varied from H, D, and 50/50 D-T and kinetic profiles were predicted. Results indicate that a weak or negligible isotope effect will be observed to arise from core turbulence in IBS conditions. The resulting energy confinement, turbulence, and density peaking, and the implications for ITER operations will be discussed.

gyrokinetics↗

DOE OSTI2026

Goated: goal-oriented tensor decompositions in python

SAND2026-20464O Goated performs goal-oriented tensor decompositions in Python, enabling efficient compression of multi-dimensional simulation data. It extends common tensor decomposition methods by incorporating domain-specific knowledge, such as conservation laws in physics, through a penalty term in the optimization process. This approach improves data compression and modeling accuracy across various applications, including physics simulations, by using specialized algorithms and structure-aware subroutines to accelerate solver performance. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy's National Nuclear Security Administration under contract DE-NA0003525.

SciDAC↗

DOE OSTI2023

Millimeter observational signatures of flares in magnetically arrested black hole accretion models

In general relativistic magnetohydrodynamic (GRMHD) simulations, accreted magnetic flux on the black hole horizon episodically decays, during which magnetic reconnection heats up the plasma near the horizon, potentially powering high-energy flares like those observed in M87* and Sgr A*. We study the mm observational counterparts of such flaring episodes in very high resolution GRMHD simulations. The change in 230 GHz flux during the expected high energy flares depends primarily on the efficiency of accelerating γ ≳ 100 (T e ≳ 10 11 K) electrons. For models in which the electrons are heated to T e ∼ 10 11 K during flares, the hot plasma produced by reconnection significantly enhances 230 GHz emission and increases the size of the 230 GHz image. By contrast, for models in which the electrons are heated to higher temperatures (which we argue are better motivated), the reconnection-heated plasma is too hot to produce significant 230 GHz synchrotron emission, and the 230 GHz flux decreases during high energy flares. We do not find a significant change in the mm polarization during flares as long as the emission is Faraday thin. We also present expectations for the ring-shaped image as observed by the Event Horizon Telescope during flares, as well as multiwavelength synchrotron spectra. Our results highlight several limitations of standard post-processing prescriptions for the electron temperature in GRMHD simulations. We also discuss the implications of our results for current and future observations of flares in Sgr A*, M87*, and related systems. Appendices contain detailed convergence studies with respect to resolution and plasma magnetization.

79 ASTRONOMY AND ASTROPHYSICS↗