Search NASA⌕ Search

SEARCH · Search NASA

Results for “momentum conservation”

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

Jet wake from linearized hydrodynamics

We explore how to improve the hybrid model description of the particles originating from the wake that a jet produced in a heavy ion collision leaves in the droplet of quark-gluon plasma (QGP) through which it propagates, using linearized hydrodynamics on a background Bjorken flow. Jet energy and momentum loss described by the hybrid model become currents sourcing linearized hydrodynamics. By solving the linearized hydrodynamic equations numerically, we investigate the development of the wake in the dynamically evolving droplet of QGP, study the effect of viscosity, scrutinize energy-momentum conservation, and check the validity of the linear approximation. We find that linearized hydrodynamics works better in the viscous case because diffusive modes damp the energy-momentum perturbation produced by the jet. We calculate the distribution of particles produced from the jet wake by using the Cooper-Frye prescription and find that both the transverse momentum spectrum and the distribution of particles in azimuthal angle are similar in shape in linearized hydrodynamics and in the hybrid model. Their normalizations are different because the momentum-rapidity distribution in the linearized hydrodynamics analysis is more spread out, due to sound modes. Since the Bjorken flow has no transverse expansion, we explore the effect of transverse flow by using local boosts to add it into the Cooper-Frye formula. After including the effects of transverse flow in this way, the transverse momentum spectrum becomes harder: more particles with transverse momenta bigger than 2 GeV are produced than in the hybrid model. Although we defer implementing this analysis in a jet Monte Carlo, as would be needed to make quantitative comparisons to data, we gain a qualitative sense of how the jet wake may modify jet observables by computing proxies for two example observables: the lost energy recovered in a cone of varying open angle, and the fragmentation function. We find that linearized hydrodynamics with transverse flow effects added improves the description of the jet wake in the hybrid model in just the way that comparison to data indicates is needed. Our study illuminates a path to improving the description of the wake in the hybrid model, highlighting the need to take into account the effects of both transverse flow and the broadening of the energy-momentum perturbation in spacetime rapidity on particle production.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

DG-IMEX method for a two-moment model for radiation transport in the $\mathscr{O}$($v$/$c$) limit

Here, we consider neutral particle systems described by moments of a phase-space density and propose a realizability-preserving numerical method to evolve a spectral two-moment model for particles interacting with a background fluid moving with nonrelativistic velocities. The system of nonlinear moment equations, with special relativistic corrections to $\mathscr{O}$($v$/$c$), expresses a balance between phase-space advection and collisions and includes velocity-dependent terms that account for spatial advection, Doppler shift, and angular aberration. The model is conservative for the correct $\mathscr{O}$($v$/$c$) Eulerian-frame number density and is consistent, to $\mathscr{O}$($v$/$c$), with Eulerian-frame energy and momentum conservation. This model is closely related to the one promoted by Lowrie et al. and similar to models currently used to study transport phenomena in large-scale simulations of astrophysical environments. The proposed numerical method is designed to preserve moment realizability, which guarantees that the moments correspond to a nonnegative phase-space density. The realizability-preserving scheme consists of the following key components: (i) a strong stability-preserving implicit-explicit (IMEX) time-integration method; (ii) a discontinuous Galerkin (DG) phase-space discretization with carefully constructed numerical uxes; (iii) a realizability-preserving implicit collision update; and(iv) a realizability-enforcing limiter. In time integration, nonlinearity of the moment model necessitates solution of nonlinear equations, which we formulate as fixed-point problems and solve with tailored iterative solvers that preserve moment realizability with guaranteed global convergence. We also analyze the simultaneous Eulerian-frame number and energy conservation properties of the semi-discrete DG scheme and propose a "spectral redistribution" scheme that promotes Eulerian-frame energy conservation. Through numerical experiments, we demonstrate the accuracy and robustness of this DG-IMEX method and investigate its Eulerian-frame energy conservation properties.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Magnetic-Field-Driven Electron Dynamics in Graphene

Graphene exhibits unique optoelectronic properties originating from the band structure at the Dirac points. It is an ideal model structure to study the electronic and optical properties under the influence of the applied magnetic field. In graphene, electric field, laser pulse, and voltage can create electron dynamics which is influenced by momentum dispersion. However, computational modeling of momentum-influenced electron dynamics under the applied magnetic field remains challenging. In this work, we perform computational modeling of the photoexcited electron dynamics achieved in graphene under an applied magnetic field. Our results show that magnetic field leads to local deviation from momentum conservation for charge carriers. With the increasing magnetic field, the delocalization of electron probability distribution increases and forms a cyclotron-like trajectory. Our work facilitates understanding of momentum resolved magnetic field effect on non-equilibrium properties of graphene, which is critical for optoelectronic and photovoltaic applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

