Search NASA⌕ Search

SEARCH · Search NASA

Results for “explicit correlation”

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 667 records · Page 37
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Stochastic generation of explicit pore structures by thresholding Gaussian random fields

We provide a description and computational investigation of an efficient method to stochastically generate realistic pore structures. Smolarkiewicz and Winter introduced this specific method in pores resolving simulation of Darcy flows (Smolarkiewicz and Winter, 2010 [1]) without giving a complete formal description or analysis of the method, or indicating how to control the parameterization of the ensemble. We address both issues in this paper. The method consists of two steps. First, a realization of a correlated Gaussian field, or topography, is produced by convolving a prescribed kernel with an initial field of independent, identically distributed random variables. The intrinsic length scales of the kernel determine the correlation structure of the topography. Next, a sample pore space is generated by applying a level threshold to the Gaussian field realization: points are assigned to the void phase or the solid phase depending on whether the topography over them is above or below the threshold. Hence, the topology and geometry of the pore space depend on the form of the kernel and the level threshold. Manipulating these two user prescribed quantities allows good control of pore space observables, in particular the Minkowski functionals. Extensions of the method to generate media with multiple pore structures and preferential flow directions are also discussed. To demonstrate its usefulness, the method is used to generate a pore space with physical and hydrological properties similar to a sample of Berea sandstone. -- Graphical abstract: -- Highlights: •An efficient method to stochastically generate realistic pore structures is provided. •Samples are generated by applying a level threshold to a Gaussian field realization. •Two user prescribed quantities determine the topology and geometry of the pore space. •Multiple pore structures and preferential flow directions can be produced. •A pore space based on Berea sandstone is generated.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Spatial and temporal agreement in climate model simulations of the Interdecadal Pacific Oscillation

Accelerated warming and hiatus periods in the long-term rise of Global Mean Surface Temperature (GMST) have, in recent decades, been associated with the Interdecadal Pacific Oscillation (IPO). Critically, decadal climate prediction relies on the skill of state-of-the-art climate models to reliably represent these low-frequency climate variations. We undertake a systematic evaluation of the simulation of the IPO in the suite of Coupled Model Intercomparison Project 5 (CMIP5) models. We track the IPO in pre-industrial (control) and all-forcings (historical) experiments using the IPO tripole index (TPI). The TPI is explicitly aligned with the observed spatial pattern of the IPO, and circumvents assumptions about the nature of global warming. We find that many models underestimate the ratio of decadal-to-total variance in sea surface temperatures (SSTs). However, the basin-wide spatial pattern of positive and negative phases of the IPO are simulated reasonably well, with spatial pattern correlation coefficients between observations and models spanning the range 0.4–0.8. Deficiencies are mainly in the extratropical Pacific. Models that better capture the spatial pattern of the IPO also tend to more realistically simulate the ratio of decadal to total variance. Of the 13% of model centuries that have a fractional bias in the decadal-to-total TPI variance of 0.2 or less, 84% also have a spatial pattern correlation coefficient with the observed pattern exceeding 0.5. This result is highly consistent across both IPO positive and negative phases. This is evidence that the IPO is related to one or more inherent dynamical mechanisms of the climate system.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Six-point functions and collisions in the black hole interior

In the eternal AdS black hole geometry, we consider two signals sent from the boundaries into the black hole interior shared between the two asymptotic regions. We compute three different out-of-time-order six-point functions to quantify various properties of the collision of these signals behind the horizons: (i) We diagnose the strength of the collision by probing the two-signal state on a late time slice with boundary operators. (ii) We quantify two-sided operator growth, which provides a dual description of the signals meeting in the black hole interior, in terms of the quantum butterfly effect and quantum circuits. (iii) We consider an explicit coupling between the left and right CFTs to make the wormhole traversable and extract information about the collision product from behind the horizon. At a technical level, our results rely on the method of eikonal resummation to obtain the relevant gravitational contributions to Lorentzian six-point functions at all orders in the GN-expansion. We observe that such correlation functions display an intriguing factorization property. We corroborate these results with geodesic computations of six-point functions in two- and three-dimensional gravity.

