Computer simulation of haemodynamic parameters changes with left ventricle assist device and mechanical ventilation

MATHEMATICAL MODELS OF HUMAN RESPIRATORY AND BLOOD CIRCULATORY SYSTEMS

Aim. To analyze modern approaches to mathematical modeling of human respiratory and blood circulatory systems. Methods. Comprehensive review of scientific literature sources extracted from domestic and international resources databases. Results. Historical information and modern data concerning mathematical modeling of human functional respiratory and blood circulatory systems were summarized and analyzed in present ¬review; current trends in approaches to the construction of these models were revealed. Conclusions. Currently, two main approaches to the mathematical modeling of respiratory and blood circulatory systems exist. One of them is the construction of models of the mechanics of respiration and blood circulation. They are based on the models of mechanics of solid deformable body, thermomechanics, hydromechanics, and continuum mechanics. This approach uses complex mathematical apparatus, including Navier-Stokes equation, which makes it possible to obtain a number of theoretical results, but it is hardly usable for real problems solutions at present time. The second approach is based on the model of F. Grodins, who represented the process of breathing as a controlled dynamic system, described by ordinary differential equations, in which the process control is carried out according to the feedback principle. There is a significant number of modifications of this model, which made it possible to simulate various disturbing influences, such as physical activity, hypoxia and hyperemia, and to predict parameters characterizing functional respiratory system under these disturbing influences.

Application of a user-friendly comprehensive circulatory model for estimation of hemodynamic and ventricular variables

Purpose

Application of a comprehensive, user-friendly, digital computer circulatory model to estimate hemodynamic and ventricular variables.

Methods

The closed-loop lumped parameter circulatory model represents the circulation at the level of large vessels. A variable elastance model reproduces ventricular ejection. The circulatory model has been modified embedding an algorithm able to adjust the model parameters reproducing specific circulatory conditions. The algorithm reads input variables: heart rate, aortic pressure, cardiac output, and left atrial pressure. After a preliminary estimate of circulatory parameters and ventricular elastance, it adjusts the amount of circulating blood, the value of the systemic peripheral resistance, left ventricular elastance, and ventricular rest volume. Input variables and the corresponding calculated variables are recursively compared: the procedure is stopped if the difference between input and calculated variables is within the set tolerance. At the procedure end, the model produces an estimate of ventricular volumes and Emaxl along with systemic and pulmonary pressures (output variables). The procedure has been tested using 4 sets of experimental data including left ventricular assist device assistance.

Results

The algorithm allows the reproduction of the circulatory conditions defined by all input variable sets, giving as well an estimate of output variables.

Conclusions

The algorithm permits application of the model in environments where the simplicity of use and velocity of execution are of primary importance. Due to its modular structure, the model can be modified adding new circulatory districts or changing the existing ones. The model could also be applied in educational applications.

Preservation of native aortic valve flow and full hemodynamic support with the TORVAD using a computational model of the cardiovascular system

This article describes the stroke volume selection and operational design for the toroidal ventricular assist device (TORVAD), a synchronous, positive-displacement ventricular assist device (VAD). A lumped parameter model was used to simulate hemodynamics with the TORVAD compared with those under continuous-flow VAD support. Results from the simulation demonstrated that a TORVAD with a 30 ml stroke volume ejecting with an early diastolic counterpulse provides comparable systemic support to the HeartMate II (HMII) (cardiac output 5.7 L/min up from 3.1 L/min in simulated heart failure). By taking the advantage of synchronous pulsatility, the TORVAD delivers full hemodynamic support with nearly half the VAD flow rate (2.7 L/min compared with 5.3 L/min for the HMII) by allowing the left ventricle to eject during systole and thus preserving native aortic valve flow (3.0 L/min compared with 0.4 L/min for the HMII, down from 3.1 L/min at baseline). The TORVAD also preserves pulse pressure (26.7 mm Hg compared with 12.8 mm Hg for the HMII, down from 29.1 mm Hg at baseline). Preservation of aortic valve flow with synchronous pulsatile support could reduce the high incidence of aortic insufficiency and valve cusp fusion reported in patients supported with continuous-flow VADs.

Improving Arterial Pulsatility by Feedback Control of a Continuous Flow Left Ventricular Assist Device via in Silico Modeling

Purpose

Continuous flow left ventricular assist devices (CF-LVADs) generally operate at a constant speed, which causes a decrease in pulse pressure and pulsatility in the arteries and allegedly may lead to late complications such as aortic insufficiency and gastrointestinal bleeding. The purpose of this study is to increase the arterial pulse pressure and pulsatility while obtaining more physiological hemodynamic signals, by controlling the CF-LVAD flow rate.

Methods

