 Research
 Open Access
 Published:
On the relationship between competitive flow and FFT analysis of the flow waves in the left internal mammary artery graft in the process of CABG
BioMedical Engineering OnLinevolume 15, Article number: 129 (2016)
Abstract
Background
During coronary artery bypass grafting (CABG), the ratio of powers of the fundamental frequency and its first harmonic (F0/H1) in fast Fourier transformation (FFT) analysis of the graft’s flow waves has been used in the field of evaluation of the patency in anastomosis. But there is no report about using the FFT method to evaluate the magnitude of competitive flow. This study is aiming at exploring the relationship between competitive flow and FFT analysis of the flow waves in left internal mammary artery (LIMA) graft, and finding a new method to evaluate the magnitude of competitive flow.
Methods
At first, establishing the CABG multiscale models of different stenosis in left anterior descending artery (LAD) to get different magnitude of competitive flows. Then, calculating the models by ANSYSCFX and getting the flow waves in LIMA. Finally, analyzing the flow waves by FFT method and comparing the FFT results with the magnitude of competitive flow.
Results
There is no relationship between competitive flow and F0/H1. As for F0/H2 and F0/H3, they both increase with the reduction of the stenosis in LAD. But the increase of F0/H3 is not obviously enough and it can’t identify the significant competitive flow clearly, so it can’t be used as the evaluation index. It is found that F0/H2 increases obviously with the increase of the competitive flow and can identify the significant competitive flow.
Conclusion
The FFT method can be used in the evaluation of competitive flow and the F0/H2 is the ideal index. High F0/H2 refers to the significant competitive flow. This method can be used during CABG to avoid the risk of competitive flow.
Background
The flow waves in graft after CABG can reflect whether the surgery is successful. The average flow (Qm), pulsatility index (PI), diastolic fraction (DF), the ratio of the diastolic peak flow and systolic peak flow (D/S) and the percentage of backflow (%Insuf) derived from the flow waves are the important indexes in evaluating the efficiency of the surgery [1–4]. Currently, there have been studies finding that the F0/H1 derived from FFT analysis of flow waves in grafts can evaluate the quality and graft patency, and it is more accurate and reliable than the above indexes [5–8]. Because of the good performance of Fast Fourier transformation (FFT) method, the FFT ratio has been used to differentiate patent grafts [3] and it is put forward as one of the modern diagnostic tools in intraoperative graft patency verification [9].
Competitive flow is the main risk factor influencing the artery graft after CABG [10, 11]. It can decline the flow and generate backflow in the graft, and cause the graft failure due to the string phenomenon [12–15]. So, how to evaluate the competitive flow is a serious problem which researchers and surgeons most concerned. The most effective method evaluating the competitive flow is to clamp the stenosis coronary artery and observe the change of the graft’s flow wave. The significant competitive flow is identified if the flow increases obviously [16]. But, clamping the patient’s coronary artery during the surgery is a dangerous action, it might cause the plaque fall off and block downstream vessel, and it will injure the myocardium. Therefore, evaluating the competitive flow by clamping the stenosis coronary artery is forbidden in clinic. It is necessary to find a method to evaluate the competitive flow without injuring the myocardium. In view of the outstanding performance of the FFT method in evaluating the patency of the anastomosis, this study is going to explore the possibility of using the FFT method to evaluate the magnitude of competitive flow.
Because of the danger of the clinical trials and the difficulty of controlling the graft flow only affected by the magnitude of competitive flow, this study chooses the method of constructing CABG multiscale models of different stenosis. Multiscale models not only guarantee the real geometry of the coronary artery, but also provide precise boundary condition [17, 18]. It has a high accuracy of simulating the physiological conditions, and it has been used in the calculation of various vessels’ hemodynamics. Multiscale models can change the stenosis ratio while other conditions keep constant, and obtain flow waves from models of different stenosis. As there is significant correlation between stenosis and competitive flow, these models can reflect different magnitude of competitive flow.
This study is aiming at exploring the relationship between competitive flow and FFT analysis of the flow waves in LIMA graft, and finding a new way to evaluate the magnitude of the competitive flow by FFT method.
Methods
The CABG multiscale model
The multiscale model is constructed by different dimensional models. In this study, the 0D/3D coupling method is used to perform a numerical simulation by coupling the lumped parameter model (0D model) and 3D model [19].
The reconstruction of 3D model
The patient data is from a male whose cardiac output is 4.6 l/min measured by Doppler ultrasound. The CT data for reconstruction has 460 slices, each of this slice is 512 * 512 pixels, and the distance between each adjacent slice is 1 mm. The patients signed informed consent forms. The threshold segmentation and manual segmentation has been used for reconstruction in Mimics software, and the Gaussian filtering is adopted for smoothing the model surface. Using the Freeform software (software of the 3D modeling system), the model is made different stenosis (100, 90, 75, 60, 50 and 40%) in LAD and completes the LIMALAD bypass. The diameter of LIMA is set up as 3 mm. Finally, all of the six models with different stenosis need to be meshed to generate the computational models. The ANSYSCFX software is used to generate hexahedral mesh in models, and the mesh is refined in the areas of interest to make the simulation results more precise. The node and element numbers for the models are listed in Table 1.
In the 3D simulation, it is assumed that the vessel wall is rigid and the blood is incompressible viscous Newtonian fluid. The density of the blood is 1050 kg/m^{3}, and the dynamic viscosity is 0.0035 Pa s.
The 0D/3D coupling model
In this study, the 0D model which act as the boundary condition for 3D model is described by Taylor [18]. The complete 0D/3D coupling model (multiscale model) is shown as Fig. 1.
The boundary condition of the 3D model are supplied by the 0D calculation and the forcing terms of the 0D model are calculated by the 3D simulation. All of the six 3D models shares the same 0D models, so the difference in the simulation can be caused only by the difference of the stenosis.
As shown in Fig. 1, the inlet of the 3D model (Aorta) is coupled with a 0D model of the heart block. The systemic outlets (A, B, C, D, E) are coupled with a 0D model of the systemic block. And the coronary outlets (an) are coupled with a 0D model of the coronary block. In each block, the resistance (R) is used to simulate the flow resistance, the capacitance (C) is used to simulate the compliance of the vessel and the inductance (L) is used to simulate the inertia of the blood flow. In the heart block, the diode is used to simulate the valve in the heart while the variable capacitors Clv is used to simulate the left ventricles. In the coronary block, the pressure of left ventricle is added at the cathode of capacitance Cim which accounts for the intramyocardial pump [18].
In the heart block, the pressure–volume relationship of the ventricle can be described as the following functions.
where E(t) is the timevarying elastance (mmHg/ml), V(t) and P(t) are the ventricle volume (ml) and pressure (mmHg) respectively, and V _{0} (ml) is the reference volume. Mathematically, the function can be approximated as:
where E _{ n }(t _{ n }) is the normalized timevarying elastance [20].
where \(t_{n} = \frac{t}{{T_{max} }}\), T _{ max } = 0.2 + 0.15t _{ c }, and t _{ c } is one cardiac cycle interval (s). In this paper, it is set that E _{ max } = 2.0, E _{ min } = 0.002458, t _{ c } = 0.8 s.
The values of parameters in the 0D model are based on data from the research of coronary artery model by Kim [21]. In Kim’s process, it adjusts the parameters of the heart block to make the cardiac output and the aortic pressure calculated by the model fit the clinical measuring values. The total flow of coronary artery is 4.0% of the cardiac output. And that the flow distribution in each coronary branch is consistent with its diameter. Adjusting the parameters of every coronary branch to fit above flow rules. In our models, we set Kim’s parameters’ data as initial values, and cardiac output and aortic pressure are set as objective functions. By using the genetic algorithm, these values are adjusted to make sure that the cardiac output, the systolic and diastolic pressures matches the clinical data, All of the values of parameters are shown in Table 2. Table 3 demonstrates the comparison of the cardiac output, the systolic and diastolic pressures between model predictions and clinical data. It is known that the model data and clinical data have a satisfactory similarity which means the values of parameters are appropriate.
Now, the complete CABG multiscale model has been constructed and can be calculated by ANSYSCFX. The flow waves are extracted from LIMA of the six models for the following analysis.
FFT analysis to flow waves
Fast Fourier transformation method is a classic signal processing method. It is valuable for clinic using this method to analyze the graft flow waves, and currently, this method has been used in evaluating the patency of anastomosis.
The powers of different frequencies will be obtained through FFT analysis, just as Fig. 2.
In the spectrogram, the first peak is the powers of fundamental frequency (F0), the second peak is the first harmonic (H1). And then, the following are the second harmonic (H2) and the third harmonic (H3), etc.
The FFT analysis of the LIMA flow waves from different stenosis models have been done to explore the possibility of using this method to evaluate the competitive flow.
To prevent the spectrum leakage, the periodic extension is adopted and we take 10 cycles for calculation. In this FFT analysis, the sampling frequency is set up as 50 Hz and the number of sampling points is 2048.
Results
The relation of stenosis ratio and competitive flow
In our CABG models, the flow in stenosis LAD and LIMA are gotten. The LIMA flow waves of different stenosis are shown in Fig. 3.
From the picture, it is known that the graft flow decreases with the reduction of the stenosis ratio. Furthermore, the backflow appears when the stenosis ratio is 50%. The presence of backflow means the increase of the oscillatory shear index (OSI) which is a dangerous hemodynamic index for the efficiency of the graft.
The ratio of the flow of the stenosis LAD and LIMA is used to evaluate the magnitude of competitive flow. The high ratio means high magnitude of competitive flow. The flow of stenosis LAD and LIMA, and the ratio of the flow of LAD and LIMA are recorded in Table 4.
As is known from the Table 4, with the reduction of the stenosis, the LAD/LIMA ratio is gradually increasing. This means the magnitude of competitive flow increases. The ratio has a significant skip when the stenosis is 50%. In this stenosis, the flow in stenosis LAD exceeds the flow in LIMA.
The FFT analysis of different stenosis
The FFT results of different stenosis have shown in Fig. 4.
Recording F0, H1, H2, H3 of different stenosis, and then calculate F0/H1, F0/H2, F0/H3. These values are shown in Table 5.
According to Fig. 4 and Table 5, it is known that F0/H1 increases when the stenosis is more than 75%, and then the ratio decreases. F0/H2 and F0/H3 both increase with the reduction of the stenosis. However, the increasing rate of F0/H2 is faster, and there is an obvious skip in 50%. Drawing the stenosisvarying line chart of F0/H1, F0/H2 and F0/H3, shown as Fig. 5.
From the Fig. 3 and Table 4, we know that the magnitude of competitive flow increases with the reduction of the stenosis. Moreover, there is a significant skip of competitive flow accompany with backflow in 50%. In the light of Fig. 5, it is known that there is no relationship between F0/H1 and the magnitude of competitive flow, while F0/H2 and F0/H3 both increase with the increase of the competitive flow. But the curve of F0/H3 is too gentle to distinguish the different magnitude of competitive flow. The F0/H2 curve has an obvious difference among different magnitude of competitive flow, and has a significant skip in 50% compared to 60% which is consistent with the actual situation. Therefore, the index of F0/H2 can reflect the magnitude of competitive flow commendably.
Discussion
In this paper, we use the 0D/3D coupling model method to construct a CABG model of different stenosis. It can ensure that LIMA flow waves only influenced by stenosis ratio. By comparing stenosis LAD flow and LIMA flow, it can be found that the magnitude of competitive flow increases with the reduction of stenosis ratio. And the competitive flow has a significant increase accompany with the backflow in 50% stenosis.
By doing FFT analysis to LIMA flow waves of different stenosis, we can get F0/H1, F0/H2 and F0/H3. Finding that F0/H1 and F0/H3 have defects in evaluating the competitive flow, while F0/H2 has an outstanding performance. The F0/H2 increases obviously with the increase of the competitive flow and can distinguish the significant competitive flow in 50% stenosis. So the F0/H2 can be used as an effective index in evaluating the competitive flow.
Therefore, to evaluate the competitive flow, the FFT analysis should be taken and F0/H2 is calculated after obtaining LIMA flow waves in clinic. High F0/H2 suggests significant competitive flow. This method avoids defects existed in the method which needs to clamp the stenosis LAD. It can evaluate competitive flow effectively while there is no damage to the heart and instruct surgeons to take measures to avoid the harm of competitive flow.
However, there is no reports published about the physiological meaning of FFT ratio. In our study, it is found that F0/H1 might be related to the ratio of systolic peak flow and diastolic peak flow, while F0/H2 and F0/H3 are lack of research. For future work, it is necessary to study the physiological meaning of FFT ratio and this will help researchers recognize and improve the method.
Conclusion
The FFT analysis could be used to evaluate the magnitude of competitive flow and F0/H2 is the index which is more effective than F0/H1 and F0/H3. High F0/H2 suggests significant competitive flow, while low F0/H2 suggests light competitive flow.
Abbreviations
 CABG:

