Fresca Stefania, Manzoni Andrea, Dedè Luca, Quarteroni Alfio
MOX, Dipartimento di Matematica, Politecnico di Milano, Milan, Italy.
Mathematics Institute, École Polytechnique Fédérale de Lausanne, Lausanne, Switzerland.
Front Physiol. 2021 Sep 22;12:679076. doi: 10.3389/fphys.2021.679076. eCollection 2021.
The numerical simulation of multiple scenarios easily becomes computationally prohibitive for cardiac electrophysiology (EP) problems if relying on usual high-fidelity, full order models (FOMs). Likewise, the use of traditional reduced order models (ROMs) for parametrized PDEs to speed up the solution of the aforementioned problems can be problematic. This is primarily due to the strong variability characterizing the solution set and to the nonlinear nature of the input-output maps that we intend to reconstruct numerically. To enhance ROM efficiency, we proposed a new generation of non-intrusive, nonlinear ROMs, based on deep learning (DL) algorithms, such as convolutional, feedforward, and autoencoder neural networks. In the proposed DL-ROM, both the nonlinear solution manifold and the nonlinear reduced dynamics used to model the system evolution on that manifold can be learnt in a non-intrusive way thanks to DL algorithms trained on a set of FOM snapshots. DL-ROMs were shown to be able to accurately capture complex front propagation processes, both in physiological and pathological cardiac EP, very rapidly once neural networks were trained, however, at the expense of huge training costs. In this study, we show that performing a prior dimensionality reduction on FOM snapshots through randomized proper orthogonal decomposition (POD) enables to speed up training times and to decrease networks complexity. Accuracy and efficiency of this strategy, which we refer to as POD-DL-ROM, are assessed in the context of cardiac EP on an idealized left atrium (LA) geometry and considering snapshots arising from a NURBS (non-uniform rational B-splines)-based isogeometric analysis (IGA) discretization. Once the ROMs have been trained, POD-DL-ROMs can efficiently solve both physiological and pathological cardiac EP problems, for any new scenario, in real-time, even in extremely challenging contexts such as those featuring circuit re-entries, that are among the factors triggering cardiac arrhythmias.
对于心脏电生理学(EP)问题,如果依赖通常的高保真全阶模型(FOM),那么多个场景的数值模拟很容易在计算上变得难以承受。同样,使用传统的降阶模型(ROM)来求解参数化偏微分方程(PDE)以加速上述问题的求解可能会出现问题。这主要是由于解集具有很强的变异性以及我们打算通过数值方法重建的输入 - 输出映射的非线性性质。为了提高ROM的效率,我们基于深度学习(DL)算法,如卷积神经网络、前馈神经网络和自动编码器神经网络,提出了新一代非侵入式非线性ROM。在所提出的DL-ROM中,由于在一组FOM快照上训练的DL算法,非线性解流形和用于在该流形上对系统演化进行建模的非线性降阶动力学都可以以非侵入的方式学习。研究表明,一旦神经网络训练完成,DL-ROM能够非常快速地准确捕捉生理和病理心脏EP中复杂的波前传播过程,然而,代价是巨大的训练成本。在本研究中,我们表明通过随机化本征正交分解(POD)对FOM快照进行先验降维能够加快训练时间并降低网络复杂度。我们将这种策略称为POD-DL-ROM,并在理想化左心房(LA)几何结构的心脏EP背景下,考虑基于非均匀有理B样条(NURBS)的等几何分析(IGA)离散化产生的快照,评估了该策略的准确性和效率。一旦ROM训练完成,POD-DL-ROM能够实时有效地解决任何新场景下的生理和病理心脏EP问题,即使在极具挑战性的情况下,如存在引发心律失常的因素之一的折返环的情况。