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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Mechanism of H 2 plasma-enabled reduction of hematite thin films

Hydrogen plasma is gaining significant interest as a promising pathway for direct iron ore reduction and for lowering process temperatures but the reduction mechanism remains poorly understood. In this work we analyzed the plasma and thermal reduction of thin-film hematite (Fe 2 O 3 ) at temperatures below 340 °C using X-ray diffraction and scanning electron microscopy with energy-dispersive X-ray spectroscopy. Plasma reduced the incubation period by an order of magnitude and increased the reduction rate by a factor of 2.6 compared to thermal reduction. Plasma-produced H-atoms facilitate the formation of numerous iron nucleation sites, bypassing the energetically unfavorable dissociative adsorption of H 2 on iron oxide. These iron nuclei can autocatalyze the reduction of the surrounding hematite via a hydrogen spillover mechanism. Our results demonstrate that plasma-derived H-atoms primarily impact the initial nucleation-limited stage. These new insights provide a mechanistic framework that can aid the implementation and optimization of hydrogen plasma-assisted iron oxide reduction at reduced temperatures.

08 HYDROGEN↗
DOE OSTI2022

Three-dimensional mapping of mineral in intact shark centra with energy dispersive x-ray diffraction

The centra of shark vertebrae consist of cartilage mineralized by a bioapatite similar to bone's carbonated hy-droxyapatite, and, without a repair mechanism analogous to remodeling in bone, these structures still survive millions of cycles of high-strain loading. The main structures of the centrum are an hourglass-shaped double cone and the intermedialia which supports the cones. Little is known about the nanostructure of shark centra, spe-cifically the relationship between bioapatite and cartilage fibers, and this study uses energy dispersive diffraction (EDD) with polychromatic synchrotron x-radiation to study the spatial organization of the mineral phase and its crystallographic texture. The unique energy-sensitive detector array at beamline 6-BM-B, the Advanced Photon Source, enables EDD to quantify the texture within each sampling volume with one exposure while constructing 3D maps via specimen translation across the sampling volume. Herein this study maps a centrum from two shark orders, a carcharhiniform and a lamniform, with different intermedialia structures. In the blue shark (Prionace glauca, Carcharhiniformes), the bioapatite's c-axes are oriented laterally within the centrum's cone walls but axially within the wide wedges of the intermedialia; the former is interpreted to resist lateral deformation, the latter to support axial loads. In the shortfin mako (Isurus oxyrinchus, Lamniformes), there is some tendency for c-axis variation with position, but the situation is unclear because one dimension of the sampling volume is consid-erably larger than the thickness and spacing of the intermedialia's radially-oriented lamellae. Because elastic modulus in collagen plus bioapatite mineralized tissues varies significantly with both volume fraction of bio-apatite and crystallographic texture, the present 3D EDD-derived maps should inform future 3D numerical models of shark centra under applied load.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Quasi-periodic travelling gravity–capillary waves

We present a numerical study of spatially quasi-periodic travelling waves on the surface of an ideal fluid of infinite depth. This is a generalization of the classic Wilton ripple problem to the case when the ratio of wavenumbers satisfying the dispersion relation is irrational. We propose a conformal mapping formulation of the water wave equations that employs a quasi-periodic variant of the Hilbert transform to compute the normal velocity of the fluid from its velocity potential on the free surface. We develop a Fourier pseudo-spectral discretization of the travelling water wave equations in which one-dimensional quasi-periodic functions are represented by two-dimensional periodic functions on the torus. This leads to an overdetermined nonlinear least-squares problem that we solve using a variant of the Levenberg–Marquardt method. We investigate various properties of quasi-periodic travelling waves, including Fourier resonances, time evolution in conformal space on the torus, asymmetric wave crests, capillary wave patterns that change from one gravity wave trough to the next without repeating and the dependence of wave speed and surface tension on the amplitude parameters that describe a two-parameter family of waves.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Three-dimensional atomic structure and local chemical order of medium- and high-entropy nanoalloys

