Integrating Patient-Specific Electrocardiogram Signals and Image-Based Computational Fluid Dynamics Method to Analyze Coronary Blood Flow in Patients during Cardiac Arrhythmias



The aim of this study was to use the computational fluid dynamics (CFD) method, patient-specific electrocardiogram (ECG) signals, and computed tomography three-dimensional image reconstruction technique to investigate the blood flow in coronary arteries during cardiac arrhythmia.


Two patients with premature ventricular contraction-type cardiac arrhythmia and one with atrial fibrillation-type cardiac arrhythmia were investigated. The inlet velocity of the coronary artery in simulation was applied with the measured velocity profile of the left ventricular outflow tract (LVOT) from the Doppler echocardiography. The measured patient central aortic blood pressure waveform was employed for the coronary artery outlet in simulation. The no-slip boundary condition was applied to the arterial wall.


For the patient with irregular cardiac rhythms (Case I), the coronary blood flow rate under the shortened and lengthened cardiac rhythms were 0.66 and 0.96 mL/s, respectively. In Case II, the maximum velocity at the LVOT under a normal heartbeat was found to be 101 cm/s, whereas the average value was 73 cm/s. In Case III, the patient was also diagnosed with a congenital stenosis problem at the myocardial bridge (MCB) at the LAD. The measured blood flow rate at the MCB of the LAD for the three heartbeats in Case III was found to be 0.68, 1.08, and 1.14 mL/s.


The integration of patient-specific ECG signals and image-based CFD methods can clearly analyze hemodynamic information for patients during cardiac arrhythmia. The cardiac arrhythmia can reduce the blood flow in the coronary arteries.



Introduction
Cardiovascular disease (CVD) is a major cause of morbidity and mortality worldwide, with the lifetime risk exceeding 60% [1]. Globally, the number of deaths due to CVD has also increased by 12.5% in the past decade [2]. Coronary heart disease is the leading cause of deaths attributable to CVD in the United States (43.8%) [1]. In particular, coronary artery disease (CAD) is ranked second among the cardiac diseases. Generally, the most common clinical methods for CAD diagnosis are noninvasive coronary computed tomography (CT) angiography (CCTA) and invasive coronary angiography (ICA) [3]. Currently, ICA is considered the gold standard for coronary artery examination. However, ICA can only provide an anatomical evaluation of cardiovascular stenosis. It cannot clearly indicate the relationship between a coronary artery and myocardial ischemia. Szu-Hsien Chou and Chun-Jung Juan contributed to this work equally. Therefore, Pijls et al. developed an indicator for the functional assessment of the coronary artery hemodynamics, called fractional flow reserve (FFR) [4][5][6][7][8][9]. However, none of the aforementioned methods, including CCTA, ICA, and FFR, can evaluate the blood flow under cardiac arrhythmia conditions. At this stage, clinical diagnosis and treatment of problems, such as cardiac arrhythmia, can focus only on the transmission direction of the cardiac electrical signals. The effects of blood flow at the coronary arteries during cardiac arrhythmia require further understanding.
Cardiac arrhythmia, characterized by heart rhythm disruption, is caused by abnormal electrical signals in the heart. Premature ventricular contractions (PVCs), which are the most common form of cardiac arrhythmias, can be detected using the electrocardiogram (ECG) signal analysis. PVC is a common event occurring in a person of any age but more frequent in elderly people where the heart beat is initiated by the heart ventricles that are independent of the pace set by the sinoatrial node. Furthermore, atrial fibrillation (AF) is the most prevalent cardiac arrhythmia, the occurrence of which increases with age. The electrocardiogram (ECG), which measures the electric activity of the heart from the body surface and reveals the heart functions, is widely used a tool by cardiologists for determining the abnormalities of the human heart [10][11][12][13][14][15]. These heart activities in the generation of P-QRS-T waves featuring specific characteristics in terms of amplitude and duration. Irregular cardiac rhythm, also known as cardiac arrhythmia, refers to the different symptoms caused by the abnormalities of the cardiac electrical signal transmission. Accurate detection and classification of life-threatening ventricular arrhythmia episodes such as ventricular fibrillation and rapid ventricular tachycardia from ECG signals is a challenging issue. However, because ECG signals are unsteady and existing signal processing methods lead to relatively high computational burden, automated and accuracy diagnoses remain challenging.
Computation fluid dynamics (CFD) is a technique that uses numerical methods to solve the governing equations of fluid dynamics in a computer and predict the flow field [16][17][18][19]. The technique has been applied in the aerospace and automobile industry for more than 50 years. However, it has only been applied to the medical field in the last 10 years to evaluate the blood flow within the human circulatory system [16,17]. CFD has also been applied in simulating the blood flow in coronary arteries [18]. Miyoshi et al. [19] used CFD techniques to compute the fractional flow reserve from CCTA images (FFR CT ), a noninvasive examination. The results indicated that FFR CT has a considerable degree of accuracy [19][20][21][22][23]. However, how the different phases of patient cardiac arrhythmias influence coronary blood flow is not yet completely understood. Furthermore, Lin et al. [24] demonstrated that the CFD can be used to simulate coronary blood flow based on patient-specific CT images with different theoretical models of blood flow. Therefore, the aim of this work is to investigate the hemodynamics of coronary blood flow for patients during cardiac arrhythmia using the patient-specific ECG signals and the image-based CFD method.

