Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Random fields”

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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
NASA NTRS2000

The Turbulent/Non-Turbulent Interface Bounding a Far-Wake

The velocity fields of a turbulent wake behind a flat plate obtained from the direct numerical simulations of Moser et al. are used to study the structure of the flow in the intermittent zone where there are, alternately, regions of fully turbulent flow and non-turbulent velocity fluctuations either side of a thin randomly moving interface. Comparisons are made with a wake that is 'forced' by amplifying initial velocity fluctuations. There is also a random temperature field T in the flow; T varies between constant values of 0.0 and 1.0 on the sides of the wake. The value of the Reynolds number based on the centreplane mean velocity defect and halfwidth b of the wake is Re approx. = 2000. It is found that the thickness of the continuous interface is about equal to 0.07b, whereas the amplitude of fluctuations of the instantaneous interface displacement y(sub I)(t) is an order of magnitude larger, being about 0.5b. This explains why the mean statistics of vorticity in the intermittent zone can be calculated in terms of the probability distribution of y(sub I) and the instantaneous discontinuity in vorticity across the interface. When plotted as functions of y - y(sub I), the conditional mean velocity (U) and temperature (T) profiles show sharp jumps Delta(U) and Delta(T) at the interface adjacent to a thick zone where (U) and (T) vary much more slowly. Statistics for the vorticity and velocity variances, available in such detail only from DNS data, show how streamwise and spanwise components of vorticity are generated by vortex stretching in the bulges of the interface. Flow fields around the interface, analyzed in terms of the local streamline pattern, confirm previous results that the advancement of the vortical interface into the irrotational flow is driven by large-scale eddy motion. It is argued that because this is an inviscid mechanism the entrainment process is not sensitive to the value of Re, and that small-scale nibbling only plays a subsidiary role. While mean Reynolds stresses decrease gradually in the intermittent zone, conditional stresses are found to decrease sharply towards zero at the interface. Using one-point turbulence models applied to either unconditional or conditional statistics for the turbulent region and then averaged, the entrainment rate E(sub b) would, if calculated exactly, be zero. But if computed with standard computational methods, E(sub b) would be non-zero because of numerical diffusion. It is concluded that the current practice in statistical models of approximating entrainment by a diffusion process is computationally arbitrary and physically incorrect. An analysis shows how E(sub b) is related to Delta(U) and the jump in shear stress at the interface, and correspondingly to Delta(T) and the heat flux.

Bisset, David K.↗
NASA NTRS2009

The Invariant Twist of Magnetic Fields in the Relativistic Jets of Active Galactic Nuclei

The origin of cosmic magnetic (B) fields remains an open question. It is generally believed that very weak primordial B fields are amplified by dynamo processes, but it appears unlikely that the amplification proceeds fast enough to account for the fields presently observed in galaxies and galaxy clusters. In an alternative scenario, cosmic B fields are generated near the inner edges of accretion disks in Active Galactic Nuclei (AGNs) by azimuthal electric currents due to the difference between the plasma electron and ion velocities that arises when the electrons are retarded by interactions with photons. While dynamo processes show no preference for the polarity of the (presumably random) seed field that they amplify, this alternative mechanism uniquely relates the polarity of the poloidal B field to the angular velocity of the accretion disk, resulting in a unique direction for the toroidal B field induced by disk rotation. Observations of the toroidal fields of 29 AGN jets revealed by parsec-scale Faraday rotation measurements show a clear asymmetry that is consistent with this model, with the probability that this asymmetry came about by chance being less than 1 %. This lends support to the hypothesis that the Universe is seeded by B fields that are generated in AGN via this mechanism

Contopoulos, Ioannis↗
DOE OSTI2022

Deconvolving the components of the sign problem

Auxiliary field quantum Monte Carlo simulations of interacting fermions require sampling over a Hubbard-Stratonovich field h introduced to decouple the interactions. The weight for a given configuration involves the products of the determinant of matrices $M_σ(h)$ where σ labels the species, and hence is typically not positive definite. Indeed, the average sign $\langle \mathscr {L} \rangle$ of the determinants goes to zero exponentially with increasing spatial size and decreasing temperature for most Hamiltonians of interest. This statement, however, does not explicitly separate two possible origins for the vanishing of $\langle \mathscr {L} \rangle$. Does $\langle \mathscr {L} \rangle$ → 0 because randomly chosen field configurations have det[M(h)] < 0, or does the sign problem arise because the specific subset of configurations chosen by the weighting function have a greater preponderance of negative values? In the latter case, the process of weighting the configurations with |det[M(h)]| might steer the simulation to a region of configuration space of h where positive and negative determinants are equally likely, even though randomly chosen h would preferentially have determinants with a single dominant sign. Here in this paper, we address the relative importance of these two mechanisms for the vanishing of $\langle \mathscr {L} \rangle$ in quantum simulations.

