SEARCH · Search NASA

Results for “Total variance”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

104 records · Page 6

All sources NASA NTRS DOE OSTI

DOE OSTI2021

A comparative study of machine learning models for predicting the state of reactive mixing

Mixing phenomena are important mechanisms controlling flow, species transport, and reaction processes in fluids and porous media. Accurate predictions of reactive mixing are critical for many Earth and environmental science problems such as contaminant fate and remediation, macroalgae growth, and plankton biomass evolution. Here, to investigate the evolution of mixing dynamics under different scenarios (e.g., anisotropy, fluctuating velocity fields), a finite-element-based numerical model was built to solve the fast, irreversible bimolecular reaction-diffusion equations to simulate a range of reactive-mixing scenarios. A total of 2,315 simulations were performed using different sets of model input parameters comprising various spatial scales of vortex structures in the velocity field, time-scales associated with velocity oscillations, the perturbation parameter for the vortex-based velocity, anisotropic dispersion contrast (i.e., ratio of longitudinal-to-transverse dispersion), and molecular diffusion. The outputs comprised concentration profiles of reactants and products. The inputs to and outputs from these simulations were concatenated into feature and label matrices, respectively, to train 20 different machine learning (ML) models intended to emulate system behavior. These 20 ML emulators, based on linear methods, Bayesian methods, ensemble learning methods, and multilayer perceptrons (MLPs), were trained to classify the state of mixing and predict three quantities of interest (QoIs) characterizing species production, decay (i.e., average concentration, square of average concentration), and degree of mixing (i.e., variances of species concentration). Unsurprisingly, linear classifiers and regressors failed to reproduce the QoIs; however, ensemble methods (classifiers and regressors) and the MLP model accurately classified the state of reactive mixing and the QoIs. Among ensemble methods, random forest and decision-tree-based AdaBoost faithfully predicted the QoIs. At run time, trained ML emulators produced results times faster than the finite-element simulations. Due to their low computational expense and high accuracy, ensemble and MLP models are excellent emulators for these numerical simulations and great utilities in uncertainty quantification exercises, which can require 1,000s of forward model runs.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

Tensor-to-scalar ratio forecasts for extended LiteBIRD frequency configurations

LiteBIRD is a planned JAXA-led cosmic microwave background (CMB) B-mode satellite experiment aiming for launch in the late 2020s, with a primary goal of detecting the imprint of primordial inflationary gravitational waves. Its current baseline focal-plane configuration includes 15 frequency bands between 40 and 402 GHz, fulfilling the mission requirements to detect the amplitude of gravitational waves with the total uncertainty on the tensor-to-scalar ratio, δr, down to δr < 0.001. A key aspect of this performance is accurate astrophysical component separation, and the ability to remove polarized thermal dust emission is particularly important. In this paper we note that the CMB frequency spectrum falls off nearly exponentially above 300 GHz relative to the thermal dust spectral energy distribution, and a relatively minor high frequency extension can therefore result in even lower uncertainties and better model reconstructions. Specifically, we compared the baseline design with five extended configurations, while varying the underlying dust modeling, in each of which the High-Frequency Telescope (HFT) frequency range was shifted logarithmically toward higher frequencies, with an upper cutoff ranging between 400 and 600 GHz. In each case, we measured the tensor-to-scalar ratio r uncertainty and bias using both parametric and minimum-variance component-separation algorithms. When the thermal dust sky model includes a spatially varying spectral index and temperature, we find that the statistical uncertainty on r after foreground cleaning may be reduced by as much as 30–50% by extending the upper limit of the frequency range from 400 to 600 GHz, with most of the improvement already gained at 500 GHz. We also note that a broader frequency range leads to higher residuals when fitting an incorrect dust model, but also it is easier to discriminate between models through higher χ 2 sensitivity. Even in the case in which the fitting procedure does not correspond to the underlying dust model in the sky, and when the highest frequency data cannot be modeled with sufficient fidelity and must be excluded from the analysis, the uncertainty on r increases by only about 5% for a 500 GHz configuration compared to the baseline.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2020

Multi-Model and Multi-Scale Global Sensitivity Analysis for Identifying Controlling Processes of Complex Systems

