Search NASA⌕ Search

SEARCH · Search NASA

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

2D Magnetohydrodynamic Simulations of the Electrothermal Instability in Metallic Liners

The Virginia Tech (VT) Plasma Dynamics Laboratory Computational (PDCL) and Lawrence Livermore National Laboratory (LLNL) are performing two dimensional (2D) simulations of the electrothermal instability (ETI) using the LLNL multi-physics code Ares. Ares is a multi-physics arbitrary–Lagrangian-Eulerian (ALE) code developed by LLNL and is of particular use in studying magnetohydrodynamic (MHD) instabilities like the ETI due to its resistive MHD, magnetic diffusion, and radiative-hydrodynamics packages. Among its capabilities, it has the ability to model material strength, perform adaptive mesh refinement (AMR), and incorporate a wide variety of equations of state models and conductivity models. The 2D Ares simulation model created by VT-LLNL for studying the development and growth of the electrothermal instability has been configured with initial conditions based on the Mykonos Electrothermal Instability II (METI-II) experiments described by this grant and conducted by team members at the University of Nevada (UNR), the University of New Mexico (UNM), and Sandia National Laboratories. Previously, preliminary 2D Ares simulations of the ETI had been run to approximately 80ns. The rods in these simulations were initiated with sinusoidal perturbations at a similar order of magnitude to those measured on the aluminum rods used for the Mykonos experiment. This model has been improved by increasing the spatial resolution of the simulations and running the simulations further in time. In addition to the simulation run-times extending, the preliminary sinusoidal perturbation has been replaced with a perturbation derived from amplitude measurements by the experimental team, thereby correlating the simulation inputs better to the experimental runs. These new runs are capable of reaching 120ns of simulated time for the uncoated cases and to 200ns the 41 μm coated cases.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Structure Formation Through Magnetohydrodynamical Instabilities in Primordial Disks

The shear flow instabilities under the presence of magnetic fields in the primordial disk can greatly facilitate the formation of density structures that serve as seeds prior to the onset of the gravitational Jeans instability. We evaluate the effects of the Parker, magnetorotational and kinematic dynamo instabilities by comparing the properties of these instabilities. We calculate the mass spectra of coagulated density structures by the above mechanism in the radial direction for an axisymmetric magnetohydrodynamic (MHD) torus equilibrium and power density profile models. Our local three-dimensional MHD simulation indicates that the coupling of the Parker and magnetorotational instability creates spiral arms and gas blobs in an accretion disk, reinforcing the theory and model. Such a mechanism for the early structure formation may be tested in a laboratory. The recent progress in experiments involving shear flows in rotating tokamak, field reversed configuration (FRC) and laser plasmas may become a key element to advance in nonlinear studies.

Noguchi, Koichi↗
DOE OSTI2020

MOCMC: Method of Characteristics Moment Closure, a Numerical Method for Covariant Radiation Magnetohydrodynamics

In this work, we present a conservative numerical method for radiation magnetohydrodynamics with frequency-dependent full transport in stationary spacetimes. This method is stable and accurate for both large and small optical depths and radiation pressures. The radiation stress–energy tensor is evolved in flux-conservative form, and closed with a swarm of samples that each transport a multigroup representation of the invariant specific intensity along a null geodesic. In each zone, the enclosed samples are used to efficiently construct a Delaunay triangulation of the unit sphere in the comoving frame, which in turn is used to calculate the Eddington tensor, average source terms, and adaptively refine the sample swarm. Furthermore, radiation four-forces are evaluated in the moment sector in a semi-implicit fashion. The radiative transfer equation is solved in invariant form deterministically for each sample. Since each sample carries a discrete representation of the full spectrum, the cost of evaluating the transport operator is independent of the number of frequency groups, representing a significant reduction of algorithmic complexity for transport in frequency-dependent problems. The major approximation we make in this work is performing scattering in an angle-averaged way. Local adaptivity in samples also makes this scheme more amenable to nonuniform meshes than a traditional Monte Carlo method. We describe the method and present results on a suite of test problems. We find that Method of Characteristics Moment Closure converges at least as ~N -1 , rather than the canonical Monte Carlo N -1/2 , where N is the number of samples per zone.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Full Transport General Relativistic Radiation Magnetohydrodynamics for Nucleosynthesis in Collapsars

