Search NASA⌕ Search

SEARCH · Search NASA

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

UQpy: A general purpose Python package and development environment for uncertainty quantification

In this paper, we present the UQpy software toolbox, an open-source Python package for general uncertainty quantification (UQ) in mathematical and physical systems. The software serves as both a user-ready toolbox that includes many of the latest methods for UQ in computational modeling and a convenient development environment for Python programmers advancing the field of UQ. The paper presents an introduction to the software's architecture and existing capabilities, divided in the code in a set of modules centered around different UQ tasks such as sampling methods, generation of random processes and random fields, probabilistic inverse modeling, reliability analysis, surrogate modeling, and active learning. The paper also highlights the importance of the RunModel module, which is used to drive simulations in the uncertainty analyses performed in UQpy. This module conveniently allows the user to define computational models directly in Python, or to run simulations from a third-party software in serial or in parallel. To illustrate the various capabilities, two examples are tracked throughout the paper and analyzed repeatedly for various UQ tasks. The first is a Python model solving a nonlinear structural dynamics problem, used to illustrate UQpy's capabilities in sampling and forward propagation of high dimensional random vectors (stochastic processes), and probabilistic inference. The second model is a third-party Abaqus finite element model solving the thermomechanical response of a beam structure. This example is used to illustrate UQpy's capabilities in variance reduction sampling techniques, reliability analysis, surrogate modeling and active learning techniques.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1994

Topological analysis of the CfA redshift survey

We study the topology of large-scale structure in the Center for Astrophysics Redshift Survey, which now includes approximately 12,000 galaxies with limiting magnitude m(sub B) is less than or equal to 15.5. The dense sampling and large volume of this survey allow us to compute the topology on smoothing scales from 6 to 20/h Mpc; we thus examine the topology of structure in both 'nonlinear' and 'linear' regimes. On smoothing scales less than or equal to 10/h Mpc this sample has 3 times the number of resolution elements of samples examined in previous studies. Isodensity surface of the smoothed galaxy density field demonstrate that coherent high-density structures and large voids dominate the galaxy distribution. We compute the genus-threshold density relation for isodensity surfaces of the CfA survey. To quantify phase correlation in these data, we compare the CfA genus with the genus of realizations of Gaussian random fields with the power spectrum measured for the CfA survey. On scales less than or equal to 10/h Mpc the observed genus amplitude is smaller than random phase (96% confidence level). This decrement reflects the degree of phase coherence in the observed galaxy distribution. In other words the genus amplitude on these scales is not good measure of the power spectrum slope. On scales greater than 10/h Mpc, where the galaxy distribution is rougly in the 'linear' regime, the genus ampitude is consistent with the random phase amplitude. The shape of the genus curve reflects the strong coherence in the observed structure; the observed genus curve appears broader than random phase (94% confidence level for smoothing scales less than or equal to 10/h Mpc) because the topolgoy is spongelike over a very large range of density threshold. This departre from random phase consistent with a distribution like a filamentary net of 'walls with holes.' On smoothing scales approaching approximately 20/h Mpc the shape of the CfA genus curve is consistent with random phase. There is very weak evidence for a shift of the genus toward a 'bubble-like' topology. To test cosmological models, we compute the genus for mock CfA surveys drawn from large (L greater than or approximately 400/h Mpc) N-body simulations of three variants of the cold dark matter (CDM) cosmogony. The genus amplitude of the 'standard' CDM model (omega h = 0.5, b = 1.5) differs from the observations (96% confidence level) on smoothing scales is less than or approximately 10/h Mpc. An open CDM model (omega h = 0.2) and a CDM model with nonzero cosmological constant (omega h = 0.24, lambda (sub 0) = 0.6) are consistent with the observed genus amplitude over the full range of smoothing scales. All of these models fail (97% confidence level) to match the broadness of the observed genus curve on smoothing scales is less than or equal to 10/h Mpc.

Vogeley, Michael S.↗
DOE OSTI2023

Unidirectional Spin Hall Magnetoresistance in Antiferromagnetic Heterostructures

Unidirectional spin Hall magnetoresistance (USMR) has been widely reported in the heavy metal/ferromagnet bilayer systems. We observe the USMR in Pt/α-Fe 2 ⁢O 3 bilayers where the α-Fe 2 ⁢O 3 is an antiferromagnetic (AFM) insulator. Systematic field and temperature dependent measurements confirm the magnonic origin of the USMR. The appearance of AFM-USMR is driven by the imbalance of creation and annihilation of AFM magnons by spin orbit torque due to the thermal random field. However, unlike its ferromagnetic counterpart, theoretical modeling reveals that the USMR in Pt/α-Fe 2 ⁢O 3 is determined by the antiferromagtic magnon number with a non-monotonic field dependence. In conclusion, our findings extend the generality of the USMR which pave the ways for the highly sensitive detection of AFM spin state.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1986

