- Research
- Open Access
- Published:

# A new CFD based non-invasive method for functional diagnosis of coronary stenosis

*BioMedical Engineering OnLine*
**volume 17**, Article number: 36 (2018)

## Abstract

### Background

Accurate functional diagnosis of coronary stenosis is vital for decision making in coronary revascularization. With recent advances in computational fluid dynamics (CFD), fractional flow reserve (FFR) can be derived non-invasively from coronary computed tomography angiography images (FFR_{CT}) for functional measurement of stenosis. However, the accuracy of FFR_{CT} is limited due to the approximate modeling approach of maximal hyperemia conditions. To overcome this problem, a new CFD based non-invasive method is proposed.

### Methods

Instead of modeling maximal hyperemia condition, a series of boundary conditions are specified and those simulated results are combined to provide a pressure-flow curve for a stenosis. Then, functional diagnosis of stenosis is assessed based on parameters derived from the obtained pressure-flow curve.

### Results

The proposed method is applied to both idealized and patient-specific models, and validated with invasive FFR in six patients. Results show that additional hemodynamic information about the flow resistances of a stenosis is provided, which cannot be directly obtained from anatomy information. Parameters derived from the simulated pressure-flow curve show a linear and significant correlations with invasive FFR (r > 0.95, P < 0.05).

### Conclusion

The proposed method can assess flow resistances by the pressure-flow curve derived parameters without modeling of maximal hyperemia condition, which is a new promising approach for non-invasive functional assessment of coronary stenosis.

## Background

Functionally significant stenosis generally causes angina symptoms, and is associated with inducible ischemia and impaired outcome. Therefore, it should be revascularized. While, if a stenosis has no functional significance, medical treatment is quite excellent [1]. Therefore, determining the functional significance of a coronary stenosis plays a pivotal role in decision making in coronary revascularization [1]. Fractional flow reserve (FFR) is considered as the gold standard for assessment of functional measurement of stenosis [2, 3]. It is a pressure derived index and defined as the ratio between the distal and proximal pressure of the stenosis at maximal hyperemia [4]. A cut-off value of 0.75–0.8 is employed to indicate the functionally significant stenosis [3]. FFR-guided strategy has been proved to be safety and also has been demonstrated to be both cost-effective and cost-saving [5,6,7]. However, FFR is an invasive procedure and requires pharmacologic intervention to induce maximal hyperemia, which limits its in-hospital utilization [2, 8]. To overcome this problem, a novel technology that combines coronary computed tomography angiography (cCTA) and computational fluid dynamics (CFD) has been developed [9,10,11]. By using the CFD method, an accurate reconstruction of coronary flow and pressure fields can be obtained from cCTA images, and then, the coronary diagnostic index (CT-derived computed FFR, FFR_{CT}) can be derived without additional medications [9, 12,13,14,15]. Multiple clinic trials have demonstrated that the performance of FFR_{CT} was superior to cCTA stenosis for diagnosing ischemic lesions [16,17,18,19,20,21]. The use of FFR_{CT} is showed to reduce the overall use of invasive angiography and more and more researchers suggest that FFR_{CT} can be the gatekeeper to the cardiac catheterization laboratory [13].

Pressure drop across a stenosis can be approximately determined by a common fluid dynamic equation [22, 23]:

where \(\Delta \overline{p}\) is the mean pressure drop, *f* the is the viscous friction, *s* is the expansion loss and \(\overline{Q}\) is the mean flow rate. Coronary stenosis will increase the viscous friction and expansion loss of the stenosis section, leading to an increase in pressure drop. The pressure derived index (FFR and FFR_{CT}) quantify the functionally significant of stenosis by using the pressure drop. However, the pressure drop relies on the flow rate, which is also determined by the behavior of distal vascular trees (micro-vascular resistance) [2, 23]. Therefore, to exclude the influence of the flow rate, FFR and FFR_{CT} should be measured during maximal hyperemia. Failure to achieve maximal hyperemia would result in an inaccuracy of FFR [1]. Additional, it will be confounded by the presence of microvascular disease [1]. Several FFR based hyperemia-free indices (basal FFR and instantaneous wave-free ratio) have been proposed for detecting coronary stenosis under basal conditions [24, 25]. However, similar to FFR, the flow rate is still an uncontrollable factor which has great impact on these indices. Since the main principle of FFR_{CT} is to model the blood flow in coronary during the hyperemia condition [9], an accurate modeling of the hyperemia condition is critical. Several physiological models are used in the modeling process: to obtain the baseline coronary flow, a fixed relationship between the baseline flow and the left ventricular myocardial mass is assumed; to get the microvascular resistance, a fixed relationship between the resistance in baseline and hyperemia condition is assumed [9]. These physiological models provide a general approach to model the hyperemia flow rate and make it possible to compute FFR_{CT} without any other medications. However, these models reflect the average behaviors in coronary circulation, ignoring the individual difference. Obviously, these assumptions would reduce the reliability in modeling the individual flow rate during hyperemia condition, which further degrades the accuracy of FFR_{CT}.