coronary artery bypass graft
 FFT:

fast Fourier transformation
 LIMA:

left internal mammary artery
 LAD:

left anterior descending artery
 F0/H1:

the ratio of powers of the fundamental frequency and its first harmonic
 F0/H2:

the ratio of powers of the fundamental frequency and its second harmonic
 F0/H3:

the ratio of powers of the fundamental frequency and its third harmonic
References
 1.
Tokuda Y, Song MH, Ueda Y, et al. Predicting early coronary artery bypass graft failure by intraoperative transit time flow measurement. Ann Thorac Surg. 2007;84(6):1928–33.
 2.
Tokuda Y, Song MH, Oshima H, et al. Predicting midterm coronary artery bypass graft failure by intraoperative transit time flow measurement. Ann Thorac Surg. 2008;86(2):532–6.
 3.
Kim KB, Kang CH, Lim C. Prediction of graft flow impairment by intraoperative transit time flow measurement in offpump coronary artery bypass using arterial grafts. Ann Thorac Surg. 2005;80(2):594–8.
 4.
Di Giammarco G, Pano M, Cirmeni S, et al. Predictive value of intraoperative transittime flow measurement for shortterm graft patency in coronary surgery. J Thorac Cardiovasc Surg. 2006;132(3):468–74.
 5.