2D Gravity↗
DOE OSTI2023

Kβ X-ray Emission Spectroscopy of Cu(I)-Lytic Polysaccharide Monooxygenase: Direct Observation of the Frontier Molecular Orbital for H 2 O 2 Activation

Lytic polysaccharide monooxygenases (LPMOs) catalyze the degradation of recalcitrant carbohydrate polysaccharide substrates. These enzymes are characterized by a mononuclear Cu(I) active site with a three-coordinate T-shaped “His-brace” configuration including the N-terminal histidine and its amine group as ligands. Here, this study explicitly investigates the electronic structure of the d 10 Cu(I) active site in a LPMO using Kβ X-ray emission spectroscopy (XES). The lack of inversion symmetry in the His-brace site enables the 3d/p mixing required for intensity in the Kβ valence-to-core (VtC) XES spectrum of Cu(I)-LPMO. These Kβ XES data are correlated to density functional theory (DFT) calculations to define the bonding, and in particular, the frontier molecular orbital (FMO) of the Cu(I) site. These experimentally validated DFT calculations are used to evaluate the reaction coordinate for homolytic cleavage of the H 2 O 2 O–O bond and understand the contribution of this FMO to the low barrier of this reaction and how the geometric and electronic structure of the Cu(I)-LPMO site is activated for rapid reactivity with H 2 O 2 .

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

A practical computational method for the anisotropic redshift-space three-point correlation function

Here, we present an algorithm enabling computation of the anisotropic redshift-space galaxy threepoint correlation function ( 3 PCF) scaling as N 2 , with N the number of galaxies. Our previous work showed how to compute the isotropic 3 PCF with this scaling by expanding the radially binned density field around each galaxy in the survey into spherical harmonics and combining these coefficients to form multipole moments. The N 2 scaling occurred because this approach never explicitly required the relative angle between a galaxy pair about the primary galaxy. Here, we generalize this work, demonstrating that in the presence of azimuthally symmetric anisotropy produced by redshift-space distortions (RSD), the 3 PCF can be described by two triangle side lengths, two independent total angular momenta, and a spin. This basis for the anisotropic 3 PCF allows its computation with negligible additional work over the isotropic 3 PCF. We also present the covariance matrix of the anisotropic 3 PCF measured in this basis. Our algorithm tracks the full 5 D redshift-space 3 PCF, uses an accurate line of sight to each triplet, is exact in angle, and easily handles edge correction. It will enable use of the anisotropic large-scale 3 PCF as a probe of RSD in current and upcoming large-scale redshift surveys.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

LA100 Equity Strategies. Chapter 12: Distribution Grid Upgrades for Equitable Resilience and Solar, Storage, and Electric Vehicle Access

The LA100 Equity Strategies project integrates community guidance with robust research, modeling, and analysis to identify strategy options that can increase equitable outcomes in Los Angeles' clean energy transition. As Los Angeles transitions toward clean energy, existing distribution grid infrastructure will need to be updated and expanded to support reliable service during routine operations, enable interconnection with distributed energy resources and electrified loads, and provide access to energy-related services during disasters. This chapter focuses on equity in distribution grid upgrades, reliability, and resilience in Los Angeles. Specifically, NREL performed grid upgrade and resilience analyses using a detailed model of the distribution grid and income-differentiated household load profiles, electric vehicle (EV) adoption patterns, distributed solar adoption, and grid reliability to explore two key questions to inform how the City of Los Angeles can ensure a resilient and reliable distribution grid for all communities during the clean energy transition: Where can distribution system upgrades can be prioritized to enable equitable access to, and adoption of, clean energy technologies and how can Los Angeles provide equitable, resilient access to electricity-related services (e.g., health care, food) during disaster events like earthquakes and flooding? The electric distribution system is the "last mile" of the grid, linking the multistate bulk power system with customers; new loads, including EVs; and distributed energy resources, such as customer and community solar and storage. This analysis focuses on the 4.8-kilovolt (kV) system, including service transformers that represent the utility-side of the grid connection for most residential customers. Chapter 17 looks at the customer-side of the grid connection with a focus on electric panel upgrade needs. The transition toward clean energy can put additional stress on the distribution system from distributed energy resources and electrification - especially EVs and increased use of electricity for heating, cooling, cooking, and hot water. This stress, measured here as the number of equipment overloads and voltage violations, correlates strongly to grid reliability and therefore is used as a proxy for understanding additional upgrades needed and to help ensure equitable access to electrification and distributed energy resources. NREL also conducted community resilience analysis to examine customer-level access to both electricity and a larger range of services, such as hospitals and grocery stores during a disaster. This analysis explicitly considers equity to understand differences in current resilience and resilience strategies to effectively improve critical services access for all Angelenos. Research was guided by input from the community engagement process, and associated equity strategies are presented in alignment with that guidance.

