SEARCH · Engineering Papers

Results for “magnetohydrodynamical 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 127 records · Page 7

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Relative alignment between magnetic fields and molecular gas structure in molecular clouds

Here, we compare the structure of synthetic dust polarization with synthetic molecular line emission from radiative transfer calculations using a three-dimensional, turbulent collapsing-cloud magnetohydrodynamics simulation. The histogram of relative orientation (HRO) technique and the projected Rayleigh statistic (PRS) are considered. In our trans-Alfvénic (more strongly magnetized) simulation, there is a transition to perpendicular alignment at densities above ~4 × 10 3 cm –3 . This transition is recovered in most of our synthetic observations of optically thin molecular tracers; however, for 12 CO it does not occur and the PRS remains in parallel alignment across the whole observer space. We calculate the physical depth of the optical depth τ = 1 surface and find that for 12 CO it is largely located in front of the cloud mid-plane, suggesting that 12 CO is too optically thick and instead mainly probes low-volume density gas. In our super-Alfvénic simulation, the magnetic field becomes significantly more tangled, and all observed tracers tend towards no preference for perpendicular or parallel alignment. An observable difference in alignment between optically thin and optically thick tracers may indicate the presence of a dynamically important magnetic field, though there is some degeneracy with viewing angle. We convolve our data with a Gaussian beam and compare it with HRO results of the Vela C molecular cloud. We find good agreement between these results and our sub-Alfvénic simulations when viewed with the magnetic field in the plane of the sky (especially when sensitivity limitations are considered), though the observations are also consistent with an intermediately inclined magnetic field.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

The three-dimensional structure of black hole accretion flows within the plunging region

ABSTRACT We analyse, using new analytical models and numerical general relativistic magnetohydrodynamic simulations, the three-dimensional properties of accretion flows inside the plunging region of black hole spacetimes (i.e. at radii smaller than the innermost stable circular orbit). These simulations are of thick discs, with aspect ratios of order unity $h/r \sim 1$, and with a magnetic field geometry given by the standard low-magnetization ‘SANE’ configuration. This work represents the first step in a wider analysis of this highly relativistic region. We show that analytical expressions derived in the ‘thin disc’ limit describe the numerical results remarkably well, despite the large aspect ratio of the flow. We further demonstrate that accretion within this region is typically mediated via spiral arms, and that the geometric properties of these spiral structures can be understood with a simple analytical model. These results highlight how accretion within the plunging region is fundamentally two-dimensional in character, which may have a number of observational implications. We derive a modified theoretical description of the pressure within the plunging region which accounts for turbulent heating and may be of use to black hole image modelling.

Astronomy & Astrophysics↗

DOE OSTI2024

The influence of the cloud virial parameter on the initial mass function

ABSTRACT Crucial for star formation is the interplay between gravity and turbulence. The observed cloud virial parameter, $\alpha _{\mathrm{vir}}$, which is the ratio of twice the turbulent kinetic energy to the gravitational energy, is found to vary significantly in different environments, where the scatter among individual star-forming clouds can exceed an order of magnitude. Therefore, a strong dependence of the initial mass function (IMF) on $\alpha _{\mathrm{vir}}$ may challenge the notion of a universal IMF. To determine the role of $\alpha _{\mathrm{vir}}$ on the IMF, we compare the star-particle mass functions obtained in high-resolution magnetohydrodynamical simulations including jet and heating feedback, with $\alpha _{\mathrm{vir}}=0.0625$, 0.125, and 0.5. We find that varying $\alpha _{\mathrm{vir}}$ from $\alpha _{\mathrm{vir}}\sim 0.5$ to $\alpha _{\mathrm{vir}}< 0.1$ shifts the peak of the IMF to lower masses by a factor of $\sim 2$ and increases the star formation rate by a similar factor. The dependence of the IMF and star formation rate on $\alpha _{\mathrm{vir}}$ is non-linear, with the dependence subsiding at $\alpha _{\mathrm{vir}}< 0.1$. Our study shows a systematic dependence of the IMF on $\alpha _{\mathrm{vir}}$. Yet, it may not be measurable easily in observations, considering the uncertainties, and the relatively weak dependence found in this study.

Mathew, Sajay Sunny (ORCID:0000000283818195)↗

DOE OSTI2026

Collapse of magnetized white dwarfs as site of heavy-element formation and kilonova signal

We present the first end-to-end calculation connecting the accretion-induced collapse (AIC) of a magnetized, rapidly rotating white dwarf to observable kilonova signatures, combining two-dimensional (2D) general-relativistic neutrino-magnetohydrodynamic simulations, followed by radiation hydrodynamics with in-situ nuclear network and 2D Monte Carlo radiative transfer with spatially resolved heating rates. Unlike all previous unmagnetized AIC models – which predicted proton-rich, $^{56}$Ni-dominated ejecta – strong magnetic fields eject ${\approx }\, 0.2\, \mathrm{ M}_\odot$ of neutron-rich material ($\langle Y_e \rangle \sim 0.24$) on dynamical time-scales, before neutrino irradiation can raise the electron fraction, enabling strong r-process nucleosynthesis up to and beyond the third peak. The resulting kilonova is lanthanide-rich ($X_{\rm lan} \approx 8~{{\ \rm per\ cent}}$) and dominated by near-infrared emission. We compute synthetic light curves in the Large Synoptic Survey Telescope and J ames Webb Space Telescope bands and find striking agreement, without parameter tuning, between the observations of AT 2023vfi/GRB 230307A and our broadband light curves for polar viewing angles. These results establish magnetized AIC as a viable channel for heavy r-process element production and a compelling progenitor candidate for long-duration gamma-ray bursts with kilonova signatures.

MHD↗

DOE OSTI2026