Materials and Methods
Three cardiac arrhythmia patients were examined. Two of them were diagnosed with premature ventricular contractions (PVCs) (i.e. PVC-type) cardiac arrhythmia and one was diagnosed with atrial fibrillation (AF) (i.e. AF-type) cardiac arrhythmia. All patient ECG signals were simultaneously recorded during CT scan acquisition. Meanwhile, recorded ECG signals of three cardiac arrhythmia patients were contained with three to four consecutive heartbeats for analysis.

CT Imaging and Geometry Reconstruction
The CT images of the patients were obtained from a 640-sliced (320-detector) CT scanner (Aquilion ONE, Toshiba Medical Systems Corporation). The CT scanner used a slice thickness of 0.5 mm, interval of 0.25 mm, rotation speed of 0.35 s/rotation, and a matrix size of 512 × 512. Three-dimensional reconstructions of the obtained coronary artery images were created using the Toshiba Vitrea workstation. The processed 3D geometrical models were output into the STL (i.e. StereoLithography) file format. The coronary artery surface geometries were imported into the ANSYS ICEM-CFD.
According to Kwon et al. [25], the computation time for the coronary artery flow simulation can be greatly reduced when the computation only includes the left and right coronary arteries but excludes the aorta. The difference in the simulated FFR and those of the model including the aorta is < 2%. Therefore, the simplified model in this study also excluded the aorta for better computation time efficiency.

Volume Meshing
Volume meshes of coronary arteries were generated with the meshing module in the ANSYS. The volume mesh was created by applying the unstructured tetrahedral mesh in this study. We performed the grid independence test using the laminar flow model and the pressure gradient difference of the right coronary artery for Case I was displayed in Table 1. Considering the balance between the numerical accuracy and the computational time, the coronary artery flow model in this study has more than 2.5 million computation grids for each case.

Boundary Conditions and Parameter Settings
In this study, the inlet velocity transient profile of the coronary artery in simulation was used with the measured velocity profile of each patient LVOT by the Doppler echocardiography. However, because the simulation of this study excluded the aorta, it is necessary to redistribute 4% of the cardiac output (CO, estimated from the output from the left ventricle to the aorta) to the coronary arteries according to the body circulation [18,26]. The redistributed flow of the left and right coronary arteries was then directly set at the entry point of the corresponding arteries [22]. The inlet velocity waveforms were determined according to the transformed periodic wave patterns from the electrocardiographic waveform obtained during the patient's CCTA examination ( Fig. 1). The periodic wave patterns were transformed into a function using MATLAB (Fig. 1) and were set at the entrance of the coronary arteries on both sides.
In order to guarantee that the boundary conditions do not affect the blood flow in the region of interest, a cylindrical extension was added to the vessel inflow and outflow (Fig. 1c). The length ( L e ) added to the vessel inflow and outflow was estimated using a semi-empirical law valid for laminar flow [27]: where D is the diameter of the vessel and Re is the dimensionless Reynolds number.
The transient waveforms of the aortic pressure of the patients were imposed at as the outlet boundary condition.
(1) L e D ≈ 0.06Re,  In general, the flow in a circular pipe is laminar for Re ≤ 2,300, turbulent for Re ≥ 4000, and transitional in between. When the Reynolds Number is less than 2300, the effect of the viscous force to the blood flow is greater than that of inertia. The Reynolds numbers in this study cases are less than 1000 that the laminar flow is appropriate. In other words, the Reynolds number is small enough that laminar flow description is valid. In this case, at small or moderate Reynolds numbers, the viscous forces are large enough to suppress these fluctuations and to keep the fluid in line. The blood flow is laminar when the Reynolds number is small. Thus, the laminar flow model was used. Furthermore, the FFR CTS [24] values obtained from a laminar flow simulation and cardiac catheterization were relatively closer to each other than the results of a turbulent flow simulation [28][29][30], even for stenosisaffected blood vessels [24]. This may be because most of blood flow within the entire blood flow field is considered laminar. The turbulence mode may underestimate or overestimate the simulated values. Therefore, laminar flow was set as the blood flow state in this study.
In this study, the density and dynamic viscosity of blood were 1060 kg/m 3 and 4 cP, respectively. On the arterial wall we imposed a no-slip boundary condition between the blood flow and the arterial wall. Furthermore, Siogkas et al. [31] reported the difference in FFR value between the rigid and elastic wall models was about 0.44% for the coronary blood flow simulation. Hence, we also assumed the arterial walls to be rigid for numerical simulations. The initial velocity and pressure were both set as 0 cm/s, and 760 mmHg (i.e., 1 atm) was taken as the reference pressure.

