Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “likelihood”

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 793 records · Page 44

OVERVIEW OF PARTICLE DEPOSITION MODELS FOR SPENT NUCLEAR FUEL STORAGE SYSTEMS

Deposition models were built to evaluate contaminant deposition on spent nuclear fuel (SNF) canisters. The primary contaminant of concern is chloride, which is dispersed in the atmosphere and then deposits onto the canisters. During dry storage, the primary degradation process is likely to be Chloride Induced Stress Corrosion Cracking (CISCC) at the heat-affected zones of the canister welds. It is known that stainless steel canisters are susceptible to CISCC; however, the rate of chloride deposition onto the canisters is poorly known, based on sparse field data from a small number of sites. The models presented in this study could be useful for determining the rate of deposition on the canisters and the likelihood of CISCC to help with SNF canister ageing management. The deposition models were developed with the commercial computational fluid dynamics (CFD) code STAR-CCM+. Various deposition mechanisms were considered and incorporated into the models, and a sensitivity study was conducted to determine the most important mechanisms for deposition within a SNF storage system. The models included both a vertical and horizontal configuration storage system: NAC International’s Modular, Advanced Generation, Nuclear All-purpose STORage System (MAGNASTOR®) and a NUHOMS® horizontal storage module respectively. The resulting canister deposition on the horizontal canister is visually compared with inspection data taken onsite at the Calvert Cliffs Nuclear Power Plant. These models are preliminary, and development of the models will continue. Future validation exercises are currently being planned, including the Canister Deposition Field Demonstration (CDFD) effort funded by U.S. Department of Energy office of Nuclear Energy office of Spent Fuel Waste Science and Technology. The goal of the modeling presented is to demonstrate a potential modeling technique that could be used to plan and inform SNF canister ageing management programs with predictive models for the timing and occurrence of canister CISCC.

Suffield, Sarah R.↗

Surface science of shape-selective metal nanocrystal synthesis from first-principles: Growth of Cu nanowires and nanocubes

The authors present results from first-principles density functional theory aimed at understanding the aqueous solution-phase growth of fivefold twinned copper nanowires and single-crystal nanocubes capped by hexadecylamine (HDA). The role of solution-phase chloride, present in the Cu salt or as an additive, is emphasized. Using ab initio thermodynamics, the authors delineate the range of solution-phase conditions, characterized by the chemical potentials of chloride and HDA, under which Cu nanowires and nanocubes can be grown. The authors discuss the likelihood of thermodynamic and/or kinetic nanostructures for various solution-phase concentrations. Their results are in good agreement with experiments and indicate that methods and insights developed for surface science in gas-phase or vacuum conditions can yield much insight into liquid-phase systems.

36 MATERIALS SCIENCE↗

Infrared thermography NDT for in-situ defect detection in sandwich composite panel manufacturing

Composite manufacturing presents numerous challenges, as defects can arise from various sources throughout the process. In sandwich composite structures, the integration of a foam core introduces additional complexity and increases the likelihood of defect formation like delamination. To mitigate these issues and reduce the risk of future structural failures, in-situ monitoring during manufacturing is essential. This study investigates infrared (IR) thermography as a non-destructive technique for detecting manufacturing defects in foam-core sandwich composite panels under thermally excited conditions representative of in-situ processing. A stationary FLIR A8590 IR camera (640 × 512 pixels, 30Hz, 17mm lens, 9 ft stand-off distance) was used to monitor prefabricated panels subjected to controlled external heating simulating compression molding and resin cure exotherm. Interlaminar delamination defects with characteristic sizes ranging from 0.25 × 0.25in² to 5 × 5in² produced measurable surface temperature depressions of approximately 4–10°C during transient cooling, exceeding the effective noise floor of the camera by more than two standard deviations. Thicker laminates exhibited prolonged defect detectability windows due to increased thermal diffusion time. In contrast, embedded Teflon inclusions generated weak thermal contrasts of ≤ 3°C, approaching the measurement noise floor, due to limited thermal property contrast with the surrounding glass fiber composite. These results establish quantitative detectability limits for stationary thermographic inspection of sandwich composite panels under manufacturing-representative thermal cycles.

