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 721 records · Page 40

thornado+FLASH-X: A Hybrid Discontinuous Galerkin–Implicit-explicit and Finite-volume Framework for Neutrino-radiation Hydrodynamics in Core-collapse Supernovae

We present neutrino-transport algorithms implemented in the toolkit for high-order neutrino-radiation hydrodynamics (thornado) and their coupling to self-gravitating hydrodynamics within the adaptive mesh refinement–based multiphysics simulation framework FLASH-X. thornado, developed primarily for simulations of core-collapse supernovae (CCSNe), employs a spectral, six-species two-moment formulation with algebraic closure and special-relativistic observer corrections accurate to $\mathcal{O}(v/c)$, and uses discontinuous Galerkin (DG) methods for phase-space discretization combined with implicit-explicit time stepping. A key development is a nonlinear neutrino–matter coupling algorithm based on nested fixed-point iteration with Anderson acceleration, enabling fully implicit treatment of collisional processes, including energy-coupling interactions such as neutrino–electron scattering and pair production. Coupling to finite-volume (FV) hydrodynamics is achieved through a hybrid DG-FV representation of the fluid variables and operator-split evolution within FLASH-X. The implementation is verified using basic transport tests with idealized opacities and relaxation and deleptonization problems with tabulated microphysics. Spherically symmetric CCSN simulations demonstrate accuracy and robustness of the coupled scheme, including close agreement with the CCSN simulation code Chimera. An axisymmetric CCSN simulation further demonstrates the viability of DG-based neutrino transport for multidimensional supernova modeling within FLASH-X. thornado’s neutrino-transport solver is GPU-enabled using OpenMP offloading or OpenACC, and all CCSN applications included in this work use the GPU implementation. Together, these results establish a foundation for future enhancements in physics fidelity, numerical algorithms, and computational performance, for increasingly realistic large-scale CCSN simulations.