Difference between using echocardiographic data obtained for LVOT and three-lead ECG for entrance waveform setting
In our previous research, we used the flow velocity waveform of the LVOT obtained from the patient ECG signals for the inlet flow velocity [24]. The focus of this study was the coronary blood flow condition during cardiac arrhythmia. Before further discussing the result of the blood flow condition during cardiac arrhythmia, we first compare the two simulation methods for any differences in the results for the first participant (Case I). Figure 2 illustrates the LVOT waveform obtained from the echocardiogram and the three-lead ECG waveform for Case I. Using the patient's CCTA image from the same heartbeat, a set of 3D models for the right coronary artery (RCA) was constructed. The same cardiac rhythm, entry flow velocity, exit pressure, and vascular wall setting were modeled. In the simulation model, the periodic waveform (i.e., LVOT) at the inlet point was used as shown in Fig. 2a and the cardiac rhythm (regular cardiac rhythm) D (i.e., ECG wave form, as displayed in Fig. 2b) was simulated in the blood flow pattern.
There was a segment of RCA affected by moderate stenosis, where S2 enters S3 for Case I. The CCTA report indicated that it was 20-30% narrower than normal arteries. By using this point of narrowing as the center, the pressure within the artery at a point 20 mm to the front end and 20 mm to the back end were measured. The ratio of these two pressures was then calculated according to the definition of FFR (i.e., the pressure at point B was divided by the pressure at point A). The ratio using the waveform obtained from LVOT as the entry flow velocity waveform was 0.872, whereas that of the ECG was 0.867. The difference was of only 0.005, illustrating the correlation between the waveforms obtained from LVOT and ECG. Hence, the simulations thereafter used the ECG waveform obtained from CCTA scanning to model the blood flow condition during the occurrence of cardiac arrhythmia.

Case I: blood flow and heartbeat simulation for PVC-type cardiac arrhythmia
Case I was of PVC-type cardiac arrhythmia. The highest flow velocity for his LVOT under regular cardiac rhythm was 88 cm/s, whereas the average value was 56 cm/s. The blood pressure during cardiac CT was found to be 109/66 mmHg. Five heartbeats were scanned. The first four were heartbeat A at 109 bpm, heartbeat B at 37 bpm, heartbeat C at 52 bpm, and heartbeat D at 52 bpm. These were analyzed as in Fig. 3. The difference in patterns can be observed from the model by comparing the blood flow at the RCA. Blood was pumped out after completing the ventricular contraction process, as indicated by heartbeat A, which should be followed by the start of the diastolic phase under normal conditions. The blood was also expected to enter the coronary arteries from the openings at the aortic root on both sides. However, because the occurrence of PVC that led to an early contraction of the ventricles, the blood at the aortic root was pushed forward again by the blood squeezed out from the left ventricle, as indicated by heartbeat B, before successfully entering the coronary arteries. This led to a decrease in blood flow. The blood flow measured at approximately 2.0 cm from the RCA entry point was 0.66 mL/s. Heartbeat B started at the point where the PVC occurred and early electrical signals were released. Because of the prior shortened diastolic phase, as indicated by heartbeat A, the blood strongly pumped out with a higher R wave directly from the left ventricle. Although the CO at this moment was relatively low, there was more time for the blood to flow into the coronary arteries, because the heartbeat that follows the early ventricular contraction will have a compensation delay to return to a normal cardiac rhythm, leading to a relatively longer diastolic period.
The blood flow measured at approximately 2.0 cm from the same opening of the RCA inlet was approximately 0.96 mL/s. Heartbeat C represented a relatively regular cardiac rhythm from the same patient. Although this cardiac rhythm was normal, owing to the extended diastolic period from the previous heartbeats, the amount of blood replenished to the ventricles was relatively larger than normal. Hence, a larger amount of blood was also pumped out during this particular cardiac contraction. The blood flow measured at approximately 2.0 cm from the RCA entry point was found to be approximately 0.89 mL/s. Heartbeat D also represented a relatively regular cardiac rhythm from the same patient.
The diastolic period before this heartbeat was also normal. Therefore, it was expected that the amount of blood replenished to the ventricles was also at the normal value. Hence, heartbeat D was considered the normal heartbeat for the patient. The blood flow at approximately 2.0 cm from the RCA entry point was approximately 0.81 mL/s.