An environmental model consists of multiple process level sub-models, and each sub-model represents a process that is key to the operation of the simulated system. Global sensitivity analysis methods have been widely used to identify important processes for system model development and improvement. The existing methods of global sensitivity analysis only consider parametric uncertainty, and are not capable of handling model uncertainty caused by multiple process models that arise from competing hypotheses about one or more processes. To address this problem, this project develops a new method to probe model output sensitivity to competing process models by integrating model averaging methods with variance-based global sensitivity analysis to address uncertainty in process models and parameters. The new method yields three process sensitivity indices. The first one is called first-order process sensitivity index, and it is derived as a single summary measure of relative process importance. Evaluating the index is computationally expensive, because it relies in a Monte Carlo scheme that requires thousands and even millions of model executions. To reduce computational cost, this project develops a computationally efficient, quasi Monte Carlo method, and this method is presented in Chapter 2 of this report with and a numerical example for demonstration. The numerical example shows that the results of the quasi Monte Carlo method are substantially close to those of the full Monte Carlo method, but the computational cost of the quasi Monte Carlo method is only 0.7% of that of the full Monte Carlo method. The second index is called total-effect process sensitivity index, and it measures interactions between different processes. Therefore, this sensitivity index includes the first-order process sensitivity index, and can be used to identify influential processes. On the other hand, the total-effect process sensitivity index can also be used to screen non-influential processes. This is demonstrated by two numerical examples using the Sobol-G* functions and groundwater flow models that consider recharge process, geological process, and snowmelt process. The numerical examples shows that the total-effect process sensitivity index is more informative than the first-order process sensitivity. The derivation of the process sensitivity index and the numerical examples are discussed in Chapter 3. Chapter 4 presents two computationally efficient methods for screening non-influential processes to exclude them from further investigation. The two methods are the multi-model difference-based sensitivity (MMDS) analysis method, which can be implemented using the Latin Hypercube Sampling. The second one is the implementation of MMDS method using a binning method. The numerical example for the Sobol-G* function indicates the two methods are capable of identifying non-influential models, and the numerical examples for the groundwater flow and reactive transport show that the two methods are effective for groundwater problems. However, it should be noted that the two methods are numerical approximations, and they can only be used for screening non-influential processes, not for ranking importance of system processes. All the sensitivity analysis methods are implemented by developing python codes, and the codes are in a software called SAMMPY: a python package for process sensitivity analysis under multiple models. The SAMMPY design and structure are discussed in Chapter 5, and the package is released to the public for free download.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Exploring the Influence of Linker Substitution and Ratios on Cooperative Framework Flexibility through the Mixed-Linker Approach

Cooperatively flexible metal–organic frameworks that exhibit step-shaped, or Type V-like, adsorption–desorption profiles can lower requisite pressure–temperature swings and thus energy input, necessary for an array of gas storage, delivery, and separations applications. However, such benefits are lost if the pressure threshold of the adsorption and desorption steps at a given temperature does not match the conditions dictated by the application, such as H 2 storage and delivery or olefin–paraffin separations. Unfortunately, the discovery of cooperatively flexible frameworks remains wholly serendipitous and de novo design remains impossible. Accordingly, there is a great need to further our understanding of flexibility such that we can intuitively derivatize known frameworks and, ultimately, design entirely new ones to meet the requisite conditions of energy-consumptive processes. In this work, we demonstrate that the mixed-linker, or multivariate, approach is a powerful tool for the derivation of a known flexible framework, with variances in linker substitution and ratio giving rise to a family wherein significant changes to the step-shaped adsorption–desorption profiles for multiple adsorbates are observed. Specifically, we report 12 isostructural mixed-linker derivatives of CdIF-13 (sod–Cd(benzimidazolate) 2 ) with six point-modified benzimidazole linkers spontaneously synthesized through prototypical solvothermal conditions. Furthermore, each is shown by PXRD to exhibit similar reversible flexibility to CdIF-13 and by TGA and DSC to be similarly thermally stable. Isothermal gas adsorption measurements with N 2 at 77 K, CO 2 at 195 K, and propane at 298 K demonstrate the wide-ranging, adsorbate-dependent effects of linker substitution and ratio, including dramatic reductions in the adsorption threshold pressure, evolution of multiple steps, and complete absence of observable adsorption. Aided by prior crystallographic characterization of CdIF-13, the observed trends were analyzed in an attempt to resolve the influence of linker substitution and ratios on structural behavior. In total, this work illustrates how the mixed-linker approach enables the synthesis and study of a large catalog of functionally modified frameworks, which can help identify the variables that influence flexibility and sorption behavior.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Evaluating Enhanced Reservoir Evaporation Losses From CMIP6‐Based Future Projections in the Contiguous United States