14 SOLAR ENERGY↗
NASA NTRS1991

Probabilistic Structures Analysis Methods (PSAM) for select space propulsion system components

The basic formulation for probabilistic finite element analysis is described and demonstrated on a few sample problems. This formulation is based on iterative perturbation that uses the factorized stiffness on the unperturbed system as the iteration preconditioner for obtaining the solution to the perturbed problem. This approach eliminates the need to compute, store and manipulate explicit partial derivatives of the element matrices and force vector, which not only reduces memory usage considerably, but also greatly simplifies the coding and validation tasks. All aspects for the proposed formulation were combined in a demonstration problem using a simplified model of a curved turbine blade discretized with 48 shell elements, and having random pressure and temperature fields with partial correlation, random uniform thickness, and random stiffness at the root.

Source record↗
NASA NTRS2002

Best Practices for Crash Modeling and Simulation

Aviation safety can be greatly enhanced by the expeditious use of computer simulations of crash impact. Unlike automotive impact testing, which is now routine, experimental crash tests of even small aircraft are expensive and complex due to the high cost of the aircraft and the myriad of crash impact conditions that must be considered. Ultimately, the goal is to utilize full-scale crash simulations of aircraft for design evaluation and certification. The objective of this publication is to describe "best practices" for modeling aircraft impact using explicit nonlinear dynamic finite element codes such as LS-DYNA, DYNA3D, and MSC.Dytran. Although "best practices" is somewhat relative, it is hoped that the authors' experience will help others to avoid some of the common pitfalls in modeling that are not documented in one single publication. In addition, a discussion of experimental data analysis, digital filtering, and test-analysis correlation is provided. Finally, some examples of aircraft crash simulations are described in several appendices following the main report.

Fasanella, Edwin L.↗
DOE OSTI2026

Correlation function metrology for warm dense matter: Recent developments and practical guidelines

X-ray Thomson scattering (XRTS) has emerged as a valuable diagnostic for matter under extreme conditions, as it captures the intricate many-body physics of the probed sample. Recent advances, such as the model-free temperature diagnostic of Dornheim et al. [Nat. Commun. 13 , 7911 (2022)], have demonstrated how much information can be extracted directly within the imaginary-time formalism. However, since the imaginary-time formalism is a concept often difficult to grasp, we provide here a systematic overview of its theoretical foundations and explicitly demonstrate its practical applications to temperature inference, including relevant subtleties. Furthermore, we present recent developments that enable the determination of the absolute normalization, Rayleigh weight, and density from XRTS measurements without reliance on uncontrolled model assumptions. Finally, we outline a unified workflow that guides the extraction of these key observables, offering a practical framework for applying the method to interpret experimental measurements.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Dynamical vertex approximation for many-electron systems with spontaneously broken SU(2) symmetry

