Search NASA⌕ Search

SEARCH · Search NASA

Results for “Function Estimation”

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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

J-PLUS: Spectral evolution of white dwarfs by PDF analysis

We estimated the spectral evolution of white dwarfs with effective temperature using the Javalambre Photometric Local Universe Survey (J-PLUS) second data release (DR2), which provides 12 photometric optical passbands over 2176 deg 2 . We analyzed 5926 white dwarfs with r ≤ 19.5 mag in common between a white dwarf catalog defined from Gaia EDR3 and J-PLUS DR2. We performed a Bayesian analysis by comparing the observed J-PLUS photometry with theoretical models of hydrogen- and helium-dominated atmospheres. We estimated the probability distribution functions for effective temperature (T eff ), surface gravity, parallax, and composition; and the probability of having a H-dominated atmosphere (p H ) for each source. We applied a prior in parallax, using Gaia EDR3 measurements as a reference, and derived a self-consistent prior for the atmospheric composition as a function of T eff . We described the fraction of white dwarfs with a He-dominated atmosphere (f He ) with a linear function of the effective temperature at 5000 < T eff < 30 000 K. We find f He = 0.24 ± 0.01 at T eff = 10 000 K, a change rate along the cooling sequence of 0.14 ± 0.02 per 10 kK, and a minimum He-dominated fraction of 0.08 ± 0.02 at the high-temperature end. We tested the obtained p H by comparison with spectroscopic classifications, finding that it is reliable. We estimated the mass distribution for the 351 sources with distance d < 100 pc, mass M > 0.45 M ⊙ , and T eff > 6000 K. The result for H-dominated white dwarfs agrees with previous studies, with a dominant M = 0.59 M ⊙ peak and the presence of an excess at M ~ 0.8 M ⊙ . This high-mass excess is absent in the He-dominated distribution, which presents a single peak. The J-PLUS optical data provide a reliable statistical classification of white dwarfs into H- and He-dominated atmospheres. We find a 21 ± 3% increase in the fraction of He-dominated white dwarfs from T eff = 20 000 K to T eff = 5000 K.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

A combined first principles study of the structural, magnetic, and phonon properties of monolayer CrI3

The first magnetic 2D material discovered, monolayer (ML) CrI3, is particularly fascinating due to its ground state ferromagnetism. However, because ML materials are difficult to probe experimentally, much remains unresolved about ML CrI3’s structural, electronic, and magnetic properties. Here, we leverage Density Functional Theory (DFT) and high-accuracy Diffusion Monte Carlo (DMC) simulations to predict lattice parameters, magnetic moments, and spin–phonon and spin–lattice coupling of ML CrI3. We exploit a recently developed surrogate Hessian DMC line search technique to determine CrI3’s ML geometry with DMC accuracy, yielding lattice parameters in good agreement with recently published STM measurements—an accomplishment given the ~10% variability in previous DFT-derived estimates depending upon the functional. Strikingly, we find that previous DFT predictions of ML CrI3’s magnetic spin moments are correct on average across a unit cell but miss critical local spatial fluctuations in the spin density revealed by more accurate DMC. DMC predicts that magnetic moments in ML CrI3 are 3.62 μB per chromium and -0.145 μB per iodine, both larger than previous DFT predictions. The large disparate moments together with the large spin–orbit coupling of CrI3’s I-p orbital suggest a ligand superexchange-dominated magnetic anisotropy in ML CrI3, corroborating recent observations of magnons in its 2D limit. We also find that ML CrI3 exhibits a substantial spin–phonon coupling of ~3.32 cm-1. Our work, thus, establishes many of ML CrI3’s key properties, while also continuing to demonstrate the pivotal role that DMC can assume in the study of magnetic and other 2D materials.

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

Atomistic simulations and machine learning of solute grain boundary segregation in Mg alloys at finite temperatures

