Search NASA⌕ Search

SEARCH · Search NASA

Results for “collisionality”

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

Pressure–Strain Interaction as the Energy Dissipation Estimate in Collisionless Plasma

The dissipative mechanism in weakly collisional plasma is a topic that pervades decades of studies without a consensus solution. We compare several energy dissipation estimates based on energy transfer processes in plasma turbulence and provide justification for the pressure–strain interaction as a direct estimate of the energy dissipation rate. The global and scale-by-scale energy balances are examined in 2.5D and 3D kinetic simulations. We show that the global internal energy increase and the temperature enhancement of each species are directly tracked by the pressure–strain interaction. The incompressive part of the pressure–strain interaction dominates over its compressive part in all simulations considered. The scale-by-scale energy balance is quantified by scale filtered Vlasov–Maxwell equations, a kinetic plasma approach, and the lag dependent von Kármán–Howarth equation, an approach based on fluid models. We find that the energy balance is exactly satisfied across all scales, but the lack of a well-defined inertial range influences the distribution of the energy budget among different terms in the inertial range. Therefore, the widespread use of the Yaglom relation in estimating the dissipation rate is questionable in some cases, especially when the scale separation in the system is not clearly defined. In contrast, the pressure–strain interaction balances exactly the dissipation rate at kinetic scales regardless of the scale separation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Resolution Study of Thermonuclear Initiation in White Dwarf Tidal Disruption Events

Abstract We study the initiation of thermonuclear detonations in tidally disrupted white dwarf stars by intermediate-mass (10 3 M ⊙ ) black holes. The length scales required to resolve the initiation mechanism are not easily reached in 3D, so instead we have devised 2D proxy models, which, together with a logarithmic gridding strategy, can adequately capture detonation wave fronts as material undergoes simultaneous compression and stretching from tidal forces. We consider 0.15 and 0.6 M ⊙ white dwarf stars parameterized by tidal strengths in the range β = 4–23. High spatial resolution elucidates the manner and conditions leading to thermonuclear detonation, linking the initiation sequence to stellar composition and tidal strength. All of our models suffer sustained detonations triggered by a combination of adiabatic compression, mild thermonuclear preconditioning, and collisional heating, in degrees depending primarily on tidal strength. We find that many diagnostics, such as temperature, total released energy, and iron-group products, are fairly well converged (better than 10%) at resolutions below 10 km along the scale height of the orbital plane. The exceptions are intermediate-mass transients like calcium, which remain uncertain up to factors of 2, even at 1 km resolution.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Quantifying the Agyrotropy of Proton and Electron Heating in Turbulent Plasmas

An important aspect of energy dissipation in weakly collisional plasmas is that of energy partitioning between different species (e.g., protons and electrons) and between different energy channels. Here we analyse pressure–strain interaction to quantify the fractions of isotropic compressive, gyrotropic, and nongyrotropic heating for each species. An analysis of kinetic turbulence simulations is compared and contrasted with corresponding observational results from Magnetospheric Multiscale Mission data in the magnetosheath. In assessing how protons and electrons respond to different ingredients of the pressure–strain interaction, we find that compressive heating is stronger than incompressive heating in the magnetosheath for both electrons and protons, while incompressive heating is stronger in kinetic plasma turbulence simulations. Concerning incompressive heating, the gyrotropic contribution for electrons is dominant over the nongyrotropic contribution, while for protons nongyrotropic heating is enhanced in both simulations and observations. Variations with plasma β are also discussed, and protons tend to gain more heating with increasing β.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

CLASSY. VI. The Density, Structure, and Size of Absorption-line Outflows in Starburst Galaxies

