Search NASA⌕ Search

SEARCH · Search NASA

Results for “Wave setup”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Passive Islanding Detection of Inverter-Based Resources in a Noisy Environment

Islanding occurs when a load is energized solely by local generators and can result in frequency and voltage instability, changes in current, and poor power quality. Poor power quality can interrupt industrial operations, damage sensitive electrical equipment, and induce outages upon the resynchronization of the island with the grid. This study proposes an islanding detection method employing a Duffing oscillator to analyze voltage fluctuations at the point of common coupling (PCC) under a high-noise environment. Unlike existing methods, which overlook the noise effect, this paper mitigates noise impact on islanding detection. Power system noise in PCC measurements arises from switching transients, harmonics, grounding issues, voltage sags and swells, electromagnetic interference, and power quality issues that affect islanding detection. Transient events like lightning-induced traveling waves to the PCC can also introduce noise levels exceeding the voltage amplitude by more than seven times, thus disturbing conventional detection techniques. The noise interferes with measurements and increases the nondetection zone (NDZ), causing failed or delayed islanding detection. The Duffing oscillator nonlinear dynamics enable detection capabilities at a high noise level. The proposed method is designed to detect the PCC voltage fluctuations based on the IEEE standard 1547 through the Duffing oscillator. For the voltages beyond the threshold, the Duffing oscillator phase trajectory changes from periodic to chaotic mode and sends an islanded operation command to the inverter. The proposed islanding detection method distinguishes switching transients and faults from an islanded operation. Experimental validation of the method is conducted using a 3.6 kW PV setup.

Energy & Fuels↗
DOE OSTI2023

Inference of the optical depth to reionization τ from Planck CMB maps with convolutional neural networks

The optical depth to reionization, τ, is the least constrained parameter of the cosmological Λ cold dark matter (ΛCDM) model. To date, its most precise value is inferred from large-scale polarized cosmic microwave background (CMB) power spectra from the High Frequency Instrument (HFI) aboard the Planck satellite. These maps are known to contain significant contamination by residual non-Gaussian systematic effects, which are hard to model analytically. Therefore, robust constraints on τ are currently obtained through an empirical cross-spectrum likelihood built from simulations. In this paper, we present a likelihood-free inference of τ from polarized Planck HFI maps which, for the first time, is fully based on neural networks (NNs). NNs have the advantage of not requiring an analytical description of the data and can be trained on state-of-the-art simulations, combining the information from multiple channels. By using Gaussian sky simulations and Planck SRoll2 simulations, including CMB, noise, and residual instrumental systematic effects, we trained, tested, and validated NN models considering different setups. We inferred the value of τ directly from Stokes Q and U maps at ~4° pixel resolution, without computing angular power spectra. On Planck data, we obtained τ NN = 0.0579 ± 0.0082, which is compatible with current EE cross-spectrum results but with a ~30% larger uncertainty, which can be assigned to the inherent nonoptimality of our estimator and to the retraining procedure applied to avoid biases. While this paper does not improve on current cosmological constraints on τ, our analysis represents a first robust application of NN-based inference on real data, and highlights its potential as a promising tool for complementary analysis of near-future CMB experiments, also in view of the ongoing challenge to achieve the first detection of primordial gravitational waves.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Protocols for x-ray transient grating pump/optical probe experiments at x-ray free electron lasers

Abstract Transient grating spectroscopy is a specialized application of the four-wave-mixing methodology and constitutes a versatile technique for investigating the dynamics of vibrational, magnetic and electronic degrees of freedom of matter in a background-free fashion. Recent developments in free-electron laser sources have enabled the extension of this technique into the extreme ultraviolet range. Ongoing efforts to expand transient grating spectroscopy into the x-ray regime promise numerous advantages: (1) substantial penetration depths that allow for probing bulk material properties, (2) element specificity via specific core-excited states, and (3) short wavelengths that allow for excitation gratings with higher momentum transfer and improved spatial resolution. In this study, we comprehensively outline the procedures for conducting x-ray transient grating pump/optical probe experiment. The process encompasses the design and alignment of the experimental setup, as well as the subsequent steps involved in data acquisition and analysis. This paper is intended as a comprehensive guide for researchers interested in implementing x-ray transient grating spectroscopy, providing valuable insights into the intricacies of the experimental workflow required for this novel technique. Furthermore, we discuss the potential for extending this methodology to an x-ray pump/x-ray probe scheme, envisioning a future direction that holds promise for enhancing the capabilities and scope of x-ray transient grating spectroscopy, opening new opportunities for studying ultrafast processes with unprecedented temporal and spatial resolutions.

