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 937 records · Page 52

Seismic Waveform Inversion Capability on Resource-Constrained Edge Devices

Seismic full wave inversion (FWI) is a widely used non-linear seismic imaging method used to reconstruct subsurface velocity images, however it is time consuming, has high computational cost and depend heavily on human interaction. Recently, deep learning has accelerated it’s use in several data-driven techniques, however most deep learning techniques suffer from overfitting and stability issues. In this work, we propose an edge computing-based data-driven inversion technique based on supervised deep convolutional neural network to accurately reconstruct the subsurface velocities. Deep learning based data-driven technique depends mostly on bulk data training. In this work, we train our deep convolutional neural network (DCN) (UNet and InversionNet) on the raw seismic data and their corresponding velocity models during the training phase to learn the non-linear mapping between the seismic data and velocity models. The trained network is then used to estimate the velocity models from new input seismic data during the prediction phase. The prediction phase is performed on a resource-constrained edge device such as Raspberry Pi. Raspberry Pi provides real-time and on-device computational power to execute the inference process. In addition, we demonstrate robustness of our models to perform inversion in the presence on noise by performing both noise-aware and no-noise training and feeding the resulting trained models with noise at different signal-to-noise (SNR) ratio values. We make great efforts to achieve very feasible inference times on the Raspberry Pi for both models. Specifically, the inference times per prediction for UNet and InversionNet models on Raspberry Pi were 22 and 4 s respectively whilst inference times for both models on the GPU were 2 and 18 s which are very comparable. Finally, we have designed a user-friendly interactive graphical user interface (GUI) to automate the model execution and inversion process on the Raspberry Pi.

Manu, Daniel (ORCID:0000000154982677)↗

A Comparison of Model- and Satellite-Derived Aerosol Optical Depth and Reflectivity

The determination of an accurate quantitative understanding of the role of tropospheric aerosols in the earth's radiation budget is extremely important because forcing by anthropogenic aerosols presently represents one of the most uncertain aspects of climate models. Here the authors present a systematic comparison of three different analyses of satellite-retrieved aerosol optical depth based on the Advanced Very High Resolution Radiometer (AVHRR)- measured radiances with optical depths derived from six different models. Also compared are the model-derived clear-sky reflected shortwave radiation with satellite-measured reflectivities derived from the Earth Radiation Budget Experiment (ERBE) satellite. The three different satellite-derived optical depths differ by between -0.10 and 0.07 optical depth units in comparison to the average of the three analyses depending on latitude and month, but the general features of the retrievals are similar. The models differ by between -0.09 and +0.16 optical depth units from the average of the models. Differences between the average of the models and the average of the satellite analyses range over -0.11 to +0.05 optical depth units. These differences are significant since the annual average clear-sky radiative forcing associated with the difference between the average of the models and the average of the satellite analyses ranges between -3.9 and 0.7 W m(exp -2) depending on latitude and is -1.7 W m (exp -2) on a global average annual basis. Variations in the source strengths of dimethylsulfide (DMS)-derived aerosols and sea salt aerosols can explain differences between the models, and between the models and satellite retrievals of up to 0.2 optical depth units. The comparison of model-generated reflected shortwave radiation and ERBE-measured shortwave radiation is similar in character as a function of latitude to the analysis of modeled and satellite-retrieved optical depths, but the differences between the modeled clear-sky reflected flux and the ERBE clear-sky reflected flux is generally larger than that inferred from the difference between the models and the AVHRR optical depths, especially at high latitudes. The difference between the mean of the models and the ERBE-analyzed clear-sky flux is 1.6 W m(exp -2). The overall comparison indicates that the model-generated aerosol optical depth is systematically lower than that inferred from measurements between the latitudes of 10 and 30 deg S. It is not likely that the shortfall is due to small values of the sea salt optical depth because increases in this component would create modeled optical depths that are larger than those from satellites in the region north of 30 deg N and near 50 deg S. Instead, the source strengths for DMS and biomass aerosols in the models may be too low. Firm conclusions, however, will require better retrieval procedures for the satellites, including better cloud screening procedures, further improvement of the model's treatment of aerosol transport and removal, and a better determination of aerosol source strengths.