Galaxy formation and evolution are regulated by the feedback from galactic winds. Absorption lines provide the most widely available probe of winds. However, since most data only provide information integrated along the line of sight, they do not directly constrain the radial structure of the outflows. In this paper, we present a method to directly measure the gas electron density in outflows (n e ), which in turn yields estimates of outflow cloud properties (e.g., density, volume filling factor, and sizes/masses). We also estimate the distance (r n ) from the starburst at which the observed densities are found. We focus on 22 local star-forming galaxies primarily from the COS Legacy Archive Spectroscopic SurveY (CLASSY). In half of them, we detect absorption lines from fine-structure excited transitions of Si $\tiny{II}$ (i.e., Si $\tiny{II}$ * ). We determine n e from relative column densities of Si $\tiny{II}$ and Si $\tiny{II}$ * , given Si $\tiny{II}$ * originates from collisional excitation by free electrons. We find that the derived n e correlates well with the galaxy's star formation rate per unit area. From photoionization models or assuming the outflow is in pressure equilibrium with the wind fluid, we get r n ~ 1–2r * or ~5r * , respectively, where r * is the starburst radius. Based on comparisons to theoretical models of multiphase outflows, nearly all of the outflows have cloud sizes large enough for the clouds to survive their interaction with the hot wind fluid. Most of these measurements are the first ever for galactic winds detected in absorption lines and, thus, will provide important constraints for future models of galactic winds.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Constraints on the Gamma-Ray Emission from Small Solar System Bodies with the Fermi Large Area Telescope Data

All known small solar system bodies have diameters between a few meters and a few thousands of kilometers. Based on the collisional evolution of solar system bodies, a larger number of asteroids with diameters down to ~2 m is thought to exist. As all solar system bodies, small bodies can be passive sources of high-energy gamma-rays, produced by the interaction of energetic cosmic rays impinging on their surfaces. Since the majority of known asteroids are in orbits between Mars and Jupiter (in a region known as the Main Belt), we expect them to produce a diffuse emission close to the ecliptic plane. In this work, we have studied the gamma-ray emission coming from the ecliptic using the data collected by the Large Area Telescope (LAT) onboard the Fermi satellite. We have fit the results with simulations of the gamma-ray intensity at the source level (calculated with the software FLUKA) to constrain the small solar system bodies population. Finally, we have proposed a model describing the distribution of asteroid sizes and we have used the LAT data to constrain the gamma-ray emission expected from this model and, in turn, on the model itself.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Expanding-box Quasilinear Model of the Solar Wind

The expanding-box model of the solar wind has been adopted in the literature within the context of magnetohydrodynamics, hybrid, and full particle-in-cell simulations to investigate the dynamic evolution of the solar wind. The present paper extends such a method to the framework of self-consistent quasilinear kinetic theory. It is shown that the expanding-box quasilinear methodology is largely equivalent to the inhomogeneous steady-state quasilinear model discussed earlier in the literature, but a distinction regarding the description of wave dynamics between the two approaches is also found. The expanding-box quasilinear formalism is further extended to include the effects of a spiraling solar-wind magnetic field as well as collisional age effects. The present finding shows that the expanding-box quasilinear approach and the steady-state global-kinetic models may be employed interchangeably in order to address other more complex problems associated with the solar-wind dynamics.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Ab-initio General-relativistic Neutrino-radiation Hydrodynamics Simulations of Long-lived Neutron Star Merger Remnants to Neutrino Cooling Timescales

Abstract We perform the first 3D ab-initio general-relativistic neutrino-radiation hydrodynamics of a long-lived neutron star merger remnant spanning a fraction of its cooling timescale. We find that neutrino cooling becomes the dominant energy loss mechanism after the gravitational-wave dominated phase (∼20 ms postmerger). Electron flavor antineutrino luminosity dominates over electron flavor neutrino luminosity at early times, resulting in a secular increase of the electron fraction in the outer layers of the remnant. However, the two luminosities become comparable ∼20–40 ms postmerger. A dense gas of electron antineutrinos is formed in the outer core of the remnant at densities ∼10 14.5 g cm −3 , corresponding to temperature hot spots. The neutrinos account for ∼10% of the lepton number in this region. Despite the negative radial temperature gradient, the radial entropy gradient remains positive, and the remnant is stably stratified according to the Ledoux criterion for convection. A massive accretion disk is formed from the material squeezed out of the collisional interface between the stars. The disk carries a large fraction of the angular momentum of the system, allowing the remnant massive neutron star to settle to a quasi-steady equilibrium within the region of possible, stable, rigidly rotating configurations. The remnant is differentially rotating, but it is stable against the magnetorotational instability. Other MHD mechanisms operating on longer timescales are likely responsible for the removal of the differential rotation. Our results indicate the remnant massive neutron star is thus qualitatively different from a protoneutron stars formed in core-collapse supernovae.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Magnetogenesis in a Collisionless Plasma: From Weibel Instability to Turbulent Dynamo

