P. F. Antonietti,L. Dedè,G. Loli,M. Montardini,G. Sangalli,P. Tesini
Abstract:This work proposes a stable and efficient space-time method for the monodomain equation coupled with the Rogers--McCulloch ionic model, which is widely used to simulate electrophysiological wave propagation in the cardiac tissue. By extending the Spline Upwind method and exploiting low-rank matrix approximations, as well as preconditioned solvers, we achieve both significant computational efficiency and accuracy. In particular, we develop a formulation that is both simple and highly effective, designed to minimize spurious oscillations and ensuring computational efficiency. We rigorously validate the method's performance through a series of numerical experiments, showing its robustness and reliability in diverse <a class="link-external link-http" href="http://scenarios.hod" rel="external noopener nofollow">this http URL</a> and also reasonable from the computational cost point of view. For these reasons we validate the method's capability with numerical experiments, focusing on accuracy and computational aspects.
What problem does this paper attempt to address?
The main problem that this paper attempts to solve is to develop a stable and efficient spatio - temporal isogeometric analysis method for simulating the single - domain equation - coupled Rogers - McCulloch ion model in cardiac electrophysiology. Specifically, the paper focuses on the following aspects:
1. **Numerical stability**: When simulating cardiac electrical activity, especially when dealing with sharp wave fronts and propagating pulses, traditional numerical methods may produce non - physical oscillations (i.e., spurious oscillations). These oscillations can seriously affect the accuracy of the results. To this end, the paper introduces the Spline Upwind (SU) method to stabilize the numerical solution and reduce or eliminate these spurious oscillations.
2. **Computational efficiency**: In order to improve computational efficiency, the paper utilizes low - rank matrix approximation techniques and pre - conditioned solvers. These techniques can significantly reduce the time and resources required to solve large - scale linear systems, making it possible to perform efficient cardiac electrophysiological simulations in three - dimensional spatial domains.
3. **Accuracy and reliability**: Through a series of numerical experiments, the paper verifies the performance of the proposed method and demonstrates its robustness and reliability in different scenarios. In particular, for both smooth solutions and cases with sharp layers, the method shows good convergence and stability.
In summary, this paper aims to provide a stable and efficient numerical tool for the research of cardiac electrophysiology, in order to more accurately simulate cardiac electrical activity and support diagnostic and predictive medicine.
### Key formulas
- The governing differential equations of the single - domain equation and the Rogers - McCulloch ion model:
\[
\begin{cases}
C_m\frac{\partial u}{\partial t}-D\Delta u + c_1u(u - a)(u - 1)+c_2uw = f & \text{in }\Omega\times(0, T),\\
\frac{\partial w}{\partial t}-b(u - d_ew)=0 & \text{in }\Omega\times(0, T),\\
\frac{\partial u}{\partial n}=0 & \text{on }\partial\Omega\times[0, T],\\
u = 0 & \text{in }\Omega\times\{0\},\\
w = 0 & \text{in }\Omega\times\{0\},
\end{cases}
\]
where \(n\) is the outer normal vector; \(a, b, c_1, c_2, d_e\) are dimensionless positive parameters; \(C_m\) and \(D\) are the local membrane capacitance and conductivity respectively.
- The definition of the spline upwind stabilization term \(S(u_h; v_h)\):
\[
S(u_h; v_h):=C_m\sum_{k = 1}^{p_t}\int_0^T\tau_k(t)\int_\Omega\theta(x, t)\frac{\partial^k u_h}{\partial t^k}\frac{\partial^k v_h}{\partial t^k}\,d\Omega\,dt
\]
- The definition of the stabilization parameter \(\theta(x, t)\):
\[
\Theta_{i_s,i_t}:=\min\left(\frac{\|C_m\frac{\partial u_h}{\partial t}-D\Delta u_h + c_1u_h(u_h - a)(u_h - 1)+c_2u_hw_h - f\|_{L^\infty(eI_{i_s,i_t})}}{C_m(T^{-1}\|u_h\|_{L^\infty(