Penner, Joyce E.↗

Physics-tailored machine learning reveals unexpected physics in dusty plasmas

Dusty plasma is a mixture of ions, electrons, and macroscopic charged particles that is commonly found in space and planetary environments. The particles interact through Coulomb forces mediated by the surrounding plasma, and as a result, the effective forces between particles can be nonconservative and nonreciprocal. Machine learning (ML) models are a promising route to learn these complex forces, yet their structure should match the underlying physical constraints to provide useful insight. Here, we demonstrate and experimentally validate an ML approach that incorporates physical intuition to infer force laws in a laboratory dusty plasma. Trained on 3D particle trajectories, the model accounts for inherent symmetries, nonidentical particles, and learns the effective nonreciprocal forces between particles with exquisite accuracy (R 2 > 0.99). We validate the model by inferring particle masses in two independent yet consistent ways. The model’s accuracy enables precise measurements of particle charge and screening length, identifying large deviations from common theoretical assumptions. Our ability to identify unknown physics from experimental data demonstrates how ML-powered approaches can guide new routes of scientific discovery in many-body systems. Furthermore, we anticipate our ML approach to be a starting point for inferring laws from dynamics in a wide range of many-body systems, from colloids to living organisms.

Science & Technology - Other Topics↗

Ringed seal demography in a changing climate

Abstract Climate change is affecting species’ distributions and abundances worldwide. Baseline population estimates, against which future observations may be compared, are necessary if we are to detect ecological change. Arctic sea ice ecosystems are changing rapidly and we lack baseline population estimates for many ice‐associated species. Provided we can detect them, changes in Arctic marine ecosystems may be signaled by changes in indicator species such as ringed seals ( Pusa hispida ). Ringed seal monitoring has provided estimates of survival and fertility rates, but these have not been used for population‐level inference. Using matrix population models, we synthesized existing demographic parameters to obtain estimates of historical ringed seal population growth and structure in Amundsen Gulf and Prince Albert Sound, Canada. We then formalized existing hypotheses about the effects of emerging environmental stressors (i.e., earlier spring ice breakup and reduced snow depth) on ringed seal pup survival. Coupling the demographic model to ice and snow forecasts available from the Coupled Model Intercomparison Project resulted in projections of ringed seal population size and structure up to the year 2100. These projections showed median declines in population size ranging from 50% to 99%. Corresponding to these projected declines were substantial changes in population structure, with increasing proportions of ringed seal pups and adults and declining proportions of juveniles. We explored if currently collected, harvest‐based data could be used to detect the projected changes in population stage structure. Our model suggests that at a present sample size of 100 seals per year, the projected changes in stage structure would only be reliably detected by mid‐century, even for the most extreme climate models. This modeling process revealed inconsistencies in existing estimates of ringed seal demographic rates. Mathematical population models such as these can contribute both to understanding past population trends as well as predicting future ones, both of which are necessary if we are to detect and interpret future observations.

Reimer, Jody R.↗

Electroweak constraints from the COHERENT experiment

We compute bounds on coefficients of effective operators in the Standard Model that can be inferred from observations of neutrino scattering by the COHERENT experiment. While many operators are bound extremely well by past experiments the full future data set from COHERENT will provide modest improvements for some operators.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Strength, deformation, and equation of state of tungsten carbide to 66 GPa

Here, strength, texture, and equation of state of hexagonal tungsten monocarbide (WC) have been determined under quasi-hydrostatic and non-hydrostatic compression to 66 GPa using angle-dispersive X-ray diffraction in the diamond anvil cell. Quasi-hydrostatic compression in a Ne pressure medium demonstrates that nanocrystalline WC is slightly less incompressible than bulk-scale WC, with respective bulk moduli of K 0 = 377 ± 7 and 397 ± 7 GPa and pressure derivatives K 0 ’ = 3.8 ± 0.3 and 3.7 ±0.3. This decrease in incompressibility with grain size is similar to behavior observed in other ceramics. Under nonhydrostatic compression, WC supports a mean differential stress of ~12-15 GPa at plastic yielding, which occurs at ~30 GPa. Strength in WC is anisotropic, with the (001) plane supporting 29-42% higher stress than stresses calculated from mean strain. Simulations using an Elasto-ViscoPlastic Self-Consistent model indicate that strength inferred from lattice strain theory may be overestimated due to effects of plastic deformation. Plastic deformation generates a texture maximum near $\langle\bar{2}110\rangle$ in the compression orientation, initially through prismatic slip on the {$10\bar{1}0$}$\langle\bar{1}2\bar{1}0\rangle$ and {$10\bar{1}0$}$\langle0001\rangle$ slip systems, followed by activation of pyramidal slip on {$10\bar{1}1$}$\langle\bar{2}113\rangle$ at ~40-50 GPa.

