Hemodynamic diagnostics of epicardial coronary stenoses: in-vitro experimental and computational study

Background The severity of epicardial coronary stenosis can be assessed by invasive measurements of trans-stenotic pressure drop and flow. A pressure or flow sensor-tipped guidewire inserted across the coronary stenosis causes an overestimation in true trans-stenotic pressure drop and reduction in coronary flow. This may mask the true severity of coronary stenosis. In order to unmask the true severity of epicardial stenosis, we evaluate a diagnostic parameter, which is obtained from fundamental fluid dynamics principles. This experimental and numerical study focuses on the characterization of the diagnostic parameter, pressure drop coefficient, and also evaluates the pressure recovery downstream of stenoses. Methods Three models of coronary stenosis namely, moderate, intermediate and severe stenosis, were manufactured and tested in the in-vitro set-up simulating the epicardial coronary network. The trans-stenotic pressure drop and flow distal to stenosis models were measured by non-invasive method, using external pressure and flow sensors, and by invasive method, following guidewire insertion across the stenosis. The viscous and momentum-change components of the pressure drop for various flow rates were evaluated from quadratic relation between pressure drop and flow. Finally, the pressure drop coefficient (CDPe) was calculated as the ratio of pressure drop and distal dynamic pressure. The pressure recovery factor (η) was calculated as the ratio of pressure recovery coefficient and the area blockage. Results The mean pressure drop-flow characteristics before and during guidewire insertion indicated that increasing stenosis causes a shift in dominance from viscous pressure to momentum forces. However, for intermediate (~80%) area stenosis, which is between moderate (~65%) and severe (~90%) area stenoses, both losses were similar in magnitude. Therefore, guidewire insertion plays a critical role in evaluating the hemodynamic severity of coronary stenosis. More importantly, mean CDPe increased (17 ± 3.3 to 287 ± 52, n = 3, p < 0.01) and mean η decreased (0.54 ± 0.04 to 0.37 ± 0.05, p < 0.01) from moderate to severe stenosis during guidewire insertion. Conclusion The wide range of CDPe is not affected that much by the presence of guidewire. CDPe can be used in clinical practice to evaluate the true severity of coronary stenosis due to its significant difference between values measured at moderate and severe stenoses.


Background
In current clinical practice, the ischemic severity of epicardial coronary stenosis can be detected by functional and anatomical measurements. The anatomical or geometric measurements, performed by imaging methods (e.g., contrast angiography, angioscopy and intravascular ultrasound), may not reveal the true severity of coronary stenosis and may lead to inappropriate clinical decisions [1]. On the other hand, functional or hemodynamic measurements of mean pressure drop across the stenosis and distal coronary flow prove to be more useful for the long term success of balloon angioplasty with or without stent placement procedure [2]. Hence, the development of simple and effective clinical diagnostic methods utilizing these functional measurements has attracted the interest of many researchers. For this purpose, the intravascular clinical procedure is performed before the angioplasty by surgically placing a 2-2.3 mm diameter guide-catheter up to the coronary ostium. The 0.35 mm diameter pressure or Doppler flow sensor-tipped guidewire is advanced through this guide-catheter and placed across the stenosed coronary artery to measure distal pressure or flow [3]. Using this invasive pressure and flow measurement, various diagnostic parameters have been developed in current clinical practice e.g. coronary flow reserve (CFR: the ratio of hyperemic flow to the basal flow) [4], trans-stenotic pressure drop [5], fractional flow reserve (FFR: the ratio of distal recovered pressure to the aortic pressure at hyperemia) [6].

