A polyconvex computational formulation for electro-activation in cardiac mechanics

Blanco, Emilio Garcia

doi:10.23889/Suthesis.52476

E-Thesis 776 views

A polyconvex computational formulation for electro-activation in cardiac mechanics / Emilio Garcia Blanco

Swansea University Author: Emilio Garcia Blanco

Full text not available from this repository: check for access using links below.

DOI (Published version): 10.23889/Suthesis.52476

Abstract

Cardiovascular diseases ( CVD) represent the main cause of death in the world with more significant prevalence in developed countries. With a constantly increasing ageing population, the burden that CVD patients put on the healthcare system has been recognised as an urgent matter in need of immediat...

Full description

Published:	Swansea 2019
Institution:	Swansea University
Degree level:	Doctoral
Degree name:	Ph.D
URI:	https://cronfa.swan.ac.uk/Record/cronfa52476
Tags:	Add Tag No Tags, Be the first to tag this record!

first_indexed	2019-10-16T20:25:39Z
last_indexed	2023-01-11T14:29:38Z
id	cronfa52476
recordtype	RisThesis
fullrecord	<?xml version="1.0"?><rfc1807><datestamp>2022-12-18T09:55:55.4009526</datestamp><bib-version>v2</bib-version><id>52476</id><entry>2019-10-16</entry><title>A polyconvex computational formulation for electro-activation in cardiac mechanics</title><swanseaauthors><author><sid>390d65783b3d69393b85ba8adfe688ec</sid><firstname>Emilio</firstname><surname>Garcia Blanco</surname><name>Emilio Garcia Blanco</name><active>true</active><ethesisStudent>false</ethesisStudent></author></swanseaauthors><date>2019-10-16</date><deptcode>FGSEN</deptcode><abstract>Cardiovascular diseases ( CVD) represent the main cause of death in the world with more significant prevalence in developed countries. With a constantly increasing ageing population, the burden that CVD patients put on the healthcare system has been recognised as an urgent matter in need of immediate attention. With the aim of addressing this menace, the computational modelling of the complex physical phenomena occurring in the human heart has become an area of increasing scientific interest over the last decades. The support that computational mechanics can provide through the design of in-silica diagnostic tools is nowadays well-acknowledged by experts in the field, especially in challenging cardiopathies such as heart infarction or dysrhythmia. Even more, great effort has been devoted to attempting the computational modelling of CVD with several objectives in mind: provide augmented simulation-based diagnosis, experiment innovative surgery techniques, offer patients personalised treatments and to better understand the very complex electro-chemobio-mechanical phenomena underpinning the behaviour of the heart. Unfortunately, most of the computational frameworks considered in the field are unable to account for the material complexity of the myocardium (heterogeneity, anisotropy, nearlyincompressibility, etc.), exhibiting spurious pressure oscillations, shear locking and volumetric locking. The main objective of this thesis is twofold: on the one hand, a novel computational framework for the numerical simulation of the electromechanical response of the myocardium during the cardiac cycle is presented; on the other hand, this novel framework is used as a reference scenario to explore alternative strategies and establish a balance between their applicability and efficiency. First, taking as a reference the mixed formulations introduced by Bonet et al. [21] in the context of nonlinear elasticity, two new mixed formulations designed for polyconvex materials have been developed and customised for active stress and active strain coupling approaches. These proposed formulations add considerable flexibility to the modelling of the heart by including the geometry and the transmembrane potential ( and possibly a Lagrange multiplier enforcing weakly the incompressibility constraint) as unknown fields as well as the deformation gradient tensor, its cofactor, its determinant, the gradient of the transmembrane potential and their respective work conjugates. Their Finite Element implementation is explained, where a static condensation procedure is presented in order to yield an extremely competitive computational approach. The superiority in accuracy of these formulations with respect to classical low order implementations is enhanced by a monolithic fully-implicit scheme for the electromechanical problem. The remarkable robustness and accuracy of these top-notch mixed formulations creates an ideal scenario to carry out a comprehensive and rigorous study of different ionic models and electromechanical activation couplings. For first time, large scale simulations have been performed in order to set up a four-way comparison among the ionic models proposed by Bueno-Orovio and Ten Tusscher, and the active strain and active stress coupling approaches under the same mesh and boundary conditions. In addition, these simulations constitute a suitable benchmark to analyse the possible loss of ellipticity and polyconvexity of the Holzapfel-Ogden constitutive model, widely used in the context of cardiac mechanics, putting forward possible polyconvexifications. Moreover, an invariant representation of Guccione's constitutive model is proposed. Secondly, despite the correct mathematical characterisation of the myocardium, the complex multi-scale interaction that takes place between the electrophysiology and electrochemistry at cellular level and the macro-scale response of the heart muscle customarily require immense computational resources. This poses an extreme challenge to the scalability of electro-mechanical solvers due to the size and conditioning of the system of equations required to obtain accurate solutions, both in terms of wall deformation and transmembrane potential propagation. In order to undertake this challenge, the results obtained by the polyconvex mixed formulations represent an ideal reference scenario to evaluate the applicability of more computationally efficient techniques and cheaper Finite Element discretisations. In the search towards an efficient modelling of electro-activation, this manuscript presents a coupled electromechanical computational framework whereby, first, the use of an efficient stabilised low order tetrahedral Finite Element methodology is explored and compared against the very accurate super enhanced mixed formulation previously introduced and, second, the use of tailor-made staggered and staggered linearised solvers is implemented in order to assess their feasibility against a fully monolithic approach. These comparisons take the global displacements and intermembrane potential as accuracy indicators as well as considering the overall computational time. Through a comprehensive set of examples, culminating in a realistic ventricular geometry, we aim to put forward some suggestions regarding the level of discretisation and coupling required to ensure sufficiently reliable results yet with an affordable computational time. Additionally, the possible locking phenomena due to shear deformations and anisotropic behaviour is assessed for our enriched mixed formulation and compared against less expensive methodologies.</abstract><type>E-Thesis</type><journal/><volume/><journalNumber/><paginationStart/><paginationEnd/><publisher/><placeOfPublication>Swansea</placeOfPublication><isbnPrint/><isbnElectronic/><issnPrint/><issnElectronic/><keywords>Cardiac electromechanics, Mixed Formulations, Polyconvexity, Finite Elements</keywords><publishedDay>31</publishedDay><publishedMonth>12</publishedMonth><publishedYear>2019</publishedYear><publishedDate>2019-12-31</publishedDate><doi>10.23889/Suthesis.52476</doi><url/><notes>A selection of third party content is redacted or is partially redacted from this thesis.</notes><college>COLLEGE NANME</college><department>Science and Engineering - Faculty</department><CollegeCode>COLLEGE CODE</CollegeCode><DepartmentCode>FGSEN</DepartmentCode><institution>Swansea University</institution><degreelevel>Doctoral</degreelevel><degreename>Ph.D</degreename><apcterm/><funders/><projectreference/><lastEdited>2022-12-18T09:55:55.4009526</lastEdited><Created>2019-10-16T14:57:16.5609500</Created><path><level id="1">Faculty of Science and Engineering</level><level id="2">School of Engineering and Applied Sciences - Uncategorised</level></path><authors><author><firstname>Emilio</firstname><surname>Garcia Blanco</surname><order>1</order></author></authors><documents/><OutputDurs/></rfc1807>
spelling	2022-12-18T09:55:55.4009526 v2 52476 2019-10-16 A polyconvex computational formulation for electro-activation in cardiac mechanics 390d65783b3d69393b85ba8adfe688ec Emilio Garcia Blanco Emilio Garcia Blanco true false 2019-10-16 FGSEN Cardiovascular diseases ( CVD) represent the main cause of death in the world with more significant prevalence in developed countries. With a constantly increasing ageing population, the burden that CVD patients put on the healthcare system has been recognised as an urgent matter in need of immediate attention. With the aim of addressing this menace, the computational modelling of the complex physical phenomena occurring in the human heart has become an area of increasing scientific interest over the last decades. The support that computational mechanics can provide through the design of in-silica diagnostic tools is nowadays well-acknowledged by experts in the field, especially in challenging cardiopathies such as heart infarction or dysrhythmia. Even more, great effort has been devoted to attempting the computational modelling of CVD with several objectives in mind: provide augmented simulation-based diagnosis, experiment innovative surgery techniques, offer patients personalised treatments and to better understand the very complex electro-chemobio-mechanical phenomena underpinning the behaviour of the heart. Unfortunately, most of the computational frameworks considered in the field are unable to account for the material complexity of the myocardium (heterogeneity, anisotropy, nearlyincompressibility, etc.), exhibiting spurious pressure oscillations, shear locking and volumetric locking. The main objective of this thesis is twofold: on the one hand, a novel computational framework for the numerical simulation of the electromechanical response of the myocardium during the cardiac cycle is presented; on the other hand, this novel framework is used as a reference scenario to explore alternative strategies and establish a balance between their applicability and efficiency. First, taking as a reference the mixed formulations introduced by Bonet et al. [21] in the context of nonlinear elasticity, two new mixed formulations designed for polyconvex materials have been developed and customised for active stress and active strain coupling approaches. These proposed formulations add considerable flexibility to the modelling of the heart by including the geometry and the transmembrane potential ( and possibly a Lagrange multiplier enforcing weakly the incompressibility constraint) as unknown fields as well as the deformation gradient tensor, its cofactor, its determinant, the gradient of the transmembrane potential and their respective work conjugates. Their Finite Element implementation is explained, where a static condensation procedure is presented in order to yield an extremely competitive computational approach. The superiority in accuracy of these formulations with respect to classical low order implementations is enhanced by a monolithic fully-implicit scheme for the electromechanical problem. The remarkable robustness and accuracy of these top-notch mixed formulations creates an ideal scenario to carry out a comprehensive and rigorous study of different ionic models and electromechanical activation couplings. For first time, large scale simulations have been performed in order to set up a four-way comparison among the ionic models proposed by Bueno-Orovio and Ten Tusscher, and the active strain and active stress coupling approaches under the same mesh and boundary conditions. In addition, these simulations constitute a suitable benchmark to analyse the possible loss of ellipticity and polyconvexity of the Holzapfel-Ogden constitutive model, widely used in the context of cardiac mechanics, putting forward possible polyconvexifications. Moreover, an invariant representation of Guccione's constitutive model is proposed. Secondly, despite the correct mathematical characterisation of the myocardium, the complex multi-scale interaction that takes place between the electrophysiology and electrochemistry at cellular level and the macro-scale response of the heart muscle customarily require immense computational resources. This poses an extreme challenge to the scalability of electro-mechanical solvers due to the size and conditioning of the system of equations required to obtain accurate solutions, both in terms of wall deformation and transmembrane potential propagation. In order to undertake this challenge, the results obtained by the polyconvex mixed formulations represent an ideal reference scenario to evaluate the applicability of more computationally efficient techniques and cheaper Finite Element discretisations. In the search towards an efficient modelling of electro-activation, this manuscript presents a coupled electromechanical computational framework whereby, first, the use of an efficient stabilised low order tetrahedral Finite Element methodology is explored and compared against the very accurate super enhanced mixed formulation previously introduced and, second, the use of tailor-made staggered and staggered linearised solvers is implemented in order to assess their feasibility against a fully monolithic approach. These comparisons take the global displacements and intermembrane potential as accuracy indicators as well as considering the overall computational time. Through a comprehensive set of examples, culminating in a realistic ventricular geometry, we aim to put forward some suggestions regarding the level of discretisation and coupling required to ensure sufficiently reliable results yet with an affordable computational time. Additionally, the possible locking phenomena due to shear deformations and anisotropic behaviour is assessed for our enriched mixed formulation and compared against less expensive methodologies. E-Thesis Swansea Cardiac electromechanics, Mixed Formulations, Polyconvexity, Finite Elements 31 12 2019 2019-12-31 10.23889/Suthesis.52476 A selection of third party content is redacted or is partially redacted from this thesis. COLLEGE NANME Science and Engineering - Faculty COLLEGE CODE FGSEN Swansea University Doctoral Ph.D 2022-12-18T09:55:55.4009526 2019-10-16T14:57:16.5609500 Faculty of Science and Engineering School of Engineering and Applied Sciences - Uncategorised Emilio Garcia Blanco 1
title	A polyconvex computational formulation for electro-activation in cardiac mechanics
spellingShingle	A polyconvex computational formulation for electro-activation in cardiac mechanics Emilio Garcia Blanco
title_short	A polyconvex computational formulation for electro-activation in cardiac mechanics
title_full	A polyconvex computational formulation for electro-activation in cardiac mechanics
title_fullStr	A polyconvex computational formulation for electro-activation in cardiac mechanics
title_full_unstemmed	A polyconvex computational formulation for electro-activation in cardiac mechanics
title_sort	A polyconvex computational formulation for electro-activation in cardiac mechanics
author_id_str_mv	390d65783b3d69393b85ba8adfe688ec
author_id_fullname_str_mv	390d65783b3d69393b85ba8adfe688ec_***_Emilio Garcia Blanco
author	Emilio Garcia Blanco
author2	Emilio Garcia Blanco
format	E-Thesis
publishDate	2019
institution	Swansea University
doi_str_mv	10.23889/Suthesis.52476
college_str	Faculty of Science and Engineering
hierarchytype
hierarchy_top_id	facultyofscienceandengineering
hierarchy_top_title	Faculty of Science and Engineering
hierarchy_parent_id	facultyofscienceandengineering
hierarchy_parent_title	Faculty of Science and Engineering
department_str	School of Engineering and Applied Sciences - Uncategorised{{{_:::_}}}Faculty of Science and Engineering{{{_:::_}}}School of Engineering and Applied Sciences - Uncategorised
document_store_str	0
active_str	0
description	Cardiovascular diseases ( CVD) represent the main cause of death in the world with more significant prevalence in developed countries. With a constantly increasing ageing population, the burden that CVD patients put on the healthcare system has been recognised as an urgent matter in need of immediate attention. With the aim of addressing this menace, the computational modelling of the complex physical phenomena occurring in the human heart has become an area of increasing scientific interest over the last decades. The support that computational mechanics can provide through the design of in-silica diagnostic tools is nowadays well-acknowledged by experts in the field, especially in challenging cardiopathies such as heart infarction or dysrhythmia. Even more, great effort has been devoted to attempting the computational modelling of CVD with several objectives in mind: provide augmented simulation-based diagnosis, experiment innovative surgery techniques, offer patients personalised treatments and to better understand the very complex electro-chemobio-mechanical phenomena underpinning the behaviour of the heart. Unfortunately, most of the computational frameworks considered in the field are unable to account for the material complexity of the myocardium (heterogeneity, anisotropy, nearlyincompressibility, etc.), exhibiting spurious pressure oscillations, shear locking and volumetric locking. The main objective of this thesis is twofold: on the one hand, a novel computational framework for the numerical simulation of the electromechanical response of the myocardium during the cardiac cycle is presented; on the other hand, this novel framework is used as a reference scenario to explore alternative strategies and establish a balance between their applicability and efficiency. First, taking as a reference the mixed formulations introduced by Bonet et al. [21] in the context of nonlinear elasticity, two new mixed formulations designed for polyconvex materials have been developed and customised for active stress and active strain coupling approaches. These proposed formulations add considerable flexibility to the modelling of the heart by including the geometry and the transmembrane potential ( and possibly a Lagrange multiplier enforcing weakly the incompressibility constraint) as unknown fields as well as the deformation gradient tensor, its cofactor, its determinant, the gradient of the transmembrane potential and their respective work conjugates. Their Finite Element implementation is explained, where a static condensation procedure is presented in order to yield an extremely competitive computational approach. The superiority in accuracy of these formulations with respect to classical low order implementations is enhanced by a monolithic fully-implicit scheme for the electromechanical problem. The remarkable robustness and accuracy of these top-notch mixed formulations creates an ideal scenario to carry out a comprehensive and rigorous study of different ionic models and electromechanical activation couplings. For first time, large scale simulations have been performed in order to set up a four-way comparison among the ionic models proposed by Bueno-Orovio and Ten Tusscher, and the active strain and active stress coupling approaches under the same mesh and boundary conditions. In addition, these simulations constitute a suitable benchmark to analyse the possible loss of ellipticity and polyconvexity of the Holzapfel-Ogden constitutive model, widely used in the context of cardiac mechanics, putting forward possible polyconvexifications. Moreover, an invariant representation of Guccione's constitutive model is proposed. Secondly, despite the correct mathematical characterisation of the myocardium, the complex multi-scale interaction that takes place between the electrophysiology and electrochemistry at cellular level and the macro-scale response of the heart muscle customarily require immense computational resources. This poses an extreme challenge to the scalability of electro-mechanical solvers due to the size and conditioning of the system of equations required to obtain accurate solutions, both in terms of wall deformation and transmembrane potential propagation. In order to undertake this challenge, the results obtained by the polyconvex mixed formulations represent an ideal reference scenario to evaluate the applicability of more computationally efficient techniques and cheaper Finite Element discretisations. In the search towards an efficient modelling of electro-activation, this manuscript presents a coupled electromechanical computational framework whereby, first, the use of an efficient stabilised low order tetrahedral Finite Element methodology is explored and compared against the very accurate super enhanced mixed formulation previously introduced and, second, the use of tailor-made staggered and staggered linearised solvers is implemented in order to assess their feasibility against a fully monolithic approach. These comparisons take the global displacements and intermembrane potential as accuracy indicators as well as considering the overall computational time. Through a comprehensive set of examples, culminating in a realistic ventricular geometry, we aim to put forward some suggestions regarding the level of discretisation and coupling required to ensure sufficiently reliable results yet with an affordable computational time. Additionally, the possible locking phenomena due to shear deformations and anisotropic behaviour is assessed for our enriched mixed formulation and compared against less expensive methodologies.
published_date	2019-12-31T04:04:52Z
_version_	1763753372983230464
score	11.036706

A polyconvex computational formulation for electro-activation in cardiac mechanics / Emilio Garcia Blanco

Similar Items