36 MATERIALS SCIENCE↗
NASA NTRS1991

Direct simulation of turbulent combustion

Understanding and modeling of turbulent combustion are key problems in the computation of numerous practical systems. Because of the lack of analytical theories in this field and of the difficulty of performing precise experiments, direct simulation appears to be one of the most attractive tools to use in addressing this problem. From Sep. 1989 - Sep. 1990, the following problems were studied using direct simulation but also experiments and theoretical models: (1) the influence of curvature on premixed flame fronts; (2) the construction of a model for the flame stretch based on direct simulation and multifractal analysis; and (3) the simulation of the interaction between random flow field and a flame front. The general objective was to improve knowledge of turbulent combustion but also to use this information for turbulent combustion models. The transfer of direct simulation results towards models is also in progress. A submodel to evaluate the flame stretch was derived from the present study and incorporated in the Coherent Flame model used in France for piston and aircraft engines.

Poinsot, Thierry J.↗
DOE OSTI2024

Disorder-induced local strain distribution in Y-substituted TmVO 4

We report an investigation of the effect of substitution of Y for Tm in Tm 1-x ⁢Y x VO 4 via low-temperature heat capacity measurements, with the yttrium content x varying from 0 to 0.997. Because the Tm ions support a local quadrupolar (nematic) moment, they act as reporters of the local strain state in the material, with the splitting of the ion's non-Kramers crystal field ground state proportional to the quadrature sum of the in-plane tetragonal symmetry-breaking transverse and longitudinal strains experienced by each ion individually. Analysis of the heat capacity, therefore, provides detailed insights into the distribution of local strains that arise as a consequence of the chemical substitution. These local strains suppress long-range quadrupole order for x > 0.22, and result in a broad Schottky-like feature for higher concentrations. Heat capacity data are compared to expectations for a distribution of uncorrelated (random) strains. For dilute Tm concentrations, the heat capacity cannot be accounted for by randomly distributed strains, demonstrating the presence of significant strain correlations between sites. For intermediate Tm concentrations, these correlations must still exist, but the data cannot be distinguished from that which would be obtained from a two-dimensional Gaussian distribution. The crossover between these limits is discussed in terms of the interplay of key lengthscales in the substituted material. Furthermore, the central result of this work, namely that local strains arising from chemical substitution are not uncorrelated, has implications for the range of validity of theoretical models based on random effective fields that are used to describe such chemically substituted materials, particularly when electronic nematic correlations are present.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Learning and meta-learning of stochastic advection–diffusion–reaction systems from sparse measurements

Physics-informed neural networks (PINNs) were recently proposed in [18] as an alternative way to solve partial differential equations (PDEs). A neural network (NN) represents the solution, while a PDE-induced NN is coupled to the solution NN, and all differential operators are treated using automatic differentiation. Here, we first employ the standard PINN and a stochastic version, sPINN, to solve forward and inverse problems governed by a non-linear advection–diffusion–reaction (ADR) equation, assuming we have some sparse measurements of the concentration field at random or pre-selected locations. Subsequently, we attempt to optimise the hyper-parameters of sPINN by using the Bayesian optimisation method (meta-learning) and compare the results with the empirically selected hyper-parameters of sPINN. In particular, for the first part in solving the inverse deterministic ADR, we assume that we only have a few high-fidelity measurements, whereas the rest of the data is of lower fidelity. Hence, the PINN is trained using a composite multi-fidelity network, first introduced in [12], that learns the correlations between the multi-fidelity data and predicts the unknown values of diffusivity, transport velocity and two reaction constants as well as the concentration field. For the stochastic ADR, we employ a Karhunen–Loève (KL) expansion to represent the stochastic diffusivity, and arbitrary polynomial chaos (aPC) to represent the stochastic solution. Correspondingly, we design multiple NNs to represent the mean of the solution and learn each aPC mode separately, whereas we employ a separate NN to represent the mean of diffusivity and another NN to learn all modes of the KL expansion. For the inverse problem, in addition to stochastic diffusivity and concentration fields, we also aim to obtain the (unknown) deterministic values of transport velocity and reaction constants. The available data correspond to 7spatial points for the diffusivity and 20 space–time points for the solution, both sampled 2000 times. We obtain good accuracy for the deterministic parameters of the order of 1–2% and excellent accuracy for the mean and variance of the stochastic fields, better than three digits of accuracy. In the second part, we consider the previous stochastic inverse problem, and we use Bayesian optimisation to find five hyper-parameters of sPINN, namely the width, depth and learning rate of two NNs for learning the modes. Here, we obtain much deeper and wider optimal NNs compared to the manual tuning, leading to even better accuracy, i.e., errors less than 1% for the deterministic values, and about an order of magnitude less for the stochastic fields.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Simulations of axion minihalos