A positivity-preserving and conservative intersection-distribution-based remapping algorithm for staggered ALE hydrodynamics on arbitrary meshes

In this study, we introduce new intersection-distribution-based remapping tools for indirect staggered arbitrary Lagrangian-Eulerian (ALE) simulations of multi-material shock hydrodynamics on arbitrary meshes. In addition to conserving momentum and total energy, the three-stage remapper proposed in this work preserves non-negativity of the internal energy. At the first stage, we construct slope-limited piecewise-linear reconstructions of all conserved quantities on zones of the source mesh and perform intersection-based remap to obtain bound-preserving zonal quantities on the target mesh. At the second stage, we define bound-preserving nodal quantities of the staggered ALE discretization as convex combinations of corner quantities. The nodal internal energy is corrected in a way which keeps it non-negative, while providing exact conservation of total energy. At the final stage, we distribute the non-negative nodal internal energy to corners, zones and materials using non-negative weights. Proofs of positivity preservation are provided for each stage. This work is a natural extension of our paper [14] in which a similar intersection-distribution-based remapping procedure was employed. The original version used a nodal kinetic energy fix which did not provably ensure positivity preservation for the zonal internal energy after the final distribution stage. The new algorithm cures this potential drawback by using ‘coordinated’ limiters for piecewise-linear reconstructions, remapping the internal energy to nodes and correcting it before redistribution. The effectiveness of the new nodal fix is illustrated by numerical examples.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Verification of an energy-conserving semi-implicit electrostatic particle-in-cell scheme for modeling high-density plasma at scale

A verification study of a semi-implicit energy-conserving electrostatic particle-in-cell algorithm is presented. The algorithm relaxes the time-step and mesh-size constraints that require resolution of the plasma period and Debye length associated with traditional explicit momentum-conserving particle-in-cell algorithms. Physical implications and applicability of using the semi-implicit scheme for modeling high-density plasmas are discussed. Where possible, numerical results are compared against analytical solutions. The simulation results indicate that the algorithm is stable at time steps larger than twice the inverse plasma frequency and cell sizes larger than the Debye length. It is found that the algorithm gives adequate results, provided that the distribution function and the spatiotemporal scales dictating the physics of the problem are resolved. As such, the algorithm may provide a robust method for kinetic modeling of high-density plasmas at scale.

Cyclotron resonance↗
DOE OSTI2023

Direct implicit and explicit energy-conserving particle-in-cell methods for modeling of capacitively coupled plasma devices

Achieving large-scale kinetic modeling is a crucial task for the development and optimization of modern plasma devices. With the trend of decreasing pressure in applications, such as plasma etching, kinetic simulations are necessary to self-consistently capture the particle dynamics. The standard, explicit, electrostatic, momentum-conserving particle-in-cell method suffers from restrictive stability constraints on spatial cell size and temporal time step, requiring resolution of the electron Debye length and electron plasma period, respectively. This results in a very high computational cost, making the technique prohibitive for large volume device modeling. We investigate the direct implicit algorithm and the explicit energy conserving algorithm as alternatives to the standard approach, both of which can reduce computational cost with a minimal (or controllable) impact on results. These algorithms are implemented into the well-tested EDIPIC-2D and LTP-PIC codes, and their performance is evaluated via 2D capacitively coupled plasma discharge simulations. The investigation reveals that both approaches enable the utilization of cell sizes larger than the Debye length, resulting in a reduced runtime, while incurring only minor inaccuracies in plasma parameters. The direct implicit method also allows for time steps larger than the electron plasma period; however, care must be taken to avoid numerical heating or cooling. It is demonstrated that by appropriately adjusting the ratio of cell size to time step, it is possible to mitigate this effect to an acceptable level.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Cross section measurements of deuteron electro-disintegration at very high recoil momenta and large 4-momentum transfers (Q^2)

