A modular computational circulatory model applicable to VAD testing and training

The hybrid (physical-computational) cardiovascular simulator to study valvular diseases

To better understand the impact of valvular heart disease (VHD) on the hemodynamics of the circulatory system, investigations can be carried out using a model of the cardiovascular system. In this study, a previously developed hybrid (hydro-numerical) simulator of the cardiovascular system (HCS) was adapted and used. In our HCS Björk-Shiley mechanical heart valves were used, playing the role of mitral and aortic ones. In order to simulate aortic stenosis (AS) and mitral regurgitation (MR), special mechanical devices have been developed and integrated with the HCS. The simulation results proved that the system works correctly. Namely, in the case of AS - the mean pulmonary arterial pressure was increased due to increased preload of the left ventricle and the decrease in right ventricular preload was caused by a decrease in systemic arterial pressure. The severity of AS was performed based on the transaortic pressure gradient as well as using the Gorlin and Aaslid equations. In the case of severe AS, when the mean gradient was above 40 mmHg, the aortic valve orifice area was 0.5 cm2, which is in line with ACC/AHA guidelines. For the case of MR - with increasing severity of MR, there was a decrease in the left ventricular pressure and an increase in left atrial pressure. Using mechanical heart valves to simulate VHD by the HCS can be a valuable tool for biomedical research, providing a safe and controlled environment to study and understand the pathophysiology of VHD.

Experimental and numerical investigation of the stenosed coronary artery taken from the clinical setting and modeled in terms of hemodynamics

The study was carried out to investigate the effect of the artery with different pulse values and stenosis rates on the pressure drop, the peristaltic pump outlet pressure, fractional flow reserve (FFR) and most importantly the amount of power consumed by the peristaltic pump. For this purpose, images taken from the clinical environment were produced as models (10 mm inlet diameter) with 0% and 70% percent areal stenosis rates (PSR) on a three-dimensional (3D) printer. In the experimental system, pure water was used as the fluid at 54, 84, 114, 132, and 168 bpm pulse values. In addition, computational fluid dynamics (CFD) analyzes of the test region were performed using experimental boundary conditions with the help of ANSYS-Fluent software. The findings showed that as PSR increases in the arteries, the pressure drop in the stenosis region increases and this amount increases dramatically with increasing effort. An increase of approximately 40% was observed in the pump outlet pressure value from 54 bpm to 168 bpm in the PSR 0% model and 51% increase in the PSR 70% model. It has been observed that the pump does more work to overcome the increased pressure difference due to increased pulse rate and PSR. With the effect of contraction, the power consumption of the pump increased from 9.2% for 54 bpm to 13.8% for 168 bpm. In both models, the Wall Shear Stress (WSS) increased significantly. WSS increased abruptly in the stenosis and arcuate regions, while sudden decreases were observed in the flow separation region.

In Vivo Material Properties of Human Common Carotid Arteries: Trends and Sex Differences

INTRODUCTION

In vivo estimation of material properties of arterial tissue can provide essential insights into the development and progression of cardiovascular diseases. Furthermore, these properties can be used as an input to finite element simulations of potential medical treatments.

MATERIALS AND METHODS

This study uses non-invasively measured pressure, diameter and wall thickness of human common carotid arteries (CCAs) acquired in 103 healthy subjects. A non-linear optimization was performed to estimate material parameters of two different constitutive models: a phenomenological, isotropic model and a structural, anisotropic model. The effect of age, sex, body mass index and blood pressure on the parameters was investigated.

RESULTS AND CONCLUSION

Although both material models were able to model in vivo arterial behaviour, the structural model provided more realistic results in the supra-physiological domain. The phenomenological model predicted very high deformations for pressures above the systolic level. However, the phenomenological model has fewer parameters that were shown to be more robust. This is an advantage when only the physiological domain is of interest. The effect of stiffening with age, BMI and blood pressure was present for women, but not always for men. In general, sex had the biggest effect on the mechanical properties of CCAs. Stiffening trends with age, BMI and blood pressure were present but not very strong. The intersubject variability was high. Therefore, it can be concluded that finding a representative set of parameters for a certain age or BMI group would be very challenging. Instead, for purposes of patient-specific modelling of surgical procedures, we currently advise the use of patient-specific parameters.

Virtual and Artificial Cardiorespiratory Patients in Medicine and Biomedical Engineering