Case II: Blood flow and heartbeat simulation for PVC-type cardiac arrhythmia
Case II was also of PVC-type cardiac arrhythmia in addition to a history of hypertension and high cholesterol for nearly 10 years. The reason for undergoing the CCTA examination was occasional chest tightness for approximately 1 year. The highest flow rate at the LVOT under a normal heartbeat was found to be 101 cm/s, whereas the average value was 73 cm/s ( Table 2). The blood pressure was found to be 137/87 during the cardiac CT. Four heartbeats were extracted from the scanning (Fig. 4). Heartbeats A, B, C, and D were found to be 44, 132, 46, and 66 bpm, respectively. In Case II, it can be observed from the CT images that there is a segment of blood vessel at the middle part of the LAD that enters the cardiac muscle, which is known as the myocardial bridge (MCB). The MCB is a common congenital anatomical deformity. From the viewpoint of the cardiac structure, the coronary arteries for normal people should be above the cardiac muscles. However, the coronary arteries for patients suffering from an MCB are within the cardiac muscles. The coronary arteries are therefore squeezed by the cardiac muscles during the systolic period and become narrower. Most patients with MCB typically show no symptoms. Some patients exhibit cardiac hypoxia during vigorous exercise. However, as the levels of hypertension, cholesterol, and blood sugar increase with age, the MCB can become increasingly detrimental.  The simulation result indicated that the LAD at the MCB had not considerably narrowed during the first extended heartbeat A. This may be a result of the relatively long diastolic period experienced by this heartbeat, which allowed the cardiac muscle to have a longer period to relax. With a larger R wave rhythm, the flow rate at that point was found to be 1.37 mL/s (Fig. 4a) ( Table 2). The LAD at the MCB was found to be considerably narrowed during the second shortened heartbeat B. This may be a result of the heartbeat image being reconstructed at the end of the contraction period and the heartbeat not entering the diastolic period. The flow rate at the point of the MCB was found to be 0.88 mL/s (Fig. 4b), as displayed in Table 2. The MCB was not found to become considerably narrower during the third heartbeat C. The reason may be similar to that for the first heartbeat. The blood flow measured at that point was found to be 1.35 mL/s (Fig. 4c). The fourth heartbeat D was found to be a normal heartbeat, where the point flow velocity at the MCB was found to be 1.21 mL/s (Fig. 4d), as shown in Table 2. Small calcification was found only at three places in the RCA of this patient. The diameter of the blood vessels was not affected, and no notable narrowing occurred. The blood flow rates at the largest calcification during the corresponding to four heartbeats were found to be 0.72, 0.52, 0.69, and 0.58 mL/s. Based on the aforementioned results, a considerable blood flow decrease was noted in heartbeat B because of the shortened diastolic time and smaller R wave. Using heartbeat D as the reference average heartbeat, the blood flow rate of heartbeat B was only 72.73% of the flow under the average heartbeat.