Uncertainties in pressure-flow diagnostics
To distinguish epicardial stenosis severity from microvascular dysfunction, simultaneous measurements [7] of the coronary flow and pressure have been recently recommended [8,9]. For this purpose, both pressure and flow sensors are simultaneously inserted into the coronary vessel. However, the insertion of guidewires across the stenosis causes a "tight fit" in the lumen of the stenosed artery and is associated with additional flow reduction and an increase in mean pressure drop [2,5] that may result in an incorrect overestimation [10,11] of hemodynamic parameters [12] and may lead to improper selection of therapeutic procedures [13,14]. Possible risks and resulting consequences to the patients because of the uncertainties in pressure drop and flow measurements due to the size and number of guidewires inside the diseased vessel were recently reported. In earlier clinical study, clinicians observed that the diastolic pressure drop measured by a balloon catheter (1.98 mm size) was 34% more than that measured by a single guidewire (0.45 mm size) [14], typically used in present clinical settings. In a study of 65 patients, Verberne et al [11] reported a 25% increase in hyperemic pressure drop and a 20% reduction in hyperaemic flow with the use of two guidewires (separate pressure-and flow sensor-tips on different guidewires; a procedure not commonly used) instead of a newly developed single guidewire having both pressure and flow sensors on its tip. Further data analysis indicated 30%, 128% and 325% increase in the misdiagnosis (false positive cases) of the true ischemic severity of coronary stenosis due to wrong measurements of CFR, FFR and hyperemic stenosis resistance (hSRv: the ratio of hyperemic pressure drop and distal blood flow velocity), respectively [11]. The uncertainty in deciding the limiting values of CFR, FFR and pressure drop measured by different sizes of guidewires, has been experimentally [15] and numerically evaluated and discussed in detail in our previous studies [16,17].
We hypothesized that the hemodynamic parameter consisting of simultaneous measurements of pressure-flow integrated with anatomical details can overcome the disadvantages of conventional diagnostic parameters. Recently, our group has developed an alternative diagnostic parameter: Lesion flow coefficient (LFC: ratio of nondimensional pressure drop at very high flow rate and at patho-physiological hyperemic flow rate) [17] to determine the true severity of epicardial coronary stenosis. Accordingly, in this study, we sought to investigate the behaviour of the non-dimensional hemodynamic or functional diagnostic parameter, pressure drop coefficient (CDP e ), over a wide range of coronary flows simulated in an in-vitro experimental set up. The CDP e utilizes simultaneous pressure-flow measurement across the stenosis, while evaluating the downstream pressure recovery factor (η) that includes simultaneous pressure-flow measurement and anatomical details.
The specific aims of this study are: 1) to compare the pressure drop-flow relation for pulsatile flow in the presence and absence of single guidewire, 2) To evaluate the behaviour of hemodynamic parameter, CDP e while evaluating η for different degrees of stenoses under variable pressureflow scenarios. where, is spatial and temporal mean blood velocity at the minimal vessel area. Therefore, we define the pressure recovery factor (η) to be the ratio of experimentally measured CPR to the quantity CPR ∞ , that represents the pressure recovery coefficient at very high (limiting) Reynolds number (Re) flow, i.e.,

Formulation
where CPR ∞ is expressed as We define Δp r = p r -p m and α is the ratio of the flow crosssectional area of the distal vessel, A r , to the minimum vessel area at the site of stenosis, A m , prior to guidewire insertion.
where, A is the cross-sectional area, subscript 'r' represents the distal vessel; 'm' represents the vessel at the site of stenosis with minimum flow area; 'e' represents the proximal healthy vessel. For simplicity, it is assumed that A r = A e i.e. proximal and distal vessel diameters are the same. Simi-larly, we define α i to be the ratio of the flow cross-sectional area of the distal vessel (A r -A i ) to the minimum vessel area at the site of stenosis (A m -A i ), during guidewire insertion.
where, subscript 'i' represents the guidewire. Further, the % occlusion of the vessel or % area stenosis is defined as follows: Here, the pressure recovery factor (η) measures the ability of distal vessel to recover the pressure downstream to the stenosis [19] and have wide application in nozzle flow fluid dynamics.

Experimental method
Simultaneously measured pressure and flow data for each stenosis test sections were compared with corresponding numerical calculations. The and were analysed for each stenosis section by comparing pressure-flow characteristics for both 'before' and 'during' guidewire insertion under pulsatile flow conditions. Finally, fluid dynamic parameters: overall pressure drop coefficient (CDP e ) and the pressure recovery factor (η) were evaluated.

Coronary Stenosis Test Sections
We created three model test sections to simulate different severity of epicardial focal coronary stenosis. These test sections are manufactured with optical grade Lexan material for following three severity conditions: moderate, intermediate and severe stenosis, having % area occlusion of 65%, 80% and 89%, respectively. Epicardial focal coronary stenosis consists of three distinct sections, converging or constricting, throat, and diverging section as shown in Fig. 1A, B [13,20]. The above manufactured test sections are as per the characteristic dimensions reported by Wilson et al. ([13], Table 1). To measure the axial pressure and time averaged pressure recovery ( ) in the distal part of the test section, 0.3 mm diameter radial pressure ports were drilled in the axial direction at an interval of ~5 mm as shown in Fig. 1. The proximal pressure ( ) was measured by pressure port located just proximal to the converging section while the throat pressure ( ) was measured by pressure port located at the middle of throat section. The mean trans-stenotic pressure drop in equation 1 can be calculated as: and the mean distal recov- ered pressure in equation 2 can be calculated as: . Details of these models are reported in our recent publication [15]. Anatomy of stenoses have curvatures and transition from different sections of the stenosis are gradual; accordingly, three custom mill bits are used to sequentially drill, first, the throat section fol-lowed by the converging section from one end and diverging section from the other end.

