Search NASA⌕ Search

SEARCH · Search NASA

Results for “Variational Bayesian”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

The chemical homogeneity of open clusters

Determining the level of chemical homogeneity in open clusters is of fundamental importance in the study of the evolution of star-forming clouds and that of the Galactic disk. Yet limiting the initial abundance spread in clusters has been hampered by difficulties in obtaining consistent spectroscopic abundances for different stellar types. Without reference to any specific model of stellar photospheres, a model for a homogeneous cluster is that it forms a one-dimensional sequence, with any differences between members due to variations in stellar mass and observational uncertainties. I present a novel method for investigating the abundance spread in open clusters that tests this one-dimensional hypothesis at the level of observed stellar spectra, rather than constraining homogeneity using derived abundances as is traditionally done. Using high-resolution APOGEE spectra for 49 giants in M67, NGC 6819, and NGC 2420 I demonstrate that these spectra form one-dimensional sequences for each cluster. With detailed forward modeling of the spectra and Approximate Bayesian Computation, I derive strong limits on the initial abundance spread of 15 elements: <0.01 (0.02) dex for C and Fe, ≲0.015 (0.03) dex for N, O, Mg, Si, and Ni, ≲0.02 (0.03) dex for Al, Ca, and Mn, and ≲0.03 (0.05) dex for Na, S, K, Ti, and V (at 68% and 95% confidence, respectively). The strong limits on C and O imply that no pollution by massive core-collapse supernovae occurred during star formation in open clusters, which, thus, need to form within ≲6 Myr. Further development of this and related techniques will bring the power of differential abundances to stars other than solar twins in large spectroscopic surveys and will help unravel the history of star formation and chemical enrichment in the Milky Way through chemical tagging.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Decadal Variability of the Indian and Pacific Walker Cells since the 1960s: Do They Covary on Decadal Time Scales?

Previous studies have investigated the centennial and multidecadal trends of the Pacific and Indian Ocean Walker cells (WCs) during the past century, but have obtained no consensus owing to data uncertainties and weak signals of the long-term trends. This paper focuses on decadal variability (periods of one to few decades) by first documenting the variability of the WCs and warm-pool convection, and their covariability since the 1960s, using in situ and satellite observations and reanalysis products. The causes for the variability and covariability are then explored using a Bayesian dynamic linear model, which can extract nonstationary effects of climate modes. The warm-pool convection exhibits apparent decadal variability, generally covarying with the Indian and Pacific Ocean WCs during winter (November–April) with enhanced convection corresponding to intensified WCs, and the Indian–Pacific WCs covary. During summer (May–October), the warm-pool convection still highly covaries with the Pacific WC but does not covary with the Indian Ocean WC, and the Indian–Pacific WCs are uncorrelated. The wintertime coherent variability results from the vital influence of ENSO decadal variation, which reduces warm-pool convection and weakens the WCs during El Niño–like conditions. During summer, while ENSO decadal variability still dominates the Pacific WC, decadal variations of ENSO, the Indian Ocean dipole, Indian summer monsoon convection, and tropical Indian Ocean SST have comparable effects on the Indian Ocean WC overall, with monsoon convection having the largest effect since the 1990s. The complex causes for the Indian Ocean WC during summer result in its poor covariability with the Pacific WC and warm-pool convection.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

Uncertainty Analysis in 3D Equilibrium Reconstruction

Reconstruction is an inverse process where a parameter space is searched to locate a set of parameters with the highest probability of describing experimental observations. Due to systematic errors and uncertainty in experimental measurements, this optimal set of parameters will contain some associated uncertainty. This uncertainty in the optimal parameters leads to uncertainty in models derived using those parameters. V3FIT is a three-dimensional (3D) equilibrium reconstruction code that propagates uncertainty from the input signals, to the reconstructed parameters, and to the final model. Here in this paper, we describe the methods used to propagate uncertainty in V3FIT. Using the results of whole shot 3D equilibrium reconstruction of the Compact Toroidal Hybrid, this propagated uncertainty is validated against the random variation in the resulting parameters. Two different model parameterizations demonstrate how the uncertainty propagation can indicate the quality of a reconstruction. As a proxy for random sampling, the whole shot reconstruction results in a time interval that will be used to validate the propagated uncertainty from a single time slice.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Bayesian characterization and forward propagation of the uncertainty in thermodynamic models for redox active materials

