Search NASA⌕ Search

SEARCH · Search NASA

Results for “inference methods”

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 1,765 records · Page 98

DeepM&Mnet: Inferring the electroconvection multiphysics fields based on operator approximation by neural networks

Electroconvection is a multiphysics problem involving coupling of the flow field with the electric field as well as the cation and anion concentration fields. Here, we use electroconvection as a benchmark problem to put forward a new data assimilation framework, the DeepM&Mnet, for simulating multiphysics and multiscale problems at speeds much faster than standard numerical methods using pre-trained neural networks. We first pre-train DeepONets that can predict independently each field, given general inputs from the rest of the fields of the coupled system. DeepONets can approximate nonlinear operators and are composed of two sub-networks, a branch net for the input fields and a trunk net for the locations of the output field. DeepONets, which are extremely fast, are used as building blocks in the DeepM&Mnet and form constraints for the multiphysics solution along with some sparse available measurements of any of the fields. We demonstrate the new methodology and document the accuracy of each individual DeepONet, and subsequently we present two different DeepM&Mnet architectures that infer accurately and efficiently 2D electroconvection fields for unseen electric potentials. Furthermore, the DeepM&Mnet framework is general and can be applied for building any complex multiphysics and multiscale models based on very few measurements using pre-trained DeepONets in a “plug-and-play” mode.

97 MATHEMATICS AND COMPUTING↗

A Bayesian Modeling Approach for Estimation of a Shape-Free Groundwater Age Distribution using Multiple Tracers

The mixing of groundwaters with different ages in aquifers, groundwater age is more appropriately represented by a distribution rather than a scalar number. To infer a groundwater age distribution from environmental tracers, a mathematical form is often assumed for the shape of the distribution and the parameters of the mathematical distribution are estimated using deterministic or stochastic inverse methods. We found that the prescription of the mathematical form limits the exploration of the age distribution to the shapes that can be described by the selected distribution. In this paper, the use of freeform histograms as groundwater age distributions is evaluated. A Bayesian Markov Chain Monte Carlo approach is used to estimate the fraction of groundwater in each histogram bin. This method was able to capture the shape of a hypothetical gamma distribution from the concentrations of four age tracers. The number of bins that can be considered in this approach is limited based on the number of tracers available. The histogram method was also tested on tracer data sets from Holten (The Netherlands; 3 H, 3 He, 85 Kr, 39 Ar) and the La Selva Biological Station (Costa-Rica; SF 6 , CFCs, 3 H, 4 He and 14 C), and compared to a number of mathematical forms. According to standard Bayesian measures of model goodness, the best mathematical distribution performs better than the histogram distributions in terms of the ability to capture the observed tracer data relative to their complexity. Among the histogram distributions, the four bin histogram performs better in most of the cases. The Monte Carlo simulations showed strong correlations in the posterior estimates of bin contributions, indicating that these bins cannot be well constrained using the available age tracers. The fact that mathematical forms overall perform better than the freeform histogram does not undermine the benefit of the freeform approach, especially for the cases where a larger amount of observed data is available and when the real groundwater distribution is more complex than can be represented by simple mathematical forms.

54 ENVIRONMENTAL SCIENCES↗

A worldwide aerosol phenomenology: Elemental and organic carbon in PM 2.5 and PM 10

