Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model 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 829 records · Page 46

A 60 kpc Galactic Wind Cone in NGC 3079

Galactic winds are associated with intense star formation and active galactic nuclei. Depending on their formation mechanism and velocity, they may remove a significant fraction of gas from their host galaxies, thus suppressing star formation, enriching the intergalactic medium, and shaping the circumgalactic gas. However, the long-term evolution of these winds remains mostly unknown. We report the detection of a wind from NGC 3079 to at least 60 kpc from the galaxy. We detect the wind in far-ultraviolet (FUV) line emission to 60 kpc (as inferred from the broad FUV filter in the Galaxy Evolution Explorer) and X-rays to at least 30 kpc. The morphology, luminosities, temperatures, and densities indicate that the emission comes from shocked material, and the O/Fe ratio implies that the X-ray-emitting gas is enriched by Type II supernovae. If so, the speed inferred from simple shock models is about 500 km s−1, which is sufficient to escape the galaxy. However, the inferred kinetic energy in the wind from visible components is substantially smaller than canonical hot superwind models.

Edmund J. Hodges-Kluck↗

Inference of bipolar neutrino flavor oscillations near a core-collapse supernova based on multiple measurements at Earth

Neutrinos in compact-object environments, such as core-collapse supernovae, can experience various kinds of collective effects in flavor space, engendered by neutrino-neutrino interactions. These include “bipolar” collective oscillations, which are exhibited by neutrino ensembles where different flavors dominate at different energies. Considering the importance of neutrinos in the dynamics and nucleosynthesis in these environments, it is desirable to ascertain whether an Earth-based detection could contain signatures of bipolar oscillations that occurred within a supernova envelope. To that end, we, in this study, continue examining a cost-function formulation of statistical data assimilation (SDA) to infer solutions to a small-scale model of neutrino flavor transformation. SDA is an inference paradigm designed to optimize a model with sparse data. Our model consists of two monoenergetic neutrino beams with different energies emanating from a source and coherently interacting with each other and with a matter background, with radially varying interaction strengths. We attempt to infer flavor transformation histories of these beams using simulated measurements of the flavor content at locations “in vacuum” (that is, far from the source), which could in principle correspond to Earth-based detectors. Within the scope of this small-scale model, we found that: (i) based on such measurements, the SDA procedure is able to infer whether bipolar oscillations had occurred within the protoneutron star envelope, and (ii) if the measurements sample the full amplitude of the neutrino oscillations in vacuum, then the amplitude of the prior bipolar oscillations is well predicted. This result intimates that the inference paradigm can well complement numerical integration codes, via its ability to infer flavor evolution at physically inaccessible locations.

79 ASTRONOMY AND ASTROPHYSICS↗

Implementation of a Binary Neural Network on a Passive Array of Magnetic Tunnel Junctions

The increasing scale of neural networks and their growing application space have produced demand for more energy- and memory-efficient artificial-intelligence-specific hardware. Avenues to mitigate the main issue, the von Neumann bottleneck, include in-memory and near-memory architectures, as well as algorithmic approaches. In this report we leverage the low-power and the inherently binary operation of magnetic tunnel junctions (MTJs) to demonstrate neural network hardware inference based on passive arrays of MTJs. In general, transferring a trained network model to hardware for inference is confronted by degradation in performance due to device-to-device variations, write errors, parasitic resistance, and nonidealities in the substrate. To quantify the effect of these hardware realities, we benchmark 300 unique weight matrix solutions of a two-layer perceptron to classify the Wine dataset for both classification accuracy and write fidelity. Despite device imperfections, we achieve software-equivalent accuracy of up to 95.3% with proper tuning of network parameters in 15 x 15 MTJ arrays having a range of device sizes. The success of this tuning process shows that new metrics are needed to characterize the performance and quality of networks reproduced in mixed signal hardware.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Bayesian differential programming for robust systems identification under uncertainty

This paper presents a machine learning framework for Bayesian systems identification from noisy, sparse and irregular observations of nonlinear dynamical systems. The proposed method takes advantage of recent developments in differentiable programming to propagate gradient information through ordinary differential equation solvers and perform Bayesian inference with respect to unknown model parameters using Hamiltonian Monte Carlo sampling. This allows an efficient inference of the posterior distributions over plausible models with quantified uncertainty, while the use of sparsity-promoting priors enables the discovery of interpretable and parsimonious representations for the underlying latent dynamics. A series of numerical studies is presented to demonstrate the effectiveness of the proposed methods, including nonlinear oscillators, predator–prey systems and examples from systems biology. Taken together, our findings put forth a flexible and robust workflow for data-driven model discovery under uncertainty. All codes and data accompanying this article are available at https://bit.ly/34FOJMj .