Barakat, Abdallah [ORNL] (ORCID:0000000296141398)↗

Evaluating the mineral commodity supply risk of the U.S. manufacturing sector

Trade tensions, resource nationalism, and various other factors are increasing concerns regarding the supply reliability of nonfuel mineral commodities. This is especially the case for commodities required for new and emerging technologies ranging from electric vehicles to wind turbines. In this analysis, we use a conventional risk-modeling framework to develop and apply a new methodology for assessing the supply risk to the U.S. manufacturing sector. Specifically, supply risk is defined as the confluence of three factors: the likelihood of a foreign supply disruption, the dependency of U.S. manufacturers on foreign supplies, and the ability of U.S. manufacturers to withstand a supply disruption. The methodology is applied to 52 commodities for the decade spanning 2007–2016. The results indicate that a subset of 23 commodities, including cobalt, niobium, rare earth elements, and tungsten, pose the greatest supply risk. This supply risk is dynamic, shifting with changes in global market conditions.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

On the role of α-alumina in the origin of life: Surface-driven assembly of amino acids

We investigate the hypothesis that mineral/water interfaces played a crucial catalytic role in peptide formation by promoting the self-assembly of amino acids. Using classical molecular dynamics simulations, we demonstrate that the α-alumina(0001) surface exhibits an affinity of 4k B Tfor individual glycine or GG dipeptide molecules due to hydrogen bonds. In simulations with multiple glycine molecules, surface-bound glycine enhances further adsorption, leading to the formation of long chains connected by hydrogen bonds between the carboxyl and amine groups of glycine molecules. We find that the likelihood of observing chains longer than 10 glycine units increases by at least five orders of magnitude at the surface compared to the bulk. This surface-driven assembly is primarily due to local high density and alignment with the alumina surface pattern. Together, these results propose a model for how mineral surfaces can induce configuration-specific assembly of amino acids, thereby promoting condensation reactions.

Science & Technology - Other Topics↗

Bayesian estimation of HIV acquisition dates for prevention trials

Accurate timing estimates of when participants acquire HIV in HIV prevention trials are necessary for determining antibody levels at acquisition. The Antibody-Mediated Prevention (AMP) Studies showed that a passively administered broadly neutralizing antibody can prevent the acquisition of HIV from a neutralization-sensitive virus. We developed a pipeline for estimating the date of detectable HIV acquisition (DDA) in AMP Study participants using diagnostic and viral sequence data. Using a Bayesian strategy that combines three streams of data (REN [rev/vpu/env/Δnef] sequence, GP [gag/Δpol] sequence, and diagnostic) where their 95% credible intervals overlap based on pre-specified criteria and decision rules. We evaluated the performance of our AMP pipeline using PacBio viral sequence data from 41 participants across two prospective acute HIV acquisition cohort studies, FRESH and RV217, with twice-weekly sampling. These cohort studies enrolled young women in South Africa and men and women in Kenya and Thailand, respectively, with a high likelihood of HIV acquisition. In evaluating performance, “true DDA” was the center of bounds between last-negative and first-positive RNA diagnostic tests (median time 4 days, range 2–7 days); bias was the mean difference between estimated and true DDA. Using diagnostic data alone yielded timing estimates with a bias of 2.4 days and root mean square error (RMSE) of 7.9 days. These results were improved using sequence + diagnostic data (bias 1.5 days, RMSE 6.9 days), as well as by restricting sequence-based estimation to samples from ≤5 weeks post-DDA (bias 0.2 days, RMSE 7.8 days).

59 BASIC BIOLOGICAL SCIENCES↗

Catabolic pathway acquisition by rhizosphere bacteria readily enables growth with a root exudate component but does not affect root colonization