Elemental carbon (EC), organic carbon (OC), and particulate matter (PM) concentrations in the inhalable (PM 10 ) and fine (PM 2.5 ) size fractions are measured worldwide, albeit with different analytical methods. These measurements from many researchers were collected and analyzed for Africa, America, Asia, and Europe for 2012–2019. EC/PM, OC/PM, and OC/EC ratios were examined based on region, site type, and season to infer potential sources and impacts. These analyses demonstrate that carbonaceous materials are important PM constituents throughout the world. Mean EC/PM ratios were lowest in PM 10 in Sahelian Africa and Europe (∼0.01), highest (>0.07) in PM 2.5 at urban sites in North America, South America, and Japan. Mean OC/PM ratios were lowest in PM 10 in the Sahel (∼0.06) and in PM 2.5 in China and Thailand (0.10), and highest in central and eastern Europe (∼0.3) and North America (∼0.4). OC/EC ratios were elevated in western and northern Europe, and at regional background sites in North America. EC/PM increased with PM 10 in Thailand, while OC/PM increased with higher PM mass in Thailand, India, and North America, highlighting the specific contribution of carbonaceous aerosols to PM pollution in these regions. At European and North American background sites, OC/EC ratios increased with PM mass. Higher OC/EC ratios in dry periods indicate influence of wildfires, prescribed burns, and secondary aerosol formation. Elevated wintertime EC/PM ratios coincide with residential heating in temperate climate zones.

54 ENVIRONMENTAL SCIENCES↗

Analysis of the neutron time-of-flight spectra from inertial confinement fusion experiments

Neutron time-of-flight diagnostics have long been used to characterize the neutron spectrum produced by inertial confinement fusion experiments. The primary diagnostic goals are to extract the d + t → n + α (DT) and d + d → n + {sup 3}He (DD) neutron yields and peak widths, and the amount DT scattering relative to its unscattered yield, also known as the down-scatter ratio (DSR). These quantities are used to infer yield weighted plasma conditions, such as ion temperature (T{sub ion}) and cold fuel areal density. We report on novel methodologies used to determine neutron yield, apparent T{sub ion}, and DSR. These methods invoke a single temperature, static fluid model to describe the neutron peaks from DD and DT reactions and a spline description of the DT spectrum to determine the DSR. Both measurements are performed using a forward modeling technique that includes corrections for line-of-sight attenuation and impulse response of the detection system. These methods produce typical uncertainties for DT T{sub ion} of 250 eV, 7% for DSR, and 9% for the DT neutron yield. For the DD values, the uncertainties are 290 eV for T{sub ion} and 10% for the neutron yield.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Bayesian tomography and integrated data analysis in fusion diagnostics

In this article, a Bayesian tomography method using non-stationary Gaussian process for a prior has been introduced. The Bayesian formalism allows quantities which bear uncertainty to be expressed in the probabilistic form so that the uncertainty of a final solution can be fully resolved from the confidence interval of a posterior probability. Moreover, a consistency check of that solution can be performed by checking whether the misfits between predicted and measured data are reasonably within an assumed data error. In particular, the accuracy of reconstructions is significantly improved by using the non-stationary Gaussian process that can adapt to the varying smoothness of emission distribution. The implementation of this method to a soft X-ray diagnostics on HL-2A has been used to explore relevant physics in equilibrium and MHD instability modes. This project is carried out within a large size inference framework, aiming at an integrated analysis of heterogeneous diagnostics.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Fast convolutional neural networks on FPGAs with hls4ml

We introduce an automated tool for deploying ultra low-latency, low-power deep neural networks with convolutional layers on field-programmable gate arrays (FPGAs). By extending the hls4ml library, we demonstrate an inference latency of 5 µs using convolutional architectures, targeting microsecond latency applications like those at the CERN Large Hadron Collider. Considering benchmark models trained on the Street View House Numbers Dataset, we demonstrate various methods for model compression in order to fit the computational constraints of a typical FPGA device used in trigger and data acquisition systems of particle detectors. In particular, we discuss pruning and quantization-aware training, and demonstrate how resource utilization can be significantly reduced with little to no loss in model accuracy. We show that the FPGA critical resource consumption can be reduced by 97% with zero loss in model accuracy, and by 99% when tolerating a 6% accuracy degradation.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Pump with Broadband Probe Experiments for Single-Shot Measurements of Plasma Conditions and Crossed-Beam Energy Transfer

