Multifidelity Monte Carlo methods rely on a hierarchy of possibly less accurate but statistically correlated simplified or reduced models, in order to accelerate the estimation of statistics of high-fidelity models without compromising the accuracy of the estimates. This approach has recently gained widespread attention in uncertainty quantification [1]. This is partly due to the availability of optimal strategies for the estimation of the expectation of scalar quantities-of-interest [2]. In practice, the optimal strategy for the expectation is also used for the estimation of variance and sensitivity indices [3]. However, a general strategy is still lacking for vector-valued problems, nonlinearly statistically-dependent models, and estimators for which a closed-form expression of the error is unavailable. The focus of the present work is to generalize the standard multifidelity estimators to the above cases. The proposed generalized estimators lead to an optimization problem that can be solved analytically and whose coefficients can be estimated numerically with few runs of the high- and low-fidelity models. We analyze the performance of the proposed approach on a selected number of experiments, with a particular focus on cardiac electrophysiology, where a hierarchy of physics-based low-fidelity models is readily available.

High-dimensional and higher-order multifidelity Monte Carlo estimators / Quaglino, A.; Pezzuto, S.; Krause, R.. - In: JOURNAL OF COMPUTATIONAL PHYSICS. - ISSN 0021-9991. - 388:(2019), pp. 300-315. [10.1016/j.jcp.2019.03.026]

High-dimensional and higher-order multifidelity Monte Carlo estimators

Pezzuto S.;
2019-01-01

Abstract

Multifidelity Monte Carlo methods rely on a hierarchy of possibly less accurate but statistically correlated simplified or reduced models, in order to accelerate the estimation of statistics of high-fidelity models without compromising the accuracy of the estimates. This approach has recently gained widespread attention in uncertainty quantification [1]. This is partly due to the availability of optimal strategies for the estimation of the expectation of scalar quantities-of-interest [2]. In practice, the optimal strategy for the expectation is also used for the estimation of variance and sensitivity indices [3]. However, a general strategy is still lacking for vector-valued problems, nonlinearly statistically-dependent models, and estimators for which a closed-form expression of the error is unavailable. The focus of the present work is to generalize the standard multifidelity estimators to the above cases. The proposed generalized estimators lead to an optimization problem that can be solved analytically and whose coefficients can be estimated numerically with few runs of the high- and low-fidelity models. We analyze the performance of the proposed approach on a selected number of experiments, with a particular focus on cardiac electrophysiology, where a hierarchy of physics-based low-fidelity models is readily available.
2019
Quaglino, A.; Pezzuto, S.; Krause, R.
High-dimensional and higher-order multifidelity Monte Carlo estimators / Quaglino, A.; Pezzuto, S.; Krause, R.. - In: JOURNAL OF COMPUTATIONAL PHYSICS. - ISSN 0021-9991. - 388:(2019), pp. 300-315. [10.1016/j.jcp.2019.03.026]
File in questo prodotto:
Non ci sono file associati a questo prodotto.

I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione

Utilizza questo identificativo per citare o creare un link a questo documento: https://hdl.handle.net/11572/360512
 Attenzione

Attenzione! I dati visualizzati non sono stati sottoposti a validazione da parte dell'ateneo

Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 11
  • ???jsp.display-item.citation.isi??? 10
social impact