Modelling mitral valvular dynamics–current trend and future directions

Dysfunction of mitral valve causes morbidity and premature mortality and remains a leading medical problem worldwide. Computational modelling aims to understand the biomechanics of human mitral valve and could lead to the development of new treatment, prevention and diagnosis of mitral valve diseases. Compared with the aortic valve, the mitral valve has been much less studied owing to its highly complex structure and strong interaction with the blood flow and the ventricles. However, the interest in mitral valve modelling is growing, and the sophistication level is increasing with the advanced development of computational technology and imaging tools. This review summarises the state-of-the-art modelling of the mitral valve, including static and dynamics models, models with fluid-structure interaction, and models with the left ventricle interaction. Challenges and future directions are also discussed

Quasi-static imaged-based immersed boundary-finite element model of human left ventricle in diastole

SUMMARY: Finite stress and strain analyses of the heart provide insight into the biomechanics of myocardial function and dysfunction. Herein, we describe progress toward dynamic patient-specific models of the left ventricle using an immersed boundary (IB) method with a finite element (FE) structural mechanics model. We use a structure-based hyperelastic strain-energy function to describe the passive mechanics of the ventricular myocardium, a realistic anatomical geometry reconstructed from clinical magnetic resonance images of a healthy human heart, and a rule-based fiber architecture. Numerical predictions of this IB/FE model are compared with results obtained by a commercial FE solver. We demonstrate that the IB/FE model yields results that are in good agreement with those of the conventional FE model under diastolic loading conditions, and the predictions of the LV model using either numerical method are shown to be consistent with previous computational and experimental data. These results are among the first to analyze the stress and strain predictions of IB models of ventricular mechanics, and they serve both to verify the IB/FE simulation framework and to validate the IB/FE model. Moreover, this work represents an important step toward using such models for fully dynamic fluid–structure interaction simulations of the heart

Towards patient-specific modelling as a pre-operative planning strategy and follow up assessment for the treatment of advanced heart failure with rotary blood pumps

