Studying the Influence of Finite Element Mesh Size on the Accuracy of Ventricular Tachycardia Simulation
Boyang Cao1,2, Nan Zhang3, Zhenyin Fu1
1College of Biomedical Engineering & Instrument Science, Zhejiang University, 310058 Hangzhou, Zhejiang, China.
Background:
Ventricular tachycardia (VT) is a life-threatening heart condition commonly seen in patients with myocardial infarction (MI). Although personalized computational modeling has been used to understand VT and its treatment noninvasively, this approach can be computationally intensive and time consuming. Therefore, finding a balance between mesh size and computational efficiency is important. This study aimed to find an optimal mesh resolution that minimizes the need for computational resources while maintaining numerical accuracy and to investigate the effect of mesh resolution variation on the simulation results.
Methods:
We constructed ventricular models from contrast-enhanced magnetic resonance imaging data from six patients with MI. We created seven different models for each patient, with average edge lengths ranging from 315 to 645 µm using commercial software, Mimics. Programmed electrical stimulation was used to assess VT inducibility from 19 sites in each heart model.
Results:
The simulation results in the slab model with adaptive tetrahedral mesh (same as in the patient-specific model) showed that the absolute and relative differences in conduction velocity (CV) were 6.1 cm/s and 7.8% between average mesh sizes of 142 and 600 µm, respectively. However, the simulation results in the six patient-specific models showed that average mesh sizes with 350 µm yielded over 85% accuracy for clinically relevant VT. Although average mesh sizes of 417 and 478 µm could also achieve approximately 80% accuracy for clinically relevant VT, the percentage of incorrectly predicted VTs increases. When conductivity was modified to match the CV in the model with the finest mesh size, the overall ratio of positively predicted VT increased.
Conclusions:
The proposed personalized heart model could achieve an optimal balance between simulation time and VT prediction accuracy when discretized with adaptive tetrahedral meshes with an average edge length about 350 µm.