Probing the axion-electron coupling with NuSTAR observations of galaxies

We search for the existence of ultralight axions coupling to electrons and photons using data from the NuSTAR telescope directed toward the galaxies M82, M87, and M31. We focus on electron bremsstrahlung and Compton scattering for axion production in stars, summing over the stellar populations found in the target galaxies when computing the axion luminosity. We then compute the hard x-ray signal that arises from the conversion of these axions to photons in each galaxy’s magnetic fields, inferred from analog galaxies in cosmological magnetohydrodynamic simulations. Analyzing NuSTAR data toward these galaxies between roughly 20 to 70 keV, we find no evidence for axions and set leading constraints on the combined axion-electron and axion-photon coupling at the level of |𝑔 𝑎⁢𝑒⁢𝑒 ×𝑔 𝑎⁢𝛾⁢𝛾 | ≲ 8.3 ×10 −27 GeV −1 for 𝑚 𝑎 ≲ 10 −10 eV at 95% confidence, with M82 providing the most stringent constraints.

Axion-like particles↗

DOE OSTI2026

Gamma-ray signatures of 𝑟-process radioactivity from the collapse of magnetized white dwarfs

We predict the gamma-ray line emission from 𝑟-process nuclei synthesized in the ejecta of the accretion-induced collapse (AIC) of a magnetized, rapidly rotating white dwarf. Using ejecta from a two-dimensional general-relativistic neutrino-magnetohydrodynamic simulation, further evolved with a radiation-hydrodynamics code coupled to an in situ nuclear reaction network, we construct angle-dependent gamma-ray spectra in the 0.01–10 MeV band via composition-dependent ray tracing through the ejecta. The emission between ∼1 and 10 d is dominated by 132 I (𝑡 1/2 = 2.3 h), continuously replenished by the decay of its parent 132 Te (𝑡 1/2 = 3.2 d), with additional contributions from 131 I, 133 Xe, and 132 Te . At 𝑡 ≳ 20 d, 56 Co (from 56 Ni decay) becomes the primary emitter. The simultaneous presence of 𝑟 process and iron-peak gamma-ray lines is distinctive of AIC ejecta and absent in binary neutron star mergers, where iron-peak nuclei are generally not synthesized. Comparing with the 3⁢𝜎 continuum sensitivities of planned MeV gamma-ray telescopes (COSI, AMEGO-X, e-ASTROGAM, GRAMS, GammaTPC), we find the brightest 𝑟-process lines detectable to ∼10 Mpc by GammaTPC and GRAMS, with the signal approaching their sensitivity threshold at 30 Mpc. As a result, the 𝑟-process spectral features survive time integration over ∼30 d exposures, demonstrating robustness against the long observation times required by gamma-ray detectors.

Nuclear reactions↗

DOE OSTI2024

Electrical conductivity of a warm neutron star crust in magnetic fields: Neutron-drip regime

We compute the anisotropic electrical conductivity tensor of the inner crust of a compact star at nonzero temperature by extending a previous work on the conductivity of the outer crust. The physical scenarios, where such crust is formed, involve protoneutron stars born in supernova explosions, binary neutron star mergers, and accreting neutron stars. The temperature-density range studied covers the transition from a semidegenerate to a highly degenerate electron gas and assumes that the nuclei form a liquid, i.e., the temperature is above the melting temperature of the lattice of nuclei. The electronic transition probabilities include (i) the screening of electron-ion interaction in the hard-thermal-loop approximation for the QED plasma, (ii) the correlations of the ionic component in a one-component plasma, and (iii) finite nuclear size effects. The conductivity tensor is obtained from the Boltzmann kinetic equation in relaxation time approximation accounting for the anisotropy introduced by a magnetic field. The sensitivity of the results towards the matter composition of the inner crust is explored by using several compositions of the inner crust, which were obtained using different nuclear interactions and methods of solving the many-body problem. The standard deviations of relaxation time and components of the conductivity tensor from the average are below ≤25% except close to crust-core transition, where nonspherical nuclear structures are expected. Finally, our results can be used in dissipative magnetohydrodynamics simulations of warm compact stars.

Physics↗

DOE OSTI2025

Influence of neutrino-electron scattering and neutrino-pair annihilation on hypermassive neutron star

Here, we investigate the influence of inelastic neutrino microphysics in general-relativistic magnetohydrodynamics simulations of a hypermassive neutron star. In particular, we include species/energy groups coupled neutrino-matter interactions, such as inelastic neutrino-electron scattering and electron-positron annihilation kernels, into simulations up to 50 ms. Neutrino-electron inelastic scattering is known to have effective neutrino-matter energy exchange. We show that, with neutrino-electron inelastic scattering, simulations predict 75% higher disc mass with slightly different mass-averaged compositions, and 18% more ejected mass with similar distributions. The enhancement of the mass of the disc and the ejecta results in stronger baryon pollution, leading to less favorable jet launching environments. Furthermore, neutrino luminosities are about 50%, 40%, and 30% higher for electron neutrino, electron antineutrino, and heavy-lepton neutrinos. In contrast, we do not see any significant impacts due to electron-positron annihilation.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2020

Magnetic Signatures of Radiation-Driven Double Ablation Fronts

In experiments performed with the OMEGA EP laser system, magnetic field generation in double ablation fronts was observed. Proton radiography measured the strength, spatial profile, and temporal dynamics of self-generated magnetic fields as the target material was varied between plastic (CH), aluminum, copper, and gold. Two distinct regions of magnetic field are generated in mid-Z targets -- one produced by gradients from electron thermal transport and the second from radiation-driven gradients. Extended magnetohydrodynamic simulations including radiation transport reproduced key aspects of the experiment, including field generation and double ablation front formation.