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 271 records · Page 15

Sequential Bayesian inference of transition rates in the hidden Markov model for multi-state system degradation

The more easily available system performance data and advances in data analytics have provided us with opportunities to optimize maintenance programs for engineered systems, for example nuclear power plants. One key task in maintenance optimization is to obtain an accurate model for system degradation. Here, in this research, we propose a Bayesian method to address this problem. Noting that systems usually exhibit multiple states and that the actual state of a system usually is not directly observable, in the method we first model the system degradation process and the observation process based on a hidden Markov model. Then we develop a sequential Bayesian inference algorithm based on importance sampling and the forward algorithm to infer the posterior distributions of the transition rates in the hidden Markov model based on available observations. The proposed Bayesian method allows us to take advantage of evidence from multiple sources, and also allows us to perform Bayesian inference sequentially, without the need to use the entire history of observations every time new observations are collected. We demonstrate the proposed method using both synthetic data for a nuclear power plant feedwater pump and realistic data for a nuclear power plant chemistry analytical device.

42 ENGINEERING↗

Quantifying discordance in the 2015 Planck CMB spectrum

We examine the internal consistency of the Planck 2015 cosmic microwave background (CMB) temperature anisotropy power spectrum. We show that tension exists between cosmological constant cold dark matter (ΛCDM) model parameters inferred from multipoles ℓ<1000 (roughly those accessible to Wilkinson Microwave Anisotropy Probe), and from ℓ⩾1000, particularly the CDM density, Ω{sub c}h{sup 2}, which is discrepant at 2.5σ for a Planck -motivated prior on the optical depth, τ=0.07±0.02. We find some parameter tensions to be larger than previously reported because of inaccuracy in the code used by the Planck Collaboration to generate model spectra. The Planck ℓ⩾1000 constraints are also in tension with low-redshift data sets, including Planck ’s own measurement of the CMB lensing power spectrum (2.4σ), and the most precise baryon acoustic oscillation scale determination (2.5σ). The Hubble constant predicted by Planck from ℓ⩾1000, H{sub 0}=64.1±1.7 km s{sup −1} Mpc{sup −1}, disagrees with the most precise local distance ladder measurement of 73.0±2.4 km s{sup −1} Mpc{sup −1} at the 3.0σ level, while the Planck value from ℓ<1000, 69.7±1.7 km s{sup −1} Mpc{sup −1}, is consistent within 1σ. A discrepancy between the Planck and South Pole Telescope high-multipole CMB spectra disfavors interpreting these tensions as evidence for new physics. We conclude that the parameters from the Planck high-multipole spectrum probably differ from the underlying values due to either an unlikely statistical fluctuation or unaccounted-for systematics persisting in the Planck data.

79 ASTRONOMY AND ASTROPHYSICS↗

Virtual to Physical: Reinforcement Learning to Optimize SNS Particle Accelerator Controls

Complex accelerators must have control systems that can handle dynamic nonlinear environments. This makes traditional control methods unsuitable as they can struggle to adapt to these uncertainties. This provides an ideal environment for reinforcement learning algorithms as they are adaptable and generalizable. We present a reinforcement learning pipeline that can effectively handle the dynamics of a complex accelerator. We test and prove our pipelines capabilities on multiple environments including the Spallation Neutron Source (SNS) and the Beam Test Facility (BTF) at Oakridge National Lab (ORNL). Due to the limited time available to train an online algorithm like reinforcement learning on a real accelerator, we utilize a virtual twin accelerator (VIRAC) developed by ORNL to pretrain the policy and show its ability to converge in the virtual environment. We then test the adaptability of the pretrained RL model by applying it on the real accelerator and comparing the results. Utilizing our Scientific Optimization and Controls Toolkit (SOCT) and open-source standards such as Gymnasium we create and solve for a MEBT orbit correction problem in the SNS and an emittance maximization problem in the BTF. We show how Twin Delayed Deep Deterministic Policy Gradient (TD3) can solve this optimization environment in the virtual accelerator and transfer this policy onto the real accelerator for inference and model retraining. We show how reinforcement learning can be utilized as a control system for complex accelerators and provide a model pipeline for how an implementation performs and can be adapted to new accelerator control problems.

