Back
Fulltext
Journal article

Space-time multilevel Monte Carlo methods and their application to cardiac electrophysiology

  • Ben Bader, Seif Center for Computational Medicine in Cardiology, Istituto di scienza computazionale (ICS), Facoltà di scienze informatiche, Università della Svizzera italiana, Svizzera
  • Benedusi, Pietro Center for Computational Medicine in Cardiology, Istituto di scienza computazionale (ICS), Facoltà di scienze informatiche, Università della Svizzera italiana, Svizzera
  • Quaglino, Alessio Center for Computational Medicine in Cardiology, Istituto di scienza computazionale (ICS), Facoltà di scienze informatiche, Università della Svizzera italiana, Svizzera - NNAISENSE SA, Lugano, Svizzera
  • Zulian, Patrick Center for Computational Medicine in Cardiology, Istituto di scienza computazionale (ICS), Facoltà di scienze informatiche, Università della Svizzera italiana, Svizzera
  • Krause, Rolf Center for Computational Medicine in Cardiology, Istituto di scienza computazionale (ICS), Facoltà di scienze informatiche, Università della Svizzera italiana, Svizzera
Show more…
    05.02.2021
Published in:
  • Journal of computational physics. - 2021, vol. 433, no. May, p. 110164
Uncertainty quantification
Multilevel methods
Space-time finite elements (3D+1)
Cardiac electrophysiology
Monodomain equation
English We present a novel approach aimed at high-performance uncertainty quantification for time-dependent problems governed by partial differential equations. In particular, we consider input uncertainties described by a Karhunen-Loève expansion and compute statistics of high-dimensional quantities-of-interest, such as the cardiac activation potential. Our methodology relies on a close integration of multilevel Monte Carlo methods, parallel iterative solvers, and a space-time discretization. This combination allows for space-time adaptivity, time-changing domains, and to take advantage of past samples to initialize the space-time solution. The resulting sequence of problems is distributed using a multilevel parallelization strategy, allocating batches of samples having different sizes to a different number of processors. We assess the performance of the proposed framework by showing in detail its application to the solution of nonlinear equations arising from cardiac electrophysiology. Specifically, we study the effect of spatially-correlated perturbations of the heart fibers' conductivities on the mean and variance of the resulting activation map. As shown by the experiments, the theoretical rates of convergence of multilevel Monte Carlo are achieved. Moreover, the total computational work for a prescribed accuracy is reduced by an order of magnitude with respect to standard Monte Carlo methods.
Language
  • English
Classification
Computer science and technology
License
CC BY-NC-ND
Open access status
hybrid
Identifiers
Persistent URL
https://n2t.net/ark:/12658/srd1319324
Statistics
Document views: 83 File downloads:
  • BenBader_2021_Else_jcp.pdf: 143