36 MATERIALS SCIENCE↗

Molecular dynamics simulations of cesium adsorption on illite nanoparticles

Here we present that the charged surfaces of micaceous minerals, especially illite, regulate the mobility of the major radioisotopes of Cs ( 134 Cs, 135 Cs, 137 Cs) in the geosphere. Despite the long history of Cs adsorption studies, the nature of the illite surface sites remains incompletely understood. To address this problem, we present atomistic simulations of Cs competition with Na for three candidate illite adsorption sites – edge, basal plane, and interlayer. Our simulation results are broadly consistent with affinities and selectivities that have been inferred from surface complexation models. Cation exchange on the basal planes is thermodynamically ideal, but exchange on edge surfaces and within interlayers shows complex, thermodynamically non-ideal behavior. The basal planes are weakly Cs-selective, while edges and interlayers have much higher affinity for Cs. The dynamics of Na-Cs exchange are rapid for both cations on the basal planes, but considerably slower for Cs localized on edge surfaces. Lastly, in addition to new insights into Cs adsorption and exchange with Na on illite, we report the development of a methodology capable of simulating fully-flexible clay mineral nanoparticles with stable edge surfaces using a well-tested interatomic potential model.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Unraveling Spatio-Temporal Chemistry Evolution in Laser Ablation Plumes and Its Relation to Initial Plasma Conditions

The chemistry evolution in a laser ablation plume depends strongly on its initial physical conditions. In this article, we investigate the impact of plasma generation conditions on the interrelated phenomena of expansion dynamics, plasma chemistry, and physical conditions. Plasmas are produced from a uranium metal target in air using nanosecond, femtosecond, and femtosecond filament-assisted laser ablation. Time-resolved two-dimensional spectral imaging was performed for evaluating the spatio-temporal evolu-tion of atoms, diatoms, polyatomic molecules, and nanoparticles in situ. Emission spectral features reveal that molecular formation occurs at early times in both femtosecond and filament ablation plumes, although with different temporal decay. In contrast, mo-lecular formation is found to occur at much later times in nanosecond plasma evolution. Spectral modeling is used to infer tem-poral behavior of plasma excitation temperature, and results highlighted that nanosecond plume temperatures decay rapidly while plumes generated by filament ablation show slowest decay. In this work, we find U atoms and UO molecules co-exist in ultrafast laser produced plasmas even at early times after plasma onset owing to favorable temperatures for molecular formation. Regardless of irradiation conditions, plume emission features showed presence of higher oxides (i.e. UxOy), although with different temporal histories. Our study provides insight into the impact of plasma generation conditions on chemistry evolution in plasmas produced from traditional focused femtosecond, nanosecond, and filament-assisted laser ablation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Arctic Cloud‐Base Ice Precipitation Properties Retrieved Using Bayesian Inference