Kasparian, Armen [Thomas Jefferson National Accele↗

Multimodal parameter spaces of a complex multi-channel neuron model

One of the most common types of models that helps us to understand neuron behavior is based on the Hodgkin–Huxley ion channel formulation (HH model). A major challenge with inferring parameters in HH models is non-uniqueness: many different sets of ion channel parameter values produce similar outputs for the same input stimulus. Such phenomena result in an objective function that exhibits multiple modes (i.e., multiple local minima). This non-uniqueness of local optimality poses challenges for parameter estimation with many algorithmic optimization techniques. HH models additionally have severe non-linearities resulting in further challenges for inferring parameters in an algorithmic fashion. To address these challenges with a tractable method in high-dimensional parameter spaces, we propose using a particular Markov chain Monte Carlo (MCMC) algorithm, which has the advantage of inferring parameters in a Bayesian framework. The Bayesian approach is designed to be suitable for multimodal solutions to inverse problems. We introduce and demonstrate the method using a three-channel HH model. We then focus on the inference of nine parameters in an eight-channel HH model, which we analyze in detail. We explore how the MCMC algorithm can uncover complex relationships between inferred parameters using five injected current levels. The MCMC method provides as a result a nine-dimensional posterior distribution, which we analyze visually with solution maps or landscapes of the possible parameter sets. The visualized solution maps show new complex structures of the multimodal posteriors, and they allow for selection of locally and globally optimal value sets, and they visually expose parameter sensitivities and regions of higher model robustness. We envision these solution maps as enabling experimentalists to improve the design of future experiments, increase scientific productivity and improve on model structure and ideation when the MCMC algorithm is applied to experimental data.

97 MATHEMATICS AND COMPUTING↗

Determinants of photochemical characteristics of the photosynthetic electron transport chain of maize

The photosynthetic electron transport chain (ETC) is the bridge that links energy harvesting during the photophysical reactions at one end and energy consumption during the biochemical reactions at the other. Its functioning is thus fundamental for the proper balance between energy supply and demand in photosynthesis. Currently, there is a lack of understanding regarding how the structural properties of the ETC are affected by nutrient availability and plant developmental stages, which is a major roadblock to comprehensive modeling of photosynthesis. Redox parameters reflect the structural controls of ETC on the photochemical reactions and electron transport. We conducted joint measurements of chlorophyll fluorescence (ChlF) and gas exchange under systematically varying environmental conditions and growth stages of maize and sampled foliar nutrient contents. We utilized the recently developed steady-state photochemical model to infer redox parameters of electron transport from these measurements. We found that the inferred values of these photochemical redox parameters varied with leaf macronutrient content. These variations may be caused either directly by these nutrients being components of protein complexes on the ETC or indirectly by their impacts on the structural integrity of the thylakoid and feedback from the biochemical reactions. Also, the redox parameters varied with plant morphology and developmental stage, reflecting seasonal changes in the structural properties of the ETC. Our findings will facilitate the parameterization and simulation of complete models of photosynthesis.

59 BASIC BIOLOGICAL SCIENCES↗

Bootstrap inversion technique for atmospheric trace gas source detection and quantification using long open-path laser measurements

Advances in natural gas extraction technology have led to increased activity in the production and transport sectors in the United States and, as a consequence, an increased need for reliable monitoring of methane leaks to the atmosphere. We present a statistical methodology in combination with an observing system for the detection and attribution of fugitive emissions of methane from distributed potential source location landscapes such as natural gas production sites. We measure long (> 500 m), integrated open-path concentrations of atmospheric methane using a dual frequency comb spectrometer and combine measurements with an atmospheric transport model to infer leak locations and strengths using a novel statistical method, the non-zero minimum bootstrap (NZMB). The new statistical method allows us to determine whether the empirical distribution of possible source strengths for a given location excludes zero. Using this information, we identify leaking source locations (i.e., natural gas wells) through rejection of the null hypothesis that the source is not leaking. The method is tested with a series of synthetic data inversions with varying measurement density and varying levels of model–data mismatch. It is also tested with field observations of (1) a non-leaking source location and (2) a source location where a controlled emission of 3.1 &times; 10 -5 kg s -1 of methane gas is released over a period of several hours. This series of synthetic data tests and outdoor field observations using a controlled methane release demonstrates the viability of the approach for the detection and sizing of very small leaks of methane across large distances (4 + km 2 in synthetic tests). The field tests demonstrate the ability to attribute small atmospheric enhancements of 17 ppb to the emitting source location against a background of combined atmospheric (e.g., background methane variability) and measurement uncertainty of 5 ppb (1 σ ), when measurements are averaged over 2 min. The results of the synthetic and field data testing show that the new observing system and statistical approach greatly decreases the incidence of false alarms (that is, wrongly identifying a well site to be leaking) compared with the same tests that do not use the NZMB approach and therefore offers increased leak detection and sizing capabilities.

03 NATURAL GAS↗

3D characterisation and quantification of an offshore freshened groundwater system in the Canterbury Bight

Although offshore freshened groundwater (OFG) systems have been documented in numerous continental margins worldwide, their geometry, controls and emplacement dynamics remain poorly constrained. Here we integrate controlled-source electromagnetic, seismic reflection and borehole data with hydrological modelling to quantitatively characterise a previously unknown OFG system near Canterbury, New Zealand. The OFG system consists of one main, and two smaller, low salinity groundwater bodies. The main body extends up to 60 km from the coast and a seawater depth of 110 m. We attribute along-shelf variability in salinity to permeability heterogeneity due to permeable conduits and normal faults, and to recharge from rivers during sea level lowstands. A meteoric origin of the OFG and active groundwater migration from onshore are inferred. However, modelling results suggest that the majority of the OFG was emplaced via topographically-driven flow during sea level lowstands in the last 300 ka. Global volumetric estimates of OFG will be significantly revised if active margins, with steep coastal topographies like the Canterbury margin, are considered.

54 ENVIRONMENTAL SCIENCES↗

Weighted greedy-optimal design of computer experiments for kernel-based and Gaussian process model emulation and calibration

This article is concerned with the approximation of high-dimensional functions by kernel-based methods. Motivated by uncertainty quantification, which often necessitates the construction of approximations that are accurate with respect to a probability density function of random variables, we aim at minimizing the approximation error with respect to a weighted $L^p$-norm. We present a greedy procedure for designing computer experiments based upon a weighted modification of the pivoted Cholesky factorization. The method successively generates nested samples with the goal of minimizing error in regions of high probability. Numerical experiments validate that this new importance sampling strategy is superior to other sampling approaches, especially when used with non-product probability density functions. We also show how to use the proposed algorithm to efficiently generate surrogates for inferring unknown model parameters from data.

97 MATHEMATICS AND COMPUTING↗

Misleading or Falsification? Inferring Deceptive Strategies and Types in Online News and Social Media

Deceptive information in online news and social media has had dramatic effect on our society in recent years. This study is the first to gain deeper insights into writers' intent behind digital misinformation by analyzing psycholinguistic signals: moral foundations and connotations extracted from different types of deceptive news ranging from strategic disinformation to propaganda and hoaxes. To ensure consistency of our findings and generalizability across domains, we experiment with data from: (1) confirmed cases of disinformation in news summaries, (2) {propaganda}, hoax, and disinformation news pages, and (3) social media news. We first contrast lexical markers of biased language, syntactic and stylistic signals, and connotations across deceptive news types including disinformation, propaganda, and hoaxes, and {deceptive} strategies including misleading or falsification. We then incorporate these insights to build machine learning and deep learning predictive models to infer deception strategies and deceptive news types. Our experimental results demonstrate that unlike earlier work on deception detection, content combined with biased language markers, moral foundations, and connotations leads to better predictive performance of deception strategies compared to syntactic and stylistic signals (as reported in earlier work on deceptive reviews). Falsification strategy is easier to identify than misleading strategy. Disinformation is more difficult to predict than to propaganda or hoaxes. Deceptive news types (disinformation, propaganda, and hoaxes), unlike deceptive strategies (falsification and misleading), are more salient, and thus easier to identify in tweets than in news reports. Finally, our novel connotation analysis across deception types provides deeper understanding of writers' perspectives and therefore reveals the intentions behind digital misinformation.

deception↗

Organic matter sulfurization and organic carbon burial in the Mesoproterozoic

Throughout the Proterozoic Era, sedimentary organic carbon burial helped set the pace of global oxygenation and acted as a major modulator of atmospheric CO 2 and climate. Although Proterozoic rocks generally contain low concentrations of organic matter (OM), there are key exceptions to this rule, including the relatively OM-rich Arctic Bay shales from Baffin Island, Canada (Bylot Supergroup, Borden Basin, ~1.05 Ga). The mechanisms driving elevated OM concentrations in these and other Proterozoic shales remain poorly understood. In the Mesozoic and Cenozoic, organic matter sulfurization can be a major driver of enhanced OM burial across a range of redox conditions comparable to those inferred for many Proterozoic environments. Therefore, in this study, we evaluate the role of sulfurization in driving OM preservation in the Mesoproterozoic Borden Basin and discuss its relevance to Proterozoic systems in general. We present the first evidence for syngenetic-to-early-diagenetic OM sulfurization in a Proterozoic basin, which begins to fill a several-billion-year gap in our record of organic S across Earth history. We find that OM sulfurization was particularly extensive in shales from a relatively shallow-water section (Alpha River) but less extensive in shales deposited in deeper water (Shale Valley), which is consistent with models that infer sulfidic ‘wedges’ or O 2 -minimum-zone-type structures on shelf margins at least intermittently at this time. At the shallower site, organic S and pyrite are similarly 34 S-depleted and thus likely formed at roughly the same time near the sediment–water interface under conditions previously interpreted to have been ferruginous to intermittently sulfidic. In contrast, at the deeper-water site, large S-isotope differences between pyrite and organic S along with low apparent OM sulfurization intensities indicate that pyrite formation was favored over OM sulfurization during early sedimentation under variable but primarily ferruginous conditions. Although Mesoproterozoic biomass can be substantially sulfurized, indicators of sulfurization intensity are not correlated with OM concentrations, and therefore sulfurization does not appear to have been the primary driver of enhanced OM concentrations in Arctic Bay Formation shales. The link between sulfurization and total OM preservation may have been modulated during the deposition of Arctic Bay Formation shales by differences in iron availability, nutrient cycling, and particle dynamics in the Mesoproterozoic.

58 GEOSCIENCES↗

Cannabis monoterpene synthases: evaluating structure–function relationships

Terpene synthases catalyze the first committed step in the biosynthesis of terpenes, a structurally diverse class of natural products that also encompasses volatiles derived from precursors in the C10 to C15 range (termed monoterpenes and sesquiterpenes, respectively). In the review section of this article, we are providing information about all functionally characterized monoterpene synthases (MTSs) and sesquiterpene synthases (STSs) of Cannabis sativa L. We are also exploring the locations of MTSs and STSs in the chromosome-level assembly of the reference chemovar CBDRx. A follow-up computational structure–function analysis focuses on MTSs, as there is already a rich literature available on the topic. More specifically, by employing sequence comparisons and homology structural modeling, we infer which amino acid residues are likely to constrain the available space in the active site of cannabis MTSs. The emphasis of these studies was to investigate why some MTSs accept only a C10 diphosphate as substrate, while mixed MTS/STS enzymes also accommodate a C15 diphosphate. Here, by combining a literature review and computational analyses in a hybrid format, we are laying the foundation for future studies to better understand the determinants of substrate and product specificity in these fascinating enzymes.

59 BASIC BIOLOGICAL SCIENCES↗

Mixture Model Framework for Traumatic Brain Injury Prognosis Using Heterogeneous Clinical and Outcome Data

Prognoses of Traumatic Brain Injury (TBI) outcomes are neither easily nor accurately determined from clinical indicators. This is due in part to the heterogeneity of damage inflicted to the brain, ultimately resulting in diverse and complex outcomes. Using a data-driven approach on many distinct data elements may be necessary to describe this large set of outcomes and thereby robustly depict the nuanced differences among TBI patients’ recovery. In this work, we develop a method for modeling large heterogeneous data types relevant to TBI. Our approach is geared toward the probabilistic representation of mixed continuous and discrete variables with missing values. The model is trained on a dataset encompassing a variety of data types, including demographics, blood-based biomarkers, and imaging findings. In addition, it includes a set of clinical outcome assessments at 3, 6, and 12 months post-injury. The model is used to stratify patients into distinct groups in an unsupervised learning setting. We use the model to infer outcomes using input data, and show that the collection of input data reduces uncertainty of outcomes over a baseline approach. In addition, we quantify the performance of a likelihood scoring technique that can be used to self-evaluate the extrapolation risk of prognosis on unseen patients.

97 MATHEMATICS AND COMPUTING↗

Multiple RGB ortho-mosaics and digital surface models in 2017 and 2018 across the Lower Montane site in the East River Watershed, Colorado

Aerial imagery was collected at the Lower Montane site (Pumphouse) in the East River Watershed, Colorado during the spring, summer, and fall seasons of 2017 and 2018 to improve the understanding of seasonal vegetation dynamics and their drivers. The datasets include Red-Green-Blue (RGB) ortho-mosaics and digital surface models (DSMs) inferred from the Unoccupied Aerial System (UAS) acquired aerial RGB imagery for June 3, June 19, July 7, and August 14, 2017, and for March 14, April 26, June 1, June 18, July 6, and August 7, 2018. Real-Time Kinematic Global Positioning System (RTK-GPS) surveyed Ground control points (GCPs) were used to increase the reconstruction accuracy. The reconstructed RGB mosaics and DSMs have been trimmed to cover a similar spatial domain. The accuracy of the RGB mosaics is considered high (~10 cm). DSM accuracy is highest (~10 cm) where sufficient GCPS are available, and more difficult to assess elsewhere (see reconstruction reports for uncertainty estimates). The dataset includes a total of 20 GeoTIFF (.tif) files, 10 PDF (.pdf) files, 3 data CSV (.csv) files, and 2 metadata CSV (.csv) files. Feel free to contact the authors with any questions or collaboration interests.This work was supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231.

54 ENVIRONMENTAL SCIENCES↗

Kinematically consistent models of viscoelastic stress evolution

Abstract Following large earthquakes, coseismic stresses at the base of the seismogenic zone may induce rapid viscoelastic deformation in the lower crust and upper mantle. As stresses diffuse away from the primary slip surface in these lower layers, the magnitudes of stress at distant locations (>1 fault length away) may slowly increase. This stress relaxation process has been used to explain delayed earthquake triggering sequences like the 1992 M w = 7.3 Landers and 1999 M w = 7.1 Hector Mine earthquakes in California. However, a conceptual difficulty associated with these models is that the magnitudes of stresses asymptote to constant values over long time scales. This effect introduces persistent perturbations to the total stress field over many earthquake cycles. Here we present a kinematically consistent viscoelastic stress transfer model where the total perturbation to the stress field at the end of the earthquake cycle is zero everywhere. With kinematically consistent models, hypotheses about the potential likelihood of viscoelastically triggered earthquakes may be based on the timing of stress maxima, rather than on any arbitrary or empirically constrained stress thresholds. Based on these models, we infer that earthquakes triggered by viscoelastic earthquake cycle effects may be most likely to occur during the first 50% of the earthquake cycle regardless of the assumed long‐term and transient viscosities.

DeVries, Phoebe M. R.↗

Tunable Microwave Conductance of Nanodomains in Ferroelectric PbZr 0.2 Ti 0.8 O 3 Thin Film

Ferroelectric materials exhibit spontaneous polarization that can be switched by electric field. Beyond traditional applications as nonvolatile capacitive elements, the interplay between polarization and electronic transport in ferroelectric thin films has enabled a path to neuromorphic device applications involving resistive switching. A fundamental challenge, however, is that finite electronic conductivity may introduce considerable power dissipation and perhaps destabilize ferroelectricity itself. In this work, tunable microwave frequency electronic response of domain walls injected into ferroelectric lead zirconate titanate (PbZr 0.2 Ti 0.8 O 3 ) on the level of a single nanodomain is revealed. Tunable microwave response is detected through first-order reversal curve spectroscopy combined with scanning microwave impedance microscopy measurements taken near 3 GHz. Contributions of film interfaces to the measured AC conduction through subtractive milling, where the film exhibited improved conduction properties after removal of surface layers, are investigated. Using statistical analysis and finite element modeling, we inferred that the mechanism of tunable microwave conductance is the variable area of the domain wall in the switching volume. These observations open the possibilities for ferroelectric memristors or volatile resistive switches, localized to several tens of nanometers and operating according to well-defined dynamics under an applied field.

36 MATERIALS SCIENCE↗

Adoption of switchgrass cultivation for biofuel under uncertainty: A discrete-time modeling approach

Production of biofuels from cellulosic sources, such as switchgrass, is being encouraged through man- dates, incentives, and subsidies. However, uncertainty in future prices coupled with large establishment costs often inhibit their cultivation. Owing to their inability to incorporate uncertainty and dynamic decision-making, standard discounted cash flow techniques are ineffective for analyzing such in- vestments. We formulate a discrete-time binomial framework to model output prices, allowing us to incorporate price uncertainty, stand age, and variable crop yields into the analytical framework. We analyze the feasibility of investments in switchgrass cultivation under varying price transition paths, evaluate the relationship between risk and profitability, and estimate the value of flexible decision- making options wherein the farmer can alter cultivation choices. Here, we find that switchgrass cultivation is only 32% likely to be profitable in the base model and infer that on-farm management could play an important role in entry and exit decisions. We also find that subsidies are important for project viability and policymakers could consider incorporating payments for ecosystem services to encourage adoption.

09 BIOMASS FUELS↗

Bedrock architecture, soil texture, and hyporheic zone characterization combining electrical resistivity and induced polarization imaging

The shallow (<15 m depth) subsurface environment of a short reach along a bedrock stream was investigated with electrical resistivity and induced polarization (IP) to map details of bedrock and soil sediments. The bedrock is mostly comprised of limestone, and is generally resistive compared to the overlying soil. The soil-bedrock interface was determined through a trial and error approach using a sharp boundary feature in the inversion model. The inferred bedrock surface determined from the inversions exhibited undulating patterns with troughs and ridges. A near continuous trough ran alongside the stream within the floodplain, and is suggestive of a paleochannel among other interpretations of this feature. The structure of the electrical resistivity above the bedrock showed small-scale elongated features. The chargeability from the IP method showed larger scale features. High values of chargeability were associated with the sediments in the floodplain, and low values were associated with bedrock, stream, and soil on the elevated banks above the stream. If the chargeability is associated with membrane polarization characteristic of clayey soils, then IP seems to highlight the mere existence of clay, while the resistivity may be more discernable of the relative proportion of clay. Vegetation differences may also explain the chargeability distribution, where parts of the survey with high chargeability had dense pine with no understory making the soils more organically rich.

54 ENVIRONMENTAL SCIENCES↗