Understanding solute segregation thermodynamics is the first step in investigating grain boundary (GB) properties, such as strong yttrium (Y) effects on grain growth and texture evolution in micro-scale polycrystalline magnesium (Mg) alloys. To estimate the average GB segregation behavior in low-solute-concentration Mg alloys (e.g., 2 at.% Y), a state-of-the-art spectral approach is applied based on a per-site segregation energy spectrum for Y solute atoms at zero K obtained from molecular statistics (MS) simulations of ~10 4 GB sites in Mg symmetric tilt GBs (STGBs). Although selected MS simulation results are consistent with verification by density functional theory (DFT) calculations, estimates of average segregation tendency based on the zero-K energy spectrum deviate from experimental observations. To resolve this problem, thermodynamic integration (TI) methods based on molecular dynamics (MD) simulations are used to determine the per-site segregation free energies of Y at representative GB sites, which show contributions beyond harmonic approximations can be important for certain GB sites at high temperatures. A surrogate model of per-site segregation free energy is constructed from a small set of TI data points using stacking cross-validation regressors and physics-informed descriptors. This model is applied to predict the Y segregation free energy spectra for thousands of GB sites in Mg STGBs with uncertainty quantification. Finally, the average segregation tendency predicted by the spectral approach based on the free energy spectra agrees well (within the uncertainty range) with experimental observations for micro-scale polycrystalline Mg alloys at typical thermomechanical processing temperatures (500 ~ 800 K), where thermodynamic equilibrium states are likely to be achieved due to fast diffusion.

Atomistic simulations↗
DOE OSTI2024

Evaluation of average leaf inclination angle quantified by indirect optical instruments in crop fields

Average leaf inclination angle ($\overline{θ}$ L ) is an important canopy structure variable that influences light regime, photosynthesis, and evapotranspiration of plants. $\overline{θ}$ L can be measured through direct methods (e.g., protractor), which are labor-intensive and time-consuming, or through indirect optical instruments, which are more efficient than the direct methods. However, uncertainties of different indirect optical instruments for quantifying $\overline{θ}$ L remain largely unquantified. In this study, we evaluated and compared the performances of three major indirect optical instruments: (1) LAI-2200, (2) 30°-tilted camera, and (3) digital hemispherical photography (DHP), in different crop fields over a growing season, benchmarked with direct measurements. LAI-2200 and 30°-tilted camera showed higher agreement with direct $\overline{θ}$ measurements (R 2 = 0.54, RMSE = 7.37°; R 2 = 0.58, RMSE = 8.08°) than DHP (R 2 = 0.14, RMSE = 13.96°). Different performances of indirect optical instruments could be attributed to the accuracy of gap fraction measurement and the performance of the $\overline{θ}$ L quantification algorithms. When using the LAI-2200 algorithm, larger gap fraction gradients over view zenith angles led to larger $\overline{θ}$ L values, and smaller gap fraction gradients led to smaller $\overline{θ}$ L values. Such error propagation was larger in sparse canopy than in dense canopy. The Wilson G function of the LAI-2200 algorithm performed better in estimating $\overline{θ}$ L than the G function based on the ellipsoidal LAD function used by the CAN_EYE algorithm. We also proposed a modification of the LAI-2200 algorithm, which further improved the performance of LAI-2200 and 30°-tilted cameras in estimating $\overline{θ}$ L . We envision that the low-cost 30°-tilted cameras provide a promising sensor solution to continuously monitor canopy structure for various ecosystems.

30°-tilted camera↗
DOE OSTI2023

Size-Dependent Energy and Adhesion of Pd Nanoparticles on Graphene on Ni(111) by Pd Vapor Adsorption Calorimetry

Carbon-supported late transition-metal nanoparticles are promising catalysts and electrocatalysts for wide-ranging applications. However, experimental investigations of the bonding energetics of metal nanoparticles on carbon supports are very limited. Here, in this work, we report heats of adsorption of Pd vapor deposited onto single-layer graphene(0001) supported on Ni(111) at 100 and 300 K as Pd grows particles of well-defined size in the range from three atom clusters to 6 nm diameter. Sizes were determined from He + low-energy ion scattering (LEIS). In this size range, the differential heat of Pd adsorption increases from 228 kJ/mol to within 10 kJ/mol of the heat of sublimation of bulk Pd (377 kJ/mol). The chemical potential of metal atoms in these nanoparticles as a function of average particle size was determined from these results. The adhesion energy at the Pd/graphene(0001)/Ni(111) interface was extracted from these data and found to be 3.5 J/m 2 for the largest Pd particles. For the three metal elements that have now been studied (Pd, Ni, and Ag), we found an increase in metal/graphene(0001)/Ni(111) adhesion energy with metal carbophilicity, which we defined here as the heat of C atom adsorption on that metal’s (111) surface estimated from published density functional theory calculations.

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