Abstract Enhanced reservoir evaporation has become an emerging concern regarding water loss, especially when compounded with the ever‐increasing water demand. In this study, we evaluated the evaporation rates and losses for 678 major reservoirs (representing nearly 90% of total storage capacity) in the Contiguous United States over historical baseline (1980–2019), near‐term (2020–2039), and mid‐term (2040–2059) future periods. The evaporation rate was estimated using the Lake Evaporation Model (LEM), an advanced lake evaporation model that addresses both heat storage and fetch effects, driven by multi‐ensemble downscaled Coupled Model Intercomparison Projects 6 (CMIP6) projections under the SSP585 emission scenario. The results project that the evaporation loss may increase by 2.5 × 10 7 m 3 /yr through the research period (1980–2059). Among all regions, the Rio Grande is projected to have the largest increasing rate in the near‐term and mid‐term future, with values of 7.11% of 10.25%, respectively. At the seasonal scale, the most significant increase in the evaporation rate is projected during the fall. The evaporation is projected to increase faster than the streamflow over many of the regions in the southwestern US during the summer/fall, suggesting that the shortage of water will be further exacerbated. The climate models contribute the most to the variance, as compared to the other components related to the projection of evaporation losses (e.g., hydrological model, downscaling method, and historical meteorological data set). These findings demonstrate the need to consider accelerated water loss through open water evaporation in long‐term water resources planning across various spatiotemporal scales.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Spatiotemporal Variability in Wind Turbine Blade Leading Edge Erosion

Wind turbine blade leading edge erosion (LEE) reduces energy production and increases wind energy operation and maintenance costs. Degradation of the blade coating and ultimately damage to the underlying blade structure are caused by collisions of falling hydrometeors with rotating blades. The selection of optimal methods to mitigate/reduce LEE are critically dependent on the rates of coating fatigue accumulation at a given location and the time variance in the accumulation of material stresses. However, no such assessment currently exists for the United States of America (USA). To address this research gap, blade coating lifetimes at 883 sites across the USA are generated based on high-frequency (5-min) estimates of material fatigue derived using a mechanistic model and robust meteorological measurements. Results indicate blade coating failure at some sites in as few as 4 years, and that the frequency and intensity of material stresses are both highly episodic and spatially varying. Time series analyses indicate that up to one-third of blade coating lifetime is exhausted in just 360 5-min periods in the Southern Great Plains (SGP). Conversely, sites in the Pacific Northwest (PNW) exhibit the same level of coating lifetime depletion in over three times as many time periods. Thus, it may be more cost-effective to use wind turbine deregulation (erosion-safe mode) for damage reduction and blade lifetime extension in the SGP, while the application of blade leading edge protective measures may be more appropriate in the PNW. Annual total precipitation and mean wind speed are shown to be poor predictors of blade coating lifetime, re-emphasizing the need for detailed modeling studies such as that presented herein.

Pryor, Sara C. (ORCID:0000000348473440)↗

DOE OSTI2023

The comparison between effects of Taichi and conventional exercise on functional mobility and balance in healthy older adults: a systematic literature review and meta-analysis

