Search NASA⌕ Search

SEARCH · Search NASA

Results for “numerical dispersion”

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

Phonon spectrum in the spin-Peierls phase of CuGeO 3

CuGeO 3 has long been studied as a prototypical example of the spin-Peierls transition in a 𝑆 = 1/2 Heisenberg chain. Despite intensive investigation of this quasi-one-dimensional material, systematic measurements and calculations of the phonon excitations in the dimerized phase have not to date been possible, leaving certain aspects of the spin-Peierls phenomenon unresolved. We perform state-of-the-art density functional theory (DFT) calculations to compute the electronic structure and phonon dynamics in the low-temperature dimerized phase. We also perform high-resolution neutron spectroscopy to measure the full phonon spectrum over multiple Brillouin zones. We find excellent agreement between our numerical and experimental results that extend to all measurement temperatures. Notable features of our phonon spectra include a number of steeply dispersive modes, nonmonotonic dispersion features, and specific phonon anticrossings, which we relate to the mode eigenvectors. By calculating the magnetic interactions within DFT and studying the effects of different phonon modes on the superexchange paths, we discuss the possibility of observing spin-phonon hybridization effects in experiments performed both in and out of equilibrium.

density functional theory↗
DOE OSTI2022

The CO Emission in the Taffy Galaxies (UGC 12914/15) at 60 pc Resolution. I. The Battle for Star Formation in the Turbulent Taffy Bridge

