Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “parametric uncertainty”

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 523 records · Page 29

Romie: A Domain-Independent Tool for Computer-Aided Robust Operations Management

Romie is a decision support tool based on AI's latest advances in the domain of robust scheduling. Unlike all its predecessors, the tool allows to (i) visually model the operational problem and context entirely (ii) optimize to find near-optimal schedules while taking uncertainty into account and deals with (iii) a combination of various {key performance indicators (KPIs). It comes with a web user interface. Part or all of the modelled activities may be associated to random variables describing their stochastic durations, in order to produce schedules that are robust w.r.t. temporal uncertainty. Hence, depending on the pursued KPIs, the schedules maximize a combination of the following terms: the probability of satisfying the problem constraints, the expected return/efficiency, the expected outcome quality, and even the operators' wellness by minimizing its expected extra-hours. Initially developed for spatial exploration and demonstration in the context of Mars analog missions, this versatile tool is here applied to operations management in both biotechnology manufacturing and robots parametrization in a cave exploration context.

Chien, Steve A.↗

Compound Continental Risk of Multiple Extreme Floods in the United States

Abstract Understanding spatially correlated floods and modeling joint hazard associated with threshold exceedances across multiple locations is crucial for accurate estimation of continental‐scale portfolio risk. This work uses a non‐parametric copula‐based spatial simulator to analyze peak floods across the United States to derive the first‐of‐its‐kind continental portfolio risk estimates at the 10‐ and 100‐year return levels. We find significant interdependence in floods across the nation, revealing the recurring pattern of extreme events affecting the Northeast, Central, West, and Northwest United States in the same year. The stochastic simulator effectively manages high‐dimensional data and offers reliable uncertainty estimates for both spatially dependent floods and the aggregated flood losses at the continental level. El Niño‐Southern Oscillation and Atlantic Multidecadal Oscillation are identified as statistically significant tele‐connectors of aggregate loss. This research aims to advance the understanding of compound continental flood hazard and the potential large‐scale climate teleconnections that lead to such compound floods.

Geology↗

Multi-modelling predictions show high uncertainty of required carbon input changes to reach a 4‰ target

Soils store vast amounts of carbon (C) on land, and increasing soil organic carbon (SOC) stocks in already managed soils such as croplands may be one way to remove C from the atmosphere, thereby limiting subsequent warming. The main objective of this study was to estimate the amount of additional C input needed to annually increase SOC stocks by 4‰ at 16 long-term agricultural experiments in Europe, including exogenous organic matter (EOM) additions. We used an ensemble of six SOC models and ran them under two configurations: (1) with default parametrization and (2) with parameters calibrated site-by-site to fit the evolution of SOC stocks in the control treatments (without EOM). We compared model simulations and analysed the factors generating variability across models. The calibrated ensemble was able to reproduce the SOC stock evolution in the unfertilised control treatments. We found that, on average, the experimental sites needed an additional 1.5 ± 1.2 Mg C ha -1 year -1 to increase SOC stocks by 4‰ per year over 30 years, compared to the C input in the control treatments (multi-model median ± median standard deviation across sites). That is, a 119% increase compared to the control. While mean annual temperature, initial SOC stocks and initial C input had a significant effect on the variability of the predicted C input in the default configuration (i.e., the relative standard deviation of the predicted C input from the mean), only water-related variables (i.e., mean annual precipitation and potential evapotranspiration) explained the divergence between models when calibrated. Here, our work highlights the challenge of increasing SOC stocks in agriculture and accentuates the need to increasingly lean on multi-model ensembles when predicting SOC stock trends and related processes. To increase the reliability of SOC models under future climate change, we suggest model developers to better constrain the effect of water-related variables on SOC decomposition.

4 per 1000 initiative↗

Bayesian Analysis of TRISO Fuel: Quantifying Model Inadequacy, Incorporating Lower-Length-Scale Effects, and Developing Parallel Active Learning Capabilities