Recently, ‘medicine in silico’ has been strongly encouraged due to ethical and legal limitations related to animal experiments and investigations conducted on patients. Computer models, particularly the very complex ones (virtual patients—VP), can be used in medical education and biomedical research as well as in clinical applications. Simpler patient-specific models may aid medical procedures. However, computer models are unfit for medical devices testing. Hybrid (i.e., numerical–physical) models do not have this disadvantage. In this review, the chosen approach to the cardiovascular system and/or respiratory system modeling was discussed with particular emphasis given to the hybrid cardiopulmonary simulator (the artificial patient), that was elaborated by the authors. The VP is useful in the education of forced spirometry, investigations of cardiopulmonary interactions (including gas exchange) and its influence on pulmonary resistance during artificial ventilation, and explanation of phenomena observed during thoracentesis. The artificial patient is useful, inter alia, in staff training and education, investigations of cardiorespiratory support and the testing of several medical devices, such as ventricular assist devices and a membrane-based artificial heart.

The use of a virtual patient to follow changes in arterial blood gases associated with therapeutic thoracentesis

PURPOSES:

Some controversies exist on the effect of therapeutic thoracentesis (TT) on arterial blood oxygen tension. The aim of this study was to evaluate this issue using a previously developed virtual patient.

METHODS:

The analysis was based and supported by clinical data collected during 36 TT. Pleural pressure and transcutaneous oxygen and carbon dioxide pressures (PtcO₂ and PtcCO₂) were measured during pleural fluid withdrawal. Arterial blood oxygen tension and arterial CO₂ tension (PaO₂ and PaCO₂) were analysed in simulations that mimicked TT. Minute ventilation was adjusted to maintain arterial CO₂ tension at a constant level unless arterial blood oxygen tension fell below 8 kPa. Specifically, the influence of hypoxic pulmonary vasoconstriction efficiency was tested.

RESULTS:

