Experimental modal analysis reveals the resonant frequencies and associated deflection patterns that exist in metallic and other solid structures with moderate-to-low material viscosity. If the structure is driven at or near one of its resonant frequencies, deflections can become extreme under low viscosity conditions. Thus, by performing a modal analysis, we can determine whether there are any resonant frequencies within the range of frequencies we expect the structure to be driven. Modal analysis was performed on a HeartMate II LVAD using a PCB Piezotronics ceramic piezoelectric accelerometer (model 352A10; sensitivity, 10.40 mV/g; mass, 0.7 g), a PCB Piezotronics model 086C01 modally tuned impulse hammer with force transducer (PCB Piezotronics, Depew, NY), and an Agilent 35670A Dynamic Signal Analyzer (Agilent Technologies, Santa Clara, CA). With the LVAD suspended by its driveline, the accelerometer was glued to a fixed point in the center of the LVAD using cyanoacrylate. Acceleration was measured radially with instrumented hammer taps at 10 points distributed circumferentially around the LVAD, 3 hammer taps per point to ensure repeatability, which was assessed by confirming near unity coherence between measurements in the frequency range of interest (0–2,000 Hz).15,16 Data were transferred to a personal computer and analyzed using custom MATLAB software.
Measurements in Human Subjects
Ten HeartMate II patients with stable LVAD function, hemodynamics, and cardiac, pulmonary, and renal function were auscultated in the outpatient setting following a consent process per an institutional review board–approved protocol. Recordings were made at the surface of the skin or through minimal clothing while the patient was sitting up and while the patient’s LVAD was connected to a pump monitor displaying instantaneous RPM, flow, pulsatility index, and power consumption information. Recordings were made at the mitral and tricuspid positions in each case. Pulse rate, blood pressure, and LVAD operation data were collected at the time of auscultation. Echocardiographic and laboratory data were collected at the closest incidence to the time of auscultation.
Two additional patients (henceforth referred to as patients A and B), separate from the group of 10 stable patients, were monitored before and after a pump exchange procedure for device thrombosis. Both patients had HeartMate II LVADs. Echocardiographic, laboratory, and hemodynamic data were collected before and after surgery. To investigate potential acoustic effects of thrombus, acoustic measurements were made before and after exchange in each patient.
Spectrographic analysis revealed the presence of pronounced banding structure with a regular low-amplitude and high-amplitude striping pattern in the spectrogram plot (graphical representation of frequency versus time) and distinct frequency peaks in the spectrogram slice (graphical representation of amplitude versus frequency). The frequency spectra of an implanted HeartMate II are shown in Figure 2A. Frequency bands occur at regular intervals and tend to decrease in amplitude as frequency increases. When plotted against harmonic order, rather than a frequency scale, it can be seen that these frequency peaks occur at harmonic intervals of the LVAD operational speed (Figure 2B). The relationship between expected harmonic intervals, calculated through Equation 1 and the measured peak frequencies at each band, was found to be consistent both in vivo and in vitro. This was demonstrated by strong correlation between measured and predicted peak frequency values (all r > 0.999).
Hλ = expected harmonic frequency for a given pump speed, n = harmonic number, and RPM = pump speed in revolutions per minute.
In Vitro Results
Frequency band localization analysis at LVAD speeds ranging from 7,000 to 12,000 RPM indicated the linear relationship between LVAD speeds and harmonic frequency. Flow (liters per minute) and power consumption (watts) increased with LVAD speed, although pulsatility index was not affected. In addition, LVAD speed-matched correlations were calculated between peak harmonic frequencies for each patient and their corresponding speed measured in vitro. All r values were greater than 0.999.
Although changes in inflow pressure, outflow pressure, and fluid viscosity affected flow rates and power consumption, no changes in peak harmonic frequency or amplitude occurred as these fluid properties were altered. Presumably, under these conditions, the fluid flow path through the LVAD conduit remained normal, and there was no significant departure from device operation typical of normal physiological use. In addition, modal analysis revealed no inherent structural resonance at the frequencies measured by the electronic stethoscope (0–2 kHz).15 The lowest frequency flexural mode of resonant vibration of the LVAD occurred at 8.5 kHz, approximately 4.25 times higher than the LVAD frequencies we recorded, indicating that all acoustic signals measured in the range of 0–2000 Hz were generated by operation of the device, not by structural resonance.
Results in Human Subjects
A total of 10 HeartMate II supported patients with mean age of 56.4 ± 12.2 years and BMI of 28.3 ± 4.8 kg/m2 were auscultated in our outpatient clinic. These patients were known to be in stable condition and relatively free of complications at the time of auscultation. The average duration on LVAD was 549.5 ± 334.7 days, and LVAD, lung, and kidney function were noted to be normal. Echocardiography revealed 3 patients (30%) with mild mitral regurgitation, 3 patients (30%) with mild aortic insufficiency, 6 patients (60%) with mild tricuspid regurgitation, 2 patients (20%) with mild pulmonic insufficiency, and no patients with mitral or tricuspid valve replacements (Table 2). In addition, hemodynamic and laboratory values did not suggest the presence of hemolysis or thrombus in any of the pumps.
Spectrogram slices for each of the 10 patients are shown in Figure 2B, and spectrographic analysis revealed that pump acoustics in vivo are very similar to pump acoustics in vitro. Left ventricular assist device speed, power consumption, flow, and patient blood pressures at the time of recording are indicated in Table 2. Occasionally, short frequency ranges and individual harmonics exhibited lower signal-to-noise ratio and obfuscated peaks in in vivo recordings compared with in in vitro recordings. Despite signal artifacts, which may be a result of a considerably more complex acoustic environment in vivo, peak harmonic frequencies occurred within a very small range of the expected values calculated by Equation 1 (r values between predicted and measured harmonic frequencies >0.999 for all patients). This was an expected result but indicates that minimal frequency distortion occurred during acoustic transmission through the human chest or the gel model used in vitro.
The location of the pump harmonics in these 10 patients again correlated well with expected values and with measured in vitro values for equivalent speeds (r > 0.999 for all). In 6 of the 10 patients, average amplitudes across all 12 recorded harmonics were significantly less than those measured in vitro at equivalent pump speeds.
By using custom MATLAB software, AUCs were measured for each patient’s spectrogram slice in 50 Hz segments over the 2,000 Hz spectrum. Comparison of in vivo AUC segments and in vitro AUC segments at each patient’s corresponding speed showed strong correlations in all incidences (r > 0.966 for all AUC segments; Figure 3).
The AUC was compared for acoustic samples taken during normal in vitro pump function and while the inflow graft was occluded to 80%. Although statistical differences in AUC were nonsignificant, the spectral slices shown in Figure 4 indicate that the two samples are closely related in behavior, but with unique amplitude values. Previous work has shown that diameter changes of less than 20% in hemodialysis grafts are detectable using a surface sensor. The presence of a stenosis is expected to increase vessel resistance and reduce fluid velocity through the LVAD, reducing vascular sounds.7 Presumably the reduction in amplitude in the occluded measurement in our study is a result of impaired pump function, specifically reduced flow through the LVAD circuit. Identical analysis for the patients who underwent pump exchange for thrombosis revealed a similar trend, suggesting that our in vitro simulation of pump dysfunction is a good representation of acoustic changes caused by thrombus or impaired flow in live patients (Figures 5 and 6).
Acoustic changes resulting from device thrombosis were investigated in patients A and B. Patient A presented with lactate dehydrogenase (LDH) levels of 4,530 U/L, plasma-free hemoglobin levels of 11.1 mg/dl, and increases in pump power to approximately 12 W. Patient B presented with LDH levels of 1,201 U/L, increasing heart failure symptoms, and continued LVAD alarms. Surgical LVAD exchange was deemed necessary in both patients.
Bronchoscopic exploration of the explanted LVADs revealed a thrombus at the inflow stator bearing, upstream of the impeller in patient A and within the device inflow cannula and pump housing in patient B. Comparison of spectra from the clotted and nonclotted pumps for both patients are shown in Figures 5 and 6. In both cases, the presence of thrombus is associated with appreciable reduction in spectral energy over the sample frequency range. When frequency is normalized to harmonic order by adjusting for device RPM, curve morphology is largely preserved, despite reduction in amplitude.
Although heart transplantation remains the most effective treatment of advanced heart failure refractory to medical therapy, because of donor shortages, LVADs have become a popular and established method to restore patient functionality and quality of life.17–19 However, there is need for noninvasive evaluation of pump function to predict complications such as thrombus. In 2005, Mansy et al.7 published a means of detecting vascular patency in hemodialysis grafts using an electronic stethoscope with custom MATLAB software for acoustic analysis. By using this system, the group was able to establish correlation between acoustic power and the degree of stenosis. In 2003, Tanishiro et al.20 showed that arterial sound can be used to detect malrotation in centrifugal LVADs in an in vitro setup, and in 2007, Slaughter et al.21 showed that acoustic changes predicted impending device malfunction in the pulsatile HeartMate XVE. Recently, Hubbert et al.22 showed significant differences in HeartMate II acoustics when inflow and outflow tubing were subject to artificial stenosis and when artificial clots were passed through the pump, and Kaufmann et al.23 showed that thrombus formation in HeartWare ventricular assist device pumps was associated with telltale changes in acoustic spectra. Similar data collection techniques were used in all studies, although particular equipment differed. It should be noted that Kaufmann et al. focused data analysis on the frequency range 0–600 Hz and that the Hubbert group assessed acoustic signatures from 0 to approximately 22 kHz. Both groups analyzed changes in acoustic amplitude resulting from device malfunction (thrombosis or simulated thrombosis). To our knowledge, this study is the first to pair clinical and laboratory investigation of HeartMate II acoustics and is the first to compare preexplantation and postexplantation spectra for patients with confirmed HeartMate II thrombosis who underwent surgical exchange.
Our results indicate that acoustic analysis is an accurate means of detecting both the location of peak harmonic frequencies and the behavior of a spectrographic slice as measured by AUC segments. The in vitro and in vivo data show strong correlations with expected values and with each other. This suggests that neither the thoracic nor the experimental acoustic environments cause frequency shifts in the measured ranges. However, both amplitude and AUC were lower in vivo than in vitro for all patients (as shown for patient 7 in Figure 3). This is likely because of amplitude attenuation in the chest where the acoustic environment is considerably more complex than the gel phantom used in vitro; transmission across the muscle, bone, fibrous tissue, and fluid may reduce the magnitude of acoustic waves. Analysis of the spectrogram shown in Figure 2A indicates that acoustics generated by the HeartMate II are not changed significantly during the cardiac cycle. The independence of HeartMate II acoustics from changes in pressure is supported by our observation that device acoustics were not significantly altered in vitro by changes in system pressures. Further, cardiac sounds, including those of valve closure, generally occur at frequencies less than 200 Hz and are largely distinct from those associated with LVAD operation.
The strong correlation between the measured peak harmonic frequencies is an expected result and indicates that this methodology is a sensitive means of detecting device rotational speed. However, analysis of 50 Hz AUC segments provides greater quantitative insight into how the spectral slice is behaving. We suggest that the areas between frequency peaks may be related to LVAD function, including flow through the device, and that calculation of AUC segments allows for quantification across all measured frequencies. The strong correlations between in vitro and in vivo AUC segments suggest that this too is a powerful means of quantifying LVAD acoustics. Further, the mirroring of in vitro and in vivo spectral slice curves shown in Figure 3 suggests that certain characteristics of the spectral slice are preserved in both laboratory and clinical measurements.
The data presented in Figures 5 and 6 support this conclusion and suggest that LVAD auscultation may be an effective means of detecting acoustic changes caused by the presence of pump thrombus in vivo. When the inflow graft was occluded in vitro, a trend of curve mirroring with reduction in amplitude was observed, as shown in Figure 4. This supports the findings of Hubbert et al.22 who detected a similar pattern when ball valves were used to occlude the inflow and outflow conduits. The similarities between these in vitro changes because of partial occlusion and the effect of thrombus in vivo may be related to the location of the thrombus at the inflow stator of the pump, upstream of the impeller in patient A, and within the inflow tract and device housing in patient B. We suggest that in the cases of partial occlusion in vitro and thrombosis in vivo, a reduction in LVAD amplitude compared with baseline measurements is related to reduced flow through the device. Notably, this trend was observed clinically in both patients with device thrombosis and may represent an early, noninvasive marker for impending device failure. This study analyzed only two patients with confirmed LVAD thrombosis. Continued study with larger sample sizes is required to determine whether reduction in amplitude may be used as a diagnostic of LVAD thrombus.
At this time, we cannot explain some of the characteristics common in many acoustic samples such as a strong attenuation at approximately 700 Hz. However, the consistently strong third harmonic amplitude is likely because of the fact that the pump impeller has three blades, causing fluctuations in flow at three times the rotation frequency of the pump. Modal analysis did not reveal any inherent resonances in the device housing, leading to the conclusion that all acoustic behavior measured in the sample range 0–2,000 Hz is resultant to LVAD operation.
Acoustic sampling and analysis is a sensitive and accurate method of detecting the presence, location, and individual amplitude of LVAD peak harmonic frequencies and AUC in spectrogram slices. Our methodology allows for characterization of pump function and may represent an important means for noninvasive detection of LVAD thrombosis. We are currently using this methodology to study the relationship of these acoustic parameters to pump complications, patient hemodynamics, and patient cardiac physiology.
The authors thank the Sejal Modi CCP, the nurses of the Advocate Christ VAD clinic, and the members of the UIC acoustics and vibrations laboratory.
G.L.Y. helped in concept/design, data collection, data analysis/interpretation, drafting, revision, and statistics; T.J.R. helped in concept/design, data interpretation, revision of article, and approval of article; G.B. helped in concept/design, data interpretation, revision of article, and approval of article; and A.J.T. helped in approval of the article.
1. Burke DJ, Burke E, Parsaei F, et al. The Heartmate II: Design and development of a fully sealed axial flow left ventricular assist system. Artific Organs. 2001;25:380–385
2. Fine NM, Topilsky Y, Oh JK, et al. Role of echocardiography in patients with intravascular hemolysis due to suspected continuous-flow LVAD thrombosis. JACC Cardiovasc Imaging. 2013;6:1129–1140
3. Acharya D, Hashim T, Kirklin JK, et al. Use of retrospectively gated CT angiography to diagnose systolic LVAD inflow obstruction. ASAIO J. 2013;59:542–546
4. Birati EY, Rame JE. Diagnosis and management of LVAD thrombosis. Curr Treat Options Cardi Med. 2015;17:361–372
5. Cowger JA, Romano MA, Shah P, et al. Hemolysis: A harbinger of adverse outcome after left ventricular assist device implant. J Heart Lung Transplant. 2014;33:35–43
6. Hasin T, Deo S, Maleszewski JJ, et al. The role of medical management for acute intravascular hemolysis in patients supported on axial flow LVAD. ASAIO J. 2014;60:9–14
7. Mansy HA, Hoxie SJ, Patel NH, Sandler RH. Computerised analysis of auscultatory sounds associated with vascular patency of haemodialysis access. Med Biol Eng Comput. 2005;43:56–62
8. Kompis M. Acoustic imaging of the human chest. Chest. 2001;120:1309–1321
9. Charif RA, Waack AM, Strickman LM Raven Pro 1.4 User’s Manual. 2010 Ithaca, NY Cornell Lab of Ornithology
10. Kay SM, Marple SL. Spectrum analysis—A modern perspective. Proc IEEE. 1981;69:1380–1481
11. Bioacoustics Research Program. Raven Pro: Interactive Sound Analysis Software (Version 1.5) [Computer Software]. 2014. Ithaca, NY The Cornell Lab of Ornithology Available at: http://www.birds.cornell.edu/raven
. Accessed November 2, 2014
12. Royston TJ, Dai Z, Chaunsali R, Liu Y, Peng Y, Magin RL. Estimating material viscoelastic properties based on surface wave measurements: A comparison of techniques and modeling assumptions. J Acous Soc Amer. 2011;130:4126–4138
13. Yasar TK, Royston TJ, Magin RL. Wideband MR elastography for viscoelasticity model identification. Magn Reson Med. 2013;70:479–489
14. Dai Z, Peng Y, Mansy HA, Sandler RH, Royston TJ. Comparison of poroviscoelastic models for sound and vibration in the lungs. J Vib Acoust. 2014;136:1–11
15. Fu HJ, He ZF. Overview of modal analysis Modal Analysis. 2001 Woburn, MA Butterworth-Heineman:15–25
16. Rossing TRossing T. Modal analysis Springer Handbook of Acoustics. 2007 New York, NY Springer-Verlag New York:1127–1138
17. Lund LH, Edwards LB, Kucheryavaya AY, et al.International Society for Heart and Lung Transplantation. The Registry of the International Society for Heart and Lung Transplantation: Thirtieth Official Adult Heart Transplant Report—2013; focus theme: Age. J Heart Lung Transplant. 2013;32:951–964
18. Pagani FD, Miller LW, Russell SD, et al.HeartMate II Investigators. Extended mechanical circulatory support with a continuous-flow rotary left ventricular assist device. J Am Coll Cardiol. 2009;54:312–321
19. Rose EA, Gelijns AC, Moskowitz AJ, et al.Randomized Evaluation of Mechanical Assistance for the Treatment of Congestive Heart Failure (REMATCH) Study Group. Long-term use of a left ventricular assist device for end-stage heart failure. N Engl J Med. 2001;345:1435–1443
20. Tanishiro H, Funakubo A, Fukui Y. Arterial sound based noninvasive malrotation detection of rotary LVAD. ASAIO J. 2004;50:306–310
21. Slaughter MS, Feldman CM, Sobieski MA, et al. Utilization of acoustic signatures to identify HeartMate XVE device end-of-life. J Heart Lung Transplant. 2007;26:579–583
22. Hubbert L, Sundbom P, Loebe M, Peterzén B, Granfeldt H, Ahn H. Acoustic Analysis of a mechanical circulatory support. Artific Organs. 2014;38:593–598
23. Kaufmann F, Hörmandinger C, Stepanenko A, et al. Acoustic spectral analysis for determining pump thrombosis in rotary blood pumps. ASAIO J. 2014;60:502–507
left ventricular assist device (LVAD); acoustics; spectral analysis; thrombosisCopyright © 2016 by the American Society for Artificial Internal Organs