We report on a first-principles numerical and theoretical study of plasma dynamo in a fully kinetic framework. By applying an external mechanical force to an initially unmagnetized plasma, we develop a self-consistent treatment of the generation of “seed” magnetic fields, the formation of turbulence, and the inductive amplification of fields by the fluctuation dynamo. Driven large-scale motions in an unmagnetized, weakly collisional plasma are subject to strong phase mixing, which leads to the development of thermal pressure anisotropy. This anisotropy triggers the Weibel instability, which produces filamentary “seed” magnetic fields on plasma-kinetic scales. The plasma is thereby magnetized, enabling efficient stretching and folding of the fields by the plasma motions and the development of Larmor-scale kinetic instabilities such as the firehose and mirror. The scattering of particles off the associated microscale magnetic fluctuations provides an effective viscosity, regulating the field morphology and turbulence. During this process, the seed field is further amplified by the fluctuation dynamo until energy equipartition with the turbulent flow is reached. By demonstrating that equipartition magnetic fields can be generated from an initially unmagnetized plasma through large-scale turbulent flows, this work has important implications for the origin and amplification of magnetic fields in the intracluster and intergalactic mediums.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Collective Neutrino Oscillations and Heavy-element Nucleosynthesis in Supernovae: Exploring Potential Effects of Many-body Neutrino Correlations

In high-energy astrophysical processes involving compact objects, such as core-collapse supernovae or binary neutron star mergers, neutrinos play an important role in the synthesis of nuclides. Neutrinos in these environments can experience collective flavor oscillations driven by neutrino–neutrino interactions, including coherent forward scattering and incoherent (collisional) effects. Recently, there has been interest in exploring potential novel behaviors in collective oscillations of neutrinos by going beyond the one-particle effective or "mean-field" treatments. Here, we seek to explore implications of collective neutrino oscillations, in the mean-field treatment and beyond, for the nucleosynthesis yields in supernova environments with different astrophysical conditions and neutrino inputs. We find that collective oscillations can impact the operation of the νp-process and r-process nucleosynthesis in supernovae. The potential impact is particularly strong in high-entropy, proton-rich conditions, where we find that neutrino interactions can nudge an initial νp-process neutron-rich, resulting in a unique combination of proton-rich low-mass nuclei as well as neutron-rich high-mass nuclei. We describe this neutrino-induced neutron-capture process as the "νi-process." In addition, nontrivial quantum correlations among neutrinos, if present significantly, could lead to different nuclide yields compared to the corresponding mean-field oscillation treatments, by virtue of modifying the evolution of the relevant one-body neutrino observables.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Modeling the Far-infrared Polarization Spectrum of a High-mass Star-forming Cloud

The polarization spectrum, or wavelength dependence of the polarization fraction, of interstellar dust emission provides important insights into the grain alignment mechanism of interstellar dust grains. We investigate the far-infrared polarization spectrum of a realistic simulated high-mass star-forming cloud under various models of grain alignment and emission. We find that neither a homogeneous grain alignment model nor a grain alignment model that includes collisional dealignment is able to produce the falling spectrum seen in observations. On the other hand, we find that a grain alignment model with grain alignment efficiency dependent on local temperature is capable of producing a falling spectrum that is in qualitative agreement with observations of OMC-1. For the model most in agreement with OMC-1, we find no correlation between the temperature and the slope of the polarization spectrum. However, we do find a positive correlation between the column density and the slope of the polarization spectrum. We suggest this latter correlation to be the result of wavelength-dependent polarization by absorption.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

An Exceptionally Powerful, Radiatively Driven Ultrafast Outflow in the Rapidly Accreting AGN RE J1034+396