A new computational approach for the solutions of generalized pantograph-delay differential equations

In this study, a new computational approach is presented to solve the generalized pantograph-delay differential equations (PDDEs). The solutions obtained by our scheme represent by a linear combination of a special kind of basis functions, and can be deduced in a straightforward manner. Firstly, using the least squares approximation method and the Lagrange-multiplier method, the given PDDE is converted to a linear system of algebraic equations, and those unknown coefficients of the solution of the problem are determined by solving this linear system. Secondly, a PDDE related to the error function of the approximate solution is constructed based on the residual error function technique, and error estimation is presented for the suggested method. The convergence of the approximate solution is proved. Several numerical examples are given to demonstrate the accuracy and efficiency. Comparisons are made between the proposed method and other existing methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Photometric redshifts probability density estimation from recurrent neural networks in the DECam local volume exploration survey data release 2

Photometric wide-field surveys are imaging the sky in unprecedented detail. These surveys face a significant challenge in efficiently estimating galactic photometric redshifts while accurately quantifying associated uncertainties. In this work, we address this challenge by exploring the estimation of Probability Density Functions (PDFs) for the photometric redshifts of galaxies across a vast area of 17,000 square degrees, encompassing objects with a median 5 σ point-source depth of g = 24.3, r = 23 . 9 , i = 23.5, and z = 22.8 mag. Our approach uses deep learning, specifically integrating a Recurrent Neural Network architecture with a Mixture Density Network, to leverage magnitudes and colors as input features for constructing photometric redshift PDFs across the whole DECam Local Volume Exploration (DELVE) survey sky footprint. Subsequently, we rigorously evaluate the reliability and robustness of our estimation methodology, gauging its performance against other well-established machine learning methods to ensure the quality of our redshift estimations. Our best results constrain photometric redshifts with the bias of − 0 . 0013 , a scatter of 0.0293, and an outlier fraction of 5.1%. These point estimates are accompanied by well-calibrated PDFs evaluated using diagnostic tools such as Probability Integral Transform and Odds distribution. We also address the problem of the accessibility of PDFs in terms of disk space storage and the time demand required to generate their corresponding parameters.We present a novel Autoencoder model that reduces the size of PDF parameter arrays to one-sixth of their original length, significantly decreasing the time required for PDF generation to one-eighth of the time needed when generating PDFs directly from the magnitudes.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Data-Centric Approach to Capture Non-Polynomial Nonlinear Dynamics

We propose an analytical construction of observable functions in the extended dynamic mode decomposition (EDMD) algorithm. EDMD is a numerical method for approximating the spectral properties of the Koopman operator. The choice of observable functions is fundamental for applying EDMD to nonlinear problems arising in systems and control. Existing methods either start from a set of dictionary functions and look for the subset that best fits the underlying nonlinear dynamics or rely on machine learning algorithms to “learn” observable functions. Conversely, in this paper, we start from the dynamical system model and lift it through the Lie derivatives, rendering it into a polynomial form. This proposed transformation into a polynomial form is exact and provides an adequate set of observable functions. The strength of the proposed approach is its applicability to a broader class of nonlinear dynamical systems, particularly those with nonpolynomial functions and compositions thereof. Moreover, it retains the physical interpretability of the underlying dynamical system and can be readily integrated into existing numerical libraries. We demonstrate the proposed approach with an application to electric power systems. The modeled system consists of a single generator connected to an infinite bus, where nonlinear terms include sine and cosine functions. The results demonstrate the effectiveness of the proposed procedure in off-attractor nonlinear dynamics for estimation and prediction; the observable functions obtained from the proposed construction outperform methods that use dictionary functions comprising monomials or radial basis functions.

