Search NASA⌕ Search

SEARCH · Search NASA

Results for “Parametric Inference”

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

A Conditional Autoencoder for Galaxy Photometric Parameter Estimation

Astronomical photometric surveys routinely image billions of galaxies, and traditionally infer the parameters of a parametric model for each galaxy. This approach has served us well, but the computational expense of deriving a full posterior probability distribution function is a challenge for increasingly ambitious surveys. In this paper, we use deep learning methods to characterize galaxy images, training a conditional autoencoder on mock data. The autoencoder can reconstruct and denoise galaxy images via a latent space engineered to include semantically meaningful parameters, such as brightness, location, size, and shape. Our model recovers galaxy fluxes and shapes on mock data with a lower variance than the Hyper Suprime-Cam photometry pipeline, and returns reasonable answers even for inputs outside the range of its training data. When applied to data in the training range, the regression errors on all extracted parameters are nearly unbiased with a variance near the Cramr-Rao bound.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Obscuration-dependent evolution of active galactic nuclei

We aim to constrain the evolution of active galactic nuclei (AGNs) as a function of obscuration using an X-ray-selected sample of ∼2000 AGNs from a multi-tiered survey including the CDFS, AEGIS-XD, COSMOS, and XMM-XXL fields. The spectra of individual X-ray sources are analyzed using a Bayesian methodology with a physically realistic model to infer the posterior distribution of the hydrogen column density and intrinsic X-ray luminosity. We develop a novel non-parametric method that allows us to robustly infer the distribution of the AGN population in X-ray luminosity, redshift, and obscuring column density, relying only on minimal smoothness assumptions. Our analysis properly incorporates uncertainties from low count spectra, photometric redshift measurements, association incompleteness, and the limited sample size. We find that obscured AGNs with N {sub H} > 10{sup 22} cm{sup –2} account for 77{sub −5}{sup +4}% of the number density and luminosity density of the accretion supermassive black hole population with L {sub X} > 10{sup 43} erg s{sup –1}, averaged over cosmic time. Compton-thick AGNs account for approximately half the number and luminosity density of the obscured population, and 38{sub −7}{sup +8}% of the total. We also find evidence that the evolution is obscuration dependent, with the strongest evolution around N {sub H} ≈ 10{sup 23} cm{sup –2}. We highlight this by measuring the obscured fraction in Compton-thin AGNs, which increases toward z ∼ 3, where it is 25% higher than the local value. In contrast, the fraction of Compton-thick AGNs is consistent with being constant at ≈35%, independent of redshift and accretion luminosity. We discuss our findings in the context of existing models and conclude that the observed evolution is, to first order, a side effect of anti-hierarchical growth.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2008

Certifying Auto-Generated Flight Code

Model-based design and automated code generation are being used increasingly at NASA. Many NASA projects now use MathWorks Simulink and Real-Time Workshop for at least some of their modeling and code development. However, there are substantial obstacles to more widespread adoption of code generators in safety-critical domains. Since code generators are typically not qualified, there is no guarantee that their output is correct, and consequently the generated code still needs to be fully tested and certified. Moreover, the regeneration of code can require complete recertification, which offsets many of the advantages of using a generator. Indeed, manual review of autocode can be more challenging than for hand-written code. Since the direct V&V of code generators is too laborious and complicated due to their complex (and often proprietary) nature, we have developed a generator plug-in to support the certification of the auto-generated code. Specifically, the AutoCert tool supports certification by formally verifying that the generated code is free of different safety violations, by constructing an independently verifiable certificate, and by explaining its analysis in a textual form suitable for code reviews. The generated documentation also contains substantial tracing information, allowing users to trace between model, code, documentation, and V&V artifacts. This enables missions to obtain assurance about the safety and reliability of the code without excessive manual V&V effort and, as a consequence, eases the acceptance of code generators in safety-critical contexts. The generation of explicit certificates and textual reports is particularly well-suited to supporting independent V&V. The primary contribution of this approach is the combination of human-friendly documentation with formal analysis. The key technical idea is to exploit the idiomatic nature of auto-generated code in order to automatically infer logical annotations. The annotation inference algorithm itself is generic, and parametrized with respect to a library of coding patterns that depend on the safety policies and the code generator. The patterns characterize the notions of definitions and uses that are specific to the given safety property. For example, for initialization safety, definitions correspond to variable initializations while uses are statements which read a variable, whereas for array bounds safety, definitions are the array declarations, while uses are statements which access an array variable. The inferred annotations are thus highly dependent on the actual program and the properties being proven. The annotations, themselves, need not be trusted, but are crucial to obtain the automatic formal verification of the safety properties without requiring access to the internals of the code generator. The approach has been applied to both in-house and commercial code generators, but is independent of the particular generator used. It is currently being adapted to flight code generated using MathWorks Real-Time Workshop, an automatic code generator that translates from Simulink/Stateflow models into embedded C code.

