Search NASA⌕ Search

SEARCH · Search NASA

Results for “KINETIC EQUATION”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Polarized QED cascades

By taking the spin and polarization of the electrons, positrons and photons into account in the strong-field QED processes of nonlinear Compton emission and pair production, we find that the growth rate of QED cascades in ultra-intense laser fields can be substantially reduced. While this means that fewer particles are produced, we also found them to be highly polarized. We further find that the high-energy tail of the particle spectra is polarized opposite to that expected from Sokolov–Ternov theory, which cannot be explained by just taking into account spin-asymmetries in the pair production process, but results significantly from 'spin-straggling'. We employ a kinetic equation approach for the electron, positron and photon distributions, each of them spin/polarization-resolved, with the QED effects of photon emission and pair production modelled by a spin/polarization dependent Boltzmann-type collision operator. For photon-seeded cascades, depending on the photon polarization, we find an excess or a shortage of particle production in the early stages of cascade development, which provides a path towards a controlled experiment. Throughout this paper we focus on rotating electric field configuration, which represent an idealized model and allows for a straightforward interpretation of the observed effects.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Three-dimensional core-collapse supernova models with phenomenological treatment of neutrino flavor conversions

Abstract We perform three-dimensional supernova simulations with a phenomenological treatment of neutrino flavor conversions. We show that the explosion energy can increase to as high as $\sim 10^{51}$ erg depending on the critical density for the onset of flavor conversions, due to a significant enhancement of the mean energy of electron antineutrinos. Our results confirm previous studies showing such energetic explosions, but for the first time in three-dimensional configurations. In addition, we predict neutrino and gravitational wave (GW) signals from a nearby supernova explosion aided by flavor conversions. We find that the neutrino event number decreases because of the reduced flux of heavy-lepton neutrinos. In order to detect GWs, next-generation GW telescopes such as Cosmic Explorer and the Einstein Telescope are needed even if the supernova event is located at the Galactic Center. These findings show that the neutrino flavor conversions can significantly change supernova dynamics and highlight the importance of further studies on the quantum kinetic equations to determine the conditions of the conversions and their asymptotic states.

Mori, Kanji↗
DOE OSTI2020

Partial equilibration of anti-Pfaffian edge modes at ν = 5 / 2

The thermal Hall conductance K of the fractional quantum Hall state at filling fraction ν = 5/2 has recently been measured to be K = 2.5π 2 $k$$^{2}_{B}$T/3h [M. Banerjee et al. Nature (London) 559, 205 (2018)]. Here, the half-integer value of this result (in units of π 2 $k$$^{2}_{B}$T/3h) provides strong evidence for the presence of a Majorana edge mode and a corresponding quantum Hall state hosting quasiparticles with non-Abelian statistics. Whether this measurement points to the realization of the PH-Pfaffian or the anti-Pfaffian state has been the subject of debate. Here we consider the implications of this measurement for anti-Pfaffian edge-state transport. We show that in the limit of strong Coulomb interaction and an approximate spin degeneracy in the lowest Landau level, the anti-Pfaffian state admits low-temperature edge phases that are consistent with the Hall conductance measurements. Furthermore, these edge phases can exhibit fully equilibrated electrical transport coexisting with partially equilibrated heat transport over a range of temperatures. Through a study of the kinetic equations describing low-temperature electrical and heat transport of these edge states, we determine regimes of parameter space, controlling the interactions between the different edge modes that agree with experiment.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Zero-field spin resonance in graphene with proximity-induced spin-orbit coupling

In this work we investigate collective spin excitations in graphene with proximity-induced spin-orbit coupling (SOC) of the Rashba and valley-Zeeman types, as it is the case, e.g., for graphene on transition-metal-dichalcogenide substrates. It is shown that, even in the absence of an external magnetic field, such a system supports collective spin modes, which correspond to coupled oscillations of the uniform and valley-staggered magnetizations. These modes can be detected via both zero-field electron spin resonance (ESR) and zero-field electric-dipole spin resonance (EDSR), with EDSR response arising solely due to Rashba SOC. We analyze the effect of electron-electron interaction within the Fermi-liquid kinetic equation and show that the interaction splits both the ESR and EDSR peaks into two. The magnitude of splitting and the relative weights of the resonances can be used to extract the spin-orbit coupling constants and many-body interaction parameters that may not be accessible by other methods.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
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 OSTI2020