Background: Ventricular Assist Devices (VADs) insertion is an established treatment for patients with end-stage heart failure waiting for a heart transplant or in need for long-term circulatory support (destination therapy). Rotary blood pumps (RBP) are the most popular devices in view of their size and performance. Pre-operative planning strategy for the insertion of a left ventricular assist device (LVAD) requires a timely discussion at a Multi-Disciplinary Team Meeting (MDT). Clinical-decision making is based according to the needs of the patient and must be processed without delays. Nevertheless, thrombus formation remains a feared complication which affects outcome. VADs operate in a flow regime which is difficult to simulate: the transitional region at the boundary of laminar and turbulent flow (low Reynolds number). Different methods have been used but the best approach remains debatable. Computational Fluid Dynamics (CFD) is an attractive and invaluable tool for the study of the interactions between VADs and the cardiovascular system. The aim of this thesis is three-fold: a) to investigate the use of pressure-volume analysis in a clinical setting through the review of six heart failure patients previously discussed at a MDT meeting with a view to predict or guide further management; b) to review the theory behind modelling approaches to VADs and their interactions with the cardiovascular system for better understanding of their clinical use. Then, an overview of computational fluid dynamics (CFD) is considered as a prelude to its application to the analysis of VADs performance. Additionally, the development of a simplified model of centrifugal pump will be used in initial simulations as preliminary analysis; c) to examine an example of a proof-of-concept pilot patient-specific model of an axial flow pump (HeartMate II) as pre-operative planning strategy in a patient-specific model with a view to identify potential critical areas that may affect pump function and outcome in a clinical setting. Material and Methods: 3D reconstruction from CT-scan images of patients who underwent the insertion of rotary blood pumps, namely HeartWare HVAD and HeartMate II. Ansys Fluent has been used for CFD analysis based on the fundamental governing equations of motion. Blood has been modelled as incompressible, Newtonian fluid with density = 1060 and viscosity = 0.0035 kg/m-s. The laminar and SST models have been used for comparison purposes. The rotational motion of the impeller has been implemented using the moving reference frame (MRF) approach. The sliding mesh method has also been used to account for unsteady interaction between stationary and moving part. The no-slip condition has been applied to all walls, which were assumed to be rigid. Boundary conditions consisting of velocity inlet and pressure outlet of the pump based on different settings and constant rotational speed for the impeller. Pressure-velocity coupling has been based on the coupled scheme. Spatial discretisation consisted of the “least square cell based” gradient for velocity and “PRESTO” or second order for pressure. Second order upwind has been set for the momentum, turbulent kinetic energy and specific dissipation rate. First order implicit has been set for transient formulation. The pseudo transient algorithm (steady state), the high order relaxation term and the warped-face gradient correction have been used to add an unsteady term to the solution equations with the aim to improve stability and enhance convergence. Specific settings have been considered for comparison purposes. Results: Pressure-volume simulation analysis in six advanced heart failure patients showed that an integrated model of the cardiovascular system based on lumped-parameter representation, modified time-varying elastance and pressure-volume analysis of ventricular function seems a feasible and suitable approach yielding a sufficiently accurate quantitative analysis in real time, therefore applicable within the time-constraints of a clinical setting. Lumped-parameter models consist of simultaneous ordinary differential equations complemented by an algebraic balance equation and are suitable for examination of global distribution of pressure, flow and volume over a range of physiological conditions with inclusion of the interaction between modelled components. Higher level lumped-parameter modelling is needed to address the interaction between the circulation and other systems based on a compromise between complexity and ability to set the required parameters to personalise an integrated lumped-parameter model for a patient-specific approach. CARDIOSIM© fulfils these requirements and does address the systems interaction with its modular approach and assembly of models with varying degree of complexity although 0-D and 1-D coupling may be required for the evaluation of long-term VAD support. The challenge remains the ability to predict outcome over a longer period of time. The preliminary CFD simulations with the HeartWare HVAD centrifugal pump demonstrated that it is possible to obtain an accurate analysis in a timely manner to complement the clinical review process. The simulations with the pilot patient-specific model of the HeartMate II axial flow pump revealed that a complex 3D reconstruction is feasible in a timely manner and can be used to generate sufficiently accurate results to be used in the context of a MDT meeting for the purposes of clinical decision-making. Overall, these three studies demonstrate that the time frame of the simulations was within hours which may fit the time constraints of the clinical environment in the context of a MDT meeting. More specifically, it was shown that the laminar model may be used for an initial evaluation of the flow development within the pump. Nonetheless, the k- model offers higher accuracy if the timeline of the clinical setting allows for a longer simulation. Conclusion: This thesis aimed at the understanding of the use of computational modelling as a pre-operative planning strategy and follow up assessment for the treatment of advanced heart failure with rotary blood pumps. The novelty lays in the use of both pressure-volume simulation analysis and 3D flow dynamics studies in VADs with a view to treatment optimisation and outcome prediction within the time constraints of a clinical setting in the context of a MDT meeting. The clinical significance and the contribution to the field is a more targeted approach for different groups of patients and a more quantitative evaluation in the clinical decision process based on a pro-active co-operation between clinicians and scientists reducing the potential for “guess work”. The results of this thesis are a proof-of-concept as a prelude to a potential future implementation of patient-specific modelling within a clinical setting on a daily basis demonstrating a clear clinical significance and contribution to the field. The proposed approach does not consider modelling and simulation as a substitute for clinical experience but an additional tool to guide therapeutic intervention and complement the clinical decision process in which the clinician remains the ultimate decision-maker. Such an approach may well add a different dimension to the problem of heart failure with potential for high return in terms of patient’s outcome and long-term surveillance. The same principles would be applicable to other cardiovascular problems in line with the current concept of “Team Approach” such as the Heart Team, the Structural Heart Team or the Aortic Team. The present work has taken this concept closer to clinical delivery and has highlighted its potential but further work remains to be done in refining the technique.Background: Ventricular Assist Devices (VADs) insertion is an established treatment for patients with end-stage heart failure waiting for a heart transplant or in need for long-term circulatory support (destination therapy). Rotary blood pumps (RBP) are the most popular devices in view of their size and performance. Pre-operative planning strategy for the insertion of a left ventricular assist device (LVAD) requires a timely discussion at a Multi-Disciplinary Team Meeting (MDT). Clinical-decision making is based according to the needs of the patient and must be processed without delays. Nevertheless, thrombus formation remains a feared complication which affects outcome. VADs operate in a flow regime which is difficult to simulate: the transitional region at the boundary of laminar and turbulent flow (low Reynolds number). Different methods have been used but the best approach remains debatable. Computational Fluid Dynamics (CFD) is an attractive and invaluable tool for the study of the interactions between VADs and the cardiovascular system. The aim of this thesis is three-fold: a) to investigate the use of pressure-volume analysis in a clinical setting through the review of six heart failure patients previously discussed at a MDT meeting with a view to predict or guide further management; b) to review the theory behind modelling approaches to VADs and their interactions with the cardiovascular system for better understanding of their clinical use. Then, an overview of computational fluid dynamics (CFD) is considered as a prelude to its application to the analysis of VADs performance. Additionally, the development of a simplified model of centrifugal pump will be used in initial simulations as preliminary analysis; c) to examine an example of a proof-of-concept pilot patient-specific model of an axial flow pump (HeartMate II) as pre-operative planning strategy in a patient-specific model with a view to identify potential critical areas that may affect pump function and outcome in a clinical setting. Material and Methods: 3D reconstruction from CT-scan images of patients who underwent the insertion of rotary blood pumps, namely HeartWare HVAD and HeartMate II. Ansys Fluent has been used for CFD analysis based on the fundamental governing equations of motion. Blood has been modelled as incompressible, Newtonian fluid with density = 1060 and viscosity = 0.0035 kg/m-s. The laminar and SST models have been used for comparison purposes. The rotational motion of the impeller has been implemented using the moving reference frame (MRF) approach. The sliding mesh method has also been used to account for unsteady interaction between stationary and moving part. The no-slip condition has been applied to all walls, which were assumed to be rigid. Boundary conditions consisting of velocity inlet and pressure outlet of the pump based on different settings and constant rotational speed for the impeller. Pressure-velocity coupling has been based on the coupled scheme. Spatial discretisation consisted of the “least square cell based” gradient for velocity and “PRESTO” or second order for pressure. Second order upwind has been set for the momentum, turbulent kinetic energy and specific dissipation rate. First order implicit has been set for transient formulation. The pseudo transient algorithm (steady state), the high order relaxation term and the warped-face gradient correction have been used to add an unsteady term to the solution equations with the aim to improve stability and enhance convergence. Specific settings have been considered for comparison purposes. Results: Pressure-volume simulation analysis in six advanced heart failure patients showed that an integrated model of the cardiovascular system based on lumped-parameter representation, modified time-varying elastance and pressure-volume analysis of ventricular function seems a feasible and suitable approach yielding a sufficiently accurate quantitative analysis in real time, therefore applicable within the time-constraints of a clinical setting. Lumped-parameter models consist of simultaneous ordinary differential equations complemented by an algebraic balance equation and are suitable for examination of global distribution of pressure, flow and volume over a range of physiological conditions with inclusion of the interaction between modelled components. Higher level lumped-parameter modelling is needed to address the interaction between the circulation and other systems based on a compromise between complexity and ability to set the required parameters to personalise an integrated lumped-parameter model for a patient-specific approach. CARDIOSIM© fulfils these requirements and does address the systems interaction with its modular approach and assembly of models with varying degree of complexity although 0-D and 1-D coupling may be required for the evaluation of long-term VAD support. The challenge remains the ability to predict outcome over a longer period of time. The preliminary CFD simulations with the HeartWare HVAD centrifugal pump demonstrated that it is possible to obtain an accurate analysis in a timely manner to complement the clinical review process. The simulations with the pilot patient-specific model of the HeartMate II axial flow pump revealed that a complex 3D reconstruction is feasible in a timely manner and can be used to generate sufficiently accurate results to be used in the context of a MDT meeting for the purposes of clinical decision-making. Overall, these three studies demonstrate that the time frame of the simulations was within hours which may fit the time constraints of the clinical environment in the context of a MDT meeting. More specifically, it was shown that the laminar model may be used for an initial evaluation of the flow development within the pump. Nonetheless, the k- model offers higher accuracy if the timeline of the clinical setting allows for a longer simulation. Conclusion: This thesis aimed at the understanding of the use of computational modelling as a pre-operative planning strategy and follow up assessment for the treatment of advanced heart failure with rotary blood pumps. The novelty lays in the use of both pressure-volume simulation analysis and 3D flow dynamics studies in VADs with a view to treatment optimisation and outcome prediction within the time constraints of a clinical setting in the context of a MDT meeting. The clinical significance and the contribution to the field is a more targeted approach for different groups of patients and a more quantitative evaluation in the clinical decision process based on a pro-active co-operation between clinicians and scientists reducing the potential for “guess work”. The results of this thesis are a proof-of-concept as a prelude to a potential future implementation of patient-specific modelling within a clinical setting on a daily basis demonstrating a clear clinical significance and contribution to the field. The proposed approach does not consider modelling and simulation as a substitute for clinical experience but an additional tool to guide therapeutic intervention and complement the clinical decision process in which the clinician remains the ultimate decision-maker. Such an approach may well add a different dimension to the problem of heart failure with potential for high return in terms of patient’s outcome and long-term surveillance. The same principles would be applicable to other cardiovascular problems in line with the current concept of “Team Approach” such as the Heart Team, the Structural Heart Team or the Aortic Team. The present work has taken this concept closer to clinical delivery and has highlighted its potential but further work remains to be done in refining the technique