The U.S. Department of Energy (DOE)’s Nuclear Energy Advanced Modeling and Simulation (NEAMS) program aims to develop predictive capabilities by applying computational methods to the analysis and design of advanced reactor and fuel-cycle systems. This program has been providing engineering-scale support for the continued development of BISON, a high-fidelity, high-resolution fuel performance tool. Fuel behavior in nuclear reactors is governed by a complex network of mechanisms that interact with various other physics aspects in the reactor system. Any model developed to represent fuel behavior will likely be idealized, resulting in uncertainties when comparing their predictions against the observed data. In Fiscal Year (FY)-23, we initiated the Uncertainty Quantification (UQ) work by using Bayesian methods to establish a level of model trustworthiness and further improve it, with a particular emphasis on TRI-Structural isOtropic (TRISO) nuclear fuel. This year, we further expanded on that UQ work by investigating an approach to quantifying model inadequacy and accounting for lower-length scale (LLS) effects in TRISO silver (Ag) release modeling. Furthermore, we are implementing parallel active learning capabilities to reduce the computational cost (i.e., required computational resources and elapsed time) of performing UQ. Specifically, we utilized The Kennedy O’Hagan framework for Bayesian uncertainty quantification (KOH) to account for model inadequacy in TRISO Ag release predictions made by BISON. The KOH framework represents an improvement over the standard Bayesian framework used in FY-23. Explicitly accounting for model inadequacy in the Bayesian framework helps establish the level of experimental noise uncertainty in the Advanced Gas Reactor (AGR) data. We compared the inverse UQ results obtained from both the standard Bayesian and KOH frameworks in light of the AGR-2/3/4 data, and also compared the predictive UQ results obtained from these two frameworks in light of the AGR-1 data. Next, we investigated the impact of considering LLS effects in the Ag release simulations. We developed an expanded database of LLS simulated effective diffusivities for Ag, covering a wide range of microstructures and temperatures. Using this database, we developed a framework for incorporating LLS effects into the engineering-scale Ag release UQ. We developed both parametric and non-parametric approaches for bridging the length scales. We then investigated the inverse UQ results in light of the AGR-2/3/4 data and the predictive UQ results in light of the AGR-1 data, and compared the LLS-informed approach and the Arrhenius equation, which does not include microstructure information. Finally, we discussed implementing parallel active learning capabilities in the Multiphysics Object Oriented Simulation Environment (MOOSE)/BISON to reduce the computational cost (i.e., computational resources and elapsed time) of Bayesian UQ. For verification purposes, we first tested these new capabil ities on a species interaction problem. We then demonstrated them on the TRISO Ag release application, showing that parallel active learning capabilities can enhance the accuracy of UQ while also substantially reducing the computational cost in comparison to the reference methods developed in FY-23.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

A Data-Driven Nonparametric Approach for Probabilistic Load-Margin Assessment Considering Wind Power Penetration

A modern power system is characterized by an increasing penetration of wind power, which results in large uncertainties in its states. These uncertainties must be quantified properly; otherwise, the system security may be threatened. Facing this challenge, here we propose a cost-effective, data-driven approach to assessing a power system's load margin probabilistically. Using actual wind data, a kernel density estimator is applied to infer the nonparametric wind speed distributions, which are further merged into the framework of a vine copula. The latter enables us to simulate complex multivariate and highly dependent model inputs with a variety of bivariate copulae that precisely represent the tail dependence in the correlated samples. Furthermore, to reduce the prohibitive computational time of traditional Monte-Carlo simulations that process a large amount of samples, we propose to use a nonparametric, Gaussian-process-emulator-based reduced-order model to replace the original complicated continuation power-flow model through a Bayesian-learning framework. To accelerate the convergence rate of this Bayesian algorithm, a truncated polynomial chaos surrogate, which serves as a highly efficient, parametric Bayesian prior, is developed. This emulator allows us to execute the time-consuming continuation power-flow solver at the sampled values with a negligible computational cost. Results of simulations that are performed on several test systems reveal the impressive performance of the proposed method in the probabilistic load-margin assessment.