Takami Y, Ina H. Relation of intraoperative flow measurement with postoperative quantitative angiographic assessment of coronary artery bypass grafting. Ann Thorac Surg. 2001;72(4):1270–4.
 6.
Une D, Chikazawa G, Karkhanis R, et al. Can Fast Fourier Transformation (FFT) analysis of graft flow predict 1 year graft failure after CABG? Can J Cardiol. 2011;27(5):S284.
 7.
Hatada A, Yoshimasu T, Kaneko M, et al. Relation of waveform of transittime flow measurement and graft patency in coronary artery bypass grafting. J Thorac Cardiovasc Surg. 2007;134(3):789–91.
 8.
Takami Y, Ina H. A simple method to determine anastomotic quality of coronary artery bypass grafting in the operating room. Vascular. 2001;9(5):499–503.
 9.
Fabio Bartolozzi MD, Michele Pilato MD. Intraoperative graft patency verification in coronary surgery: modern diagnostic tools. Intraoperative Flow Meas Coron Artery Surg. 2009;23(2):129.
 10.
Sabik JF, Lytle BW, Blackstone EH, et al. Does competitive flow reduce internal thoracic artery graft patency? Ann Thorac Surg. 2003;76(5):1490–7.
 11.
Pagni S, Storey J, Ballen J, et al. Factors affecting internal mammary artery graft survival: how is competitive flow from a patent native coronary vessel a risk factor? J Surg Res. 1997;71(2):172–8.
 12.