Statistical characterization of sea surface geometry for a wave slope field discontinuous in the mean square

Statistics of two-dimensional wave groups, of steep wave events, and of a cascade pattern manifested in the surface geometry in a developed sea state are derived. A systematic view is presented of the spectral moment problem, highlighting its physical meaning and showing its relationship with the theory of random fields. A solution of this problem is suggested that is appropriate to the case of wind-generated surface waves. The solution method treats the surface elevation field as specified on a spatial (temporal) running grid, developing an averaging procedure which employs the Taylor microscale as the mesh size. The technique is illustrated by first exposing errors in direct calculations of the effective surface impedance for a coherently reflected L-band radio wave. The technique is then used to study wave groups and steep waves for a Gaussian, two-dimensional, time-varying surface. Finally, the theory is applied to estimate breaking wave statistics,

Glazman, R. E.↗
DOE OSTI2019

Heterodyne Detection of Radio-Frequency Electric Fields Using Point Defects in Silicon Carbide.

Sensing electric fields with high sensitivity, high spatial resolution, and at radio frequencies can be challenging to realize. Recently, point defects in silicon carbide have shown their ability to measure local electric fields by optical conversion of their charge state. Here, we report the combination of heterodyne detection with charge-based electric field sensing, solving many of the previous limitations of this technique. Owing to the nonlinear response of the charge conversion to electric fields, the application of a separate pump electric field results in a detection sensitivity as low as 1.1(V/cm)/ Hz , with a near-diffraction limited spatial resolution and tunable control of the sensor dynamic range. In addition, we show both incoherent and coherent heterodyne detection, allowing measurements of either unknown random fields or synchronized fields with higher sensitivities. Finally, we demonstrate in-plane vector measurements of the electric field by combining orthogonal pump electric fields. Overall, this work establishes charge-based measurements as highly relevant for solid-state defect sensing.

Wolfowicz, Gary↗
DOE OSTI2025

Multifrequency Models of Black Hole Photon Rings from Low-luminosity Accretion Disks

Images of black holes encode both astrophysical and gravitational properties. Detecting highly lensed features in images can differentiate between these two effects. We present an accretion disk emission model coupled to the Adaptive Analytical Ray Tracing (AART) code that allows a fast parameter space exploration of black hole photon ring images produced from synchrotron emission from 10 to 670 GHz. As an application, we systematically study several disk models and compute their total flux density, average radii, and optical depth. The model parameters are chosen around fiducial values calibrated to general relativistic magnetohydrodynamic (GRMHD) simulations and observations of M87*. For the parameter space studied, we characterize the transition between optically thin and thick regimes and the frequency at which the first photon ring is observable. Our results highlight the need for careful definitions of photon ring radius in the image domain, as in certain models the highly lensed photon ring is dimmer than the direct emission at certain angles. We find that at low frequencies the ring radii are set by the electron temperature, while at higher frequencies the magnetic field strength plays a more significant role, demonstrating how multifrequency analysis can also be used to infer plasma parameters. Lastly, we show how our implementation can qualitatively reproduce multifrequency black hole images from GRMHD simulations when adding time variability to our disk model through Gaussian random fields. This approach provides a new method for simulating observations from the Event Horizon Telescope and the proposed Black Hole Explorer space mission.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1993

Probabilistic nonlinear finite element analysis of composite structures

A probabilistic finite element analysis procedure for laminated composite shells is developed. A total Lagrangian finite element formulation, employing a degenerated three-dimensional laminated composite shell element with the full Green-Lagrange strains and first-order shear deformable kinematics, is used. The first-order second-moment technique for probabilistic finite element analysis of random fields is employed, and results are presented in the form of mean and variance of the structural response. Reliability calculations are made by using the first-order reliability method combined with sensitivity derivatives from the finite element analysis. Both ply-level and micromechanics-level random variables are incorporated, the latter by means of the Aboudi micromechanics model. Two sample problems are solved to verify the accuracy of the procedures developed and to quantify the variability of certain material type/structure combinations. In general, the procedure is quite effective in determining the response statistics and reliability for linear and geometric nonlinear behavior of laminated composite shells.

Engelstad, S. P.↗
DOE OSTI2022

Interfacial Exchange Phenomena Driven by Ferromagnetic Domains

