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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Nonmagnetic fractional conductance in high mobility InAs quantum point contacts

In this paper, we report the magnetoelectronic properties of high mobility InAs quantum point contacts grown on InP substrates. The InAs quantum well is embedded between In 0.72 ⁢Ga 0.28 ⁢As cladding layers and In 0.81 ⁢Al 0.19 ⁢As barrier layers, and is populated via self-accumulation. The one-dimensional (1D) conductance reaches a maximum value of 17 plateaus, quantized in units of 2⁢𝑒 2 /ℎ, where 𝑒 is the fundamental unit of charge and ℎ is Planck's constant. The in-plane effective 𝑔 factor was estimated to be −10.9 ± 1.5 for subband 𝑁 = 1 and −10.8 ± 1.6 for subband 𝑁 = 2. Furthermore, a study of the nonmagnetic fractional conductance states at 0.2 (𝑒 2 /ℎ) and 0.1 (𝑒 2 /ℎ) is provided. While their origin remains under discussion, evidence suggests that they arise from strong electron-electron interactions and momentum-conserving backscattering between electrons in two distinct channels within the 1D region. This phenomenon may also be interpreted as an entanglement between the two channel directions facilitated by momentum-conserving backscattering.

Ballistic transport↗
DOE OSTI2020

Finite spatial-grid effects in energy-conserving particle-in-cell algorithms

Finite-grid (or aliasing) instabilities are pervasive in particle-in-cell (PIC) plasma simulation algorithms, and force the modeler to resolve the smallest (Debye) length scale in the problem regardless of dynamical relevance. These instabilities originate in the aliasing of interpolation errors between mesh quantities and particles (which live in the space–time continuum). Recently, strictly energy-conserving PIC (EC-PIC) algorithms have been developed that promise enhanced robustness against aliasing instabilities. In this study, we confirm by analysis that EC-PIC is stable against aliasing instabilities for stationary plasmas. For drifting plasmas, we demonstrate by analysis and numerical experiments that, while EC-PIC algorithms are not free from these instabilities in principle, they feature a benign stability threshold for finite-temperature plasmas that make them usable in practice for a large class of problems (featuring ambipolarity and realistic ion-electron mass ratios) without the need to consider the size of the Debye length. We also demonstrate that this threshold is absent for the popular momentum-conserving PIC algorithms, which are therefore unstable for both drifting and stationary plasmas beyond a threshold in cell size compared to Debye length. Finite-grid (or aliasing) instabilities are pervasive in particle-in-cell (PIC) plasma simulation algorithms, and force the modeler to resolve the smallest (Debye) length scale in the problem regardless of dynamical relevance. These instabilities originate in the aliasing of interpolation errors between mesh quantities and particles (which live in the space–time continuum). Recently, strictly energy-conserving PIC (EC-PIC) algorithms have been developed that promise enhanced robustness against aliasing instabilities. In this study, we confirm by analysis that EC-PIC is stable against aliasing instabilities for stationary plasmas. For drifting plasmas, we demonstrate by analysis and numerical experiments that, while EC-PIC algorithms are not free from these instabilities in principle, they feature a benign stability threshold for finite-temperature plasmas that make them usable in practice for a large class of problems (featuring ambipolarity and realistic ion-electron mass ratios) without the need to consider the size of the Debye length. Finally, we also demonstrate that this threshold is absent for the popular momentum-conserving PIC algorithms, which are therefore unstable for both drifting and stationary plasmas beyond a threshold in cell size compared to Debye length.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Landau diagrams in AdS and S-matrices from conformal correlators

Quantum field theories in AdS generate conformal correlation functions on the boundary, and in the limit where AdS is nearly flat one should be able to extract an S-matrix from such correlators. We discuss a particularly simple position-space procedure to do so. It features a direct map from boundary positions to (on-shell) momenta and thereby relates cross ratios to Mandelstam invariants. This recipe succeeds in several examples, includes the momentum-conserving delta functions, and can be shown to imply the two proposals in based on Mellin space and on the OPE data. Interestingly the procedure does not always work: the Landau singularities of a Feynman diagram are shown to be part of larger regions, to be called ‘bad regions’, where the flat-space limit of the Witten diagram diverges. To capture these divergences we introduce the notion of Landau diagrams in AdS. As in flat space, these describe on-shell particles propagating over large distances in a complexified space, with a form of momentum conservation holding at each bulk vertex. As an application we recover the anomalous threshold of the four-point triangle diagram at the boundary of a bad region.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