extended dynamic mode decomposition↗
DOE OSTI2019

Satellite-based monitoring of groundwater depletion in California’s Central Valley

Abstract Range change data, obtained from Synthetic Aperture Radar satellites, form the basis for estimates of aquifer volume change in California’s Central Valley. The estimation algorithm incorporates a function penalizing changes far from known well locations, linking the aquifer volume changes to agricultural, industrial, and municipal pumping within the Tulare basin. We show that the range changes are compatible with the hypothesis that the source of aquifer volume changes are variations in effective pressure around documented wells. Specifically, inclusion of the well distance penalty does not degrade the fit to the observations, inversions with and without it both give variance reductions of 99.6%. The patterns of aquifer volume change vary significantly from the drought year, between October 2015 and October 2016, to a wet year in 2017, and into 2018, a year with near average rainfall. The 2.3 million acre-feet of estimated volume reduction, a lower bound on the amount of water extracted from the basin between October 2015 and 2016, agrees with independent estimates of 1.8 and 2.3 million acre-feet. The aquifer volume reduction is also compatible with a loss of 3.1 km 3 (2.5 million acre-feet) in groundwater volume derived from Gravity Recovery and Climate Experiment (GRACE) satellite data.

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

Density estimation via measure transport: Outlook for applications in the biological sciences

Abstract One among several advantages of measure transport methods is that they allow or a unified framework for processing and analysis of data distributed according to a wide class of probability measures. Within this context, we present results from computational studies aimed at assessing the potential of measure transport techniques, specifically, the use of triangular transport maps, as part of a workflow intended to support research in the biological sciences. Scenarios characterized by the availability of limited amount of sample data, which are common in domains such as radiation biology, are of particular interest. We find that when estimating a distribution density function given limited amount of sample data, adaptive transport maps are advantageous. In particular, statistics gathered from computing series of adaptive transport maps, trained on a series of randomly chosen subsets of the set of available data samples, leads to uncovering information hidden in the data. As a result, in the radiation biology application considered here, this approach provides a tool for generating hypotheses about gene relationships and their dynamics under radiation exposure.

gene expression data↗
DOE OSTI2022

Catalyst design to direct high-octane gasoline fuel properties for improved engine efficiency

The paraffin-to-olefin (P/O) ratio in gasoline fuel is a critical metric affecting fuel properties and engine efficiency. In the conversion of dimethyl ether (DME) to high-octane hydrocarbons over BEA zeolite catalysts, the P/O ratio can be controlled through catalyst design. Here, we report bimetallic catalysts that balance the net hydrogenation and dehydrogenation activity during DME homologation. The Cu-Zn/BEA catalyst exhibited greater relative dehydrogenation activity attributed to higher ionic site density, resulting in a lower P/O ratio (6.6) versus the benchmark Cu/BEA (9.4). The Cu-Ni/BEA catalyst exhibited increased hydrogenation due to reduced Ni species, resulting in a higher P/O ratio (19). The product fuel properties were estimated with an efficiency merit function and compared against finished gasolines and a typical alkylate blendstock. Merit values for the hydrocarbon product from all three BEA catalysts exceeded those of the comparison fuels (0–5.3), with the product from Cu-Zn/BEA exhibiting the highest merit value (9.7).

Catalyst design↗
DOE OSTI2022

Cross Inference of Throughput Profiles Using Micro Kernel Network Method

Dedicated network connections are being increasingly deployed in cloud, centralized and edge computing and data infrastructures, whose throughput profiles are critical indicators of the underlying data transfer performance. Due to the cost and disruptions to physical infrastructures, network emulators, such as Mininet, are often used to generate measurements needed to estimate throughput profiles, typically expressed as a function of the connection round trip time. The profiles estimated using measurements from such emulated networks are usually inaccurate for high bandwidth and high latency connections, since they do not accurately reflect the critical network transport dynamics mainly due to computing and memory constraints of the host. We present a machine learning (ML) method to estimate the throughput profiles using emulation measurements to closely match the testbed and production network profiles. In particular, we propose a micro Kernel Network (mKN) that provides baseline throughput measurements on the host running Mininet emulations, which are used to learn a regression map that converts them to the corresponding testbed measurement estimates. Once initially learned, this map is applied to measurements from subsequent network emulations on the same host. We present experimental measurements to illustrate this approach, and derive generalization equations for the proposed mKN-ML method. Using a four-site scenario emulation, we show the effectiveness of this method in providing accurate concave throughput profiles from inaccurate convex or non-smooth ones indicated by Mininet emulation.