Abstract Interfacial proximity effects in antiferromagnetic/ferromagnetic (AFM/FM) bilayers control the exchange‐bias (EB) phenomena exploited in most spintronic devices, although still is lack of full understanding. Discordant results, including different exchange‐bias field ( H E ), coercivity ( H C ), or blocking temperature ( T B ) found even in similar systems, are usually ascribed to uncontrolled parameters, namely dissimilar interfacial defects, structure, and thicknesses. Here, it is shown in the very same sample that the magnetic domain structure during the magnetization reversal of the FM layer controls those mentioned effects. Simultaneous transport and vectorial‐resolved magnetic measurements performed in a V 2 O 3 /Co system during warming after different field cooling (FC) procedures exhibit a strong dependence on the FC angle and the domain structure of the FM layer. Remarkably, magnetization reversal analysis reveals 35 K of variation in T B and up to a factor of two in H E . These observations can be explained within the random‐field model for the interfacial exchange coupling with a fixed AFM domain structure in contact with a variable (angle‐dependent) FM domain structure. The results highlight the importance of the domain structure and magnetization reversal of the FM layer (not previously considered) in the EB phenomena, with potential to tailor interfacial effects in future spintronic devices.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Far-field thermal imaging below diffraction limit

Non-uniform self-heating and temperature hotspots are major concerns compromising the performance and reliability of submicron electronic and optoelectronic devices. At deep submicron scales where effects such as contact-related artifacts and diffraction limits accurate measurements of temperature hotspots, non-contact thermal characterization can be extremely valuable. In this work, we use a Bayesian optimization framework with generalized Gaussian Markov random field (GGMRF) prior model to obtain accurate full-field temperature distribution of self-heated metal interconnects from their thermoreflectance thermal images (TRI) with spatial resolution 2.5 times below Rayleigh limit for 530nm illumination. Finite element simulations along with TRI experimental data were used to characterize the point spread function of the optical imaging system. In addition, unlike iterative reconstruction algorithms that use ad hoc regularization parameters in their prior models to obtain the best quality image, we used numerical experiments and finite element modeling to estimate the regularization parameter for solving a real experimental inverse problem.

42 ENGINEERING↗
NASA NTRS1993

Nonlinear probabilistic finite element models of laminated composite shells

A probabilistic finite element analysis procedure for laminated composite shells has been developed. A total Lagrangian finite element formulation, employing a degenerated 3-D laminated composite shell with the full Green-Lagrange strains and first-order shear deformable kinematics, forms the modeling foundation. The first-order second-moment technique for probabilistic finite element analysis of random fields is employed and results are presented in the form of mean and variance of the structural response. The effects of material nonlinearity are included through the use of a rate-independent anisotropic plasticity formulation with the macroscopic point of view. Both ply-level and micromechanics-level random variables can be selected, the latter by means of the Aboudi micromechanics model. A number of sample problems are solved to verify the accuracy of the procedures developed and to quantify the variability of certain material type/structure combinations. Experimental data is compared in many cases, and the Monte Carlo simulation method is used to check the probabilistic results. In general, the procedure is quite effective in modeling the mean and variance response of the linear and nonlinear behavior of laminated composite shells.

Engelstad, S. P.↗
DOE OSTI2023

Deep learning Hamiltonians from disordered image data in quantum materials

The capabilities of image probe experiments are rapidly expanding, providing new information about quantum materials on unprecedented length- and timescales. Many such materials feature inhomogeneous electronic properties with intricate pattern formation on the observable surface. This rich spatial structure contains information about interactions, dimensionality, and disorder—a spatial encoding of the Hamiltonian driving the pattern formation. Image recognition techniques from machine learning are an excellent tool for interpreting information encoded in the spatial relationships in such images. Here, we develop a deep learning framework for using the rich information available in these spatial correlations in order to discover the underlying Hamiltonian driving the patterns. We first vet the method on a known case, scanning near-field optical microscopy on a thin film of V⁢O 2 . We then apply our trained convolutional neural network architecture to new optical microscope images of a different V⁢O 2 film as it goes through the metal-insulator transition. We find that a two-dimensional Hamiltonian with both interactions and random field disorder is required to explain the intricate, fractal intertwining of metal and insulator domains during the transition. This detailed knowledge about the underlying Hamiltonian paves the way for using the model to control the pattern formation via, e.g., tailored hysteresis protocols. Finally, we also introduce a distribution-based confidence measure on the results of a multilabel classifier, which does not rely on adversarial training. In addition, we propose a machine-learning-based criterion for diagnosing a physical system's proximity to criticality.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2017

Probabilistic Risk Assessment of CO2 Trapping Mechanisms in a Sandstone CO2-EOR Field in Northern Texas, USA