Experimental Setup
The diastole-dominated coronary arterial flow waveform was generated by compliance-resistance method as  explained in our previous in-vitro study [15]. The schematic diagram of experimental set up is shown in Fig. 2.
The basic pulsatile flow waveform (T = 0.8 sec, = 350 ml/min) was generated by the pulsatile pump (Harvard Apparatus, MA), which is similar to the aortic flow waveform. This flow was then bifurcated into two flow conduits: A and B in order to divert very high flow rate from the pulsatile pump. The compliances (C1 and C2), and resistances (R1 and R2) were connected and adjusted within the conduits to generate the physiological coronary arterial flow waveform (Fig. 3). The mean flow rate was increased from basal (~50 ml/min) to hyperemic flow (moderate: 180 ml/min, intermediate: 165 ml/min, severe: 115 ml/min) in four steps. The pulsatile pressure and flow were measured simultaneously with the help of a reference trigger square pulse. Figure 4 explains the method of simultaneous pressure-flow measurement. A typical left anterior descending (LAD) coronary flow pulse was reported by many researchers [21,22]. Generally, the systolic dominance in the aortic flow profile is reversed to diastolic dominance in the normal to moderately stenosed vessels. The experimentally obtained flow pulse ( Fig.  3) was almost similar to the pulse velocity measured by a Doppler catheter in the LAD of patients undergoing angioplasty [21,22]. the stenotic test section to measure trans-stenotic pressure and flow.

Pressure and Flow Measurements
The axial pressure was measured by a DSA3207 digital sensor array (Scanivalve corp., WA) with 12 ms time interval. The pulsatile flow was measured with an in-line Doppler flow sensor (Transonic Inc., NY) with a time interval of 1 ms. Simultaneous pressure-flow values were measured for the following two conditions: 1) Before guidewire insertion: This method can be considered as "non-invasive or patho-physiological" measurements, since pressure was measured without guidewire insertion across the stenotic models.
2) During guidewire insertion: This method can be considered as "invasive" measurements, since pressure was measured after guidewire insertion across the stenotic model. The guidewire was connected to the ComboMap system (Volcano Therapeutic, CA) for continuous pressure recording during its insertion across the stenotic models.

Blood Analog Fluid (BAF)
The experiments were performed with a fluid exhibiting the shear thinning, non-Newtonian viscous property of blood. The BAF was prepared by mixing 65% water, 35% glycerine, and 0.02% Xanthum gum (by weight) [23]. The viscosity was measured by concentric cylinder viscometer (DV-II+ PRO Digital Viscometer, Brookfield, MA). The experimental data, showing shear-thinning behaviour of BAF, were fitted to Carreau model (Eq. in Fig 5), where the Carreau coefficients are as follows: zero shear rate viscosity, μ 0 = 55 cP, infinite shear rate viscosity, μ ∞ = 3.39 cP, time constant, λ = 9.56 s, and power index, n = 0.2. The comparison of BAF viscosity and real blood viscosity [24] is provided in Fig. 5. The density of BAF was measured as 1.05 g/cm 3 . The comparison of results between BAF and blood viscosities is discussed in Additional file 1.  Table 1. Two numerical models were generated for each stenotic test section; namely, before and during guidewire insertion. Mesh diagrams 'during' and 'before' guidewire insertion are shown in Fig. 6A and Fig. 6B. It is assumed that: (a) the arterial wall has a smooth, rigid and round concentric shape; (b) the stenosis geometry remains unchanged during basal and hyperemic flow due to failure of the flow-dependent dilation mechanisms in atherosclerotic coronary arteries; and (c) the rigid arterial wall has insignificant effect of pressure pulsation on its dimensions.