Rao, Nageswara↗
DOE OSTI2023

Residual estimation for grid modification in wall-modeled large eddy simulation using unstructured high-order methods

Here, the accuracy and computational cost of a large eddy simulation are highly dependent on the computational grid. Building optimal grids manually from a priori knowledge is not feasible in most practical use cases; instead, solution-adaptive strategies can provide a robust and cost-efficient method to generate a grid with the desired accuracy. We adapt the residual estimation algorithm developed by Toosi and Larsson for Discontinuous Galerkin Spectral Elements Methods (DGSEM) to guide the grid-adaptation process. The core of the method is the computation of the estimated modeling residual using the polynomial basis functions used in DGSEM, and the averaging of the estimated residual over each element. The final method is assessed in multiple channel flow test cases and for the transonic flow over an airfoil, in both cases making use of mortar interfaces between elements with hanging nodes. The method is found to be robust and reliable, and to provide solutions on grids with significantly fewer elements at comparable accuracy compared to when using human-generated grids.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Sequential Bayesian Methods for Analyzing Computer Models

Efficient analysis of computer models is essential for the validation and uncertainty quantification of those models. Surrogate models, and Gaussian processes in particular, are a common and powerful approach to analyzing computer models that treat computer models as a black-box function. Gaussian processes form a Bayesian model over a space of functions that gives a measure of uncertainty about the computer model output at unobserved locations and a framework for sequential sampling. This tutorial will show how to fit Gaussian processes on a series of test functions and apply sequential design techniques to estimate extrema, level sets, and reliabilities of those test functions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Adaptive Reweighted Variance Estimation for Monte Carlo Eigenvalue Simulations

Monte Carlo (MC) simulation is used to solve the eigenvalue form of the Boltzmann transport equation to estimate various parameters such as fuel pin flux distributions that are crucial for the safe and efficient operation of nuclear systems (e.g., a nuclear reactor). Monte Carlo eigenvalue simulation uses a sample mean over many stationary cycles (iterations) to estimate various parameters important to nuclear systems. A variance estimate of the sample mean is often used for calculating the confidence intervals. However, MC eigenvalue simulation variance estimators that ignore the intercycle correlation underestimate the true variance of the estimated quantity. This paper presents novel data-adaptive approaches based on a simple autoregressive (AR) model and sigmoid functions to improve MC variance estimation. The standard MC sample-based variance estimator (or naïve estimator) and the spectral density–based MC variance estimator are enhanced by adding data-adaptive components that reduce their bias and improve performance. By investigating the frequency pattern of the AR(1) (order 1) model, two adaptive spectral estimators and one adaptive naïve estimator are proposed. The proposed estimators manifest superior performance when applied to three test problems compared to the standard spectral density–based estimator previously introduced by the authors. In conclusion, these new estimators are straightforward, as they use online algorithms and do not require storage of tallies from all active cycles.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Pragmatic Uncertainty Quantification and Propagation in Inverse Estimation of Structural Dynamics Parameters given Material Property Uncertainties and Limited Sensor Data