17 WIND ENERGY↗

The diphoton $$q_T$$ spectrum at N$$^3$$LL$$^\prime $$ + NNLO

Abstract We present a $$q_T$$ q T -resummed calculation of diphoton production at order N $$^3$$ 3 LL $$^\prime $$ ′ + NNLO. To reach the primed level of accuracy we have implemented the recently published three-loop $${\mathcal {O}}(\alpha _s^3)$$ O ( α s 3 ) virtual corrections in the $$q\bar{q}$$ q q ¯ channel and the three-loop transverse momentum dependent beam functions and combined them with the existing infrastructure of , a code performing resummation at order N $${}^3$$ 3 LL. While the primed predictions are parametrically not more accurate, one typically observes from lower orders and other processes that they are the dominant effect of the next order. We include in both the $$q\bar{q}$$ q q ¯ and loop-induced gg channel the hard contributions consistently together at order $$\alpha _s^3$$ α s 3 and find that the resummed $$q\bar{q}$$ q q ¯ channel without matching stabilizes indeed. Due to large matching corrections and large contributions and uncertainties from the gg channel, the overall improvements are small though. We furthermore study the effect of hybrid-cone photon isolation and hard-scale choice on our fully matched results to describe the ATLAS $$\text {8}$$ 8 TeV data and find that the hybrid-cone isolation worsens agreement at small $$q_T$$ q T compared to smooth-cone isolation.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Semileptonic form factors for $$B\rightarrow D^*\ell \nu $$ at nonzero recoil from $$2+1$$-flavor lattice QCD: Fermilab Lattice and MILC Collaborations

Abstract We present the first unquenched lattice-QCD calculation of the form factors for the decay $$B\rightarrow D^*\ell \nu $$ B → D ∗ ℓ ν at nonzero recoil. Our analysis includes 15 MILC ensembles with $$N_f=2+1$$ N f = 2 + 1 flavors of asqtad sea quarks, with a strange quark mass close to its physical mass. The lattice spacings range from $$a\approx 0.15$$ a ≈ 0.15 fm down to 0.045 fm, while the ratio between the light- and the strange-quark masses ranges from 0.05 to 0.4. The valence b and c quarks are treated using the Wilson-clover action with the Fermilab interpretation, whereas the light sector employs asqtad staggered fermions. We extrapolate our results to the physical point in the continuum limit using rooted staggered heavy-light meson chiral perturbation theory. Then we apply a model-independent parametrization to extend the form factors to the full kinematic range. With this parametrization we perform a joint lattice-QCD/experiment fit using several experimental datasets to determine the CKM matrix element $$|V_{cb}|$$ | V cb | . We obtain $$\left| V_{cb}\right| = (38.40 \pm 0.68_{\text {th}} \pm 0.34_{\text {exp}} \pm 0.18_{\text {EM}})\times 10^{-3}$$ V cb = ( 38.40 ± 0 . 68 th ± 0 . 34 exp ± 0 . 18 EM ) × 10 - 3 . The first error is theoretical, the second comes from experiment and the last one includes electromagnetic and electroweak uncertainties, with an overall $$\chi ^2\text {/dof} = 126/84$$ χ 2 /dof = 126 / 84 , which illustrates the tensions between the experimental data sets, and between theory and experiment. This result is in agreement with previous exclusive determinations, but the tension with the inclusive determination remains. Finally, we integrate the differential decay rate obtained solely from lattice data to predict $$R(D^*) = 0.265 \pm 0.013$$ R ( D ∗ ) = 0.265 ± 0.013 , which confirms the current tension between theory and experiment.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Precision measurement of the electron energy-loss function in tritium and deuterium gas for the KATRIN experiment