Boundary and Initial Conditions
The 'no slip' boundary condition, u i = 0, was specified on the arterial wall and on the guidewire wall (in case of guidewire inserted conditions). A stress free boundary condition (σ i = 0) was applied at the flow outlet. The symmetric boundary condition (u r = 0) was applied at the central axis in the case of 'before' guidewire insertion. The coronary flow waveform, Q(t), used in the computational analyses, was obtained from the corresponding in-vitro experimental data. The spatially averaged velocity along the vessel cross-section, , needed for the computational analysis, was obtained from the mass balance equation: = Q(t)/A e . The with analogous poiseuille flow relation was applied at 3 cm (at the point where l/d = 10) proximal to the converging section of the numerical stenotic models in case of 'before' guidewire insertion case. Thus, the inlet boundary conditions are [25,26]: Simultaneous pressure-flow measurement  Because of long entry length, the spatial velocity profiles at any instance of time adjusted to a profile that is consistent for non-Newtonian BAF fluid.
The temporal mean pressures were calculated along the arterial wall to find the mean pressure drop, and mean pressure recovery, , where as , and represent the mean pressure measured proximal, distal (at the site of pressure recovery) and at the throat region of the stenosis. The experimentally observed non-Newtonian viscosity and density of BAF were used for this calculation.

Solution Strategy
The unsteady problem was solved by finite element method utilizing the Galerkin scheme (FIDAP, ANSYS, Inc., NH). The pressure was discredited by mixed formulation and approximated as discontinuous across the element boundaries for this incompressible flow problem. The 2 nd order implicit trapezoidal time integration scheme was used to control local truncation error. A successive substitution type of fully coupled iterative solver was used to obtain the solution at each time step of this non-linear time-dependent problem. This method solves the linearized system of governing equations by direct Gaussian elimination approach. The convergence criteria for the velocity and residual vector were tightened to a value of 0.01% which is two orders of magnitude lower than the recommended value [27]. In order to optimize the convergence time at each time step, a relaxation factor was used (0.5). To achieve smooth converged results, a relatively small value of streamline upwinding (a value less than 0.45) that adds numerical diffusion along the flow direction was used to calculate primary flow varia-

Mesh refinement study
A series of meshes were created, with each one 20% higher in density than the previous one, in order to check the overall convergence accuracy of the numerical calculations [25][26][27]. When the improvement with the refined mesh was less than 1% in velocity vectors, wall shear stress, and pressure, the numerical calculation was considered to be converged. For pulsatile flow calculations, the convergence is not only dependent on the mesh resolution but also on the selection of the time step. Depending on the velocity pulse shape and stenosis geometry, the calculation time steps varied between 1 × 10 -4 to 1 × 10 -5 sec. These calculations were started from a time where the velocities were near zero in order to make the stiffness matrix well-balanced and stable. Considering the refinements in time-steps, and appropriate starting point of numerical calculations, one may not need much finer mesh for an unsteady calculation [25][26][27].

Trans-stenotic mean axial pressure drop
The experimentally measured and numerically calculated temporal mean axial pressure data were analyzed to deter-mine the overall mean trans-stenotic pressure drop ( ), mean distal pressure recovery ( ). The and were quantified based on the experimentally measured mean pressure drop and mean flow rate values having quadratic relation. Figure 7A shows the trans-stenotic mean axial pressure as a function of axial distance before guidewire insertion.  Table 2. The axial pressure drop profiles for above flow rates are plotted in Fig. 7B. The CPR values increased from 0.09 to 0.39 as Re e increased from 109 to 367; however, they were less than the corresponding values for 'before guidewire insertion'. Moreover, for all flow rates, the mean distal pressure recovered within a shorter axial distance as compared to 'before guidewire insertion'. This

Moderate (65% area) stenosis
Numerical models Figure 6 Numerical models. A -Numerical model for moderate stenosis (during guidewire insertion). B -Numerical model for moderate stenosis (before guidewire insertion).
Mean pressure drop vs. Axial distance for moderate stenosis Intermediate (80% area) stenosis Before guidewire insertion (Fig. 8A), sudden pressure drop in the converging and throat region was observed with increased overall pressure drop as compared with moderate stenosis. The and are summarized in Mean pressure drop vs. Axial distance for intermediate stenosis Severe (89% area) stenosis Before guidewire insertion, the overall pressure drop ( ) was much higher with sharp pressure drop in the converging section as compared with that for moderate and intermediate stenosis (Fig. 9A), primarily because of an area constriction effect. The distal pressure recovery was highly unstable because of enhanced and organized vortical cells. The and are summarized in  Figure 9B shows the axial pressure drop profile during guidewire insertion for 34, 61, and 81 ml/min. The guidewire insertion significantly reduced the blood flow with increase in . The CPR reduced from 0.46 to 0.30 for Re e range of 68-161. The pressure recovery values were observed at 17 and 13 mm distal to the stenosis before and during guidewire insertion, respectively. This type of significant stenosis is also known to cause ischemia in the subendocardium and coincides with symptomatic angina [13].
This study also suggests that the tip of the guidewire sensor should be positioned at sufficient distance distal to the stenosis. The pressure sensor tip should be positioned after the pressure recovery in order to avoid inaccuracy in distal pressure measurement. From axial pressure drop plots, shown in Figs. 7, 8 and 9, it is evident that the pressure tip should be at least 14 mm distal to the stenosis.