A lumped parameter model was used to simulate the cardiovascular system including the heart chambers, heart valves, systemic and pulmonary arteries and veins. A baroreflex model was used to regulate the heart rate and a model of the Micromed DeBakey CF-LVAD (Micromed Technology, Houston, TX, USA) to simulate the pump dynamics at different operating speeds. A model simulating the flow rate through the aortic valve served as reference model. CF-LVAD operating speed was regulated by applying proportional-integral (PI) control to the pump flow rate. For comparison, the CF-LVAD was also operated at a constant speed, equaling the mean CF-LVAD speed as applied in pulsatile mode.

Results

In different operating modes, at the same mean operating speeds, mean pump output, mean arterial pressure, end-systolic and end-diastolic volume and heart rate were the same over the cardiac cycle. However, the arterial pulse pressure and index of pulsatility increased by 50% in the pulsatile CF-LVAD support mode with respect to constant speed pump support.

Conclusions

This study shows the possibility of obtaining more physiological pulsatile hemodynamics in the arteries by applying output-driven varying speed control to a CF-LVAD.

Evaluating the hemodynamical response of a cardiovascular system under support of a continuous flow left ventricular assist device via numerical modeling and simulations

Dilated cardiomyopathy is the most common type of the heart failure which can be characterized by impaired ventricular contractility. Mechanical circulatory support devices were introduced into practice for the heart failure patients to bridge the time between the decision to transplant and the actual transplantation which is not sufficient due to the state of donor organ supply. In this study, the hemodynamic response of a cardiovascular system that includes a dilated cardiomyopathic heart under support of a newly developed continuous flow left ventricular assist device—Heart Turcica Axial—was evaluated employing computer simulations. For the evaluation, a numerical model which describes the pressure-flow rate relations of Heart Turcica Axial, a cardiovascular system model describing the healthy and pathological hemodynamics, and a baroreflex model regulating the heart rate were used. Heart Turcica Axial was operated between 8000 rpm and 11000 rpm speeds with 1000 rpm increments for assessing the pump performance and response of the cardiovascular system. The results also give an insight about the range of the possible operating speeds of Heart Turcica Axial in a clinical application. Based on the findings, operating speed of Heart Turcica Axial should be between 10000 rpm and 11000 rpm.

Verification of a computational cardiovascular system model comparing the hemodynamics of a continuous flow to a synchronous valveless pulsatile flow left ventricular assist device

The purpose of this investigation is to use a computational model to compare a synchronized valveless pulsatile left ventricular assist device with continuous flow left ventricular assist devices at the same level of device flow, and to verify the model with in vivo porcine data. A dynamic system model of the human cardiovascular system was developed to simulate the support of a healthy or failing native heart from a continuous flow left ventricular assist device or a synchronous pulsatile valveless dual-piston positive displacement pump. These results were compared with measurements made during in vivo porcine experiments. Results from the simulation model and from the in vivo counterpart show that the pulsatile pump provides higher cardiac output, left ventricular unloading, cardiac pulsatility, and aortic valve flow as compared with the continuous flow model at the same level of support. The dynamic system model developed for this investigation can effectively simulate human cardiovascular support by a synchronous pulsatile or continuous flow ventricular assist device.

Study of systolic pressure variation (SPV) in presence of mechanical ventilation

Systolic pressure variation (SPV) and its components (dUp and dDown) have been demonstrated to be of interest in assessing preload in mechanically ventilated patients. The aim of this paper is to analyse the sensitivity of these variables to preload and volemic changes during mechanical ventilation in different conditions of the cardiovascular system. Computer simulation experiments have been done using a modular lumped parameter model of the cardiovascular system. The effect of mechanical ventilation has been reproduced operating on intrathoracic pressure. Experiments have been performed varying preload through filling pressure. Sensitivity of SVP dUp and dDown is described varying separately left ventricular elastance (Ev), systemic arterial resistance (Ras) and systemic arterial compliance (Cas). The sensitivity of SPV and dDown to preload and filling pressure is appreciable for high values of Ev and for a wide variation of Ras. Preliminary clinical data concerning the three parameters show good correlation with simulation results.

A virtual instrument (VI) for haemodynamic management in ICU and during surgery

Systolic pressure variations (SPV) during mechanical ventilation and its single components, related to short apnea, reflect changes of the volemic condition of the patient. To introduce their determination during clinical monitoring for different fluid states and for different tidal volumes, they must be computed on-line without introducing interference with standard activities. A system computing on-line systolic pressure variation during mechanical ventilation, connected to standard monitoring devices, has been proposed. It is based on a notebook PC implemented with graphical software comprising a user panel in the form of a virtual instrument and is able to acquire, process and present signals from different instruments utilized in ICU and during surgery. It can be used as a base to assess the ability of computed parameters helpful in clinical decision. The use of a notebook PC and open software allows operators, even if non-expert in computer science, to test and implement this, as well as other innovative tools in clinical practice.