Horizontal gene transfer (HGT) is a fundamental evolutionary process that plays a key role in bacterial evolution. The likelihood of a successful transfer event is expected to depend on the precise balance of costs and benefits resulting from pathway acquisition. Most experimental analyses of HGT have focused on phenotypes that have large fitness benefits under appropriate selective conditions, such as antibiotic resistance. However, many examples of HGT involve phenotypes that are predicted to provide smaller benefits, such as the ability to catabolize additional carbon sources. We have experimentally simulated the consequences of one such HGT event in the laboratory, studying the effects of transferring a pathway for catabolism of the plant-derived aromatic compound salicyl alcohol between rhizosphere isolates from the Pseudomonas genus. We find that pathway acquisition enables rapid catabolism of salicyl alcohol with only minor disruptions to the existing metabolic and regulatory networks of the new host. However, this new catabolic potential does not confer a measurable fitness advantage during competitive growth in the rhizosphere. We conclude that the phenotype of salicyl alcohol catabolism is readily transferable but is selectively neutral under environmentally relevant conditions. We propose that this condition is common and that HGT of many pathways will be self-limiting because the selective benefits are small.

59 BASIC BIOLOGICAL SCIENCES↗

Genetic Predictive Factors for Nonsusceptible Phenotypes and Multidrug Resistance in Expanded-Spectrum Cephalosporin-Resistant Uropathogenic Escherichia coli from a Multicenter Cohort: Insights into the Phenotypic and Genetic Basis of Coresistance

Antimicrobial resistance in urinary tract infections (UTIs) is a major public health concern. This study aims to characterize the phenotypic and genetic basis of multidrug resistance (MDR) among expanded-spectrum cephalosporin-resistant (ESCR) uropathogenic Escherichia coli (UPEC) causing UTIs in California patient populations. Between February and October 2019, 577 ESCR UPEC isolates were collected from patients at 6 clinical laboratory sites across California. Lineage and antibiotic resistance genes were determined by analysis of whole-genome sequence data. The lineages ST131, ST1193, ST648, and ST69 were predominant, representing 46%, 5.5%, 4.5%, and 4.5% of the collection, respectively. Overall, 527 (91%) isolates had an expanded-spectrum β-lactamase (ESBL) phenotype, with bla CTX-M-15 , bla CTX-M-27 , bla CTX-M-55 , and bla CTX-M-14 being the most prevalent ESBL genes. In the 50 non-ESBL phenotype isolates, 40 (62%) contained bla CMY-2 , which was the predominant plasmid-mediated AmpC (pAmpC) gene. Narrow-spectrum β-lactamases, bla TEM-1B and bla OXA-1 , were also found in 44.9% and 32.1% of isolates, respectively. Among ESCR UPEC isolates, isolates with an ESBL phenotype had a 1.7-times-greater likelihood of being MDR than non-ESBL phenotype isolates (P < 0.001). The cooccurrence of bla CTX-M-15 , bla OXA-1 , and aac(6')-Ib-cr within ESCR UPEC isolates was strongly correlated. Cooccurrence of bla CTX-M-15 , bla OXA-1 , and aac(6')-Ib-cr was associated with an increased risk of nonsusceptibility to piperacillin-tazobactam, cefepime, fluoroquinolones, and amikacin as well as MDR. Multivariate regression revealed the presence of bla CTX-M-55 , bla TEM-1B , and the ST131 genotype as predictors of MDR.

59 BASIC BIOLOGICAL SCIENCES↗

The Microbial Community and Functional Potential in the Midland Basin Reveal a Community Dominated by Both Thiosulfate and Sulfate-Reducing Microorganisms

The Permian Basin is the highest producing oil and gas reservoir in the United States. Hydrocarbon resources in this region are often accessed by unconventional extraction methods, including horizontal drilling and hydraulic fracturing. Despite the importance of the Permian Basin, there is no publicly available microbiological data from this region. We completed an analysis of Permian produced water samples to understand the dynamics present in hydraulically fractured wells in this region. We analyzed produced water samples taken from 10 wells in the Permian region of the Midland Basin using geochemical measurements, 16S rRNA gene sequencing, and metagenomic sequencing. Compared to other regions, we found that Permian Basin produced water was characterized by higher sulfate and lower total dissolved solids (TDS) concentrations, with a median of 1,110 mg/L and 107,000 mg/L. Additionally, geochemical measurements revealed the presence of frac hits, or interwell communication events where an established well is affected by the pumping of fracturing fluid into a new well. The occurrence of frac hits was supported by correlations between the microbiome and the geochemical parameters. Our 16S rRNA gene sequencing identified a produced water microbiome characterized by anaerobic, halophilic, and sulfur reducing taxa. Interestingly, sulfate and thiosulfate reducing taxa including Halanaerobium, Orenia, Marinobacter, and Desulfohalobium were the most prevalent microbiota in most wells. We further investigated the metabolic potential of microorganisms in the Permian Basin with metagenomic sequencing. We recovered 15 metagenome assembled genomes (MAGs) from seven different samples representing 6 unique well sites. These MAGs corroborated the high presence of sulfate and thiosulfate reducing genes across all wells, especially from key taxa including Halanaerobium and Orenia. The observed microbiome composition and metabolic capabilities in conjunction with the high sulfate concentrations demonstrate a high potential for hydrogen sulfide production in the Permian Basin. Additionally, evidence of frac hits suggests the possibility for the exchange of microbial cells and/or genetic information between wells. This exchange would increase the likelihood of hydrogen sulfide production and has implications for the oil and gas industry.