Optics↗
DOE OSTI2025

Engineering long-lived entanglement through dissipation in quantum hybrid solid-state platforms

Spin squeezing, a form of many-body entanglement, is a crucial resource in quantum metrology and information processing. While experimentally viable protocols for generating stable spin squeezing have been proposed in quantum optics setups, there is growing interest in quantum hybrid solid-state systems as alternative platforms for both engineering and exploring many-body quantum phenomena. In this work, we propose a scheme to generate long-lived spin squeezing in an ensemble of solid-state qubits interacting with electromagnetic noise emitted by a squeezed solid-state bath. We identify the conditions under which quantum correlations within the bath can be transferred to the qubit array, driving it into an entangled state independently of its initial configuration. To assess the experimental feasibility of our approach, we analyze the dynamics of an array of solid-state spin defects coupled to a common ferromagnetic bath, which is driven into a non-equilibrium squeezed state through its interaction with a surface acoustic wave mode. Our results demonstrate that the ensemble can exhibit steady-state spin squeezing under suitable conditions, opening new pathways for the generation of robust many-body entanglement in solid-state spin ensembles.

NV centers↗
DOE OSTI2022

Investigating growth to detonation in vapor-deposited hexanitrostilbene and pentaerythritol tetranitrate films using high-throughput methods

In this work, a high-throughput experimental setup was used to characterize initiation threshold and growth to detonation in the explosives hexanitrostilbene (HNS) and pentaerythritol tetranitrate (PETN). The experiment sequentially launched an array of laser-driven flyers to shock samples arranged in a 96-well microplate geometry, with photonic Doppler velocimetry diagnostics to characterize flyer velocity and particle velocity at the explosive–substrate interface. Vapor-deposited films of HNS and PETN were used to provide numerous samples with various thicknesses, enabling characterization of the evolution of growth to detonation. One-dimensional hydrocode simulations were performed with reactions disabled to illustrate where the experimental data deviate from the predicted inert response. Prompt initiation was observed in 144 μm thick HNS films at flyer velocities near 3000 m/s and in 125 μm thick PETN films at flyer velocities near 2400 m/s. This experimental setup enables rapid quantification of the growth of reactions in explosive materials that can reach detonation at sub-millimeter length scales. These data can subsequently be used for parameterizing reactive burn models in hydrocode simulations, as discussed in Paper II [D. E. Kittell, R. Knepper, and A. S. Tappan, J. Appl. Phys. 131, 154902 (2022)].

36 MATERIALS SCIENCE↗
DOE OSTI2021

Design and Testing of a Water-Cooled Rotating Detonation Combustor at Elevated Operating Pressures - Abstract

This is an extended abstract being submitted for consideration that briefly describes the material that will be included in the final paper. This study will detail the design and testing of a water-cooled Rotating Detonation Combustor (RDC) that permits operation for extended periods of time ensuring that the device has reach a stable operating temperature. Extended run times also provides an opportunity to consider transient behaviors that occur as a result of altering the operating conditions such as equivalence ratio. The experimental setup is also unique in that it consists of a ducted exhaust with a downstream high-temperature valve that can control the pre-combustion pressure in the RDC independent of the combustor annulus and exit geometry. This paper will examine the results from the extended operation of the water-cooled RDC over a range of both transient and steady state equivalence ratios (0.5 - 1.0), pre-combustion operating pressures (0 - 207 kPa), mass flow rates (0.42 – 0.5 kg/sec) and air inlet throat area to combustor channel area ratios (0.09 – 0.32). All tests were performed while operating on hydrogen in air at a nominal air inlet temperature of 340 – 350 K. Limited data obtained at an elevated inlet air temperature of 480 K will also be discussed. Results suggest that under certain conditions the operating mode of the detonation wave can vary as it approaches a stable operating temperature. For the RDC tested, a stable operating temperature is not achieved until approximately 8-10 seconds of run time.

Ferguson, Donald↗
DOE OSTI2020