The 2H(e,e'p)n cross sections have been measured at negative 4-momentum transfers of Q^2 = 4.5 +/- 0.5 (GeV/c)^2 and Q^2 = 3.5 +/- 0.5 (GeV/c)^2 reaching neutron recoil (missing) momenta up to p_r ?1.0 GeV/c. The data have been obtained at fixed neutron recoil angles 5 <= theta_nq <= 95 degrees with respect to the 3-momentum transfer q. The new data agree well with the previous data which reached p_r ?550 MeV/c. At theta_nq = 35 and 45 degrees, final state interactions (FSI), meson exchange currents (MEC) and isobar configurations (IC) are suppressed and the plane wave impulse approximation (PWIA) provides the dominant cross section contribution. The new data are compared to recent theoretical calculations, and a significant disagreement for recoil momenta p_r > 700 MeV/c is observed. The experiment was carried out in experimental Hall C at the Thomas Jefferson National Accelerator Facility (TJNAF) and formed part of a group of four experiments that were used to commission the new Super High Momentum Spectrometer (SHMS). The experiment consisted of a 10.6 GeV electron beam incident on a liquid deuterium target which resulted in the break-up of the deuteron into a proton and neutron. The scattered electrons were detected by the SHMS in coincidence with the knocked-out protons detected in the previously existing High Momentum Spectrometer (HMS) and the recoiling neutrons were reconstructed from energy-momentum conservation laws. To ensure that the 2H(e,e'p)n reaction channel was selected, we required the missing energy of the system to be the binding energy of the deuteron (?2.22 MeV). The spectrometers? central angles and momenta were set to measure three central missing momentum settings of the neutron corresponding to p_r = 80, 580 and 750 MeV/c, which required the SHMS central angle and momentum to be fixed and the HMS to be rotated from smaller to larger angles corresponding to the lower and higher missing momentum settings, respectively. The experiment was carried out in a time period of six days with typical electron beam currents of 45-60 uA at about 50% beam efficiency.

Yero, Carlos↗
DOE OSTI2021

(U) A Theoretical Study of Asay Foil Trajectories (V.02)

We consider the trajectory of a generic Asay foil ejecta momentum diagnostic sensor, for scenarios where ejecta are produced at a planar surface and fly ballistically through a perfect vacuum to the sensor. To do so, we build upon a previously established mathematical formalism derived for the analytic study of stationary sensors (i.e., piezopins). First, we derive the momentum conservation equation for the problem, in a form amenable to accelerating sensors, in terms of a generic ejecta source areal mass function (“source model”). This defines an integro-differential equation (IDE) for the foil trajectory. When ejecta production is instantaneous - as is generally assumed in momentum diagnostic data analyses - the IDE leads to an implicit and easily calculable closed-form solution for the foil trajectory in a perfect system, as long as the ejecta particle velocity distribution is twice-integrable. General properties of the instant-production trajectory solution indicate the existence of a boundary condition the particle velocity distribution must satisfy in order for the analytically predicted foil trajectories to be compatible with certain features commonly observed in foil data. This boundary condition is identical to one derived previously from a consideration of piezopin data. Armed with the analytic solution for instant production, we also consider various techniques used to extract time-dependent cumulative ejecta masses from foil trajectories, and derive an expression for the error imposed by using an approximated equation of motion. This analytic trajectory solution furthermore makes it possible to examine the common practice of presenting inferred cumulative ejecta masses as a function of a normalized implied velocity; we derive conditions under which this methodology is and is not meaningful. We also propose a strategy for extending the instant-production trajectory solution to time-dependent source functions.

42 ENGINEERING↗
DOE OSTI2026

Generative models on phase space

Deep generative models such as diffusion and flow matching are powerful machine learning tools capable of learning and sampling from high-dimensional distributions. They are particularly useful when the training data appears to be concentrated on a submanifold of the data embedding space. For high-energy physics data, consisting of collections of relativistic energy-momentum 4-vectors, this submanifold can enforce extremely strong physically-motivated priors, such as energy and momentum conservation. If these constraints are learned only approximately, rather than exactly, this can inhibit the interpretability and reliability of such generative models. To remedy this deficiency, we introduce generative models which are, by construction, confined at every step of their sampling trajectory to the manifold of massless N-particle Lorentz-invariant phase space in the center-of-momentum frame. In the case of diffusion models, the "pure noise" forward process endpoint corresponds to the uniform distribution on phase space, which provides a clear starting point from which to identify how correlations among the particles emerge during the reverse (de-noising) process. We demonstrate that our models are able to learn both few-particle and many-particle distributions with various singularity structures, paving the way for future interpretability studies using generative models trained on simulated jet data.

Bogorad, Zachary [Fermilab]↗
DOE OSTI2022

Spinning up the Surface: Evidence for Planetary Engulfment or Unexpected Angular Momentum Transport?

