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 253 records · Page 14

Shear thinning rheology in a concentrated suspension of fibers: the role of attractive interactions

Fiber-reinforced composites are ubiquitously encountered in the engineering, automobile, and aerospace industries. Fabrication of these composites requires the mixing of fibers dispersed in a liquid, and the final structure is affected by fiber properties, interactions, and flow fields. A good understanding of the rheology of fiber suspensions can aid in the design and optimization of the fabrication processes. To this end, we simulate the simple shear flow of fiber suspensions by accounting for short range lubrication, van der Waals attractive, electrostatic repulsion interactions, as well as friction between fibers. Direct numerical simulations are performed using the Immersed Boundary Method where the fibers are modeled as flexible slender bodies governed by the Euler-Bernoulli beam theory and the fluid flow is resolved using the Navier-Stokes equations. The simulation results for the suspension viscosity and yield stress are consistent with the experimental data from the literature for polyamide fiber suspensions. The shear-thinning behavior becomes stronger as we increase the magnitude of attractive interactions. Furthermore, we perform a parametric study varying fiber flexibility, volume fraction, aspect ratio, inertia and examine suspension viscosity and yield stress.

Khan, Monsurul↗

Observation of maxon-like ultrasound in elastic metabeam

We observe maxon-like dispersion of ultrasonic guided waves in elastic metamaterials consisting of a rectangular beam and an array of cylindrical resonators. The pillars act as asymmetric resonators that induce a strong modal hybridization. We experimentally observe the strongly localized maxon mode with zero group velocity. Our study also demonstrates a unique feature of the maxon with a down-shifting peak frequency in space. To reveal the fundamental mechanism, we conduct comprehensive numerical studies on all frieze group symmetries and key geometric parameters.

Materials Science↗

Scalable branching on dual decomposition of stochastic mixed-integer programming problems

In this work, we present a scalable branching method for the dual decomposition of stochastic mixed-integer programming. Our new branching method is based on the branching method proposed by Caroe and Schultz that creates branching disjunctions on first-stage variables only. We propose improvements to the process for creating branching disjunctions, including (1) branching on the optimal solutions of the Dantzig-Wolfe reformulation of the restricted master problem and (2) using a more comprehensive (yet simple) measure for the dispersions associated with subproblem solution infeasibility. We prove that the proposed branching process leads to an algorithm that terminates finitely, and we provide conditions under which globally optimal solutions can be identified after termination. We have implemented our new branching method, as well as the Caroe-Schultz method and a branch-and-price method, in the open-source software package DSP. Using SIPLIB test instances, we present extensive numerical results to demonstrate that the proposed branching method significantly reduces the number of node subproblems and solution times.

97 MATHEMATICS AND COMPUTING↗

Dynamically controlled random lasing with colloidal titanium carbide MXene

Control of lasing properties through tailorable and dynamically tunable materials and reconfigurable compositions can augment the performance of random lasers for a wide range of applications. Here, a colloid of randomly dispersed weakly scattering single-layer titanium carbide (Ti 3 C 2 T x ) MXene flakes embedded within rhodamine 101 gain medium is experimentally shown to provide feedback for random lasing. Additionally, in contrast to previously reported random laser systems where the optical properties of scatterers are static, the relative permittivity of Ti 3 C 2 T x MXene flakes can be varied under optical pumping due to the saturable absorption properties. Numerical simulations indicate that the observed nonlinear response of Ti 3 C 2 T x MXene flakes enables dynamically tunable random lasing. Thus, pumping the Ti 3 C 2 T x MXene flakes with a second optical source decreases the gain threshold required to obtain random lasing. Also, using numerical simulations, it is shown that the control over the intensity of the second pump enables tuning the field distribution of the random lasing modes. Considering the diversity of the MXenes family, the proposed MXene colloidal metamaterial design opens up a new avenue to advanced control of lasing properties for photonic applications.

36 MATERIALS SCIENCE↗

Planck 2018 results XII. Galactic astrophysics using polarized dust emission