Denney, Ewen↗
DOE OSTI2011

Global Neutrino Data Analysis and the Quest to Pin Down $\sin\theta_{13}$ in Different Mixing Matrix Parametrizations

Hints for sizable $\sin^2 \theta_{13}$ have been reported in earlier global neutrino oscillation data analyses as well as will be reported in this work, and quite recently by the Double Chooz experiment. However, as we enter the era of precision neutrino oscillation experiments, terms linear in $\sin\theta_{13}$ will no longer be negligible, and its sign would affect the extraction of other oscillation parameters. The sign of $\sin\theta_{13}$ also plays a crucial role in the determination of the CP-violating phase. In this work we show that by adopting an alternative parametrization for the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) mixing matrix, one already has a chance to infer the sign of each mixing angle in the conventional parametrization using existing global neutrino data. A weak preference for negative $\sin \theta_{13}$ is found. In particular, the solar data suggest that $\sin\theta_{13} > 0$ while all other data the opposite. This leads to the speculation on whether the Mikheyev-Smirnov-Wolfenstein (MSW) effect is responsible. In this work we found that in the new mixing matrix parametrization, the 68% CL constraints on the three mixing angles are comparable to those estimated in the conventional parametrization adopted in the literature. Owing to the strong correlations among the three mixing angles in the new parametrization, the advantages of doing the global neutrino oscillation analysis using data from past, current, and near future neutrino oscillation experiments shall become manifest.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Bayesian uncertainty quantification of tristructural isotropic particle fuel silver release: Decomposing model inadequacy plus experimental noise and parametric uncertainties

Tristructural isotropic (TRISO) particle fuel is one of the most promising fuel concepts enabling high temperature and high burnup reactor operation. One dominant source of radioactivity released from the TRISO particles is silver (Ag), which is subject to a high release fraction and long decay life compared to other fission products. Previous modeling efforts using the fuel performance code BISON indicated nonnegligible uncertainties in modeling the diffusion process of fission products in TRISO compared to the Advanced Gas Reactor experiments. The overall uncertainties observed when modeling the fission product diffusion can result from uncertainties in model parameters, noisy experimental measurements, and deficiencies in the developed models. The three types of underlying uncertainties have not yet been properly quantified in open literature. Here, this paper presents the Bayesian uncertainty quantification (UQ) using massively parallelizable Markov chain Monte Carlo samplers. The uncertainties due to model parameters, model inadequacy, and experimental measurement noise are quantified, with the σ term used to represent the sum of the model inadequacy and measurement noise uncertainties. It is worth noting that this is the first time the σ term is inferred for nuclear fuel experiments, as compared to using prescribed values for uncertainty quantification in previous work. The parallelizable Markov chain Monte Carlo samplers efficiently infer the model parameters and the σ term, giving insight into physical parameters like diffusion coefficients and the combined model discrepancy and measurement noise. A subsequent forward uncertainty quantification (UQ) is also performed based on the calibration results to generate more accurate predictions of the Ag release. The model inadequacy plus experimental noise is the most dominant source of uncertainty compared to the parametric uncertainty. All the UQ analyses presented in this work are based on the second series of the irradiation experiments in the Advanced Gas Reactor program.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS2012

Parametric 3D Atmospheric Reconstruction in Highly Variable Terrain with Recycled Monte Carlo Paths and an Adapted Bayesian Inference Engine

We describe a method for accelerating a 3D Monte Carlo forward radiative transfer model to the point where it can be used in a new kind of Bayesian retrieval framework. The remote sensing challenge is to detect and quantify a chemical effluent of a known absorbing gas produced by an industrial facility in a deep valley. The available data is a single low resolution noisy image of the scene in the near IR at an absorbing wavelength for the gas of interest. The detected sunlight has been multiply reflected by the variable terrain and/or scattered by an aerosol that is assumed partially known and partially unknown. We thus introduce a new class of remote sensing algorithms best described as "multi-pixel" techniques that call necessarily for a 3D radaitive transfer model (but demonstrated here in 2D); they can be added to conventional ones that exploit typically multi- or hyper-spectral data, sometimes with multi-angle capability, with or without information about polarization. The novel Bayesian inference methodology uses adaptively, with efficiency in mind, the fact that a Monte Carlo forward model has a known and controllable uncertainty depending on the number of sun-to-detector paths used.