Science & Technology - Other Topics↗

INFERRING THE CORONAL DENSITY IRREGULARITY FROM EUV SPECTRA

Understanding the density structure of the solar corona is important for modeling both coronal heating and the solar wind. Direct measurements are difficult because of line-of-sight integration and possible unresolved structures. We present a new method for quantifying such structures using density-sensitive extreme ultraviolet line intensities to derive a density irregularity parameter, a relative measure of the amount of structure along the line of sight. We also present a simple model to relate the inferred irregularities to physical quantities, such as the filling factor and density contrast. For quiet-Sun regions and interplume regions of coronal holes, we find a density contrast of at least a factor of 3–10 and corresponding filling factors of about 10%–20%. Our results are in rough agreement with other estimates of the density structures in these regions. The irregularity diagnostic provides a useful relative measure of unresolved structure in various regions of the corona.

79 ASTRONOMY AND ASTROPHYSICS↗

Computer Model Calibration with Time Series Data Using Deep Learning and Quantile Regression

Computer models play a key role in many scientific and engineering problems. One major source of uncertainty in computer model experiments is input parameter uncertainty. Computer model calibration is a formal statistical procedure to infer input parameters by combining information from model runs and observational data. The existing standard calibration framework suffers from inferential issues when the model output and observational data are high-dimensional dependent data, such as large time series, due to the difficulty in building an emulator and the nonidentifiability between effects from input parameters and data-model discrepancy. To overcome these challenges, we, in this work, propose a new calibration framework based on a deep neural network (DNN) with long short-term memory layers that directly emulates the inverse relationship between the model output and input parameters. Adopting the “learning with noise” idea, we train our DNN model to filter out the effects from data-model discrepancy on input parameter inference. We also formulate a new way to construct interval predictions for DNN using quantile regression to quantify the uncertainty in input parameter estimates. Through a simulation study and real data application with the Weather Research and Forecasting Model Hydrological modeling system (WRF-Hydro), we show our approach can yield accurate point estimates and well-calibrated interval estimates for input parameters.

97 MATHEMATICS AND COMPUTING↗

TRAJECTORIES AND DISTRIBUTION OF INTERSTELLAR DUST GRAINS IN THE HELIOSPHERE

The solar wind carves a bubble in the surrounding interstellar medium (ISM) known as the heliosphere. Charged interstellar dust grains (ISDG) encountering the heliosphere may be diverted around the heliopause or penetrate it depending on their charge-to-mass ratio. We present new calculations of trajectories of ISDG in the heliosphere, and the dust density distributions that result. We include up-to-date grain charging calculations using a realistic UV radiation field and full three-dimensional magnetohydrodynamic fluid + kinetic models for the heliosphere. Models with two different (constant) polarities for the solar wind magnetic field (SWMF) are used, with the grain trajectory calculations done separately for each polarity. Small grains a {sub gr} {approx}< 0.01 {mu}m are completely excluded from the inner heliosphere. Large grains, a {sub gr} {approx}> 1.0 {mu}m, pass into the inner solar system and are concentrated near the Sun by its gravity. Trajectories of intermediate size grains depend strongly on the SWMF polarity. When the field has magnetic north pointing to ecliptic north, the field de-focuses the grains resulting in low densities in the inner heliosphere, while for the opposite polarity the dust is focused near the Sun. The ISDG density outside the heliosphere inferred from applying the model results to in situ dust measurements is inconsistent with local ISM depletion data for both SWMF polarities but is bracketed by them. This result points to the need to include the time variation in the SWMF polarity during grain propagation. Our results provide valuable insights for interpretation of the in situ dust observations from Ulysses.

79 ASTRONOMY AND ASTROPHYSICS↗

Uncertainty quantification and propagation in nuclear density functional theory

