Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dispersion relation”

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

Finite-volume quantization condition from the N/D representation

We propose a new model-independent method for determining hadronic resonances from lattice QCD. The formalism is derived from the general principles of unitarity and analyticity, as encoded in the N/D representation of a partial-wave two-body amplitude. The associated quantization condition relates the finite-volume spectrum to the infinite-volume numerator $\mathscr{N}$, used to reconstruct the scattering amplitude from dispersive relations. Unlike the original Lüscher condition, this new formalism is valid for energies coinciding with the left-hand cuts from arbitrary one- and multiparticle exchanges.

Dawid, Sebastian M. [Univ. of Washington, Seattle,↗
DOE OSTI2021

Constraints on quasinormal modes and bounds for critical points from pole-skipping

We consider a holographic thermal state and perturb it by a scalar operator whose associated real-time Green’s function has only gapped poles. These gapped poles correspond to the non-hydrodynamic quasinormal modes of a massive scalar perturbation around a Schwarzschild black brane. Relations between pole-skipping points, critical points and quasinormal modes in general emerge when the mass of the scalar and hence the dual operator dimension is varied. First, this novel analysis reveals a relation between the location of a mode in the infinite tower of quasinormal modes and the number of pole-skipping points constraining its dispersion relation at imaginary momenta. Second, for the first time, we consider the radii of convergence of the derivative expansions about the gapped quasinormal modes. These convergence radii turn out to be bounded from above by the set of all pole-skipping points. Furthermore, a transition between two distinct classes of critical points occurs at a particular value for the conformal dimension, implying close relations between critical points and pole-skipping points in one of those two classes. We show numerically that all of our results are also true for gapped modes of vector and tensor operators.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

An FFT-based approach for Bloch wave analysis: application to polycrystals

A method based on the Fast Fourier Transform is proposed to obtain the dispersion relation of acoustic waves in heterogeneous periodic media with arbitrary microstructures. The microstructure is explicitly considered using a voxelized Representative Volume Element (RVE). The dispersion diagram is obtained solving an eigenvalue problem for Bloch waves in Fourier space. To this aim, two linear operators representing stiffness and mass are defined through the use of differential operators in Fourier space. The smallest eigenvalues are obtained using the implicitly restarted Lanczos and the subspace iteration methods, and the required inverse of the stiffness operator is done using the conjugate gradient with a preconditioner. The method is used to study the propagation of acoustic waves in elastic polycrystals, showing the strong effect of crystal anistropy and polycrystaline texture on the propagation. It is shown that the method combines the simplicity of classical Fourier series analysis with the versatility of Finite Elements to account for complex geometries proving an efficient and general approach which allows the use of large RVEs in 3D.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Terahertz-field activation of polar skyrons

Unraveling collective modes arising from coupled degrees of freedom is crucial for understanding complex interactions in solids and developing new functionalities. Unique collective behaviors emerge when two degrees of freedom, ordered on distinct length scales, interact. Polar skyrmions, three-dimensional electric polarization textures in ferroelectric superlattices, disrupt the lattice continuity at the nanometer scale with nontrivial topology, leading to previously unexplored collective modes. Here, using terahertz-field excitation and femtosecond x-ray diffraction, we discover subterahertz collective modes, dubbed “skyrons”, which appear as swirling patterns of atomic displacements functioning as atomic-scale gearsets. The key to activating skyrons is the use of the THz field that couples primarily to skyrmion domain walls. Momentum-resolved time-domain measurements of diffuse scattering reveal an avoided crossing in the dispersion relation of skyrons. Atomistic simulations and dynamical phase-field modeling provide microscopic insights into the three-dimensional crystallographic and polarization dynamics. The amplitude and dispersion of skyrons are demonstrated to be controlled by sample temperature and electric-field bias. The discovery of skyrons and their coupling with terahertz fields opens avenues for ultrafast control of topological polar structures.

ferroelectrics↗
DOE OSTI2020

Distributions in CFT. Part I. Cross-ratio space

We show that the four-point functions in conformal field theory are defined as distributions on the boundary of the region of convergence of the conformal block expansion. The conformal block expansion converges in the sense of distributions on this boundary, i.e. it can be integrated term by term against appropriate test functions. This can be interpreted as a giving a new class of functionals that satisfy the swapping property when applied to the crossing equation, and we comment on the relation of our construction to other types of functionals. Our language is useful in all considerations involving the boundary of the region of convergence, e.g. for deriving the dispersion relations. We establish our results by elementary methods, relying only on crossing symmetry and the standard convergence properties of the conformal block expansion. This is the first in a series of papers on distributional properties of correlation functions in conformal field theory.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Quantum field theory and the Bieberbach conjecture

An intriguing correspondence between ingredients in geometric function theory related to the famous Bieberbach conjecture (de Branges' theorem) and the non-perturbative crossing symmetric representation of 2-2 scattering amplitudes of identical scalars is pointed out. Using the dispersion relation and unitarity, we are able to derive several inequalities, analogous to those which arise in the discussions of the Bieberbach conjecture. We derive new and strong bounds on the ratio of certain Wilson coefficients and demonstrate that these are obeyed in one-loop Φ 4 theory, tree level string theory as well as in the S-matrix bootstrap. Further, we find two sided bounds on the magnitude of the scattering amplitude, which are shown to be respected in all the contexts mentioned above. Translated to the usual Mandelstam variables, for large |s|, fixed t, the upper bound reads |M(s,t)| ≲ |s 2 |. We discuss how Szegö's theorem corresponds to a check of univalence in an EFT expansion, while how the Grunsky inequalities translate into nontrivial, nonlinear inequalities on the Wilson coefficients.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

The Kelvin–Helmholtz Instability at the Boundary of Relativistic Magnetized Jets

We study the linear stability of a planar interface separating two fluids in relative motion, focusing on conditions appropriate for the boundaries of relativistic jets. The jet is magnetically dominated, whereas the ambient wind is gas-pressure-dominated. We derive the most general form of the dispersion relation and provide an analytical approximation of its solution for an ambient sound speed much smaller than the jet Alfvén speed v A , as appropriate for realistic systems. The stability properties are chiefly determined by the angle ψ between the wavevector and the jet magnetic field. For ψ = π/2, magnetic tension plays no role, and our solution resembles the one of a gas-pressure-dominated jet. Here, only sub-Alfvénic jets are unstable ($0\lt {M}_{e}\equiv (v/{v}_{{\rm{A}}})\cos \theta \lt 1$ , where v is the shear velocity and θ the angle between the velocity and the wavevector). For ψ = 0, the free energy in the velocity shear needs to overcome the magnetic tension, and only super-Alfvénic jets are unstable ($1\lt {M}_{e}\lt \sqrt{{(1+{{\rm{\Gamma }}}_{w}^{2})/[1+({v}_{{\rm{A}}}/c)}^{2}{{\rm{\Gamma }}}_{w}^{2}]}$, with Γ w the wind adiabatic index). Our results have important implications for the propagation and emission of relativistic magnetized jets.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Observed and Simulated Variability of Droplet Spectral Dispersion in Convective Clouds Over the Amazon

Abstract In this study, the variability of the spectral dispersion of droplet size distributions (DSDs) in convective clouds is investigated. Analyses are based on aircraft measurements of growing cumuli near the Amazon basin, and on numerical simulations of an idealized ice‐free cumulus. In cleaner clouds, the relative dispersion , defined as the ratio of the standard deviation to the mean value of the droplet diameter, is negatively correlated with the ratio of the cloud water content ( ) to the adiabatic liquid water content ( ), while no strong correlation between and is seen in polluted clouds. Bin microphysics numerical simulations suggest that these contrasting behaviors are associated with the effect of collision‐coalescence in cleaner clouds, and secondary droplet activation in polluted clouds, in addition to the turbulent mixing of parcels that experienced different paths within the cloud. Collision‐coalescence simultaneously broadens the DSDs and decreases , explaining the inverse relationship between and in cleaner clouds. Secondary droplet activation broadens the DSDs but has little direct impact on . The combination of a rather modest DSD broadening due to weak collision‐coalescence with enhanced droplet activation in both diluted and highly undiluted cloud regions may contribute to maintain a relatively uniform within polluted clouds. These findings can be useful for parameterizing the shape parameter ( ) of gamma DSDs in bulk microphysics cloud‐resolving models. It is shown that emulating the observed relationship improves the estimation of the collision‐coalescence rate in bulk microphysics simulations compared to the bin simulations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

The anomalous magnetic moment of the muon in the Standard Model

We review the present status of the Standard Model calculation of the anomalous magnetic moment of the muon. This is performed in a perturbative expansion in the fine-structure constant α and is broken down into pure QED, electroweak, and hadronic contributions. The pure QED contribution is by far the largest and has been evaluated up to and including O(α 5 ) with negligible numerical uncertainty. The electroweak contribution is suppressed by (m μ / M W ) 2 and only shows up at the level of the seventh significant digit. It has been evaluated up to two loops and is known to better than one percent. Hadronic contributions are the most difficult to calculate and are responsible for almost all of the theoretical uncertainty. The leading hadronic contribution appears at O(α 2 ) and is due to hadronic vacuum polarization, whereas at O(α 3 ) the hadronic light-by-light scattering contribution appears. Given the low characteristic scale of this observable, these contributions have to be calculated with nonperturbative methods, in particular, dispersion relations and the lattice approach to QCD. The largest part of this review is dedicated to a detailed account of recent efforts to improve the calculation of these two contributions with either a data-driven, dispersive approach, or a first- principle, lattice-QCD approach. The final result reads $a$$^{SM}_{μ}$ = 116 591 810(43) × 10 −11 and is smaller than the Brookhaven measurement by 3.7σ. The experimental uncertainty will soon be reduced by up to a factor four by the new experiment currently running at Fermilab, and also by the future J-PARC experiment. Furthermore, this and the prospects to further reduce the theoretical uncertainty in the near future – which are also discussed here – make this quantity one of the most promising places to look for evidence of new physics.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Multifield ultralight dark matter

Ultralight dark matter (ULDM) is usually taken to be a single scalar field. Here we explore the possibility that ULDM consists of N light scalar fields with only gravitational interactions. This configuration is more consistent with the underlying particle physics motivations for these scenarios than a single ultralight field. ULDM halos have a characteristic granular structure that increases stellar velocity dispersion and can be used as observational constraints on ULDM models. In multifield simulations, we find that inside a halo the amplitude of the total density fluctuations decreases as $1/\sqrt{N}$ and that the fields do not become significantly correlated over cosmological timescales. Smoother halos heat stellar orbits less efficiently, reducing the velocity dispersion relative to the single field case and thus weakening the observational constraints on the field mass. Analytically, we show that for N equal-mass fields with mass m the ULDM contribution to the stellar velocity dispersion scales as 1/(Nm 3 ). Finally, lighter fields heat the most efficiently and if the smallest mass mL is significantly below the other field masses the dispersion scales as $1/(N^{2}m^3_L$).

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Stability and Characteristics of Lower-hybrid Drift Waves: Dependence on Electron Beta and Cross-field Relative Drift

Lower-hybrid drift waves (LHDWs) are frequently observed microinstabilities in both space and laboratory plasmas. Despite decades of study, the relationship between electrostatic (ES-LHDW) and electromagnetic (EM-LHDW) variants and the plasma parameters controlling their stability remains unclear. Here, we systematically examine LHDW behavior by solving the local linear dispersion relation over a wide range of plasma and field conditions. Our results demonstrate that ES-LHDWs and EM-LHDWs are not distinct modes but rather two different regimes of the same drift wave whose character evolves smoothly with electron beta (β e ) and the cross-field electron drift velocity relative to ions, normalized to the ion sound speed (u 0x /C s ). The nature of the waves changes from electrostatic to electromagnetic when β e increases. Growth rates increase with u 0x /C s but decrease with β e , while the most unstable wavelength remains nearly universal, with kρ e ∼ 0.8 (k is the magnitude of the wave vector and ρ e is the electron gyroradius). We further present quasi-linear estimates of nonlinear saturation properties, including energy partition among electric fields, magnetic fields, and particle kinetic responses. We show that ES-LHDWs reach higher electric-field saturation amplitudes, whereas EM-LHDWs generate strong magnetic perturbations and parallel electric fields that may enable efficient particle heating. Comparisons with the classical model reveal that retaining electromagnetic effects is essential for accurate predictions of frequency, growth rate, and the propagation angle. These findings provide a unified framework for understanding LHDWs across diverse collisionless plasma environments, including current sheets of magnetic reconnection, shear layers, collisionless shocks, and boundary regions.

Solar coronal waves↗
DOE OSTI2022

Modified gravitational wave propagation with higher modes and its degeneracies with lensing

Low-energy alternatives to General Relativity (GR) generically modify the phase of gravitational waves (GWs) during their propagation. As detector sensitivities increase, it becomes key to understand how these modifications affect the GW higher modes and to disentangle possible degeneracies with astrophysical phenomena. We apply a general formalism — the WKB approach — for solving analytically wave propagation in the spatial domain with a modified dispersion relation (MDR). We compare this WKB approach to applying a stationary phase approximation (SPA) in the temporal domain with time delays associated to the group or particle velocity. To this end, we extend the SPA to generic signals with higher modes, keeping careful track of reference phases and arrival times. We find that the WKB approach coincides with the SPA using the group velocity, in agreement with the principles of wave propagation. We then explore the degeneracies between a GW propagation with an MDR and a strongly-lensed GW in GR, since the latter can introduce a frequency-independent phase shift which is not degenerate with source parameters in the presence of higher modes. We find that for a particular MDR there is an exact degeneracy for wave propagation, unlike with the SPA for particle propagation. For the other cases, we search for the values of the MDR parameters that minimize the χ 2 and conclude that strongly-lensed GR GWs could be misinterpreted as GWs in modified gravity. As a result, future MDR constraints with higher mode GWs should include the possibility of frequency-independent phase shifts, allowing for the identification of modified gravity and strong lensing distortions at the same time.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Bent Laue crystal anatomy: new insights into focusing and energy-dispersion properties

X-ray Laue-type monochromators are common and essential optical components at many high-power X-ray facilities, e.g. synchrotron facilities. The X-ray optics of bent Laue crystals is a well developed area. An incident X-ray beam penetrating a bent Laue crystal will result in a diffracted beam with different angles and energies. There is a need for a way of organizing the rays that allows one to sort out the energy and spatial properties of the diffracted beam. The present work introduces a new approach for describing the general behaviour of bent Laue crystals from a ray-tracing point of view. This quasi-monochromatic beam approach provides an intuitive view of bent-crystal diffraction and leads to deeper understanding. It explains the energy and spatial properties of common and special cases of bent Laue optics, predicts phenomena that can improve energy-dispersion-related X-ray imaging techniques and provides a theoretical framework that makes ray-tracing simulation easier to realize.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Evaluation of Turbulence and Dispersion in Multiscale Atmospheric Simulations over Complex Urban Terrain during the Joint Urban 2003 Field Campaign

Abstract This paper evaluates the representation of turbulence and its effect on transport and dispersion within multiscale and microscale-only simulations in an urban environment. These simulations, run using the Weather Research and Forecasting Model with the addition of an immersed boundary method, predict transport and mixing during a controlled tracer release from the Joint Urban 2003 field campaign in Oklahoma City, Oklahoma. This work extends the results of a recent study through analysis of turbulence kinetic energy and turbulence spectra and their role in accurately simulating wind speed, direction, and tracer concentration. The significance and role of surface heat fluxes and use of the cell perturbation method in the numerical simulation setup are also examined. Our previous study detailed the model development necessary for our multiscale simulations, examined model skill at predicting wind speeds and tracer concentrations, and demonstrated that dynamic downscaling from mesoscale to microscale through a sequence of nested simulations can improve predictions of transport and dispersion relative to a microscale-only simulation forced by idealized meteorology. Here, predictions are compared with observations to assess qualitative agreement and statistical model skill at predicting wind speed, wind direction, tracer concentration, and turbulent kinetic energy at locations throughout the city. We also investigate the scale distribution of turbulence and the associated impact on model skill, particularly for predictions of transport and dispersion. Our results show that downscaled large-scale turbulence, which is unique to the multiscale simulations, significantly improves predictions of tracer concentrations in this complex urban environment. Significance Statement Simulations of atmospheric transport and mixing in urban environments have many applications, including pollution modeling for urban planning or informing emergency response following a hazardous release. These applications include phenomena with spatial scales spanning from millimeters to kilometers. Most simulations resolve flow only within the urban area of interest, omitting larger scales of turbulence and regional influences. This study examines a method that resolves both the small and large-scale flow features. We evaluate simulation accuracy by comparing predictions with observations from an experiment involving the release of a tracer gas in Oklahoma City, Oklahoma, with emphasis on correctly modeling turbulent fluctuations. Our results demonstrate the importance of resolving large-scale flow features when predicting transport and dispersion in urban environments.

42 ENGINEERING↗
DOE OSTI2023

Dirichlet-type absorbing boundary conditions for peridynamic scalar waves in two-dimensional viscous media

Construction of absorbing boundary conditions (ABCs) for nonlocal models is generally challenging, primarily due to the fact that nonlocal operators are commonly associated with volume constrained boundary conditions. Moreover, application of Fourier and Laplace transforms, which are essential for the majority of available methods for ABCs, to nonlocal models is complicated. In this paper, we propose a simple method to construct accurate ABCs for peridynamic scalar wave-type problems in viscous media. The proposed ABCs are constructed in the time and space domains and are of Dirichlet type. Consequently, their implementation is relatively simple, since no derivatives of the wave field are required. The proposed ABCs are derived at the continuum level, from a semi-analytical solution of the exterior domain using harmonic exponential basis functions in space and time (plane-wave modes). The numerical implementation is done using a meshfree collocation approach employed within a boundary layer adjacent to the interior domain boundary. The modes satisfy the peridynamic numerical dispersion relation, resulting in a compatible solution of the interior region (near-field) with that of the exterior region (far-field). The accuracy and stability of the proposed ABCs are demonstrated with several numerical examples in two-dimensional unbounded domains.

42 ENGINEERING↗
DOE OSTI2024

Nonlinear wavepacket dynamics in proximity to a stationary inflection point

A stationary inflection point (SIP) in the Bloch dispersion relation of a periodic waveguide is an exceptional point degeneracy where three Bloch eigenmodes coalesce, forming the so-called frozen mode with a divergent amplitude and vanishing group velocity of its propagating component. For this work, we have developed a theoretical framework to study the time evolution of wavepackets centered at an SIP. Analysis of the evolution of the statistical moment distribution of linear pulses shows a strong deviation from the conventional ballistic wavepacket dynamics in dispersive media. The presence of nonlinear interactions dramatically changes the situation, resulting in a mostly ballistic propagation of nonlinear wavepackets with the speed and even the direction of propagation essentially dependent on the wavepacket amplitude. Such a behavior is unique to nonlinear wavepackets centered at an SIP and can be used for the realization of a novel family of beam power routers for classical waves.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Viscous dissipation in a gas of one-dimensional fermions with generic dispersion

A well-known feature of the classical monoatomic gas is that its bulk viscosity is strongly suppressed because the single-particle dispersion is quadratic. On the other hand, in condensed matter systems the effective single -particle dispersion is altered by lattice effects and interactions. In this work, we study the bulk viscosity of one-dimensional Fermi gases with generic energy-momentum dispersion relations. As an application, viscous dissipation arising from lattice effects is analyzed for the tight-binding model. In addition, we investigate how weak interactions affect the bulk viscosity. Finally, we discuss viscous dissipation in the regime in which the Fermi gas is not fully equilibrated, as can occur when the system is driven at frequencies that exceed the rate of fermion backscattering. In this case, the Fermi gas is described by three bulk viscosities, which we obtain for a generic single-particle dispersion.

1-dimensional systems↗
DOE OSTI2022

Theory of gradient drift instabilities in low-temperature, partially magnetised plasmas

A fluid dispersion theory in partially magnetised plasmas is analysed to examine the conditions under which large-wavelength modes develop in Penning-type configurations, that is, where an electric field is imposed perpendicular to a homogeneous magnetic field. The fluid dispersion relation assuming a slab geometry shows that two types of low-frequency, gradient drift instabilities occur in the direction of the E×B and diamagnetic drifts. One type of instability, observed when the equilibrium electric field and plasma density gradient are in the same direction, is similar to the classic modified Simon–Hoh instability. A second instability is found for conditions in which (i) the diamagnetic drift is in the direction opposite to the E×B drift and (ii) the magnitude of the diamagnetic drift is sufficiently larger than the electron thermal speed. The present fluid dispersion theory suggests that the rotating spokes driven by such fluid instabilities propagate in the same direction as the diamagnetic drift, which can be in the same direction as or opposite to the E×B drift, depending on the plasma conditions. This finding may account for the observation, in some plasma devices, of the rotation of large-scale structures in both the E×B and -E×B directions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