In this paper, we report the potential detection of a nonmonotonic radial rotation profile in a low-mass lower-luminosity giant star. For most low- and intermediate-mass stars, the rotation on the main sequence seems to be close to rigid. As these stars evolve into giants, the core contracts and the envelope expands, which should suggest a radial rotation profile with a fast core and a slower envelope and surface. KIC 9267654, however, seems to show a surface rotation rate that is faster than its bulk envelope rotation rate, in conflict with this simple angular momentum conservation argument. We improve the spectroscopic surface constraint, show that the pulsation frequencies are consistent with the previously published core and envelope rotation rates, and demonstrate that the star does not show strong chemical peculiarities. We discuss the evidence against any tidally interacting stellar companion. Finally, we discuss the possible origin of this unusual rotation profile, including the potential ingestion of a giant planet or unusual angular momentum transport by tidal inertial waves triggered by a close substellar companion, and encourage further observational and theoretical efforts.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Measurement of phonon angular momentum

In condensed matter, angular momentum is intimately related to the emergence of topological quantum states, including chiral superconductivity, quantum spin liquids and various chiral quasiparticles. Recently, it has been predicted that microscopic lattice excitations, known as phonons, can carry finite angular momentum, leading to specific physical properties of materials. However, phonon angular momentum has not yet been observed directly. Here, in this study, we demonstrate that angular momentum conservation results in a macroscopic mechanical torque when applying a time-reversal symmetry-breaking thermal gradient along the chiral axis of single-crystal tellurium. We probe this torque using a cantilever-based device and establish that it changes sign by flipping the chirality or thermal gradient. This behavior disappears in polycrystalline samples that lack a preferred chirality. Our experimental results align well with theoretical calculations. We provide compelling evidence for phonon angular momentum, which might enable quantum states with potential applications in microelectronics.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Spitzer Variability Properties of Low-gravity L Dwarfs

We present Spitzer Space Telescope variability monitoring observations of three low-gravity L dwarfs with previous detections of variability in the near-IR: 2MASS J0045+16, 2MASS J0501–00, and 2MASS J1425–36. We detect significant periodic variability in two of our targets, 2MASS J0045+16 and 2MASS J0501–00. We do not detect variability in 2MASS J1425–36. Combining our new rotation periods with rotational velocities, we calculate inclination angles of 22° ± 1°, 60{sub -8}{sup +13°}, and 52{sub -13}{sup +19°} for 2MASS J0045+16, 2MASS J0501–00, and 2MASS J1425–36, respectively. Our three new objects are consistent with the tentative relations between inclination, amplitude, and color anomaly previously reported. Objects with the highest variability amplitudes are inclined equator on, while the maximum observed amplitude decreases as the inclination angle decreases. We also find a correlation between the inclination angle and (J–K){sub 2MASS} color anomaly for the sample of objects with measured inclinations. Compiling the entire sample of brown dwarfs with Spitzer variability detections, we find no enhancement in amplitude for young, early-L dwarfs compared to the field dwarf population. We find a possible enhancement in amplitude of low-gravity late-L dwarfs at 4.5 μm. We do not find a correlation between amplitude ratio and spectral type for field dwarfs or for the young population. Finally, we compile the rotation periods of a large sample of brown dwarfs with ages 1 Myr–1 Gyr and compare the rotation rates predicted by evolutionary models assuming angular momentum conservation. We find that the rotation rates of the current sample of brown dwarfs fall within the expected range set by evolutionary models and breakup limits.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Examination of flow and nonflow factorization methods in small collision systems

Two-particle correlations have been used extensively to study hydrodynamic flow patterns in heavy-ion collisions. In small collision systems, such as $p + p$ and $p + A$, where particle multiplicities are much smaller than in $A + A$ collisions, nonflow effects from jet correlations, momentum conservation, particle decays, etc. can be significant, even when imposing a large pseudorapidity gap between the particles. A number of techniques to subtract the nonflow contribution in two-particle correlations have been developed by experiments at the CERN Large Hadron Collider (LHC) and then used to measure particle flow in $p + p$ and $p + Pb$ collisions. Recently, experiments at the BNL Relativistic Heavy Ion Collider (RHIC) have explored the possibility of adopting these techniques for small collision systems at lower energies. In this paper, we test these techniques using Monte Carlo generators pythia and hijing, which do not include any collective flow, and ampt, which does. We find that it is crucial to examine the results of such tests both for correlations integrated over particle transverse momentum $p_T$ and differentially as a function of $p_T$. Our results indicate reasonable nonflow subtraction for $p + p$ collisions at the highest LHC energies, while failing if applied to $p + p$ collisions at RHIC. In the case of $p + Au$ collisions at RHIC, both hijing and ampt results indicate a substantial over-estimate of nonflow for $p_T ≳ GeV/c$ and hence an underestimate of elliptic flow and overestimate of triangular flow.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Completely dark matter from rapid-turn multifield inflation

