An adaptive hybridizable discontinuous Galerkin approach for cardiac electrophysiology

Multivariable Technique for the Evaluation of the Trans-stenotic Pressure Gradient

PURPOSE

This study establishes a reliable image-based multivariable technique for measuring the trans-stenotic pressure gradient.

METHODS

A self-made in vitro steady flow model based on adjustable velocities and stenotic properties were used as the experimental subject. The pre-stenotic flow velocity, severity, and length of the stenosis were used as the input variables. Based on equations used to fit the plots of the physically measured pressure gradient values versus each input variable, a multivariable formula for the pressure gradient measurement could then be derived. The flow model was scanned using velocity-encoded phase-contrast magnetic resonance imaging (PC-MRI) to validate the derived formula while simultaneously measuring the trans-stenotic pressure gradient. The correlation between the physically-measured pressure gradient values and the pressure gradient values calculated using the new formula were subsequently analyzed.

RESULTS

The results of linear regression analysis using the physically measured pressure gradient values for the new method were compared to values obtained using the simplified Bernoulli equation (R², 0.991, and 0.975, respectively). In a paired t-test, no statistically significant difference was found between the new method and the physical measurements.

CONCLUSIONS

The derived multivariable technique was found to reliably measure the trans-stenotic pressure gradient, with better performance than a traditional procedure based on the simplified Bernoulli equation.

Efficient time splitting schemes for the monodomain equation in cardiac electrophysiology

Approximating the fast dynamics of depolarization waves in the human heart described by the monodomain model is numerically challenging. Splitting methods for the PDE-ODE coupling enable the computation with very fine space and time discretizations. Here, we compare different splitting approaches regarding convergence, accuracy, and efficiency. Simulations were performed for a benchmark problem with the Beeler-Reuter cell model on a truncated ellipsoid approximating the left ventricle including a localized stimulation. For this configuration, we provide a reference solution for the transmembrane potential. We found a semi-implicit approach with state variable interpolation to be the most efficient scheme. The results are transferred to a more physiological setup using a bi-ventricular domain with a complex external stimulation pattern to evaluate the accuracy of the activation time for different resolutions in space and time.

Balancing conduction velocity error in cardiac electrophysiology using a modified quadrature approach

Conduction velocity error is often the main culprit behind the need for very fine spatial discretizations and high computational effort in cardiac electrophysiology problems. In light of this, a novel approach for simulating an accurate conduction velocity in coarse meshes with linear elements is suggested based on a modified quadrature approach. In this approach, the quadrature points are placed at arbitrary offsets of the isoparametric coordinates. A numerical study illustrates the dependence of the conduction velocity on the spatial discretization and the conductivity when using different quadrature rules and calculation approaches. Additionally, examples using the modified quadrature in coarse meshes for wave propagation demonstrate the improved accuracy of the conduction velocity with this method. This novel approach possesses great potential in reducing the computational effort required but remains limited to specific linear elements and experiences a reduction in accuracy for irregular meshes and heterogeneous conductivities. Further research can focus on developing an adaptive quadrature and extending the approach to other element formulations in order to make the approach more generally applicable.

A numerical study on the effects of spatial and temporal discretization in cardiac electrophysiology

Millions of degrees of freedom are often required to accurately represent the electrophysiology of the myocardium due to the presence of discretization effects. This study seeks to explore the influence of temporal and spatial discretization on the simulation of cardiac electrophysiology in conjunction with changes in modeling choices. Several finite element analyses are performed to examine how discretization affects solution time, conduction velocity and electrical excitation. Discretization effects are considered along with changes in the electrophysiology model and solution approach. Two action potential models are considered: the Aliev-Panfilov model and the ten Tusscher-Noble-Noble-Panfilov model. The solution approaches consist of two time integration schemes and different treatments for solving the local system of ordinary differential equations. The efficiency and stability of the calculation approaches are demonstrated to be dependent on the action potential model. The dependency of the conduction velocity on the element size and time step is shown to be different for changes in material parameters. Finally, the discrepancies between the wave propagation in coarse and fine meshes are analyzed based on the temporal evolution of the transmembrane potential at a node and its neighboring Gauss points. Insight obtained from this study can be used to suggest new methods to improve the efficiency of simulations in cardiac electrophysiology.

Stabilized hybrid discontinuous Galerkin finite element method for the cardiac monodomain equation

Numerical methods for solving the cardiac electrophysiology model, which describes the electrical activity in the heart, are proposed. The model problem consists of a nonlinear reaction-diffusion partial differential equation coupled to systems of ordinary differential equations that describes electrochemical reactions in cardiac cells. The proposed methods combine an operator splitting technique for the reaction-diffusion equation with primal hybrid methods for spatial discretization considering continuous or discontinuous approximations for the Lagrange multiplier. A static condensation is adopted to form a reduced global system in terms of the multiplier only. Convergence studies exhibit optimal rates of convergence and numerical experiments show that the proposed schemes can be more efficient than standard numerical techniques commonly used in this context when preconditioned iterative methods are used for the solution of linear systems.

Automatic mapping of atrial fiber orientations for patient-specific modeling of cardiac electromechanics using image registration

Knowledge of appropriate local fiber architecture is necessary to simulate patient-specific electromechanics in the human heart. However, it is not yet possible to reliably measure in vivo fiber directions especially in human atria. Thus, we present a method that defines the fiber architecture in arbitrarily shaped atria using image registration and reorientation methods based on atlas atria with fibers predefined from detailed histological observations. Thereby, it is possible to generate detailed fiber families in every new patient-specific geometry in an automated, time-efficient process. We demonstrate the good performance of the image registration and fiber definition on 10 differently shaped human atria. Additionally, we show that characteristics of the electrophysiological activation pattern that appear in the atlas atria also appear in the patients' atria. We arrive to analogous conclusions for coupled electro-mechano-hemodynamical computations.