Background Taichi is beneficial for functional mobility and balance in older adults. However, such benefits of Taichi when comparing to conventional exercise (CE) are not well understood due to large variance in study protocols and observations. Methods We reviewed publications in five databases. Eligible studies that examined the effects of Taichi on the outcomes of functional mobility and balance in healthy older adults as compared to CE were included. Subgroup analyses compared the effects of different types of CE (e.g., single and multiple-type exercise) and different intervention designs (e.g., Taichi types) on those outcomes (Registration number: CRD42022331956). Results Twelve studies consisting of 2,901 participants were included. Generally, compared to CE, Taichi induced greater improvements in the performance of Timed-Up-and-Go (SMD = −0.18, [−0.33 to −0.03], p = 0.040, I 2 = 59.57%), 50-foot walking (MD = −1.84 s, [−2.62 to −1.07], p < 0.001, I 2 = 0%), one-leg stance with eyes open (MD = 6.00s, [2.97 to 9.02], p < 0.001, I 2 = 83.19%), one-leg stance with eyes closed (MD = 1.65 s, [1.35 to 1.96], p < 0.001, I 2 = 36.2%), and functional reach (SMD = 0.7, [0.32 to 1.08], p < 0.001, I 2 = 86.79%) tests. Subgroup analyses revealed that Taichi with relatively short duration (<20 weeks), low total time (≤24 h), and/or using Yang-style, can induce significantly greater benefits for functional mobility and balance as compared to CE. Uniquely, Taichi only induced significantly greater improvements in Timed-Up-and-Go compared to single- (SMD = −0.40, [−0.55 to −0.24], p < 0.001, I 2 = 6.14%), but not multiple-type exercise. A significant difference between the effects of Taichi was observed on the performance of one-leg stance with eyes open when compared to CE without balance (MD = 3.63 s, [1.02 to 6.24], p = 0.006, I 2 = 74.93%) and CE with balance (MD = 13.90s, [10.32 to 17.48], p < 0.001, I 2 = 6.1%). No other significant difference was shown between the influences of different CE types on the observations. Conclusion Taichi can induce greater improvement in functional mobility and balance in older adults compared to CE in a more efficient fashion, especially compared to single-type CE. Future studies with more rigorous design are needed to confirm the observations here.

Li, Yiting↗

DOE OSTI2025

Constraints on Local Primordial Non-Gaussianity with 3D Velocity Reconstruction from the Kinetic Sunyaev-Zeldovich Effect

The cosmic velocity field is an unbiased probe of the total matter distribution but is challenging to measure directly at intermediate and high redshifts. The large-scale velocity field imprints a signal in the cosmic microwave background (CMB) through the kinetic Sunyaev-Zeldovich (kSZ) effect. We perform the first 3D reconstruction of the large-scale velocity field from the kSZ effect by applying a quadratic estimator to CMB temperature maps and the 3D positions of galaxies. We do so by combining CMB data from the fifth data release of the Atacama Cosmology Telescope (in combination with Planck) and a spectroscopic galaxy sample from the Sloan Digital Sky Survey. We then measure the galaxy-velocity cross-power spectrum and detect the presence of the kSZ signal at a signal-to-noise ratio of 7.2⁢𝜎. Using this galaxy-velocity cross-correlation alone, we constrain the amplitude of local primordial non-Gaussianity finding 𝑓 NL =−9⁢0$^{+210}_{−350}$. In conclusion, this pathfinder measurement sets the stage for joint galaxy-CMB kSZ constraints to significantly enhance the 𝑓 NL information obtained from galaxy surveys through sample variance cancellation.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

Evaluating Sea Breezes and Associated Convective Cloud Evolution in the Model Gray Zone

We characterize convective clouds associated with sea‐breeze circulations (SBC) using multi‐agency observations and multi‐case ensemble model simulations. The focus is on assessing convective cloud lifecycle properties and their merging behavior, as well as the environmental conditions they are embedded in, particularly SBC features. In total, 46 SBC days over the Houston‐Galveston region are selected and simulated using the Weather Research and Forecasting (WRF) model at a gray zone scale with a forecast‐like parameterization setup. Advanced techniques, including change‐point detection, a Lagrangian cloud tracking method, and a newly developed cell merging and splitting detection algorithm, are applied and/or developed for this study. Our findings indicate that the WRF model at 1 km grid spacing well represents the thermodynamic conditions over the region, as well as SBC timing and intensity. However, for the associated convective cells, WRF overestimates the 30‐dBZ echo top height, cell area, and maximum radar reflectivity compared to radar observations. This overestimation is potentially due to under‐resolved entrainment processes, an overestimated merging frequency, and the overestimation of updraft intensity. Furthermore, the model exhibits a deficiency in simulating congestus clouds, showing a more rapid transition from shallow to deep convection compared to observed behavior. Moreover, observations indicate stronger, deeper, and wider clouds when merging happens. Conversely, in simulations, the merging process does not necessarily lead to higher or longer‐lived cells, as many cases experience rapid and frequent merging and splitting which may result in more variance in convective updraft velocity during the convection lifetime.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Mixing Pathways in Simple Box Models