Abstract The KATRIN experiment is designed for a direct and model-independent determination of the effective electron anti-neutrino mass via a high-precision measurement of the tritium $$\upbeta $$ β -decay endpoint region with a sensitivity on $$m_\nu $$ m ν of 0.2 $$\hbox {eV}/\hbox {c}^2$$ eV / c 2 (90% CL). For this purpose, the $$\upbeta $$ β -electrons from a high-luminosity windowless gaseous tritium source traversing an electrostatic retarding spectrometer are counted to obtain an integral spectrum around the endpoint energy of 18.6 keV. A dominant systematic effect of the response of the experimental setup is the energy loss of $$\upbeta $$ β -electrons from elastic and inelastic scattering off tritium molecules within the source. We determined the energy-loss function in-situ with a pulsed angular-selective and monoenergetic photoelectron source at various tritium-source densities. The data was recorded in integral and differential modes; the latter was achieved by using a novel time-of-flight technique. We developed a semi-empirical parametrization for the energy-loss function for the scattering of 18.6-keV electrons from hydrogen isotopologs. This model was fit to measurement data with a 95% $$\hbox {T}_2$$ T 2 gas mixture at 30 K, as used in the first KATRIN neutrino-mass analyses, as well as a $$\hbox {D}_2$$ D 2 gas mixture of 96% purity used in KATRIN commissioning runs. The achieved precision on the energy-loss function has abated the corresponding uncertainty of $$\sigma (m_\nu ^2)< {{10}^{-2}}{\hbox {eV}^{2}}$$ σ ( m ν 2 ) < 10 - 2 eV 2 [1] in the KATRIN neutrino-mass measurement to a subdominant level.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Theory uncertainties of the irreducible background to VBF Higgs production

Higgs boson production through gluon fusion in association with two jets is an irreducible background to Higgs boson production through vector boson fusion, one of the most important channels for analyzing and understanding the Higgs boson properties at the Large Hadron Collider. Despite a range of available simulation tools, precise predictions for the corresponding final states are notoriously hard to achieve. Using state-of-the-art fixed-order calculations as the baseline for a comparison, we perform a detailed study of similarities and differences in existing event generators. We provide consistent setups for the simulations that can be used to obtain identical parametric precision in various programs used by experiments. We find that NLO calculations for the two-jet final state are essential to achieve reliable predictions.

Chen, Xuan [Shandong U.] (ORCID:0000000167722196)↗

Jülich-Bonn-Washington model for pion electroproduction multipoles

Pion electroproduction off the proton is analyzed in a new framework based on a general parametrization of transition amplitudes, including constraints from gauge invariance and threshold behavior. Data with energies 1.13 GeV < W < 1.6 GeV and Q 2 below 6 Ge V 2 are included. The model is an extension of the latest Jülich-Bonn solution incorporating constraints from pion-induced and photoproduction data. Here, performing large-scale fits (≈10 5 data) we find a set of solutions with $χ^2_{\text{dof}}$ = 1.69–1.81 which allows us to assess the systematic uncertainty of the approach.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

The Problem of Size in Robust Design

To facilitate the effective solution of multidisciplinary, multiobjective complex design problems, a departure from the traditional parametric design analysis and single objective optimization approaches is necessary in the preliminary stages of design. A necessary tradeoff becomes one of efficiency vs. accuracy as approximate models are sought to allow fast analysis and effective exploration of a preliminary design space. In this paper we apply a general robust design approach for efficient and comprehensive preliminary design to a large complex system: a high speed civil transport (HSCT) aircraft. Specifically, we investigate the HSCT wing configuration design, incorporating life cycle economic uncertainties to identify economically robust solutions. The approach is built on the foundation of statistical experimentation and modeling techniques and robust design principles, and is specialized through incorporation of the compromise Decision Support Problem for multiobjective design. For large problems however, as in the HSCT example, this robust design approach developed for efficient and comprehensive design breaks down with the problem of size - combinatorial explosion in experimentation and model building with number of variables -and both efficiency and accuracy are sacrificed. Our focus in this paper is on identifying and discussing the implications and open issues associated with the problem of size for the preliminary design of large complex systems.