16S RNA↗

Optimal Experimental Design for Inverse Problems in the Presence of Observation Correlations

Optimal experimental design (OED) is the general formalism of sensor placement and decisions on the data collection strategy for engineered or natural experiments. This problem is prevalent in many critical fields such as battery design, numerical weather prediction, geosciences, environmental and urban studies. State-of-the-art computational methods for experimental design do not accommodate correlation structure in observational errors produced by many expensive-to-operate devices such as X-ray machines, radars, and satellites. Discarding evident data correlations leads to biased results, higher expenses, and waste of valuable resources. We present a general formulation of the OED formalism for model-constrained large-scale Bayesian linear inverse problems, where measurement errors are generally correlated. The proposed approach utilizes the Hadamard product of matrices to formulate the weighted-likelihood, and is valid for both finite as well as infinite-dimensional Bayesian inverse problems. Extensive numerical experiments are carried out for empirical verification of the proposed approach using an advection-diffusion model, where the objective is to optimally place a small set of sensors, under a limited budget, to predict the concentration of a contaminant in a closed and bounded domain.

Correlated observations↗

Enabling Hyper-Differential Sensitivity Analysis for Ill-Posed Inverse Problems

Inverse problems constrained by partial differential equations (PDEs) play a critical role in model development and calibration. In many applications, there are multiple uncertain parameters in a model that must be estimated. However, high dimensionality of the parameters and computational complexity of the PDE solves make such problems challenging. A common approach is to reduce the dimension by fixing some parameters (which we will call auxiliary parameters) to a best estimate and use techniques from PDE-constrained optimization to estimate the other parameters. In this article, hyper-differential sensitivity analysis (HDSA) is used to assess the sensitivity of the solution of the PDE-constrained optimization problem to changes in the auxiliary parameters. Foundational assumptions for HDSA require satisfaction of the optimality conditions which are not always practically feasible as a result of ill-posedness in the inverse problem. Here we introduce novel theoretical and computational approaches to justify and enable HDSA for ill-posed inverse problems by projecting the sensitivities on likelihood informed subspaces and defining a posteriori updates. Our proposed framework is demonstrated on a nonlinear multiphysics inverse problem motivated by estimation of spatially heterogeneous material properties in the presence of spatially distributed parametric modeling uncertainties.

97 MATHEMATICS AND COMPUTING↗

Generalized Bayesian MARS: Tools for Stochastic Computer Model Emulation

The multivariate adaptive regression spline (MARS) approach of Friedman and its Bayesian counterpart are effective approaches for the emulation of computer models. The traditional assumption of Gaussian errors limits the usefulness of MARS, and many popular alternatives, when dealing with stochastic computer models. Here, we propose a generalized Bayesian MARS (GBMARS) framework which admits the broad class of generalized hyperbolic distributions as the induced likelihood function. This allows us to develop tools for the emulation of stochastic simulators which are parsimonious, scalable, and interpretable and require minimal tuning, while providing powerful predictive and uncertainty quantification capabilities. GBMARS is capable of robust regression with t distributions, quantile regression with asymmetric Laplace distributions, and a general form of “Normal-Wald” regression in which the shape of the error distribution and the structure of the mean function are learned simultaneously. We demonstrate the effectiveness of GBMARS on various stochastic computer models, and we show that it compares favorably to several popular alternatives.