Monte Carlo↗
DOE OSTI2014

Observational implications of mattergenesis during inflation

The observed baryon asymmetry, as well as potentially an asymmetry in the dark matter sector, can be produced through dissipative particle production during inflation. A distinctive feature of this mechanism is the generation of matter isocurvature perturbations that are fully (anti-)correlated with the dominant adiabatic curvature perturbations. We show that chaotic warm inflation models yield anti-correlated isocurvature modes that may partially or even completely screen the contribution of primordial gravity waves to the CMB temperature power spectrum. The tensor-to-scalar ratio inferred from the latter may thus be parametrically smaller than the one deduced from B-mode polarization maps, which is particularly relevant in the light of the recently announced results of the BICEP2 experiment.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1978

Nitric oxide formation in gas turbine engines: A theoretical and experimental study

A modified Zeldovich kinetic scheme was used to predict nitric oxide formation in the burned gases. Nonuniformities in fuel-air ratio in the primary zone were accounted for by a distribution of fuel-air ratios. This was followed by one or more dilution zones in which a Monte Carlo calculation was employed to follow the mixing and dilution processes. Predictions of NOX emissions were compared with various available experimental data, and satisfactory agreement was achieved. In particular, the model is applied to the NASA swirl-can modular combustor. The operating characteristics of this combustor which can be inferred from the modeling predictions are described. Parametric studies are presented which examine the influence of the modeling parameters on the NOX emission level. A series of flow visualization experiments demonstrates the fuel droplet breakup and turbulent recirculation processes. A tracer experiment quantitatively follows the jets from the swirler as they move downstream and entrain surrounding gases. Techniques were developed for calculating both fuel-air ratio and degree of nonuniformity from measurements of CO2, CO, O2, and hydrocarbons. A burning experiment made use of these techniques to map out the flow field in terms of local equivalence ratio and mixture nonuniformity.

Mikus, T.↗
DOE OSTI2020

Bayesian Inference of the Symmetry Energy of Superdense Neutron-rich Matter from Future Radius Measurements of Massive Neutron Stars

In this work, using as references the posterior probability distribution functions of the equation of state (EOS) parameters inferred from the radii of canonical neutron stars (NSs) reported by the LIGO/VIRGO and NICER Collaborations based on their observations of GW170817 and PSR J0030+0451, we investigate how future radius measurements of more massive NSs will improve our current knowledge about the EOS of superdense neutron-rich nuclear matter, especially its symmetry energy term. Within the Bayesian statistical approach using an explicitly isospin-dependent parametric EOS for the core of NSs, we infer the EOS parameters of superdense neutron-rich nuclear matter from three sets of imagined mass–radius correlation data representing typical predictions by various nuclear many-body theories, that is, the radius stays the same, decreases, or increases with increasing NS mass within ±15% between 1.4 and 2.0 M$_\odot$. The corresponding NS average density increases quickly or slowly or slightly decreases as the NS mass increases from 1.4 to 2.0 M$_\odot$. While the EOSs of symmetric nuclear matter (SNM) inferred from the three data sets are approximately the same, the corresponding symmetry energies above about twice the saturation density of nuclear matter are very different, indicating that the radii of massive NSs carry important information about the high-density behavior of nuclear symmetry energy with little influence from the remaining uncertainties of the SNM EOS at suprasaturation densities.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

A multishock model for the density variance of anisotropic, highly magnetized, supersonic turbulence

ABSTRACT Shocks form the basis of our understanding for the density and velocity statistics of supersonic turbulent flows, such as those found in the cool interstellar medium (ISM). The variance of the density field, $\sigma ^2_{\rho /\rho _0}$, is of particular interest for molecular clouds (MCs), the birthplaces of stars in the Universe. The density variance may be used to infer underlying physical processes in an MC, and parametrizes the star formation (SF) rate of a cloud. However, models for $\sigma ^2_{\rho /\rho _0}$ all share a common feature – the variance is assumed to be isotropic. This assumption does not hold when a trans-/sub-Alfvénic mean magnetic field, ${B}_0$, is present in the cloud, which observations suggest is relevant for some MCs. We develop an anisotropic model for $\sigma _{\rho /\rho _0}^2$, using contributions from hydrodynamical and fast magnetosonic shocks that propagate orthogonal to each other. Our model predicts an upper bound for $\sigma _{\rho /\rho _0}^2$ in the high Mach number $(\mathcal {M})$ limit as small-scale density fluctuations become suppressed by the strong ${B}_0$. The model reduces to the isotropic $\sigma _{\rho /\rho _0}^2\!-\!\mathcal {M}$ relation in the hydrodynamical limit. To validate our model, we calculate $\sigma _{\rho /\rho _0}^2$ from 12 high-resolution, three-dimensional, supersonic, sub-Alfvénic magnetohydrodynamical (MHD) turbulence simulations and find good agreement with our theory. We discuss how the two MHD shocks may be the bimodally oriented overdensities observed in some MCs and the implications for SF theory in the presence of a sub-Alfvénic ${B}_0$. By creating an anisotropic, supersonic density fluctuation model, this study paves the way for SF theory in the highly anisotropic regime of interstellar turbulence.