Observations of the submillimetre emission from Galactic dust, in both total intensity I and polarization, have received tremendous interest thanks to the Planck full-sky maps. In this paper we make use of such full-sky maps of dust polarized emission produced from the third public release of Planck data. As the basis for expanding on astrophysical studies of the polarized thermal emission from Galactic dust, we present full-sky maps of the dust polarization fraction p, polarization angle ψ, and dispersion function of polarization angles S . The joint distribution (one-point statistics) of p and NH confirms that the mean and maximum polarization fractions decrease with increasing NH. The uncertainty on the maximum observed polarization fraction, pmax = 22.0 -1.4 +3.5% at 353 GHz and 80' resolution, is dominated by the uncertainty on the Galactic emission zero level in total intensity, in particular towards diffuse lines of sight at high Galactic latitudes. Furthermore, the inverse behaviour between p and S found earlier is seen to be present at high latitudes. This follows the S ∝ p -1 relationship expected from models of the polarized sky (including numerical simulations of magnetohydrodynamical turbulence) that include effects from only the topology of the turbulent magnetic field, but otherwise have uniform alignment and dust properties. Thus, the statistical properties of p, ψ, and S for the most part reflect the structure of the Galactic magnetic field. Nevertheless, we search for potential signatures of varying grain alignment and dust properties. First, we analyse the product map S × p, looking for residual trends. While the polarization fraction p decreases by a factor of 3-4 between NH = 10 20 cm -2 and NH = 2 x 10 22 cm -2 , out of the Galactic plane, this product S × p only decreases by about 25%. Because S is independent of the grain alignment efficiency, this demonstrates that the systematic decrease in p with N H is determined mostly by the magnetic-field structure and not by a drop in grain alignment. This systematic trend is observed both in the diffuse interstellar medium (ISM) and in molecular clouds of the Gould Belt. Second, we look for a dependence of polarization properties on the dust temperature, as we would expect from the radiative alignment torque (RAT) theory. We find no systematic trend of S × p with the dust temperature T d , whether in the diffuse ISM or in the molecular clouds of the Gould Belt. In the diffuse ISM, lines of sight with high polarization fraction p and low polarization angle dispersion S tend, on the contrary, to have colder dust than lines of sight with low p and high S . We also compare the Planck thermal dust polarization with starlight polarization data in the visible at high Galactic latitudes. The agreement in polarization angles is remarkable, and is consistent with what we expect from the noise and the observed dispersion of polarization angles in the visible on the scale of the Planck beam. The two polarization emission-to-extinction ratios, R P/p and R S/V , which primarily characterize dust optical properties, have only a weak dependence on the column density, and converge towards the values previously determined for translucent lines of sight. We also determine an upper limit for the polarization fraction in extinction, p V /E(B - V), of 13% at high Galactic latitude, compatible with the polarization fraction p ≈ 20% observed at 353 GHz. Taken together, these results provide strong constraints for models of Galactic dust in diffuse gas.

79 ASTRONOMY AND ASTROPHYSICS↗

Multispectral Quantitative Phase Imaging Using a Diffractive Optical Network

As a label‐free imaging technique, quantitative phase imaging (QPI) provides optical path length information of transparent specimens for various applications in biology, materials science, and engineering. Multispectral QPI measures quantitative phase information across multiple spectral bands, permitting the examination of wavelength‐specific phase and dispersion characteristics of samples. Herein, the design of a diffractive processor is presented that can all‐optically perform multispectral quantitative phase imaging of transparent phase‐only objects within a snapshot. The design utilizes spatially engineered diffractive layers, optimized through deep learning, to encode the phase profile of the input object at a predetermined set of wavelengths into spatial intensity variations at the output plane, allowing multispectral QPI using a monochrome focal plane array. Through numerical simulations, diffractive multispectral processors are demonstrated to simultaneously perform quantitative phase imaging at 9 and 16 target spectral bands in the visible spectrum. The generalization of these diffractive processor designs is validated through numerical tests on unseen objects, including thin Pap smear images. Due to its all‐optical processing capability using passive dielectric diffractive materials, this diffractive multispectral QPI processor offers a compact and power‐efficient solution for high‐throughput quantitative phase microscopy and spectroscopy.