97 MATHEMATICS AND COMPUTING↗

Score-Based Physics-Informed Neural Networks for High-Dimensional Fokker–Planck Equations

The Fokker-Planck (FP) equation is a foundational partial differential equation (PDE) in stochastic processes involving Brownian motions. However, the curse of dimensionality (CoD) poses a formidable challenge when dealing with high-dimensional FP equations. Although Monte Carlo simulation and (vanilla) Physics-Informed Neural Networks (PINNs) have shown the potential to tackle CoD, both methods exhibit significant numerical errors in high dimensions when dealing with the probability density function (PDF) associated with Brownian motion. The point-wise PDF values tend to decrease exponentially as dimensionality increases, surpassing the precision of numerical simulations and resulting in substantial errors. In addition, due to its massive sampling, Monte Carlo fails to offer fast sampling. Modeling the logarithm likelihood (LL) via vanilla PINNs transforms the FP equation into a notoriously difficult Hamilton-Jacobi-Bellman (HJB) equation, which is impractical for PINN learning, whose error grows rapidly with dimension. To this end, we propose a novel approach utilizing a score-based solver to fit the score function in stochastic differential equations (SDEs). The score function, defined as the gradient of the LL, plays a fundamental role in inferring LL and PDF and enables fast SDE sampling, offering an effective means to overcome the CoD. Three fitting methods, Score Matching (SM), Sliced Score Matching (SSM), and Score-PINN, are introduced, each contributing unique advantages in computational complexity, accuracy, and generality. The proposed score-based SDE solver operates in two stages: first, employing score matching or Score-PINN to acquire the score function; and second, solving the LL via an ordinary differential equation (ODE) using the obtained score function. Comparative evaluations across these methods showcase varying trade-offs. The proposed methodology is evaluated across diverse SDEs, including anisotropic Ornstein-Uhlenbeck processes, geometric Brownian motion, and Brownian motion with varying eigenspace. We also test various distributions, including Gaussian, Log-normal, Laplace, and Cauchy distributions. The numerical results demonstrate the score-based SDE solver’s stability, speed, and performance across different experimental settings, solidifying its potential as a solution to CoD for high-dimensional FP equations.

97 MATHEMATICS AND COMPUTING↗

Enhancing Gaussian Process Surrogates for Optimization and Posterior Approximation via Random Exploration

This paper proposes novel noise-free Bayesian optimization strategies that rely on a random exploration step to enhance the accuracy of Gaussian process surrogate models. The new algorithms retain the ease of implementation of the classical GP-UCB algorithm, but the additional random exploration step accelerates their convergence, nearly achieving the optimal convergence rate. Furthermore, to facilitate Bayesian inference with intractable likelihoods, we propose to utilize optimization iterates for maximum a posteriori estimation to build a Gaussian process surrogate model for the unnormalized log-posterior density. We provide bounds for the Hellinger distance between the true and the approximate posterior distributions in terms of the number of design points. We demonstrate the effectiveness of our Bayesian optimization algorithms in nonconvex benchmark objective functions, in a machine learning hyperparameter tuning problem, and in a black-box engineering design problem. The effectiveness of our posterior approximation approach is demonstrated in two Bayesian inference problems for parameters of dynamical systems.

Bayesian inference↗

Efficient Neural Network Approaches for Conditional Optimal Transport with Applications in Bayesian Inference

In this work, we present two neural network approaches that approximate the solutions of static and dynamic conditional optimal transport (COT) problems. Both approaches enable conditional sampling and conditional density estimation, which are core tasks in Bayesian inference—particularly in the simulation-based (“likelihood-free”) setting. Our methods represent the target conditional distribution as a transformation of a tractable reference distribution. Obtaining such a transformation, chosen here to be an approximation of the COT map, is computationally challenging even in moderate dimensions. To improve scalability, our numerical algorithms use neural networks to parameterize candidate maps and further exploit the structure of the COT problem. Our static approach approximates the map as the gradient of a partially input convex neural network. It uses a novel numerical implementation to increase computational efficiency compared to state-of-the-art alternatives. Our dynamic approach approximates the conditional optimal transport via the flow map of a regularized neural ODE; compared to the static approach, it is slower to train but offers more modeling choices and can lead to faster sampling. We demonstrate both algorithms numerically, comparing them with competing state-of-the-art approaches, using benchmark datasets and simulation-based Bayesian inverse problems.

