<p>Computer models for simulating cardiac electrophysiology are valuable tools for research and clinical applications. Traditional reaction–diffusion (RD) models used for these purposes are computationally expensive. While eikonal models offer a faster alternative, they are not well-suited to study cardiac arrhythmias driven by reentrant activity. The present work extends the diffusion–reaction eikonal alternant model (DREAM), incorporating conduction velocity (CV) restitution for simulating complex cardiac arrhythmias. The DREAM modifies the fast iterative method to model cyclical behavior, dynamic boundary conditions, and frequency-dependent anisotropic CV. Additionally, the model alternates with an approximated RD model, using a detailed ionic model for the reaction term and a triple-Gaussian to approximate the diffusion term. The DREAM and monodomain models were compared, simulating reentries in 2D manifolds with different resolutions. The DREAM produced similar results across all resolutions, while experiments with the monodomain model failed at lower resolutions. CV restitution curves obtained using the DREAM closely approximated those produced by the monodomain simulations. Reentry in 2D sheets yielded similar results in vulnerable window and mean reentry duration for low CV in both models. In the left atrium, most inducing points identified by the DREAM were also present in the high-resolution monodomain model. DREAM’s reentry simulations on meshes with an average edge length of <InlineEquation ID="IEq1"> <InlineMediaObject> <ImageObject Color="BlackWhite" FileRef="366_2024_2094_Article_IEq1.gif" Format="GIF" Height="16" Rendition="HTML" Resolution="72" Type="Linedraw" Width="60" /> </InlineMediaObject> <EquationSource Format="TEX">\(1600\,\upmu \hbox {m}\)</EquationSource> <EquationSource Format="MATHML"><math> <mrow> <mn>1600</mn> <mspace width="0.166667em" /> <mi mathvariant="normal">μ</mi> <mtext>m</mtext> </mrow> </math></EquationSource> </InlineEquation> were 87x faster than monodomain simulations at <InlineEquation ID="IEq2"> <InlineMediaObject> <ImageObject Color="BlackWhite" FileRef="366_2024_2094_Article_IEq2.gif" Format="GIF" Height="16" Rendition="HTML" Resolution="72" Type="Linedraw" Width="52" /> </InlineMediaObject> <EquationSource Format="TEX">\(200\,\upmu \hbox {m}\)</EquationSource> <EquationSource Format="MATHML"><math> <mrow> <mn>200</mn> <mspace width="0.166667em" /> <mi mathvariant="normal">μ</mi> <mtext>m</mtext> </mrow> </math></EquationSource> </InlineEquation>. This work establishes the mathematical foundation for using the accelerated DREAM simulation method for cardiac electrophysiology. Cardiac research applications are enabled by a publicly available implementation in the openCARP simulator.</p>

错误:搜索内容不能为空,请输入英文关键词
错误:关键词超出字数限制,请精简
高级检索

A cyclical fast iterative method for simulating reentries in cardiac electrophysiology using an eikonal-based model

  • Cristian Barrios Espinosa,
  • Jorge Sánchez,
  • Stephanie Appel,
  • Silvia Becker,
  • Jonathan Krauß,
  • Patricia Martínez Díaz,
  • Laura Unger,
  • Marie Houillon,
  • Axel Loewe

摘要

Computer models for simulating cardiac electrophysiology are valuable tools for research and clinical applications. Traditional reaction–diffusion (RD) models used for these purposes are computationally expensive. While eikonal models offer a faster alternative, they are not well-suited to study cardiac arrhythmias driven by reentrant activity. The present work extends the diffusion–reaction eikonal alternant model (DREAM), incorporating conduction velocity (CV) restitution for simulating complex cardiac arrhythmias. The DREAM modifies the fast iterative method to model cyclical behavior, dynamic boundary conditions, and frequency-dependent anisotropic CV. Additionally, the model alternates with an approximated RD model, using a detailed ionic model for the reaction term and a triple-Gaussian to approximate the diffusion term. The DREAM and monodomain models were compared, simulating reentries in 2D manifolds with different resolutions. The DREAM produced similar results across all resolutions, while experiments with the monodomain model failed at lower resolutions. CV restitution curves obtained using the DREAM closely approximated those produced by the monodomain simulations. Reentry in 2D sheets yielded similar results in vulnerable window and mean reentry duration for low CV in both models. In the left atrium, most inducing points identified by the DREAM were also present in the high-resolution monodomain model. DREAM’s reentry simulations on meshes with an average edge length of \(1600\,\upmu \hbox {m}\) 1600 μ m were 87x faster than monodomain simulations at \(200\,\upmu \hbox {m}\) 200 μ m . This work establishes the mathematical foundation for using the accelerated DREAM simulation method for cardiac electrophysiology. Cardiac research applications are enabled by a publicly available implementation in the openCARP simulator.