We model a compact black hole-accretion disk system in the collapsar scenario with full transport, frequency dependent, general relativistic radiation magnetohydrodynamics. We examine whether or not winds from a collapsar disk can undergo rapid neutron capture (r-process) nucleosynthesis and significantly contribute to solar r-process abundances. We find the inclusion of accurate transport has significant effects on outflows, raising the electron fraction above ${Y}_{{\rm{e}}}\sim 0.3$ and preventing third-peak r-process material from being synthesized. We analyze the time evolution of neutrino processes and electron fraction in the disk and present a simple one-dimensional model for the vertical structure that emerges. Finally, we compare our simulation to semi-analytic expectations and argue that accurate neutrino transport and realistic initial and boundary conditions are required to capture the dynamics and nucleosynthetic outcome of a collapsar.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Small-amplitude Compressible Magnetohydrodynamic Turbulence Modulated by Collisionless Damping in Earth’s Magnetosheath: Observation Matches Theory

Plasma turbulence is a ubiquitous dynamical process that transfers energy across many spatial and temporal scales and affects energetic particle transport. Recent advances in the understanding of compressible magnetohydrodynamic (MHD) turbulence demonstrate the important role of damping in shaping energy distributions on small scales, yet its observational evidence is still lacking. This study provides the first observational evidence of substantial collisionless damping (CD) modulation on the small-amplitude compressible MHD turbulence cascade in Earth’s magnetosheath using four Cluster spacecraft. Based on an improved compressible MHD decomposition algorithm, turbulence is decomposed into three eigenmodes: incompressible Alfvén modes and compressible slow and fast (magnetosonic) modes. Our observations demonstrate that CD enhances the anisotropy of compressible MHD modes because CD has a strong dependence on wave propagation angle. The wavenumber distributions of slow modes are mainly stretched perpendicular to the background magnetic field (B 0 ) and weakly modulated by CD. In contrast, fast modes are subjected to a more significant CD modulation. Fast modes exhibit a weak, scale-independent anisotropy above the CD truncation scale. Below the CD truncation scale, the anisotropy of fast modes enhances as wavenumbers increase. As a result, fast-mode fractions in the total energy of compressible modes decrease with the increase of perpendicular wavenumber (to B 0 ) or wave propagation angle. Our findings reveal how the turbulence cascade is shaped by CD and its consequences for anisotropies in the space environment.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Magnetic Flux Plays an Important Role during a Black Hole X-Ray Binary Outburst in Radiative Two-temperature General Relativistic Magnetohydrodynamic Simulations

Abstract Black hole (Bh) X-ray binaries cycle through different spectral states of accretion over the course of months to years. Although persistent changes in the Bh mass accretion rate are generally recognized as the most important component of state transitions, it is becoming increasingly evident that magnetic fields play a similarly important role. In this article, we present the first radiative two-temperature general relativistic magnetohydrodynamics simulations in which an accretion disk transitions from a quiescent state at an accretion rate of M ̇ ∼ 10 − 10 M ̇ Edd to a hard-intermediate state at an accretion rate of M ̇ ∼ 10 − 2 M ̇ Edd . This huge parameter space in mass accretion rate is bridged by artificially rescaling the gas density scale of the simulations. We present two jetted BH models with varying degrees of magnetic flux saturation. We demonstrate that in “standard and normal evolution” models, which are unsaturated with magnetic flux, the hot torus collapses into a thin and cold accretion disk when M ̇ ≳ 5 × 10 − 3 M ̇ Edd . On the other hand, in “magnetically arrested disk” models, which are fully saturated with vertical magnetic flux, the plasma remains mostly hot with substructures that condense into cold clumps of gas when M ̇ ≳ 1 × 10 − 2 M ̇ Edd . This suggests that the spectral signatures observed during state transitions are closely tied to the level of magnetic flux saturation.

