- Research
- Open Access

# Numerical analysis of the pressure drop across highly-eccentric coronary stenoses: application to the calculation of the fractional flow reserve

- R. Agujetas
^{1}, - M. R. González-Fernández
^{2}, - J. M. Nogales-Asensio
^{2}and - J. M. Montanero
^{1}Email authorView ORCID ID profile

**Received:**18 December 2017**Accepted:**23 May 2018**Published:**30 May 2018

## Abstract

### Background

Fractional flow reverse (FFR) is the gold standard assessment of the hemodynamic significance of coronary stenoses. However, it requires the catheterization of the coronary artery to determine the pressure waveforms proximal and distal to the stenosis. On the contrary, computational fluid dynamics enables the calculation of the FFR value from relatively non-invasive computed tomography angiography (CTA).

### Methods

We analyze the flow across idealized highly-eccentric coronary stenoses by solving the Navier–Stokes equations. We examine the influence of several aspects (approximations) of the simulation method on the calculation of the FFR value. We study the effects on the FFR value of errors made in the segmentation of clinical images. For this purpose, we compare the FFR value for the nominal geometry with that calculated for other shapes that slightly deviate from that geometry. This analysis is conducted for a range of stenosis severities and different inlet velocity and pressure waveforms.

### Results and conclusions

The errors made in assuming a uniform velocity profile in front of the stenosis, as well as those due to the Newtonian and laminar approximations, are negligible for stenosis severities leading to FFR values around the threshold 0.8. The limited resolution of the stenosis geometry reconstruction is the major source of error when predicting the FFR value. Both systematic errors in the contour detection of just 1-pixel size in the CTA images and a low-quality representation of the stenosis surface (coarse faceted geometry) may yield wrong outcomes of the FFR assessment for an important set of eccentric stenoses. On the contrary, the spatial resolution of images acquired with optical coherence tomography may be sufficient to ensure accurate predictions for the FFR value.

## Keywords

- CFD
- Stenosis
- Fractional flow reserve

## Background

A stenosis is an abnormal narrowing in a blood vessel. Coronary arteries are known to be common sites for stenoses, which frequently lead to myocardial infarction and even sudden cardiac death. The vast majority of these stenoses are caused by atheromatous lesions, which give rise to highly eccentric narrowing. They are frequently characterized by a lumen lying in the outer region of the coronary artery, and delimited by an arc of nearly normal wall [1].

Coronary revascularization is performed to treat stenoses with hemodynamic significance. The decision to carry out this procedure is frequently based on the stenosis severity, defined as the percentage of narrowing of the coronary artery. Thus, coronary revascularization is often conducted for severe stenoses, i.e., those with severities above the threshold value 75%. While the stenosis severity can be measured with relatively non-invasive computed tomography angiography (CTA), there is no simple relationship between that parameter and the stenosis hemodynamic significance and the plaque vulnerability [2]. In fact, more than 50% of the stenoses judged severe by CTA do not cause ischemia [3]. Therefore, criteria based on measured stenosis hemodynamic significance seem to be more appropriate.

Fractional flow reverse (FFR) is defined as the ratio of maximal (i.e., on maximum hyperaemia conditions) coronary blood flow through a stenotic artery to the maximal blood flow in the hypothetical case that the artery were normal. Because this definition is not operational, it is frequently substituted by \(\text{FFR}\equiv \overline{P}_1/\overline{P}_0,\) where \(\overline{P}_0\) and \(\overline{P}_1\) are the averages over the cardiac cycle of the pressure waveforms \(P_0(t)\) and \(P_1(t)\) proximal and distal to the stenosis, respectively, both measured at maximum hyperaemia. Assuming that the loss of pressure in a normal epicardial artery is negligible as compared to that in the coronary microcirculation, the two definitions would be equivalent if the distal pressure were directly proportional to the coronary blood flow [4]. This last condition does not strictly hold due to the colateral flow, among other factors, as indicated by the non-zero pressure intercept for flow in the coronary circulation [5]. In any case, the above fraction is adopted as the operational definition of FFR.

FFR is the gold standard assessment of the hemodynamic significance of coronary stenoses. Revascularization is typically conducted for lesions with FFR < 0.8 [5]. FFR-guided revascularization has proved to be long-lived and cost-saving [6]. As compared with other assessments, such as angiography for multivessel evaluation (AME), FFR-guided revascularization exhibits lower rates of major adverse cardiac events. In addition, it correctly classifies as hemodynamically insignificant many lesions judged “obstructive” by AME [7]. The main drawback of the FFR-based assessment is its invasive character. In fact, this method uses pressure wires to determine the pressure waveforms proximal and distal to the stenosis, which requires the catheterization of the coronary arteries. Less invasive methods to determine the FFR value of coronary stenoses are very desirable.

Computational fluid dynamics (CFD) is a method with high diagnostic performance to determine coronary lesions that cause ischemia [8–13]. In particular, it enables the calculation of the FFR value from relatively non-invasive measurements [10, 11, 14]. CFD has been frequently used to analyze idealized 2D (axisymmetric) stenoses and understand relevant aspects of hemodynamics in these geometries [15–20], especially turbulence. There are much fewer studies dealing with 3D (asymmetric) shapes, most of them restricted to small eccentricities [21–23]. As mentioned above, many of the coronary stenoses exhibit highly eccentric shapes whose necks are delimited by an arc of the normal artery wall and a atheromatous lesion. Analyzing the influence of different hemodynamics aspects and model approximations on the pressure loss across these stenoses is of great importance, and has practical consequences in calculating the FFR value.

