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

- Rupak K Banerjee
^{1, 4}Email author, - Koustubh D Ashtekar
^{1}, - Tarek A Helmy
^{2}, - Mohamed A Effat
^{2}, - Lloyd H Back†
^{3}and - Saeb F Khoury
^{2}

**7**:24

https://doi.org/10.1186/1475-925X-7-24

© Banerjee et al; licensee BioMed Central Ltd. 2008

**Received: **26 September 2007

**Accepted: **27 August 2008

**Published: **27 August 2008

## Abstract

### 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 (CDP_{e}) 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 CDP_{e} 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 CDP_{e} is not affected that much by the presence of guidewire. CDP_{e} 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.

## Keywords

## 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 non-dimensional 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 pressure-flow scenarios.

## Methods

### Formulation

#### Hemodynamic diagnostic parameter: pressure drop coefficient (CDP_{e})

_{e}to be the ratio of mean trans-stenotic pressure drop, $\Delta \tilde{p}$, (superscript '~' indicates temporal average quantity) to the proximal dynamic pressure, i.e.,

where, *ρ* is blood density; ${\tilde{\overline{u}}}_{e}$ is spatial and temporal mean blood velocity in the proximal vessel (superscript '-' indicates a spatially average quantity). The $\Delta \tilde{p}$, in equation 1, is composed of viscous- ($\Delta {\tilde{p}}_{vis}=B\times \tilde{Q}$) and momentum-change ($\Delta {\tilde{p}}_{mom}=A\times {\tilde{Q}}^{2}$) related pressure losses, where A and B are momentum and viscous pressure loss coefficients, respectively. These pressure loss coefficients are the functions of flow rate (*Q*) and anatomical (geometrical) details of a stenosis [18]. Thus, the total pressure drop exhibits a quadratic relation with flow rate ($\Delta \tilde{p}=A\times {\tilde{Q}}^{2}+B\times \tilde{Q}$), which, in turn, influences CDP_{e} for different size of stenotic geometries.

#### Pressure recovery factor (*η*)

*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.,

*CPR*

_{∞}is expressed as