The axion, motivated as a solution to the strong CP problem, is also a viable dark matter candidate. The axion field takes random values in causally disconnected regions if the symmetry breaking that establishes the particle occurs after inflation, leading to white-noise density fluctuations at low wave numbers and forming dense minihalos with sub-planetary masses subsequently. There have been two recent proposals that appear capable of testing this scenario, namely using pulsar timing arrays and studying cosmological microlensing caustics. Motivated by these proposals, we use N-body simulations to study the formation of substructures from white-noise density fluctuations. The density profiles of our relaxed axion minihalos can be described by the Navarro–Frenk–White profile, and the minihalos’ concentration number agrees well with a simple, physically-motivated model. We develop a semi-analytic formula to fit the mass function from our simulation, which agrees broadly at different redshifts and only differs at factor of two level from classic halo mass functions. This analytic mass function allows us to consider uncertainties in the post-inflation axion scenario, as well as extrapolate our high-redshift simulations results to the present. Our work estimates the present-day abundance of axion substructures, as is necessary for predicting their effect on cosmological microlensing caustics and pulsar timing. Our calculations suggest that if pulsar timing and microlensing probes can reach recent sensitivity forecasts, they may be sensitive to the post-inflation axion dark matter scenario, even when accounting for uncertainties pertaining to axion strings. For pulsar timing, the most significant caveat is whether axion minihalos are disrupted by stars, which our estimates show is mildly important at the most relevant masses. In conclusion, as our gravitational simulations are scale invariant, the results can be extended to models where the dark matter is comprised of other axion-like particles and even clusters of primordial black holes.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Entanglement and confinement in coupled quantum systems

We study some general properties of coupled quantum systems. We consider simple interactions between two copies of identical Hamiltonians such as the SYK model, Pauli spin chains with random magnetic field and harmonic oscillators. Such couplings make the ground states close to the thermofield double states of the uncoupled Hamiltonians. For the coupled SYK model, we push the numerical computation further towards the thermodynamic limit so that an extrapolation in the size of the system is possible. We find good agreement between the extrapolated numerical result and the analytic result in the large-q limit. We also consider the coupled gauged matrix model and vector model, and argue that the deconfinement is associated with the loss of the entanglement, similarly to the previous observation for the coupled SYK model. The understanding of the microscopic mechanism of the confinement/deconfinement transition enables us to estimate the quantum entanglement precisely, and backs up the dual gravity interpretation which relates the deconfinement to the disappearance of the wormhole. Our results demonstrate the importance of the entanglement between the color degrees of freedom in the emergence of the bulk geometry from quantum field theory via holography.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Localization dynamics in a centrally coupled system

In systems in which interactions couple a central degree of freedom and a bath, one would expect signatures of the bath's phase to be reflected in the dynamics of the central degree of freedom. This has been recently explored in connection with many-body localized baths coupled with a central qubit or a single-cavity mode - systems with growing experimental relevance in various platforms. Such models also have an interesting connection with Floquet many-body localization via quantizing the external drive, although this has been relatively unexplored. Here we adapt the multilayer multiconfigurational time-dependent Hartree (ML-MCTDH) method, a well-known tree tensor network algorithm, to numerically simulate the dynamics of a central degree of freedom, represented by a d-level system (qudit), coupled to a disordered interacting one-dimensional spin bath. ML-MCTDH allows us to reach ≈10 2 lattice sites, a far larger system size than what is feasible with exact diagonalization or kernel polynomial methods. From the intermediate time dynamics, we find a well-defined thermodynamic limit for the qudit dynamics upon appropriate rescaling of the system-bath coupling. The spin system shows similar scaling collapse in the Edward-Anderson spin-glass order parameter or entanglement entropy at relatively short times. At longer timescales, we see slow growth of the entanglement, which may arise from dephasing mechanisms in the localized system or long-range interactions mediated by the central degree of freedom. Similar signs of localization are shown to appear as well with unscaled system-bath coupling.

1-dimensional spin chains↗
DOE OSTI2020

Characterization of quantum chaos by two-point correlation functions

In this work, we propose a characterization of quantum many-body chaos: given a collection of simple operators, the set of all possible pair correlations between these operators can be organized into a matrix with a random-matrix-like spectrum. This approach is particularly useful for locally interacting systems, which do not generically show exponential Lyapunov growth of out-of-time-ordered correlators. We demonstrate the validity of this characterization by numerically studying the Sachdev-Ye-Kitaev model and a one-dimensional spin chain with random magnetic field (XXZ model).

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Lattice QCD Calculation of Electroweak Box Contributions to Superallowed Nuclear and Neutron Beta Decays