We report the analysis of ∼1 Ms of XMM-Newton observations of the rapidly accreting active galactic nucleus RE J1034+396. The 0.3–9 keV EPIC-pn spectra are well described by a model consisting of steep continuum emission from the corona accompanied by relativistically blurred reflection from a highly ionized accretion disk. The source is known to exhibit strong excess soft X-ray emission, which we show is well represented by thermal disk photons Comptonized by a warm plasma spanning the inner accretion flow. Additionally, the EPIC-pn data provide compelling evidence (ΔC ∼ 60 for four additional parameters) for the presence of an ultrafast outflow (UFO) with a line-of-sight velocity v/c = 0.307$^{+0.001}_{-0.005}$, and an emission signature consistent with reflection of the corona from modestly ionized, outflowing gas. The simultaneous 0.5–2.5 keV RGS spectra show clear absorption lines. Modelling of these data confirms the presence of the UFO and constrains its equivalent hydrogen column density, log N H /(atom cm −2 ) = 21.7$^{+0.1}_{-0.2}$. The RGS data also reveal at least two warm absorber components with a modest outflow velocity (1680$^{+40}_{-50}$ km s -1 ). The measured properties and time evolution of the UFO in RE J1034+396 suggest that it is formed from collisionally ionized plasma, launched from the disk surface and accelerated by radiation pressure. The high terminal velocity and substantial absorbing column density imply that the outflow carries sufficient momentum and energy to transform its environment, being capable of driving out essentially all dust and gas it interacts with along the line of sight, even if the AGN were initially surrounded by a Compton-thick absorber.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Stability and convergence of nuclear detonations in white dwarf collisions

We investigate the numerical stability of thermonuclear detonations in 1D accelerated reactive shocks and 2D binary collisions of equal-mass magnetized and unmagnetized white dwarf stars. To achieve high resolution at initiation sites, we devised geometric gridding and mesh velocity strategies specially adapted to the unique requirements of head-on collisional geometries, scenarios in which one expects maximum production of iron-group products. We study the effects of grid resolution and the limiting of temperature, energy, and reactants for different stellar masses, separations, magnetic fields, inigenerationtial compositions, detonation mechanisms, and limiter parameters across a range of cell sizes from 1 to 100 km. Our results set bounds on the parameter space of limiter amplitudes for which both temperature- and energy-limiting procedures yield consistent and monotonically convergent solutions. Within these bounds, we find that grid resolutions of 5 km or better are necessary for uncertainties in total released energy and iron-group products to drop below 10%. Intermediate-mass products (e.g., calcium) exhibit similar convergence trends but with somewhat greater uncertainty. These conclusions apply equally to pure C/O white dwarfs, multispecies compositions (including helium shells), magnetized and unmagnetized cores, and either single or multiple detonation scenarios.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Laboratory Measurements of n >= 3 K-shell Transition Energies of Sulfur Ions from F-like S viii to Li-like S xiv

Inner-shell transitions are ubiquitous in nonequilibrium collisionally ionized plasmas, such as supernova remnants, and in photoionized plasmas, such as outflows from active galactic nuclei and X-ray binaries. Inner-shell X-ray emission can help determine key parameters of these systems, such as ionization time, τ, and ionization parameter, ξ. Despite their importance, only theoretical inner-shell transition energies are available for many ions. To provide experimental benchmarks, we have measured the dominant n → 1 K-shell transitions of sulfur ions where n ≥ 3 from Li-like S xiv to F-like S viii using LLNL’s SuperEBIT electron beam ion trap and the NASA/GSFC EBIT Calorimeter Spectrometer (ECS). We identify over 30 spectral features and measure their energies with uncertainties in the ∼0.1–1 eV range. We compare these results to Flexible Atomic Code (FAC) and multireference Møller–Plesset (MR-MP) calculations and find differences between theory and experiment of ∼1 eV for FAC and <0.5 eV for most MR-MP calculations. We also compare these results to two widely used atomic databases, AtomDB and CHIANTI, and find discrepancies as high as 7 eV. Furthermore, many transitions are missing from these databases despite being prominent in our data.

Atomic data benchmarking↗
DOE OSTI2025

Turbulent Heating in Collisionless Low-beta Plasmas: Imbalance, Landau Damping, and Electron–Ion Energy Partition

An understanding of how turbulent energy is partitioned between ions and electrons in weakly collisional plasmas is crucial for modeling many astrophysical systems. Using theory and simulations of a four-dimensional reduced model of low-beta gyrokinetics (the “Kinetic Reduced Electron Heating Model”), we investigate the dependence of collisionless heating processes on plasma beta and imbalance (normalized cross-helicity). These parameters are important because they control the helicity barrier, the formation of which divides the parameter space into two distinct regimes with remarkably different properties. In the first, at lower beta and/or imbalance, the absence of a helicity barrier allows the cascade of injected power to proceed to small (perpendicular) scales, but its slow cascade rate makes it susceptible to significant electron Landau damping, in some cases leading to a marked steepening of the magnetic spectra on scales above the ion Larmor radius. In the second, at higher beta and/or imbalance, the helicity barrier halts the cascade, confining electron Landau damping to scales above the steep “transition-range” spectral break, resulting in dominant ion heating. We formulate quantitative models of these processes that compare well to simulations in each regime, and combine them with results of previous studies to construct a simple formula for the electron–ion heating ratio as a function of beta and imbalance. This model predicts a “winner takes all” picture of low-beta plasma heating, where a small change in the fluctuations' properties at large scales (the imbalance) can cause a sudden switch between electron and ion heating.