Relaxation times for Bose-Einstein condensation in axion miniclusters

We study the Bose condensation of scalar dark matter in the presence of both gravitational and self-interactions. Axions and other scalar dark matter in gravitationally bound miniclusters or dark matter halos are expected to condense into Bose-Einstein condensates called Bose stars. This process has been shown to occur through attractive self-interactions of the axionlike particles or through the field’s self gravitation. We show that in the high-occupancy regime of scalar dark matter, the Boltzmann collision integral does not describe either gravitational or self-interactions, and derive kinetic equations valid for these interactions. We use this formalism to compute relaxation times for the Bose-Einstein condensation, and find that condensation into Bose stars could occur within the lifetime of the Universe. The self-interactions reduce the condensation time only when they are very strong.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Once-in-a-lifetime encounter models for neutrino media: From coherent oscillations to flavor equilibration

Collective neutrino oscillations are typically studied using the lowest-order quantum kinetic equation, also known as the mean-field approximation. However, some recent quantum many-body simulations suggest that quantum entanglement among neutrinos may be important and may result in flavor equilibration of the neutrino gas. In this work, we develop new quantum models for neutrino gases in which any pair of neutrinos can interact at most once in their lifetimes. A key parameter of our models is γ = μ Δ z , where μ is the neutrino coupling strength, which is proportional to the neutrino density, and Δ z is the duration over which a pair of neutrinos can interact each time. Our models reduce to the mean-field approach in the limit γ → 0 and achieve flavor equilibration in time t ≫ ( γ μ ) − 1 . These models demonstrate the emergence of coherent flavor oscillations from the particle perspective and may help elucidate the role of quantum entanglement in collective neutrino oscillations. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Chaos in inhomogeneous neutrino fast flavor instability

In dense neutrino gases, the neutrino-neutrino coherent forward scattering gives rise to a complex flavor oscillation phenomenon not fully incorporated in simulations of neutron star mergers (NSM) and core collapse supernovae (CCSNe). Moreover, it has been proposed to be chaotic, potentially limiting our ability to predict neutrino flavor transformations in simulations. To address this issue, we explore how small flavor perturbations evolve in the nonlinear regime of the neutrino quantum kinetic equation within a narrow centimeter-scale region inside a NSM and a toy neutrino distribution. Our findings reveal that paths in the flavor state space of solutions with similar initial conditions diverge exponentially, exhibiting chaos. This inherent chaos makes the microscopic scales of neutrino flavor transformations unpredictable. However, the domain-averaged neutrino density matrix remains relatively stable, with chaos minimally affecting it. This particular property suggests that domain-averaged quantities remain reliable despite the exponential amplification of errors. Published by the American Physical Society 2024

Astronomy & Astrophysics↗
DOE OSTI2024

Three-flavor, full momentum space neutrino spin oscillations in neutron star mergers

In the presence of anisotropic neutrino and antineutrino fluxes, the quantum kinetic equations drive coherent oscillations in neutrino helicity, frequently referred to as spin oscillations. These oscillations depend directly on the absolute mass scale and Majorana phase, but are usually too transient to produce important effects. In this paper we present a full momentum-space analysis of Majorana neutrino spin oscillations in a snapshot of a three-dimensional neutron star merger simulation. We find an interesting angular dependence that allows for that resonant and adiabatic oscillations to occur along specific directions in a large volume of the merger remnant. The solid angle spanned by these directions is extremely narrow in general. We then analyze spin transformation in the presence of flavor transformation by characterizing how the effect's resonance and timescale change during a fast flavor instability. For this analysis, we derive a generalized resonance condition that poses a restrictive requirement for resonance to exist in any flavor channel. Finally, we determine that spin oscillations at all locations in the merger snapshot have a length scale that is too large for significant oscillations to be expected even where there exist resonant and adiabatic directions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Ionization waves in low-current dc discharges in noble gases obtained with a hybrid kinetic-fluid model