We generalize the formalism of the dynamical vertex approximation (DΓA)—a diagrammatic extension of the dynamical mean-field theory (DMFT)—to treat magnetically ordered phases. To this aim, we start by concisely illustrating the many-electron formalism for performing ladder resummations of Feynman diagrams in systems with broken SU(2) symmetry associated to ferromagnetic (FM) or antiferromagnetic (AF) order. We then analyze the algorithmic simplifications introduced by taking the local approximation of the two-particle irreducible vertex functions in the Bethe-Salpeter equations, which defines the ladder implementation of DΓA for magnetic systems. The relation of this assumption with the DMFT limit of large coordination-number/high dimensions is explicitly discussed. As a last step, we derive the expression for the ladder DΓA self-energy in the FM- and AF-ordered phases of the Hubbard model. The physics emerging in the AF-ordered case is explicitly illustrated by means of approximated calculations based on a static mean-field input for DΓA equations. The results obtained capture fundamental aspects of both metallic and insulating ground states of two-dimensional antiferromagnets, providing a reliable compass for future, more extensive applications of our approach. Furthermore, possible routes to further develop diagrammatic-based treatments of magnetic phases in correlated electron systems are briefly outlined in the Conclusions.

2-dimensional systems↗
DOE OSTI2021

Geodyn Material Library: Pseudocap models for dry porous tocks

This report describes the second edition of the Pseudocap Strength models for porous rocks implemented in GEODYN material library. The first model was developed in 2007 and calibrated for concrete. Then, the model parameters were calibrated based on triaxial tests reported for limestones and sandstones of various porosities. In these models some key parameters were chosen as functions of the reference porosity,Φ. Since then multiple modifications were implemented in the model, therefore, it has been recalibrated for some common porous materials such as limestones, sandstones, alluvium, tuffs and granite. Two types of models are described in this repot. The first type (called Pseudocap Model or PM) is for rocks from a specific location. Parameters were calibrated for several specific geologic materials. The second type (called Generic Pseudocap Model or GPM) is useful for the sites where only basic information (rock type, porosity) is available. Generic models include built-in correlations between porosities and other mechanical properties observed for certain rock types. The models of both types were validated by comparing not only to quasi-static triaxial tests for these materials but also to shock Hugoniot data and spherical explosion data for some materials. All models were derived in the frame of isotropic plasticity. They are designed to be used in explicit finite element/difference codes. Tangent stiffness tensor is not provided but can be calculated numerically for the model to be used in implicit finite element codes. For an isotropic material the stress can be decomposed into volumetric and deviatoric parts. The volumetric part is modeled using an Equation of state (EOS) which calculates the pressure and the bulk sound speed as functions of the internal specific energy and density. Here a simple, Mie-Gruneisen EOS is presented, but tabulated EOS (LEOS) provided by the library can be used as well. The stress is limited by the yield surface which depends on three invariants of the stress tensor and specific internal energy. In addition, to capture the strain-rate dependence a simple multiplier is used for the yield surface which depends on the equivalent plastic strain rate. The failure surface (the ultimate yield, Y f , defined later) is chosen in the Hoek-Brown form, commonly used in rock mechanics. It includes measurable parameters such as Unconfined Compressive Strength (UCS) as well as scale parameters characterizing the quality of the rock such as GSI (Geologic Strength Index). Thus, even though the model is calibrated for small samples it offers a way to extrapolate the strength to the field scale using geological characterization of the rock mass. The model captures effects of brittle-ductile transition in rocks by introducing a cap multiplier to the yield function. The rate of dilatancy (bulking) is proportional to the slope of the yield surface affected by the cap. Therefore, it takes place only at low confinements when the pressure is less than the brittle-ductile transition pressure, P BD . On the contrary, the porous compaction takes place at pressures higher than P BD . The cap moves as the porosity is compacted or new porosity is generated due to dilatancy. The porous compaction is modeled using an evolution equation which includes deviatoric stress so that the onset of compaction corresponds to the cap surface. The model captures effects of shear-enhanced compaction which is an important for porous rocks. Section 2 describes the modeling framework and Section 3 presents the model calibration procedure. Section 4 compares experimental data for various rocks versus model predictions. The model parameters used for this comparison are given in Appendix. The files with material constants are available with the latest GEODYN material library distribution.