Here, a novel technique for measuring plasma conditions using monochromatic pump-broadband-probe laser interactions has been experimentally demonstrated. Originally proposed in Ludwig et al. [Phys. Plasmas 26, 113108 (2019)], this method utilizes crossed-beam energy transfer between the broadband probe and the pump, mediated by plasma ion acoustic waves. The complete energy transfer spectrum can be captured in a single shot, enabling the inference of plasma parameters such as density, electron and ion temperatures, and flow velocity. Compared to Thomson scattering, this technique offers signal enhancements typically larger than 9 orders of magnitude, significantly reducing the required probe laser intensity and facilitating interactions that are linear and measurements that are nonperturbative of the plasma. Furthermore, it provides a powerful tool for advancing studies of crossed-beam energy transfer under conditions relevant to inertial confinement fusion experiments.

Physics - Plasma physics↗

Neural network methods for radiation detectors and imaging

Recent advances in image data proccesing through deep learning allow for new optimization and performance-enhancement schemes for radiation detectors and imaging hardware. This enables radiation experiments, which includes photon sciences in synchrotron and X-ray free electron lasers as a subclass, through data-endowed artificial intelligence. We give an overview of data generation at photon sources, deep learning-based methods for image processing tasks, and hardware solutions for deep learning acceleration. Most existing deep learning approaches are trained offline, typically using large amounts of computational resources. However, once trained, DNNs can achieve fast inference speeds and can be deployed to edge devices. A new trend is edge computing with less energy consumption (hundreds of watts or less) and real-time analysis potential. While popularly used for edge computing, electronic-based hardware accelerators ranging from general purpose processors such as central processing units (CPUs) to application-specific integrated circuits (ASICs) are constantly reaching performance limits in latency, energy consumption, and other physical constraints. These limits give rise to next-generation analog neuromorhpic hardware platforms, such as optical neural networks (ONNs), for high parallel, low latency, and low energy computing to boost deep learning acceleration (LA-UR-23-32395).

edge computing↗

Determination of Coronal Mass Ejection physical parameters from a combination of polarized visible light and UV Lyα observations

Visible-light observations of Coronal Mass Ejections (CMEs) performed with coronagraphs and heliospheric imagers (in primis on board the Solar and Heliospheric Observatory and STEREO missions) have offered the best way to study the kinematics and geometrical structure of these fundamental events so far. Nevertheless, it has been widely demonstrated that only combination of multi-wavelength data (including X-ray spectra, EUV images, EUV-UV spectra, and radio dynamic spectra) can provide complete information on the plasma temperature and density distributions, non-thermal motions, magnetic fields, and other physical parameters, for both CMEs and CME-related phenomena. In this work, we analyze three CMEs by combining simultaneous data acquired in the polarized visible light by the LASCO-C2 coronagraph and in the UV H i Lyα line (1216 Å) by the UVCS spectrometer, in order to estimate the CME plasma electron density (using the polarization-ratio technique to infer the 3D structure of the CME) and temperature (from the comparison between the expected and measured Lyα intensities) along the UVCS field of view. This analysis is primarily aimed at testing the diagnostic methods that will be applied to coronagraphic observations of CMEs delivered by the Metis instrument on board the next ESA-Solar Orbiter mission. We find that CME cores are usually associated with cooler plasma (T∼10{sup 6} K), and that a significant increase of the electron temperatures is observed from the core to the front of the CME (where T>10{sup 6.3} K), which seems to be correlated, in all cases, with the morphological structure of the CME as derived from visible-light images.

79 ASTRONOMY AND ASTROPHYSICS↗

Experimental investigation of the excess charge and time constant of minority carriers in the thin diffused layer of 0.1 ohm-cm silicon solar cells

An experimental method is presented that can be used to interpret the relative roles of bandgap narrowing and recombination processes in the diffused layer. This method involves measuring the device time constant by open-circuit voltage decay and the base region diffusion length by X-ray excitation. A unique illuminated diode method is used to obtain the diode saturation current. These data are interpreted using a simple model to determine individually the minority carrier lifetime and the excess charge. These parameters are then used to infer the relative importance of bandgap narrowing and recombination processes in the diffused layer.

Godlewski, M. P.↗

Remote sensing of snow and ice