As mentioned above, FFR is measured at maximum hyperaemia, which greatly reduces the artery expansion/contraction during the systole/diastole phase. Therefore, one can safely adopt the rigid wall approximation in the CFD simulations, which avoids modelling the fluid solid interaction [24, 25]. FFR is a measure of the mechanical energy dissipated by viscosity in the stenosis averaged over the cardiac cycle. In some phases of this cycle, the blood speed decreases and so does the shear rate over the entire fluid domain. Therefore, it is natural to wonder whether using the Newtonian approximation with the viscosity at infinite shear rate is a good approximation to get reliable predictions for the FFR value. Here, we will answer this question by calculating the pressure drop across a stenosis with both the Newtonian and Carreau models.

It is well known that, due to the stabilizing effect of the blood acceleration, turbulence in arteries arises at Reynolds numbers much larger than those for steady flow. The peak Reynolds number in a coronary stenosis is typically much smaller than that leading to the laminar-to-turbulent transition in an oscillatory flow inside a cylindrical duct, and consequently the flow is usually assumed to be laminar. However, blood greatly accelerates in severe stenoses (like the ones studied here), which gives rise to complex post-stenotic flows, including separation, recirculation, as well as localized transition to turbulence and relaminarization [15]. Both the inlet flow perturbations and the stenosis asymmetry enhance the post-stenotic transition to turbulence [21], and thus considerable turbulent outbreaks can appear in that region even for inlet Reynolds numbers as low as 500 [17, 26]. CFD simulations of the coronary system are conducted in the laminar regime because most of the flow is truly laminar, and turbulence (if there is any) is restricted to the post-stenotic region. However, the flow in this very region essentially determines the FFR value, and therefore significant errors may be made in neglecting turbulence. Here, we will evaluate the influence of turbulence on the pressure drop across the stenosis by using the large eddy simulation (LES) model, which has shown to provide good results for asymmetric shapes [21].

Typically, the only (if there is any) available kinematic information at the stenosis inlet section is the instantaneous flow rate crossing that section [27, 28]. To complete the velocity inlet boundary condition, one frequently assumes the uniform velocity profile compatible with the measured flow rate. In fact, the pulsatility of the flow in arteries hinders the inwards viscous diffusion of momentum from the wall, which explains why the uniform velocity approximation is commonly preferred. However, significant extra loss of pressure may occur in coronary stenoses due to the growth of the boundary layer when that approximation is made [29]. In this work, we will examine the influence of this choice on the pressure drop across the stenosis.

Simulating the entire cardiac cycle increases very considerably the computing time. It implies running the simulation over several cycles to ensure that the periodic regime has been reached. The Womersley number (the ratio between local acceleration and viscous force per unit mass) takes relatively low values in the epicardial arteries. One may wonder whether the flow pulsatile character can be neglected in calculating the FFR value. In that case, the unsteady simulation could be replaced with a stationary one with an average (effective) inlet velocity. Here, we will assess to what extent the instantaneous pressure drop across the stenosis is affected by the blood acceleration.

Ideally, the FFR quantification should be conducted by simulating the flow in the stenotic region exclusively. This would allow one to increase the spatial resolution of the simulation, capturing the critical influence of the intricate stenosis shape on the FFR value. However, such a simulation requires the precise knowledge of the boundary conditions in the pre- and post-stenotic segments, which cannot be easily determined. For this reason, the numerical simulation normally involves not only the coronary artery where the stenosis is located, but also the peripheral vasculature, including the aorta. The inlet and outlet boundary conditions for this numerical domain are relatively well established, which allows CFD to produce FFR predictions in an autonomous way, i.e., just from the knowledge of the vasculature geometry and some patient-specific parameters.

The above mentioned approach does not lack important disadvantages. Even if the boundary conditions are correctly prescribed [30–32], the disparity between the sizes of the coronary stenosis and the whole numerical domain constitutes an important barrier to obtain accurate results. It is obvious that the pressure drop across the stenosis critically depends on the complicated shape of the latter. The vasculature geometrical reconstruction has limited accuracy. The stenosis shape is determined from segmentation of CTA images of artery sections. This segmentation is carried out by manually selecting a range of values in the Hounsfield scale corresponding to the area enclosed by the artery inner wall. Pixels that have been excluded (included) erroneously in this process are added (eliminated) also manually. The contours are the boundaries that separate the selected pixels from the rest. Finally, the detected contours are used to construct a 3D representation of the real vasculature geometry in terms of triangles connected to each other (the faceted geometry). One can identify at least two sources of errors in this process. Firstly, the contour detection produces errors in the artery wall position at least on the order of 1-pixel size. Lastly, both hardware and software characteristics limit the total number of triangles used in the geometrical representation of the stenosis. It must be pointed out that a coarse faceted geometry hinders refining the simulation mesh in that critical region, and prevents generating the grid necessary to accurately calculate thin boundary layers. This factor may become important when those layers separate from the artery surface under the action of adverse pressure gradients.