58 GEOSCIENCES↗
NASA NTRS1990

Nonlinear damping model for flexible structures

The study of nonlinear damping problem of flexible structures is addressed. Both passive and active damping, both finite dimensional and infinite dimensional models are studied. In the first part, the spectral density and the correlation function of a single DOF nonlinear damping model is investigated. A formula for the spectral density is established with O(Gamma(sub 2)) accuracy based upon Fokker-Planck technique and perturbation. The spectral density depends upon certain first order statistics which could be obtained if the stationary density is known. A method is proposed to find the approximate stationary density explicitly. In the second part, the spectral density of a multi-DOF nonlinear damping model is investigated. In the third part, energy type nonlinear damping model in an infinite dimensional setting is studied.

Zang, Weijian↗
DOE OSTI2011

Nyquist-Shannon sampling theorem applied to refinements of the atomic pair distribution function

We have systematically studied the optimal real-space sampling of atomic pair distribution (PDF) data by comparing refinement results from oversampled and resampled data. Based on nickel and a complex perovskite system, we show that not only is the optimal sampling bounded by the Nyquist interval described by the Nyquist-Shannon (NS) sampling theorem as expected, but near this sampling interval, the data points in the PDF are minimally correlated, which results in more reliable uncertainty estimates in the modeling. Surprisingly, we find that PDF refinements quickly become unstable for data on coarser grids. Although the Nyquist-Shannon sampling theorem is well known, it has not been applied to PDF refinements, despite the growing popularity of the PDF method and its adoption in a growing number of communities. Here, we give explicit expressions for the application of NS sampling theorem to the PDF case, and establish through modeling that it is working in practice, which lays the groundwork for this to become more widely adopted. This has implications for the speed and complexity of possible refinements that can be carried out many times faster than currently with no loss of information, and it establishes a theoretically sound limit on the amount of information contained in the PDF that will prevent over-parametrization during modeling.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

Point model equations for neutron correlation counting: Extension of Böhnel's equations to any order

Various methods of autocorrelation neutron analysis may be used to extract information about a measurement item containing spontaneously fissioning material. The two predominant approaches being the time correlation analysis (that make use of a coincidence gate) methods of multiplicity shift register logic and Feynman sampling. The common feature is that the correlated nature of the pulse train can be described by a vector of reduced factorial multiplet rates. We call these singlets, doublets, triplets etc. Within the point reactor model the multiplet rates may be related to the properties of the item, the parameters of the detector, and basic nuclear data constants by a series of coupled algebraic equations – the so called point model equations. Solving, or inverting, the point model equations using experimental calibration model parameters is how assays of unknown items is performed. Currently only the first three multiplets are routinely used. In this work we develop the point model equations to higher order multiplets using the probability generating functions approach combined with the general derivative chain rule, the so called Faà di Bruno Formula. Explicit expression up to 5th order are provided, as well the general iterative formula to calculate any order. This study represents the first necessary step towards determining if higher order multiplets can add value to nondestructive measurement practice for nuclear materials control and accountancy.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2014

CONNECTING STAR FORMATION QUENCHING WITH GALAXY STRUCTURE AND SUPERMASSIVE BLACK HOLES THROUGH GRAVITATIONAL HEATING OF COOLING FLOWS

Recent observations suggested that star formation quenching in galaxies is related to galaxy structure. Here we propose a new mechanism to explain the physical origin of this correlation. We assume that while quenching is maintained in quiescent galaxies by a feedback mechanism, cooling flows in the hot halo gas can still develop intermittently. We study cooling flows in a large suite of around 90 hydrodynamic simulations of an isolated galaxy group, and find that the flow development depends significantly on the gravitational potential well in the central galaxy. If the galaxy's gravity is not strong enough, cooling flows result in a central cooling catastrophe, supplying cold gas and feeding star formation to galactic bulges. When the bulge grows prominent enough, compressional heating starts to offset radiative cooling and maintains cooling flows in a long-term hot mode without producing a cooling catastrophe. Our model thus describes a self-limited growth channel for galaxy bulges and naturally explains the connection between quenching and bulge prominence. In particular, we explicitly demonstrate that M{sub ∗}/R{sub eff}{sup 1.5} is a good structural predictor of quenching. We further find that the gravity from the central supermassive black hole also affects the bimodal fate of cooling flows, and we predict a more general quenching predictor to be M{sub bh}{sup 1.6}M{sub ∗}/R{sub eff}{sup 1.5}, which may be tested in future observational studies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