Abstract Tracer variance budgets can be used to estimate bulk mixing in a control volume. For example, simple, analytical, bulk formulations of salt mixing, defined here as the destruction of salinity variance, can be found for estuaries with a riverine source of freshwater and a two-layer exchange flow at the mouth using salinity as a representative tracer. For a steady case, the bulk salt mixing M can be calculated as , where S in and S out are the representative salinities in the estuarine exchange flow, and Q R and Q in are the river and landward volume fluxes, respectively. The bulk salt mixing M can be considered as the sum of mixing pathways, where each pathway has a mixing of Q (Δ S ) 2 , where Q is the volume transport and Δ S is the salinity difference across the pathway. For the estuary case, one mixing path is associated with the river inflow, and the other is associated with the inflow of salty, oceanic water. This concept of linking mixing to input–output pathways is extended, in simple box models, from estuaries to scenarios with multiple inputs/outputs, as might be found in a complex estuarine/fjord network, in a region on a continental shelf, or any other control volume with multiple exchanges. This approach allows for the estimation of the relative contributions of each input–output pathway to the total mixing within a control volume.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2026

Computationally inexpensive part-scale thermal history of additive friction-stir deposition

This study presents an analytical model for steady-state power generation and tool heat loss in additive friction-stir deposition (AFSD), developed to enable part-scale thermal simulation while remaining computationally inexpensive. The model predicts total generated power, yielding 3.7–4.7 kW across deposition temperature setpoints of 400–460 °C for the deposition of AA6061 with a Be-Cu tool. This corresponds to 90–95% of the reported spindle power. Tool heat loss is experimentally determined by calibrating a steady-state energy balance between the generated power, the substrate-deposition thermal gradient, and a temperature dependent tool heat loss term: q tool (T) = a + b (T - 400°C) with a = 2.7 x 10 6 Wm -2 and b = 9.5 x 10 3 Wm -2 K -1 . The calibration indicates that about 69% of the generated heat is conducted into the tool for this configuration, which is much higher than previously reported. The calibrated heat-source is implemented in finite element software (Adamantine) to simulate the transient thermal history of a 100 cm 3 representative build in 8 min on a standard desktop (at 0.635 mm build-height resolution). For the first three layers, the substrate temperatures between simulation and experiment are within 10% mean absolute percentage error. Sensitivity analysis indicates that uncertainties in average deposition temperature and deformation localization (stir-zone geometry, depth, and spatial dependance of strain-rate and flow stress) dominate model variance, motivating additional experimental verification.

Additive Friction-Stir Deposition↗

DOE OSTI2020

Nitrogen Fertilization Restructured Spatial Patterns of Soil Organic Carbon and Total Nitrogen in Switchgrass and Gamagrass Croplands in Tennessee USA

Nitrogen (N) fertilizers can potentially alter spatial distribution of soil organic carbon (SOC) and total nitrogen (TN) concentrations in croplands such as switchgrass (SG: Panicum virgatum L.) and gamagrass (GG: Tripsacum dactyloides L.), but it remains unclear whether these effects are the same between crops and under different rates of fertilization. 13 C and 15 N are two important proxy measures of soil biogeochemistry, but they were rarely examined as to their spatial distributions in soil. Based on a three-year long fertilization experiment in Middle Tennessee, USA, the top mineral horizon soils (0–15 cm) were collected using a spatially explicit design within two 15-m 2 plots under three fertilization treatments in SG and GG croplands. A total of 288 samples were collected based on 12 plots and 24 samples in each plot. The fertilization treatments were no N input (NN), low N input (LN: 84 kg N ha -1 in urea) and high N input (HN: 168 kg N ha -1 in urea). The SOC, TN, SOC/TN (C: N), δ 13 C and δ 15 N were quantified and their within-plot variations and spatial distributions were achieved via descriptive and geostatistical methods. Results showed that SG generally displayed 10~120% higher plot-level variations in all variables than GG, and the plot-level variations were 20~77% higher in NN plots than LN and HN plots in SG but they were comparable in unfertilized and fertilized plots in GG. Relative to NN, LN and HN showed more significant surface trends and spatial structures in SOC and TN in both croplands, and the fertilization effect appeared more pronounced in SG. Spatial patterns in C: N, δ 13 C and δ 15 N were comparable among different fertilization treatments in both croplands. The descending within-plot variations were also identified among variables (SOC > TN > δ 15 N > C: N > δ 13 C). This study demonstrated that N fertilizations generally reduced the plot-level variance and simultaneously re-established spatial structures of SOC and TN in bioenergy croplands, which little varied with fertilization rate but was more responsive in switchgrass cropland.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2024