Demonstrating the thermodynamic efficiency of hydrogen conversion processes using various materials is a critical step in developing new technologies for storing concentrated solar energy, and is largely accomplished by using a thermodynamic model derived from experimental data. A main goal of this project is to calculate the uncertainty of the thermodynamic efficiency by calculating the uncertainty of the components that feed into the efficiency. Many different models and data sets were used to test the workflow. First, the models were fit to the data using a Bayesian Inference and a method called Markov Chain Monte Carlo (MCMC), which found the maximum a priori parameters, and a posterior probability distribution of the parameters. Next, the different models were compared to each other using model evidence values. It was found that for cleaner data sets, overfitting had not yet been reached, and the most complicated model was ideal, but on the noisier data sets, the less complex models were favored because the more complicated models resulted in overfitting. Next, forward propagation was used to calculate the enthalpy change and its associated uncertainty. A few variations on the models were tried, such as fitting in a different variable, producing negligible or negative effects on the fits of the models. Thus, the original models were used. A sensitivity analysis was performed, and used to calculate the model error. On the cleaner data sets, there was very minimal experimental noise, and thus, all resulting error was from the model. With consideration of the model error, the models fit the data very well, and the simpler model had a high model error, as expected. All these components will then be used to calculate the thermodynamic efficiency of the different materials.

14 SOLAR ENERGY↗
DOE OSTI2019

Modeling the Formation of the Family of the Dwarf Planet Haumea

The dwarf planet (136108) Haumea has an intriguing combination of unique physical properties: near-breakup spin, two regular satellites, and an unexpectedly compact family. While these properties point toward formation by a collision, there is no self-consistent and reasonably probable formation hypothesis that can connect the unusually rapid spin and low relative velocities of Haumea family members (“Haumeans”). We explore and test the proposed formation hypotheses (catastrophic collision, graze-and-merge, and satellite collision) in detail. We flexibly parameterize the properties of the collision (e.g., the collision location) and use simple models for the unique three-dimensional velocity ejection field expected from each model to generate simulated families. These are then compared to the observed Kuiper Belt objects using Bayesian parameter inference, including a mixture model that robustly allows for interlopers from the background population. After testing our methodology, we find that the best match to the observed Haumeans is an essentially isotropic ejection field with a typical velocity of 150 m s{sup −1}. The graze-and-merge formation hypothesis—in which Haumeans are shed due to excess angular momentum—is clearly disfavored because the observed Haumeans are not oriented in a plane. The satellite collision model is also disfavored. Including these new constraints, we present a detailed discussion of the formation hypotheses, including variations, some of which are tested. Some new hypotheses are proposed (a cratering collision and a collision where Haumea’s upper layers are “missing”) and scrutinized. We do not identify a satisfactory formation hypothesis, but we do propose several avenues of additional investigation. In the process of these analyses, we identify many new candidate Haumeans and dynamically confirm seven of them as consistent with the observed family. We also confirm that Haumeans have a shallow size distribution and discuss implications for the discovery and identification of new Haumeans.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

BEYONDPLANCK II. CMB mapmaking through Gibbs sampling

We present a Gibbs sampling solution to the mapmaking problem for cosmic microwave background (CMB) measurements that builds on existing destriping methodology. Gibbs sampling breaks the computationally heavy destriping problem into two separate steps: noise filtering and map binning. Considered as two separate steps, both are computationally much cheaper than solving the combined problem. This provides a huge performance benefit as compared to traditional methods and it allows us, for the first time, to bring the destriping baseline length to a single sample. Here, we applied the Gibbs procedure to simulated Planck 30 GHz data. We find that gaps in the time-ordered data are handled efficiently by filling them in with simulated noise as part of the Gibbs process. The Gibbs procedure yields a chain of map samples, from which we are able to compute the posterior mean as a best-estimate map. The variation in the chain provides information on the correlated residual noise, without the need to construct a full noise covariance matrix. However, if only a single maximum-likelihood frequency map estimate is required, we find that traditional conjugate gradient solvers converge much faster than a Gibbs sampler in terms of the total number of iterations. The conceptual advantages of the Gibbs sampling approach lies in statistically well-defined error propagation and systematic error correction. This methodology thus forms the conceptual basis for the mapmaking algorithm employed in the BEYONDPLANCK framework, which implements the first end-to-end Bayesian analysis pipeline for CMB observations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Parkes Pulsar Timing Array constraints on ultralight scalar-field dark matter