Determination of acoustic nonlinearity parameters using thermal modulation of ultrasonic waves

This study presents a test method and its theoretical framework to determine the acoustic nonlinearity parameters (α, β, δ) of material using thermal modulation of ultrasonic waves. Temperature change-induced thermal strain excites the nonlinear response of the material and modulates the ultrasonic wave propagating in it. Experimental results showed a strong correlation between the relative wave velocity change and the temperature change. With a quadratic polynomial model, the acoustic nonlinearity parameters were obtained from the polynomial coefficients by curve fitting the experimental data. Their effects on thermal-induced velocity change were discussed. The parameters α, β, and δ govern the hysteretic gap, average slope, and curvature of the correlation curve, respectively. The proposed theory was validated on aluminum, steel, intact and damaged concrete samples. The obtained nonlinear parameters show reasonable agreement with values reported in the literature. Compared to other nonlinear acoustic methods using vibration or acoustic excitation, the thermal modulation method generates more uniform, slow changing, and larger strain field in the test sample. Employing the thermal effect as the driving force for nonlinearity instead of an undesired influencing factor, this method can measure the absolute values of α, β, and δ with good accuracy using a simple ultrasonic test setup.

Physics↗
DOE OSTI2022

Dynamics in a stellar convective layer and at its boundary: Comparison of five 3D hydrodynamics codes

Our ability to predict the structure and evolution of stars is in part limited by complex, 3D hydrodynamic processes such as convective boundary mixing. Hydrodynamic simulations help us understand the dynamics of stellar convection and convective boundaries. However, the codes used to compute such simulations are usually tested on extremely simple problems and the reliability and reproducibility of their predictions for turbulent flows is unclear. We define a test problem involving turbulent convection in a plane-parallel box, which leads to mass entrainment from, and internal-wave generation in, a stably stratified layer. We compare the outputs from the codes FLASH, MUSIC, PPMSTAR, PROMPI, and SLH, which have been widely employed to study hydrodynamic problems in stellar interiors. The convection is dominated by the largest scales that fit into the simulation box. All time-averaged profiles of velocity components, fluctuation amplitudes, and fluxes of enthalpy and kinetic energy are within ≲3σ of the mean of all simulations on a given grid (128 3 and 256 3 grid cells), where σ describes the statistical variation due to the flow’s time dependence. They also agree well with a 512 3 reference run. The 128 3 and 256 3 simulations agree within 9% and 4%, respectively, on the total mass entrained into the convective layer. The entrainment rate appears to be set by the amount of energy that can be converted to work in our setup and details of the small-scale flows in the boundary layer seem to be largely irrelevant. Our results lend credence to hydrodynamic simulations of flows in stellar interiors. We provide in electronic form all outputs of our simulations as well as all information needed to reproduce or extend our study.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Power transport efficiency during O-X-B 2nd harmonic electron cyclotron heating in a helicon linear plasma device

The principal objective of this work is to report on the power coupled to a tungsten target in the prototype-material plasma exposure experiment device during oblique injection of a microwave beam (<70 kW at 28 GHz) into a high-power (~100 kW at 13.56 MHz) over-dense (${n_e} > 1 \times {10^{19}}{\text{ }}{{\text{m}}^{ - 3}}$) deuterium helicon plasma column. The experimental setup, electron heating system, electron heating scheme, and IR thermographic diagnostic for quantifying the power transport is described in detail. It is demonstrated that the power transported to the target can be effectively controlled by adjusting the magnetic field profile. Using this method, heat fluxes up to 22 MW m -2 and power transport efficiencies in the range of 17%–20% have been achieved using 70 kW of microwave power. It is observed that most of the heat flux is confined to a narrow region at the plasma periphery. Ray-tracing calculations are presented which indicate that the power is coupled to the plasma electrons via an O-X-B mode conversion process. Calculations indicate that the microwave power is absorbed in a single pass at the plasma periphery via collisions and in the over-dense region via 2nd harmonic cyclotron resonance of the electron Bernstein wave. The impact of our results is discussed in the context of MPEX.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

GR-Athena++: Puncture Evolutions on Vertex-centered Oct-tree Adaptive Mesh Refinement