36 MATERIALS SCIENCE↗

Evaluation of several relevant fractionation processes as possible explanation for radioxenon isotopic activity ratios in samples taken near underground nuclear explosions in shafts and tunnels

Gas samples taken from two historic underground nuclear tests done in 1989 at the Nevada National Security Site (NNSS), formerly the Nevada Test Site (NTS), were examined to determine how xenon isotopes fractionate because of early-time cavity processes, transport through the rock, or dispersal through tunnels. Xenon isotopes are currently being used to distinguish civilian sources of xenon in the atmosphere from sources associated with underground nuclear explosions (UNEs). The two nuclear tests included (1) BARNWELL, a test conducted in a vertical shaft approximately 600 m below ground surface at Pahute Mesa, and (2) DISKO ELM, a horizontal line-of-sight test done in P-tunnel approximately 261 m below the surface of Aqueduct Mesa. Numerical flow and transport models developed for the two sites had mixed success when attempting to match the observed xenon isotope ratios. At the BARNWELL site, the simulated xenon isotope ratios were consistent with measurements from the chimney and ground surface and appeared to have been affected primarily by fractionation during subsurface transport. At the DISKO ELM site, samples taken from two elevations in the chimney failed to show the degree of fractionation predicted by the models during transport and did not show evidence for significant fractionation due to early-time condensation of refractory xenon-precursor radionuclides into the melt glass. Gas samples taken from the adjacent tunnels in the days following the test showed mixed evidence for early-time separation of xenon isotopes from their iodine precursors.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Toward more mechanistic representations of biogeochemical processes in river networks: Implementation and demonstration of a multiscale model

A recently introduced multiscale model for representing the effects of small-scale hyporheic-zone biogeochemical processes is extended from the reach scale to river network scales. The model uses advection-dispersion-reaction equations for the channel network and one-dimensional advection-reaction subgrid models for the hyporheic zone. We summarize implementation in the integrated surface/subsurface hydrology modeling system Amanzi-ATS. The extension uses topologically defined meshes to represent stream/river networks and associates a hyporheic subgrid model with each channel grid cell. Biogeochemical reaction modeling capability residing in community software is accessed through an application programming interface. The implementation is verified against independent numerical solutions on a single reach. Mesh convergence studies show that commonly used semi-distributed representations can introduce significant spatial discretization error. Denitrification of farm runoff in a subbasin of the Portage River Basin in Ohio, USA is used to demonstrate the general-purpose reactive transport capability.

54 ENVIRONMENTAL SCIENCES↗

Direct observation of the formation and stabilization of metallic nanoparticles on carbon supports

Direct formation of ultra-small nanoparticles on carbon supports by rapid high temperature synthesis method offers new opportunities for scalable nanomanufacturing and the synthesis of stable multi-elemental nanoparticles. However, the underlying mechanisms affecting the dispersion and stability of nanoparticles on the supports during high temperature processing remain enigmatic. In this work, we report the observation of metallic nanoparticles formation and stabilization on carbon supports through in situ Joule heating method. We find that the formation of metallic nanoparticles is associated with the simultaneous phase transition of amorphous carbon to a highly defective turbostratic graphite (T-graphite). Molecular dynamic (MD) simulations suggest that the defective T-graphite provide numerous nucleation sites for the nanoparticles to form. Furthermore, the nanoparticles partially intercalate and take root on edge planes, leading to high binding energy on support. This interaction between nanoparticles and T-graphite substrate strengthens the anchoring and provides excellent thermal stability to the nanoparticles. These findings provide mechanistic understanding of rapid high temperature synthesis of metal nanoparticles on carbon supports and the origin of their stability.

36 MATERIALS SCIENCE↗

Crosstalk-robust quantum control in multimode bosonic systems