It is widely accepted that dark matter contributes about a quarter of the critical mass-energy density in our Universe. The nature of dark matter is currently unknown, with the mass of possible constituents spanning nearly one hundred orders of magnitude. The ultralight scalar field dark matter, consisting of extremely light bosons with m~10 –22 eV and often called “fuzzy” dark matter, provides intriguing solutions to some challenges at sub-Galactic scales for the standard cold dark matter model. As shown by Khmelnitsky and Rubakov, such a scalar field in the Galaxy would produce an oscillating gravitational potential with nanohertz frequencies, resulting in periodic variations in the times of arrival of radio pulses from pulsars. The Parkes Pulsar Timing Array (PPTA) has been monitoring 20 millisecond pulsars at two- to three-week intervals for more than a decade. In addition to the detection of nanohertz gravitational waves, PPTA offers the opportunity for direct searches for fuzzy dark matter in an astrophysically feasible range of masses. We analyze the latest PPTA data set which includes timing observations for 26 pulsars made between 2004 and 2016. Here, we perform a search in this data set for evidence of ultralight dark matter in the Galaxy using Bayesian and Frequentist methods. No statistically significant detection has been made. We, therefore, place upper limits on the local dark matter density. Our limits, improving on previous searches by a factor of 2 to 5, constrain the dark matter density of ultralight bosons with m ≤ 10 –23 eV to be below 6 GeV cm –3 with 95% confidence in the Earth neighborhood. Finally, we discuss the prospect of probing the astrophysically favored mass range m ≳ 10 –22 eV with next-generation pulsar timing facilities.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Can we use antipredator behavior theory to predict wildlife responses to high-speed vehicles?

Animals seem to rely on antipredator behavior to avoid vehicle collisions. There is an extensive body of antipredator behavior theory that have been used to predict the distance/time animals should escape from predators. These models have also been used to guide empirical research on escape behavior from vehicles. However, little is known as to whether antipredator behavior models are appropriate to apply to an approaching high-speed vehicle scenario. We addressed this gap by (a) providing an overview of the main hypotheses and predictions of different antipredator behavior models via a literature review, (b) exploring whether these models can generate quantitative predictions on escape distance when parameterized with empirical data from the literature, and (c) evaluating their sensitivity to vehicle approach speed using a simulation approach wherein we assessed model performance based on changes in effect size with variations in the slope of the flight initiation distance (FID) vs. approach speed relationship. The slope of the FID vs. approach speed relationship was then related back to three different behavioral rules animals may rely on to avoid approaching threats: the spatial, temporal, or delayed margin of safety. We used literature on birds for goals (b) and (c). Our review considered the following eight models: the economic escape model, Blumstein’s economic escape model, the optimal escape model, the perceptual limit hypothesis, the visual cue model, the flush early and avoid the rush (FEAR) hypothesis, the looming stimulus hypothesis, and the Bayesian model of escape behavior. We were able to generate quantitative predictions about escape distance with the last five models. However, we were only able to assess sensitivity to vehicle approach speed for the last three models. The FEAR hypothesis is most sensitive to high-speed vehicles when the species follows the spatial (FID remains constant as speed increases) and the temporal margin of safety (FID increases with an increase in speed) rules of escape. The looming stimulus effect hypothesis reached small to intermediate levels of sensitivity to high-speed vehicles when a species follows the delayed margin of safety (FID decreases with an increase in speed). The Bayesian optimal escape model reached intermediate levels of sensitivity to approach speed across all escape rules (spatial, temporal, delayed margins of safety) but only for larger (> 1 kg) species, but was not sensitive to speed for smaller species. Overall, no single antipredator behavior model could characterize all different types of escape responses relative to vehicle approach speed but some models showed some levels of sensitivity for certain rules of escape behavior. We derive some applied applications of our findings by suggesting the estimation of critical vehicle approach speeds for managing populations that are especially susceptible to road mortality. Overall, we recommend that new escape behavior models specifically tailored to high-speeds vehicles should be developed to better predict quantitatively the responses of animals to an increase in the frequency of cars, airplanes, drones, etc. they will face in the next decade.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Precision mass measurement of lightweight self-conjugate nucleus 80 Zr