This paper reviews remote sensing of snow and ice, techniques for improved monitoring, and incorporation of the new data into forecasting and management systems. The snowcover interpretation of visible and infrared data from satellites, automated digital methods, radiative transfer modeling to calculate the solar reflectance of snow, and models using snowcover input data and elevation zones for calculating snowmelt are discussed. The use of visible and near infrared techniques for inferring snow properties, microwave monitoring of snowpack characteristics, use of Landsat images for collecting glacier data, monitoring of river ice with visible imagery from NOAA satellites, use of sequential imagery for tracking ice flow movement, and microwave studies of sea ice are described. Applications of snow and ice research to commercial use are examined, and it is concluded that a major problem to be solved is characterization of snow and ice in nature, since assigning of the correct properties to a real system to be modeled has been difficult.

Rango, A.↗

Building distributed rule-based systems using the AI Bus

The AI Bus software architecture was designed to support the construction of large-scale, production-quality applications in areas of high technology flux, running heterogeneous distributed environments, utilizing a mix of knowledge-based and conventional components. These goals led to its current development as a layered, object-oriented library for cooperative systems. This paper describes the concepts and design of the AI Bus and its implementation status as a library of reusable and customizable objects, structured by layers from operating system interfaces up to high-level knowledge-based agents. Each agent is a semi-autonomous process with specialized expertise, and consists of a number of knowledge sources (a knowledge base and inference engine). Inter-agent communication mechanisms are based on blackboards and Actors-style acquaintances. As a conservative first implementation, we used C++ on top of Unix, and wrapped an embedded Clips with methods for the knowledge source class. This involved designing standard protocols for communication and functions which use these protocols in rules. Embedding several CLIPS objects within a single process was an unexpected problem because of global variables, whose solution involved constructing and recompiling a C++ version of CLIPS. We are currently working on a more radical approach to incorporating CLIPS, by separating out its pattern matcher, rule and fact representations and other components as true object oriented modules.

Schultz, Roger D.↗

Bayesian Approach to the Joint Inversion of Gravity and Magnetic Data, with Application to the Ismenius Area of Mars

This viewgraph presentation reviews a Bayesian approach to the inversion of gravity and magnetic data with specific application to the Ismenius Area of Mars. Many inverse problems encountered in geophysics and planetary science are well known to be non-unique (i.e. inversion of gravity the density structure of a body). In hopes of reducing the non-uniqueness of solutions, there has been interest in the joint analysis of data. An example is the joint inversion of gravity and magnetic data, with the assumption that the same physical anomalies generate both the observed magnetic and gravitational anomalies. In this talk, we formulate the joint analysis of different types of data in a Bayesian framework and apply the formalism to the inference of the density and remanent magnetization structure for a local region in the Ismenius area of Mars. The Bayesian approach allows prior information or constraints in the solutions to be incorporated in the inversion, with the "best" solutions those whose forward predictions most closely match the data while remaining consistent with assumed constraints. The application of this framework to the inversion of gravity and magnetic data on Mars reveals two typical challenges - the forward predictions of the data have a linear dependence on some of the quantities of interest, and non-linear dependence on others (termed the "linear" and "non-linear" variables, respectively). For observations with Gaussian noise, a Bayesian approach to inversion for "linear" variables reduces to a linear filtering problem, with an explicitly computable "error" matrix. However, for models whose forward predictions have non-linear dependencies, inference is no longer given by such a simple linear problem, and moreover, the uncertainty in the solution is no longer completely specified by a computable "error matrix". It is therefore important to develop methods for sampling from the full Bayesian posterior to provide a complete and statistically consistent picture of model uncertainty, and what has been learned from observations. We will discuss advanced numerical techniques, including Monte Carlo Markov

data analysis↗

Evolution of Geometric Sensitivity Derivatives from Computer Aided Design Models