Beattie, James R.↗
DOE OSTI2021

On the Stochastic Stability of Deep Markov Models

Deep Markov models (DMM) are generative models which are scalable and expressive generalization of Markov models for representation, learning, and inference problems. DMMs using deep neural networks to parametrize the transition of Markov probability distributions have recently been shown to provide more expressiveness in modeling sequential data and dynamical system responses. However, the fundamental stochastic stability guarantees of such models have not been thoroughly investigated. In this paper, we present a rigorous analytical method to prove the necessary and sufficient conditions of DMM's stochastic stability. This task is achieved by spectral analysis of the efficiently computed Jacobians of probabilistic maps modeled by deep neural networks. We make theoretical connections between the eigenvalues of neural network's weights and the different activation function types used on the stability and overall dynamic behavior of DMMs with Gaussian distributions. We empirically substantiate our theoretical results on stochastic stability and eigenvalue spectra via several numerical experiments. Formal stability guarantees of DMMs can substantially improve their robustness and trustworthiness, necessary for reliable use in safety-critical real-world applications.

Drgona, Jan↗
NASA NTRS1999

Development of Corrections for Biomass Burning Effects in Version 2 of GEWEX/SRB Algorithm

The objectives of this project were: (1) To incorporate into an existing version of the University of Maryland Surface Radiation Budget (SRB) model, optical parameters of forest fire aerosols, using best available information, as well as optical properties of other aerosols, identified as significant. (2) To run the model on regional scales with the new parametrization and information on forest fire occurrence and plume advection, as available from NASA LARC, and test improvements in inferring surface fluxes against daily values of measured fluxes. (3) Develop strategy how to incorporate the new parametrization on global scale and how to transfer modified model to NASA LARC.

Pinker, Rachel T.↗
NASA NTRS1982

Quantiles, parametric-select density estimation, and bi-information parameter estimators

A quantile-based approach to statistical analysis and probability modeling of data is presented which formulates statistical inference problems as functional inference problems in which the parameters to be estimated are density functions. Density estimators can be non-parametric (computed independently of model identified) or parametric-select (approximated by finite parametric models that can provide standard models whose fit can be tested). Exponential models and autoregressive models are approximating densities which can be justified as maximum entropy for respectively the entropy of a probability density and the entropy of a quantile density. Applications of these ideas are outlined to the problems of modeling: (1) univariate data; (2) bivariate data and tests for independence; and (3) two samples and likelihood ratios. It is proposed that bi-information estimation of a density function can be developed by analogy to the problem of identification of regression models.

Parzen, E.↗
DOE OSTI2025

Discriminative versus generative approaches to simulation-based inference

Most of the fundamental, emergent, and phenomenological parameters of particle and nuclear physics are determined through parametric template fits. Simulations are used to populate histograms which are then matched to data. This approach is inherently lossy, since histograms are binned and low-dimensional. Deep learning has enabled unbinned and high-dimensional parameter estimation through neural likelihood(-ratio) estimation. We compare two approaches for neural simulation-based inference (NSBI): one based on discriminative learning (classification) and one based on generative modeling. These two approaches are directly evaluated on the same datasets, with a similar level of hyperparameter optimization in both cases. In addition to a Gaussian dataset, we study NSBI using a Higgs boson dataset from the FAIR Universe Challenge. We find that both the direct likelihood and likelihood ratio estimation are able to effectively extract parameters with reasonable uncertainties. For the numerical examples and within the set of hyperparameters studied, we found that the likelihood ratio method is more accurate and/or precise. Both methods have a significant spread from the network training and would require ensembling or other mitigation strategies in practice.

high energy physics↗
DOE OSTI2020

Turbulence model reduction by deep learning