Protons and neutrons in the atomic nucleus move in shells analogous to the electronic shell structures of atoms. The nuclear shell structure varies due to changes of the nuclear mean field with the number of neutrons N and protons Z, and these variations can be probed by measuring mass differences between nuclei. The N=Z=40 self-conjugate nucleus 80 Zr is of particular interest as its proton and neutron shell structures are expected to be very similar, and its ground state is highly deformed. Here, we provide evidence for the existence of a deformed double shell closure in 80Zr through high precision Penning trap mass measurements of 80-83 Zr. Our mass values show that 80Zr is significantly lighter, and thus more strongly bound than predicted. This can be attributed to the deformed shell closure at N=Z=40 and the large Wigner energy. A statistical Bayesian model mixing analysis employing several global nuclear mass models demonstrates difficulties with reproducing the observed mass anomaly using current theory.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2016

Gaussian processes with built-in dimensionality reduction: Applications to high-dimensional uncertainty propagation

Uncertainty quantification (UQ) tasks, such as model calibration, uncertainty propagation, and optimization under uncertainty, typically require several thousand evaluations of the underlying computer codes. To cope with the cost of simulations, one replaces the real response surface with a cheap surrogate based, e.g., on polynomial chaos expansions, neural networks, support vector machines, or Gaussian processes (GP). However, the number of simulations required to learn a generic multivariate response grows exponentially as the input dimension increases. This curse of dimensionality can only be addressed, if the response exhibits some special structure that can be discovered and exploited. A wide range of physical responses exhibit a special structure known as an active subspace (AS). An AS is a linear manifold of the stochastic space characterized by maximal response variation. The idea is that one should first identify this low dimensional manifold, project the high-dimensional input onto it, and then link the projection to the output. If the dimensionality of the AS is low enough, then learning the link function is a much easier problem than the original problem of learning a high-dimensional function. The classic approach to discovering the AS requires gradient information, a fact that severely limits its applicability. Furthermore, and partly because of its reliance to gradients, it is not able to handle noisy observations. The latter is an essential trait if one wants to be able to propagate uncertainty through stochastic simulators, e.g., through molecular dynamics codes. In this work, we develop a probabilistic version of AS which is gradient-free and robust to observational noise. Our approach relies on a novel Gaussian process regression with built-in dimensionality reduction. In particular, the AS is represented as an orthogonal projection matrix that serves as yet another covariance function hyper-parameter to be estimated from the data. To train the model, we design a two-step maximum likelihood optimization procedure that ensures the orthogonality of the projection matrix by exploiting recent results on the Stiefel manifold, i.e., the manifold of matrices with orthogonal columns. The additional benefit of our probabilistic formulation, is that it allows us to select the dimensionality of the AS via the Bayesian information criterion. We validate our approach by showing that it can discover the right AS in synthetic examples without gradient information using both noiseless and noisy observations. We demonstrate that our method is able to discover the same AS as the classical approach in a challenging one-hundred-dimensional problem involving an elliptic stochastic partial differential equation with random conductivity. Finally, we use our approach to study the effect of geometric and material uncertainties in the propagation of solitary waves in a one dimensional granular system.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

BEYONDPLANCK VIII. Efficient sidelobe convolution and corrections through spin harmonics

We introduce a new formulation of the Conviqt convolution algorithm in terms of spin harmonics, and apply this to the problem of sidelobe correction for BEYONDPLANCK, the first end-to-end Bayesian Gibbs sampling framework for CMB analysis. We compare our implementation to the previous Planck LevelS implementation, and find good agreement between the two codes in terms of accuracy, but with a speed-up reaching a factor of 3–10, depending on the frequency bandlimits, l max and m max . The new algorithm is significantly simpler to implement and maintain, since all low-level calculations are handled through an external spherical harmonic transform library. We find that our mean sidelobe estimates for Planck LFI are in good agreement with previous efforts. Additionally, we present novel sidelobe rms maps that quantify the uncertainty in the sidelobe corrections due to variations in the sky model.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Phylogenetic relationships and phylogeography of relevant lineages within the complex Campanulaceae family in Macaronesia