To overcome this problem, a new CFD based non-invasive approach for functional measurement of stenosis is proposed. It is based on an invasive non-dimensional index, pressure drop coefficient (CDP), which combining both the pressure and flow velocity information. CDP is defined as the ratio of trans-stenotic pressure drop to distal dynamic pressure during maximal hyperemia [26]. Compared with pressure-derived indices (FFR), CDP directly quantifies the hemodynamic behavior of a stenosis section [26]. Thus, the presence of microvascular disease has limited impact on CDP [27,28,29]. However, due to the limitation of the invasive procedure (the viscous loss is not considered in order to get a result with a single measurement; the flow velocity is measured, but not the flow rate), CDP is still dependent on the flow rate across a stenosis [30]. Although previous studies reported that CDP were useful during both basal and hyperemia condition, the cutoff values were quite different in those two conditions [27, 30]. Additional, guide wire insertion has great impact on flow patterns, and reliable measurement of flow velocity is also technically difficult [23]. Fortunately, those limitations of CDP can be overcame by combining with CFD method. The viscous loss can be included with multiple measurements, simply by specifying different boundary conditions; flow rates also can be accurately obtained without the influence of guide wire insertion. Unlike FFR_{CT}, the proposed approach do not need to model the hyperemia condition. Instead, the stenosis flow with a series of flow rates will be simulated and the pressure-flow relationships is obtained. Then, *f* and *s* in Eq. 1 can be estimated and further employed to functional measurement of stenosis.

## Methods

### Vascular geometry reconstruction

In total, 20 3D coronary models with stenosis (one right coronary artery (RCA), one left circumflex coronary artery (LCX) and 18 left anterior descending coronary arteries (LAD)) have been reconstructed from the cCTA images of 19 patients. The cCTA images are acquired in mid-diastole (60–70% of R–R interval) or systole (40–50% of R–R interval) with 0.4-mm slice interval and 0.28 mm/pixel in-plane image resolution, using a second-generation DSCT (SOMATOM Definition Flash, Siemens Healthcare, Forchheim, Germany) scanner. Coronary arteries are semi-automatically segmented and reconstructed in Mimics (Mimics, Materialise, Leuven, Belgium). And then, only the stenosis sections with nearby branches are retained for further CFD analysis (as shown in Fig. 1a). Area stenosis (AS%) and minimal lumen area (MLA) are calculated based on the reconstructed coronary artery model. AS% is defined as 100% minus the percentage of minimal lumen area to the reference lumen area. The reference site is in close proximity to the lesion, without intervening branch vessels.

Besides the reconstructed patient-specific models, several idealized coronary stenosis models are also created. As shown in Fig. 1b, two types of stenosis are modeled: concentric shaped stenosis and eccentric shaped stenosis. For each types, seven models with the AS% of 65, 70, 75, 80, 85, 90 and 95% are created. These idealized stenosis sections are modeled as 3D pipes with a diameter of 3 mm and the stricture length is kept to be two diameters for all models.

### Boundary conditions