Koch, Patrick N.↗

Radiative Corrections to Superallowed β Decays in Effective Field Theory

The accuracy of V u d determinations from superallowed β decays critically hinges on control over radiative corrections. Recently, substantial progress has been made on the single-nucleon, universal corrections, while nucleus-dependent effects, typically parametrized by a quantity δ NS , are much less well constrained. Here, we lay out a program to evaluate this correction from effective field theory (EFT), highlighting the dominant terms as predicted by the EFT power counting. Moreover, we compare the results to a dispersive representation of δ NS and show that the expected momentum scaling applies even in the case of low-lying intermediate states. Our EFT framework paves the way toward calculations of δ NS and thereby addresses the dominant uncertainty in V u d . Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

A Reference Ocean Surface Emission and Backscatter Model from Microwaves to Infrared

Satellite observations are vital for the initialization of Numerical Weather Prediction models, and very important for climate monitoring and prediction, as well as other applications such as hydrology and flood awareness prediction. Knowledge of radiative contributions from the Earth's surface is needed to sound the lower troposphere from space. The lack of a reference quality ocean emission and backscatter model is a major gap in our ability to provide absolute calibration of the satellite based observing system. Uncertainty in emissivity models is not well characterized and different models are used for different spectral bands, for active and passive instruments. An International Space Science Institute (ISSI) team was put together [4] to address these issues. The objectives of the team are to provide a reference model as a community software (i.e., documented and freely available code), that is maintained and supported, has traceable uncertainty estimations, and that enables new science from microwaves to infrared with bidirectional reflectance distribution function (BRDF) capability. We will present the model and its various components, discussing the choices between various parameterizations, building on the LOCEAN model of [2]. The model predictions will be evaluated at various frequencies, including comparisons to radiometric observations by SMAP, AMSR2 and GMI (e.g., [5]). We will discuss early model evaluation in the infrared and for active microwave sensors. Areas of ongoing research include improving the foam parametrization (coverage and emissivity) to provide consistent performances across frequencies, building on [1], and the azimuthal dependence of the active and passive signals. The model will be used to generate training data for fast models e.g., Fastem, [3], that are used in operational data assimilation and climate re-analysis.

Emmanuel Dinnat↗

A composite likelihood approach for inference under photometric redshift uncertainty

ABSTRACT Obtaining accurately calibrated redshift distributions of photometric samples is one of the great challenges in photometric surveys like LSST, Euclid, HSC, KiDS, and DES. We present an inference methodology that combines the redshift information from the galaxy photometry with constraints from two-point functions, utilizing cross-correlations with spatially overlapping spectroscopic samples, and illustrate the approach on CosmoDC2 simulations. Our likelihood framework is designed to integrate directly into a typical large-scale structure and weak lensing analysis based on two-point functions. We discuss efficient and accurate inference techniques that allow us to scale the method to the large samples of galaxies to be expected in LSST. We consider statistical challenges like the parametrization of redshift systematics, discuss and evaluate techniques to regularize the sample redshift distributions, and investigate techniques that can help to detect and calibrate sources of systematic error using posterior predictive checks. We evaluate and forecast photometric redshift performance using data from the CosmoDC2 simulations, within which we mimic a DESI-like spectroscopic calibration sample for cross-correlations. Using a combination of spatial cross-correlations and photometry, we show that we can provide calibration of the mean of the sample redshift distribution to an accuracy of at least 0.002(1 + z), consistent with the LSST-Y1 science requirements for weak lensing and large-scale structure probes.

(cosmology:) large-scale structure of Universe↗

Single nucleon emission in relativistic nucleus-nucleus reactions