High-coherence superconducting cavities offer a hardware-efficient platform for quantum information processing. To achieve universal operations of these bosonic modes, the requisite nonlinearity is realized by coupling them to a transmon ancilla. However, this configuration is susceptible to crosstalk errors in the dispersive regime, where the ancilla frequency is Stark shifted by the state of each coupled bosonic mode. This leads to a frequency mismatch of the ancilla drive, lowering the gate fidelities. To mitigate such coherent errors, we employ quantum optimal control to engineer ancilla pulses that are robust to the frequency shifts. These optimized pulses are subsequently integrated into a recently developed echoed conditional displacement protocol for executing single- and two-mode operations. Through numerical simulations, we examine two representative scenarios: the preparation of single-mode Fock states in the presence of spectator modes and the generation of two-mode entangled Bell-cat states. Our approach markedly suppresses crosstalk errors, outperforming conventional ancilla control methods by orders of magnitude. These results provide guidance for experimentally achieving high-fidelity multimode operations and pave the way for developing high-performance bosonic quantum information processors.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Crossing Symmetric Dispersion Relations in Quantum Field Theories

For 2-2 scattering in quantum field theories, the usual fixed t dispersion relation exhibits only two-channel symmetry. This Letter considers a crossing symmetric dispersion relation, reviving certain old ideas from the 1970s. Rather than the fixed t dispersion relation, this needs a dispersion relation in a different variable z, which is related to the Mandelstam invariants s, t, u via a parametric cubic relation making the crossing symmetry in the complex z plane a geometric rotation. The resulting dispersion is manifestly three-channel crossing symmetric. We give simple derivations of certain known positivity conditions for effective field theories, including the null constraints, which lead to two sided bounds and derive a general set of new nonperturbative inequalities. We show how these inequalities enable us to locate the first massive string state from a low energy expansion of the four dilaton amplitude in type II string theory. We also show how a generalized (numerical) Froissart bound, valid for all energies, is obtained from this approach.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

“Gray Zone” Simulations Using a Three-Dimensional Planetary Boundary Layer Parameterization in the Weather Research and Forecasting Model

Abstract Generating accurate weather forecasts of planetary boundary layer (PBL) properties is challenging in many geographical regions, oftentimes due to complex topography or horizontal variability in, for example, land characteristics. While recent advances in high-performance computing platforms have led to an increase in the spatial resolution of numerical weather prediction (NWP) models, the horizontal gridcell spacing (Δ x ) of many regional-scale NWP models currently fall within or are beginning to approach the gray zone (i.e., Δ x ≈ 100–1000 m). At these gridcell spacings, three-dimensional (3D) effects are important, as the most energetic turbulent eddies are neither fully parameterized (as in traditional mesoscale simulations) nor fully resolved [as in traditional large-eddy simulations (LES)]. In light of this modeling challenge, we have implemented a 3D PBL parameterization for high-resolution mesoscale simulations using the Weather Research and Forecasting Model. The PBL scheme, which is based on the algebraic model developed by Mellor and Yamada, accounts for the 3D effects of turbulence by calculating explicitly the momentum, heat, and moisture flux divergences in addition to the turbulent kinetic energy. In this study, we present results from idealized simulations in the gray zone that illustrate the benefit of using a fully consistent turbulence closure framework under convective conditions. While the 3D PBL scheme reproduces the evolution of convective features more appropriately than the traditional 1D PBL scheme, we highlight the need to improve the turbulent length scale formulation. Significance Statement The spatial resolution of weather models continues to increase at a rapid rate in accordance with the enhancement of computing power. As a result, smaller-scale atmospheric features become more explicitly resolved. However, most numerical models still ignore the impact of horizontal weather variations on boundary layer flows, which becomes more important at these smaller spatial scales. To address this issue, we have implemented a new modeling approach, using fundamental principles, which accounts for horizontal variability. Our results show that including three-dimensional effects of turbulence is necessary to achieve realistic boundary layer characteristics. This novel technique may be useful for many applications including complex terrain flows, pollutant dispersion, and surface–atmosphere interaction studies.

Juliano, Timothy W.↗

Parameterizations of Entrainment-Mixing Mechanisms and their Effects on Cloud Droplet Spectral Width Based on Numerical Simulations