Villareal RP, Mathur VS. The string phenomenon: an important cause of internal mammary artery graft failure. Tex Heart Inst J. 2000;27(4):346.
 13.
Beijk MA, Harskamp RE. Treatment of coronary artery bypass graft failure. In: Aronow WS (ed) Artery bypass. InTech; 2013. p. 193–237.
 14.
Barner HB. Double internal mammarycoronary artery bypass. Arch Surg. 1974;109(5):627–30.
 15.
Seki T, Kitamura S, Kawachi K, et al. A quantitative study of postoperative luminal narrowing of the internal thoracic artery graft in coronary artery bypass surgery. J Thorac Cardiovasc Surg. 1992;104(6):1532–8.
 16.
Bolotin G, Kypson AP, Nifong LW, et al. A technique for evaluating competitive flow for intraoperative decision making in coronary artery surgery. Ann Thorac Surg. 2003;76(6):2118–20.
 17.
Sankaran S, Moghadam ME, Kahn AM, et al. Patientspecific multiscale modeling of blood flow for coronary artery bypass graft surgery. Ann Biomed Eng. 2012;40(10):2228–42.
 18.
Taylor CA, Fonte TA, Min JK. Computational fluid dynamics applied to cardiac computed tomography for noninvasive quantification of fractional flow reserve: scientific basis. J Am Coll Cardiol. 2013;61(22):2233–41.
 19.
Zhao X, Liu Y, Li L, et al. Hemodynamics of the string phenomenon in the internal thoracic artery grafted to the left anterior descending artery with moderate stenosis. J Biomech. 2016;49(7):983–91.
 20.
Stergiopulos N, Meister JJ, Westerhof N. Determinants of stroke volume and systolic and diastolic aortic pressure. Am J Physiol Heart Circ Physiol. 1996;270(6):H2050–9.
 21.
Kim HJ, VignonClementel IE, Coogan JS, et al. Patientspecific modeling of blood flow and pressure in human coronary arteries. Ann Biomed Eng. 2010;38(10):3195–209.
Declarations
Authors’ contributions
BM was responsible for the design and signal processing. WW and XZ were responsible for model construction. LL and HZ were responsible for the parameter adjusting. ZZ was responsible for the data collection. YL instructed all works. All authors have made substantial contributions to conception and design, or acquisition of data, or analysis and interpretation of data and have been involved in drafting the manuscript or revising it critically for important intellectual content. Each author has participated sufficiently in the work to take public responsibility for appropriate portions of the content. All authors read and approved the final manuscript.
Competing interests
The authors declare that they have no competing interests.
About this supplement
This article has been published as part of BioMedical Engineering OnLine Volume 15 Supplement 2, 2016. Computational and experimental methods for biological research: cardiovascular diseases and beyond. The full contents of the supplement are available online http://biomedicalengineeringonline.biomedcentral.com/articles/supplements/volume15supplement2.
Availability of data and materials
Data available upon request.
Ethics approval and consent to participate
The study was approved by the Ethic Committee of Beijing University of Technology. Witten sign consent was obtained from the participant.
Funding
Publication of this article was funded by the National Natural Science Foundation of China (11172016, 11472022), the Advanced Medical Instruments Development project of Beijing University of Technology.
Author information
Rights and permissions
Open Access This 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.