Instead of providing approximate boundary conditions for pulsatile flow simulation (as FFR_{CT}), a series of boundary conditions with specified pressures and resistances are applied for steady flow simulation. In the proposed method, a static pressure (88 mmHg) is applied to the inlet. For each outlets, a lumped parameter model with only one resistance is coupled (Fig. 1c). The resistance represents the microvascular resistance. The “form–function” relationship (which relates the resistance of a downstream vessel to the vessel size at each outlet) is assumed to derive the values of microvascular resistance for each outlets with a given total distal microvascular resistance [9, 31]. Finally, a series values of total distal microvascular resistance are specified, providing different outlet boundary conditions for each steady flow simulation. The total distal microvascular resistance is set to be 240 (mmHg s/cm^{3}) at first, and then reduced to 87.5, 75.0, 62.5, 50.0, 37.5 and 25.0%. Thus, for one case, steady flow is simulated seven times with seven different total distal microvascular resistances. Since only one outlet exists in each idealized models, the outlet may not be coupled to the lumped parameter model to enforce the “form–function” relationship. Thus, a more simple boundary conditions can be applied to get the pressure-flow curves. For idealized models, a relative zero pressure is assumed at the outlet and a series of static pressure is applied at the inlet. The inlet static pressure is set to be 0.5 mmHg at first, and then increased by 0.5 mmHg each time, until it reaches 11 mmHg.

### Computational method

ANSYS FLUENT V14 (ANSYS Inc.) was used to perform flow simulations. The blood was modeled as an incompressible Newtonian fluid. The dynamic viscosity is set to be 3.5 cP and the density to be 1050 kg/m^{3}. The k–ω shear stress transport turbulence model was adopted for turbulence modeling of the low Reynolds number flow in stenosis arteries [32]. Mesh independence was judged by comparing both the computed velocities and the pressure. For each case, further grid refinement (doubled mesh resolution) led to < 1% relative error in velocity and pressure profiles.

### Post-processing

The start and end of a lesion in curved multiplanar reformats are defined by an experienced observer. Then, for each case, the mean pressures at these two planars and the mean flow rate across the stenosis are extracted from each steady flow simulation results. The pressure drop is obtained by subtracting the mean pressure at the end of a lesion from that at the start of a lesion. Finally, *f* and *s* in Eq. 1 are estimated from the Pressure-Flow data by using iterative least squares estimation for nonlinear regression [33]. Two additional parameters are further extracted based on the estimated *f* and *s*. The area (*S*) below the pressure-flow curve (area between flow rate of 0 ml/s to *q* ml/s) can be obtained from the following equation:

with setting *q* to be 1 and 2 ml/s, respectively, two areas, S_{1} and S_{2}, can be calculated from Eq. 2.

### Validation with FFR

To further validate the proposed method, six patients who undergo both FFR and cCTA scans are employed. FFRs are measured with pressure wire (St. Jude Medical, Inc.) by experienced invasive cardiologists. cCTA images are acquired and six 3D coronary models with stenosis in LAD are reconstructed as described in “Vascular Geometry Reconstruction”. Parameters derived from the simulated pressure-flow curve are compared with the invasive FFR.

## Results

The pressure-flow relationships for idealized models are shown in Fig. 2. As the flow rate increases, the pressure drops across the stenosis section increases for both concentric and eccentric models, and a larger AS% will lead to a more rapid increase in pressure drop. Compared with concentric models, the eccentric models with a same AS% always have a larger pressure drop at given flow rate, and this phenomena will be more pronounced in models with smaller AS %. For each idealized models, the *f* and *s* in Eq. 1 are estimated from the simulated pressure drops and flow rates. From the results, the model-predicted pressure drops (by using Eq. 1 with estimated *f* and *s*) are consistent well with the simulated ones.

The estimated parameters (*f* and *s*) for idealized models with different AS% are shown in Fig. 3. Both *f* and *s* will increase as AS% increases, and a more rapid increase is observed as AS% become greater than 80%. Besides the degree of the stenosis, the shapes also have great impact on the values of the estimated parameters. For models with the same AS%, the estimated *f* and *s* for eccentric stenosis are always greater than that for concentric stenosis.

The pressure-flow relationships for reconstructed patient-specific models are shown in Fig. 4. The 20 models are separated into three groups based on AS%: AS% > 70% (severe), 50% < AS% < 69% (middle) and AS% < 49% (mild). However, no obvious differences in pressure-flow relationships are observed between severe stenosis and middle stenosis. Flat pressure-flow curves are observed in two cases in severe group, while sharp pressure-flow curves are observed in two cases in middle group. Five cases in middle group have a similar pressure-flow curves as those in mild group. For all cases, the model-predicted pressure drops (using seven steady flow simulation results) are consistent well with the simulated ones.