This study employs a polynomial chaos expansion (PCE) technique to develop the reduced order models (ROMs) integrated with the Monte Carlo simulations to quantify uncertainties of CO 2 storage associated with enhanced oil recovery (EOR) in a commercial-scale CO 2 -EOR field. The Morrow B sandstone at the Farnsworth Unit (FWU) EOR field in northern Texas was selected as a case study. A 3-D reservoir model, including simulation of a five-spot injection/production well pattern, was constructed based on the data from well logs and cores, and an updated seismic interpretation in the FWU. Porosity and permeability were considered as sources of uncertainty. The 1000-realization heterogeneous random fields of porosity and permeability were generated based on the variogram analysis on the core measurements of the FWU. Three model outputs of interest include the amount of CO 2 trapped by three trapping mechanisms, i.e. hydrodynamic trapping, oil dissolution trapping, and aqueous dissolution trapping. The total of 25 runs were conducted with the simulations of a 20-year CO 2 -EOR period, and a 20-year monitoring period to develop the ROMs between uncertain model input parameters (i.e., porosity and permeability in this study) and model outputs at the end of simulation for each grid cell using PCE approach. The high coefficient of determination (R 2 ) and small normalized root mean square error (NRMSE) between the simulated (using reservoir modeling) and predicted (using ROMs) CO 2 storage indicate that the ROMs are acceptable and reliable to be used for the predictions. Given the 1000 Monte Carlo samplings of the model input parameters, cumulative distribution functions (CDFs) and uncertainty bounds (10 th and 90 th percentiles) of model outputs were estimated based on ROMs. At the end of simulation, the most injected CO 2 at the FWU was stored by hydrodynamic trapping (between 121,400 and 166,860 tonnes), followed by oil dissolution trapping (between 55,303 and 76,816 tonnes), and aqueous dissolution trapping (between 1979 and 2751 tonnes).

58 GEOSCIENCES↗
DOE OSTI2021

Bayesian Cross-matching of High Proper-motion Stars in Gaia DR2 and Photometric Metallicities for ∼1.7 million K and M Dwarfs

We present a Bayesian method to cross-match 5,827,988 high proper-motion Gaia sources (μ > 40 mas yr{sup −1}) to various photometric surveys: Two Micron All Sky Survey, AllWISE data release from the Wide-field Infrared Explorer (WISE) mission, Galaxy Evolution Explorer, Radial Velocity Experiment, Sloan Digital Sky Survey, and Panoramic Survey Telescope and Rapid Response System (Pan-STARRS). To efficiently associate these objects across catalogs, we develop a technique that compares the multidimensional distribution of all sources in the vicinity of each Gaia star to a reference distribution of random field stars obtained by extracting all sources in a region on the sky displaced 2′. This offset preserves the local field stellar density and magnitude distribution, allowing us to characterize the frequency of chance alignments. The resulting catalog with Bayesian probabilities >95% has a marginally higher match rate than current internal Gaia data release 2 (DR2) matches for most catalogs. However, a significant improvement is found with Pan-STARRS, where ∼99.8% of the sample within the Pan-STARRS footprint is recovered, as compared to a low ∼20.8% in Gaia DR2. Using these results, we train a Gaussian process regressor to calibrate two photometric metallicity relationships. For dwarfs of 3500 < T{sub eff} < 5280 K, we use metallicity values of 4378 stars from the Apache Point Observatory Galactic Evolution Experiment and Hejazi et al. to calibrate the relationship, producing results with a 1σ precision of 0.12 dex and few systematic errors. We then indirectly infer the metallicity of 4018 stars with 2850 < T{sub eff} < 3500 K, which are wide companions of primaries whose metallicities are estimated with our first regressor, to produce a relationship with a 1σ precision of 0.21 dex and significant systematic errors. Additional work is needed to better remove unresolved binaries from this second sample to reduce these systematic errors.

47 OTHER INSTRUMENTATION↗
NASA NTRS1986

Spatial characteristics of observed precipitation fields: A catalog of summer storms in Arizona, Volume 1

Eight years of summer raingage observations are analyzed for a dense, 93 gage, network operated by the U. S. Department of Agriculture, Agricultural Research Service, in their 150 sq km Walnut Gulch catchment near Tucson, Arizona. Storms are defined by the total depths collected at each raingage during the noon to noon period for which there was depth recorded at any of the gages. For each of the resulting 428 storms, the 93 gage depths are interpolated onto a dense grid and the resulting random field is anlyzed. Presented are: storm depth isohyets at 2 mm contour intervals, first three moments of point storm depth, spatial correlation function, spatial variance function, and the spatial distribution of total rainstorm depth.