Abstract Cloud‐climate feedbacks are still the greatest source of uncertainty in current climate projections. Arctic clouds, which are predominantly stratiform and supercooled, often long‐lived, and nearly continuously precipitate ice particles, contribute roughly 10% of the uncertainty attributed to the global cloud feedback. This Arctic cloud uncertainty is driven by incomplete observational and theoretical knowledge required to estimate and explain the state and active processes occurring in those clouds. A focus on ice precipitation properties at Arctic cloud base rather than the surface deconfounds the product of cloud condensate sink processes from the influence of the atmospheric thermodynamic state below cloud base, rendering cloud‐base properties a more appealing target for inference and evaluation of model simulations. Here I describe an inverse model for the estimation of cloud base ice precipitation properties over Utqiagvik, North Slope of Alaska, using the synthesis of ground‐based radar and lidar measurements. By leveraging a Markov Chain Monte Carlo algorithm as the core of the inverse model, a wide range of particle size distributions are sampled, and different combinations of ice habit models are examined, both of which are typically fixed in other retrieval methods. Results show intriguing links between different cloud base thermodynamic and ice precipitation properties. Apparent ice number concentration enhancements at temperatures of −5 and −15°C suggest possible secondary ice production (SIP). The analysis alludes to an overestimation of SIP occurrence and intensity, especially in studies relying only on radar or lidar measurements. Finally, reflectivity‐dependent ice precipitation rate and ice water content parameterizations are presented.

54 ENVIRONMENTAL SCIENCES↗

Constraining Bedrock Groundwater Residence Times in a Mountain System With Environmental Tracer Observations and Bayesian Uncertainty Quantification

Groundwater residence time distributions provide fundamental insights on the hydrological processes within watersheds. Yet, observations that can constrain groundwater residence times over broad timescales remain scarce in mountain catchment studies. We use environmental tracers (CFC-12, SF 6 , 3 H, and 4 He) to investigate groundwater residence times along a hillslope in the East River Watershed, Colorado, USA. We develop a Bayesian inference framework that applies a Markov-chain Monte Carlo (MCMC) approach to estimate noble gas recharge temperature, elevation, and excess-air parameters and the resulting environmental tracer concentrations. MCMC is then used to propagate the environmental tracer uncertainties to estimates of groundwater mean residence times inferred with lumped parameter models. All samples contain 3 H, CFC-12, and SF 6 in addition to terrigenic 4 He, suggesting a mixture of water characterized by modern and premodern residence times. 4He exponential mean residence times range from hundreds of years at the upslope well to thousands of years at the toe-slope well assuming average crustal production rates. We find that binary mixing residence time distributions with separate young and old mixing fractions are needed to predict the 4 He, CFC-12, SF 6 , and 3 H observations, supporting the importance of flow path mixing in this bedrock system. Our findings that the fractured bedrock hosts groundwater with a mixture of residence times ranging from decades to millennia suggest variable recharge dynamics and flow path mixing along the hillslope and highlight the importance of characterizing groundwater systems with observations that are sensitive to transport over a broad range of residence times.

54 ENVIRONMENTAL SCIENCES↗

Early aqueous activity on the ordinary and carbonaceous chondrite parent bodies recorded by fayalite

Here, chronology of aqueous activity on chondrite parent bodies constrains their accretion times and thermal histories. Radiometric 53 Mn– 53 Cr dating has been successfully applied to aqueously formed carbonates in CM carbonaceous chondrites. Owing to the absence of carbonates in ordinary (H, L and LL), and CV and CO carbonaceous chondrites, and the lack of proper standards, there are no reliable ages of aqueous activity on their parent bodies. Here we report the first 53 Mn– 53 Cr ages of aqueously formed fayalite in the L3 chondrite Elephant Moraine 90161 as 2.4 +1.8 -1.3 Myr after calcium–aluminium-rich inclusions (CAIs), the oldest Solar System solids. In addition, measurements using our synthesized fayalite standard show that fayalite in the CV3 chondrite Asuka 881317 and CO3-like chondrite MacAlpine Hills 88107 formed and 4.2 +0.8 -0.7 Myr after CAIs, respectively. Thermal modelling, combined with the inferred conditions (temperature and water/rock ratio) and 53 Mn– 53 Cr ages of aqueous alteration, suggests accretion of the L, CV and CO parent bodies ~1.8–2.5 Myr after CAIs.

79 ASTRONOMY AND ASTROPHYSICS↗

Hardware-accelerated inference for real-time gravitational-wave astronomy

