- Research
- Open Access
Windkessel model of hemodynamic state supported by a pulsatile ventricular assist device in premature ventricle contraction
- Keun Her1,
- Joon Yeong Kim2,
- Ki Moo Lim3 and
- Seong Wook Choi2Email authorView ORCID ID profile
- Received: 21 July 2017
- Accepted: 10 January 2018
- Published: 2 February 2018
Abstract
Background
Counter-pulsation control (CPC) by ventricular assist devices (VADs) is believed to reduce cardiac load and increase coronary perfusion. However, patients with VADs have a higher risk of arrhythmia, which may cause the CPC to fail. Consequently, CPC has not been applied by VADs in clinical practice. The phase-locked loop (PLL) algorithm for CPC is readily implemented in VADs; however, it requires a normal, consistent heartbeat for adequate performance. When an arrhythmia occurs, the algorithm maintains a constant pumping rate despite the unstable heartbeat. Therefore, to apply the PLL algorithm to CPC, the hemodynamic effects of abnormal heartbeats must be analyzed.
Objectives
This study sought to predict the hemodynamic effects in patients undergoing CPC using VADs, based on electrocardiogram (ECG) data, including a wide range of heart rate (HR) changes caused by premature ventricular contraction (PVC) or other reasons.
Methods
A four-element Windkessel hemodynamic model was used to reproduce the patient’s aortic blood pressure in this study. ECG data from 15 patients with severe congestive heart failure were used to assess the effect of the CPC on the patients’ hemodynamic state. The input and output flow characteristics of the pulsatile VAD (LibraHeart I, Cervika, Korea) were measured using an ultrasound blood flow meter (TS410, Transonic, USA), with the aortic pressure maintained at 80–120 mmHg. All other patient conditions were also reproduced.
Results
In patients with PVCs or normal heartbeats, CPC controlled by a VAD reduced the cardiac load by 20 and 40%, respectively. When the HR was greater for other reasons, such as sinus tachycardia, simultaneous ejection from the heart and VAD was observed; however, the cardiac load was not increased by rapid cardiac contractions resulting from decreased left ventricle volume. These data suggest that the PLL algorithm reduces the cardiac load and maintains consistent hemodynamic changes.
Keywords
- Windkessel model
- Pulsatile ventricular assist device
- Counter-pulsation control
- Phase-locked loop
- Arrhythmia
Background
Use of a ventricular assist device (VAD) is the most effective way to improve the survival rates of patients with end-stage heart failure [1–5]. However, to avoid additional heart problems or thrombus, new approaches such as pulsatile flow and counter-pulsation control (CPC) are required [5, 6]. Although CPC is expected to increase cardiovascular circulation and reduce cardiac load, it has not been implemented clinically because of the possible adverse effects on the heart in the event of CPC failure [7–9]. Patients with VADs are at a higher risk of developing arrhythmia, which can induce CPC failure [10, 11]. Therefore, prior to engaging in CPC, it should be confirmed that failure of the VAD will not cause the patient’s hemodynamic state or cardiac load to worsen.
Multiple studies of CPC using VADs have assumed normal HRs without considering the possible occurrence of arrhythmias. Lim et al. reported that CPC increased coronary perfusion and decreased cardiac load [7]. Koshimoto et al. assessed such factors using pulsatile rotary VADs in animal studies [8, 9]. Some CPC algorithms have considered tachycardia and bradycardia; however, this approach has not been tested clinically or in animals due to the difficulties associated with inducing arrhythmias [12]. Clinical studies involving CPC are complicated despite the increasing use of VAD applications. This is because many reports have suggested that VADs should not be used for patients with arrhythmias [13–15]. Moreover, most VADs used to date did not have a CPC function [12, 16]. Intra-aortic balloon pumps (IABPs) have successfully engaged in CPC and beneficially affected heart treatment [17, 18]. However, in cases involving VADs, it is difficult to apply the CPC algorithm for IABPs. IABPs inject helium gas through a catheter into the balloon in the aorta rapidly, whereas the outflow pulse of a VAD as it passes through the blood stream is slow [9, 18, 19]. Therefore, the CPC algorithm of VADs should predict future heartbeats based on prior heartbeats, and preemptively provide subsequent pulses [12].
The LibraHeart I (Cervika, Korea) is a pulsatile VAD that analyzes prior heart beats and initiates its next pulse flow to keep the time interval between the pulses of the patient’s heart and the VAD constant according to the phase locked loop (PLL) algorithm [6]. The PLL is an algorithm that determines the pumping rate of a VAD whose output pulses are related to the QRS complex of the patient’s electrocardiogram (ECG) [20]. However, the PLL algorithm for a VAD can fail at CPC when arrhythmias occur. Therefore, the likelihood of CPC failure, and the effects it would have on cardiac load and the hemodynamic state, should be investigated before determining whether to implement CPC using VADs. This study analyzed ECG data from 15 patients with congestive heart failure and observed 35 arrhythmias, including premature ventricular contractions (PVCs) and sinus tachycardia (ST), in those ECGs [21].
The hemodynamic model was based on the four-element Windkessel model, which reproduces patient blood pressure and blood flow data [22, 23]. Because the conventional Windkessel model does not consider the reflection wave in the aorta, the actual blood pressure and blood flow data of a patient cannot be reproduced [24]; whereas the four-element model considers the reflection wave and reproduces the same aortic pressure and blood flow as measured in the patient [25]. The four-element model has been used to evaluate the accuracy of medical devices that analyze blood pressure and blood flow [22].
The purpose of this study was to predict the effects of a VAD in CPC mode on the cardiac and hemodynamic states using a four-element Windkessel model that includes the data for a VAD and arrhythmias that induced CPC failure.
Methods
Counter pulsation control by the ventricular assist device using the phase-locked loop method
The conventional method of determining when to initiate VADs following heartbeat detection misses the appropriate counter-pulsation point. If the delay from the time of heart rate (HR) detection to actuator operation is not reduced to within 0.1 s [9], a normal heartbeat of over 70 bpm can be missed. Following heartbeat detection, delay occurs due to the time required to establish the pulse flow by increasing the actuator speed, thereby generating the blood pressure pulse at the outlet of the VAD and delivering the pulse to the aortic valve through vessels and blood [18]. Comparatively, an IABP successfully applies CPC and reduces delays by injecting helium gas through a stiff, narrow catheter to a balloon located near the heart [19]. To overcome the slow response of a pulsatile VAD, the controller should predict the starting point by analyzing the regularity of the heartbeat and increasing motor speed prior to detecting the following beat [6, 12].
Phase-locked loop (PLL) algorithm for counter-pulsation control (CPC)
When an arrhythmia occurred, TR-R changed temporarily from 30 to 170%, which caused temporary CPC failure. In most cases, however, the VAD slowly adjusted its PR to the average HR according to the PLL algorithm, and it seemed to maintain a constant PR as opposed to changing the PR of the IABP for each heartbeat.
Hemodynamic model of the left ventricle and aorta
Parameters related to reflection wave in the 4-element Windkessel model
Symbol | Values | Units | Meaning of parameter |
---|---|---|---|
CAorta | 0.15 | cc/mmHg | Compliance of the aorta |
ZAorta | |||
LAorta | 0.0015 | cc/s2 | Impedance of the aorta composed of inertia and resistance of blood flow |
RAorta | 0.08 | mmHg s/cc | |
CPA | 1.5 | cc/mmHg | Compliance of the peripheral arteries |
RPA | 1.35 | mmHg s/cc | Resistance of the peripheral vessels |
Auxiliary parameters used to estimate the cardiac load supported by the VAD
Symbol | Values | Units | Meaning of parameter |
---|---|---|---|
CLV | 12.5 | cc/mmHg | Compliance of the LV during LV relaxation |
0.66–12.5 | Compliance of the LV during LV contraction | ||
RAV | 0.002 | mmHg s/cc | Resistance of the aortic valve when it opens |
1000 | Resistance of the aortic valve when it closes | ||
RMV | 0.01 | mmHg s/cc | Resistance of the mitral valve when it opens |
1000 | Resistance of the mitral valve when it closes | ||
PLA | 10 | mmHg | Blood pressure in the left atrium |
VAD | 2 (average) | L/min | Blood flow through the VAD |
a 4-element Windkessel model, including the ventricular assist device (VAD), b aortic blood pressure (AoP), and c pressure–volume curves before and during VAD perfusion
a Representative electrocardiogram (ECG) data showing temporary premature ventricle contraction (PVC) obtained from a congestive heart failure patient, and b the patient’s AoP, reproduced by the Windkessel model
ECG data were obtained from 15 patients (11 men, aged 22–71 years; 4 women, aged 54–63 years) with severe congestive heart failure (NYHA class 3–4) at a sampling rate of 250 per second. The data included 1260 heart beats and 28 arrhythmias [21]. Nine patients had tachycardia with a rapid HR over 100 bpm. The primary arrhythmias were PVCs (Fig. 3a), which animal studies suggest are more likely to occur when using a pulsatile VAD [6]. Although the contraction force of the left ventricle was set to be constant, the Windkessel model predicted that the AoP would decrease during PVC (Fig. 3b), as shown in animal studies. In total, 25 PVCs were analyzed among 15 patients and were found to increase the HR by 15.5–34.5%. LV contractions following a PVC occurred regularly, so the TR-R following PVCs was temporarily extended, changing the hemodynamic states in the model. Such hemodynamic changes at PVC and at the next beat were observed and compared as the LV stroke volume fluctuated. The HR changes caused by PVCs were applied to two Windkessel models, one with a VAD and one without a VAD. Three episodes of ST occurred in two patients and resulted in simultaneous ejection of the heart and VAD. The effect of co-pulsations was compared with the results for normal counter-pulsation. In addition, the relationship between temporary heart rate changes and the cardiac load was investigated by applying HR changes over a wide range to the model.
a Experimental setup for the flow measurement of the pulsatile VAD, and b the inflow and outflow waveform of the VAD
Results
Representative curves for a patient ECGs, b AoP without VAD perfusion, c AoP with VAD perfusion, d pressure–volume (PV) curves of LV without VAD, and e PV curves of LV with VAD. At first, the patient’s heart beat was normal (①, ②). However, when PVC occurred, the heart rate increased abruptly (③) and then immedicately decreased (④)
The volume of the ventricle was reduced by the inflow to the VAD (Fig. 5e). When PVC occurred and the ratio of Tp-d/TR-R deviated significantly from the target value (45%), the cardiac output and the output of the VAD did not occur simultaneously, because the blood flow output period of the VAD was too short, so the VAD output finished before the start of the following PVC. As shown in Fig. 5b, c, the diastolic AoP increased due to the output of the VAD, while the systolic pressure decreased due to decreases in the cardiac output caused by the flow from the LV to the VAD. The VAD pulses caused the area of the PV curve to decrease (Fig. 5e). Even when the PVC was generated, the VAD reduced the volume of the ventricle, which subsequently reverted to its prior size following the PVC.
Representations of a patient ECGs, b LV volumes, c LV pressures, d AoP and e PV curves during normal (①) and subsequent abnormal heart beat caused by sinus tachycardia (②, ③)
a Predicted cardiac load reduction ratio according to the ratio of Tp-d to TR-R, and b histogram of the ratios of Tp-d to TR-R for observed patients’ heartbeats
Discussion and conclusion
Although PVCs are not dangerous arrhythmias, their frequency affects the occurrence of dangerous arrhythmias such as ventricular fibrillation [28, 29]. PVC data enable analysis of the effect of various types of arrhythmia, because PVC shows fast and slow heart rate changes within a short period. In animal studies using pigs and a PLL-controlled VAD, CPC by the VAD induced clear changes and the systolic blood pressure decreased by 23%, while the diastolic BP increased by 25% [6]; this resulted from the different characteristics of animals compared with humans and edema that was observed at thoracotomy. In vivo, PVCs occurred, when the VAD operated not in CPC mode, but in asynchronous mode. Changes in ventricular load caused by VADs are thought to cause PVCs [30]. As a result, the animal studies did not show the effects of CPC on the hemodynamic state during arrhythmias, because arrhythmias did not occur in the presence of CPC. As the causes of PVCs include alcohol consumption and stress [28], PVC can still occur, even when a VAD imparts CPC; however, they are difficult to detect due to their intermittency. Thus, PVCs can be the prevailing arrhythmia, despite all efforts to avoid the use of VADs in arrhythmia patients.
The results of this study were obtained from simulations; however, the input and output blood flow data of the VAD were measured under the same conditions based on patients and in vitro [26]. The ECG data were measured from patients who required the application of a VAD [21]. The Windkessel model was designed using accepted mathematical and physical formulas relating to hemodynamic states [22]. The parameters in the model, including blood pressure, cardiac output, and systemic vascular resistance, were identical to those in another study and were measured in patients [22, 25]. Since the model includes the AoP waveform, which is similar to that measured in patients, the model is considered to include hemodynamic data [23]. Therefore, the results of this study are expected to be realistic for patients with similar hemodynamic characteristics.
The model used in this study predicted the effects of the VAD on events that were difficult to reproduce in clinical or animal experiments. However, a limitation of this study is that predictions using other input conditions require a substantial volume of data to be collected beforehand. Because the cardiovascular characteristics of patients vary due to multiple factors, including the influence of sympathetic and parasympathetic nervous systems and the patient’s movement and posture, predictions must encompass a wide range of conditions. As the technologies for acquiring data become more encompassing, the limitations of numerical simulation will be minimized [31, 32].
Previous studies of CPC were designed to detect heart rhythms and to maintain the VAD at constant intervals [9, 12]. However, in those studies, the VAD was expected to alter the blood flow irregularly, when irregular heartbeats occurred [33]. It is difficult to predict the outcome of CPC by a VAD under arrhythmia conditions, given that the algorithm accounts for complex scenarios to avoid co-pulsation [12, 19]. Because the Novacor left ventricular assist system (LVAS) (Novacor, Oakland, CA, USA) also has a PLL algorithm, our results should be useful for predicting the hemodynamic effects of the device during an arrhythmia [34].
In this study, the PLL algorithm showed a 97.3% success rate at CPC when tachycardia and arrhythmias appeared in patients and it could reduce the heart load. It predicted that CPC and temporary CPC failure would not induce an increased cardiac load, even when co-pulsation occurred.
Declarations
Authors’ contributions
HK wrote this manuscript. KJY participated in vitro experiments. LKM did the simulation of Windkessel model. CSW analyzed obtained data. All authors read and approved the final manuscript.
Acknowledgements
This work was supported by the National Research Foundation of Korea (NRF) (C1013660-01-01) and was also supported by the Ministry of Trade, Industry & Energy (MOTIE), Korea Institute for Advancement of Technology (KIAT) through the Encouragement Program for The Industries of Economic Cooperation Region (No. R0004491).
Competing interests
The authors declare that they have no competing interests.
Availability of data and materials
All data and material are obtained or produced through legal procedures. They have not been published elsewhere and they are not under consideration by another journal.
Consent for publication
All authors consent for the publication of this manuscript.
Ethics approval and consent to participate
Not applicable.
Funding
National Research Foundation of Korea (NRF) (C1013660-01-01). Ministry of Trade, Industry & Energy (MOTIE), Korea Institute for Advancement of Technology (KIAT) (No. R0004491).
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
Authors’ Affiliations
References
- Slaughter MS, Sihgn R. The role of ventricular assist devices in advanced heart failure. Rev Esp Cardiol. 2012;65(11):983–5.View ArticleGoogle Scholar
- Calvaruso DF, Ocello S. Implantation of a Berlin heart as single ventricle by-pass on Fontan circulation in univentricular heart failure. ASAIO. 2007;53:e1–2.View ArticleGoogle Scholar
- Trivedi JR, Cheng A, Singh R, Williams ML, Slaughter MS. Survival on the heart transplant waiting list: impact of continuous flow left ventricular assist device as bridge to transplant. Ann Thorac Surg. 2009;98:830–4.View ArticleGoogle Scholar
- Frazier OH, Myers TJ, Westaby S, Gregoric ID. Use of the Jarvik 2000 left ventricular assist system as a bridge to heart transplantation or as destination therapy for patients with chronic heart failure. Ann Surg. 2003;237:631–7.Google Scholar
- Schmitto JD, Hanke JS, Rojas SV, Avsar M, Haverich A. First implantation in man of a new magnetically levitated left ventricular assist device. J Heart Lung Transplant. 2015;34(6):858–60.View ArticleGoogle Scholar
- Choi SW, Nam KW, Lim KM, Shim EB, Won YS, Woo HM, Kwak HH, Noh MR, Kim IY, Park SM. Effect of counter-pulsation control of a pulsatile left ventricular assist device on working load variations of the native heart. Biomed Eng Online. 2014;13:35.View ArticleGoogle Scholar
- Lim KM, Kim IS, Choi SW, Min BG, Won YS, Kim HY, Shim EB. Computational analysis of the effect of the type of LVAD flow on coronary perfusion and ventricular afterload. J Physiol Sci. 2009;59:307–16.View ArticleGoogle Scholar
- Shi Y, Korakianitis T, Bowles C. Numerical simulation of cardiovascular dynamics with different types of VAD assistance. J Biomech. 2007;40(13):2919–33.View ArticleGoogle Scholar
- Koshimoto Y, Takewa Y, Arakawa M, Umeki A, Ando M, Nishimura T, Fujii Y, Mizuno T, Nishimura M, Tatsumi E. Development of a novel drive mode to prevent aortic insufficiency during continuous-flow LVAD support by synchronizing rotational speed with heartbeat. J Artif Organs. 2013;16:129–37.View ArticleGoogle Scholar
- Rosenbaum AN, Kremers WK, Duval S, Sakaguchi S, John R, Eckman PM. Arrhythmias in patients with cardiac implantable devices after implantation of a left ventricular assist device. ASAIO. 2016;62(3):274–80.View ArticleGoogle Scholar
- Refaat M, Chemaly E, Lebeche D, Gwathmey JK, Hajjar RJ. Ventricular arrhythmias after left ventricular assist device implantation. Pacing Clin Electrophysiol. 2008;31:1246–52.View ArticleGoogle Scholar
- Nam KW, Chung J, Choi SW, Kim WE, Sun K, Ming BG. Development of counterpulsation algorithm for a moving-actuator type pulsatile LVAD. Int J Artif Organs. 2004;27(10):882–90.View ArticleGoogle Scholar
- Harding JD, Piacentino V 3rd, Gauhan JP, Houser SR, Margulies KB. Electrophysiological alterations after mechanical circulatory support in patients with advanced cardiac failure. Circulation. 2001;104(11):1241–7.View ArticleGoogle Scholar
- Fasseas P, Kutalesk SP, Samuels FL, Holmes EC, Samuels LE. Ventricular assist device support for management of sustained ventricular arrhythmias. Tex Heart Inst J. 2002;29:33–6.Google Scholar
- Ziv O, Dizon J, Thosni A, Naka Y, Magnano AR, Garan H. Effects of left ventricular assist device therapy on ventricular arrhythmias. J Am Coll Cardiol. 2005;45:1428–34.View ArticleGoogle Scholar
- Harding JD, Piacentino V 3rd, Rothman S, Chambers S, Jessup M, Margulies KB. Prolonged repolarization after ventricular assist device support is associated with arrhythmias in humans with congestive heart failure. J Card Fail. 2005;11(3):227–32.View ArticleGoogle Scholar
- Annamalai SK, Buiten L, Esposito ML, Paruchuri V, Mullin A, Breton C, Redicini R, O’kelly R, Morine K, Wessler B, Patel AR, Kieman MS, Karas RH, Kapur NK. Acute hemodynamic effects of intra-aortic balloon counterpulsation pumps in advanced heart failure. J Card Fail. 2017. https://doi.org/10.1016/j.cardfail.2017.05.015.Google Scholar
- Parissis H, Graham V, Lampridis S, Lau M, Hooks G, Mhandu C. IABP: history-evolution-pathophysiology-indications: what we need to know. J Cardiothorac Surg. 2016;11:122.View ArticleGoogle Scholar
- Farrar DJ, Compton PG, Lawson JH, Hershon JJ, Hill JD. Control modes of a clinical ventricular assist device. IEEE Eng Med Biol Mag. 1986;5:19–24.View ArticleGoogle Scholar
- Miura H, Arai S, Sato F, Matsuki H, Sato T. A synchronous rectification using a digital PLL technique for contactless power supplies. IEE Trans Magn. 2005;41:3997–9.View ArticleGoogle Scholar
- Goldberger AL, Amaral LAN, Glass L, Hausdorff JM, Ivanov PCh, Mark RG, Mietus JE, Moody GB, Peng CK, Stanley HE. PhysioBank, PhysioToolkit, and PhysioNet: components of a new research resource for complex physiologic signals. Circulation. 2000;101(23):e215–20.View ArticleGoogle Scholar
- Choi SW. Estimation of blood pressure diagnostic methods by using the four elements blood pressure model simulating aortic wave reflection. J Biomed Eng Res. 2015;36:183–90.View ArticleGoogle Scholar
- Linton NWF, Linton RAF. Estimation of changes in cardiac output from the arterial blood pressure waveform in the upper limb. Br J Anaesth. 2001;86:486–96.View ArticleGoogle Scholar
- Wesseling KH, Jansen JRC, Settels JJ, Schreuder JJ. Computation of aortic flow from pressure in humans using a nonlinear, three-element model. J Appl Physiol. 1993;74:2566–73.View ArticleGoogle Scholar
- Stergiopulos N, Westerhof BE, Westerhof N. Total arterial inertance as the fourth element of the windkessel model. Am J Physiol. 1999;276:H81–8.Google Scholar
- Kang SM, Choi SW. Blood flow and pressure evaluation for a pulsatile conduit-shaped ventricular assist device with structural characteristic of conduit shape. Trans Korean Soc Mech Eng B. 2009;35:1191–8.View ArticleGoogle Scholar
- Hamilton PS, Tompkins WJ. Quantitative investigation of QRS detection rules using the MIT/BIH arrhythmia database. IEEE Trans Biomed Eng. 1986;BME-33(9):1157–65.View ArticleGoogle Scholar
- Yokokawa M, Kim HM, Eric G, et al. Relation of symptoms and symptom duration to premature ventricular complex induced cardiomyopathy. Heart Rhythm. 2012;9:92–5.View ArticleGoogle Scholar
- Baman TS, Lange DC, Ilg KJ, et al. Relationship between burden of premature ventricular complexes and left ventricular function. Heart Rhythm. 2010;7:865–9.View ArticleGoogle Scholar
- Stacy G Jr, Jobe RL, Taylor LK, Hansen DE. Stretch-induced depolarizations as a trigger of arrhythmias in isolated canine left ventricles. Am J Physiol. 1992;263:613–21.Google Scholar
- Wong KKL, Cheung SCP, Yang W, Tu J. Numerical simulation and experimental validation of swirling flow in spiral vortex ventricular assist device. Int J Artif Organs. 2010;33(12):856–67.Google Scholar
- Wong KKL, Tu J, Kelso RM. Vortical flow analysis. J Mech Med Biol. 2010;10(02):191–212.View ArticleGoogle Scholar
- Michael RF. Mechano-electrical feedback in ventricular myocardium. Cardiovasc Res. 1996;32(15–24):30.Google Scholar
- Robbins RC, Kown MH, Portner PM, Oyer PE. The totally implantable novacor left ventricular assist system. Ann Thorac Surg. 2001;71:S162–5.View ArticleGoogle Scholar