Empirical Validation of UBEM: An Assessment of Bias in Urban Building Energy Modeling for Chicago

Residential and commercial buildings currently account for 30% of total global final energy consumption. Urban-scale building energy modeling (UBEM) can enable scalable investments and unlock building improvements by quantifying energy, demand, emissions, and cost reductions of specific measures or packages for building-specific technologies in large geographic regions. While the sophistication of UBEM data sources and technologies have increased dramatically in the past decade, there remains a knowledge gap for empirical validation and sources of bias between building-specific energy models and measured data at varying geographic scales.As UBEM continues to develop, systemic analysis of accuracy, bias, and limitations of the resulting models is necessary to inform best practices and move toward standardization. These are characterized for the Automatic Building Energy Modeling (AutoBEM) software suite with an initial case study involving metered electricity consumption data from 247,188 buildings in Chicago, Illinois, USA - averaged across years 2019-2021 - compared to the following datasets: (1) the AutoBEM-generated nation-scale Model America version 2 (MAv2) data for 596,064 buildings, (2) tax assessor data for 579,829 buildings, (3) tax assessor data filled with MAv2, and (4) 102 representative dynamic archetypes. The accuracy is reported for every building type and vintage combination, along with multiple sources of bias for unique building descriptors. The AutoBEM simulation workflow produced energy consumption estimates that closely match aggregated metered electricity consumption data for different types of buildings constructed during various time periods at the city scale - with initial normalized mean bias error of 10.9%, and 1.1% after removing outliers. Contribution of statistically significant factors including building type, land use, age, and size to variance in UBEM bias is quantified.

Garg, Ankur↗

DOE OSTI2025

Personalized and uncertainty-aware coronary hemodynamics simulations: From Bayesian estimation to improved multi-fidelity uncertainty quantification

Non-invasive simulations of coronary hemodynamics have improved clinical risk stratification and treatment outcomes for coronary artery disease, compared to relying on anatomical imaging alone. However, simulations typically use empirical approaches to distribute total coronary flow amongst the arteries in the coronary tree, which ignores patient variability, the presence of disease, and other clinical factors. Further, uncertainty in the clinical data often remains unaccounted for in the modeling pipeline. We present an end-to-end uncertainty-aware pipeline to (1) personalize coronary flow simulations by incorporating vessel-specific coronary flows as well as cardiac function; and (2) predict clinical and biomechanical quantities of interest with improved precision, while accounting for uncertainty in the clinical data. We assimilate patient-specific measurements of myocardial blood flow from clinical CT myocardial perfusion imaging to estimate branch-specific coronary artery flows. Simulated noise in the clinical data is used to estimate the joint posterior distributions of the model parameters using adaptive Markov Chain Monte Carlo sampling. Additionally, the posterior predictive distribution for the relevant quantities of interest is determined using a new approach combining multi-fidelity Monte Carlo estimation with non-linear, data-driven dimensionality reduction. This leads to improved correlations between high- and low-fidelity model outputs. Our framework accurately recapitulates clinically measured cardiac function as well as branch-specific coronary flows under measurement noise uncertainty. We observe substantial reductions in confidence intervals for estimated quantities of interest compared to single-fidelity Monte Carlo estimation and state-of-the-art multi-fidelity Monte Carlo methods. This holds especially true for quantities of interest that showed limited correlation between the low- and high-fidelity model predictions. In addition, the proposed multi-fidelity Monte Carlo estimators are significantly cheaper to compute than traditional estimators, under a specified confidence level or variance. The proposed pipeline for personalized and uncertainty-aware predictions of coronary hemodynamics is based on routine clinical measurements and recently developed techniques for CT myocardial perfusion imaging. The proposed pipeline offers significant improvements in precision and reduction in computational cost.

Bayesian parameter estimation↗