Figure 5 shows the linear correlation of the estimated parameters (*f* and *s*) with AS% and MLA. When correlated with AS%, only *s* show a moderate but significant correlation (r = 0.66, P < 0.05). However, when correlated with MLA, both *f* and *s* have linear and significant correlations (*f*: r = 0.84, P < 0.05; *s*: r = 0.88, P < 0.05).

Figure 6 shows the linear correlation of the extracted parameters (S_{1} and S_{2}) with AS% and MLA. For S_{1}, it has a moderate but significant correlation with AS% (r = 0.56, P < 0.05) and a linear and significant correlation with MLA (r = 0.90, P < 0.05). Similar results are observed for S_{2} (correlated with AS%: r = 0.61, P < 0.05; correlated with MLA: r = 0.90, P < 0.05).

Figure 7 shows the linear correlation of the estimated parameters with invasive FFR. All the four parameters show a linear and significant correlations with invasive FFR (*f*: r = 0.96, P < 0.05; *s*: r = 0.95, P < 0.05; *S*_{1}: r = 0.96, P < 0.05; *S*_{2}: r = 0.96, P < 0.05).

## Discussion

A promising noninvasive method was proposed for assessing of hemodynamic significance of coronary stenosis. Anatomy information of a stenosis is first derived from the cCTA images; and then, the pressure-flow relationship of a stenosis is obtained by using CFD method. Unlike previous invasive (FFR and CDP) or noninvasive method (FFR_{CT}) which evaluating the stenosis at a specific flow rate (maximum hyperemia), the proposed method characters the hemodynamic stenosis at a series of flow rate. The main advantages of the proposed method are: (1) it is a non-invasive functional diagnosis method; (2) do not need to model the maximum hyperemia condition; (3) fully characters the pressure-flow behavior of a stenosis within normal physiological range. The proposed method is validate in six patients who undergo both FFR and cCTA scans, and all the proposed four parameters have linear and significant correlations with invasive FFR.

CFD is always combined with cardiovascular imaging for blood flow analysis in various cardiac disease [34,35,36,37,38,39]. However, it is usually quite time-consuming for clinical use due to the complex computational models. To overcome this problem, steady flow simulation is used in our proposed method. And the two parameters, *f* and *s*, which are estimated with seven steady flow simulation results, can be estimated even with only two steady flow simulation results. The root-mean-square errors (RMSE) of the predicted pressure drops for models fitted with different numbers of selected steady flow simulation results are shown in Fig. 8. It is demonstrated that the pressure-flow curves can be fitted well with much less computation time (only needs two steady flow simulation results).

### Compared with cCTA

It was reported that anatomic assessments of stenosis evaluated by cCTA had significant correlation with those evaluated by digital subtraction angiography [40]. Although cCTA alone can accurately detect anatomically obstructive coronary stenosis, it cannot define the hemodynamic significance of a stenosis due to the complex relationship between stenosis and flow [1]. The spatial characteristics of stenosis have great impact on flow [41,42,43]. Even for stenosis with similar AS%, the pressure-flow behaviors still can be different due to the different shapes of stenosis (as shown in Figs. 2 and 4). With the anatomic models obtained by cCTA, a further CFD analysis is performed to fill the gap between anatomically obstructive and hemodynamic significance. In the proposed method, pressure-flow curves are obtained and further employed to define the hemodynamic significance. It characters the nonlinear relationship between the flow rate and the pressure drop across a stenosis, reflecting the flow resistances of a given stenosis. A sharper curve represents a larger flow resistance, while more flat curve means smaller flow resistance. As shown in Fig. 4, severe stenosis defined by cCTA can have a relative smaller flow resistance, but it cannot be known before the pressure-flow curve obtained. Several parameters are derived from the curves to quantify the flow resistance. As compared with AS%, flow resistance has a more direct relationship with the radius of a stenosis [44]. Thus, from the results, moderate correlations are observed when correlated these parameters with AS%, but linear and significant correlations are obtained when correlated with MLA.

### Compared with FFR_{CT}