In this paper, we will analyze the influence on the pressure drop across a highly-eccentric coronary stenosis of several aspects (approximations) of the simulation method. Specifically, we will consider the accuracy and spatial resolution achieved in reconstructing the stenosis geometry, the velocity distribution at the stenosis inlet, the post-stenotic turbulence, the blood non-Newtonian character, and the values taken by the Reynolds and Womersley numbers. The stenosis geometry will be parameterized by adopting an idealized model that represents a wide class of atheromatous lesions. First, we will simulate the flow through stenoses with varied severities. Then, we will focus on a stenosis with FFR around the critical value 0.8 to determine to what extent the factors mentioned above can influence the outcome of the FFR assessment. We will conclude that this is a geometry-dominated problem, and therefore the limited accuracy and spatial resolution of the stenosis reconstruction is the major source of error in predicting the FFR value. To the best of our knowledge, this aspect of the problem has not been quantitatively analyzed yet.

## Formulation of the problem

*L*and diameter

*D*. This subtraction is carried out at a distance \(L_s\) from the cylinder inlet. The pre-stenotic \(P_0\) and post-stenotic \(P_1\) pressures are measured at that distance form the stenosis neck. For fixed values of

*L*/

*D*, \(D_s/D\) and \(L_s/D,\) the fluid domain is univocally characterized by the stenosis severity \({\mathcal {S}}=1-4A_{\text{min}}/(\pi D^2),\) where \(A_{\text{min}}\) is the minimum cross section area Fig. 3b. This asymmetric stenosis is similar to the real cases described in the “Background” section and considered in, e.g., Refs. [33] and [34].

As mentioned in the “Background” section, the segmentation of the CTA images provides a representation of the real geometry in terms of triangles connected to each other. In order to quantify the quality of that representation, we define the stenosis Surface Resolution Index SRI as \(\text{SRI}\equiv N/(S/D^2),\) where *N* is the number of triangles in the stenosis and *S* the area of the surface represented by those triangles. The ratio \(S/D^2\) indicates the “complexity” of the surface analyzed: the larger \(S/D^2\) the more complex that surface. Therefore, \(N/(S/D^2)\) can be regarded as a measure of the surface representation resolution in terms of its complexity. We assign the value \(\text{SRI}=\infty\) to the perfectly smooth surface resulting from a mathematical model like the one sketched in Fig. 3.

*r*is the distance to the cylinder axis, and \(v_0(t)\) is the velocity averaged over the inlet section. The velocity distribution (3) is the solution for an infinite pipe subject to a harmonic pressure gradient characterized by a Womersley number value equal to \(\eta^{-1}\) [35]. This distribution is considered here just as a family of axisymmetric velocity profiles, where \(\eta\) is a dimensionless parameter that measures the deviation from the uniform distribution. In fact, the uniform and parabolic velocity profiles correspond to \(\eta=0\) and \(\infty,\) respectively (Fig. 4).

^{−1}respectively.

## Numerical method

In this work, the Navier–Stokes equations were integrated in the incompressible regime with the finite volume method [40] implemented in the commercial software Fluent [41]. Both steady and pulsatile flow simulations were conducted. When turbulence was accounted for, the LES equations were integrated with implicit filtering by calculating the subgrid scale Reynolds stress tensor with the Smagorinsky model. The Smagorinsky coefficient was set to the value 0.13, which leads to accurate predictions in stenosed arteries [22].