Mean pressure drop and flow relation
The time averaged transient pressure drop-flow ( ) characteristic for each stenosis geometry is compared and validated numerically (Fig. 10). The experimentally measured time averaged trans-stenotic pressure drop, pressure recovery, flow rate and Reynolds number are summarized in Table 2. The corresponding quadratic correlation, as provided in method section: formulation, is presented in Fig 10 for each stenosis geometry for pulsatile flow. Then the overall and for each stenosis model can be compared based on the coefficients 'A' and 'B'. These coefficients are summarized in Table 3.
For all stenotic models, the guidewire insertion increased the trans-stenotic pressure drop. The mean pressure drop for moderate stenosis (65% area stenosis) was governed by the viscous component, , which was further augmented by guidewire insertion (Fig. 10A; coefficient 'B' in Table 3). Due to dominance of , increased noticeably during guidewire insertion for severe stenosis (89% area stenosis) as shown in Fig. 10C. Before guidewire inserted condition, coefficient 'A' for severe stenosis increased by factors of 16 (= 32 × 10 -4 /2.01 × 10 -4 ) and 2.4 (= 32 × 10 -4 /13.5 × 10 -4 ), as shown in Table 3 could not be clearly distinguished (Fig. 10B). Thus, the guidewire insertion acts as an enhancer of viscous-losses in moderate stenosis, whereas it is the momentum-change pressure loss that dominates the severe stenosis and this distinction is unpredictable for intermediate stenosis.
In addition, comparison of coefficients 'A' and 'B' for reference steady flow (not shown in the figure; shown in Table 3) and pulsatile flow showed that the pulsatile flow increased momentum-change effects significantly. For a selected mean flow rate, pressure drop for steady flow was lower than the time-averaged pulsatile flow cases. The comparison of steady and pulsatile flow coefficients shows that pulsatile flow increased considerably than . The pulsatile flow enhances momentum-change pressure loss due to the significant effect of convective term: .
Mean pressure drop vs. Axial distance for severe stenosis blood flow resistance. Figure 11 shows mean CDP e vs. Re e characteristic for each stenosis and flow condition before and during guidewire insertion. The CDP e was primarily dependent on stenosis geometry, shape of flow pulse, presence of guidewire and flow rate. The mean CDP e increased from 17 ± 3.3 for moderate stenosis to 287 ± 52 for severe stenosis (n = 3, p < 0.01 using unpaired t-test). For a given stenosis, values of CDP e : decreased with increase in Re e ; increased after guidewire insertion; and are higher for pulsatile flow than those for steady flow cases.
Mean pressure drop vs. Mean flow rate  profile in conjunction with the area measurement by diagnostic imaging techniques (e.g. IVUS, angiograph) in clinical settings.