Numerical relativity is central to the investigation of astrophysical sources in the dynamical and strong-field gravity regime, such as binary black hole and neutron star coalescences. Current challenges set by gravitational-wave and multimessenger astronomy call for highly performant and scalable codes on modern massively parallel architectures. We present GR-Athena++, a general-relativistic, high-order, vertex-centered solver that extends the oct-tree, adaptive mesh refinement capabilities of the astrophysical (radiation) magnetohydrodynamics code Athena++. To simulate dynamical spacetimes, GR-Athena++ uses the Z4c evolution scheme of numerical relativity coupled to the moving puncture gauge. We demonstrate stable and accurate binary black hole merger evolutions via extensive convergence testing, cross-code validation, and verification against state-of-the-art effective-one-body waveforms. GR-Athena++ leverages the task-based parallelism paradigm of Athena++ to achieve excellent scalability. We measure strong-scaling efficiencies above 95% for up to ~1.2 × 10 4 CPUs and excellent weak scaling is shown up to ~10 5 CPUs in a production binary black hole setup with adaptive mesh refinement. GR-Athena++ thus allows for the robust simulation of compact binary coalescences and offers a viable path toward numerical relativity at exascale.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Simulink Model of a New Energy EVG-005 5kW Hydropower Turbine

Simulink model for a New Energy 5kW hydropower turbine. ADCP data ("ds_streamwise_7_13.nc") and DC voltage, DC current, and rotor rotation observed from the New Energy EVG-005 Current Energy Converter (CEC) ("electrical_7_13_10ohms.nc") were collected at the Tanana River Test Site (TRTS) near Nenana Alaska. - Data was collected on July 13th, 2023 with a constant 10 ohms resistance applied with a DC load bank. - Simulink model is meant to resemble the electrical setup at the TRTS. - Model is initialized by running the "NewEnergy_2023_10hz.m" Matlab script. Then the Simulink model ("New_Energy_Model_PMSM.slx") can be run. - Results are processed with the "NewEnergy_2023_processing.m" Matlab code. The Matlab results are also saved in the "New_Energy_7_13_model_validation_results.mat" Matlab file. This can be directly loaded into the Matlab Workspace using the "Load()" command. - The timetable variables "electrical_model_downsampled" is the model results and the "electrical_limited" is the experimental data from the TRTS.

16 TIDAL AND WAVE POWER↗
DOE OSTI2021

Learning the ground state of a non-stoquastic quantum Hamiltonian in a rugged neural network landscape

Strongly interacting quantum systems described by non-stoquastic Hamiltonians exhibit rich low-temperature physics. Yet, their study poses a formidable challenge, even for state-of-the-art numerical techniques. Here, we investigate systematically the performance of a class of universal variational wave-functions based on artificial neural networks, by considering the frustrated spin- 1/2 1 / 2 J_1-J_2 J 1 − J 2 Heisenberg model on the square lattice. Focusing on neural network architectures without physics-informed input, we argue in favor of using an ansatz consisting of two decoupled real-valued networks, one for the amplitude and the other for the phase of the variational wavefunction. By introducing concrete mitigation strategies against inherent numerical instabilities in the stochastic reconfiguration algorithm we obtain a variational energy comparable to that reported recently with neural networks that incorporate knowledge about the physical system. Through a detailed analysis of the individual components of the algorithm, we conclude that the rugged nature of the energy landscape constitutes the major obstacle in finding a satisfactory approximation to the ground state wavefunction, and prevents learning the correct sign structure. In particular, we show that in the present setup the neural network expressivity and Monte Carlo sampling are not primary limiting factors.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Toward UV models of kinetic mixing and portal matter. VI. A more complex dark matter sector?