CFD is applied in both FFR_{CT} and proposed method to provide the hemodynamic information of a stenosis. For FFR_{CT}, pressure drop during maximal hyperemia is employed [9]; while for our method, pressure-flow relationship is derived. With a given flow rate, the pressure drop across the stenosis can be determined from its pressure-flow relationship. Thus, FFR_{CT} can be derived from the obtained pressure-flow relationship if the exact flow rate during maximal hyperemia is known. However, the main problem for CFD based method is that there is no additional measurement about the flow rate during maximal hyperemia [17]. Although clinical trials (DISCOVER-FLOW, DeFACTO and NXT) demonstrated that FFR_{CT} improved diagnostic accuracy in differentiating ischemic stenosis as compared with cCTA alone [19,20,21], the performance is still unsatisfactory when compared with invasive FFR. One possible reason is that FFR_{CT} has to employ several physiological models to estimate an approximate maximal hyperemia condition. Instead of modeling maximal hyperemia condition, a series of boundary conditions are specified and those simulated results are combined to provide a pressure-flow curve in proposed method. By this way, the accuracy do not rely on the modeling of maximal hyperemia condition. Additionally, pressure-flow relationship represents the flow resistances of a stenosis, which is determined mainly by its anatomical morphology. A functionally significant stenosis will increase the flow resistances and lead to a decrease in distal perfusion pressure, further cause ischemia [23, 43]. Thus, flow resistance maybe a more essential index for functional measurement of stenosis, as compared with pressure derived index.

### Compared with CDP

CDP is a parameter which incorporates both pressure and flow variations [26]. By ignoring the viscous friction effect and replacing the flow rate as flow velocity in Eq. (1), the expansion loss (*s*) can be approximately obtained with a single measurement of pressure drop and flow velocity. In the presence of microvascular dysfunction and submaximal hyperemia, both the pressure drop and the flow velocity will decreased, leading to a significant increase of pressure-derived parameters (FFR and FFR_{CT}), but it has a limited impact on CDP [27]. Ideally, the flow rate should have little impact on CDP. However, due to the assumptions made by CDP, the measured values in basal and hyperemic conditions are still quite different [30]. Additionally, it is also difficult to measure the flow velocity accurately. Fortunately, those problem can be overcome by the proposed CFD based method. By combining several measurements, both the viscous friction (*f*) and expansion loss (*s*) can be estimated and the pressure-flow relationship can be obtained. With the CFD method, the flow rate also can be measured without the impact of guide wire insertion. It is reported that CDP has the ability to localize the differential culpability for flow impediment: translesional or microvascular [27,28,29]. But in proposed method, it is noninvasive and only the anatomy information derived from cCTA is used, thus the microvascular disease is excluded and it is concentrated to evaluate the culpability for flow impediment of translesional stenosis.

## Limitations

In this paper, a CFD based non-invasive method is proposed for functional measurement of coronary stenosis, which do not need to model the maximal hyperemia condition. However, only six patients are employed to validate with current “gold standard” invasive FFR. Besides, due to lack of FFR data, currently there is no cut-off value for those parameters derived from pressure-flow curves. Even though, this is a first approach to functional measurement of coronary stenosis non-invasively without modeling the maximal hyperemia conditions. Given this advantage over FFR_{CT}, we believe that it will have a promising role in clinical practice. Still, further validations need to be performed, especially corrected with FFR. To limit the computational time, steady state flow is simulated to get the pressure-flow curves. However, blood flow is pulsating, which is inconsistence with steady state flow assumption. Fortunately, current “gold standard” invasive FFR is based on time-averaged pressure measured over several cardiac cycles [4], and previous studies also demonstrated that FFR_{CT} obtained by using steady state flow simulation still performed well [45]. Additionally, the boundary conditions used for pulsating flow simulation in FFR_{CT} is derived from lumped parameter models, which is also inconsistence with real conditions [11, 46]. Thus, steady state flow is simulated in our method.

## References

- 1.
Pijls NHJ, Sels J-WEM. Functional measurement of coronary stenosis. J Am Coll Cardiol. 2012;59(12):1045–57.

- 2.
Koo B-K. The present and future of fractional flow reserve. Circ J. 2014;78(5):1048–54.