Interplanetary turbulence↗
DOE OSTI2025

On the Stability Analysis of Astrophysical Cooling Functions

To model the temperature evolution of optically thin astrophysical environments at MHD scales, radiative and collisional cooling rates are typically either pretabulated or fit into a functional form and then input into MHD codes as a radiative loss function. Thermal balance requires estimates of the analogous heating rates, which are harder to calculate, and due to uncertainties in the underlying dissipative heating processes these rates are often simply parameterized. The resulting net cooling function defines an equilibrium curve that varies with density and temperature. Such cooling functions can make the gas prone to thermal instability (TI), which will cause departures from equilibrium. There has been no systematic study of thermally unstable parameter space for nonequilibrium states. Motivated by our recent finding that there is a related linear instability, catastrophic cooling instability, that can dominate over TI, here we carry out such a study. We show that Balbus instability criteria for TI can be used to define a critical cooling rate, Λ c , that permits a nonequilibrium analysis of cooling functions through the mapping of TI zones. We furthermore illustrate how thermal conduction modifies the shape of TI zones. Upon applying a Λ c -based stability analysis to coronal loop simulations, we find that loops undergoing periodic episodes of coronal rain formation are linearly unstable to catastrophic cooling instability, while TI is stabilized by thermal conduction.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Enigmatic Line Broadening during Solar Flares: Magnetic Field Broadening?

The origin of the extreme broadening observed in chromospheric metal lines during solar and stellar flares—particularly Mg II h and k and Ca II H and K—remains poorly understood. These lines often display Lorentzian-like wings whose widths exceed standard Stark broadening predictions by factors of ∼30, with no known collisional mechanism capable of producing such enhancements. We posit that magnetic fields are responsible for this additional broadening, due to the increase of magnetic activity during flares. A magnetic field distribution of the form P(B) ∝ B −3 reproduces the observed Mg II profile wings while leaving H I Balmer lines and optically thin transitions largely unaffected. To explain the broadening using magnetic fields, the high-B tail can extend up to 10 6 G with extremely low probabilities where the filling factors are ≲10 −6 . We propose that observations of flares using spectropolarimetry can verify whether the anomalous broadening is from magnetic structures in flare ribbons.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

A DG-IMEX Method for Two-moment Neutrino Transport: Nonlinear Solvers for Neutrino–Matter Coupling

Neutrino-matter interactions play an important role in core-collapse supernova (CCSN) explosions as they contribute to both lepton number and/or four-momentum exchange between neutrinos and matter, and thus act as the agent for neutrino-driven explosions. Due to the multiscale nature of neutrino transport in CCSN simulations, an implicit treatment of neutrino-matter interactions is desired, which requires solutions of coupled nonlinear systems in each step of the time integration scheme. In this paper we design and compare nonlinear iterative solvers for implicit systems with energy coupling neutrino-matter interactions commonly used in CCSN simulations. Specifically, we consider electron neutrinos and antineutrinos, which interact with static matter configurations through the Bruenn 85 opacity set. The implicit systems arise from the discretization of a nonrelativistic two-moment model for neutrino transport, which employs the discontinuous Galerkin (DG) method for phase-space discretization and an implicit-explicit (IMEX) time integration scheme. In the context of this DG-IMEX scheme, we propose two approaches to formulate the nonlinear systems — a coupled approach and a nested approach. For each approach, the resulting systems are solved with Anderson-accelerated fixedpoint iteration and Newton’s method. The performance of these four iterative solvers has been compared on relaxation problems with various degree of collisionality, as well as proto-neutron star deleptonization problems with several matter profiles adopted from spherically symmetric CCSN simulations. Here, numerical results suggest that the nested Anderson-accelerated fixed-point solver is more efficient than other tested solvers for solving implicit nonlinear systems with energy coupling neutrino-matter interactions.

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↗