Abstract Macaronesia has long been recognized as a natural model for studying evolutionary processes in plant diversification. Several studies have attempted to focus on single lineages, and few have covered the diversification of a family across all the archipelagos. We used a comprehensive sample to clarify the phylogenetic relationships and the biogeographic history of the Macaronesian Campanulaceae. Hypotheses related to the colonization of these archipelagos will be used to examine the diversification patterns of different lineages. We sequenced the ITS region and six cp DNA markers ( atpB , matK , petD , rbcL , trnL‐F , and psbA‐trnH ) from 10 Campanulaceae species, including seven endemic species in Macaronesia. The phylogeny of these taxa was reconstructed using maximum parsimony, maximum likelihood, and Bayesian inference. To study the relationships within each lineage, haplotype networks were calculated using NeighborNet and TCS algorithms. Moreover, data were combined with fossil information to construct time‐calibrated trees for the Macaronesian Campanulaceae species. The phylogenetic analyses are largely congruent with current taxon circumscriptions, and all the endemic genera formed monophyletic clades, namely Azorina in Azores; Musschia in Madeira; and Campanula in Cape Verde. The Azorina clade and the Cape Verde endemic Campanula may share a common ancestor in North Africa, and the divergence was dated ca. 12.3 million years ago (Mya). The divergence of the Musschia clade began in the Pliocene ca. 3.4 Mya. Moreover, several examples of intraspecific variation were revealed among the native species with a clear geographic structured patterns, suggesting that cryptic diversity might exist within the native Macaronesian Campanulaceae when compared to the close mainland taxa (e.g., Campanula erinus , Trachelium caeruleum ), but additional studies are needed to support the molecular data. This study highlights the power of combining data (e.g., phylogeny and divergence times, with species distribution data) for testing diversification hypotheses within the unique Macaronesian flora, providing useful information for future conservation efforts.

Menezes, Tiago↗
DOE OSTI2022

A mixed, unified forward/inverse framework for earthquake problems: fault implementation and coseismic slip estimate

SUMMARY We introduce a new finite-element (FE) based computational framework to solve forward and inverse elastic deformation problems for earthquake faulting via the adjoint method. Based on two advanced computational libraries, FEniCS and hIPPYlib for the forward and inverse problems, respectively, this framework is flexible, transparent and easily extensible. We represent a fault discontinuity through a mixed FE elasticity formulation, which approximates the stress with higher order accuracy and exposes the prescribed slip explicitly in the variational form without using conventional split node and decomposition discrete approaches. This also allows the first order optimality condition, that is the vanishing of the gradient, to be expressed in continuous form, which leads to consistent discretizations of all field variables, including the slip. We show comparisons with the standard, pure displacement formulation and a model containing an in-plane mode II crack, whose slip is prescribed via the split node technique. We demonstrate the potential of this new computational framework by performing a linear coseismic slip inversion through adjoint-based optimization methods, without requiring computation of elastic Green’s functions. Specifically, we consider a penalized least squares formulation, which in a Bayesian setting—under the assumption of Gaussian noise and prior—reflects the negative log of the posterior distribution. The comparison of the inversion results with a standard, linear inverse theory approach based on Okada’s solutions shows analogous results. Preliminary uncertainties are estimated via eigenvalue analysis of the Hessian of the penalized least squares objective function. Our implementation is fully open-source and Jupyter notebooks to reproduce our results are provided. The extension to a fully Bayesian framework for detailed uncertainty quantification and non-linear inversions, including for heterogeneous media earthquake problems, will be analysed in a forthcoming paper.

58 GEOSCIENCES↗
DOE OSTI2021

A Scalable Gradient Free Method for Bayesian Experimental Design with Implicit Models

Bayesian experimental design (BED) is to answer the question that how to choose designs that maximize the information gathering. For implicit models, where the likelihood is intractable but sampling is possible, conventional BED methods have difficulties in efficiently estimating the posterior distribution and maximizing the mutual information (MI) between data and parameters. Recent work proposed the use of gradient ascent to maximize a lower bound on MI to deal with these issues. However, the approach requires a sampling path to compute the pathwise gradient of the MI lower bound with respect to the design variables, and such a pathwise gradient is usually inaccessible for implicit models. In this paper, we propose a novel approach that leverages recent advances in stochastic approximate gradient ascent incorporated with a smoothed variational MI estimator for efficient and robust BED. Without the necessity of pathwise gradients, our approach allows the design process to be achieved through a unified procedure with an approximate gradient for implicit models. Several experiments show that our approach outperforms baseline methods, and significantly improves the scalability of BED in high-dimensional problems.

Zhang, Jiaxin↗
DOE OSTI2021

A fast multi-fidelity method with uncertainty quantification for complex data correlations: Application to vortex-induced vibrations of marine risers