A hybrid kinetic-fluid model is used to study ionization waves (striations) in a low-current plasma column of dc discharges in noble gases. Coupled solutions of a kinetic equation for electrons, a drift-diffusion equation for ions, and a Poisson equation for the electric field are obtained to clarify the nature of plasma stratification in the positive column. A simplified two-level excitation-ionization model is used for the conditions when the nonlinear effects due to stepwise ionization, gas heating, and Coulomb interactions among electrons are negligible. It is confirmed that the nonlocal effects are responsible for the formation of moving striations in dc discharges at low plasma densities and low values of pR (the product of gas pressure and tube radius). Here, the calculated properties of self-excited waves of S–, P–, and R types in neon and S type in argon agree with available experimental data. The reason for helium plasma stability to stratification is clarified. It is shown that sustaining stratified plasma is more efficient than striation-free plasma when the ionization rate is a nonlinear function of the electric field. However, the nonlinear dependence of the ionization rate on the electric field is not required for plasma stratification. Striations of S–, P–, and R types in neon exist with minimal or no ionization enhancement. Effects of the column length and plasma density on the wave properties are demonstrated.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Does a Quark-Gluon Plasma Feature an Extended Hydrodynamic Regime?

Here we investigate the response of the near-equilibrium quark-gluon plasma (QGP) to perturbation at nonhydrodynamic gradients. We propose a conceivable scenario under which sound mode continues to dominate the medium response in this regime. Such a scenario has been observed experimentally for various liquids and liquid metals. We further show that this extended hydrodynamic regime (EHR) indeed exists for both the weakly coupled kinetic equation in the relaxation time approximation (RTA) and the strongly coupled $\mathcal{N}$ = 4 supersymmetric Yang-Mills (SYM) theory. We construct a simple but nontrivial extension of Müller-Israel-Stewart (MIS) theory—namely MIS*—and demonstrate that it describes the EHR response for both the RTA and SYM theory. This indicates that MIS* equations can potentially be employed to search for QGP EHR via heavy-ion collisions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

High Order Strong Stability Preserving MultiDerivative Implicit and IMEX Runge--Kutta Methods with Asymptotic Preserving Properties

In this article we present a class of high order unconditionally strong stability preserving (SSP) implicit two-derivative Runge--Kutta schemes and SSP implicit-explicit (IMEX) multi-derivative Runge--Kutta schemes where the time-step restriction is independent of the stiff term. The unconditional SSP property for a method of order $p>2$ is unique among SSP methods and depends on a backward-in-time assumption on the derivative of the operator. We show that this backward derivative condition is satisfied in many relevant cases where SSP IMEX schemes are desired. We devise unconditionally SSP implicit Runge--Kutta schemes of order up to $p=4$ and IMEX Runge--Kutta schemes of order up to $p=3$. For the multiderivative IMEX schemes, we also derive and present the order conditions, which have not appeared previously. The unconditional SSP condition ensures that these methods are positivity preserving, and we present sufficient conditions under which such methods are also asymptotic preserving when applied to a range of problems, including a hyperbolic relaxation system, the Broadwell model, and the Bhatnagar--Gross--Krook kinetic equation. We present numerical results to support the theoretical results on a variety of problems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Wave-kinetic approach to zonal-flow dynamics: recent advances

Basic physics of drift-wave turbulence and zonal flows has long been studied within the framework of wave-kinetic theory. Recently, this framework has been re-examined from first principles, which has led to more accurate yet still tractable "improved" wave-kinetic equations. In particular, these equations reveal an important effect of the zonal-flow "curvature" (the second radial derivative of the flow velocity) on dynamics and stability of drift waves and zonal flows. We overview these recent findings and present a consolidated high-level picture of (mostly quasilinear) zonal-flow physics within reduced models of drift-wave turbulence.

fusion plasmas↗
DOE OSTI2021

Emu v1.1

Emu is a particle-in-cell code for solving the neutrino quantum kinetic equations in 1, 2, or 3 spatial dimensions with arbitrary angular resolution in order to simulate neutrino flavor transformation in neutron star merger and core collapse supernova environments. Emu represents the neutrino distribution function as a set of particles, each of which represent a collection of neutrinos and antineutrinos with unique position and momentum. Each particle carries two density matrices to define the flavor state of the neutrinos and antineutrinos it represents. Emu includes the vacuum, matter, and neutrino self-interaction potentials. Emu calculates the self-interaction potential using PIC deposition and interpolation algorithms that efficiently compute the local number density and flux at particle locations. Emu is implemented in C++ and is based on the AMReX library for high-performance, block-structured adaptive mesh refinement. Emu is parallelized with MPI + OpenMP for CPUs and MPI + CUDA for GPUs.

Willcox, Donald↗
DOE OSTI2022

Development of a three-dimensional APOLLO3 neutrons deterministic scheme for the CABRI reactor