We define Δ*p*
_{
r
}= *p*
_{
r
}- *p*
_{
m
}and *α* is the ratio of the flow cross-sectional area of the distal vessel, A_{r}, to the minimum vessel area at the site of stenosis, A_{m}, prior to guidewire insertion.
*α* = *u*
_{
m
}/*u*
_{
r
}= *A*
_{
r
}/*A*
_{
m
}= *A*
_{
e
}/*A*
_{
m
}

(4)

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. Similarly, 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.
*α*
_{
i
}= (*A*
_{
r
}- *A*
_{
i
})/(*A*
_{
m
}- *A*
_{
i
}) = (*A*
_{
e
}- *A*
_{
i
})/(*A*
_{
m
}- *A*
_{
i
})

(5)

where, subscript '*i*' represents the guidewire. Further, the % occlusion of the vessel or % area stenosis is defined as follows:
Percentage area stenosis = (*A*
_{
e
}- *A*
_{
m
})/(*A*
_{
e
})

(6)

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 $\Delta {\tilde{p}}_{mom}$ and $\Delta {\tilde{p}}_{vis}$ 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

Dimensions for three stenotic test sections. All dimensions are measured with ImageJ software^{®}

Dimensions (mm) | 65% area (moderate) stenosis | 80% area (intermediate) stenosis | 89% area (severe) stenosis |
---|---|---|---|

proximal vessel diameter ( | 2.95 | 2.95 | 2.96 |

length of converging section ( | 6.96 | 6.35 | 6.28 |

throat minimal diameter ( | 1.75 | 1.32 | 0.98 |

throat length ( | 3.15 | 0.95 | 0.39 |

length of diverging section ( | 1.79 | 1.62 | 1.59 |

distal vessel diameter ( | 2.95 | 2.98 | 3.00 |

#### Experimental Setup

*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, $\tilde{\text{Q}}$ = 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].

For reference steady flow experiments, gravity induced constant flow was supplied through the stenotic test section. The steady flow (*Q*) was increased from basal to hyperemic flow in four steps, while trans-stenotic pressure and flow were measured simultaneously. During each step of flow increment, the guidewire was inserted across the stenotic test section to measure trans-stenotic pressure and flow.

#### Pressure and Flow Measurements

- 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)

*μ*

_{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.

### Numerical method

#### Geometry

*in-vitro*stenotic test sections as provided in 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

*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, $\overline{u}\left(t\right)$, needed for the computational analysis, was obtained from the mass balance equation: $\overline{u}\left(t\right)$ =

*Q*(

*t*)/

*A*

_{ e }. The $\overline{u}\left(t\right)$ 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]:

during guidewire insertion for annular flow

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, $\Delta \tilde{p}={\tilde{p}}_{r}-{\tilde{p}}_{e}$ and mean pressure recovery, $\Delta {\tilde{p}}_{r}={\tilde{p}}_{r}-{\tilde{p}}_{m}$, where as ${\tilde{p}}_{e}$, ${\tilde{p}}_{r}$ and ${\tilde{p}}_{m}$ 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 variables (*u*). For these analyses, a Compaq Linux machine with dual processor (2.4 GHz, 1 GB RAM, 80 GB hard disk) was utilized so that CPU time for each time step was ~1.5 s.

#### 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–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–27].

## Results

### Trans-stenotic mean axial pressure drop

The experimentally measured and numerically calculated temporal mean axial pressure data were analyzed to determine the overall mean trans-stenotic pressure drop ($\Delta \tilde{p}$), mean distal pressure recovery ($\Delta {\tilde{p}}_{r}$). The $\Delta {\tilde{p}}_{vis}$ and $\Delta {\tilde{p}}_{mom}$ were quantified based on the experimentally measured mean pressure drop and mean flow rate values having quadratic relation.

#### Moderate (65% area) stenosis

_{e}) increased from 125 to 414. After inserting the guidewire, mean flow rates were reduced to 54, 132, and 184 ml/min. Corresponding $\Delta \tilde{p}$ and $\Delta \tilde{p}$ are summarized in 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 signify the dominance of viscous pressure loss (*$\Delta {\tilde{p}}_{vis}$

*) over momentum change-related pressure loss (*$\Delta {\tilde{p}}_{mom}$

*) for moderate stenosis*.

Mean Pressure drop and mean pressure recovery values before and during guidewire insertion

Before guidewire insertion | |||||||
---|---|---|---|---|---|---|---|

% area stenosis | Flow ( | pressure drop ($\Delta \tilde{p}$) | pressure recovery ($\Delta {\tilde{p}}_{r}$) | Pressure recovery coefficient (CPR) | Pressure recovery factor ( | CDP | Re |

65% area (moderate) stenosis | 57 | 1.32 | 0.26 | 0.42 | 0.48 | 17 | 125 |

136 | 5.77 | 1.67 | 0.48 | 0.55 | 13 | 298 | |

189 | 9.80 | 3.50 | 0.52 | 0.59 | 11 | 414 | |

80% area (intermediate) stenosis | 54 | 6.29 | 0.80 | 0.45 | 0.47 | 89 | 122 |

82 | 12.46 | 2.10 | 0.52 | 0.54 | 77 | 195 | |

132 | 27.97 | 4.66 | 0.46 | 0.48 | 68 | 294 | |

89% area (severe) stenosis | 44 | 10.17 | 1.80 | 0.46 | 0.47 | 215 | 99 |

74 | 24.79 | 4.04 | 0.38 | 0.38 | 190 | 164 | |

93 | 37.85 | 6.05 | 0.36 | 0.36 | 184 | 207 | |

During guidewire insertion | |||||||

% area stenosis | Flow ( | pressure drop ($\Delta \tilde{p}$) | pressure recovery ($\Delta {\tilde{p}}_{r}$) | Pressure recovery coefficient (CPR) | Pressure recovery factor ( | CDP | Re |

65% area (moderate) stenosis | 54 | 1.62 | 0.05 | 0.09 | 0.11 | 22 | 109 |

132 | 7.29 | 1.25 | 0.35 | 0.41 | 17 | 262 | |

184 | 12.41 | 2.76 | 0.39 | 0.45 | 15 | 367 | |

80% area (intermediate) stenosis | 51 | 7.93 | 0.62 | 0.35 | 0.36 | 125 | 101 |

75 | 15.07 | 1.35 | 0.35 | 0.36 | 109 | 151 | |

121 | 32.18 | 4.07 | 0.41 | 0.42 | 90 | 244 | |

89% area (severe) stenosis | 34 | 10.90 | 1.41 | 0.46 | 0.46 | 372 | 68 |

61 | 26.76 | 3.57 | 0.37 | 0.37 | 291 | 122 | |

81 | 39.44 | 5.06 | 0.30 | 0.31 | 247 | 161 |

#### Intermediate (80% area) stenosis

_{e}(122–294). Because of increased shear layer instabilities with flow separation and reattachment, the distal pressure was recovered with vortical cells formations [25, 26]. The location of distal pressure recovery was approximately 14 mm downstream of the throat region for all flows. Generally, guidewire insertion adds viscous effects by increasing surface area and momentum-change effect by blocking more of the throat area. However, for intermediate stenosis, guidewire insertion relatively increased $\Delta {\tilde{p}}_{mom}$ by constricting more throat area than with $\Delta {\tilde{p}}_{vis}$. Table 2 summarizes the $\Delta \tilde{p}$ and $\Delta {\tilde{p}}_{r}$ and Fig. 8B shows the axial pressure drop profiles for 'during guidewire insertion condition'. The CPR was changed from 0.35 to 0.41 for Re

_{e}range of 101 to 150. The distal pressure recovery occurred approximately 13 mm distal to the diverging section. For larger mean flow rates, $\Delta {\tilde{p}}_{mom}$ was more dominant than $\Delta {\tilde{p}}_{vis}$, due to reduction in the throat area, as compared to moderate stenosis. The clinically measured fractional flow reserve (FFR) falls below the limiting value of 0.75 (leading to angioplasty or stent placement) [6] because of guidewire insertion; hence, true severity of intermediate stenosis is very difficult to diagnose.

#### Severe (89% area) stenosis

_{e}range. However, CPR was reduced from 0.46 to 0.36 for Re

_{e}range of 99–207. Before guidewire insertion, the vortical flow cells are much stronger than that during guidewire insertion. 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 $\Delta {\tilde{p}}_{mom}$. 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

Viscous (A) and momentum-change (B) coefficients; Pulsatile flow values are obtained from Figs 7A, B, and 7C

Flow status | Steady flow | Pulsatile flow | ||||||
---|---|---|---|---|---|---|---|---|

Stenosis severity | Before guidewire | During Guidewire | Before guidewire | During Guidewire | ||||

A | B | A | B | A | B | A | B | |

Moderate | 1.46E-04 | 1.00E-02 | 1.87E-04 | 1.85E-02 | 2.01E-04 | 1.44E-02 | 2.34E-04 | 2.39E-02 |

Intermediate | 4.70E-04 | 2.65E-02 | 8.00E-04 | 4.34E-02 | 1.35E-03 | 3.71E-02 | 1.55E-03 | 7.60E-02 |

Severe | 2.50E-03 | 2.79E-02 | 4.00E-03 | 6.01E-02 | 3.20E-03 | 9.64E-02 | 5.00E-03 | 1.13E-01 |

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, $\Delta {\tilde{p}}_{vis}$, which was further augmented by guidewire insertion (Fig. 10A; coefficient 'B' in Table 3). Due to dominance of $\Delta {\tilde{p}}_{mom}$, $\Delta \tilde{p}$ 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, from moderate and intermediate stenosis (80% area stenosis), respectively. During guidewire insertion, coefficient 'A' for severe stenosis increased by factors of 21 (= 50 × 10^{-4}/2.34 × 10^{-4}) and 3.2 (= 50 × 10^{-4}/15.5 × 10^{-4}), as shown in Table 3, from moderate and intermediate stenosis, respectively. *This confirms the dominance of* $\Delta {\tilde{p}}_{mom}$
*for severe stenosis*. Interestingly, the intermediate stenosis was a trade-off point for which the dominance of $\Delta {\tilde{p}}_{vis}$ and $\Delta {\tilde{p}}_{mom}$ 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 $\Delta {\tilde{p}}_{mom}$ considerably than $\Delta {\tilde{p}}_{vis}$. *The pulsatile flow enhances momentum-change pressure loss due to the significant effect of convective term:* $\partial \overline{u}/\partial t$.

### Overall pressure drop coefficient (CDP_{e})

_{e}, a non-dimensional quantity, represents the fluid flow resistance added by any arterial conduit having both momentum-change- and viscous- related pressure losses. The coronary stenosis adds resistance to flow by blocking more lumen area and its severity ranges from moderate stenosis offering negligible blood flow resistance to the severe stenosis contributing considerable 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.

Before guidewire insertion, from basal to hyperemic flow, CDP_{e} decreased from 215 to 184 for severe stenosis and from 17 to 12 for moderate stenosis. Similarly, during guidewire insertion, CDP_{e} changed from 372 to 247 for severe stenosis and from 22 to15 for moderate stenosis. The statistically significant difference in the range of CDP_{e} values between moderate and severe stenosis can be used in the clinical practice in diagnosing the severity of stenosis.

_{e}correlation showing a linear variation. The CDP

_{e}increased by a factor of 15 from moderate to severe stenosis model. This wide range of CDP

_{e}can be used to delineate the severity of coronary stenosis. Unlike FFR, CDP

_{e}is not limited between small range; therefore a better and much accurate cut-off value can be established for CDP

_{e}after human clinical trials. The linear correlation indicates CDP

_{e}could be a viable diagnostic parameter under clinical setting.

### Pressure recovery factor (*η*)

Pressure recovery, downstream of the stenosis, can be quantified by the pressure recovery factor (*η*). For given flow rate, as stenotic severity increases from moderate to intermediate, the CPR increases, signifying increased pressure recovery, while approaching CPR_{8}. As a result, *η* could also increase to a limiting value of unity. Similarly, for a particular stenosis, as flow increases, *η* increases due to increase in CPR. Table 2 compares the *η* values, calculated before and during guidewire insertion. From basal to hyperemic flow, during guidewire insertion, *η* increased from 0.11 to 0.45 for moderate stenosis and from 0.36 to 0.42 for intermediate stenosis. However, it decreased from 0.46 to 0.31 for severe stenosis. Similarly, before guidewire insertion, *η* increased from 0.48 to 0.59 for moderate stenosis, remained nearly constant (0.47 to 0.48) for intermediate stenosis and reduced from 0.47 to 0.36 for severe stenosis. This indicates that the guidewire insertion reduces the ability of the diverging section to recover the downstream pressure. In future, the CPR and *η* values can be obtained by guidewire pressure pullback 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 $\Delta {\tilde{p}}_{vis}$ and $\Delta {\tilde{p}}_{mom}$ 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, having different $\partial \overline{u}/\partial t$ 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} (a limiting value of pressure drop coefficient at very high flow where viscous losses can be neglected over momentum-change losses) and square root of CDP_{m} (pressure drop coefficient evaluated at the site of constriction at hyperemic flow where viscous and momentum-change pressure losses are equally important). As the flow increases, the CDP_{m} approaches CDP_{8}, leading to an LFC value of unity. The relation between CDP_{m}, CDP_{e}, CDP_{8} and LFC are summarized in Additional file 2. *CDP*
_{∞} can be expressed as (1-1/*α*)^{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 change-related 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 *in-vitro* 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 pressure-flow data

*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 ($CFR-{\tilde{p}}_{rh}$) characteristic (where subscript '

*h*' indicates hyperemic flow). Therefore, the $CFR-{\tilde{p}}_{rh}$ line was constructed based on physiological limits, without microvascular impairment, reported in Table 4. Figure 13 shows that the intersection of $CFR-{\tilde{p}}_{rh}$ line and $\Delta {\tilde{p}}_{h}$ curves (for ${\tilde{p}}_{rh}$, $\tilde{Q}/{\tilde{Q}}_{b}-{\tilde{p}}_{rh}$ = CFR) for each stenosis, before and during guidewire insertion, are related by: $\tilde{Q}={\tilde{Q}}_{h}$, where ${\tilde{Q}}_{h}/{\tilde{Q}}_{b}$ is hyperemic flows and ${\tilde{p}}_{rh}={\tilde{p}}_{a}-A\times {\tilde{Q}}_{h}^{2}-B\times {\tilde{Q}}_{h}$. 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.

Mean Hemodynamic end points in terms of CFR and ${\tilde{Q}}_{h}$

Stenosis | Coronary flow reserve (CFR) | Proximal pressure (${\tilde{p}}_{a}$) mmHg | Hyperemic trans-stenotic pressure drop ($\Delta \tilde{p}={\tilde{p}}_{rh}-{\tilde{p}}_{a}$) mmHg | Hyperemic distal recovered pressure (${\tilde{Q}}_{h}$) mmHg | Fractional flow reserve (FFR) |
---|---|---|---|---|---|

Moderate [13] | 3.6 ± 0.3 | 84 ± 3 | 7.4 | 75.2 | 0.89 |

Intermediate [16] | 3.3 | 86 | 14.3 | 70.4 | 0.82 |

2.3 ± 0.1 | 89 ± 3 | 34 | 55 | 0.62 |

Estimated hyperemic diagnostic parameters during guidewire insertion

% area stenosis | Estimated Hyperemic flow (from Fig. 10) | Estimated pressure drop (mmHg) $\Delta \tilde{p}=A\times {\tilde{Q}}^{2}+B\times \tilde{Q}$ | Pressure recovery $\Delta {\tilde{p}}_{r}$ (mmHg) | Pressure recovery coefficient (CPR) | Pressure recovery factor ( | CDP |
---|---|---|---|---|---|---|

65% area (moderate) stenosis | 172 | 11 | 2.5 | 0.41 | 0.46 | 15 |

80% area (Intermediate) stenosis | 116 | 31 | 3.7 | 0.40 | 0.42 | 96 |

89% area (Severe) stenosis | 80 | 45 | 5.0 | 0.31 | 0.31 | 292 |

### Eccentricity effect

## Conclusion

The guidewire insertion caused reduction in hyperemic flow and increased the trans-stenotic pressure drop, that affected the viscous- and momentum change-related pressure losses. The pressure drop coefficient (CDPe) increased by a factor of 15 from moderate to severe stenosis model and thus diminishes the overall effect of guidewire flow obstruction and related misdiagnosis during invasive functional measurements. This wide range of 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.

## Notes

## List of Abbreviations used

- Subscripts
*a*– aorta: -
*c*– converging section of the stenosis,*e*– proximal section to the stenosis,*h*– hyperemic condition for flow and pressure,*i*– guidewire,*m*–*t*– minimum temporal value,*m*– throat section,*o*– mean diameter,*p*–*t*– peak temporal value,*r*– distal section/diverging region of the stenosis.

- Superscripts
^{~}– time average (mean) quantity over cardiac cycle: -
^{-}– spatial average quantity across the vessel

## Declarations

### Acknowledgements

This work is supported by American Heart Association: National-Scientific Development Grant (#0335270N). Authors would like to thank Dr. Bernard J. Dardzinski of Cincinnati Children's hospital, for providing the MicroCT facility used for the accurate dimensional measurements of the test sections. We are also grateful to Tristate cutting tool Inc, for manufacturing the stenosis test models to highest precision level. Authors want to acknowledge the help of Dr. Abhijit Sinha Roy during the experimental setup and test section design. We are thankful to all members of the Bio-Fluid, Heat, and Mass Transport Lab (Bio-FHM), Mr. Doug Hurd of Mechanical engineering machine shop at University of Cincinnati for their assistance during these *in-vitro* experiments.

## Authors’ Affiliations

## References

- Kern MJ, Lerman A, Bech JW, De Bruyne B, Eeckhout E, Fearon WF, Higano ST, Lim MJ, Meuwissen M, Piek JJ, Pijls NH, Siebes M, Spaan JA:
**Physiological assessment of coronary artery disease in the cardiac catheterization laboratory: a scientific statement from the American Heart Association Committee on Diagnostic and Interventional Cardiac Catheterization.***Circulation*2006,**114**(12):1321–1341. 10.1161/CIRCULATIONAHA.106.177276View ArticleGoogle Scholar - Back LH, Kwack EY, Back MR:
**Flow rate-pressure drop relation in coronary angioplasty: catheter obstruction effect.***J Biomech Eng*1996,**118**(1):83–89. 10.1115/1.2795949View ArticleGoogle Scholar - Pijls NH, De Bruyne B:
*Coronary pressure.*2nd edition. Dordrecht, The Netherlands , Kluwer Academic; 2000.View ArticleGoogle Scholar - Gould KL, Lipscomb K, Hamilton GW:
**Physiologic basis for assessing critical coronary stenosis. Instantaneous flow response and regional distribution during coronary hyperemia as measures of coronary flow reserve.***Am J Cardiol*1974/01/01 edition. 1974,**33**(1):87–94. 10.1016/0002-9149(74)90743-7View ArticleGoogle Scholar - Anderson HV, Roubin GS, Leimgruber PP, Cox WR, Douglas JS Jr., King SB 3rd, Gruentzig AR:
**Measurement of transstenotic pressure gradient during percutaneous transluminal coronary angioplasty.***Circulation*1986,**73**(6):1223–1230.View ArticleGoogle Scholar - Pijls NH, De Bruyne B, Peels K, Van Der Voort PH, Bonnier HJ, Bartunek JKJJ, Koolen JJ:
**Measurement of fractional flow reserve to assess the functional severity of coronary-artery stenoses.***N Engl J Med*1996/06/27 edition. 1996,**334**(26):1703–1708. 10.1056/NEJM199606273342604View ArticleGoogle Scholar - de Bruyne B, Bartunek J, Sys SU, Pijls NH, Heyndrickx GR, Wijns W:
**Simultaneous coronary pressure and flow velocity measurements in humans.***Circulation*1996,**94**(8):1842–1849.View ArticleGoogle Scholar - Fearon WF, Nakamura M, Lee DP, Rezaee M, Vagelos RH, Hunt SA, Fitzgerald PJ, Yock PG, Yeung AC:
**Simultaneous assessment of fractional and coronary flow reserves in cardiac transplant recipients: Physiologic Investigation for Transplant Arteriopathy .***Circulation*2003,**108**(13):1605–1610. 10.1161/01.CIR.0000091116.84926.6FView ArticleGoogle Scholar - Siebes M, Verhoeff BJ, Meuwissen M, de Winter RJ, Spaan JA, Piek JJ:
**Single-wire pressure and flow velocity measurement to quantify coronary stenosis hemodynamics and effects of percutaneous interventions.***Circulation*2004,**109**(6):756–762. 10.1161/01.CIR.0000112571.06979.B2View ArticleGoogle Scholar - Jenni R, Buchi M, Zweifel HJ, Ritter M:
**Impact of Doppler guidewire size and flow rates on intravascular velocity profiles.***Cathet Cardiovasc Diagn*1998,**45**(1):96–100. 10.1002/(SICI)1097-0304(199809)45:1<96::AID-CCD23>3.0.CO;2-IView ArticleGoogle Scholar - Verberne HJ, Meuwissen M, Chamuleau SA, Verhoeff BJ, van Eck-Smit BL, Spaan JA, Piek JJ, Siebes M:
**Effect of simultaneous intracoronary guidewires on the predictive accuracy of functional parameters of coronary lesion severity.***Am J Physiol Heart Circ Physiol*2007,**292**(5):H2349-H2355. 10.1152/ajpheart.01042.2006View ArticleGoogle Scholar - De Bruyne B, Pijls NH, Paulus WJ, Vantrimpont PJ, Sys SU, Heyndrickx GR:
**Transstenotic coronary pressure gradient measurement in humans: in vitro and in vivo evaluation of a new pressure monitoring angioplasty guide wire.***J Am Coll Cardiol*1993,**22**(1):119–126.View ArticleGoogle Scholar - Wilson RF, Johnson MR, Marcus ML, Aylward PE, Skorton DJ, Collins S, White CW:
**The effect of coronary angioplasty on coronary flow reserve.***Circulation*1988,**77**(4):873–885.View ArticleGoogle Scholar - Emanuelsson H, Dohnal M, Lamm C, Tenerz L:
**Initial experiences with a miniaturized pressure transducer during coronary angioplasty.***Cathet Cardiovasc Diagn*1991,**24**(2):137–143. 10.1002/ccd.1810240213View ArticleGoogle Scholar - Ashtekar KD, Back LH, Khoury SF, Banerjee RK:
**In-vitro quantification of guidewire flow obstruction effect in model coronary stenoses for interventional diagnostic procedure.***ASME J Medical Devices*2007,**1**(3):185–196. 10.1115/1.2776336View ArticleGoogle Scholar - Roy AS, Banerjee RK, Back LH, Back MR, Khoury S, Millard RW:
**Delineating the guide-wire flow obstruction effect in assessment of fractional flow reserve and coronary flow reserve measurements.***Am J Physiol Heart Circ Physiol*2005,**289**(1):H392-H397. 10.1152/ajpheart.00798.2004View ArticleGoogle Scholar - Banerjee RK, Sinha Roy A, Back LH, Back MR, Khoury SF, Millard RW:
**Characterizing momentum change and viscous loss of a hemodynamic endpoint in assessment of coronary lesions.***J Biomech*2007,**40**(3):652–662. 10.1016/j.jbiomech.2006.01.014View ArticleGoogle Scholar - Young DF, Tsai FY:
**Flow characteristics in models of arterial stenoses. II. Unsteady flow.***J Biomech*1973/09/01 edition. 1973,**6**(5):547–559. 10.1016/0021-9290(73)90012-2View ArticleGoogle Scholar - Lefebvre AH:
*Gas turbine combustion.*2nd edition. Philadelphia, PA , Taylor & Francis; 1999.Google Scholar - Back LH, Denton TA:
**Some Arterial Wall Shear Stress Estimates in Coronary Angioplasty.***Angiology*1997,**48**(2):99–109. 10.1177/000331979704800201View ArticleGoogle Scholar - Cho YI, Back LH, Crawford DW, Cuffel RF:
**Experimental study of pulsatile and steady flow through a smooth tube and an atherosclerotic coronary artery casting of man.***J Biomech*1983,**16**(11):933–946. 10.1016/0021-9290(83)90057-XView ArticleGoogle Scholar - Segal J, Kern MJ, Scott NA, King SB 3rd, Doucette JW, Heuser RR, Ofili E, Siegel R:
**Alterations of phasic coronary artery flow velocity in humans during percutaneous coronary angioplasty.***J Am Coll Cardiol*1992,**20**(2):276–286.View ArticleGoogle Scholar - Brookshier KA, Tarbell JM:
**Evaluation of a transparent blood analog fluid: aqueous xanthan gum/glycerin.***Biorheology*1993/03/01 edition. 1993,**30**(2):107–116.Google Scholar - Cho YI, Kensey KR:
**Effects of the non-Newtonian viscosity of blood on flows in a diseased arterial vessel. Part 1: Steady flows.***Biorheology*1991/01/01 edition. 1991,**28**(3–4):241–262.Google Scholar - Banerjee RK, Back LH, Back MR, Cho YI:
**Physiological flow simulation in residual human stenoses after coronary angioplasty.***J Biomech Eng*2000,**122**(4):310–320. 10.1115/1.1287157View ArticleGoogle Scholar - Banerjee RK, Back LH, Back MR, Cho YI:
**Physiological flow analysis in significant human coronary artery stenoses.***Biorheology*2003,**40**(4):451–476.Google Scholar - FIDAP manual 10 Cavendish Court, Lebanon, NH 03766, USA , Fluent Inc; 1998.Google Scholar
- Brown BG, Bolson EL, Dodge HT:
**Dynamic mechanisms in human coronary stenosis.***Circulation*1984,**70**(6):917–922.View ArticleGoogle Scholar - Dapira P, BDash RK:
**A Numerical Study of Pulsatile Blood Flow in an Eccentric Catheterized Artery Using a Fast Algorithm.***J Engg Mathematics*2002,**42:**1–22. 10.1023/A:1014332225766View ArticleGoogle Scholar - Back LH:
**Estimated mean flow resistance increase during coronary artery catheterization.***J Biomech*1994/02/01 edition. 1994,**27**(2):169–175. 10.1016/0021-9290(94)90205-4View ArticleGoogle Scholar

## Copyright

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.