Astronomy & Astrophysics↗
DOE OSTI2024

Cholla-MHD: An Exascale-capable Magnetohydrodynamic Extension to the Cholla Astrophysical Simulation Code

Abstract We present an extension of the massively parallel, GPU native, astrophysical hydrodynamics code Cholla to magnetohydrodynamics (MHD). Cholla solves the ideal MHD equations in their Eulerian form on a static Cartesian mesh utilizing the Van Leer + constrained transport integrator, the HLLD Riemann solver, and reconstruction methods at second and third order. Cholla’s MHD module can perform ≈260 million cell updates per GPU-second on an NVIDIA A100 while using the HLLD Riemann solver and second order reconstruction. The inherently parallel nature of GPUs combined with increased memory in new hardware allows Cholla’s MHD module to perform simulations with resolutions ∼500 3 cells on a single high-end GPU (e.g., an NVIDIA A100 with 80 GB of memory). We employ GPU direct Message Passing Interface to attain excellent weak scaling on the exascale supercomputer Frontier, while using 74,088 GPUs and simulating a total grid size of over 7.2 trillion cells. A suite of test problems highlights the accuracy of Cholla’s MHD module and demonstrates that zero magnetic divergence in solutions is maintained to round off error. We also present new testing and CI tools using GoogleTest, GitHub Actions, and Jenkins that have made development more robust and accurate and ensure reliability in the future.

Astronomy & Astrophysics↗
DOE OSTI2024

How to Turn Jets into Cylinders near Supermassive Black Holes in 3D General Relativistic Magnetohydrodynamic Simulations

Abstract Accreting supermassive black holes (SMBHs) produce highly magnetized relativistic jets that tend to collimate gradually as they propagate outward. However, recent radio interferometric observations of the 3C 84 galaxy reveal a stunning, cylindrical jet already at several hundred SMBH gravitational radii,r≳ 350r g . We explore how such extreme collimation emerges via a suite of 3D general relativistic magnetohydrodynamic simulations. We consider an SMBH surrounded by a magnetized torus immersed in a constant-density ambient medium that starts at the edge of the SMBH sphere of influence, chosen to be much larger than the SMBH gravitational radius,r B = 10 3 r g . We find that radiatively inefficient accretion flows (e.g., M87) produce winds that collimate the jets into parabolas near the black hole. After the disk winds stop collimating the jets atr≲r B , they turn conical. Once outsider B , the jets run into the ambient medium and form backflows that collimate the jets into cylinders some distance beyondr B . Interestingly, for radiatively efficient accretion, as in 3C 84, the radiative cooling saps the energy out of the disk winds; at early times, they cannot efficiently collimate the jets, which skip the initial parabolic collimation stage, start out conical near the SMBH, and turn into cylinders already atr≃ 300r g , as observed in 3C 84. Over time, the jet power remains approximately constant, whereas the mass accretion rate increases; the winds grow in strength and start to collimate the jets, which become quasi-parabolic near the base, and the transition point to a nearly cylindrical jet profile moves outward while remaining insider B .

Astronomy & Astrophysics↗
DOE OSTI2025

Magnetically Driven Neutron-rich Ejecta Unleashed: Global 3D Neutrino–General Relativistic Magnetohydrodynamic Simulations of Collapsars Probe the Conditions for r -process Nucleosynthesis