- 3.
Pijls NH, van Schaardenburgh P, Manoharan G, Boersma E, Bech JW, van’t Veer M, Bar F, Hoorntje J, Koolen J, Wijns W, de Bruyne B. Percutaneous coronary intervention of functionally nonsignificant stenosis: 5-year follow-up of the DEFER Study. J Am Coll Cardiol. 2007;49(21):2105–11.

- 4.
Pijls NH, van Son JA, Kirkeeide RL, De Bruyne B, 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(4):1354–67.

- 5.
Elgendy IY, Conti CR, Bavry AA. Fractional flow reserve: an updated review. Clin Cardiol. 2014;37(6):371–80.

- 6.
De Bruyne BPN, Kalesan B, Barbato E, Tonino PA, Piroth Z, et al. Fractional flow reserve-guided pci versus medical therapy in stable coronary disease. N Engl J Med. 2012;367(18):1768.

- 7.
Tonino PAL, De Bruyne B, Pijls NHJ, Siebert U, Ikeno F, Veer MVT, Klauss V, Manoharan G, Engstrom T, Oldroyd KG. Fractional flow reserve versus angiography for guiding percutaneous coronary intervention. N Engl J Med. 2009;360(3):213–24.

- 8.
Pothineni NV, Shah NS, Rochlani Y, Nairooz R, Raina S, Leesar MA, Uretsky BF, Hakeem A. US trends in inpatient utilization of fractional flow reserve and percutaneous coronary intervention. J Am Coll Cardiol. 2016;67(6):732–3.

- 9.
Taylor CA, Fonte TA, Min JK. Computational fluid dynamics applied to cardiac computed tomography for noninvasive quantification of fractional flow reserve: scientific basis. J Am Coll Cardiol. 2013;61(22):2233–41.

- 10.
Serruys PW, Girasis C, Papadopoulou SL, Onuma Y. Non-invasive fractional flow reserve: scientific basis, methods and perspectives. EuroIntervention. 2012;8(4):511–9.

- 11.
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(10):3195–209.

- 12.
Norgaard BL, Leipsic J. From Newton to the coronaries computational fluid dynamics has entered the clinical scene. JACC Cardiovasc Imaging. 2016;9(6):700–2.

- 13.
Hecht HS, Narula J, Fearon WF. Fractional flow reserve and coronary computed tomographic angiography a review and critical analysis. Circ Res. 2016;119(2):300–16.

- 14.
Min JK, Taylor CA, Achenbach S, Koo BK, Leipsic J, Norgaard BL, Pijls NJ, De Bruyne B. Noninvasive fractional flow reserve derived from coronary ct angiography: clinical data and scientific principles. JACC Cardiovasc Imaging. 2015;8(10):1209–22.

- 15.
Zarins C, Taylor C, Min J. Computed fractional flow reserve (FFTCT) derived from coronary CT angiography. J Cardiovasc Trans Res. 2013;6(5):708–14.

- 16.
Kruk M, Wardziak L, Demkow M, Pleban W, Pregowski J, Dzielinska Z, Witulski M, Witkowski A, Ruzyllo W, Kepka C. Workstation-based calculation of CTA-based FFR for intermediate stenosis. JACC Cardiovasc Imaging. 2016;9(6):690–9.

- 17.
De Geer J, Sandstedt M, Bjorkholm A, Alfredsson J, Janzon M, Engvall J, Persson A. Software-based on-site estimation of fractional flow reserve using standard coronary CT angiography data. Acta Radiol. 2016;57(10):1186–92.

- 18.
Renker M, Schoepf UJ, Wang R, Meinel FG, Rier JD, Bayer RR 2nd, Mollmann H, Hamm CW, Steinberg DH, Baumann S. Comparison of diagnostic value of a novel noninvasive coronary computed tomography angiography method versus standard coronary angiography for assessing fractional flow reserve. Am J Cardiol. 2014;114(9):1303–8.

- 19.
Norgaard BL, Leipsic J, Gaur S, Seneviratne S, Ko BS, Ito H, Jensen JM, Mauri L, De Bruyne B, Bezerra H, et al. Diagnostic performance of noninvasive fractional flow reserve derived from coronary computed tomography angiography in suspected coronary artery disease: the NXT trial (Analysis of Coronary Blood Flow Using CT Angiography: Next Steps). J Am Coll Cardiol. 2014;63(12):1145–55.