Endeve, Eirik [Oak Ridge National Laboratory (ORNL↗

On the Injection Scale of the Turbulence in the Partially Ionized Very Local Interstellar Medium

The cascade of magnetohydrodynamic (MHD) turbulence is subject to ion–neutral collisional damping and neutral viscous damping in the partially ionized interstellar medium. By examining the damping effects in the warm and partially ionized local interstellar medium, we find that the interstellar turbulence is damped by neutral viscosity at ~261 au and cannot account for the turbulent magnetic fluctuations detected by Voyager 1 and 2. The MHD turbulence measured by Voyager in the very local interstellar medium (VLISM) should be locally injected in the regime where ions are decoupled from neutrals for its cascade to survive the damping effects. With the imposed ion–neutral decoupling condition and the strong turbulence condition for the observed Kolmogorov magnetic energy spectrum, we find that the turbulence in the VLISM is sub-Alfvénic, and its largest possible injection scale is ~194 au.

79 ASTRONOMY AND ASTROPHYSICS↗

XRISM Finds the Changing-look Active Galactic Nucleus NGC 1365 in an Extended Low State: A Dense, Highly Ionized Outflow Obscures the Central Source

We present the first XRISM/Resolve observations of the active galactic nucleus NGC 1365, obtained in 2024 February and July. NGC 1365 is known for rapid transitions between Compton-thick and Compton-thin states, along with strong absorption from a highly ionized wind. During our observations, the source was found in a persistent low-flux state, characterized by a decrease in hard-X-ray luminosity and significant line-of-sight obscuration. In this state, XRISM/Resolve reveals clear Fe XXV and Fe XXVI absorption lines together with, for the first time in this source, corresponding emission lines. These features may arise either from reemission from a photoionized wind (P Cygni profile) or from collisionally ionized gas associated with outflow-driven shocks in the interstellar medium. We estimate the wind launch radius to be approximately 10 16 cm (∼10 4 R g ), consistent with the location of the X-ray broadline region. We also resolve a broadened Fe Kα line by σ ∼ 1300 km s −1 , placing it at similar scales to the wind, consistent with radii inferred from disk-broadening models and the variability of the Fe Kα broad line. The similarity of the Fe Kα profile to the Hβ wing and broad Paα width indicates that the X-ray-emitting region is likely cospatial with the optical/infrared broadline region and originates from the same gas.

Astronomy and AstroPhysics↗

Impacts of secondary ice production on Arctic mixed-phase clouds based on ARM observations and CAM6 single-column model simulations

For decades, measured ice crystal number concentrations have been found to be orders of magnitude higher than measured ice-nucleating particle number concentrations in moderately cold clouds. This observed discrepancy reveals the existence of secondary ice production (SIP) in addition to the primary ice nucleation. However, the importance of SIP relative to primary ice nucleation remains highly unclear. Furthermore, most weather and climate models do not represent SIP processes well, leading to large biases in simulated cloud properties. This study demonstrates a first attempt to represent different SIP mechanisms (frozen raindrop shattering, ice–ice collisional breakup, and rime splintering) in a global climate model (GCM). The model is run in the single column mode to facilitate comparisons with the Department of Energy (DOE)'s Atmospheric Radiation Measurement (ARM) Mixed-Phase Arctic Cloud Experiment (M-PACE) observations. We show the important role of SIP in four types of clouds during M-PACE (i.e., multilayer, single-layer stratus, transition, and frontal clouds), with the maximum enhancement in ice crystal number concentrations up to 4 orders of magnitude in moderately supercooled clouds. We reveal that SIP is the dominant source of ice crystals near the cloud base for the long-lived Arctic single-layer mixed-phase clouds. The model with SIP improves the occurrence and phase partitioning of the mixed-phase clouds, reverses the vertical distribution pattern of ice number concentrations, and provides a better agreement with observations. The findings of this study highlight the importance of considering SIP in GCMs.

54 ENVIRONMENTAL SCIENCES↗

Primary and secondary ice production: interactions and their relative importance

A discrepancy of up to 5 orders of magnitude between ice crystal and ice nucleating particle (INP) number concentrations was found in the measurements, indicating the potentially important role of secondary ice production (SIP) in the clouds. However, the interactions between primary and SIP processes and their relative importance remain unexplored. In this study, we implemented five different ice nucleation schemes as well as physical representations of SIP processes (i.e., droplet shattering during rain freezing, ice-ice collisional break-up, and rime splintering) in the Community Earth System Model version 2 (CESM2). We ran CESM2 in the single column mode for model comparisons with the DOE Atmospheric Radiation Measurement (ARM) Mixed-Phase Arctic Cloud Experiment (M-PACE) observations. We found that the model experiments with aerosol-aware ice nucleation schemes and SIP processes yield the best simulation results for the M-PACE single-layer mixed-phase clouds. We further investigated the relative importance of ice nucleation and SIP to ice number and cloud phase as well as interactions between ice nucleation and SIP in the M-PACE single-layer mixed-phase clouds. Our results show that SIP contributes 80 % to the total ice formation and transforms ~30 % of pure liquid-phase clouds simulated in the model experiments without considering SIP into mixed-phase clouds. The SIP is not only a result of ice crystals produced from ice nucleation, but also competes with the ice nucleation by reducing the number concentrations of cloud droplets and cloud-borne dust INPs. Conversely, strong ice nucleation also suppresses SIP by glaciating mixed-phase clouds and thereby reducing the amount of precipitation particles (rain and graupel).

54 ENVIRONMENTAL SCIENCES↗

Dynamic Phase Alignment in Inertial Alfven Turbulence

In weakly-collisional plasma environments with sufficiently low electron beta, Alfvenic turbulence transforms into inertial Alfvenic turbulence at scales below the electron skin-depth. We argue that, in inertial Alfvenic turbulence, both energy and generalized kinetic helicity exhibit direct cascades. We demonstrate that the two cascades are compatible due to the existence of a strong scale-dependence of the phase alignment angle between velocity and magnetic field fluctuations, with the phase alignment angle scaling as proportional to the perpendicular wavenumber to the power of -1. The kinetic and magnetic energy spectra scale as the perpendicular wavenumber to the power of -5/3 and -11/3, respectively. As a result of the dual direct cascade, the generalized-helicity spectrum scales as the perpendicular wavenumber to the power of -5/3, implying progressive balancing of the turbulence as the cascade proceeds to smaller scales in the sub-electron-skin-depth range. Turbulent eddies exhibit a phase-space anisotropy, with the parallel wavenumber scaling as the perpendicular wavenumber to the power of 5/3, consistent with critically-balanced inertial Alfven fluctuations. Our results may be applicable to a variety of geophysical, space, and astrophysical environments, including the Earth's magnetosheath and ionosphere, solar corona, non-relativistic pair plasmas, as well as to strongly rotating non-ionized fluids.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Multi-device Studies of Pedestal Physics and Confinement in the I-mode Regime

This paper describes joint ITPA studies of the I-mode regime, which features an edge thermal barrier but with L-mode-like particle and impurity transport, and without ELMs. The regime has been demonstrated on Alcator C-Mod, ASDEX Upgrade and DIII-D tokamaks, over a wide range of device parameters and pedestal conditions, including low collisionality. Access to I-mode is independent of heating method (neutral beam injection, ICRH and/or ECRH). Normalized energy confinement H98y,2 ≥ 1 has been achieved for a range of 3≤q95≤4.9 and scales favourably with power. Changes in turbulence in the pedestal region accompany the transition from L-mode to I-mode. The L-I threshold increases with plasma density and current, and with device size, but weakly with BT. The upper limit of power for I-modes increases with toroidal field, and is largest on Alcator C-Mod at B > 5 T. Issues for extrapolation to ITER and other future fusion devices are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Alcator C-Mod: Research in support of ITER and steps beyond

This paper presents an overview of recent highlights from research on Alcator C-Mod. Significant progress has been made across all research areas over the last two years, with particular emphasis on divertor physics and power handling, plasma-material-interaction studies, ELM-suppressed pedestal dynamics, core transport and turbulence, and RF heating and current drive utilizing Ion Cyclotron and Lower Hybrid tools. Specific results of particular relevance to ITER include: inner wall SOL transport studies that have led, together with results from other experiments, to the change of the detailed shape of the inner wall in ITER; runaway electron studies showing that the critical electric field required for runaway generation is much higher than predicted from collisional theory; core tungsten impurity transport studies reveal that tungsten accumulation is naturally avoided in typical C-Mod conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Multi-scale transport in the DIII-D ITER baseline scenario with direct electron heating and projection to ITER

Multi-scale fluctuations measured by turbulence diagnostics spanning long and short wavelength spatial scales impact energy confinement and the scale-lengths of plasma kinetic profiles in the DIII-D ITER baseline scenario with direct electron heating. Contrasting discharge phases with ECH + neutral beam injection (NBI) and NBI only at similar rotation reveal higher energy confinement and lower fluctuations when only NBI heating is used. Modeling of the core transport with TGYRO using the TGLF turbulent transport model and NEO neoclassical transport reproduces the experimental profile changes upon application of direct electron heating and indicates that multi-scale transport mechanisms are responsible for changes in the temperature and density profiles. Intermediate and high-k fluctuations appear responsible for the enhanced electron thermal flux, and intermediate-k electron modes produce an inward particle pinch that increases the inverse density scale length. Projection to ITER is performed with TGLF and indicates a density profile that has a finite scale length due to intermediate-k electron modes at low collisionality and increases the fusion gain. For a range of E×B shear, the dominant mechanism that increases fusion performance is suppression of outward low-k particle flux and increased density peaking.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Investigations of LHW-plasma coupling and current drive at high density related to H-mode experiments in EAST

Two important issues in achieving lower hybrid current drive (LHCD) high confinement plasma in EAST are to improve lower hybrid wave (LHW)-plasma coupling and to drive the plasma current at a high density. Studies in different configurations with different directions of toroidal magnetic field (Bt) show that the density near the antenna is affected by both the radial electric field induced by plasma without a LHW (Er-plasma) in the scrape off layer (SOL), and the radial electric field induced by LHW power (Er-LH) near the grill. Investigations indicate that Er× Btin the SOL leads to a different effect of configuration on the LHW-plasma coupling and Er-LH× Btaccounts for the asymmetric density behaviour in the SOL observed in the experiments, where Eris the total radial electric field in the SOL. Modelling of parametric instability (PI), collisional absorption (CA) and scattering from density fluctuations (SDF) in the edge region, performed considering the parameters of high density LHCD experiments in EAST, has shown that these mechanisms could be responsible for the low current drive (CD) efficiency at high density. Radiofrequency probe spectra, useful for documenting PI occurrence, show sidebands whose amplitude in the case of the lithiated vacuum chamber is smaller than in the case of poor lithiation, consistently with growth rates from PI modeling of the respective reference discharges. Since strong lithiation is also expected to diminish the parasitic effect on the LHCD of the remaining possible mechanisms, this appears to be a useful method for improving LHCD efficiency at a high density.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

One and Two Dimensional Quantum Lattice Algorithms for Maxwell Equations in Inhomogeneous Scalar Dielectric Media I: Theory

A quantum lattice algorithm (QLA) is developed for Maxwell equations in scalar dielectric media using the Riemann-Silberstein representation on a Cartesian grid. For x-dependent and y-dependent dielectric inhomogeneities, the corresponding QLA requires a minimum of 8 qubits/spatial lattice site. This is because the corresponding Pauli spin matrices have off-diagonal components which permit the local collisional entanglement of these qubits. However, z-dependent inhomogeneities require a QLA with a minimum of 16 qubits/lattice site since the Pauli spin matrix σz is diagonal. For 2 dimensional inhomogeneities, one can readily couple the 8-8 qubit schemes for x-y variations. z-x and y-z variations can be treated by either a 16-8 qubit scheme or a 16-16 qubit representation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Investigating profile stiffness and critical gradients in shaped TCV discharges using local gyrokinetic simulations of turbulent transport

The experimental observation made on the TCV tokamak of a significant confinement improvement in plasmas with negative triangularity ( δ < 0) compared to those with standard positive triangularity has been interpreted in terms of different degrees of profile stiffness (Sauter et al. Phys. Plasmas 21 055906, 2014) and/or different critical gradients. Employing the Eulerian gyrokinetic code GENE (Jenko et al. Phys. Plasmas 7 1904, 2000), profile stiffness and critical gradients are studied under TCV relevant conditions. For the considered experimental discharges, trapped electron modes (TEMs) and electron temperature gradient (ETG) modes are the dominant microinstabilities, with the latter providing a significant contribution to the non-linear electron heat fluxes near the plasma edge. Two series of simulations with different levels of realism are performed, addressing the question of profile stiffness at various radial locations. Retaining finite collisionality, impurities and electromagnetic effects, as well as the physical electron-to-ion mass ratio are all necessary in order to approach the experimental flux measurements. However, flux-tube simulations are unable to fully reproduce the TCV results, pointing towards the need to carry out radially nonlocal (global) simulations, i.e. retaining finite machine size effects, in a future study. Some conclusions about the effect of triangularity can nevertheless be drawn based on the flux-tube results. In particular, the importance of considering the sensitivity to both temperature and density gradient is shown. The flux tube results show an increase of the critical gradients towards the edge, further enhanced when δ < 0, and they also appear to indicate a reduction of profile stiffness towards plasma edge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Deep modelling of plasma and neutral fluctuations from gas puff turbulence imaging

The role of turbulence in setting boundary plasma conditions is presently a key uncertainty in projecting to fusion energy reactors. To robustly diagnose edge turbulence, we develop and demonstrate a technique to translate brightness measurements of HeI line radiation into local plasma fluctuations via a novel integrated deep learning framework that combines neutral transport physics and collisional radiative theory for the $3^3 D - 2^3 P$ transition in atomic helium. The tenets for experimental validity are reviewed, illustrating that this turbulence analysis for ionized gases is transferable to both magnetized and unmagnetized environments with arbitrary geometries. Based upon fast camera data on the Alcator C-Mod tokamak, we present the first 2-dimensional time-dependent experimental measurements of the turbulent electron density, electron temperature, and neutral density revealing shadowing effects in a fusion plasma using a single spectral line.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Final Technical Report for the SciDAC Center for the Study of Plasma Microturbulence

This final technical report for the SciDAC Center for the Study of Plasma Microturbulence documents contributions from MIT to the project over the period of the award (February 15, 2011 to February 14, 2020). This project focused on the development of neoclassical and gyrokinetic simulation capabilities relevant for the description of particle, momentum and energy transport in tokamak fusion plasmas. Extensive validation against experimental results from Alcator C-Mod and DIII-D experiments was performed. Overall areas of emphasis included gyrokinetic simulations, experiments, and observation of density gradient driven Trapped Electron Mode (TEM) turbulence relevant to burning plasmas; theory and gyrokinetic simulations of the parametric variation of the nonlinear upshift of the TEM critical density gradient; identification of a new mechanism for the isotope effect due to density gradient driven TEM turbulence; new synthetic diagnostics based on full-wave simulations which enabled unprecedented agreement between gyrokinetic simulations and measured fluctuation spectra with and without strong electron heating; development of new continuum pedestal neoclassical and turbulence codes; formulation and implementation of the first gyrokinetic exact Fokker-Planck collision operator in gyrokinetic turbulence simulations; development of new spectral approaches in velocity space for highly efficient treatments of collisions; new model collision operators implemented in the discontinuous Galerkin gyrokinetic edge turbulence code GKEYLL; and a new model and spectral code to study the interaction of ion- and electron-scale turbulence, which will serve as an efficient test-bed for multi-rate and multi-scale methods. These fundamental advances contribute to the development of a validated, predictive capability for turbulent and collisional particle, momentum and energy transport in both the core and edge of magnetically confined fusion plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Nanostructure Accelerators: Novel concept and path to its realization

TeV/m acceleration gradients using crystals as originally envisioned by R. Hofstadter, an early pioneer of HEP, have remained unrealizable. Fundamental obstacles that have hampered efforts on particle acceleration using bulk-crystals arise from collisional energy loss and emittance degradation in addition to severe beam disruption despite the favorable effect of particle channeling along interatomic planes in bulk. We aspire for the union of nanoscience with accelerator science to not only overcome these problems using nanostructured tubes to avoid direct impact of the beam on bulk ion-lattice but also to utilize the highly tunable characteristics of nanomaterials. We pioneer a novel surface wave mechanism in nanostructured materials with a strong electrostatic component which not only attains tens of TeV/m gradients but also has focusing fields. Under our initiative, the proof-of-principle demonstration of tens of TeV/m gradients and beam nanomodulation is underway. Realizable nanostructure accelerators naturally promise new horizons in HEP as well as in a wide range of areas of research that utilize beams of high-energy particles or photons.

43 PARTICLE ACCELERATORS↗

The Impact of Anomalous Resistivity in Vacuum Contaminant Plasmas on the Electrothermal Instability

This manuscript presents an assessment of the electrothermal instability (ETI) in the presence of anomalous resistivity (AR) in vacuum contaminant plasmas (VCP) when applied to a magnetized liner inertial fusion (MagLIF)-like load.Pulsed-power driven dielectrically coated metallic liners, like in MagLIF, experience the current-driven electrothermal instability which occurs when a material’s resistivity changes with temperature and is subject to ohmic heating. Large scale pulsed-power facilities that use magnetically insulated transmission lines (MITL) have been shown to generate low-density plasma which enters the target chamber and coalesces around the load. The low-density high-temperature vacuum contaminant plasmas (VCP) can parasitically divert current from the load through causing a short in the anode-cathode gap inside the target chamber. Resistive magnetohydrodynamic (MHD) simulations of these VCP experience unphysical runaway ohmic heating due to under predicting the resistivity by using a purely collisional resistivity model.AR provides a physics-based way to address this runaway heating through increasing the resistivity in a proportional way with the drift speed. In this work, 1D simulations probe the effect that AR in VCP has on the magnetic diffusion rate, and 2D simulations show how this effect manifests in the nonlinear striation form of the ETI for a MagLIF-like load. Beryllium and aluminum dielectrically coated liners are used for the 1D and 2D simulations in this work. The1D simulations show that a VCP causes a delay in the current delivery to the load by upwards of 8 ns at 60 ns into atypical current pulse for an 1 × 10 7 A scale pulsed-power accelerator. The 2D simulations show the delay observed in the 1D simulations is reduced substantially (4 ns) by 2D turbulence that disrupts the VCP layer early in time. Here, the 2Dsimulated ETI growth varies across AR models, more so for the beryllium liner than the aluminum liner because the beryllium liner shows an enhanced rate of penetration for the magnetic diffusion wave in comparison to the aluminum liner. The 2D simulations show the bulk dielectric thickness varies across AR models with the Davidson AR model beingthe largest and the Buneman AR model being the smallest, and in connection with the thickness the Rayleigh-Taylor bubble-spike distances varies correspondingly.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Quantifying the Relevance of Prompt Dissociations of Resonance Stabilized Free Radicals in Flames

In prior studies from our laboratory we have shown that weakly-bound free radicals have a propensity to promptly dissociate at elevated temperatures relevant to flames. Dissociation can occur so rapidly at high temperatures in such weakly-bound radicals that it interferes with collisional relaxation of internal energies leading to prompt dissociation. While prompt dissociation is a ubiquitous feature in such free radical dissociations, we have hypothesized that its relevance to global combustion observables such as flame speeds depends on satisfying specific criteria. Two criteria that motivate the present study are that the free radicals must play a central role in the oxidation mechanism and that their thermal dissociation rate coefficients should be small enough. Resonance stabilized radicals have bond strengths ~ 50 kcal/mol (or larger) and consequently are known to dissociate slowly and persist even in high temperature environments. In this work we have calculated prompt dissociation probabilities for some simple allylic (allyl, methyl allyl, and dimethyl allyl) radicals. Flame simulations were performed to quantify the relevance of including prompt dissociations of these allylic radicals.

Allyl Radicals↗

Recycling Measurements for the Lithium Tokamak Experiment

The lithium tokamak experiment (LTX) is a spherical torus designed to investigate a low-recycling wall regime for magnetically confined plasmas. Traditional measurements of recycling using Dα emission are not feasible in LTX due to reflections off the plasma-facing-components: stainless-steel shell, liquid lithium layer, and molybdenum limiters. Since lithium has low VUV reflectivity, LTX will use Lyman-α emission to measure particle recycling. LTX will use 10- and 20-channel AXUV diode arrays with Lyman-α filters to view recycling off the inboard and outboard shells. A single detector views recycling off a limiter. Measurements will be used with neutral modeling to calculate recycling and neutral deposition. To investigate the performance of an AXUV diode with directly-deposited filter, measurements were obtained with this detector and a visible detector viewing a calibration plasma then compared with emission calculated using a collisional-radiative model.

Physics - Plasma physics↗