Medium- and high-entropy alloys (M/HEAs) mix several principal elements with near-equiatomic composition and represent a model-shift strategy for designing previously unknown materials in metallurgy, catalysis and other fields. One of the core hypotheses of M/HEAs is lattice distortion, which has been investigated by different numerical and experimental techniques. However, determining the three-dimensional (3D) lattice distortion in M/HEAs remains a challenge. Moreover, the presumed random elemental mixing in M/HEAs has been questioned by X-ray and neutron studies, atomistic simulations, energy dispersive spectroscopy and electron diffraction, which suggest the existence of local chemical order in M/HEAs. However, direct experimental observation of the 3D local chemical order has been difficult because energy dispersive spectroscopy integrates the composition of atomic columns along the zone axes and diffuse electron reflections may originate from planar defects instead of local chemical order. Here, in this work, we determine the 3D atomic positions of M/HEA nanoparticles using atomic electron tomography and quantitatively characterize the local lattice distortion, strain tensor, twin boundaries, dislocation cores and chemical short-range order (CSRO). We find that the high-entropy alloys have larger local lattice distortion and more heterogeneous strain than the medium-entropy alloys and that strain is correlated to CSRO. We also observe CSRO-mediated twinning in the medium-entropy alloys, that is, twinning occurs in energetically unfavoured CSRO regions but not in energetically favoured CSRO ones, which represents, to our knowledge, the first experimental observation of correlating local chemical order with structural defects in any material. We expect that this work will not only expand our fundamental understanding of this important class of materials but also provide the foundation for tailoring M/HEA properties through engineering lattice distortion and local chemical order.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Shallow Crustal Structure and Site Response at Fort Greely, Alaska

Fort Greely, Alaska is located in proximity to numerous faults capable of generating large (MW > 6.5) and damaging earthquakes. However, little knowledge exists about the noise properties, shallow crustal structure and site response of the area. We address this by examining ambient noise conditions, Rayleigh wave group wavespeed dispersion, and site response measurements using new and existing seismic data from December 2020 – April 2021 for four sites along a north-south line—stations GREN, GRES, PS09, K24K operated by the Alaska Earthquake Center. We find that ambient noise levels are relatively high on base compared to other places in Central Alaska, especially above 1 Hz, while noise levels in the microseismic band are typical of the region. Short-period (0.5 – 2Hz) Rayleigh wave group wavespeeds appear to decrease with increasing sediment thickness and increase when going southward from GREN to K24K. The predominant site frequency related to sediment thickness is lowest beneath PS09 and highest beneath K24K. We use the group velocities and site response values to propose a rough north-south layered model for the region.

58 GEOSCIENCES↗
DOE OSTI2026

Excitation and tunneling spectra of a fractional quantum Hall system in the thin-cylinder limit

Here, the excitations of fractional quantum Hall effect (FQHE) states have been largely inaccessible to experimental probes until recently. New electron scanning tunneling microscopy (STM) results from Hu et al. [Nat. Phys. 21, 716 (2025)] show promise in detecting and identifying these excited states via the local density of states (LDOS) spectrum. On a torus, there exists a mapping from the lowest Landau level states to a 1D lattice with a Hamiltonian that features dipole moment conservation. In this work, we apply perturbation theory starting from the thin-cylinder limit (𝐿 𝑥 → ∞, 𝐿 𝑦 < 𝑙 𝐵 for torus dimensions 𝐿 𝑥 and 𝐿 𝑦 and magnetic length 𝑙 𝐵 ) to obtain an analytical approach to the low-lying neutral and charged excitations of the 𝜈 = 1/3 FQHE state. Notably, in the thin cylinder, we can systematically enumerate all the low-lying excitations by the patterns of “dipoles” formed by the electron occupation pattern on the 1D lattice. We find that the thin-cylinder limit predicts a significant dispersion of the low-lying neutral excitations but sharpness of the LDOS spectra, which measure charged excitations. We also discuss connections between our work and several different approaches to the FQHE STM spectra, including those using the composite fermion theory. Numerical exact diagonalization beyond the thin-cylinder limit suggests that the energies of charged excitations remain largely confined to a narrow range of energies, which in experiments might appear as a single peak.

