GPU-accelerated N-body simulation for arrhythmia research and ablation training
The Interactive Left Atrium Model is a GPU-accelerated N-body computational simulation that enables real-time, interactive exploration of cardiac arrhythmias. Implemented in CUDA C++ with OpenGL visualization, this project combines physics-based cardiac muscle modeling with interactive ablation capabilities for both research and medical education.
Published: Journal of Electrocardiology, Volume 86 (2024), Article 153762 | DOI: 10.1016/j.jelectrocard.2024.153762
CUDA implementation on NVIDIA GPUs (Dual RTX A6000) for real-time interactive simulation of thousands of cardiac myocytes
Validated progression from 1D strands to patient-specific 3D geometries from medical imaging
Reproduces real arrhythmias (micro-reentry, flutter, SVT) and demonstrates ablation effectiveness
Irregular heartbeats affect millions worldwide and are leading contributors to stroke and heart failure. The left atrium is the primary site of complex arrhythmias, particularly atrial fibrillation (AF) and supraventricular tachycardia (SVT).
Catheter ablation is the gold standard intervention, but success rates vary widely (70-80%). Procedural planning relies heavily on operator experience rather than patient-specific computational modeling.
An interactive digital twin enables researchers and clinicians to understand arrhythmia mechanisms, test ablation strategies, and plan procedures without patient risk—all in real time.
Cardiac muscle has a unique property that distinguishes it from other muscle types: a long refractory period. This means that once the heart muscle contracts, it must wait before it can contract again. This built-in delay is essential for proper heart function, ensuring that the heart chambers fill with blood between beats. Our model simulates this critical physiological behavior, allowing researchers to study how disruptions in this timing can lead to arrhythmias.
Our model uses an N-body simulation approach where individual nodes represent cardiac muscle cells (myocytes) connected by virtual muscle fibers. These connections allow electrical signals to propagate through the tissue, mimicking how real heart muscle conducts impulses. The model captures both electrical activation and mechanical contraction, enabling users to observe and manipulate cardiac behavior in real-time. This interactive framework makes it possible to induce arrhythmias and test virtual ablation strategies.
The model was systematically validated through dimensional progression, starting with simple geometries and advancing to anatomically accurate structures:
Validated electrical wavefront propagation and conduction blockage. The fundamental unit for testing signal transmission.
Demonstrated macro-reentrant circuits and reentry tachycardia. Tests boundary conditions and circular wave propagation.
Stable 3D spiral wave dynamics and multidirectional propagation. Validates 3D behavior before anatomical complexity.
Anatomical features: pulmonary veins, mitral valve, Bachmann's bundle. Tests realistic geometry before patient data.
Real left atrial geometries from medical imaging (CTA, MRI). Clinically relevant models for arrhythmia reproduction and ablation planning.
The digital twin successfully reproduced multiple clinically relevant arrhythmia mechanisms observed in patients. By adjusting conduction velocities and applying precisely timed ectopic events, we induced various forms of reentrant activity including micro-reentry, macro-reentry, and atrial flutter. These arrhythmias demonstrated self-sustaining propagation patterns consistent with clinical observations.
We successfully induced left atrial flutter by slowing conduction velocity between pulmonary vein openings and triggering ectopic events at specific locations and times. Reentrant circuits and rotors emerged naturally, producing sustained activation patterns that resembled those seen in real patients.
Simulated catheter ablation reliably terminated reentrant activity and restored organized conduction patterns. We successfully eliminated induced flutters using virtual ablations, allowing the system to return to normal rhythm. The model demonstrated effectiveness on both idealized and patient-specific atrial geometries.
The interactive GUI lets you trigger single ectopic events, schedule recurring ectopic beats, perform virtual ablations, and adjust beat period and visualization—all in real time without digging through menus.
Principal Investigator:
Dr. Bryant Wyatt
Tarleton State University, Department of Mathematics
Co-Investigator:
Dr. Charles Puelz
University of Houston, Department of Mathematics
Accepted Poster Presentation
1st Place Poster
1st Place Undergraduate Poster
Heart Rhythm O2, 6(9):S2
Published: September 2025
Journal of Electrocardiology, 86 (2024), 153762
Published: July 22, 2024 | Open Access
Scroll through full-size media without squinting. Use the arrows or thumbnails to jump between videos and images; everything plays inline.