Development and characterization of the arterial windkessel and its role during left ventricular assist device assistance

Modeling of the cardiovascular system is challenging, but it has the potential to further advance our understanding of normal and pathological conditions. Morphology and function are closely related. The arterial system provides steady blood flow to each organ and damps out wave fluctuations as a consequence of intermittent ventricular ejection. These actions can be approached separately in terms of mathematical relationships between pressure and flow. Lumped parameter models are helpful for the study of the interactions between the heart and the arterial system. The arterial windkessel model still plays a significant role despite its limitations. This review aims to discuss the model and its modifications and derive the fundamental equations by applying electric circuits theory. In addition, its role during left ventricular assist device assistance is explored and discussed in relation to rotary blood pumps

A Multi-Scale CFD Analysis of Patient-Specific Geometries to Tailor LVAD Cannula Implantation Under Pulsatile Flow Conditions: an investigation aimed at reducing stroke incidence in LVADs

A Left Ventricular Assist Device (LVAD) is a mechanical pump that provides temporary circulatory support when used as bridge-to-transplantation and relieves workload demand placed on a failing heart allowing for myocardia recovery when used as destination therapy. Stroke is the most devastating complication after ventricular assist device (VAD) implantation, with an incidence of 14-47% over 3-6 months. This complication due to thrombus formation and subsequent transport through the vasculature to cerebral vessels continues to limit the widespread implementation of VAD therapy. Patient-specific computational fluid dynamics (CFD) analysis may elucidate ways to reduce this risk. We employed a multi-scale model of the aortic circulation in order to examine the effects on flow conditions resulting from varying the VAD cannula implantation location and angle of incidence of the anastomosis to the ascending aorta based on a patient-specific geometry obtained from CT scans. The multi-scale computation consists of a 0D lumped parameter model (LPM) of the circulation modeled via a 50 degree of freedom (DOF) electrical circuit analogy that includes an LVAD model coupled to a 3D computational fluid dynamics model of the circulation. An in-house adaptive Runge-Kutta method is utilized to solve the 50 DOF LPM, and the Starccm+ CFD code is utilized to solve the flowfield. This 0D-3D coupling for the flow is accomplished iteratively with the 0D LPM providing the pulsatile boundary conditions that drive the 3D CFD time-accurate computations of the flowfield. Investigated angle configurations include cannula implantations at 30°, 60° and 90° to the right lateral wall of the ascending aorta. We also considered placements of the VAD cannula along the ascending aorta in which distances of the VAD anastomosis is varied relative to the take-off of the innominate artery. We implemented a mixed Eulerian-Lagrangian particle-tracking scheme to quantify the number of stroke-inducing particles reaching cerebral vessel outlets and included flow visualization through streamlines to identify regions of strong vorticity and flow stagnation, which can promote thrombus formation. Thrombi were modeled as spheres with perfectly elastic interactions numerically released randomly in time and space at cannula inlet plane. Based on clinical observation of the range of thrombus sizes encountered in such cases, particle diameters of 2.5mm and 3.5mm were investigated in our numerical computations. Pulsatile flow results for aforementioned angles suggest that a 90° cannula implementation causes flow impingement on the left lateral aortic wall and appears to be highly thrombogenic due to large momentum losses and zones of large re-circulation and that shallow and intermediate cannula angles promote more regular flow carrying particles towards the lower body potentially reducing stroke risk. Indications from this pulsatile numerical study suggest that up to a 50% reduction in stroke rate can be achieve with tailoring of cannula implantation. Results are consistent with significant reduction in stroke incidence achieved by tailoring cannula implantation as reported in previous steady flow computations carried out by our group. As such, results of this study suggest that a simple surgical maneuver in the process of VAD implantation may significantly improve patient life

