Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “magnetohydrodynamics”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

An Extension of the Athena++ Code Framework for Radiation-magnetohydrodynamics in General Relativity Using a Finite-solid-angle Discretization

We extend the general-relativistic magnetohydrodynamics (GRMHD) capabilities of Athena++ to incorporate radiation. The intensity field in each finite-volume cell is discretized in angle, with explicit transport in both space and angle properly accounting for the effects of gravity on null geodesics, and with matter and radiation coupled in a locally implicit fashion. Here we describe the numerical procedure in detail, verifying its correctness with a suite of tests. Motivated in particular by black hole accretion in the high-accretion-rate, thin-disk regime, we demonstrate the application of the method to this problem. With excellent scaling on flagship computing clusters, the port of the algorithm to the GPU-enabled AthenaK code now allows the simulation of many previously intractable radiation-GRMHD systems.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Multifrequency General Relativistic Radiation Magnetohydrodynamic Simulations of Thin Disks

We present a set of six general relativistic, multifrequency, radiation magnetohydrodynamic simulations of thin accretion disks with different target mass accretion rates around black holes with spins ranging from nonrotating to rapidly spinning. The simulations use the M 1 closure scheme with 12 independent frequency (or energy) bins ranging logarithmically from 5 × 10− 3 keV to 5 × 10 3 keV. The multifrequency capability allows us to generate crude spectra and energy-dependent light curves directly from the simulations without a need for special postprocessing. While we generally find roughly thermal spectra with peaks around 1–4 keV, our high-spin cases showed harder-than-expected tails for the soft or thermally dominant state. This leads to radiative efficiencies that are up to five times higher than expected for a Novikov–Thorne disk at the same spin. We attribute these high efficiencies to the high-energy, coronal emission. These coronae mostly occupy the effectively optically thin regions near the inner edges of the disks and also cover or sandwich the inner ∼15GM/c 2 of the disks.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

X-Ray Spectra from General Relativistic Radiation Magnetohydrodynamic Simulations of Thin Disks

We compare X-ray emission from several general relativistic, multifrequency, radiation magnetohydrodynamic simulations of thin black hole accretion disks with different accretion rates and spins. The simulations were performed using the M1 closure scheme, resolved with 12 frequency (energy) bins logarithmically spaced from 5 × 10 −3 to 5 × 10 3 ​​​ keV. We apply a general relativistic Monte Carlo transport code to postprocess the simulation data with greater fidelity in frequency resolution and Compton scattering treatment. Despite the relatively few energy bins and Kompaneets approximation to Compton scattering utilized in the M1 method, we find generally good agreement between the methods. Both produce prominent thermal profiles with peaks around 2–2.5 keV, where agreement is particularly strong and representative of the soft state. Both also find weaker (lower-luminosity) thermally sourced emission extending out to 100 keV due to the hotter innermost regions of the disks. Inverse Compton scattering becomes increasingly effective at hardening spectral outputs with increasing black hole spin, and becomes the dominant mechanism for photons that escape with energies between 10 to several hundred keV. At very high rates of spin, the radiation flux in this upscattered component becomes comparable to the thermal flux, a phenomenon typically associated with intermediate states. Beyond 10 4 keV, we observe faint, free–free emission from hot, optically thin coronal regions developing near the horizon, common to both spinning and nonspinning black holes.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Temporal Properties of Compressible Magnetohydrodynamic Turbulence

Describing the temporal properties of compressible magnetohydrodynamic (MHD) turbulence is a fundamental problem that has important implications for particle acceleration and transport in astrophysical plasmas. Here, by carefully analyzing the spatial and temporal properties of compressible MHD turbulence, we derive a new spectral power density function that is supported by simulations. This new function reveals that the low-frequency fluctuations are dominated by modes with small parallel wavenumbers with respect to the mean background magnetic field. Furthermore, for fluctuations with dynamically significant parallel wavenumbers, broadening around their eigenfrequencies is described by this function, which is in close agreement with simulations. We use this formalism to present the scaling properties of individual MHD modes. Such broadening is a direct consequence of nonlinear processes and is different for the three fundamental MHD modes. Our results provide a new window to investigate the temporal properties of turbulence and will enable further studies on the interaction between compressible MHD turbulence and energetic plasmas.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

The Athena++ Adaptive Mesh Refinement Framework: Design and Magnetohydrodynamic Solvers

The design and implementation of a new framework for adaptive mesh refinement calculations are described. It is intended primarily for applications in astrophysical fluid dynamics, but its flexible and modular design enables its use for a wide variety of physics. The framework works with both uniform and nonuniform grids in Cartesian and curvilinear coordinate systems. It adopts a dynamic execution model based on a simple design called a "task list" that improves parallel performance by overlapping communication and computation, simplifies the inclusion of a diverse range of physics, and even enables multiphysics models involving different physics in different regions of the calculation. We describe physics modules implemented in this framework for both nonrelativistic and relativistic magnetohydrodynamics (MHD). These modules adopt mature and robust algorithms originally developed for the Athena MHD code and incorporate new extensions: support for curvilinear coordinates, higher-order time integrators, more realistic physics such as a general equation of state, and diffusion terms that can be integrated with super-time-stepping algorithms. The modules show excellent performance and scaling, with well over 80% parallel efficiency on over half a million threads. The source code has been made publicly available.

79 ASTRONOMY AND ASTROPHYSICS↗
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↗
DOE OSTI2026

Relativistic Magnetohydrodynamic Simulations of Giant Magnetar Bursts

Gradual crustal deformation can generate strongly twisted magnetic fields around magnetars, potentially triggering giant flares with total energies exceeding 10 44 erg. In this letter, we present the first relativistic magnetohydrodynamic simulation of a surface shear-driven magnetar eruption, capturing reconnection-driven plasma heating, the ejection of relativistically hot plasma, and the formation of a hot fireball confined within the inner magnetosphere. We find that magnetic reconnection in the equatorial current sheet launches a hot trailing outflow capable of powering the initial spike observed in giant flares, while simultaneously leaving behind a thermally stratified fireball with sufficient thermal energy to produce the pulsating, decaying tail. Together, these features provide a self-consistent physical framework for understanding the observed energetics of magnetar giant flares. The eruption also expels a magnetically dominated giant plasmoid carrying up to ∼9% of the magnetosphere’s total magnetic energy. Furthermore, our simulation demonstrates how the plasmoid drives the formation of a blast wave—an important ingredient in models linking magnetar eruptions to fast radio bursts.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1962

The Analysis of a Vortex Type Magnetohydrodynamic Induction Generator

Consideration it is given to the performance to the characteristics of an AC magnetohydrodynamic power generator, A rotating magnetic field is imposed on the vortex flow of an electrically conducting fluid, which is injected tangentially into an annulus formed by two nonconducting concentric cylinders and two nonconducting end plates. A perturbation technique is used to determine the two dimensional velocity and three dimensional electromagnetic field and current distributions. Finally, the generated power, the ohmic losses, the effective power and the electrical efficiency of the converter system are calculated.

MAGNETOHYDRODYNAMIC GENERATOR↗