Fennessey, N. M.↗
DOE OSTI2017

Insight into 3D micro-CT data: exploring segmentation algorithms through performance metrics

Three-dimensional (3D) micro-tomography (µ-CT) has proven to be an important imaging modality in industry and scientific domains. Understanding the properties of material structure and behavior has produced many scientific advances. An important component of the 3D µ-CT pipeline is image partitioning (or image segmentation), a step that is used to separate various phases or components in an image. Image partitioning schemes require specific rules for different scientific fields, but a common strategy consists of devising metrics to quantify performance and accuracy. The present article proposes a set of protocols to systematically analyze and compare the results of unsupervised classification methods used for segmentation of synchrotron-based data. Additionally, the proposed dataflow for Materials Segmentation and Metrics (MSM) provides 3D micro-tomography image segmentation algorithms, such as statistical region merging (SRM), k-means algorithm and parallel Markov random field (PMRF), while offering different metrics to evaluate segmentation quality, confidence and conformity with standards. Both experimental and synthetic data are assessed, illustrating quantitative results through the MSM dashboard, which can return sample information such as media porosity and permeability. The main contributions of this work are: (i) to deliver tools to improve material design and quality control; (ii) to provide datasets for benchmarking and reproducibility; (iii) to yield good practices in the absence of standards or ground-truth for ceramic composite analysis.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Local structure graph models with higher-order dependence

Local structure graph models (LSGMs) describe random graphs and networks as a Markov random field (MRF)—each graph edge has a specified conditional distribution dependent on explicit neighbourhoods of other graph edges. Centered parameterizations of LSGMs allow for direct control and interpretation of parameters for large- and small-scale structures (e.g., marginal means vs. dependence). Here, we extend this parameterization to account for triples of dependent edges and illustrate the importance of centered parameterizations for incorporating covariates and interpreting parameters. Using a MRF framework, common exponential random graph models are also shown to induce conditional distributions without centered parameterizations and thereby have undesirable features. This work attempts to advance graph models through conditional model specifications with modern parameterizations, covariates and higher-order dependencies.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1991

Non-Gaussian microwave background fluctuations from nonlinear gravitational effects

Whether the statistics of primordial fluctuations for structure formation are Gaussian or otherwise may be determined if the Cosmic Background Explorer (COBE) Satellite makes a detection of the cosmic microwave-background temperature anisotropy delta T(sub CMB)/T(sub CMB). Non-Gaussian fluctuations may be generated in the chaotic inflationary model if two scalar fields interact nonlinearly with gravity. Theoretical contour maps are calculated for the resulting Sachs-Wolfe temperature fluctuations at large angular scales (greater than 3 degrees). In the long-wavelength approximation, one can confidently determine the nonlinear evolution of quantum noise with gravity during the inflationary epoch because: (1) different spatial points are no longer in causal contact; and (2) quantum gravity corrections are typically small-- it is sufficient to model the system using classical random fields. If the potential for two scalar fields V(phi sub 1, phi sub 2) possesses a sharp feature, then non-Gaussian fluctuations may arise. An explicit model is given where cold spots in delta T(sub CMB)/T(sub CMB) maps are suppressed as compared to the Gaussian case. The fluctuations are essentially scale-invariant.

Salopek, D. S.↗
DOE OSTI2024

Classical field approximation of ultralight dark matter: Quantum break times, corrections, and decoherence

The classical field approximation is widely used to better understand the predictions of ultralight dark matter. Here, in this work, we use the truncated Wigner approximation method to test the classical field approximation of ultralight dark matter. This method approximates a quantum state as an ensemble of independently evolving realizations drawn from its Wigner function. The method is highly parallelizable and allows the direct simulation of quantum corrections and decoherence times in systems many times larger than have been previously studied in reference to ultralight dark matter. Our study involves simulation of systems in 1, 2, and 3 spatial dimensions. We simulate three systems, the condensation of a Gaussian random field in three spatial dimensions, a stable collapsed object in three spatial dimensions, and the merging of two stable objects in two spatial dimensions. We study the quantum corrections to the classical field theory in each case. We find that quantum corrections grow exponentially during nonlinear growth with the timescale being approximately equal to the system dynamical time. In stable systems the corrections grow quadratically. We also find that the primary effect of quantum corrections is to reduce the amplitude of fluctuations on the de Broglie scale in the spatial density. Finally, we find that the timescale associated with decoherence due to gravitational coupling to baryonic matter is at least as fast as the quantum corrections due to gravitational interactions. These results are consistent with the predictions of the classical field theory being accurate.

79 ASTRONOMY AND ASTROPHYSICS↗