Portal matter (PM), having both Standard Model (SM) and dark sector charges, can induce kinetic mixing between the 𝑈⁢(1) 𝐷 dark photon and the SM gauge fields at the 1-loop level offering an attractive mechanism by which light (≲1 GeV) thermal dark matter (DM) can interact with visible matter and obtain its observed relic density. In doing so, if the DM is fermionic, the CMB and other astrophysical observations inform us that it must be Majorana/pseudo-Dirac in nature to avoid velocity/temperature-independent 𝑠-wave annihilation to SM final states. How does this idea fit into a more UV-complete picture also including the SM interactions? There are some reasons to believe that at least a first step along this path may not lie too far away in energy due to the renormalization group equations running of the dark gauge coupling, which for a significant range of parameters, becomes nonperturbative at/before the ∼10’s of TeV energy range. This implies that 𝑈⁢(1) 𝐷 must become embedded in an asymptotically free, non-Abelian group, 𝐺 𝐷 , before this can occur. The breaking of this larger group then produces the masses for the PM and the additional gauge fields associated with 𝐺 𝐷 then can lead to new interactions between the SM and the dark sector. Following several bottom-up approaches, we have examined a set of distinctive and testable phenomenological features associated with this general setup, based upon a number of simplifying assumptions. Clearly, it behooves us to explore the impact of these specific assumptions on these predictions for the array of possible experimental tests of this class of models. In most past analyses it has been assumed that DM is a vectorlike, complex singlet under the group 𝐺 𝐷 . If this assumption is relaxed, the dark sector must be augmented by additional fermion(s) and the associated scalar fields needed to break the gauge symmetries while generating the needed Majorana-like mass terms for the DM. In this paper, we analyze the simplest extension of this kind wherein the DM lies in a vectorlike doublet of 𝐺 𝐷 , which we take to have the structure 𝑆⁢𝑈⁢(2) 𝐼 ×𝑈⁢(1) 𝑌 𝐼 as in earlier work, leading to new phenomenological implications. We find, for example, that given the current LHC search constraints on the masses of heavy gauge bosons, the production of these new dark states with large rates is unlikely to occur at colliders unless they are produced singly in 𝑔⁡𝑔 fusion or their pair production cross sections are resonantly enhanced. Here, we also find that an additional mechanism arises to generate hierarchal neutrino masses in such a setup.

Extensions of Higgs sector↗
DOE OSTI2021

Study of the Effect of Reflections on High-Power, 110-GHz Pulsed Gyrotron Operation

The effect of reflection is studied experimentally and theoretically on a high-power 110-GHz gyrotron operating in the TE 22,6 mode in 3 μs pulses at 96 kV, 40 A. The experimental setup allows variation of the reflected power from 0 to 33% over a range of gyrotron operating conditions. The phase of the reflection is varied by translating the reflector along the axis. Operating at a higher efficiency point, at 4.40 T with 940kW of output power, reflected power exceeding 11% causes a switch from operation in the TE 22,6 to simultaneous operation in the TE 22,6 and TE 21,6 modes with a large decrease of the total gyrotron output power. This switching effect is in good agreement with simulations using the MAGY code. Operating at a more stable point, 4.44 T with 580kW of output power, when the reflection is increased, the output power remains in the TE 22,6 mode but it decreases monotonically with increasing reflection, dropping to 200kW at 33% reflection. Furthermore, at a reflection above 22%, a power modulation at 25 to 30MHz is observed, independent of the phase of the reflected wave. Such a modulated signal may be useful in spectroscopic and other applications.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Two-nucleon 𝑆-wave interactions at the SU(3) flavor-symmetric point with 𝑚 𝑢⁢𝑑 ≈ 𝑚$^{phys}_𝑠$: A first lattice QCD calculation with the stochastic Laplacian Heaviside method

We report on the first application of the stochastic Laplacian Heaviside method for computing multiparticle interactions with lattice QCD to the two-nucleon system. Like the Laplacian Heaviside method, this method allows for the construction of interpolating operators which can be used to construct a set of positive-definite two-nucleon correlation functions, unlike nearly all other applications of lattice QCD to two nucleons in the literature. It also allows for a variational analysis in which optimal linear combinations of the interpolating operators are formed that couple predominantly to the eigenstates of the system. Utilizing such methods has become of paramount importance to help resolve the discrepancy in the literature on whether two nucleons in either isospin channel form a bound state at pion masses heavier than physical, with the discrepancy persisting even in the SU(3)-flavor-symmetric point with all quark masses near the physical strange quark mass. This is the first in a series of papers aimed at resolving this discrepancy. In the present work, we employ the stochastic Laplacian Heaviside method without a hexaquark operator in the basis at a lattice spacing of 𝑎 ≈0.086 fm, lattice volume of 𝐿 = 48⁢𝑎 ≈ 4.1 fm and pion mass 𝑚 𝜋 ≈ 714 MeV. With this setup, the observed spectrum of two-nucleon energy levels strongly disfavors the presence of a bound state in either the deuteron or dineutron channel.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Soft x-ray high-harmonic generation in an anti-resonant hollow core fiber driven by a 3 μ m ultrafast laser

