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-Loeeve 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.
翻译:我们提出了一种新颖的方法,旨在对部分差异方程式所制约的、时间依赖时间的问题进行高性能不确定性量化,特别是,我们考虑了Karhunen-Loeeve扩展和计算高维利益数量统计(如心脏激活潜能)所描述的投入不确定性,我们的方法依靠的是将多层次的蒙特卡洛方法、平行迭接求解答器和时空分解方法紧密结合,这种结合可以允许空间时间适应性、时间变化域,并利用过去的样本来启动时空解决方案。由此产生的问题序列采用多级平行战略进行分配,将不同尺寸的样本分到不同数量的处理器。我们评估拟议框架的绩效,详细展示其对由心电物理学产生的非线性方程式解决方案的应用情况。具体地说,我们研究心脏纤维的与空间变化性扰动对启动图的平均值和差异的影响。实验表明,多层次蒙特卡洛的合并理论速度是按不同处理器数量分配的。此外,我们通过详细显示拟议框架在对心脏电物理学产生的非线性方程式解决方案的应用情况。我们研究了与计算方法的总体精确度要求的计算方法。