In this report we demonstrate some relatively simple and inexpensive methods to effectively account for various sources of epistemic lack-of-knowledge type uncertainty in inverse problems. The demonstration problem involves inverse estimation of six parameters of a bolted joint that attaches a kettlebell shaped object to a thick plate. The parameters are efficiently inverted in a modal-based model calibration using gradient-based optimization. Two material properties of the kettlebell are treated as uncertain to within given epistemic uncertainty bounds. We apply and test interval and sparse-sample probabilistic approaches to account for uncertainty in the estimated parameters (and various scalar functionals of the parameters as generic quantities of interest, QOIs) due to uncertainties in the material properties. We also investigate the error effects of limited numbers of vibration sensors (accelerometers) on the kettlebell and plate, and therefore abbreviated excitation/response information in the parameter inversions. We propose and demonstrate a Leave-K-Sensors-Out “cross-prediction” UQ approach to estimate related uncertainties on the parameters and QOI functionals. We indicate how uncertainties from material properties and limited sensors are treated in a combined manner. The economical combined UQ approach involves just three to five samples (i.e. three to five inverse simulations), with no added complication or error/uncertainty from use of surrogate models for affordability. Finally, we describe a related economical UQ approach for handling potential parameter solution non-uniqueness and numerical optimization related precision uncertainties in the estimated parameter values. Indicated further research is identified.

36 MATERIALS SCIENCE↗
DOE OSTI2025

The Triple Catalytic Action of Tertiary Nitrogen Catalysts in Recyclable Epoxy-Anhydride Thermosets

The thermosetting polymer matrix in fiber reinforced composites is an important component for energy related applications, such as the lightweighting of vehicles or their use in wind and waterpower turbine blades, due to their ability to provide superior adhesion, stiffness, and applicability to a wide range of manufacturing processes. Despite these benefits, today's thermosets are widely considered to be unrecyclable; thus, there is a large interest in redesigning these materials to be inherently recyclable so that energy intensive production of fibers and monomers can be circumvented, bolstering composite manufacture supply chains. Polyester covalent adaptable networks (PECANs) are one such promising alternative to the incumbent, nonrecyclable epoxy-amine thermosets. PECANs can be formed from the ring-opening co-polymerization (ROCOP) of epoxy-anhydride monomer mixtures and subsequent curing at mild temperatures to exhibit similar performance to conventional epoxies while also possessing unique dynamic chemistries along the ester-hydroxyl backbone that are capable of transesterification and thus reprocessability. While significant advancements have been made in formulating these materials for improved mechanical properties or optimizing solvolysis and reprocessing strategies, less attention has been placed on the impact of the residing amine catalyst used to generate the polyester network. In this work, we evaluated the triple-catalytic efficacy of 12 tertiary amines that act as a curing (bulk ROCOP), a transesterification (internal bond exchange), and a deconstruction (methanolysis) catalyst for PECAN thermosets. Specifically, we first distinguish between chain-growth and step-growth polymerization mechanisms for epoxy-amine and epoxy-anhydride mechanisms. We also utilized density functional theory (DFT) to estimate the basicity (pKb) of each catalyst. Of the tested catalysts, the ROCOP of the studied PECAN network can be completed between 95 and 247 min (at 80 degrees C), with variable gelation phenomena. Additionally, the stress relaxation (transesterification metric) efficiency of the tested PECAN networks with alternative embedded catalysts ranged from 95% to 15% reduction in stress after 5 h at 200 degrees C, and the depolymerization efficacy ranged from 2.5% to 9.8% deconstruction after 36 h at 130 degrees C. Overall, the nitrogen-based moieties were demonstrated to influence polymerization kinetics, catalyze the dynamic transesterification exchange mechanism, and aid in the solvolysis of the thermosets at end-of-life.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Non-uniform active learning for Gaussian process models with applications to trajectory informed aerodynamic databases

The ability to non-uniformly weight the input space is desirable for many applications, and has been explored for space-filling approaches. Increased interests in linking models, such as in a digital twinning framework, increases the need for sampling emulators where they are most likely to be evaluated. In particular, here we apply non-uniform sampling methods for the construction of aerodynamic databases. This paper combines non-uniform weighting with active learning for Gaussian Processes (GPs) to develop a closed-form solution to a non-uniform active learning criterion. We accomplish this by utilizing a kernel density estimator as the weight function. We demonstrate the need and efficacy of this approach with an atmospheric entry example that accounts for both model uncertainty as well as the practical state space of the vehicle, as determined by forward modeling within the active learning loop.

42 ENGINEERING↗