Case III: Blood flow and heartbeat simulation for AF-type cardiac arrhythmia
Case III was of AF-type cardiac arrhythmia, with a 20-30% stenosis at the MCB in the CT cardiac examination (Fig. 5).
The largest flow velocity at the LVOT under a normal cardiac rhythm was 93 cm/s, whereas the average value was 72 cm/s. The blood pressure during cardiac CT was 130/84 mmHg. Four heartbeats were extracted from the scanning. Heartbeats A, B, C, and D were found to be 80, 46, 67, and 68 bpm, respectively. The first three heartbeats were selected for further discussion focusing on the MCB of the LAD. The blood flow measured at the MCB of the LAD for the three heartbeats in Case III was found to be 0.68, 1.08, and 1.14 mL/s. The results showed a good fit with the Case I and II simulation results, where a decrease in blood flow was also observed in the shortened cardiac rhythm. Notably, a greater decrease in blood flow was found to be possible in Case III, as the occurrence of AF could lead to a longer MCB. With heartbeat C as the reference average heartbeat, the blood flow for heartbeat A was only 59.65% of the normal value. Thus, AF can lead to a greater decrease in blood flow at the coronary arteries when compared with PVC. Given that patients in both Cases II and III had a narrowing point at the MCB of the LAD,  27.27% in Case II when PVC-type cardiac arrhythmia occurred, whereas in Case III, it decreased by 40.35% when AF-type cardiac arrhythmia occurred. This indicated that AF-type malignant cardiac arrhythmia had a stronger effect on the blood flow rate than did the PVC type.

Discussion
The ECG is noninvasive and widely used in clinical practice. As a recording of the body surface electrical activity, it provides information about heart rhythm abnormalities and helps detect CADs. However, visual inspection of the ECG provides discrete clinically interpreted features which cannot reveal the coronary blood flow information of ECG abnormalities and morphologies. Thus, this is why CFD method are required, as it can make sense of coronary blood flow and provide the ECG-based blood flow of patients during cardiac arrhythmia. Treatment for cardiac arrhythmia differs according to the type of cardiac arrhythmia and cause of the associated cardiac diseases. From the three simulation results mentioned previously, this study found that blood flow decreased regardless of whether the heartbeat was (1) an incomplete or compensatory-delayed because of early PVC or (2) shortened or lengthened because of AF where incorrect electric signals were released. In addition to studying more about the blood flow changes caused by PVC-and AF-type cardiac arrhythmia, this study aimed to extend its application to understanding other types of malignant cardiac arrhythmia, such as those related to left bundle branch block and paroxysmal supraventricular tachycardia. The preliminary results from our study demonstrated the capability of integrating electrocardiogram and image-based CFD method to demonstrate the coronary blood flow in patients during cardiac arrhythmia.
The feasibility of the integrating of the patient-specific image-based CFD method and ECG signals was tested in three arrhythmia patients. In addition, there were some limitations in the current study. First, we did not consider the motion of the coronary arteries during the cardiac cycle. In reality, the coronary arteries move significantly to contract and relax during the cardiac cycle. Second, we assumed the rigid coronary arteries for the whole computational model although we know that coronary arteries are elastic. Third, this study had a limited number of subjects recruited and should be extended by further experience. A follow-up study with a larger number of subjects may provide a better valid assessment of coronary blood flow in patients with cardiac arrhythmia.

Conclusions
In this study, we proposed a novel approach of integrating patient-specific ECG signals into CT image-based CFD simulation to analyze the blood flow in coronary arteries in patients with two types (two having premature ventricular contraction and one having atrial fibrillation) of cardiac arrhythmias. This study successfully demonstrated that cardiac arrhythmia (i.e., irregular heartbeats) can reduce the blood flow rate in the coronary arteries, especially in the presence of the AF-type cardiac arrhythmia patient.
Author contributions All authors were actively involved in the conception, implementation, and data collection for the study. All authors met to discuss the findings of the study with respect to manuscript preparation prior to it being written. SHC was the principal author of the manuscript's first draft, produced the figures and tables, and edited subsequent revisions. CJJ clarified the idea, critically and strengthened the method of this study and edited the entire revised manuscript in terms of scientific and English writing. TCS was co-author of the manuscript's first draft. He also was the principal author involved in data analysis and helped produce many of the figures and tables. SHC, KYL, CYH, CJJ, and ZYC collected the clinical data and analyzed. TCS reviewed and edited the manuscript first draft and made significant editorial contributions to the final version which has been submitted.

Compliance with ethical standards
Conflict of interest None of the listed authors has any financial or personal relationships with other people or organizations that could potentially/inappropriately influence their work and conclusions with respect to this study or the submitted manuscript.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/ by/4.0/.