Adhidewata, Jyesta M. [University of California, B↗
DOE OSTI2024

Design of experiments to spectroscopically characterize radiation flow in stochastic media

Precise characterization of experimental radiation flow is required to validate the high energy density physics models, numerical methods, and codes that are used to simulate radiation-hydrodynamics phenomena such as thermal radiation transport in stochastic media. The Cassio code is used to simulate thermal radiation flow through inhomogeneous, stochastic-media-foam configurations containing optically thick clumps dispersed within an optically thin background aerogel. Cassio can model small inhomogeneous problems directly, but most problems require approximations to meet computer limitations on run-times and memory usage. Various examples of these approximations are methods that produce, in one calculation, an ensemble-averaged solution and associated standard deviation; reduced spatial dimensionality with approximate geometries; and full material homogenization with no geometric detail. Cassio simulations are used to design experiments at the OMEGA-60 Laser Facility that can measure the radiation flow using the spatially resolved COAX absorption spectroscopy diagnostic. The experimental platforms flow radiation through foam targets ranging from a background-only aerogel, to a single configuration of a specified stochastic medium, to a fully homogenized foam of the background and clump materials. Under constant total clump mass, larger clumps (here, larger than 10 μm diameter) will mix more slowly with the background such that the bulk radiation flow is faster than it would be in a fully homogenized material. The COAX platform can be used to infer temperature and density profiles in both the background material and clumps, simultaneously, and therefore to differentiate radiation flow in a range of stochastic and homogeneous media.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Characterizing the magnetic fields of nearby molecular clouds using submillimeter polarization observations

Of all the factors that influence star formation, magnetic fields are perhaps the least well understood. The goal of this paper is to characterize the 3D magnetic field properties of nearby molecular clouds through various methods of statistically analysing maps of polarized dust emission. Our study focuses on nine clouds, with data taken from the Planck Sky Survey as well as data from the Balloon-borne Large Aperture Submillimeter Telescope for Polarimetry observations of Vela C. We compare the distributions of polarization fraction (p), dispersion in polarization angles ($\mathcal {S}$), and hydrogen column density (N H ) for each of our targeted clouds. To broaden the scope of our analysis, we compare the distributions of our clouds’ polarization observables with measurements from synthetic polarization maps generated from numerical simulations. We also use the distribution of polarization fraction measurements to estimate the inclination angle of each cloud’s cloud-scale magnetic field. We obtain a range of inclination angles associated with our clouds, varying from 16° to 69°. We establish inverse correlations between p and both $\mathcal {S}$ and NH in almost every cloud, but we are unable to establish a statistically robust $\mathcal {S}$ versus NH trend. By comparing the results of these different statistical analysis techniques, we are able to propose a more comprehensive view of each cloud’s 3D magnetic field properties. These detailed cloud analyses will be useful in the continued studies of cloud-scale magnetic fields and the ways in which they affect star formation within these molecular clouds.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Emergent pair density wave order across a Lifshitz transition

We numerically investigate the telltale signs of pair-density-wave order (PDW) in the Kondo-Heisenberg chain by focusing on the momentum resolved spectrum in different parameter regimes. Density matrix renormalization group calculations reveal that this phase is characterized by a dispersion with two minima and four Fermi points, indicating the emergence of an effective next-nearest-neighbor hopping that arises as a second-order effect to avoid magnetic frustration. The pairs appear in the spectrum as in-gap bound states with weight concentrated in the hole pockets. The low-energy physics can be understood by means of a generalized t−J model with next-nearest-neighbor hopping. Our results offer a guide for searching for experimental signatures, and for other models that can realize PDW physics.

Yang, Luhang [University of Tennessee Knoxville]↗
DOE OSTI2020

Strings in Ramond-Ramond backgrounds from the Neveu-Schwarz-Ramond formalism

We treat RR flux backgrounds of type II string theory in the framework of closed superstring field theory based on the NSR formalism, focusing on two examples: (1) the pp-wave background supported by 5-form flux, and (2) AdS 3 × S 3 × M 4 supported by mixed 3-form fluxes. In both cases, we analyze the classical string field solution perturbatively, and compute the correction to the dispersion relation of string states to quadratic order in the RR flux. In the first example, our result is in a delicate way consistent with that obtained from lightcone quantization of the Green-Schwarz string. In the second example, we will obtain numerically the mass corrections to pulsating type IIB strings in AdS 3 × S 3 × M 4 . Our results, valid at finite AdS radius, agree with previously known answers in the semiclassical limit and in the BMN limit respectively.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Suppression of one-dimensional weak localization by band asymmetry

Herein we investigate disorder-induced localization in metals that break time-reversal and inversion symmetries through their energy dispersion, ε k ≠ε– k , but lack Berry phases. In the perturbative regime of disorder, we show that weak localization is suppressed due to a mismatch of the Fermi velocities of left and right movers. To substantiate this analytical result, we perform quench numerics on chains shorter than the Anderson localization length ξ—the latter computed and verified to be finite using the recursive Green's function method—and find a sharp rise in the saturation value of the participation ratio due to band asymmetry, indicating a tendency to delocalize. Interestingly, for weak disorder strength η, we see a better fit to the scaling behavior ξ∝1/η 2 for asymmetric bands than conventional symmetric ones.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Multimode description of self-mode locking in a single-section quantum-dot laser

This paper describes a theory for mode locking and frequency comb generation by four-wave mixing in a semiconductor quantum-dot active medium. The derivation uses a multimode semiclassical laser theory that accounts for fast carrier collisions within an inhomogeneous distribution of quantum dots. Numerical simulations are presented to illustrate the role of active medium nonlinearities in mode competition, gain saturation, carrier-induced refractive index and creation of combination tones that lead to locking of beat frequencies among lasing modes in the presence of cavity material dispersion.

Chow, Weng W.↗
DOE OSTI2022

Shock-driven dispersal of a corrugated finite-thickness particle layer

A research area emerging in the multiphase flow community is the study of shock-driven multiphase instability (SDMI), a gas–particle analog of the traditional fluid-fluid Richtmyer–Meshkov instability (RMI). In this work, we study the interaction of planar air shocks with corrugated glass particle curtains through the use of numerical simulations with an Eulerian–Lagrangian approach. One objective of this study is to compare the simulated particle curtains to a comparable set of shock tube experiments performed to analyze traditional RMI of a gas curtain. The simulations are set to match the experimental shock Mach numbers and perturbation wavelengths (3.6 and 7.2 mm) while also matching the Atwood number of the experiments to the multiphase Atwood number of the simulations. Varying particle diameters are tested in the simulations to explore the impact of particle diameter on the evolution of the particle curtain. This simulation setup allows for a one-to-one comparison between RMI and SDMI under comparable conditions while also allowing for a separate study into the validity of the use of the multiphase Atwood number to compare the single-phase and multiphase instabilities. In particular, we show that the comparison depends on the diameter of the particles (thus, dependent on the Stokes number of the flow). A second objective of this study is to analyze the effect of the initial particle volume fraction on the evolution of the curtain and the behavior of the instability. This is done through analyzing the effect of the multiphase terms of the vorticity evolution equation on the vorticity deposition in SDMI. Also discussed is the effect of the particle diameter on the multiphase generation terms as well as in the baroclinic vorticity generation term in SDMI as the shock passes over the curtain.

42 ENGINEERING↗
DOE OSTI2023

A conservative Galerkin solver for the quasilinear diffusion model in magnetized plasmas

We propose a conservative Galerkin scheme for the quasilinear model in three-dimensional momentum space and three-dimensional spectral space, with cylindrical symmetry. We construct an unconditionally conservative weak form and use a discretization that preserves conservation properties independent of the wave emission probability. The discrete operators, combined with a consistent quadrature rule, preserve all the conservation laws rigorously. The proposed scheme is quite general: it works for both relativistic and non-relativistic systems, for both magnetized and unmagnetized plasmas, and even for problems with time-dependent dispersion relations. We represent the particle distribution by continuous basis functions and use discontinuous basis functions for the wave spectral energy density, which enables the application of a positivity-preserving technique. We adopt the marching simplex algorithm, designed initially for computer graphics, for numerical integration on the resonance manifold. Furthermore, the numerical examples with a bump-on-tail initial configuration show how the unstable waves produce strong momentum space diffusion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Complete quasilinear model for the acceleration-driven lower hybrid drift instability and a computational assessment of its validity

A complete quasilinear model is derived for the electrostatic acceleration-driven lower hybrid drift instability in a uniform two-species low-beta plasma in which current is perpendicular to the background magnetic field. The model consists of coupled nonlinear velocity space diffusion equations for the volume-averaged ion and electron distribution functions. Each species' diffusion coefficient depends on a time-evolving spectral density of the electric-field energy per unit volume and a time-evolving dispersion relation. The dispersion relation is expressed analytically in integral form without the use of asymptotic limits and applies to arbitrary distribution functions, so long as they can be expressed as a function of one velocity coordinate, e.g., f⁡(vy) or f⁡(v⊥). The quasilinear model conserves energy and is complete in that it fully describes the evolution of the distribution functions, including resonant and nonresonant particle-wave interactions, while accounting for distribution-function-dependent mixed-complex frequencies. Further, the quasilinear diffusion model is solved numerically and self-consistently using a Crank-Nicolson temporal discretization and a second-order finite-volume velocity-space discretization. Numerical solutions are compared to nonlinear fourth-order accurate continuum kinetic Vlasov-Poisson simulations. Evolution of electric-field energy, growth rates, distribution functions, and diffusion coefficients are shown to be in agreement with Vlasov simulations. The quasilinear model is shown to predict anomalous transport terms, like resistivity and heating, to within a factor of order unity. Discrepancies between the quasilinear model and Vlasov simulations are assessed and attributed primarily to lack of damping in the quasilinear description and to the use of unperturbed-orbit susceptibilities in the linear theory dispersion relation. The results illuminate the predictive accuracy of the quasilinear model, place approximate bounds on its validity, and provide much needed vetting of quasilinear theory's ability to predict the nonlinear state of a microturbulent plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

From Intermolecular Interaction Energies and Observable Shifts to Component Contributions and Back Again: A Tale of Variational Energy Decomposition Analysis

Quantum chemistry in the form of density functional theory (DFT) calculations is a powerful numerical experiment for predicting intermolecular interaction energies. However, no chemical insight is gained in this way beyond predictions of observables. Energy decomposition analysis (EDA) can quantitatively bridge this gap by providing values for the chemical drivers of the interactions, such as permanent electrostatics, Pauli repulsion, dispersion, and charge transfer. These energetic contributions are identified by performing DFT calculations with constraints that disable components of the interaction. This review describes the second-generation version of the absolutely localized molecular orbital EDA (ALMO-EDA-II). The effects of different physical contributions on changes in observables such as structure and vibrational frequencies upon complex formation are characterized via the adiabatic EDA. Example applications include red- versus blue-shifting hydrogen bonds; the bonding and frequency shifts of CO, N 2 , and BF bound to a [Ru(II)(NH 3 ) 5 ] 2 + moiety; and the nature of the strongly bound complexes between pyridine and the benzene and naphthalene radical cations. Additionally, the use of ALMO-EDA-II to benchmark and guide the development of advanced force fields for molecular simulation is illustrated with the recent, very promising, MB-UCB potential.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Interaction of Inner Heliosheath Ions with Electromagnetic Ion Cyclotron Waves

Abstract Inner heliosheath (IHS) ions are expected to be exposed to various waves, shocks, and turbulence, which can affect ion distributions and thus their charge-exchange rates with interstellar neutral atoms. This work addresses the potential significance of electromagnetic ion cyclotron (EMIC) waves under expected IHS conditions. From a kinetic dispersion relation, we find the possibility of frequent triggering of EMIC instability in the IHS. The threshold anisotropy of proton temperatures required for the instability is small, ( T ⊥ – T ∣∣ )/ T ∣∣ ≈ 0.1 or less, mainly due to high plasma β (ratio of the plasma to magnetic pressures). Numerical calculations on the scattering of ions (protons, He + , and He 2+ with energy of 0.01–50 keV) based on two models for EMIC waves with a moderate intensity indicate significant scattering in the pitch angle (mostly a few tens of degrees) and energy (mostly a few tens of percent) although details depend on the energy and pitch angle of each species and adopted EMIC wave models. This occurs on a short timescale (<100 times the gyro-period of each ion species). Resonant scattering in a few to a few tens of keV (corresponding to the expected pickup ion energy in the IHS) is easily expected unless the wavenumber is too large. The scattering effect is distinguished among different species such that for lower gyrofrequency ions (He + versus He 2+ and He 2+ versus protons), the main scattering effect moves toward a lower energy domain. All these results imply continuous disturbance of ion distributions by possibly prevailing EMIC waves in the IHS.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Multiscale modeling of microwave-heated multiphase systems

Modeling of microwave heating of multiphase systems requires resolving the electromagnetic field from the single particle to the entire cavity. We introduce a multiscale methodology for computationally affordable simulations of microwave heating of multiphase systems consisting of one phase dispersed in a continuum phase. The methodology homogenizes the original multiphase system, by taking advantage of the large separation of length scales, to calculate the effective permittivity, effective thermal conductivity, and volumetric power absorbed in each phase. The methodology is rigorously assessed against particle-resolved detailed numerical simulations of model systems. We demonstrate a significant reduction in computational cost while retaining the accuracy of detailed simulations. Lastly, the present methodology enables high throughput exploration of novel multiphase system designs utilizing microwave heating.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