BSE versus StarTrack: Implementations of new wind, remnant-formation, and natal-kick schemes in NBODY7 and their astrophysical consequences

As a result of their formation via massive single and binary stellar evolution, the masses of stellar-remnant black holes (BH) are subjects of great interest in this era of gravitational-wave detection from binary black hole (BBH) and binary neutron star merger events. In this work, we present new developments in the stellar-remnant formation and related schemes of the current N-body evolution program NBODY7. We demonstrate that the newly implemented stellar-wind and remnant-formation schemes in the stellar-evolutionary sector or BSE of the NBODY7 code, such as the “rapid” and the “delayed” supernova (SN) schemes along with an implementation of pulsational-pair-instability and pair-instability supernova (PPSN/PSN), now produce neutron star (NS) and BH masses that agree nearly perfectly, over large ranges of zero-age-main-sequence (ZAMS) mass and metallicity, with those from the widely recognised StarTrack population-synthesis program. We also demonstrate the new, recipe-based implementations of various widely debated mechanisms of natal kicks on NSs and BHs, such as “convection-asymmetry-driven”, “collapse-asymmetry-driven”, and “neutrino-emission-driven” kicks, in addition to a fully consistent implementation of the standard, fallback-dependent, momentum-conserving natal kick. All the above newly implemented schemes are also shared with the standalone versions of SSE and BSE. Furthermore, all these demonstrations are performed with both the updated standalone BSE and the updated NBODY7/BSE. When convolved with stellar and primordial-binary populations as observed in young massive clusters, such remnant-formation and natal-kick mechanisms crucially determine the accumulated number, mass, and mass distribution of the BHs retained in young massive, open, and globular clusters (GCs); these BHs would eventually become available for long-term dynamical processing. Among other conclusions, we find that although the newer, delayed SN remnant formation model gives birth to the largest number (mass) of BHs, the older remnant-formation schemes cause the largest number (mass) of BHs to survive in clusters, when incorporating SN material fallback onto the BHs. The SN material fallback also causes the convection-asymmetry-driven SN kick to effectively retain similar numbers and masses of BHs in clusters as for the standard, momentum-conserving kick. The collapse-asymmetry-driven SN kick would cause nearly all BHs to be retained in clusters irrespective of their mass, remnant-formation model, and metallicity, whereas the inference of a large population of BHs in GCs would potentially rule out the neutrino-driven SN kick mechanism. Pre-SN mergers of massive primordial binaries would potentially cause BH masses to deviate from the theoretical, single-star ZAMS to mass-remnant mass relation unless a substantial of the total merging stellar mass of up to ≈40% is lost during a merger process. In particular, such mergers, at low metallicities, have the potential to produce low-spinning BHs within the PSN mass gap that can be retained in a stellar cluster and be available for subsequent dynamical interactions. As recent studies indicate, the new remnant-formation modelling reassures us that young massive and open clusters would potentially contribute to the dynamical BBH merger detection rate to a similar extent as their more massive GC counterparts.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Local spin polarizations in relativistic heavy-ion collisions

Here, based on a generalized side-jump formalism for massless chiral fermions, which naturally takes into account the spin-orbit coupling in the scattering of two chiral fermions and the chiral vortical effect in a rotating chiral fermion matter, we have developed a covariant and total-angular-momentum-conserved chiral transport model to study both the global and local polarizations of this matter. For a system of massless quarks of random spin orientations and finite vorticity in a box, we have demonstrated that the model can exactly conserve the total angular momentum of the system and dynamically generate the quark spin polarization expected from a thermally equilibrated quark matter. Using this model to study the spin polarization in relativistic heavy-ion collision, we have found that the local quark spin polarizations depend strongly on the reference frame where they are evaluated as a result of the nontrivial axial charge distribution caused by the chiral vortical effect. We have further shown that because of the anomalous orbital or side-jump contribution to the quark spin polarization, the local quark polarizations calculated in the medium rest frame are qualitatively consistent with the local polarizations of Lambda hyperons measured in experiments.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

X-ray Parametric Down Conversion at an XFEL

Understanding and harnessing X-ray quantum effects could open new frontiers in imaging and quantum optics. In this study, we measured the process of X-ray Parametric Down Conversion, where a single high-energy X-ray photon splits into two lower-energy photons. Using the SACLA X-ray Free Electron Laser in Japan at 9.83 keV, we found clear evidence that pairs of photons were produced along the energy-angle relationship that conserved both energy and momentum, as predicted for down-conversion, and consistent with quantum entanglement of X-ray photons. By matching specific photon pairs for energy and momentum conservation, we observe a signal rate of 1,250 pairs per hour, confirming that correlated photon pairs can be generated and observed in the absence of explicit time correlations. Our results show that with further refinement, the number of entangled photons produced per laser pulse could increase by an order of magnitude. This paves the way for demonstrating quantum-enhanced X-ray imaging, and confirmation of X-ray photon entanglement.

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