Entrainment-mixing mechanisms significantly affect cloud droplet number concentration, radius, and spectral shape. Quantitative examination of entrainment-mixing effects on cloud droplet spectral width is lacking. Here, we examine the effects of entrainment-mixing processes on cloud microphysics by 12,218 different setups, each simulated 10 times using the Explicit Mixing Parcel Model (EMPM) driven by the observational data from the Third Tibetan Plateau Atmospheric Scientific Experiment (TIPEX-III) campaign. Parameterizations of entrainment-mixing mechanisms are developed by relating homogeneous mixing degree to transition scale number that depends on the dissipation rate and droplet evaporation time scale. The correlation between relative dispersion of cloud droplet size distribution and homogeneous mixing degree changes from negative to positive with the decreasing homogeneous mixing degree. The different relationships are closely related to the competition between complete and partial droplet evaporation and the number concentration of small droplets, which are quantitatively described by two newly introduced dimensionless numbers. The competition is significantly affected by relative humidity and mixing fraction of entrained air as well as turbulence dissipation rate, but not much by cloud droplet number concentration. Especially, when relative humidity and dissipation rate are high, there is only a negative correlation. This study sheds new light on generalizing the homogeneous/inhomogeneous concept by considering relative dispersion, and also provides parameterizations of entrainment-mixing processes and relative dispersion for atmospheric models.

54 ENVIRONMENTAL SCIENCES↗

Turbulent transport and reactions of plant-emitted hydrocarbons in an Amazonian rain forest

The processes governing the temporal and spatial patterns of isoprene and monoterpenes emitted by a rainforest in the central Amazon region of Brazil is investigated using a combination of field experiments and numerical simulations. Specifically, Large Eddy Simulations (LES) were used to resolve emissions of isoprene and monoterpenes, turbulent transport, and air chemistry. The coupled chemistry-transport LES included the effects of isoprene and monoterpenes reactivity due to reactions with hydroxyl radical and ozone. The LES results are used to compute vertically resolved budgets of isoprene and monoterpenes in the rainforest canopy in response to emissions, turbulent transport, surface deposition, and air chemistry. Our results indicated that emission and dispersion dominated the isoprene budget as the gases were transported out of the canopy space. In a region limited by nitrogen oxides (with prevailing nitric oxide levels of < 0.5 parts per billion), the in-canopy chemical destruction removed approximately 10% of locally emitted monoterpenes. Hydroxyl radical production rates from the ozonolysis of monoterpenes amounted to ≈ 2 x 10 6 radicals cm -3 s -1 and had similar magnitude to the light-dependent hydroxyl radical formation. One key conclusion was that the Amazonia rainforest abundantly emitted monoterpenes whose in-canopy ozonolysis yielded hydroxyl radicals in amounts similar to the magnitude of light-dependent formation. Reactions of monoterpenes and isoprene with hydroxyl radical and ozone were necessary for the maintenance of the Amazon rainforest canopy as a photochemically active environment suitable to generate oxidants and secondary organic aerosols.

54 ENVIRONMENTAL SCIENCES↗

Pass-efficient methods for compression of high-dimensional turbulent flow data

The future of high-performance computing, specifically on future Exascale computers, will presumably see memory capacity and bandwidth fail to keep pace with data generated, for instance, from massively parallel partial differential equation (PDE) systems. Current strategies proposed to address this bottleneck entail the omission of large fractions of data, as well as the incorporation of in situ compression algorithms to avoid overuse of memory. To ensure that post-processing operations are successful, this must be done in a way that a sufficiently accurate representation of the solution is stored. Moreover, in situations where the input/output system becomes a bottleneck in analysis, visualization, etc., or the execution of the PDE solver is expensive, the number of passes made over the data must be minimized. In the interest of addressing this problem, this work focuses on the utility of pass-efficient, parallelizable, low-rank, matrix decomposition methods in compressing high-dimensional simulation data from turbulent flows. Additionally, a particular emphasis is placed on using coarse representation of the data – compatible with the PDE discretization grid – to accelerate the construction of the low-rank factorization. This includes the presentation of a novel single-pass matrix decomposition algorithm for computing the so-called interpolative decomposition. The methods are described extensively and numerical experiments on two turbulent channel flow data are performed. In the first (unladen) channel flow case, compression factors exceeding 400 are achieved while maintaining accuracy with respect to first- and second-order flow statistics. In the particle-laden case, compression factors of 100 are achieved and the compressed data is used to recover particle velocities. These results show that these compression methods can enable efficient computation of various quantities of interest in both the carrier and disperse phases.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Quantum Drude oscillators coupled with Coulomb potential as an efficient model for bonded and non-covalent interactions in atomic dimers