Review of Zero-D and 1-D Models of Blood Flow in the Cardiovascular System

<p>Abstract</p> <p>Background</p> <p>Zero-dimensional (lumped parameter) and one dimensional models, based on simplified representations of the components of the cardiovascular system, can contribute strongly to our understanding of circulatory physiology. Zero-D models provide a concise way to evaluate the haemodynamic interactions among the cardiovascular organs, whilst one-D (distributed parameter) models add the facility to represent efficiently the effects of pulse wave transmission in the arterial network at greatly reduced computational expense compared to higher dimensional computational fluid dynamics studies. There is extensive literature on both types of models.</p> <p>Method and Results</p> <p>The purpose of this review article is to summarise published 0D and 1D models of the cardiovascular system, to explore their limitations and range of application, and to provide an indication of the physiological phenomena that can be included in these representations. The review on 0D models collects together in one place a description of the range of models that have been used to describe the various characteristics of cardiovascular response, together with the factors that influence it. Such models generally feature the major components of the system, such as the heart, the heart valves and the vasculature. The models are categorised in terms of the features of the system that they are able to represent, their complexity and range of application: representations of effects including pressure-dependent vessel properties, interaction between the heart chambers, neuro-regulation and auto-regulation are explored. The examination on 1D models covers various methods for the assembly, discretisation and solution of the governing equations, in conjunction with a report of the definition and treatment of boundary conditions. Increasingly, 0D and 1D models are used in multi-scale models, in which their primary role is to provide boundary conditions for sophisticate, and often patient-specific, 2D and 3D models, and this application is also addressed. As an example of 0D cardiovascular modelling, a small selection of simple models have been represented in the CellML mark-up language and uploaded to the CellML model repository <url>http://models.cellml.org/</url>. They are freely available to the research and education communities.</p> <p>Conclusion</p> <p>Each published cardiovascular model has merit for particular applications. This review categorises 0D and 1D models, highlights their advantages and disadvantages, and thus provides guidance on the selection of models to assist various cardiovascular modelling studies. It also identifies directions for further development, as well as current challenges in the wider use of these models including service to represent boundary conditions for local 3D models and translation to clinical application.</p