- 20.
Leipsic J, Yang T-H, Thompson A, Koo B-K, Mancini GBJ, Taylor C, Budoff MJ, Park H-B, Berman DS, Min JK. CT angiography (CTA) and diagnostic performance of noninvasive fractional flow reserve: results from the determination of fractional flow reserve by anatomic CTA (DeFACTO) study. Am J Roentgenol. 2014;202(5):989–94.

- 21.
Min JK, Koo B-K, Erglis A, Doh J-H, Daniels DV, Jegere S, Kim H-S, Dunning A, Defrance T, Leipsic J. Effect of image quality on diagnostic accuracy of noninvasive fractional flow reserve: results from the prospective multicenter international DISCOVER-FLOW study. J Cardiovasc Comput Tomogr. 2012;6(3):191–9.

- 22.
Gould KL. Pressure-flow characteristics of coronary stenoses in unsedated dogs at rest and during coronary vasodilation. Circ Res. 1978;43(2):242–53.

- 23.
Banerjee RK, Roy AS, 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–62.

- 24.
Sen S, Escaned J, Malik IS, Mikhail GW, Foale RA, Mila R, Tarkin J, Petraco R, Broyd C, Jabbour RJ. Development and validation of a new adenosine-independent index of stenosis severity from coronary wave-intensity analysis: results of the ADVISE (ADenosine Vasodilator Independent Stenosis Evaluation) study. J Am Coll Cardiol. 2012;59(15):1392–402.

- 25.
Kim JS, Lee HD, Suh YK, Kim JH, Chun KJ, Park YH, Kim J, Han DC, Sohn CB. Prediction of fractional flow reserve without hyperemic induction based on resting baseline Pd/Pa. Korean Circ J. 2013;43(5):309–15.

- 26.
Banerjee RK, Ashtekar KD, Helmy T, Effat M, Back LH, Khoury SF. Hemodynamic diagnostics of epicardial coronary stenoses: in vitro experimental and computational study. Biomed Eng Online. 2008;7(1):24.

- 27.
Kolli KK, van de Hoef TP, Effat MA, Banerjee RK, Peelukhana SV, Succop P, Leesar MA, Imran A, Piek JJ, Helmy TA. Diagnostic cutoff for pressure drop coefficient in relation to fractional flow reserve and coronary flow reserve: a patient-level analysis. Catheter Cardiovasc Interv. 2016;87(2):273–82.

- 28.
Kolli KK, Helmy TA, Peelukhana SV, Arif I, Leesar MA, Back LH, Banerjee RK, Effat MA. Functional diagnosis of coronary stenoses using pressure drop coefficient: a pilot study in humans. Catheter Cardiovasc Interv. 2014;83(3):377–85.

- 29.
Kolli KK, Arif I, Peelukhana SV, Succop P, Back LH, Helmy TA, Leesar MA, Effat MA, Banerjee RK. Diagnostic performance of pressure drop coefficient in relation to fractional flow reserve and coronary flow reserve. J Invasive Cardiol. 2014;26(5):188–95.

- 30.
Kolli KK, Effat MA, Peelukhana SV, Succop P, Back LH, Leesar MA, Helmy TA, Imran A, Banerjee RK. Hyperemia-free delineation of epicardial and microvascular impairments using a basal index. Ann Biomed Eng. 2014;42(8):1681–90.

- 31.
Murray CD. The physiological principle of minimum work: I. The vascular system and the cost of blood volume. Proc Natl Acad Sci USA. 1926;12(3):207–14.

- 32.
Jozwik K, Obidowski D. Numerical simulations of the blood flow through vertebral arteries. J Biomech. 2010;43(2):177–85.

- 33.
George AF, Seber CJW. Nonlinear regression. NJ: Wiley-Interscience; 2003.

- 34.
Liu G, Wu J, Huang W, Wu W, Zhang H, Wong K, Ghista D. 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(4):605–11.

- 35.
Wong KKL, Thavornpattanapong P, Cheung SCP, Tu J. Biomechanical investigation of pulsatile flow in a three-dimensional atherosclerotic carotid bifurcation model. J Mech Med Biol. 2013;13(1):1–21.