The field of transient astronomy has seen a revolution with the first gravitational-wave detections and the arrival of multi-messenger observations they enabled. Transformed by the first detection of binary black hole and binary neutron star mergers, computational demands in gravitational-wave astronomy are expected to grow by at least a factor of two over the next five years as the global network of kilometer-scale interferometers are brought to design sensitivity. With the increase in detector sensitivity, real-time delivery of gravitational-wave alerts will become increasingly important as an enabler of multi-messenger followup. In this work, we report a novel implementation and deployment of deep learning inference for real-time gravitational-wave data denoising and astrophysical source identification. This is accomplished using a generic Inference-as-a-Service model that is capable of adapting to the future needs of gravitational-wave data analysis. Overall, our implementation allows seamless incorporation of hardware accelerators and also enables the use of commercial or private (dedicated) as-a-service computing. Based on our results, we propose a paradigm shift in low-latency and offline computing in gravitational-wave astronomy. Such a shift can address key challenges in peak-usage, scalability and reliability, and provide a data analysis platform particularly optimized for deep learning applications. The achieved sub-millisecond scale latency will also be relevant for any machine learning-based real-time control systems that may be invoked in the operation of near-future and next generation ground-based laser interferometers, as well as the front-end collection, distribution and processing of data from such instruments.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Electron-beam energy reconstruction for neutrino oscillation measurements

Neutrinos exist in one of three types or ‘flavours’—electron, muon and tau neutrinos—and oscillate from one flavour to another when propagating through space. This phenomena is one of the few that cannot be described using the standard model of particle physics (reviewed in ref. 1), and so its experimental study can provide new insight into the nature of our Universe (reviewed in ref. 2). Neutrinos oscillate as a function of their propagation distance (L) divided by their energy (E). Therefore, experiments extract oscillation parameters by measuring their energy distribution at different locations. As accelerator-based oscillation experiments cannot directly measure E, the interpretation of these experiments relies heavily on phenomenological models of neutrino–nucleus interactions to infer E. Here we exploit the similarity of electron–nucleus and neutrino–nucleus interactions, and use electron scattering data with known beam energies to test energy reconstruction methods and interaction models. We find that even in simple interactions where no pions are detected, only a small fraction of events reconstruct to the correct incident energy. More importantly, widely used interaction models reproduce the reconstructed energy distribution only qualitatively and the quality of the reproduction varies strongly with beam energy. This shows both the need and the pathway to improve current models to meet the requirements of next-generation, high-precision experiments such as Hyper-Kamiokande (Japan) and DUNE (USA).

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Multi-fidelity information fusion with concatenated neural networks

Recently, computational modeling has shifted towards the use of statistical inference, deep learning, and other data-driven modeling frameworks. Although this shift in modeling holds promise in many applications like design optimization and real-time control by lowering the computational burden, training deep learning models needs a huge amount of data. This big data is not always available for scientific problems and leads to poorly generalizable data-driven models. This gap can be furnished by leveraging information from physics-based models. Exploiting prior knowledge about the problem at hand, this study puts forth a physics-guided machine learning (PGML) approach to build more tailored, effective, and efficient surrogate models. For our analysis, without losing its generalizability and modularity, we focus on the development of predictive models for laminar and turbulent boundary layer flows. In particular, we combine the self-similarity solution and power-law velocity profile (low-fidelity models) with the noisy data obtained either from experiments or computational fluid dynamics simulations (high-fidelity models) through a concatenated neural network. We illustrate how the knowledge from these simplified models results in reducing uncertainties associated with deep learning models applied to boundary layer flow prediction problems. The proposed multi-fidelity information fusion framework produces physically consistent models that attempt to achieve better generalization than data-driven models obtained purely based on data. While we demonstrate our framework for a problem relevant to fluid mechanics, its workflow and principles can be adopted for many scientific problems where empirical, analytical, or simplified models are prevalent. In line with grand demands in novel PGML principles, this work builds a bridge between extensive physics-based theories and data-driven modeling paradigms and paves the way for using hybrid physics and machine learning modeling approaches for next-generation digital twin technologies.

42 ENGINEERING↗

Three-dimensional simulations of turbulent convective mixing in ONe and CO classical nova explosions