CABRI is an experimental reactor to study the fuel behavior during reactivity injection transients. These transients being highly multiphysics, the development of suitable modeling and simulation tools to simulate them is important for the optimization of the tests and the control of the experimental conditions. This paper focuses on the development of an APOLLO3 deterministic core calculation tool dedicated to the CABRI transient analysis. It represents the first stage of the incremental process for the implementation of a multiphysics time-dependent modeling of the CABRI transient. The neutron calculation scheme is based on a classical two-step approach. The first step consists of a 281-energy group calculation flux with the TDT-MOC (Method Of Characteristics) solver for cross-section space and energy (23 groups) collapsing for the CABRI different assembly clusters. The bias on a 2D core neutron calculation due to the self-shielding calculation and collapsing on a restricted pattern are investigated thanks to a comparison with a direct full 2D calculation on a quarter of core. The second step relies on a pin-resolved transport 3D transport core calculation with the SN solver MINARET. A progressive numerical validation process is followed to quantify the calculation biases on reactivity and reaction rates at each step using reference calculations with the stochastic code TRIPOLI4. The next development stage toward a multiphysics scheme will be the implementation of the 3D-kinetics equation resolution and the coupling with a core thermal-hydraulics model. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2020

An expanded kinetics model describing injector region desorption

During the initial stages of electron beam injection, absorbed or entrained material can be released into the nascent beam pulse by thermal desorption or electron beam impact. These processes, and subsequent ionization and dissociation, are examined assuming a prototypical temporal beam pulse. An expanded set of kinetic equations that includes thermal desorption, electron-stimulated ion emission, electron impact dissociation, electron impact ionization, and thermal decomposition of adsorbed water molecules is developed. Estimates of the three most important ionized species H 2 O + , H + and OH + are provided for typical injector conditions. Electron impact dissociative ionization of neutral H 2 O is the dominant source of H + under these conditions. Since the abundance of thermally desorbed neutral H 2 O has an exponential rate, the production of H+ is likewise controlled by this exponential desorption rate. Thus thermal management of the confining walls could suppress proton production in the injector region. Thermal decomposition of H 2 O will not be a contributing physical process in the injector region on time scales less than a μs. The formalism is sufficiently flexible to include entrained hydrogen.

42 ENGINEERING↗
DOE OSTI2021

Institutional Computing: Annual Progress Report, Precision early universe simulations to constrain nuclear reactions and beyond standard model physics

The current one-year allocation (w20 qburst) is the first year in a two year allocation of the same name. The project has supported two publications in peer-review journals, another manuscript in preparation for publication, a white paper for the Astro2020 Decadal Survey, and a number of talks (see below). Code development continues on implementing a neutrino quantum kinetic equation solver into BURST, now called QBURST. While implementing coherent and forward-scattering processes into the dynamical neutrino evolution equations, multiple numerical issues have appeared. We have remedied a number of these issues through extensive debugging and testing, which include streamlining interpolation and implementing new ODE integration techniques. Further testing is required to solve the problem at a precision of better than one part in 10 6 . We expect that we will understand the numerical difficulties and implement solutions by end of year 2021.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

INITIAL EXPLORATION OF A NOVEL TRANSIENT ARREST SYSTEM INVOLVING FUEL HEATING

A preliminary analysis on a novel accident response system to diminish the severity of supercritical transients was conducted. The novel accident response system, called the instant shock arrest system, involves using electricity to heat the nuclear fuel at the onset of a large accidental reactivity insertion. This system is specifically designed for reactors with metallic fuel, such that the fuel is capable of conducting electricity, and being resistively heated. A reactor dynamics model of the advanced test reactor was created using the point kinetics equations and a linear reactivity feedback model to simulate how the system would effect the maximum fuel temperatures experienced during the transient. Transients with the instant shock arrest system were compare to those without it. It was found that the instant shock arrest system initially heated the fuel more than the unaffected transient but the negative reactivity inserted from such heating was enough to lower the maximum fuel temperature experienced during the transient. After simulating six different accident scenarios with reactivity insertions ranging from 0.5 \$ to 1.3 \$, it was found that the an optimal system response could reduce peak fuel temperatures during the transient by 3.5\% to 5\%. Furthermore, discussion was given on how the optimal system response could be obtained using relatively simple numerical optimization algorithms due to the smoothness of the optimization problem.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