An improved numerical method for strong coupling of excitation and contraction models in the heart
Quantifying the interactions between excitation and contraction is fundamental to furthering our understanding of cardiac physiology. To date simulating these effects in strongly coupled excitation and contraction tissue models has proved computationally challenging. This is in part due to the numer...
| Main Authors: | , |
|---|---|
| Format: | Journal article |
| Published: |
2008
|
| _version_ | 1797052128117653504 |
|---|---|
| author | Niederer, S Smith, N |
| author_facet | Niederer, S Smith, N |
| author_sort | Niederer, S |
| collection | OXFORD |
| description | Quantifying the interactions between excitation and contraction is fundamental to furthering our understanding of cardiac physiology. To date simulating these effects in strongly coupled excitation and contraction tissue models has proved computationally challenging. This is in part due to the numerical methods implemented to maintain numerical stability in previous simulations, which produced computationally intensive problems. In this study, we analytically identify and quantify the velocity and length dependent sources of instability in the current established coupling method and propose a new method which addresses these issues. Specifically, we account for the length and velocity dependence of active tension within the finite deformation equations, such that the active tension is updated at each intermediate Newton iteration, within the mechanics solution step. We then demonstrate that the model is stable and converges in a three-dimensional rod under isometric contraction. Subsequently, we show that the coupling method can produce stable solutions in a cube with many of the attributes present in the heart, including asymmetrical activation, an inhomogeneous fibre field and a nonlinear constitutive law. The results show no instabilities and quantify the error introduced by discrete length updates. This validates our proposed coupling framework, demonstrating significant improvement in the stability of excitation and contraction simulations. |
| first_indexed | 2024-03-06T18:28:22Z |
| format | Journal article |
| id | oxford-uuid:08c4b7a4-cba1-40c4-8a94-01e986e7ffd3 |
| institution | University of Oxford |
| last_indexed | 2024-03-06T18:28:22Z |
| publishDate | 2008 |
| record_format | dspace |
| spelling | oxford-uuid:08c4b7a4-cba1-40c4-8a94-01e986e7ffd32022-03-26T09:14:41ZAn improved numerical method for strong coupling of excitation and contraction models in the heartJournal articlehttp://purl.org/coar/resource_type/c_dcae04bcuuid:08c4b7a4-cba1-40c4-8a94-01e986e7ffd3Department of Computer Science2008Niederer, SSmith, NQuantifying the interactions between excitation and contraction is fundamental to furthering our understanding of cardiac physiology. To date simulating these effects in strongly coupled excitation and contraction tissue models has proved computationally challenging. This is in part due to the numerical methods implemented to maintain numerical stability in previous simulations, which produced computationally intensive problems. In this study, we analytically identify and quantify the velocity and length dependent sources of instability in the current established coupling method and propose a new method which addresses these issues. Specifically, we account for the length and velocity dependence of active tension within the finite deformation equations, such that the active tension is updated at each intermediate Newton iteration, within the mechanics solution step. We then demonstrate that the model is stable and converges in a three-dimensional rod under isometric contraction. Subsequently, we show that the coupling method can produce stable solutions in a cube with many of the attributes present in the heart, including asymmetrical activation, an inhomogeneous fibre field and a nonlinear constitutive law. The results show no instabilities and quantify the error introduced by discrete length updates. This validates our proposed coupling framework, demonstrating significant improvement in the stability of excitation and contraction simulations. |
| spellingShingle | Niederer, S Smith, N An improved numerical method for strong coupling of excitation and contraction models in the heart |
| title | An improved numerical method for strong coupling of excitation and contraction models in the heart |
| title_full | An improved numerical method for strong coupling of excitation and contraction models in the heart |
| title_fullStr | An improved numerical method for strong coupling of excitation and contraction models in the heart |
| title_full_unstemmed | An improved numerical method for strong coupling of excitation and contraction models in the heart |
| title_short | An improved numerical method for strong coupling of excitation and contraction models in the heart |
| title_sort | improved numerical method for strong coupling of excitation and contraction models in the heart |
| work_keys_str_mv | AT niederers animprovednumericalmethodforstrongcouplingofexcitationandcontractionmodelsintheheart AT smithn animprovednumericalmethodforstrongcouplingofexcitationandcontractionmodelsintheheart AT niederers improvednumericalmethodforstrongcouplingofexcitationandcontractionmodelsintheheart AT smithn improvednumericalmethodforstrongcouplingofexcitationandcontractionmodelsintheheart |