We study cosmological gravitational particle production as applied to “rapid-turn” models of inflation involving two scalar fields. We are interested in the production of massive spin-0 particles that only interact gravitationally and provide a candidate for the dark matter. Specifically, we study two models of rapid-turn multifield inflation, motivated in part by the de Sitter swampland conjecture, that are distinguished by the curvature of field space and the presence or absence of field space ‘angular momentum’ conservation. We find that one of these models leads to insufficient particle production and cannot explain the observed dark matter relic abundance. The second model is able to explain the origin of spin-0 dark matter via gravitational production, and we identify the relevant region of parameter space that is consistent with measurements of the dark-matter relic abundance, the dark-matter-photon isocurvature perturbations, and the spectrum of curvature perturbations that is probed by cosmological observations. Our work demonstrates the compatibility of the de Sitter swampland conjecture with completely dark matter.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

BCJ amplitude relations for Anti-de Sitter boundary correlators in embedding space

We generalize the color/kinematics duality of flat-space scattering amplitudes to the embedding space formulation of AdS boundary correlators. Kinematic numerators and propagators are replaced with differential operators acting on a scalar contact diagram that is the AdS generalization of the momentum conserving delta function of flat space scattering amplitudes. We show that color/kinematics duality implies differential relations among AdS boundary correlators that naturally generalize the flat space BCJ amplitude relations and verify them for the correlators of Yang-Mills theory and of the Nonlinear Sigma Model through four- and six-points, respectively. For the latter we also find representations of the four- and six-point correlator that manifest the duality. Possible double-copy procedures in AdS space are also discussed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2019

Activation of Water by Thorium Cation: A Guided Ion Beam and Quantum Chemical Study

The reaction of atomic thorium cations with deuterated water as a function of kinetic energy from thermal to 10 eV was studied using guided ion beam tandem mass spectrometry. At thermal energies, both ThO{sup +} + D{sub 2} and DThO{sup +} + D are formed in barrierless exothermic processes and reproduce results in the literature obtained using ion cyclotron resonance mass spectrometry. As the energy is increased, the branching ratio between these two channels changes such that the dominant product changes from ThO{sup +} to DThO{sup +} and back to ThO{sup +}, until ThD{sup +} + OD is energetically available and is the dominant product channel. To help understand these experimental results, a variety of theoretical approaches were tried and used to establish a potential energy surface, which compares well with previous theoretical studies. Utilizing the theoretical results, the kinetic energy dependent branching ratio between the ThO{sup +} + D{sub 2} and DThO{sup +} + D channels was calculated using both RRKM and phase space theory (PST). The results indicate that consideration of angular momentum conservation (as in PST) and spin-orbit corrected energies are needed to reproduce experimental results quantitatively. The PST modeling also provides relative energies for the two competing transition states that lead to the primary products, for which theory provides reasonable agreement. .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Action principles for quantum automata and Lorentz invariance of discrete time quantum walks

A discrete action principle for general quantum automata is proposed. This action principle is particularized to Discrete Time Quantum Walks (DTQWs) and then extended into an energy and momentum preserving, manifestly covariant formulation. Space–time coordinates are introduced as new variables of the action and their equations of motion enforce energy and momentum conservation. This guarantees that the proposed action can be used to build future, DTQW-based self-consistent models of spinors interacting with gauge fields. A discrete stress–energy tensor for the DTQW is also obtained. The manifest covariance of the formulation highlights the special role played by the grid reference frame in the DTQW dynamics. The main discussion is complemented by three appendices.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

The hierarchy of excitation lifetimes in two-dimensional Fermi gases

Momentum-conserving quasiparticle collisions in two-dimensional Fermi gases give rise to a large family of exceptionally long-lived excitation modes. The lifetimes of these modes exceed by a factor (T{sub F}∕T){sup 2}≫1 the conventional Landau Fermi-liquid lifetimes τ∼T{sub F}∕T{sup 2}. The long-lived modes have a distinct angular structure, taking the form of cosmθ and sinmθ with odd m values for a circular Fermi surface, with relaxation rate dependence on m of the form m{sup 4}logm, valid at not-too-large m. In contrast, the even-m harmonics feature conventional lifetimes with a weak m dependence. The long-time dynamics, governed by these modes, occurs at “non-Fermi-liquid” timescales, taking values as large as τ∼T{sub F}{sup 3}∕T{sup 4}. The hierarchy of timescales, short and long, leads to unusual long-time memory effects, defining an intriguing transport regime that lies between the conventional ballistic and hydrodynamic regimes.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