We present the first lattice QCD calculation of the universal axial γ W -box contribution □ γ W V A to both superallowed nuclear and neutron beta decays. This contribution emerges as a significant component within the theoretical uncertainties surrounding the extraction of | V u d | from superallowed decays. Our calculation is conducted using two domain wall fermion ensembles at the physical pion mass. To construct the nucleon four-point correlation functions, we employ the random sparsening field technique. Furthermore, we incorporate long-distance contributions to the hadronic function using the infinite-volume reconstruction method. Upon performing the continuum extrapolation, we arrive at □ γ W V A = 3.65 ( 7 ) lat ( 1 ) PT × 10 − 3 . Consequently, this yields a slightly higher value of | V u d | = 0.973 86 ( 11 ) exp ( 9 ) RC ( 27 ) NS , reducing the previous 2.1 σ tension with the CKM unitarity to 1.8 σ . Additionally, we calculate the vector γ W -box contribution to the axial charge g A , denoted as □ γ W V V , and explore its potential implications. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Length scales and scale-free dynamics of dislocations in dense solid solutions

The fundamental interactions between an edge dislocation and a random solid solution are studied by analyzing dislocation line roughness profiles obtained from molecular dynamics simulations of Fe 0.70 Ni 0.11 Cr 0.19 over a range of stresses and temperatures. These roughness profiles reveal the hallmark features of a depinning transition. Namely, below a temperature-dependent critical stress, the dislocation line exhibits roughness in two different length scale regimes which are divided by a so-called correlation length. This correlation length increases with applied stress and at the critical stress (depinning transition or yield stress) formally goes to infinity. Above the critical stress, the line roughness profile converges to that of a random noise field. Motivated by these results, a physical model is developed based on the notion of coherent line bowing over all length scales below the correlation length. Above the correlation length, the solute field prohibits such coherent line bow outs. Using this model, we identify potential gaps in existing theories of solid solution strengthening and show that recent observations of length-dependent dislocation mobilities can be rationalized.

36 MATERIALS SCIENCE↗
DOE OSTI2021

PastoralScape: An Environment-Driven Model of Vaccination Decision Making Within Pastoralist Groups in East Africa

Economic and cultural resilience among pastoralists in East Africa is threatened by the interconnected forces of climate change, contagious diseases spread and evolving national and international trade. A key factor in the resilience of livestock that communities depend on is human decision making regarding vaccination against prevalent diseases such as Rift Valley fever and Contagious Bovine Pleuropneumonia. This paper describes an agent-based model that couples models of disease propagation, animal health, human decision making, and external GIS data sources capturing measures of foraging condition. We describe the design of the sub-models, their coupling, and demonstrate the sensitivity of the model to parameters that relate to controllable factors such as government and NGO information sources that can influence human decision making patterns. This model is intended to form the basis upon which richer economic and human factor models can be built.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Galaxy Cluster Contribution to the Diffuse Extragalactic Ultraviolet Background

The diffuse ultraviolet background radiation has been mapped over most of the sky with 2′ resolution using data from the Galaxy Evolution Explorer survey. We utilize this map to study the correlation between the UV background and clusters of galaxies discovered via the Sunyaev–Zeldovich effect in the Planck survey. We use only high Galactic latitude (|b|>60{sup ∘}) galaxy clusters to avoid contamination by Galactic foregrounds, and we only analyze clusters with a measured redshift. This leaves us with a sample of 142 clusters over the redshift range of 0.02 ≤ z ≤ 0.72, which we further subdivide into four redshift bins. In analyzing our stacked samples binned by redshift, we find evidence for a central excess of UV background light compared to local backgrounds for clusters with z < 0.3. We then stacked these z < 0.3 clusters to find a statistically significant excess of 12 ± 2.3 photon cm{sup −2} s{sup −1} sr{sup −1} Å{sup −1} over the median of ∼380 photon cm{sup −2} s{sup −1} sr{sup −1} Å{sup −1} measured around random blank fields. We measure the stacked radial profile of these clusters, and find that the excess UV radiation decays to the level of the background at a radius of ∼1 Mpc, roughly consistent with the maximum radial extent of the clusters. Analysis of possible physical processes contributing to the excess UV brightness indicates that non-thermal emission from relativistic electrons in the intracluster medium and faint, unresolved UV emission from cluster member galaxies and intracluster light are likely the dominant contributors.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1966

Cosmic-ray propagation. I.

Charged particle motion in random magnetic field, describing time evolution of particle distribution in pitch angle and position in terms of Fokker- Planck coefficients

INTERPLANETARY MAGNETIC FIELD↗
NASA NTRS1968

Generalized stochastic equations.

Generalized stochastic equations of probability density distribution for particle system in random electric field, noting relationship to Fokker-Planck equation

So, K. C.↗