Effect of Cathode Cooling in Three-Dimensional Simulations of an Atmospheric Pressure Glow Discharge

The Atmospheric Pressure Glow Discharge (APGD) is a relatively simple and versatile plasma source used in a wide range of applications. Active cooling of the cathode can effectively mitigate instabilities, leading to glow-to-arc transitions. This study investigates the effect of varying the degree of cathode cooling in APGD with a planar cathode in helium. The plasma flow model incorporates mass conservation, chemical species transport, momentum conservation, conservation of thermal energy of heavy species and of electrons, and electrostatics. The model is applied to time-dependent simulations through a three-dimensional computational domain describing the whole discharge, without geometric symmetry or steady-state assumptions. Simulations of an experimentally characterized APGD explore the effects of electric current and cathode cooling—ranging from thermally insulated to extreme convective cooling. Results show the formation of an annular region with high electric field over the cathode surface under conditions of high current and low cooling.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Conservation of angular momentum in ultrafast spin dynamics

The total angular momentum of a closed system is a conserved quantity, which should remain constant in time for any excitation experiment once the pumping signal has been extinguished. Such conservation, however, is never satisfied in practice in any real-time first-principles description of the demagnetization process. Furthermore, there is growing experimental evidence that the same takes place in experiments. The missing angular momentum is usually associated to lattice vibrations, which are not measured experimentally and are never considered in real-time simulations. Here, in this study, we critically analyze the issue and conclude that current state-of-the-art simulations violate angular momentum conservation already at the electronic level of description. This shortcoming originates from an oversimplified description of the spin-orbit coupling, which includes atomic contributions but neglects completely that of itinerant electrons. We corroborate our findings with time-dependent simulations using model tight-binding Hamiltonians, and show that indeed such conservation can be reintroduced by an appropriate choice of spin-orbit coupling. The consequences of our findings on recent experiments are also discussed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Axisymmetric hydrodynamics in numerical relativity: treating coordinate singularity, artificial heating and modeling MHD instabilities

Two-dimensional axisymmetric simulations of binary neutron star (BNS) merger remnant are a cheap alternative to 3D simulations. To maintain realism for secular timescales, simulations must avoid accumulated errors from drifts in conserved quantities and artificial heating, and they must model turbulent transport in a way that remains plausible throughout the evolution. It is also crucial to avoid numerical artifacts due to the polar coordinate axis singularity. Methods that behave well near the axis often break flux-conservative form of the hydrodynamic equations, resulting in significant drifts in conserved quantities. We present a flux-conservative scheme that maintains smoothness near the axis without sacrificing conservative formulation of the equations or incurring drifts in conserved global quantities. We compare the numerical performance of different treatments of the hydrodynamic equations when evolving a hypermassive neutron star resembling the remnant of a BNS merger. These simulations demonstrate that the new scheme combines the axis smoothness of non-conservative methods with the mass and angular momentum conservation of other conservative methods on $\sim 10^2$ ms timescales of viscous and neutrino-driven evolution. Because fluid profiles remain smooth in the remnant interior, it is possible to remove artificial heating by evolving the entropy density. We show how physical heating and cooling terms can be easily calculated from source terms of the conservative evolution variables and demonstrate our implementation. Finally, we discuss and implement improvements to the effective viscosity scheme to better model the effect of magnetohydrodynamic instabilities as the remnant evolves.

axisymmetry↗
DOE OSTI2021

Conservation of angular momentum in an elastic medium with spins

Exact conservation of the angular momentum is worked out for an elastic medium with spins. The intrinsic anharmonicity of the elastic theory is shown to be crucial for conserving the total momentum. As a result, any spin-lattice dynamics inevitably involves multiphonon processes and interaction between phonons. Furthermore, this makes transitions between spin states in a solid fundamentally different from transitions between atomic states in a vacuum governed by linear electrodynamics. Consequences for using solid-state spins as qubits are discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Conservation of angular momentum in second harmonic generation from under-dense plasmas