The role of elastic and plastic anisotropy in intergranular spall failure

Recent mesoscale experimental observations of dynamic ductile failure have demonstrated a strong relationship between grain boundary (GB) misorientation and the likelihood of failure initiation along said GB. This correlation has been attributed to inherent GB weakness of particular misorientation. Here we discuss the role played by mechanics, i.e. elastic and plastic anisotropy, on the experimental observation. We make use of a recently developed framework for modeling dislocation-based crystal plasticity and ductile failure of single crystals under dynamic loading (CPD-FE). Polycrystals are studied at the mesoscale level through the explicit resolution of individual grains, i.e. resolving each individual grain's size, shape, and orientation. In our simulations, failure naturally localizes along the GBs with no necessity for ad hoc rules governing damage nucleation. We carry out a few thousand mesoscale calculations, systematically varying the misorientation angles of the GB in the computational microstructure. Despite the fact that we neglect the possibility of variations in inherent GB weakness, our simulations agree favorably with the experimental observations, implying that stress concentration generated by elastic and plastic anisotropies across GBs is a dominant governing factor in this phenomenon. Lastly, we find that misorientation angle is an insufficient GB descriptor to predict the likelihood of intergranular spall failure, which is better understood through the consideration of additional GB degrees of freedom.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Planck 2015 results: XVIII. Background geometry and topology of the Universe

We report that maps of cosmic microwave background (CMB) temperature and polarization from the 2015 release of Planck data provide the highestquality full-sky view of the surface of last scattering available to date. This enables us to detect possible departures from a globally isotropic cosmology. We present the first searches using CMB polarization for correlations induced by a possible non-trivial topology with a fundamental domain that intersects, or nearly intersects, the last-scattering surface (at comoving distance χ rec ), both via a direct scan for matched circular patterns at the intersections and by an optimal likelihood calculation for specific topologies. We specialize to flat spaces with cubic toroidal (T3) and slab (T1) topologies, finding that explicit searches for the latter are sensitive to other topologies with antipodal symmetry. These searches yield no detection of a compact topology with a scale below the diameter of the last-scattering surface. The limits on the radius R i of the largest sphere inscribed in the fundamental domain (at log-likelihood ratio ΔlnL > -5 relative to a simply-connected flat Planck best-fit model) are: R i > 0.97 χ rec for the T3 cubic torus; and R i > 0.56 χ rec for the T1 slab. The limit for the T3 cubic torus from the matched-circles search is numerically equivalent, R i > 0.97 χ rec at 99% confidence level from polarization data alone. We also perform a Bayesian search for an anisotropic global Bianchi VII h geometry. In the non-physical setting, where the Bianchi cosmology is decoupled from the standard cosmology, Planck temperature data favour the inclusion of a Bianchi component with a Bayes factor of at least 2.3 units of log-evidence. However, the cosmological parameters that generate this pattern are in strong disagreement with those found from CMB anisotropy data alone. Fitting the induced polarization pattern for this model to the Planck data requires an amplitude of -0.10 ± 0.04 compared to the value of + 1 if the model were to be correct. In the physically motivated setting, where the Bianchi parameters are coupled and fitted simultaneously with the standard cosmological parameters, we find no evidence for a Bianchi VII h cosmology and constrain the vorticity of such models to (ω/H) 0 < 7.6 × 10 -10 (95% CL).

79 ASTRONOMY AND ASTROPHYSICS↗