Discussion
The present study examines the pressure drop-flow relation for the range of epicardial coronary stenoses with quantification of and before and during guidewire insertion. More importantly, this research provides further hemodynamic insights on clinical data [13,28] and corroborates our earlier numerical studies [25,26]. The effect of guidewire eccentricity in relation to lumen diameter has also been performed, which is discussed at the end of this section.
The large difference between pressure drop values for the same mean flow implies that the contribution of viscous and momentum-change related pressure drops are different for steady and pulsatile flows. The pulsatile flows, hav-ing different values, show diverse pressure drop values. Hence, it is necessary to distinguish the steady and pulsatility flow effects on viscous-and momentum-change pressure losses. Further, the effect of physiological parameters such as pulse shape, heart rate, and systole to diastole ratio are needed to be studied in the future.
Depending upon blood disorder, the effect of blood viscosity on diagnostic parameters should be analyzed since it may vary from patient to patient. The blood viscosity may affect the diagnostic procedure for intermediate and moderate stenosis as pressure drop for these geometries are dominated by viscous losses. In contrast, for the severe stenosis cases only momentum-change effects are significant and thus, blood viscosity may not have significant effect on diagnostic parameters.
A new diagnostic parameter 'lesion flow coefficient' (LFC; see Additional file 2) was introduced to assess the severity of epicardial stenosis [15,17]. LFC is the ratio of square root of CDP 8 2 where, α is area ratio term obtained from equation 5 [17]. Thus, LFC is a function of flow rate, the severity of stenosis, guidewire size and the pulse shape. The nature of viscous-and momentum changerelated pressure losses for CDP e will help in understanding the characteristic of the LFC. The detailed axial pressure drop profiles for each stenosis severity obtained from this study provide additional insights for determining viscous-and momentum change-related pressure losses.
A wide range of CDP e reveals that it can itself be used as a diagnostic index and thus, it has its own clinical importance. Since calculation of CDP e requires simultaneous pressure and flow measurement, it may be used to detect epicardial stenosis under presence of microvascular disease. The wide and non-overlapping range of CDP e , e.g. 14-22 for moderate stenosis, 87-132 for intermediate stenosis and 251-375 for severe stenosis may be useful for the estimation of severity of coronary stenosis. This invitro data will be helpful to design future pre-clinical and clinical trials.
The distal pressure recovery and η signifies the influence of pressure recovery, distal vascular-bed pressure and position of epicardial stenosis in relation to the location of guidewire on diagnostic parameters. The pull-back image obtained by IVUS catheter and pressure pull-back profile obtained by pressure guidewire can be recorded in a catheterization laboratory and may be combined to determine distal pressure recovery and η.

Physiological limitations to the distal coronary pressureflow data
The curves can be used to evaluate non-invasive diagnostic parameters [15]. The current in-vitro experimental data only simulate epicardial coronary hemodynamic and does not simulate the coronary microvasculature. The hyperemic flow range observed in human, governed by microvasculature, can be obtained by constructing coronary flow reserve vs. hyperemic distal recovered pressure ( ) characteristic (where subscript 'h' indicates hyperemic flow). Therefore, the line was constructed based on physiological limits, without microvascular impairment, reported in Table 4. Figure 13 shows that the intersection of line and curves (for , = CFR) for each stenosis, before and during guidewire insertion, are related by: , where is hyperemic flows and . However, in clinical scenario, hemodynamic measurements are performed only by inserting guidewire or by invasive method. Consequently, hyperemic trans-stenotic pressure drop was estimated based on Fig. 13 for a given hyperemic flow for each stenosis severity. Further, estimated hyperemic pressure recovery, CDP e and η are provided in Table 5.

Eccentricity effect
Maximum difference between the experimentally measured and numerically calculated was 20% during guidewire insertion. This difference is expected due to the assumption of concentric guidewire position in the numerical formulation. (The guidewire insertion was simulated by placing a guidewire of diameter 0.35 mm in concentric configuration with the stenosis axis as shown in Fig. 6A). In contrast, the guidewire was observed to lie eccentrically and fluctuates with the pulsatile flow inside the stenotic test section during the experiments. Thus, position of guidewire inside the stenosis test section was unpredictable. The effect of guidewire eccentricity was well explained in the literatures for straight artery under pulsatile flow [29] and under steady flow [30]. To validate the eccentricity effect in moderate stenosis, 3 dimensional half plane-symmetric models were constructed ( Figure  14). The extreme conditions of guidewire placements namely, concentric and maximum eccentric positions were simulated and steady state pressure drop calculations were performed using the procedure explained in the method section. Figure 15 compares the steady state experimental pressure data with that obtained numerically for these extreme conditions. The pressure drops between above stated extreme conditions differ by ~2 mmHg which matches closely with previously reported numerical data [30]. The eccentricity effect under pulsatile flow condition should be quantified in future.
CDPe can be used to delineate the severity of coronary stenosis. Unlike FFR, CDPe is not limited between small range; therefore a better and much accurate cut-off value can be established after human clinical trials. The pressure recovery factor (η) quantifies the recovered distal vascular-bed pressure downstream of the stenosis. The geometrical details of coronary stenosis included in the expression for η combines the functional and anatomical information. CDPe and η need to be further evaluated for pre-clinical studies in order to distinguish between epicardial and microvascular dysfunction.

List of Abbreviations used
A -flow cross-sectional area, d -diameter of the vessel at specified section, l -length of vessel at specified section, p line for mean pressure and flow data Figure 13 line for mean pressure and flow data.
CFR p rh − % CFR p rh − %  Three dimensional half plane-symmetric model showing eccentricity in moderate stenosis Figure 14 Three dimensional half plane-symmetric model showing eccentricity in moderate stenosis.
and numerical data analysis and critically revised the draft. All authors read and approved the final manuscript.