Spin and orbital angular momentum of an optical beam are two independent parameters that exhibit distinct effects on mechanical objects. However, when laser beams with angular momentum interact with plasmas, one can observe the interplay between the spin and the orbital angular momentum. Here, by measuring the helical phase of the second harmonic 2ω radiation generated in an underdense plasma using a known spin and orbital angular momentum pump beam, we verify that the total angular momentum of photons is conserved and observe the conversion of spin to orbital angular momentum. We further determine the source of the 2ω photons by analyzing near field intensity distributions of the 2ω light. The 2ω images are consistent with these photons being generated near the largest intensity gradients of the pump beam in the plasma as predicted by the combined effect of spin and orbital angular momentum when Laguerre-Gaussian beams are used.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Linearised Fokker–Planck collision model for gyrokinetic simulations

We introduce a gyrokinetic, linearised Fokker–Planck collision model that satisfies conservation laws and is accurate at arbitrary collisionalities. The differential test-particle component of the operator is exact; the integral field-particle component is approximated using a spherical harmonic and a modified Laguerre polynomial expansion developed by Hirshman and Sigmar (1976 Phys. Fluids 19 1532). The numerical methods of the implementation in the δf-gyrokinetic code stella (Barnes et al 2019 J. Comput. Phys. 391 365–80) are discussed, and conservation properties of the operator are demonstrated. The collision model is then benchmarked against the collision model of the gyrokinetic solver GS2 in the limiting cases of a reduced test-particle collision operator and energy- and momentum-conserving operator. The accuracy of the full collision model is investigated by solving the parallel Spitzer-Härm problem for the transport coefficients. It is shown that retaining collisional energy flux and higher-order terms in the field-particle operator reduces errors in the transport coefficients from 10%–25% for a simple momentum- and energy-conserving model to under 1%.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Conservative discontinuous Galerkin schemes for nonlinear Dougherty–Fokker–Planck collision operators

Here, we present a novel discontinuous Galerkin algorithm for the solution of a class of Fokker–Planck collision operators. These operators arise in many fields of physics, and our particular application is for kinetic plasma simulations. In particular, we focus on an operator often known as the ‘Lenard–Bernstein’ or ‘Dougherty’ operator. Several novel algorithmic innovations, based on the concept of weak equality, are reported. These weak equalities are used to define weak operators that compute primitive moments, and are also used to determine a reconstruction procedure that allows an efficient and accurate discretization of the diffusion term. We show that when two integrations by parts are used to construct the discrete weak form, and finite velocity-space extents are accounted for, a scheme that conserves density, momentum and energy exactly is obtained. One novel feature is that the requirements of momentum and energy conservation lead to unique formulas to compute primitive moments. Careful definition of discretized moments also ensure that energy is conserved in the piecewise linear case, even though the kinetic-energy term, $v^2$ is not included in the basis set used in the discretization. A series of benchmark problems is presented and shows that the scheme conserves momentum and energy to machine precision. Empirical evidence also indicates that entropy is a non-decreasing function. The collision terms are combined with the Vlasov equation to study collisional Landau damping and plasma heating via magnetic pumping.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Momentum-exchange current drive by electrostatic waves in an unmagnetized collisionless plasma

For a planar electrostatic wave interacting with a single species in a collisionless plasma, momentum conservation implies current conservation. However, when multiple species interact with the wave, they can exchange momentum, leading to current drive. A simple, general formula for this driven current is derived. As examples, we show how currents can be driven for Langmuir waves in electron–positron–ion plasmas, and for ion-acoustic waves in electron–ion plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Probing early-time longitudinal dynamics with the Λ hyperon's spin polarization in relativistic heavy-ion collisions

Here, we systematically study the hyperon global polarization's sensitivity to the collision systems' initial longitudinal flow velocity in hydrodynamic simulations. By explicitly imposing local energy-momentum conservation when mapping the initial collision geometry to macroscopic hydrodynamic fields, we study the evolution of systems' orbital angular momentum (OAM) and fluid vorticity. We find that a simultaneous description of the $\Lambda$ hyperons' global polarization and the slope of pion's directed flow can strongly constrain the size of longitudinal flow at the beginning of hydrodynamic evolution. We extract the size of the initial longitudinal flow and the fraction of orbital angular momentum in the produced QGP fluid as a function of collision energy with the STAR measurements in the RHIC Beam Energy Scan program. We find that there is about 100-200 $\hbar$ OAM that remains in the mid-rapidity fluid at the beginning of hydrodynamic evolution. We further exam the effects of different hydrodynamic gradients on the spin polarization of $\Lambda$ and $\bar{\Lambda}$. The gradients of $\mu_B/T$ can change the ordering between $\Lambda$'s and $\bar{\Lambda}$'s polarization.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Realizability-preserving discontinuous Galerkin method for spectral two-moment radiation transport in special relativity