High-harmonic upconversion driven by a mid-infrared femtosecond laser can generate coherent soft x-ray beams in a tabletop-scale setup. Here, we report on a compact ytterbium-pumped optical parametric chirped pulse amplifier (OPCPA) laser system seeded by an all-fiber front-end and employing periodically poled lithium niobate (PPLN) nonlinear media operated near the pulse fluence limits of current commercially available PPLN crystals. The OPCPA delivers 3 µm wavelength pulses with 775 µJ energy at 1 kHz repetition rate, with transform-limited 120 fs pulse duration, diffraction-limited beam quality, and ultrahigh 0.33% rms energy stability over >18 h. Using this laser, we generate soft x-ray high harmonics (HHG) in argon gas by focusing into a low-loss, high-pressure gas-filled anti-resonant hollow core fiber (ARHCF), generating coherent light at photon energies up to the argon L-edge (250 eV) and carbon K-edge (284 eV), with high beam quality and ∼1% rms energy stability. This work demonstrates soft x-ray HHG in a high-efficiency guided-wave phase matched geometry, overcoming the high losses inherent to mid-IR propagation in unstructured waveguides, or the short interaction lengths of gas cells or jets. The ARHCF can operate in the long term without damage and with the repetition rate, stability, and robustness required for demanding applications in spectromicroscopy and imaging. Finally, we discuss routes for further optimizing the soft x-ray HHG flux by driving He at higher laser intensities using either the signal (1.5 μm) or idler wavelengths (3 μm).

Femtosecond lasers↗
DOE OSTI2022

Embedded, Real-Time, and Distributed Traveling Wave Fault Location Method Using Graph Convolutional Neural Networks

This work proposes and develops an implementation of a fault location method to provide a fast and resilient protection scheme for power distribution systems. The method analyzes the transient dynamics of traveling waves (TWs) to generate features using the discrete wavelet transform (DWT), which are then used to train several graph convolutional network (GCN) models. Faults are simulated in the IEEE 34-node system, which is divided into three protection zones (PZs). The goal is to identify the PZ in which the fault occurs. The GCN models create a distributed protection scheme, as all nodes are able to retrieve a prediction. Given that message-passing between nodes occurs both during training and in the execution of the model, the resiliency of such schemes to communication losses was analyzed and demonstrated. One of the models, which only uses voltage measurements, was implemented on a Texas Instruments F28379D development board. The execution times were monitored to assess the speed of the protection scheme. It is shown that the proposed method can be executed in approximately a millisecond, which is comparable to existing TW protection in the transmission system. For experimental purposes, a DWT-based detection method is employed. A design of a setup to playback TWs using two development boards is also addressed.

Jiménez-Aparicio, Miguel (ORCID:000000016864461X)↗
DOE OSTI2025

GR-Athena++: General-relativistic Magnetohydrodynamics Simulations of Neutron Star Spacetimes

We present the extension of GR-Athena++ to general-relativistic magnetohydrodynamics (GRMHD) for applications to neutron star spacetimes. The new solver couples the constrained transport implementation of Athena++ to the Z4c formulation of the Einstein equations to simulate dynamical spacetimes with GRMHD using oct-tree adaptive mesh refinement. We consider benchmark problems for isolated and binary neutron star spacetimes demonstrating stable and convergent results at relatively low resolutions and without grid symmetries imposed. The code correctly captures magnetic field instabilities in nonrotating stars with total relative violation of the divergence-free constraint of 10 –16 . It handles evolutions with a microphysical equation of state and black hole formation in the gravitational collapse of a rapidly rotating star. For binaries, we demonstrate correctness of the evolution under the gravitational radiation reaction and show convergence of gravitational waveforms. We showcase the use of adaptive mesh refinement to resolve the Kelvin–Helmholtz instability at the collisional interface in a merger of magnetised binary neutron stars. GR-Athena++ shows strong scaling efficiencies above 80% in excess of 10 5 CPU cores and excellent weak scaling is shown up to ~5 × 10 5 CPU cores in a realistic production setup. GR-Athena++ allows for the robust simulation of GRMHD flows in strong and dynamical gravity with exa-scale computers.

79 ASTRONOMY AND ASTROPHYSICS↗