Here we develop a fast multi-fidelity modeling method for very complex correlations between high- and low-fidelity data by working in modal space to extract the proper correlation function. We apply this method to infer the amplitude of motion of a flexible marine riser in cross-flow, subject to vortex-induced vibrations (VIV). VIV are driven by an absolute instability in the flow, which imposes a frequency (Strouhal) law that requires a matching with the impedance of the structure; this matching is easily achieved because of the rapid parametric variation of the added mass force. As a result, the wavenumber of the riser spatial response is within narrow bands of uncertainty. Hence, an error in wavenumber prediction can cause significant phase-related errors in the shape of the amplitude of response along the riser, rendering correlation between low- and high-fidelity data very complex. Working in modal space as outlined herein, dense data from low-fidelity data, provided by the semi-empirical computer code VIVA, can correlate in modal space with few high-fidelity data, obtained from experiments or fully-resolved CFD simulations, to correct both phase and amplitude and provide predictions that agree very well overall with the correct shape of the amplitude response. We also quantify the uncertainty in the prediction using Bayesian modeling and exploit this uncertainty to formulate an active learning strategy for the best possible location of the sensors providing the high fidelity measurements.

42 ENGINEERING↗
DOE OSTI2020

A Hybrid Gradient Method to Designing Bayesian Experiments for Implicit Models

Bayesian experimental design (BED) aims at designing an experiment to maximize the information gathering from the collected data. The optimal design is usually achieved by maximizing the mutual information (MI) between the data and the model parameters. When the analytical expression of the MI is unavailable, e.g.,having implicit models with intractable data distributions, a neural network-based lower bound of the MI was recently proposed and a gradient ascent method was used to maximize the lower bound [1]. However, the approach in [1] requires a pathwise sampling path to compute the gradient of the MI lower bound with respect to the design variables, and such a pathwise sampling path is usually inaccessible for implicit models. In this work, we propose a hybrid gradient approach that leverages recent advances in variational MI estimator and evolution strategies (ES)combined with black-box stochastic gradient ascent (SGA) to maximize the MI lower bound. This allows the design process to be achieved through a unified scalable procedure for implicit models without sampling path gradients. Several experiments demonstrate that our approach significantly improves the scalability of BED for implicit models in high-dimensional design space.

Zhang, Jiaxin↗
DOE OSTI2019

Ecological memory of daily carbon exchange across the globe and its importance in drylands

How do antecedent (past) conditions influence land-carbon dynamics after those conditions no longer persist? In particular, quantifying such memory effects associated with the influence of past environmental (exogenous) and biological (endogenous) conditions is crucial for understanding and predicting the carbon cycle. Here we show, using data from 42 eddy covariance sites across six major biomes, that ecological memory—decomposed into environmental and biological memory components—of daily net carbon exchange (NEE) is critical for understanding the land-carbon metabolism, especially in drylands for which memory explains ~ 32% of the variation in NEE. The strong environmental memory in drylands was primarily driven by short- and long-term moisture status. Moreover, the strength of environmental memory scales with increasing water stress. This universal scaling relationship, emerging within and among major biomes, suggests a potential adaptive response to water limitation. Our findings underscore the necessity of considering ecological memory in experiments, observations and modelling.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Rapid assessment of the creep rupture life of metals: A model enabling experimental design

Prediction of the creep rupture life of engineering metals is critical for qualification and design of new materials. The use of long-term creep tests and the need to quantify the performance variability in a priori similar systems hinder the rapid creep assessment of a given material. Therefore, it is essential to develop methods that can extrapolate the long-term performance of alloys and the associated variability from short-term experiments. To this end, this study introduces a new model which enables the estimation of the rupture life of a material for a given stress and temperature. This model relies on two components. First, a new relation for the minimum creep rate (MCR) of materials is introduced. It includes a stress dependent stress exponent allowing the model to capture the variation of MCR across a wide range of temperatures and stresses. Second, employing the Monkman-Grant (MG) law, we establish a relation between stress, temperature and creep rupture life. Together, these two elements yield a new closed-form mathematical expression for the Larson Miller parameter as a function of stress and temperature. This expression captures the creep rupture time for many metals (Gr91, Copper, Gr122 and 347H) and compares favorably with alternate empirical approaches. The model is then used to assess the minimum duration of creep rates necessary to qualify the material up to 100000h. Furthermore, it is found that depending on the material system, creep tests as few as five limited to 5000 h for steels (Gr91, Gr122, 347H) and 100 h for copper are sufficient to model creep lifetimes. Finally, using a Bayesian inference-based approach to calibrate the model, we demonstrate that variability in rupture life can be captured via the quantification of the uncertainty in the model parameters and extrapolated from a limited number of short to moderately short creep tests; thereby paving the way for accelerated creep testing.

36 MATERIALS SCIENCE↗