Abstract Collapsars—rapidly rotating stellar cores that form black holes—can power gamma-ray bursts and are proposed to be key contributors to the production of heavy elements in the Universe via the rapid neutron capture process ( r -process). Previous neutrino-transport collapsar simulations have been unable to unbind neutron-rich material from the disk. However, these simulations have not included sufficiently strong magnetic fields and the black hole (BH), both of which are essential for launching mass outflows. We present ν h-amr , a novel neutrino-transport general relativistic magnetohydrodynamic ( ν GRMHD) code, which we use to perform the first 3D global ν GRMHD collapsar simulations. We find a self-consistent formation of a weakly magnetized dense accretion disk, which has sufficient time to neutronize. Eventually, substantial magnetic flux accumulates near the BH, becomes dynamically important, leads to a magnetically arrested disk (MAD), and unbinds some of the neutron-rich material. However, the strong flux also hinders accretion, lowers density, and increases neutrino-cooling timescale, which prevents further disk neutronization. Typical collapsar progenitors with mass accretion rates, M ̇ ∼ 0.1 − 1 M ⊙ s - 1 , do not produce significant neutron-rich ( Y e < 0.25) ejecta. However, we find that MADs at higher mass accretion rates, M ̇ ≳ few M ⊙ s - 1 (e.g., for more centrally concentrated progenitors), can unbind M ej ≲ M ⊙ of neutron-rich ejecta. The outflows inflate a shocked cocoon that mixes with the infalling neutron-poor stellar gas and raises the final outflow Y e ; however, the final r -process yield may be determined earlier at the point of neutron capture freeze-out. Future work will explore under what conditions more typical collapsar engines become r -process factories.

Issa, Danat (ORCID:0009000524787631)↗
DOE OSTI2025

Self-similar Cosmic-Ray Transport in High-resolution Magnetohydrodynamic Turbulence

We study the propagation of cosmic rays (CRs) through a simulation of magnetohydrodynamic (MHD) turbulence at an unprecedented resolution of 10,240 3 . We drive turbulence that is subsonic and super-Alfvénic, characterized by δB rms /B 0 = 2. The high resolution enables an extended inertial range such that the Alfvén scale l A , where δB(l A ) ≈ B 0 , is well resolved. This allows us to properly capture how the cascade transitions from large amplitudes on large scales to small amplitudes on small scales. We find that sharp bends in the magnetic field are key mediators of particle transport even on small scales via resonant curvature scattering. We further find that particle scattering in the turbulence shows strong hints of self-similarity: (1) the diffusion has a weak energy dependence over almost two decades in particle energy, and (2) the particles’ random walk exhibits a broad power-law distribution of collision times such that the diffusion is dominated by the rarest, long-distance excursions. Our results suggest that large-amplitude MHD turbulence can provide efficient scattering over a wide range of CR energies and may help explain many CR observations above a ∼TeV: the flattening of the B/C spectrum, the hardening of CR primary spectra, and the weak dependence of arrival anisotropy on CR energy.

Kempski, Philipp [Princeton Univ., NJ (United Stat↗
NASA NTRS1968

Study of convective magnetohydrodynamic channel flow

Study involves the effects of the interactions of electromagnetic, velocity, and temperature fields to aid in the design of a magnetohydrodynamic device. It concerns a theoretical analysis of the convective flow of an electrically conducting gas in a channel composed of conducting walls.

Singer, R. M.↗
NASA NTRS1969

Magnetohydrodynamic generators using two-phase liquid-metal flows

Two-phase flow generator cycle of a magnetohydrodynamic /MHD/ generator uses a working fluid which is compressible and treated as an expanding gas. The two-phase mixture passes from the heat source through the MHD generator, where the expansion process takes place and the electrical energy is extracted.

Petrick, M.↗
NASA NTRS1969

Studies of cycles for liquid-metal magnetohydrodynamic generation of power

Studies of liquid-metal magnetohydrodynamic power cycles indicate that the overall efficiency of a binary cycle, employing a liquid-metal topping cycle and a bottoming steam cycle, may reach 60 percent. Details of analyses and data on cycles are presented, and the commercial potential of the binary cycle is discussed.

Lee, K.↗