The quantum Drude oscillator (QDO) model has been widely used as an efficient surrogate to describe the electric response properties of matter as well as long-range interactions in molecules and materials. Most commonly, QDOs are coupled within the dipole approximation so that the Hamiltonian can be exactly diagonalized, which forms the basis for the many-body dispersion method [Phys. Rev. Lett. 108, 236402 (2012)]. The dipole coupling is efficient and allows us to study non-covalent many-body effects in systems with thousands of atoms. However, there are two limitations: (i) the need to regularize the interaction at short distances with empirical damping functions and (ii) the lack of multipolar effects in the coupling potential. In this work, we convincingly address both limitations of the dipole-coupled QDO model by presenting a numerically exact solution of the Coulomb-coupled QDO model by means of quantum Monte Carlo methods. We calculate the potential-energy surfaces of homogeneous QDO dimers, analyzing their properties as a function of the three tunable parameters: frequency, reduced mass, and charge. We study the coupled-QDO model behavior at short distances and show how to parameterize this model to enable an effective description of chemical bonds, such as the covalent bond in the H2 molecule.

Chemistry↗

Crystalline symmetry-protected non-trivial topology in prototype compound BaAl 4

The BaAl 4 prototype crystal structure is the most populous of all structure types, and is the building block for a diverse set of sub-structures including the famous ThCr 2 Si 2 family that hosts high-temperature superconductivity and numerous magnetic and strongly correlated electron systems. The MA 4 family of materials (M = Sr, Ba, Eu; A = Al, Ga, In) themselves present an intriguing set of ground states including charge and spin orders, but have largely been considered as uninteresting metals. We predict the exemplary compound BaAl 4 to harbor a three-dimensional Dirac spectrum with non-trivial topology and possible nodal lines crossing the Brillouin zone, wherein one pair of semi-Dirac points with linear dispersion along the k z direction and quadratic dispersion along the k x / k y direction resides on the rotational axis with C 4 v point group symmetry. An extremely large, unsaturating positive magnetoresistance in BaAl 4 despite an uncompensated band structure is revealed, and quantum oscillations and angle-resolved photoemission spectroscopy measurements confirm the predicted multiband semimetal structure with pockets of Dirac holes and a Van Hove singularity (VHS) remarkably consistent with the theoretical prediction. We thus present BaAl 4 as a topological semimetal, casting its prototype status into a role as a building block for a vast array of topological materials.

36 MATERIALS SCIENCE↗

Estimating perturbed stress from 3-D borehole displacements induced by fluid injection in fractured or faulted shales

SUMMARY Hydrofracturing stress measurements in fractured and anisotropic shales are notoriously difficult, because opening of existing geological features tends to prevent the creation of a pure hydraulic fracture perpendicular to the least compressive principal stress. Here we show how adding 3-D borehole-displacement measurements while conducting the hydraulic injection test helps to better constrain the principal stress orientations and magnitudes. We developed a 3-D fully coupled hydromechanical numerical model to analyse the displacement, fluid pressure and injection flow-rate data measured during an injection pressure-step-rate test conducted to activate a faulted borehole interval in the Mont Terri Opalinus Clay (Switzerland). We find that injected fluids can only penetrate the fault when it is at or above the Coulomb failure pressure. Borehole displacement orientations are sensitive to a ∼15° variation in the stress–tensor orientation and a 1 MPa stress magnitude variation. Although some dispersion occurs while rupture is propagating along the fault plane ∼4 m away from the borehole, the maximum density of displacement orientations consistently informs about the stress orientation. Thus, an extended injection step-rate approach coupled with an accurate in situ measurement of the borehole wall displacements can be used to better constrain the local stress field perturbations in fractured shales and in heterogeneous rock in general.

Guglielmi, Yves↗