Classical novae are thermonuclear explosions that take place in the envelopes of accreting white dwarfs in binary systems. The material piles up under degenerate conditions, driving a thermonuclear runaway. The energy released by the suite of nuclear processes operating at the envelope heats the material up to peak temperatures of ~(1-4) × 10 8 K. During these events, about 10 -3 -10 -7 M ⊙ , enriched in CNO and, sometimes, other intermediate-mass elements (e.g., Ne, Na, Mg, Al) are ejected into the interstellar medium. To account for the gross observational properties of classical novae (in particular, the high concentrations of metals spectroscopically inferred in the ejecta), models require mixing between the (solar-like) material transferred from the secondary and the outermost layers (CO- or ONe-rich) of the underlying white dwarf. Recent multidimensional simulations have demonstrated that Kelvin-Helmholtz instabilities can naturally produce self-enrichment of the accreted envelope with material from the underlying white dwarf at levels that agree with observations. However, the feasibility of this mechanism has been explored in the framework of CO white dwarfs, while mixing with different substrates still needs to be properly addressed. We performed three-dimensional simulations of mixing at the core-envelope interface during nova outbursts with the multidimensional code FLASH, for two types of substrates: CO- and ONe-rich. We also show that the presence of an ONe-rich substrate, as in “neon novae”, yields higher metallicity enhancements in the ejecta than CO-rich substrates (i.e., non-neon novae). Finally, a number of requirements and constraints for such 3D simulations (e.g., minimum resolution, size of the computational domain) are also outlined.

79 ASTRONOMY AND ASTROPHYSICS↗

Quantum learning of classical stochastic processes: The completely positive realization problem

Among several tasks in Machine Learning, a specially important one is the problem of inferring the latent variables of a system and their causal relations with the observed behavior. A paradigmatic instance of this is the task of inferring the hidden Markov model underlying a given stochastic process. This is known as the positive realization problem (PRP), [L. Benvenuti and L. Farina, IEEE Trans. Autom. Control 49(5), 651–664 (2004)] and constitutes a central problem in machine learning. The PRP and its solutions have far-reaching consequences in many areas of systems and control theory, and is nowadays an important piece in the broad field of positive systems theory. We consider the scenario where the latent variables are quantum (i.e., quantum states of a finite-dimensional system) and the system dynamics is constrained only by physical transformations on the quantum system. The observable dynamics is then described by a quantum instrument, and the task is to determine which quantum instrument — if any — yields the process at hand by iterative application. We take as a starting point the theory of quasi-realizations, whence a description of the dynamics of the process is given in terms of linear maps on state vectors and probabilities are given by linear functionals on the state vectors. This description, despite its remarkable resemblance with the hidden Markov model, or the iterated quantum instrument, is however devoid of any stochastic or quantum mechanical interpretation, as said maps fail to satisfy any positivity conditions. The completely positive realization problem then consists in determining whether an equivalent quantum mechanical description of the same process exists. We generalize some key results of stochastic realization theory, and show that the problem has deep connections with operator systems theory, giving possible insight to the lifting problem in quotient operator systems. Our results have potential applications in quantum machine learning, device-independent characterization and reverse-engineering of stochastic processes and quantum processors, and more generally, of dynamical processes with quantum memory [M. Guţă, Phys. Rev. A 83(6), 062324 (2011); M. Guţă and N. Yamamoto, e-print http://arxiv.org/abs/1303.3771 (2013)].

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Effects of preheat and mix on the fuel adiabat of an imploding capsule

We demonstrate the effect of preheat, hydrodynamic mix and vorticity on the adiabat of the deuterium-tritium (DT) fuel in fusion capsule experiments. We show that the adiabat of the DT fuel increases resulting from hydrodynamic mixing due to the phenomenon of entropy of mixture. An upper limit of mix, M clean =M DT ≥ 0:98 is found necessary to keep the DT fuel on a low adiabat. We demonstrate in this study that the use of a high adiabat for the DT fuel in theoretical analysis and with the aid of 1D code simulations could explain some aspects of 3D effects and mix in capsule implosion. Furthermore, we can infer from our physics model and the observed neutron images the adiabat of the DT fuel in the capsule and the amount of mix produced on the hot spot.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