We present Atacama Large Millimeter/submillimeter Array observations at a spatial resolution of $0^{"}_{.}2$ (60 pc) of CO emission from the Taffy galaxies (UGC 12914/5). The observations are compared with narrowband Paα, mid-IR, radio continuum and X-ray imaging, plus optical spectroscopy. The galaxies have undergone a recent head-on collision, creating a massive gaseous bridge that is known to be highly turbulent. The bridge contains a complex web of narrow molecular filaments and clumps. The majority of the filaments are devoid of star formation, and fall significantly below the Kennicutt–Schmidt relationship for normal galaxies, especially for the numerous regions undetected in Paα emission. Within the loosely connected filaments and clumps of gas we find regions of high velocity dispersion that appear gravitationally unbound for a wide range of likely values of X CO . Like the "Firecracker" region in the Antennae system, they would require extremely high external dynamical or thermal pressure to stop them dissipating rapidly on short crossing timescales of 2–5 Myr. We suggest that the clouds may be transient structures within a highly turbulent multiphase medium that is strongly suppressing star formation. Despite the overall turbulence in the system, stars seem to have formed in compact hotspots within a kiloparsec-sized extragalactic H ii region, where the molecular gas has a lower velocity dispersion than elsewhere, and shows evidence for a collision with an ionized gas cloud. Like the shocked gas in the Stephan's Quintet group, the conditions in the Taffy bridge shows how difficult it is to form stars within a turbulent, multiphase, gas.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Magneto-Stokes flow in a shallow free-surface annulus

In this study, we analyse ‘magneto-Stokes’ flow, a fundamental magnetohydrodynamic (MHD) flow that shares the cylindrical-annular geometry of the Taylor–Couette cell but uses applied electromagnetic forces to circulate a free-surface layer of electrolyte at low Reynolds numbers. The first complete, analytical solution for time-dependent magneto-Stokes flow is presented and validated with coupled laboratory and numerical experiments. Three regimes are distinguished (shallow-layer, transitional and deep-layer flow regimes), and their influence on the efficiency of microscale mixing is clarified. The solution in the shallow-layer limit belongs to a newly identified class of MHD potential flows, and thus induces mixing without the aid of axial vorticity. We show that these shallow-layer magneto-Stokes flows can still augment mixing in distinct Taylor dispersion and advection-dominated mixing regimes. The existence of enhanced mixing across all three distinguished flow regimes is predicted by asymptotic scaling laws and supported by three-dimensional numerical simulations. Mixing enhancement is initiated with the least electromagnetic forcing in channels with order-unity depth-to-gap-width ratios. If the strength of the electromagnetic forcing is not a constraint, then shallow-layer flows can still yield the shortest mixing times in the advection-dominated limit. Our robust description of momentum evolution and mixing of passive tracers makes the annular magneto-Stokes system fit for use as an MHD reference flow.

58 GEOSCIENCES↗
DOE OSTI2023

CFL Optimized Forward–Backward Runge–Kutta Schemes for the Shallow-Water Equations

Abstract We present the formulation and optimization of a Runge–Kutta-type time-stepping scheme for solving the shallow-water equations, aimed at substantially increasing the effective allowable time step over that of comparable methods. This scheme, called FB-RK(3,2), uses weighted forward–backward averaging of thickness data to advance the momentum equation. The weights for this averaging are chosen with an optimization process that employs a von Neumann–type analysis, ensuring that the weights maximize the admittable Courant number. Through a simplified local truncation error analysis and numerical experiments, we show that the method is at least second-order in time for any choice of weights and exhibits low dispersion and dissipation errors for well-resolved waves. Further, we show that an optimized FB-RK(3,2) can take time steps up to 2.8 times as large as a popular three-stage, third-order strong stability-preserving Runge–Kutta method in a quasi-linear test case. In fully nonlinear shallow-water test cases relevant to oceanic and atmospheric flows, FB-RK(3,2) outperforms SSPRK3 in admittable time step by factors roughly between 1.6 and 2.2, making the scheme approximately twice as computationally efficient with little to no effect on solution quality. Significance Statement The purpose of this work is to develop and optimize time-stepping schemes for models relevant to oceanic and atmospheric flows. Specifically, for the shallow-water equations we optimize for schemes that can take time steps as large as possible while retaining solution quality. We find that our optimized schemes can take time steps between 1.6 and 2.2 times larger than schemes that cost the same number of floating point operations, translating directly to a corresponding speedup. Our ultimate goal is to use these schemes in climate-scale simulations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Numerical Study of Coupled Water and Vapor Flow, Heat Transfer, and Solute Transport in Variably-Saturated Deformable Soil During Freeze-Thaw Cycles

As climate change intensifies, soil water flow, heat transfer, and solute transport in the active, unfrozen zones within permafrost and seasonally frozen ground exhibit progressively more complex interactions that are difficult to elucidate with measurements alone. For example, frozen conditions impede water flow and solute transport in soil, while heat and mass transfer are significantly affected by high thermal inertia generated from water-ice phase change during the freeze-thaw cycle. To assist in understanding these subsurface processes, the current study presents a coupled two-dimensional model, which examines heat conduction-convection with water-ice phase change, soil water (liquid water and vapor) and groundwater flow, advective-dispersive solute transport with sorption, and soil deformation (frost heave and thaw settlement) in variably saturated soils subjected to freeze-thaw actions. This coupled multiphysics problem is numerically solved using the finite element method. The model's performance is first verified by comparison to a well-documented freezing test on unsaturated soil in a laboratory environment obtained from the literature. Then based on the proposed model, we quantify the impacts of freeze-thaw cycles on the distribution of temperature, water content, displacement history, and solute concentration in three distinct soil types, including sand, silt and clay textures. The influence of fluctuations in the air temperature, groundwater level, hydraulic conductivity, and solute transport parameters was also comparatively studied. The results show that (a) there is a significant bidirectional exchange between groundwater in the saturated zone and soil water in the vadose zone during freeze-thaw periods, and its magnitude increases with the combined influence of higher hydraulic conductivity and higher capillarity; (b) the rapid dewatering ahead of the freezing front causes local volume shrinkage within the non-frozen region when the freezing front propagates downward during the freezing stage and this volume shrinkage reduces the impact of frost heave due to ice formation. This gradually recovers when the thawed water replenishes the water loss zone during the thawing stage; and (c) the profiles of soil moisture, temperature, displacement, and solute concentration during freeze-thaw cycles are sensitive to the changes in amplitude and freeze-thaw period of the sinusoidal varying air temperature near the ground surface, hydraulic conductivity of soil texture, and the initial groundwater levels. Our modeling framework and simulation results highlight the need to account for coupled thermal-hydraulic-mechanical-chemical behaviors to better understand soil water and groundwater dynamics during freeze-thaw cycles and further help explain the observed changes in water cycles and landscape evolution in cold regions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Distinction of Electron Dispersion in Time-Resolved Photoemission Spectroscopy

While recent experiments provided compelling evidence for an intricate dependence of attosecond photoemission-time delays on the solid’s electronic band structure, the extent to which electronic transport and dispersion in solids can be imaged in time-resolved photoelectron (PE) spectra remains poorly understood. Emphasizing the distinction between photoemission time delays measured with two-photon, two-color interferometric spectroscopy, and transport times, we demonstrate how the effect of energy dispersion in the solid on photoemission delays can, in principle, be observed in interferometric photoemission. In this work we reveal analytically a scaling relation between the PE transport time in the solid and the observable photoemission delay and confirm this relation in numerical simulations for a model system. We trace photoemission delays to the phase difference the PE accumulates inside the solid and, in particular, predict negative photoemission delays. Based on these findings, we suggest a novel time-domain interferometric solid-state energy-momentum-dispersion imaging method.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Stochastic properties of ultralight scalar field gradients

Ultralight axion-like particles are well-motivated dark matter candidates that are the target of numerous direct detection efforts. In the vicinity of the Solar System, such particles can be treated as oscillating scalar fields. The velocity dispersion of the Milky Way determines a coherence time of about 10 6 oscillations, beyond which the amplitude of the axion field fluctuates stochastically. Any analysis of data from an axion direct detection experiment must carefully account for this stochastic behavior to properly interpret the results. This is especially true for experiments sensitive to the gradient of the axion field that are unable to collect data for many coherence times. Indeed, the direction, in addition to the amplitude, of the axion field gradient fluctuates stochastically. We present the first complete stochastic treatment for the gradient of the axion field, including multiple computationally efficient methods for performing likelihood-based data analysis, which can be applied to any axion signal, regardless of coherence time. Additionally, we demonstrate that ignoring the stochastic behavior of the gradient of the axion field can potentially result in failure to discover a true axion signal

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

What is the iε for the S-matrix?

Can the S-matrix be complexified in a way consistent with causality? Since the 1960's, the affirmative answer to this question has been well-understood for 2→2 scattering of the lightest particle in theories with a mass gap at low momentum transfer, where the S-matrix is analytic everywhere except at normal-threshold branch cuts. We ask whether an analogous picture extends to realistic theories, such as the Standard Model, that include massless fields, UV/IR divergences, and unstable particles. Especially in the presence of light states running in the loops, the traditional iε prescription for approaching physical regions might break down, because causality requirements for the individual Feynman diagrams can be mutually incompatible. We demonstrate that such analyticity problems are not in contradiction with unitarity. Instead, they should be thought of as finite-width effects that disappear in the idealized 2→2 scattering amplitudes with no unstable particles, but might persist at higher multiplicity. To fix these issues, we propose an iε -like prescription for deforming branch cuts in the space of Mandelstam invariants without modifying the analytic properties. This procedure results in a complex strip around the real part of the kinematic space, where the S-matrix remains causal. In addition to giving a pedagogical introduction to the analytic properties of the perturbative S-matrix from a modern point of view, we illustrate all the points on explicit examples, both symbolically and numerically. Furthermore, to help with the investigation of related questions, we introduce a number of tools, including holomorphic cutting rules, new approaches to dispersion relations, as well as formulae for local behavior of Feynman integrals near branch points.

Anomalous thresholds↗
DOE OSTI2022

A Laser Frequency Transverse Modulation Might Compensate for the Spectral Broadening Due to Large Electron Energy Spread in Thomson Sources

Compact laser plasma accelerators generate high-energy electron beams with increasing quality. When used in inverse Compton backscattering, however, the relatively large electron energy spread jeopardizes potential applications requiring small bandwidths. We present here a novel interaction scheme that allows us to compensate for the negative effects of the electron energy spread on the spectrum, by introducing a transverse spatial frequency modulation in the laser pulse. Such a laser chirp, together with a properly dispersed electron beam, can substantially reduce the broadening of the Compton bandwidth due to the electron energy spread. We show theoretical analysis and numerical simulations for hard X-ray Thomson sources based on laser plasma accelerators.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Modal analysis of the behavior of inertial particles in turbulence subjected to Stokes drag

The clustering of small heavy inertial particles subjected to Stokes drag in turbulence, which is the concern of present study, is known to be minimal at small and large Stokes number and substantial at St = $\textit{O}$(1). This nonmonotonic trend, which has been shown computationally and experimentally, is yet to be explained analytically. In this study, we obtain an analytical expression for the Lyapunov exponents that quantitatively predicts this trend. The sum of the exponents, which is the normalized rate of change of the signed volume of a small cloud of particles, is correctly predicted to be negative and positive at small and large Stokes numbers, respectively, asymptoting to $\textit{τQ}$ as $\textit{τ}$ → 0 and $τ^{–1/2}|Q|^{1/4}$ as $\textit{τ} → ∞$, where $\textit{τ}$ is the particle relaxation time and $\textit{Q(τ)}$ is the difference between the norm of the rotation- and strain-rate tensors computed along the particle trajectory, which remains negative at all $\textit{τ}$ in turbulent flows. Additionally, the trajectory crossing is predicted only in hyperbolic flows where Q < 0 for sufficiently inertial particles with a $\textit{τ}$ that scales with $|Q|^{–1/2}$. Following the onset of crossovers, a transition from clustering to dispersion is predicted correctly. We show these behaviors are not unique to three-dimensional isotropic turbulence and can be reproduced closely by a one-dimensional mono-harmonic flow, which appears as a fundamental canonical problem in the study of particle clustering. Analysis of this one-dimensional canonical flow shows that the rate of clustering, quantified as the product of the Lyapunov exponent and particle relaxation time, is bounded by –1/2, behaving with extreme nonlinearity in the hyperbolic flows and always remaining positive in the elliptic flows. Finally, this model problem also confirms that dispersion in a hyperbolic flow will occur only if particles' trajectory cross. These findings, which stem from our analysis, are corroborated by the direct numerical simulations.

42 ENGINEERING↗
DOE OSTI2023

Ground-state order in magic-angle graphene at filling ν = - 3 : A full-scale density matrix renormalization group study

Here we investigate twisted bilayer graphene (TBG) at filling ν = -3 in the presence of realistic heterostrain. Strain amplifies the band dispersion and drives the system beyond the strong-coupling regime of previous theoretical studies. We use DMRG to conduct an unbiased, large-scale numerical calculations that include all spin and valley degrees of freedom, up to bond dimension χ = 24576. We establish a global phase diagram that unifies a number of theoretical and experimental results. Near zero strain we find an intervalley-coherent quantized anomalous Hall (QAH-IVC) state, a competitive strong-coupling order that evaded past numerical studies. A tiny strain around 0.05% drives a transition into an incommensurate Kekulé spiral (IKS) phase, supporting the mean-field prediction in [Kwan, Phys. Rev. X 11, 041063 (2021)2160-330810.1103/PhysRevX.11.041063]. Even higher strains above 0.2% favor a flavor-symmetric metallic order, which may explain metals found at ν = -3 in many experiments.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Eigenstate Entanglement: Crossover from the Ground State to Volume Laws

For the typical quantum many-body systems that obey the eigenstate thermalization hypothesis (ETH), we argue that the entanglement entropy of (almost) all energy eigenstates is described by a single crossover function. The ETH implies that the crossover functions can be deduced from subsystem entropies of thermal ensembles and have universal properties. These functions capture the full crossover from the ground-state entanglement regime at low energies and small subsystem size (area or log-area law) to the extensive volume-law regime at high energies or large subsystem size. For critical one-dimensional systems, a universal scaling function follows from conformal field theory and can be adapted for nonlinear dispersions. We use it to also deduce the crossover scaling function for Fermi liquids in d > 1 dimensions. Here, the analytical results are complemented by numerics for large noninteracting systems of fermions in d ≤ 3 dimensions and have also been confirmed for bosonic systems and nonintegrable spin chains.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Compact, All-PM Fiber Integrated and Alignment-Free Ultrafast Yb:Fiber NALM Laser With Sub-Femtosecond Timing Jitter

We report a simple and compact design of a dispersion compensated mode-locked Yb:fiber oscillator based on a nonlinear amplifying loop mirror (NALM). The fully polarization maintaining (PM) fiber integrated laser features a chirped fiber Bragg grating (CFBG) for dispersion compensation and a fiber integrated compact non-reciprocal phase bias device, which is alignment-free. The main design parameters were determined by numerically simulating the pulse evolution in the oscillator and by analyzing their impact on the laser performance. Experimentally, we achieved an 88 fs compressed pulse duration with sub-fs timing jitter at 54 MHz repetition rate and 51 mW of output power with 5.5×10 -5 [20 Hz, 1 MHz] integrated relative intensity noise (RIN). Furthermore, we demonstrate tight phase-locking of the laser's carrier-envelope offset frequency ( fceo ) to a stable radio frequency (RF) reference and of one frequency comb tooth to a stable optical reference at 291 THz.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Noncolinear Synthesis of Picosecond Flat-Top Laser Pulses to Increase Electron Photoinjector Brightness

Next-generation electron photoinjector accelerators, such as the LCLS-II photoinjector, have increasingly tight requirements on the excitation lasers, often calling for tens of picosecond, temporally flat-top, ultraviolet (UV) pulse trains to be delivered at up to 1 MHz[1]. We present a numerical and experimental implementation of a non-colinear sum frequency generation scheme wherein the incident optical pulses have equal and opposite amounts of spectral dispersion resulting in a spectrally narrowband pulse with flat-top temporal profile[2]. In preliminary experiments we achieve upwards of 40% conversion efficiency with this process. We additionally demonstrate the ability for the narrowband pulses to be directly upconverted to the UV for driving a photocathode.

Lemons, R.↗
DOE OSTI2022

Supernova-driven Turbulent Metal Mixing in High-redshift Galactic Disks: Metallicity Fluctuations in the Interstellar Medium and its Imprints on Metal-poor Stars in the Milky Way

The extent to which turbulence mixes gas in the face of recurrent infusions of fresh metals by supernovae (SN) could help provide important constraints on the local star formation conditions. This includes predictions of the metallicity dispersion among metal-poor stars, which suggests that the interstellar medium was not very well mixed at these early times. The purpose of this Letter is to help isolate, via a series of numerical experiments, some of the key processes that regulate turbulent mixing of SN elements in galactic disks. We study the gas interactions in small simulated patches of a galaxy disk with the goal of resolving the small-scale mixing effects of metals at parsec scales, which enables us to measure the turbulent diffusion coefficient in various galaxy environments. By investigating the statistics of variations of α elements in these simulations, we are able to derive constraints not only on the allowed range of intrinsic yield variations in SN explosions but also on the star formation history of the Milky Way. We argue that the observed dispersion of [Mg/Fe] in metal-poor halo stars is compatible with the star-forming conditions expected in dwarf satellites or in an early low-star-forming Milky Way progenitor. In particular, metal variations in stars that have not been phase-mixed can be used to infer the star-forming conditions of disrupted dwarf satellites.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Simulation studies for beam commissioning at FRIB Advanced Rare Isotope Separator

The Facility for Rare Isotope Beams (FRIB) includes a powerful superconducting driver accelerator and an Advanced Rare Isotope Separator (ARIS). The ARIS collects and purifies the rare isotope fragments of interest for experiments in nuclear physics, nuclear astrophysics, fundamental symmetries, etc. ARIS consists of a vertical pre-separator and downstream horizontal separator section (C-Bend). Each section can provide a high-resolution separation alternatively. The resolution reduction due to the emittance induced by momentum compression can be avoided by isotope separation in different dispersive planes. Beam commissioning of ARIS for the first experiments was completed and demonstrated particle identification of fragments. The beam tuning in ARIS largely relies on numerical simulations since the limited space for diagnostics. Here, we report the result of the beam trajectory correction, transverse matching, and beam-based misalignment studies at ARIS.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

Observation and Simulation of Methane Plumes During the Morning Boundary Layer Transition

Abstract Methane (CH 4 ) contributes significantly to global warming. However, accurate identification of CH 4 sources for reducing CH 4 emissions is often hampered by inadequate accuracy and spatiotemporal coverage of CH 4 detection, and lack of accurate CH 4 forward modeling used in top‐down inversion systems. In this study, a field experiment was conducted in Pampa, Texas using two CH 4 sensors (LI‐COR and OGI camera) to detect CH 4 releases. We investigated whether high‐resolution simulations using the Weather Research and Forecasting (WRF) model with greenhouse gases (WRF‐GHG) could accurately simulate the CH 4 plumes in the presence of evolving atmospheric boundary layer from sunrise to noon. CH 4 plumes showed substantial variation in time. At a release rate of ∼17.5 kg hr −1 , the maximum enhancement of CH 4 measured by LI‐COR was 2.6 ppm at sunrise (7:36 a.m.), 250 m from the release location. Within half an hour after sunrise, this enhancement decreased to 0.3–0.4 ppm. The enhancement was 0.2 ppm by 10:00 a.m. and further dropped to less than 0.1 ppm after 11:30 a.m. Due to the low temperature at sunrise, the OGI camera failed to detect the CH 4 plume. The WRF‐GHG large‐eddy simulation (LES) with 32 m grid spacing successfully reproduced these CH 4 enhancements. In situ measurements together with numerical simulations illustrate the impact of the transition from a stable boundary layer in the early morning to a convective boundary layer at noon on the dispersion of CH 4 plumes. Additionally, CH 4 plumes from a cattle farm in Oklahoma are briefly examined using the same modeling approach.

Hu, Xiao‐Ming [Center for Analysis and Prediction ↗