In patients, PtcCO₂ remained at an approximately constant level (average amplitude: 0.63 ± 0.29 kPa), while some fluctuations of PtcO2 were observed (amplitude: (1.65 ± 1.18 kPa) were observed. In 42% of patients, TT was associated with decrease in PtcCO₂. Simulations showed the following: (a) there were similar PaO₂ fluctuations in the virtual patient; (b) the lower the hypoxic pulmonary vasoconstriction efficiency, the more pronounced the PaO₂ fall during fluid withdrawal; and (c) the lower the atelectatic lung areas recruitment rate, the slower the PaO₂ normalization. The decrease in PaO₂ was caused by an increase of pulmonary shunt.

CONCLUSION:

Therapeutic thoracentesis may cause both an increase and a decrease in PaO₂ during the procedure. Pleural pressure decrease, caused by pleural fluid withdrawal, improves the perfusion of atelectatic lung areas. If the rate of recruitment of these areas is low, a lack of ventilation causes the arterial blood oxygen tension to fall. Effective hypoxic pulmonary vasoconstriction may protect against the pulmonary shunt.

A mock heart engineered with helical aramid fibers for in vitro cardiovascular device testing

Mock heart circulation loops (MHCLs) serve as in-vitro platforms to investigate the physiological interaction between circulatory systems and cardiovascular devices. A mock heart (MH) engineered with silicone walls and helical aramid fibers, to mimic the complex contraction of a natural heart, has been developed to advance the MHCL previously developed in our group. A mock aorta with an anatomical shape enables the evaluation of a cannulation method for ventricular assist devices (VADs) and investigation of the usage of clinical measurement systems like pressure-volume catheters. Ventricle and aorta molds were produced based on MRI data and cast with silicone. Aramid fibers were layered in the silicone ventricle to reproduce ventricle torsion. A rotating hollow shaft was connected to the apex enabling the rotation of the MH and the connection of a VAD. Silicone wall thickness, aramid fiber angle and fiber pitch were varied to generate different MH models. All MH models were placed in a tank filled with variable amounts of water and air simulating the compliance. In this work, physiological ventricular torsion angles (15°-26°) and physiological pressure-volume loops were achieved. This MHCL can serve as a comprehensive testing platform for cardiovascular devices, such as artificial heart valves and cannulation of VADs.

Hemodynamic analysis of outflow grafting positions of a ventricular assist device using closed-loop multiscale CFD simulations: Preliminary results

Subclavian arteries are a possible alternate location for left ventricular assist device (LVAD) outflow grafts due to easier surgical access and application in high risk patients. As vascular blood flow mechanics strongly influence the clinical outcome, insights into the hemodynamics during LVAD support can be used to evaluate different grafting locations. In this study, the feasibility of left and right subclavian artery (SA) grafting was investigated for the HeartWare HVAD with a numerical multiscale model. A 3-D CFD model of the aortic arch was coupled to a lumped parameter model of the cardiovascular system under LVAD support. Grafts in the left and right SA were placed at three different anastomoses angles (90°, 60° and 30°). Additionally, standard grafting of the ascending and descending aorta was modelled. Full support LVAD (5l/min) and partial support LVAD (3l/min) in co-pulsation and counter-pulsation mode were analysed. The grafting positions were investigated regarding coronary and cerebral perfusion. Furthermore, the influence of the anastomosis angle on wall shear stress (WSS) was evaluated. Grafting of left or right subclavian arteries has similar hemodynamic performance in comparison to standard cannula positions. Angularity change of the graft anastomosis from 90° to 30° slightly increases the coronary and cerebral blood flow by 6-9% while significantly reduces the WSS by 35%. Cannulation of the SA is a feasible anastomosis location for the HVAD in the investigated vessel geometry.

Hemodynamic Effects of Ventricular Assist Device Implantation on Norwood, Glenn, and Fontan Circulation: A Simulation Study

The growing population of failing single-ventricle (SV) patients might benefit from ventricular assist device (VAD) support as a bridge to heart transplantation. However, the documented experience is limited to isolated case reports. Considering the complex and different physiopathology of Norwood, Glenn, and Fontan patients and the lack of established experience, the aim of this work is to realize and test a lumped parameter model of the cardiovascular system able to simulate SV hemodynamics and VAD implantation effects to support clinical decision. Hemodynamic and echocardiographic data of 30 SV patients (10 Norwood, 10 Glenn, and 10 Fontan) were retrospectively collected and used to simulate patients' baseline. Then, the effects of VAD implantation were simulated. Simulation results suggest that the implantation of VAD: (i) increases the cardiac output and the mean arterial systemic pressure in all the three palliation conditions (Norwood 77.2 and 19.7%, Glenn 38.6 and 32.2%, and Fontan 17.2 and 14.2%); (ii) decreases the SV external work (Norwood 55%, Glenn 35.6%, and Fontan 41%); (iii) decreases the pressure pulsatility index (Norwood 65.2%, Glenn 81.3%, and Fontan 64.8%); (iv) increases the pulmonary arterial pressure in particular in the Norwood circulation (Norwood 39.7%, Glenn 12.1% and Fontan 3%); and (v) decreases the atrial pressure (Norwood 2%, Glenn 10.6%, and Fontan 8.6%). Finally, the VAD work is lower in the Norwood circulation (30.4 mL·mm Hg) in comparison with Fontan (40.3 mL·mm Hg) and to Glenn (64.5 mL·mm Hg) circulations. The use of VAD in SV physiology could be helpful to bridge patients to heart transplantations by increasing the CO and unloading the SV with a decrement of the atrial pressure and the SV external work. The regulation of the pulmonary flow is challenging because the Pap is increased by the presence of VAD. The hemodynamic changes are different in the different SV palliation step. The use of numerical models could be helpful to support patient and VAD selection to optimize the clinical outcome.

Simulation of Ventricular, Cavo-Pulmonary, and Biventricular Ventricular Assist Devices in Failing Fontan

Considering the lack of donors, ventricular assist devices (VADs) could be an alternative to heart transplantation for failing Fontan patients, in spite of the lack of experience and the complex anatomy and physiopathology of these patients. Considering the high number of variables that play an important role such as type of Fontan failure, type of VAD connection, and setting (right VAD [RVAD], left VAD [LVAD], or biventricular VAD [BIVAD]), a numerical model could be useful to support clinical decisions. The aim of this article is to develop and test a lumped parameter model of the cardiovascular system simulating and comparing the VAD effects on failing Fontan. Hemodynamic and echocardiographic data of 10 Fontan patients were used to simulate the baseline patients' condition using a dedicated lumped parameter model. Starting from the simulated baseline and for each patient, a systolic dysfunction, a diastolic dysfunction, and an increment of the pulmonary vascular resistance were simulated. Then, for each patient and for each pathology, the RVAD, LVAD, and BIVAD implantations were simulated. The model can reproduce patients' baseline well. In the case of systolic dysfunction, the LVAD unloads the single ventricle and increases the cardiac output (CO) (35%) and the arterial systemic pressure (Pas) (25%). With RVAD, a decrement of inferior vena cava pressure (Pvci) (39%) was observed with 34% increment of CO, but an increment of the single ventricle external work (SVEW). With the BIVAD, an increment of Pas (29%) and CO (37%) was observed. In the case of diastolic dysfunction, the LVAD increases CO (42%) and the RVAD decreases the Pvci, while both increase the SVEW. In the case of pulmonary vascular resistance increment, the highest CO (50%) and Pas (28%) increment is obtained with an RVAD with the highest decrement of Pvci (53%) and an increment of the SVEW but with the lowest VAD power consumption. The use of numerical models could be helpful in this innovative field to evaluate the effect of VAD implantation on Fontan patients to support patient and VAD type selection personalizing the assistance.

A Modeling Tool to Study the Combined Effects of Drug Administration and Lvad Assistance in Pathophysiological Circulatory Conditions

The aim of this work is to develop a tool to study the effect of sodium nitroprusside (SNP) on hemodynamics in conjunction with baroreflex and mechanical circulatory assistance. To this aim, a numerical model of the pharmacodynamic effect of SNP was developed and inserted into a cardiovascular circulatory model integrated with baroreflex and LVAD (continuous flow pump with atrio-aortic connection) sub-models.

The experiments were carried out in two steps. In the first step the model was verified comparing simulations with experimental data acquired from mongrel dogs on mean arterial pressure (MAP), cardiac output (CO), heart rate (HR), peripheral resistance, and left ventricular properties.

In the second step, the combined action of SNP and mechanical circulatory assistance was studied. Data were measured at pump off and at pump on (20000 rpm and 24000 rpm). At pump off, with a 2.5 μg/kg per min SNP infusion in heart failure condition, the MAP was reduced by approximately 8%, CO and HR increased by about 16% and 18%, respectively. In contrast, during assistance (24000 rpm) the changes in MAP, CO and HR were around −9%, +12%, and +20%, respectively. Furthermore, the effects of the drug on hemodynamic parameters at different heart conditions were significantly different. Thus, the model provides insight into the complex interactions between baroreflex, drug infusion, and LVAD and could be a support for clinical decision-making in cardiovascular pathologies.

Simulation of Apical and Atrio-aortic VAD in Patients with Transposition or Congenitally Corrected Transposition of the Great Arteries

Purpose

VADs could be used for transportation of the great arteries (TGA) and for congenitally corrected transposition (ccTGA) treatment. A cardiovascular numerical model (NM) may offer a useful clinical support in these complex physiopathologies. This work aims at developing and preliminarily verifying a NM of ccTGA and TGA interacting with VADs.

Methods

Hemodynamic data were collected at the baseline (BL) and three months (FUP) after apical (atrio-aortic) VAD implantation in a TGA (ccTGA) patient and used in a lumped parameter NM to simulate the patient's physiopathology. Measured (MS) and simulated (SIM) data were compared.

Results

MS and SIM data are in accordance at the BL and at FUP. Cardiac output (l/min): BL_m = 2.9 ± 0.4, BL_s = 3.0 ± 0.3; FUP_m = 4.2 ± 0.2, FUP_s = 4.1 ± 0.1. Right atrial pressure (mmHg): BL_m = 21.4 ± 4.1, BL_s = 18.5 ± 4.5; FUP_m = 13 ± 4, FUP_s = 14.8 ± 3.6. Pulmonary arterial pressure (mmHg): BL_m = 56 ± 6.3, BL_s = 57 ± 2, FUP_m = 37.5 ± 7.5, FUP_s = 35.5 ± 5.9. Systemic arterial pressure (mmHg): BL_m = 71 ± 2, BL_s = 74.6 ± 2.1; FUP_m = 84 ± 9, FUP_s = 81.9 ± 9.8.

Conclusions

NM can simulate the effect of a VAD in complex physiopathologies, with the inclusion of changes in circulatory parameters during the acute phase and at FUP. The simulation of differently assisted physiopathologies offers a useful support for clinicians.

Continuous-flow pump model study: the effect on pump performance of pump characteristics and cardiovascular conditions

This model study evaluates the effect of pump characteristics and cardiovascular data on hemodynamics in atrio-aortic VAD assistance. The model includes a computational circulatory sub-model and an electrical sub-model representing two rotary blood pumps through their pressure-flow characteristics. The first is close to a pressure generator-PG (average flow sensitivity to pressure variations, -0.047 l mmHg(-1)); the second is closer to a flow generator-FG (average flow sensitivity to pressure variations, -0.0097 l mmHg(-1)). Interaction with VAD was achieved by means of two interfaces, behaving as impedance transformers. The model was verified by use of literature data and VAD onset conditions were used as a control for the experiments. Tests compared the two pumps, at constant pump speed, in different ventricular and circulatory conditions: maximum ventricular elastance (0.44-0.9 mmHg cm(-3)), systemic peripheral resistance (781-1200 g cm(-4) s(-1)), ventricular diastolic compliance C p (5-10-50 cm(3) mmHg(-1)), systemic arterial compliance (0.9-1.8 cm(3) mmHg(-1)). Analyzed variables were: arterial and venous pressures, flows, ventricular volume, external work, and surplus hemodynamic energy (SHE). The PG pump generated the highest SHE under almost all conditions, in particular for higher C p (+50 %). PG pump flow is also the most sensitive to E max and C p changes (-26 and -33 %, respectively). The FG pump generally guarantees higher external work reduction (54 %) and flow less dependent on circulatory and ventricular conditions. The results are evidence of the importance of pump speed regulation with changing ventricular conditions. The computational sub-model will be part of a hydro-numerical model, including autonomic controls, designed to test different VADs.