A defining problem of turbulence theory is to produce a predictive model for turbulent fluxes. These have profound implications for virtually all aspects of the turbulence dynamics. In magnetic confinement devices, drift-wave turbulence produces anomalous fluxes via cross-correlations between fluctuations. In this work, we introduce an alternative, data-driven method for parametrizing these fluxes. The method uses deep supervised learning to infer a reduced mean-field model from a set of numerical simulations. We apply the method to a simple drift-wave turbulence system and find a significant new effect which couples the particle flux to the local gradient of vorticity. Notably, here, this effect is much stronger than the oft-invoked shear suppression effect. We also recover the result via a simple calculation. The vorticity gradient effect tends to modulate the density profile. In addition, our method recovers a model for spontaneous zonal flow generation by negative viscosity, stabilized by nonlinear and hyperviscous terms. We highlight the important role of symmetry to implementation of the new method.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Application of Sequential Design of Experiments (SDoE) to Large Pilot-Scale Solvent-Based CO2 Capture Process at Technology Centre Mongstad (TCM)

The United States Department of Energy’s Carbon Capture Simulation for Industry Impact (CCSI2) program has developed a framework for sequential design of experiments (SDoE) that aims to maximize knowledge gained from budget- and schedule-limited pilot scale testing. SDoE was applied to the planning and execution of campaigns for testing CO2 capture systems at pilot-scale in order to optimally allocate resources available for the testing. In this methodology, a stochastic process model is developed by quantifying the parametric uncertainty in submodels of interest; for a solvent-based CO2 capture system, these may include physical properties and equipment performance submodels (e.g., mass transfer, interfacial area). This uncertainty is propagated through the full process model, over variable operating conditions, for estimating the resulting uncertainty in key model outputs (e.g., percentage of CO2 capture, solvent regeneration energy requirement). In developing a data collection plan, the predicted output uncertainty is incorporated into an algorithm that seeks simultaneously to select process operating conditions for which the predicted uncertainty is relatively high and to ensure that the entire space of operation is well represented. This test plan is then used to guide operation of the pilot plant at varying steady-state conditions, with resulting process data incorporated into the existing model using Bayesian inference to refine parameter distributions. The updated stochastic model, with reduced parametric uncertainty from data collected, is then used to guide additional data collection, thus the sequential nature of the experimental design. The SDoE process was implemented at the pilot test unit (12 MWe in scale) at Norway’s Technology Centre Mongstad (TCM) in a summer 2018 test campaign with aqueous monoethanolamine (MEA). During the test campaign, the varied operating conditions included the flowrates of circulated solvent, flue gas, and reboiler steam and the CO2 concentration in the flue gas. The process data were used to update probability distributions of mass transfer and interfacial area parameters of a stochastic process model developed by the CCSI2 team. Two iterations of the SDoE process were executed, resulting in the uncertainty in model predicted CO2 capture percentage decreasing by an average of 58.0 ± 4.7% over the full input space of interest. This work demonstrates the potential of the SDoE process for model refinement through reduction in process model parametric uncertainty, and ultimately risk in scale-up, in CO2 capture technology performance.

carbon capture↗
DOE OSTI2023

Imaging electron angular distributions to assess a full-power petawatt-class laser focus

Here we present a technique to assess the focal volume of petawatt-class lasers at full power. Our approach exploits quantitative measurement of the angular distribution of electrons born in the focus via ionization of rarefied gas, which are accelerated forward and ejected ponderomotively by the field. We show that a bivariate (θ,Φ) angular distribution, which was obtained with image plates, not only enables the peak intensity to be extracted, but also reflects nonideality of the focal-spot intensity distribution. In our prototype demonstration at intensities of a few ×10 19 to a few ×10 20 W/cm 2 , an f/10 optic produced a focal spot in the paraxial regime. This allows a plane-wave parametrization of the peak intensity given by tanθ c =2/a 0 (a 0 being the normalized vector potential and θ c the minimum ejection angle) to be compared with our measurements. Qualitative agreement was found using an a 0 inferred from the pulse energy, pulse duration, and focal spot distribution with a modified parametrization, tanθ c =2η/a 0 ($η=2.02^{+0.26}_{–0.22}$). This highlights the need for (i) better understanding of intensity degradation due to focal-spot distortions and (ii) more robust modeling of the ejection dynamics. Using single-shot detection of electrons, we showed that while there is significant shot-to-shot variation in the number of electrons ejected at a given angular position, the average distribution scales with the pulse energy in a way that is consistent with that seen with the image plates. Finally, we note that the asymptotic behavior as θ→0° limits the usability of angular measurement. For 800 nm, this limit is at an intensity ~10 21 W/cm 2 .

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