97 MATHEMATICS AND COMPUTING↗

Discovery of Probabilistic Dirichlet-to-Neumann Maps on Graphs

Dirichlet-to-Neumann maps enable the coupling of multiphysics simulations across computational subdomains by ensuring continuity of state variables and fluxes at artificial interfaces. We present a novel method for learning Dirichlet-to-Neumann maps on graphs using Gaussian processes, specifically for problems where the data obey a conservation law arising from an underlying partial differential equation. Our approach combines discrete exterior calculus and nonlinear optimal recovery to infer relationships between vertex and edge values. This framework yields data-driven predictions with uncertainty quantification across the entire graph, even when observations are limited to a subset of vertices and edges. By minimizing the reproducing kernel Hilbert space norm while penalizing kernel complexity through maximum likelihood estimation, our method ensures that the resulting surrogate strictly enforces conservation laws without overfitting. We demonstrate our method on two representative applications: subsurface flow in fracture networks and arterial blood flow. Finally, the results demonstrate that the method maintains high accuracy and well-calibrated uncertainty estimates even under severe data scarcity, highlighting its potential for scientific applications where limited data and reliable uncertainty quantification are critical.

Dirichlet-to-Neumann map↗

Evidence for $\text {W}\text {W}$ production from double-parton interactions in proton–proton collisions at $\sqrt{s} = 13 \,\text {TeV} $ = 13TeV

A search for WW production from double-parton scattering processes using same-charge electron-muon and dimuon events is reported, based on proton-proton collision data collected at a center-of-mass energy of 13 TeV. The analyzed data set corresponds to an integrated luminosity of 77.4 fb$^{-1}$, collected using the CMS detector at the LHC in 2016 and 2017. Multivariate classifiers are used to discriminate between the signal and the dominant background processes. A maximum likelihood fit is performed to extract the signal cross section. This leads to the first evidence for WW production via double-parton scattering, with a significance of 3.9 standard deviations. The measured inclusive cross section is 1.41 $\pm$ 0.28 (stat) $\pm$ 0.28 (syst) pb.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

A multi-dimensional search for new heavy resonances decaying to boosted $\text{W}$ $\text{W}$, $\text{W}$ $\text{Z}$, or $\text{Z}$ $\text{Z}$ boson pairs in the dijet final state at 13 $\text {Te}\text{V}$

A search in an all-jet final state for new massive resonances decaying to $\text{ W }{}{}\text{ W }{}{}$, $\text{ W }{}{}\text{ Z }{}{}$, or $\text{ Z }{}{}\text{ Z }{}{}$ boson pairs using a novel analysis method is presented. The analysis is performed on data corresponding to an integrated luminosity of 77.3 $\,\text {fb}^{-1}$ recorded with the CMS experiment at the LHC at a centre-of-mass energy of 13 $\text {Te}\text {V}$. The search is focussed on potential narrow-width resonances with masses above 1.2 $\text {Te}\text {V}$, where the decay products of each $\text{ W }{}{}$ or $\text{ Z }{}{}$ boson are expected to be collimated into a single, large-radius jet. The signal is extracted using a three-dimensional maximum likelihood fit of the two jet masses and the dijet invariant mass, yielding an improvement in sensitivity of up to 30% relative to previous search methods. No excess is observed above the estimated standard model background. In a heavy vector triplet model, spin-1 ${\text {Z}}^{\prime }$ and ${\text {W}}^{\prime }$ resonances with masses below 3.5 and 3.8 $\text {Te}\text {V}$, respectively, are excluded at 95% confidence level. In a bulk graviton model, upper limits on cross sections are set between 27 and 0.2 $\,\text {fb}$ for resonance masses between 1.2 and 5.2 $\text {Te}\text {V}$, respectively. The limits presented in this paper are the best to date in the dijet final state.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