Nuclear density functional theory (DFT) is one of the main theoretical tools used to study the properties of heavy and superheavy elements, or to describe the structure of nuclei far from stability. While on-going eff orts seek to better root nuclear DFT in the theory of nuclear forces, energy functionals remain semi-phenomenological constructions that depend on a set of parameters adjusted to experimental data in fi nite nuclei. In this study, we review recent eff orts to quantify the related uncertainties, and propagate them to model predictions. In particular, we cover the topics of parameter estimation for inverse problems, statistical analysis of model uncertainties and Bayesian inference methods. Illustrative examples are taken from the literature.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DeepLensSBI: Deep inference of simulated strong lenses in ground-based surveys

This code is used to train and test machine learning models and generate results and plots presented in 2501.08524 [astro-ph.IM]. The code is written in python. The goal of this work is to train ML models trained on simulated images of strong gravitational lenses. The trained model can then quickly infer properties of the lensed objects with uncertainty quantification.

Poh, Jason [Univ. of Chicago, IL (United States)] ↗

The Effect of the Prior and the Experimental Design on the Inference of the Precision Matrix in Gaussian Chain Graph Models

Here, we investigate whether (and how) experimental design could aid in the estimation of the precision matrix in a Gaussian chain graph model, especially the interplay between the design, the effect of the experiment and prior knowledge about the effect. Estimation of the precision matrix is a fundamental task to infer biological graphical structures like microbial networks. We compare the marginal posterior precision of the precision matrix under four priors: flat, conjugate Normal-Wishart, Normal-MGIG and a general independent. Under the flat and conjugate priors, the Laplace-approximated posterior precision is not a function of the design matrix rendering useless any efforts to find an optimal experimental design to infer the precision matrix. In contrast, the Normal-MGIG and general independent priors do allow for the search of optimal experimental designs, yet there is a sharp upper bound on the information that can be extracted from a given experiment. We confirm our theoretical findings via a simulation study comparing (i) the KL divergence between prior and posterior and (ii) the Stein’s loss difference of MAPs between random and no experiment. Our findings provide practical advice for domain scientists conducting experiments to better infer the precision matrix as a representation of a biological network.

54 ENVIRONMENTAL SCIENCES↗

Modeling Harris Current Sheets with Themis Observations

Current sheets are ubiquitous in nature. occurring in such varied locations as the solar atmosphere. the heliosphere, and the Earth's magnetosphere. The simplest current sheet is the one-dimensional Harris neutral sheet, with the lobe field strength and scale-height the only free parameters. Despite its simplicity, confirmation of the Harris sheet as a reasonable description of the Earth's current sheet has remained elusive. In early 2009 the orbits of the 5 THEMIS probes fortuitously aligned such that profiles of the Earth's current sheet could be modeled in a time dependent manner. For the few hours of alignment we have calculated the time history of the current sheet parameters (scale height and current) in the near-Earth region. during both quiet and active times. For one particular substorm. we further demonstrate good quantitative agreement with the diversion of cross tail current inferred from the Harris modeling with the ionospheric current inferred from ground magnetometer data.

Kepko, L.↗

Variability of moisture in the upper troposphere as inferred from TOVS satellite observations and the ECMWF model analyses in 1989

Large-scale variability of moisture in the upper troposphere is examined using TIROS (Television Infrared Observation Satellite) Operational Vertical Sounder (TOVS) satellite observations and European Centre for Medium-Range Weather Forecasts (ECMWF) model analyses from 1989 in the latitude band from 40 deg N to 40 S. To compare these dissimilar datasets, upwelling radiances were computed for the 6-7 micrometer water vapor band from the ECMWF temperature and moisture analyses, and these computed radiances were compared to the corresponding TOVS satellite observations. The ECMWF-based radiances reproduce the general locations and seasonal cycle of the TOVS-observed moisture features, particularly after an improved convective parameterization scheme was adopted by ECMWF in May 1989. However, the ECMWF analysis scheme still results in much milder lateral moisture gradients and seasonal contrasts than indicated by the TOVS observations. Seasonally, the upper troposphere in each hemisphere dries in winter and moistens in summer, but there are regions in each hemisphere that run counter to this seasonal trend, apparently depending on continental- and monsoon-scale dynamics. Dynamically, the TOVS-observed regions of significant subtropical dryness are correlated with persistent subsidence indicated by ECMWF 300-mb vertical velocity analyses. The TOVS radiance observations indicate large variations in space and time of the upper-tropospheric moisture field, which are not fully captured by the ECMWF analyses.

Salathe, Eric P., Jr.↗