Computer simulation of platelet activation in a pulsatile ventricular assist device, through finite elements and a simplified geometry

Ventricular assist devices are a technological solution for patients who suffer from cardiac insufficiency and await a transplant. In this work, a new design of an implantable pulsatile blood pump is analyzed in terms of blood damage, by means of finite elements on a simplified geometry. It is a double effect volumetric pump, which has a noncontact driven piston and four active valves. The analysis is done by means of blood flow simulation into the pump and the prediction of platelets activation. The last is a measure of the pump compatibility with human life.The platelet activation state is evaluated by an equivalent or representative shear stress and compared with bibliographic data corresponding to other VAD kinds and cardiac prosthetic valves. The results show that, for the complementary blood flow rate supplied by the simulated VAD, the predicted platelet damage is in the same levels of current cardiac devices, particularly other VADs.Fil: Fries, Exequiel. Universidad Nacional de Entre Ríos. Facultad de Ingeniería. Departamento de Física Química. Laboratorio de Biomecánica Computacional; ArgentinaFil: Berli, Marcelo Eduardo. Universidad Nacional de Entre Ríos. Facultad de Ingeniería. Departamento de Física Química. Laboratorio de Biomecánica Computacional; ArgentinaFil: Campana, Diego Martin. Universidad Nacional de Entre Ríos. Facultad de Ingeniería. Departamento de Física Química. Laboratorio de Biomecánica Computacional; ArgentinaFil: Ubal, Sebastian. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Santa Fe. Instituto de Desarrollo Tecnológico para la Industria Química. Universidad Nacional del Litoral. Instituto de Desarrollo Tecnológico para la Industria Química; ArgentinaFil: Di Paolo, José. Universidad Nacional de Entre Ríos. Facultad de Ingeniería. Departamento de Física Química. Laboratorio de Biomecánica Computacional; Argentin