Here we present a realizability-preserving numerical method for solving a spectral two-moment model to simulate the transport of massless, neutral particles interacting with a steady background material moving with relativistic velocities. The model is obtained as the special relativistic limit of a four-momentum-conservative general relativistic two-moment model. Using a maximum-entropy closure, we solve for the Eulerian-frame energy and momentum. The proposed numerical method is designed to preserve moment realizability, which corresponds to moments defined by a nonnegative phase-space density. The realizability-preserving method is achieved with the following key components: (i) a discontinuous Galerkin phase-space discretization with specially constructed numerical fluxes in the spatial and energy dimensions; (ii) a strong stability-preserving implicit-explicit time-integration method; (iii) a realizability-preserving conserved to primitive moment solver; (iv) a realizability-preserving implicit collision solver; and (v) a realizability-enforcing limiter. Component (iii) is necessitated by the closure procedure, which closes higher order moments nonlinearly in terms of primitive moments. The nonlinear conserved to primitive and the implicit collision solves are formulated as fixed-point problems, which are solved with custom iterative solvers designed to preserve the realizability of each iterate. With a series of numerical tests, we demonstrate the accuracy and robustness of this discontinuous-Galerkin-implicit-explicit method.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Gauge-symmetrization method for energy-momentum tensors in high-order electromagnetic field theories

For electromagnetic field theories, canonical energy-momentum conservation laws can be derived from the underpinning spacetime translation symmetry according to the Noether procedure. However, the canonical energy-momentum tensors (EMTs) are neither symmetric nor gauge-symmetric (gauge invariant). The Belinfante-Rosenfeld (BR) method is a well-known procedure to symmetrize the EMTs, which also renders them gauge symmetric for first-order field theories. High-order electromagnetic field theories appear in the study of gyrokinetic systems for magnetized plasmas and the Podolsky system for the radiation reaction of classical charged particles. For these high-order field theories, gauge-symmetric EMTs are not necessarily symmetric and vice versa. In the present study, we develop a new gauge-symmetrization method for EMTs in high-order electromagnetic field theories. The Noether procedure is carried out using the Faraday tensor $F_{μν}$, instead of the 4-potential $A_{μ}$, to derive a canonical EMT $T^{μν}_{N}$. We show that the gauge-dependent part of $T^{μν}_{N}$ can be removed using the displacement-potential tensor $F^{σμν}$ ≡ $D^{σμ}A^{ν}/4π$, where $D^{σμ}$ is the antisymmetric electric displacement tensor. This method gauge-symmetrizes the EMT without necessarily making it symmetric, which is adequate for applications not involving general relativity. For first-order electromagnetic field theories, such as the standard Maxwell system, $F^{σμν}$ reduces to the familiar BR superpotential $S^{σμν}$, and the method developed can be used as a simpler procedure to calculate $S^{σμν}$ without employing the angular momentum tensor in 4D spacetime. When the electromagnetic system is coupled to classical charged particles, the gauge-symmetrization method for EMTs is shown to be effective as well.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Investigating the physical mechanisms that modify wind plant blockage in stable boundary layers

Abstract. Wind plants slow down the approaching wind, a phenomenon known as blockage. Wind plant blockage undermines turbine performance for front-row turbines and potentially for turbines deeper into the array. We use large-eddy simulations to characterize blockage upstream of a finite-size wind plant in flat terrain for different atmospheric stability conditions and investigate the physical mechanisms modifying the flow upstream of the turbines. To examine the influence of atmospheric stability, we compare simulations of two stably stratified boundary layers using the Weather Research and Forecasting model in large-eddy simulation mode, representing wind turbines using the generalized actuator disk approach. For a wind plant, a faster cooling rate at the surface, which produces stronger stably stratified flow in the boundary layer, amplifies blockage. As a novelty, we investigate the physical mechanisms amplifying blockage by evaluating the different terms in the momentum conservation equation within the turbine rotor layer. The velocity deceleration upstream of a wind plant is caused by an adverse pressure gradient and momentum advection out of the turbine rotor layer. The cumulative deceleration of the flow upstream of the front-row turbines instigates vertical motions. The horizontal flow is diverted vertically, reducing momentum availability in the turbine rotor layer. Although the adverse pressure gradient upstream of the wind plant remains unchanged with atmospheric stability, vertical advection of horizontal momentum is amplified in the more strongly stable boundary layer, mainly by larger shear of the horizontal velocity, thus increasing the blockage effect.

17 WIND ENERGY↗