Significant discrepancies between theory and experiment have previously been noted for nucleon emission via electromagnetic processes in relativistic nucleus-nucleus collisions. The present work investigates the hypothesis that these discrepancies have arisen due to uncertainties about how to deduce the experimental electromagnetic cross section from the total measured cross section. An optical-model calculation of single neutron removal is added to electromagnetic cross sections and compared to the total experimental cross sections. Good agreement is found thereby resolving some of the earlier noted discrepancies. A detailed comparison to the recent work of Benesh, Cook, and Vary is made for both the impact parameter and the nuclear cross section. Good agreement is obtained giving an independent confirmation of the parametrized formulas developed by those authors.

NASA Discipline Radiation Health↗

A Generic Nonlinear Aerodynamic Model for Aircraft

A generic model of the aerodynamic coefficients was developed using wind tunnel databases for eight different aircraft and multivariate orthogonal functions. For each database and each coefficient, models were determined using polynomials expanded about the state and control variables, and an othgonalization procedure. A predicted squared-error criterion was used to automatically select the model terms. Modeling terms picked in at least half of the analyses, which totalled 45 terms, were retained to form the generic nonlinear aerodynamic (GNA) model. Least squares was then used to estimate the model parameters and associated uncertainty that best fit the GNA model to each database. Nonlinear flight simulations were used to demonstrate that the GNA model produces accurate trim solutions, local behavior (modal frequencies and damping ratios), and global dynamic behavior (91% accurate state histories and 80% accurate aerodynamic coefficient histories) under large-amplitude excitation. This compact aerodynamics model can be used to decrease on-board memory storage requirements, quickly change conceptual aircraft models, provide smooth analytical functions for control and optimization applications, and facilitate real-time parametric system identification.

Grauer, Jared A.↗

Mind the gap: the discrepancy between simulation and reality drives interpretations of the Galactic Center Excess

The Galactic Center Excess (GCE) in GeV gamma rays has been debated for over a decade, with the possibility that it might be due to dark matter annihilation or undetected point sources such as millisecond pulsars (MSPs). This study investigates how the gamma-ray emission model (γEM) used in Galactic center analyses affects the interpretation of the GCE's nature. To address this issue, we construct an ultra-fast and powerful inference pipeline based on convolutional Deep Ensemble Networks. We explore the two main competing hypotheses for the GCE using a set of γEMs with increasing parametric freedom. We calculate the fractional contribution (f src ) of a dim population of MSPs to the total luminosity of the GCE and analyze its dependence on the complexity of the γEM. For the simplest γEM, we obtain fsrc = 0.10 ± 0.07, while the most complex model yields f src = 0.79 ± 0.24. In conclusion, we find that the statement about the nature of the GCE (dark matter or not) strongly depends on the assumed γEM. The quoted results for f src do not account for the additional uncertainty arising from the fact that the observed gamma-ray sky is out-of-distribution concerning the investigated γEM iterations. We quantify the reality gap between our γEMs using deep-learning-based One-Class Deep Support Vector Data Description networks, revealing that all employed γEMs have gaps to reality. Finally, our study casts doubt on the validity of previous conclusions regarding the GCE and dark matter, and underscores the urgent need to account for the reality gap and consider previously overlooked "out of domain" uncertainties in future interpretations.

79 ASTRONOMY AND ASTROPHYSICS↗

Dual adaptive control: Design principles and applications

The design of an actively adaptive dual controller based on an approximation of the stochastic dynamic programming equation for a multi-step horizon is presented. A dual controller that can enhance identification of the system while controlling it at the same time is derived for multi-dimensional problems. This dual controller uses sensitivity functions of the expected future cost with respect to the parameter uncertainties. A passively adaptive cautious controller and the actively adaptive dual controller are examined. In many instances, the cautious controller is seen to turn off while the latter avoids the turn-off of the control and the slow convergence of the parameter estimates, characteristic of the cautious controller. The algorithms have been applied to a multi-variable static model which represents a simplified linear version of the relationship between the vibration output and the higher harmonic control input for a helicopter. Monte Carlo comparisons based on parametric and nonparametric statistical analysis indicate the superiority of the dual controller over the baseline controller.

Mookerjee, Purusottam↗