- 36.
Tu J, Wong KKL, Cheung SCP, Beare R, Phan TG. Analysis of patient-specific carotid bifurcation models using computational fluid dynamics. J M Imaging Health Inform. 2011;1(2):116–25.

- 37.
Wong KKL, Thavornpattanapong P, Tu J, Chaichana T, Sun Z. Analysis on intra-aneurysmal flow influence by stenting. Biomed Eng Inform. 2010;1:113–7.

- 38.
Wong KKL, Kelso RM, Worthley SG, Sanders P, Mazumdar J, Abbott D. Cardiac flow analysis applied to phase contrast magnetic resonance imaging of the heart. Ann Biomed Eng. 2009;37(8):1495–515.

- 39.
Wong KKL, Kelso RM, Worthley S, Sanders P, Mazumdar J, Abbott D. Medical imaging and processing methods for cardiac flow reconstruction. J Mech Med Biol. 2009;9(1):1–20.

- 40.
Guo W, Liu X, Gao Z, Pirbhulal S, Huang W, Lin WH, Zhang H, Tan N, Zhang YT. Quantification of three-dimensional computed tomography angiography for evaluating coronary luminal stenosis using digital subtraction angiography as the standard of reference. Biomed Eng Online. 2015;14:50.

- 41.
Yang Y, Liu X, Xia Y, Liu X, Wu W, Xiong H, Zhang H, Xu L, Wong KKL, Ouyang H, Huang W. Impact of spatial characteristics in the left stenotic coronary artery on the hemodynamics and visualization of 3D replica models. Sci Rep. 2017;7(1):15452.

- 42.
Liu X, Peng C, Xia Y, Gao Z, Xu P, Wang X, Xian Z, Yin Y, Jiao L, Wang D, et al. Hemodynamics analysis of the serial stenotic coronary arteries. Biomed Eng Online. 2017;16(1):127.

- 43.
Peng C, Wang X, Xian Z, Liu X, Huang W, Xu P, Wang J. The impact of the geometric characteristics on the hemodynamics in the stenotic coronary artery. PLoS ONE. 2016;11(6):e0157490.

- 44.
Sutera PS, Skalak R. The history of Poiseuille’s law. Annu Rev Fluid Mech. 2003;25(1):1–19.

- 45.
De Rosa S, Zhang J-M, Zhong L, Luo T, Lomarda AM, Huo Y, Yap J, Lim ST, Tan RS, Wong ASL, et al. Simplified models of non-invasive fractional flow reserve based on CT images. PLoS ONE. 2016;11(5):e0153070.

- 46.
Kim HJ, Vignon-Clementel IE, Figueroa CA, LaDisa JF, Jansen KE, Feinstein JA, Taylor CA. On coupling a lumped parameter heart model and a three-dimensional finite element aorta model. Ann Biomed Eng. 2009;37(11):2153–69.

## Authors’ contributions

XX contributed to the experimental design, analysis of data. MZ and DW contributed to the cCTA data analysis. YL and SX carried out the experiments and the statistical analysis. XX provided final approval of the version of submitted manuscript. All authors read and approved the final manuscript.

### Acknowledgements

The authors gratefully acknowledge the support provided by the National Natural Science Foundation of China (Grant Nos. 61601368, 81471663), the Natural Science Basic Research Plan in Shaanxi Province of China (2017JQ6012) and the Fundamental Research Funds for the Central Universities (No. 3102017jc11002).

### Competing interests

The authors declare that they have no competing interests.

### Ethical approval and consent to participate

Written informed consent was obtained from the patient for the publication of this report and any accompanying images.

### Publisher’s Note

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

## Author information

## Rights and permissions

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

## About this article

### Cite this article

Xie, X., Zheng, M., Wen, D. *et al.* A new CFD based non-invasive method for functional diagnosis of coronary stenosis.
*BioMed Eng OnLine* **17, **36 (2018). https://doi.org/10.1186/s12938-018-0468-6

Received:

Accepted:

Published:

### Keywords

- Functional coronary stenosis
- Non-invasive assessment
- Pressure-flow curve
- CFD