The velocity profile (3) was imposed at the inlet section. In the unsteady laminar simulations conducted to calculate the FFR value, we set \(\eta=0\) in Eq. (3), which corresponds to the common simplification of uniform flow at the inlet section. The flow was assumed to be fully developed at the outlet section given the large value of *L*/*D*. For this reason, outflow boundary conditions (\(\partial a/\partial n=0,\) where *a* is any hydrodynamic quantity and *n* the direction normal to the surface) were prescribed at that section. We verified that this choice leads to virtually the same steady pressure drops as those calculated for constant outlet pressure. On the other side, the non-slip boundary condition was prescribed at the artery inner surface. In the LES simulations, the turbulence intensity was set to zero at the inlet section.

The numerical integration of the hydrodynamic equations was carried out with the pressure-based solver. The gradients in the cell center and faces were evaluated with the green-Gauss node-based scheme and a multidimensional Taylor expansion, respectively. The spatial discretization of the pressure equation was conducted with the second-order approximation, while the momentum equations were discretized with the second-order upwind scheme. The velocity-pressure coupling was conducted with the SIMPLE procedure. In the unsteady laminar simulations, the time step \(\Delta t^*=0.005\) was much smaller than the period of the cardiac cycle. In the LES simulations, the time step was reduced to \(\Delta t^*=10^{-4}.\) The Courant number \(C=v\Delta t /\Delta x\) based on the cell size \(\Delta x,\) the velocity *v* in the cell, and the time step \(\Delta t^*,\) took values in the range 0.1–10. The average value was around 0.6. To check the solution convergence, the pressure drop across the stenosis was monitored in the course of the simulation. In order to reach the quasi-periodic regime in the laminar unsteady case, the simulation was run over three cardiac cycles. The simulation results presented in “Results and discussion” section correspond to the last cycle. The turbulent simulations were run until the flow became statistically steady.

## Results and discussion

### Pressure drop across the stenosis

Figure 7 shows the values of \(\Pi\) calculated for different stenosis severities. There is a strong influence of this parameter on the pressure drop for severities around 0.6. As will be seen in “Influence of the segmentation error on the fractional flow reverse” section, this value approximately corresponds to the threshold FFR = 0.8 for clinical intervention when the inlet velocity waveform in Fig. 5 is assumed. Therefore, one expects that small variations of the stenosis shape may lead to different outcomes of the FFR assessment. The results presented in Fig. 7 allows one to estimate the errors in \(\Pi\) due to the image segmentation process for our stenosis geometrical model. Consider a coronary artery of diameter *D* with a stenosis like that represented in Fig. 3. Assume that the image segmentation produces an error \(\delta D_s\) when determining \(D_s\) from CTA. That error translates into a variation of the area \(A_{\mathrm{min}},\) and therefore of the severity degree \({\mathcal {S}}.\) The resulting deviation of the pressure drop \(\Pi\) from its true value can be calculated from the curve \(\Pi ({\mathcal {S}}).\) As an example, consider the following values: \(D=5.15\) mm (estimated from Fig. 2), \(D_s=10\) mm, \({\mathcal {S}}=0.583\) and \(\delta D_s\) = + 300 μm (the positive/negative sign stands for an increase/decrease of \(D_s\)). The latter corresponds to 1-pixel size in Fig. 2b. Due to the error \(\delta D_s,\) the stenosis severity and corresponding pressure drop increase from (\(S=0.583,\) \(\Pi =4.03\)) up to \((S=0.63,\;\Pi =4.93).\) To get variations of \(\Pi\) smaller than 5%, the segmentation error should be smaller than 100 μm, around one third of 1-pixel size in Fig. 2b. This example shows the importance of accurate image segmentation in the analysis of highly-eccentric stenosis with FFR close to the critical value.

As mentioned in the “Background” section, the surface of a real stenosis is represented in a Fluent simulation in terms of triangles connected to each other, and is stored in the so-called STL format [41]. Another potential source of error in calculating the FFR value is the limited resolution of this faceted geometry. To analyze this aspect of the problem, we first considered the perfectly smooth surface resulting from the mathematical model sketched in Fig. 3. We assign the surface resolution index \(\text{SRI}=\infty\) to this surface. Then, several faceted surfaces consisting of different numbers of triangles were generated from the mathematical model. This process was automatically conducted by Fluent when exporting the smooth geometry to the STL format. Finally, we calculated the pressure drop \(\Pi\) and the corresponding SRI value for each of those surfaces. It must be noted that the results may significantly depend not only on the SRI value but also on the specific method used to build the faceted geometry. Figure 8 shows \(\Pi\) as a function of SRI for the characteristic case \({\mathcal {S}}=0.583.\) The SRI value of the point marked with an arrow coincides with that of the stenosis shown in Fig. 2c, and therefore it can be regarded as a realistic value. The maximum value of SRI corresponds to the maximum number of triangles that can be generated from the procedure described above. The horizontal line indicates the value of \(\Pi\) obtained for \(\text{SRI}=\infty.\) The pressure drop of the marked point deviates around 19% from that value. Therefore, this aspect of the simulation may constitute a considerable source of error as well.

Figure 7 shows the pressure drops for the uniform \((\eta=0)\) and parabolic (\(\eta=\infty\)) inlet velocity distributions considering different stenosis degrees \({\mathcal {S}}.\) As can be observed, the uniform and parabolic profiles lead to significantly different results for small \({\mathcal {S}}.\) In fact, the growth of the boundary layer for the plug flow case causes an extra loss of pressure of order unity, and therefore increases the value of \(\Pi\) with respect to that of the parabolic profile. This effect was predicted analytically in Ref. [29]. For higher values of \({\mathcal {S}},\) the drop of pressure takes place essentially in the stenotic region, where the flow loses memory of the inlet velocity distribution. For this reason, the values of \(\Pi\) obtained in the two cases practically coincide. It must be noted that the numerical simulations provide useful clinical information for our stenosis model when the severity is about 60%, because this case corresponds to \(\text{FFR}\simeq\) 0.8 (using the inlet velocity waveform in Fig. 5), and therefore the need of revascularization would be debatable. As can be seen, the pre-stenotic velocity distribution does not play a significant role in this case. Figure 9 shows the slight influence of the inlet velocity distribution for \({\mathcal {S}}=0.583.\)

The stenosis eccentricity enhances turbulent motion in the post-stenotic region [21]. It is natural to wonder whether turbulence affects the pressure drop across the stenosis for the Reynolds numbers typically found in the coronary arteries. To answer this question, we have solved the LES model for the average Reynolds number Re = 1268. Figure 10 shows the instantaneous vorticity field measured in the artery symmetry plane. Vorticity vanishes at the entrance because the uniform velocity profile is imposed in that section. It takes small values in the nearly-irrotational core of both the stream in front of the stenosis and the post-stenotic jet. High vorticity values are confined within the thin boundary layer attached to the wall in the pre-stenotic region. The boundary layer separation taking place in the stenosis neck triggers the turbulent motion. The turbulent domain occupies a significant portion of the artery in the post-stenotic region. Relaminarization takes place relatively far away from the outlet. Turbulence slightly enhances energy dissipation. In fact, \(\Pi\) increases only in about 3% with respect to the laminar case. As will be explained above, increments of \(\Pi\) of this order of magnitude translate into negligible increases of the FFR value. Therefore, the laminar approximation is fully justified when it comes to the calculation of the FFR value for our stenosis model.

Figure 10 also shows the complex flow pattern behind the stenosis. The blood jet originated from the stenosis neck destabilizes, rotates and slows down in the post-stenotic region where the flow reattaches to the artery wall (Fig. 11). That blood deceleration is associated with a significant recovery of static pressure in a relatively short distance. This implies that the FFR value measured in this configuration may considerably depend on the section where the post-stenotic pressure is measured if that section is not located sufficiently far away from the stenosis.

In order to analyze the influence of the Womerslay number \(\alpha\) on the pressure drop across the stenosis, we compared the results for the pulsatile flow with those of the steady regime for the same Reynolds numbers. For this purpose, we defined the Reynolds number \(\text{Re}^*=\rho v_{0}(t) D/{\mu }_{\infty }\) and pressure drop \(\Pi^*=[P_0(t)-P_1(t)]/\rho v_{0}^2(t),\) which are the instantaneous counterparts of the corresponding dimensionless numbers in the steady regime. Figure 13 shows the comparison between the curves \(\Pi^*(\mathrm{Re}^*)\) and \(\Pi(\text{Re}).\) The results indicate that \(\Pi^*\) varies over the cardiac cycle due to both the variation of \(\text{Re}^*\) and the existence of acceleration and deceleration phases. The pressure drop takes larger values in the acceleration phase due to the extra force exerted to increasing liquid momentum, while the opposite occurs in the deceleration phase. The values for the steady regime approximately coincide with those of the pulsatile flow at the Reynolds numbers corresponding to zero acceleration at the inlet section \((dv_0/dt=0).\) The pulsatile character of the flow produces variations of around 5% in the pressure loss across the stenosis with respect to the corresponding steady regime value.

The results presented in the previous section and the order-of-magnitude analysis described above allow one to conclude that errors in the calculation of FFR made when assuming a certain velocity profile in front of the stenosis, as well as those associated with the Newtonian and laminar approximations, must be smaller than 1% for our stenosis model with \({\mathcal {S}}=0.583\) (which corresponds to \(\text{FFR}\simeq 0.8\) if the inlet velocity waveform in Fig. 5 is assumed). On the contrary, errors in the segmentation of the CTA images and the surface representation may translate into considerable deviations of the FFR value.

### Influence of the segmentation error on the fractional flow reverse

In order to show the validity of the above conclusions, we quantified the influence on the FFR value of one of the two major sources of error: the segmentation of the CTA images. Specifically, we simulated the pulsatile flow across the nominal and distorted geometries considered in “Pressure drop across the stenosis” section. In all the cases, we took \(D=5.15\) mm and \(D_s=10\) mm. The nominal geometry and distorted geometries correspond to \(\delta D_s=0\) and ± 300 μm, respectively. The latter value corresponds to 1-pixel size in Fig. 2b. The severity of the nominal geometry is 58.3%, while \({\mathcal {S}}\) = 63 and 53% for \(\delta D_s\) = + 300 and − 300 μm, respectively.

The simulations were conducted by imposing velocity inlet and outflow boundary conditions. In this case, the simulation calculates pressures relative to that of a reference point of the numerical domain. In order to obtain absolute values of the pressure field, the value of this quantity was prescribed at the inlet section. We calculated the FFR values by considering two inlet conditions: (i) the pre-stenotic velocity and pressure waveforms taken from Ref. [36], measured in a mildly stenotic coronary artery during hyperaemia, and (ii) the waveforms measured also on hyperaemia conditions when treating a coronary artery with a severe eccentric stenosis in the Hospital Infanta Cristina in Badajoz (Spain).

Based on our analysis, one can propose a two-step method to improve the CFD calculation of the FFR value for highly-eccentric coronary stenoses. In the first step, one determines the internal geometry of the large-scale numerical domain (including the aorta and the rest of peripheral vasculature) from CTA images of the cardiac vasculature, and conduct the corresponding numerical simulation. This calculation allows one to determine both the time-dependent velocity profile and pressure waveform at a given section of the pre-stenotic segment. If possible, that velocity profile might be re-calculated by scaling all the velocities so that the flow rate coincides with that measured from any relatively non-invasive technique, like MR phase-contrast flow measurements [27, 28]. In the second step, images of the stenosed coronary artery are acquired to increasing the accuracy of the image segmentation and the resolution of the surface representation in that region. For this purpose, OCT might be used [43, 44]. Then, the FFR value is determined by simulating the pulsatile flow across that region using the inlet pressure and velocity obtained in the first step. A similar idea was proposed in Ref. [19]. In that case, the pre-simulation was conducted to calculate just the pre-stenotic pressure, and the uniform velocity profile was assumed in that region.

## Conclusions

FFR-guided revascularization offers important advantages over other assessments such as CTA or AME. For instance, it proves to be long-lived and cost-saving, exhibits lower rates of major adverse cardiac events, and correctly classifies as hemodynamically insignificant many lesions judged “obstructive” by CTA or AME. However, FFR assessment requires the catheterization of the coronary arteries to measure the pressure waveforms both proximal and distal to the stenosis. CFD has arisen as a useful alternative to this invasive technique [10, 11]. It can predict the FFR value just from CTA images of the patient cardiac vasculature when some approximations are taken. The major concern when conducting CFD simulations is typically the validity of the inlet and outlet boundary conditions necessary to correctly pose the mathematical problem. This difficulty has been addressed by enlarging the numerical domain, so that it ends where certain boundary conditions can be safely imposed [10, 11]. Much less (if any) attention has been paid to the critical role played by the geometrical reconstruction and representation of the stenosis surface. It must be pointed out that the disparity between the scales of the fluid domain considered in the simulations and the critical stenotic region considerably hinders an accurate geometrical reconstruction of the latter. Our results have shown that both the limited resolution of the CTA images and a coarse representation of the artery surface constitute an important obstacle to determine the FFR value of highly-eccentric stenoses with sufficient accuracy. It must be noted that eccentric coronary plaques are prone to dynamic changes in stenosis geometry during the cardiac cycle, especially partial collapse with decreasing distending pressure during higher hyperemic velocity in diastole. These changes in geometry may well exceed the effect of imaging inaccuracy. The spatial resolution of images acquired with optical coherence tomography (OCT) may be sufficient to ensure accurate predictions for the FFR value.

We have extended our study to analyze the influence of a number of approximations on the pressure drop across highly-eccentric coronary stenoses. For this purpose, we have parameterized the stenosis geometry by adopting an idealized model that represents a wide class of stenoses caused by atheromatous lesions. We have considered realistic values of the rest of parameters governing the problem (artery diameter, velocity and pressure at the inlet section,…). The conclusions derived from this analysis are expected to be qualitatively valid for real stenoses in curved, flexible and moving vessels with side branches. Our results on the pressure drop across the stenosis confirm the common belief that the flow can be regarded as Newtonian and laminar. In addition, the inlet velocity distribution plays a secondary role in severe stenoses because of the very localized character of the pressure drop in the artery narrowing. Finally, the flow pulsatility has little influence on the pressure drop due to both the moderate value of the Womerslay number and the smooth temporal evolution over the cardiac cycle of the velocity in the coronary artery inlet.

We have restricted our analysis to the geometry described in “Formulation of the problem” section because the majority of coronary artery stenoses are caused by atheromatous lesions producing an eccentric narrowing similar that geometry [1]. A natural question is whether the above conclusions obtained for highly-eccentric coronary stenoses can also be applied to concentric ones. Most of the features analyzed here are not expected to be specific of eccentric stenoses; in particular, the fact that both the limited resolution of the CTA images and a coarse representation of the artery surface constitute important obstacles to calculate accurately the FFR value. On the contrary, the destabilization of the blood jet associated with a significant adverse pressure gradient may considerably depend on the stenosis eccentricity degree. This is an important phenomenon because it can make the FFR value very sensitive to the location where distal pressure is measured. The sample of analyzed geometries must be enlarged to complete the present study. It is also interesting to assess quantitatively the validity of the above conclusions for real patient-specific geometries and their corresponding hemodynamic conditions.

The present work raises the question of whether the image spatial resolution obtained in standard CTA is high enough for CFD to provide reliable predictions in FFR assessment. To answer this question, one must conduct a systematic study of the effects produced by choices made in the segmentation process on the resulting FFR value. In particular, attention must be paid to the range of values in the Hounsfield scale corresponding to the area delimited by the artery inner wall.

In principle, the FFR value could be also determined by conducting experiments on models fabricated with 3D printing [46, 47]. Increasing the size of the model does not reduce the error coming from image segmentation, but does decrease that associated with a poor representation of the artery surface. In fact, the use of big models allows one to increase the number *N* of triangles of the faceted geometry while keeping the ratio \(S/D^2\) constant, which results in higher values of SRI. In this sense, high-resolution printers have proved to be capable of forming precise replicas of big models of coronary arteries [47].

## Declarations

### Authors’ contributions

RA conducted the numerical simulations. MRG-F and JMN-A participated in the analysis and elaboration of the results. JMM participated in the elaboration of the results and drafted the manuscript. All authors read and approved the final manuscript.

### Acknowledgements

Partial support from the Junta de Extremadura through Grant No. IB16119 (partially financed by FEDER funds) is gratefully acknowledged.

### Competing interests

The authors declare that they have no competing interests.

### Ethics approval and consent to participate

Not applicable.

### Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

**Open Access**This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

## Authors’ Affiliations

## References

- Klein L. Coronary stenosis morphology: analysis and implication. New York: Springer; 1997.View ArticleGoogle Scholar
- Xin L, Zhenzhou L, Lina P, Zhifan G, WenHua H, Heye Z, Yuan-Ting Z. Numerical study of the mechanical interaction between the carotid atherosclerosis plaque and the blood flow field distribution: a geometric assessment of the plaque vulnerability. J Med Imaging Health Inform. 2015;5:1–12.View ArticleGoogle Scholar
- Meijboom WB. Comprehensive assessment of coronary artery stenoses: computed tomography coronary angiography versus conventional coronary angiography and correlation with fractional flow reserve in patients with stable angina. J Am Coll Cardiol. 2008;52:636–43.View ArticleGoogle Scholar
- Pijls NH, van Son JA, Kirkeeide RL, Bruyne BD, Gould KL. Experimental basis of determining maximum coronary, myocardial, and collateral blood flow by pressure measurements for assessing functional stenosis severity before and after percutaneous transluminal coronary angioplasty. Circulation. 1993;87:1354–67.View ArticleGoogle Scholar
- van de Hoef TP, Meuwissen M, Escaned J, Davies JE, Siebes M, Spaan JE, Piek JJ. Fractional flow reserve as a surrogate for inducible myocardial ischaemia. Nat Rev Cardiol. 2013;10:439–52.View ArticleGoogle Scholar
- Fearon WF, Bornschein B, Tonino PA, Gothe RM, Bruyne BD, Pijls NH, Siebert U. Economic evaluation of fractional flow reserve-guided percutaneous coronary intervention in patients with multivessel disease. Circulation. 2010;122:2545–50.View ArticleGoogle Scholar
- Tonino PAL. Angiographic versus functional severity of coronary artery stenoses in the FAME study. J Am Coll Cardiol. 2010;55:2816–21.View ArticleGoogle Scholar
- Taylor CA, Figueroa CA. Patient-specific modeling of cardiovascular mechanics. Annu Rev Biomed Eng. 2009;11:109–34.View ArticleGoogle Scholar
- Kim HJ, Vignon-Clementel IE, Coogan JS, Figueroa CA, Jansen KE, Taylor CA. Patient-specific modeling of blood flow and pressure in human coronary arteries. Ann Biomed Eng. 2010;38:3195–209.View ArticleGoogle Scholar
- Koo BK. Diagnosis of ischemia-causing coronary stenoses by noninvasive fractional flow reserve computed from coronary computed tomographic angiograms. J Am Coll Cardiol. 2011;58:1989–97.View ArticleGoogle Scholar
- Taylor CA, Fonte TA, Min JK. Computational fluid dynamics applied to cardiac computed tomography for noninvasive quantification of fractional flow reserve. J Am Coll Cardiol. 2013;61:2233–41.View ArticleGoogle Scholar
- Chiastra C, Morlacchi S, Gallo D, Morbiducci U, Cárdenes R, Larrabide I, Migliavacca F. Computational fluid dynamic simulations of image-based stented coronary bifurcation models. J R Soc Interface. 2013;10:20130193.View ArticleGoogle Scholar
- Vorobstova N, Chiastra C, Stremler MA, Sane DC, Migliavacc F, Vlachos P. Effects of vessel tortuosity on coronary hemodynamics: an idealized and patient-specific computational study. Ann Biomed Eng. 2016;44:2228–39.View ArticleGoogle Scholar
- Morris PD, Ryan D, Morton AC, Lycett R, Lawford PV, Hose DR, Gunn JP. Virtual fractional flow reserve from coronary angiography: modeling the significance of coronary lesions: results from the VIRTU-1 (VIRTUal fractional flow reserve from coronary angiography) study. JACC Cardiovasc Interv. 2013;6:149–57.View ArticleGoogle Scholar
- Ghalichi F, Deng X, Champlain A, Douville Y, King M, Guidoin R. Low Reynolds number turbulence modeling of blood flow in arterial stenoses. Biorheology. 1998;35:281–94.View ArticleGoogle Scholar
- Varghese SS, Frankel SH. Numerical modeling of pulsatile turbulent flow in stenotic vessels. J Biomech Eng. 2003;125:445–60.View ArticleGoogle Scholar
- Varghese SS, Frankel SH, Fischer PF. Direct numerical simulation of stenotic flows. Part 1. Steady flow. J Fluid Mech. 2007;582:253.MathSciNetView ArticleMATHGoogle Scholar
- Varghese SS, Frankel SH, Fischer PF. Direct numerical simulation of stenotic flows. Part 2. Pulsatile flow. J Fluid Mech. 2007;582:281.MathSciNetView ArticleMATHGoogle Scholar
- Sommer K, Bernat D, Schmidt R, Breit H, Schreiber LM. Resting myocardial blood flow quantification using contrast-enhanced magnetic resonance imaging in the presence of stenosis: a computational fluid dynamics study. Med Phys. 2015;42:4375–84.View ArticleGoogle Scholar
- Wong KKL, Wang D, Ko JKL, Mazumdar J, Le TT, Ghista D. Computational medical imaging and hemodynamics framework for functional analysis and assessment of cardiovascular structures. BioMed Eng Online. 2017;16(35):1–25.Google Scholar
- Varghese SS, Frankel SH, Fischer PF. Modeling transition to turbulence in eccentric stenotic flows. J Biomech Eng. 2008;130:014503.View ArticleGoogle Scholar
- Tan FPP, Wood NB, Tabor G, Xu XY. Comparison of LES of steady transitional flow in an idealized stenosed axisymmetric artery model with a RANS transitional model. J Biomech Eng. 2011;133:051001.View ArticleGoogle Scholar
- Bahramian F, Mohammadi H. Modeling blood flow in an eccentric stenosed artery using large eddy simulation and parallel computing. J Mech Med Biol. 2015;15:1550086.View ArticleGoogle Scholar
- Wong KKL, Thavornpattanapong P, Cheung SCP, Tu JY. Biomechanical investigation of pulsatile flow in a three-dimensional atherosclerotic carotid bifurcation model. J Mech Med Biol. 2013;13:1350001.View ArticleGoogle Scholar
- Liu G, Wu J, Huang W, Wu W, Zhang H, Wong KKL, Ghista DN. Numerical simulation of flow in curved coronary arteries with progressive amounts of stenosis using fluid-structure interaction modelling. J Med Imaging Health Inform. 2014;4:605–11.View ArticleGoogle Scholar
- Khair A, Wang B, Kuhn DCS. Study of laminar-turbulent flow transition under pulsatile conditions in a constricted channel. Int. J Comput Fluid Dyn. 2015;29:447–63.MathSciNetView ArticleGoogle Scholar
- Sakuma H, Saeed M, Takeda K, Wendland MF, Schwitter J, Szolar DH, Derugin N, Shimakawa A, Foo TKE, Higgins CB. Quantification of coronary artery volume flow rate using fast velocity-encoded cine mr imaging. AJR. 1997;168:1363–7.View ArticleGoogle Scholar
- Nagel E, Thouet T, Klein C, Schalla S, Bornstedt A, Schnackenburg B, Hug J, Wellnhofer E, Fleck E. Noninvasive determination of coronary blood flow velocity with cardiovascular magnetic resonance in patients after stent deployment. Circulation. 2003;107:1738–43.View ArticleGoogle Scholar
- Huo Y, Svendsen M, Choy JS, Zhang Z-D, Kassab GS. A validated predictive model of coronary fractional flow reserve. J R Soc Interface. 2012;9:1325–38.View ArticleGoogle Scholar
- Pietrabissa R, Mantero S, Marottat T, Menicantit L. A lumped parameter model to evaluate the fluid dynamics of different coronary bypasses. Med Eng Phys. 1996;18:477–84.View ArticleGoogle Scholar
- Mohammadi H, Bahramian F. Boundary conditions in simulation of stenosed coronary arteries. Cardiovasc Eng. 2009;9:83–91.View ArticleGoogle Scholar
- van der Giessen AG, et al. The influence of boundary conditions on wall shear stress distribution in patients specific coronary trees. J Biomech. 2011;44:1089–95.View ArticleGoogle Scholar
- Mates RE, Gupta RL, Bell AC, Klocke FJ. Fluid dynamics of coronary artery stenosis. Circ Res. 1978;42:152–62.View ArticleGoogle Scholar
- Kim Y, Park Y, Lim S. 3D simulations of blood flow dynamics in compliant vessels: normal, aneurysmal, and stenotic arteries. Commun Comput Phys. 2016;19:1167–90.MathSciNetView ArticleMATHGoogle Scholar
- Womersley JR. Method for the calculation of velocity, rate of flow and viscous drag in arteries when the pressure gradient is known. J Physiol. 1955;127:553–63.View ArticleGoogle Scholar
- Spaan J, Kolyva C, van Wijngaard J, Wee RT, van Horssen P, Piek J, Siebes M. Coronary structure and perfusion in health and disease. Philos Trans R Soc A. 2008;366:3137–53.View ArticleGoogle Scholar
- Siebes M, Verhoeff BJ, Meuwissen M, de Winter RJ, Spaan JAE, Piek JJ. Single-wire pressure and flow velocity measurement to quantify coronary stenosis hemodynamics and effects of percutaneous interventions. Circulation. 2004;109:756–62.View ArticleGoogle Scholar
- de Bruyne B, Bartunek J, Sys SU, Pijls NHJ, Heyndrickx GR, Wijns W. Simultaneous coronary pressure and flow velocity measurements in humans. Circulation. 1996;94:1842–9.View ArticleGoogle Scholar
- Marques KMJ, van Eenige MJ, Spruijt HJ, Westerhof N, Twisk J, Visser CA, Visser FC. The diastolic flow velocity-pressure gradient relation and dpv50 to assess the hemodynamic significance of coronary stenoses. Am J Physiol Heart Circ Physiol. 2006;291:2630–5.View ArticleGoogle Scholar
- Versteeg HK, Malalasekera W. An introduction to computational fluid dynamics. England: Pearson Education Limited; 2007.Google Scholar
- Ansys. Ansys Fluent 16.2. User’s guide. Canonsburg: Ansys; 2015.Google Scholar
- Lupi A, Buffon A, Finocchiaro ML, Conti E, Maseri A, Crea F. Mechanisms of adenosine-induced epicardial coronary artery dilatation. Eur Heart J. 1997;18:614–7.View ArticleGoogle Scholar
- Kousera CA, Nijjer S, Torii R, Petraco R, Sen S, Foin N, Hughes AD, Francis DPP, Xu XY, Davies JE. Patient-specific coronary stenoses can be modeled using a combination of OCT and flow velocities to accurately predict hyperemic pressure gradients. IEEE Trans Biomed Eng. 2014;61:1902–13.View ArticleGoogle Scholar
- Migliori S, Chiastra C, Bologna M, Montin E, Dubini G, Aurigemma C, Fedele R, Burzotta F, Mainardi L, Migliavacca F. A framework for computational fluid dynamics analyses of patient-specific stented coronary arteries from optical coherence tomography images. Med Eng Phys. 2017;47:105–16.View ArticleGoogle Scholar
- Park J et al. Clinical relevance of functionally insignificant moderate coronary artery stenosis assessed by 3-vessel fractional flow reserve measurement. J Am Heart Assoc. 2018. https://doi.org/10.1161/JAHA.117.008055.Google Scholar
- Agujetas R, Ferrera C, Marcos AC, Alejo JP, Montanero JM. Numerical and experimental analysis of the transitional flow across a real stenosis. Biomech Model Mechanobiol. 2017;16:1447–58.View ArticleGoogle Scholar
- Yang Y, Liu X, Xia Y, Liu X, Wu W, Xiong H, Zhang H, Xu L, Wong KKL, Ouyang H. Impact of spatial characteristics in the left stenotic coronary artery on the hemodynamics and visualization of 3d replica models. Sci Rep. 2017;7:15452.View ArticleGoogle Scholar