The generation of design parameter sensitivity derivatives is required for gradient-based optimization. Such sensitivity derivatives are elusive at best when working with geometry defined within the solid modeling context of Computer-Aided Design (CAD) systems. Solid modeling CAD systems are often proprietary and always complex, thereby necessitating ad hoc procedures to infer parameter sensitivity. A new perspective is presented that makes direct use of the hierarchical associativity of CAD features to trace their evolution and thereby track design parameter sensitivity. In contrast to ad hoc methods, this method provides a more concise procedure following the model design intent and determining the sensitivity of CAD geometry directly to its respective defining parameters.

Jones, William T.↗

Constraining physical models at gigabar pressures

High-energy-density (HED) experiments in convergent geometry are able to test physical models at pressures beyond hundreds of millions of atmospheres. The measurements from these experiments are generally highly integrated and require unique analysis techniques to procure quantitative infor mation. This work describes a methodology to constrain the physics in convergent HED experiments by adapting the methods common to many other fields of physics. As an example, a mechanical model of an imploding shell is constrained by data from a thin-shelled direct-drive exploding-pusher experiment on the OMEGA Laser System using Bayesian inference, resulting in the reconstruction of the shell dynamics and energy transfer during the implosion. The model is tested by analyzing synthetic data from a 1-D hydrodynamics code and is sampled using a Markov chain Monte Carlo to generate the posterior distributions of the model parameters. As a result, the goal of this work is to demonstrate a general methodology that can be used to draw conclusions from a wide variety of HED experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Physics-Informed Neural Networks for Heat Transfer Problems

Abstract Physics-informed neural networks (PINNs) have gained popularity across different engineering fields due to their effectiveness in solving realistic problems with noisy data and often partially missing physics. In PINNs, automatic differentiation is leveraged to evaluate differential operators without discretization errors, and a multitask learning problem is defined in order to simultaneously fit observed data while respecting the underlying governing laws of physics. Here, we present applications of PINNs to various prototype heat transfer problems, targeting in particular realistic conditions not readily tackled with traditional computational methods. To this end, we first consider forced and mixed convection with unknown thermal boundary conditions on the heated surfaces and aim to obtain the temperature and velocity fields everywhere in the domain, including the boundaries, given some sparse temperature measurements. We also consider the prototype Stefan problem for two-phase flow, aiming to infer the moving interface, the velocity and temperature fields everywhere as well as the different conductivities of a solid and a liquid phase, given a few temperature measurements inside the domain. Finally, we present some realistic industrial applications related to power electronics to highlight the practicality of PINNs as well as the effective use of neural networks in solving general heat transfer problems of industrial complexity. Taken together, the results presented herein demonstrate that PINNs not only can solve ill-posed problems, which are beyond the reach of traditional computational methods, but they can also bridge the gap between computational and experimental heat transfer.

Engineering↗

The Pantheon+ Analysis: Forward Modeling the Dust and Intrinsic Color Distributions of Type Ia Supernovae, and Quantifying Their Impact on Cosmological Inferences

Abstract Recent studies have shown that the observed color distributions of Type Ia supernovae (SNe Ia) can be well described by a combination of a dust distribution and an intrinsic color distribution. Using the Pantheon+ sample of 1701 SN Ia, we apply a new forward-modeling fitting method (Dust2Dust) to measure the parent dust and color distributions, including their dependence on host-galaxy mass. At each fit step, the SN Ia selection efficiency is determined from a large simulated sample that is reweighted to reflect the proposed distributions. We use five separate metrics to describe the goodness of fit: distribution of fitted light-curve color c , cosmological residual trends with c , cosmological residual scatter with c , fitted color–luminosity relationship β SALT2 , and intrinsic scatter σ int . We present the results and the uncertainty in 12-dimensional space. Furthermore, we measure that the uncertainty on this modeling propagates to an upper threshold uncertainty in the equation of state of dark energy w of 0.014(1) for the Pantheon+ cosmology analysis and contributes negligible uncertainty to the Hubble constant H 0 . The Dust2Dust code is made publicly available at https://github.com/djbrout/dustdriver .

79 ASTRONOMY AND ASTROPHYSICS↗