Skip to main content

Clinical validation and assessment of aortic hemodynamics using computational fluid dynamics simulations from computed tomography angiography



Hemodynamic information including peak systolic pressure (PSP) and peak systolic velocity (PSV) carry an important role in evaluation and diagnosis of congenital heart disease (CHD). Since MDCTA cannot evaluate hemodynamic information directly, the aim of this study is to provide a noninvasive method based on a computational fluid dynamics (CFD) model, derived from multi-detector computed tomography angiography (MDCTA) raw data, to analyze the aortic hemodynamics in infants with CHD, and validate these results against echocardiography and cardiac catheter measurements.


This study included 25 patients (17 males, and 8 females; a median age of 2 years, range: 4 months–4 years) with CHD. All patients underwent both transthoracic echocardiography (TTE) and MDCTA within 2 weeks prior to cardiac catheterization. CFD models were created from MDCTA raw data. Boundary conditions were confirmed by lumped parameter model and transthoracic echocardiography (TTE). Peak systolic velocity derived from CFD models (PSVCFD) was compared to TTE measurements (PSVTTE), while the peak systolic pressure derived from CFD (PSPCFD) was compared to catheterization (PSPCC). Regions with low and high peak systolic wall shear stress (PSWSS) were also evaluated.


PSVCFD and PSPCFD showed good agreements between PSVTTE (r = 0.968, p < 0.001; mean bias = − 7.68 cm/s) and PSPCC (r = 0.918, p < 0.001; mean bias = 1.405 mmHg). Regions with low and high PSWSS) can also be visualized. Skewing of velocity or helical blood flow was also observed at aortic arch in patients.


Our result demonstrated that CFD scheme based on MDCTA raw data is an accurate and convenient method in obtaining the velocity and pressure from aorta and displaying the distribution of PSWSS and flow pattern of aorta. The preliminary results from our study demonstrate the capability in combining clinical imaging data and novel CFD tools in infants with CHD and provide a noninvasive approach for diagnose of CHD such as coarctation of aorta in future.


Congenital heart disease (CHD) is a common malformation affecting approximately six per 1000 live births, occurring as an isolated trait or related to multiple congenital anomalies [1]. Despite anatomical evaluation has a great contribution to the diagnosis and treatments of CHD, the hemodynamic evaluation is also indispensable. Peak systolic pressure (PSP) and peak systolic velocity (PSV) were widely used in diagnosis or grading of CHD such as pulmonary arterial hypertension (PAH), aortic valve stenosis (AS), and coarctation of aorta (CoA) [2,3,4]. An accurate and noninvasive approach to evaluate such hemodynamic information may carry an important role to benefit the clinical diagnosis or grading of CHD.

Cardiac catheterization is the reference standard in obtaining patient-specific hemodynamic analysis, but is associated with patient discomfort and carries potential peri-procedural risks [5]. Transthoracic echocardiography (TTE) is the effective first-line technology for obtaining PSV, but is limited in the ability to provide accurate PSP information. Cardiac magnetic resonance (CMR) is another noninvasive method to evaluate hemodynamics in patients [6, 7], but it is expensive and more technically demanding. Currently, multi-detector computed tomography angiography (MDCTA) has been widely accepted as an accurate imaging modality to evaluate cardiovascular anatomy for its conveniences, noninvasive procedure and low cost; however, MDCTA cannot provide hemodynamic information directly. Recent studies showed that reliable hemodynamic results of carotid artery, coronary artery and other artery can be acquired using MDCTA raw data alone [8,9,10,11,12], but only few studies had focused on aortic hemodynamic schemes based MDCTA [13].

The aim of this study is to provide a noninvasive method based on a computational fluid dynamics (CFD) model, derived from multi-detector computed tomography angiography (MDCTA) raw data, to analyze the aortic hemodynamics in infants with CHD, and validate these results against echocardiography and cardiac catheter measurements.


Study population

We conducted a single center search of all patients with CHD from July 2015 to October 2016 in our department. Patients were included if they: (1) infants under 4 years old. (2) Had MDCTA and TTE studies, with flow measurement by TTE. (3) Cardiac catheterization within 2 weeks from MDCTA and TTE, with PSP measurement. (4) No malformation and lesion of aorta; Patients were excluded from this study if the studies had poor image quality, or if the above inclusion criteria were not met. There are twenty-five patients (17 males, 8 females) with known CHD were enrolled finally, with a median age of 2 years (range, 4 months–4 years).

Transthoracic echocardiography, MDCTA and cardiac catheterization protocol

Before undergoing MDCTA, TTE (Philips iE 33 imaging system, Philips Healthcare, Andover, Massachusetts, USA) was performed.

A second-generation dual source CT scanner (Somatom Definition Flash, Siemens Health-care, Forchheim, Germany) was performed using electrocardiographic-gated step and shoot protocol. Short-term sedation was achieved with 0.1 mg/ml of oral chloral hydrate solution. The scans were performed in cranio-caudal direction from the thoracic inlet to the bottom of the heart. MDCTA parameters were as follows: 0.28 s gantry rotation time, 2 × 64 × 0.6 mm detector collimation, CARE kV (weight adapted setting for tube voltage and tube current). In the sequential mode, the acquisition window was set at 35–45% of the R–R interval. Safire (strength 3) is selected as the iterative reconstruction algorithm, the kernel is I26, the slice thickness is 0.75 mm, and increment is 0.5 mm. Iodinated contrast medium (Iopamidol, 300 mg I/ml, BRACCO, Italy) was injected intravenously at a volume of 1.5–2.0 ml/kg body weight, followed by 1.0 ml/kg body weight saline chaser with injection rate 1–2 ml/s. The acquisition delay was determined by the time of contrast medium entering both ventricles.

CC was performed by using Philips Allura Xper FD10 system (Philips Medical Systems, Best, the Netherlands) to get hemodynamics including PSP of the aortic isthmus (AI).

Construction of aorta and mesh generation

Computational representations of the aorta were created using Mimics 17 (Materialise, Leuven, Belgium) software that facilitates volume visualization and conversion of the MDCTA raw data into geometrically representative computer models, as demonstrated in Fig. 1. Models were discretized using an automatic mesh generation software (Ansys ICEM 14.5, ANSYS, Inc., Canonsburg, Pennsylvania, USA). The mesh generation followed a custom standard protocol with unstructured, formatted tetrahedral. In addition, prism mesh was added to the wall boundary (height = 1, ratio = 1.2, numbers = 3) in order to improve the accuracy of hemodynamic parameters adjacent to aortic wall. The total elements in each mesh ranging from 1,000,000 to 2,000,000.

Fig. 1

The patient-specific aortic geometry. An inlet boundary and four outlet boundaries were defined in left image. O, R, L and I represented the four partition of the aortic wall demonstrated in right image. The region identified by ‘AI’ was aortic isthmus

Boundary conditions

Inflow and outflow boundaries were defined in Fig. 1, and the flow domain was defined as cavity of the reconstructed geometry. Velocity information in TTE data was mapped to the inlet of CFD models, while the blood pressure (BP) was used to prescribe the outlet boundaries. Lumped parameter model (LPM) was applied to confirm the outflow boundary conditions if BP of outlet was inaccessible. The parameter of LPM was shown in Table 1, and the schematic illustration of LPM was shown in Fig. 2. The unknown pressure (P) of outlet for CFD models was calculated from LPM, the equation of LPM was given as follow:

Table 1 Parameter of LPMs
Fig. 2

Scheme of LPM. R1 represented the characteristic resistance, R2 represented the peripheral impedance, and C represented the compliance of artery. Q was flow rate of artery, and pressure on outlet was calculated by solving the differential equation

$$P = (R_{1} + R_{2} )\,Q - R_{2} C\frac{dP}{dt} + R_{1} R_{2} \frac{dQ}{dt}$$

where Q was the volume flow rate through brachiocephalic artery (BA), left common carotid artery (LCCA), left subclavian artery (LSA), or descending aorta (DAo) respectively. R1 was the characteristic resistance of the artery, while R2 was Peripheral impedance of the artery, and C was the compliance of artery. The value of R1, R2, C were calculated according to previous work [13, 14]. The derivative items can be calculated using the backward Euler method:

$$P_{n + 1} = \frac{{(R_{1} + R_{2} + R_{1} \beta ) - R_{1} \beta Q_{n} + \beta P_{n} }}{1 + \beta }$$
$$\beta = R_{2} C/\Delta t$$

where Δt was time interval between Pn and Pn+1.

Computational fluid dynamics simulation

Simulation was performed by using Ansys FLUENT 14.5 (Fluent Inc., Lebanon, New Hamshire, USA). In our calculation, there were some assumptions of physical properties of blood according to previous work: (1) Blood was assumed to be incompressible, viscous, Newtonian fluid. (2) The density (ρ) and viscosity (μ) of blood was assumed to be constant and equals to 1050 kg/m3 and 0.004 Pa s [14,15,16,17]. Calculation solved the Navier–Stokes equations as follow [18]:

$$\rho \left( {\frac{du}{dt} + u \bullet \nabla u} \right) = - \nabla P + \mu \nabla^{2} + F$$
$$- \nabla \bullet u = 0$$

where u was the blood velocity, F was body force equals to zero. Fluid structure interaction simulation was wildly used to illustrate the interaction between blood and vessels [19,20,21]. Since the hemodynamics of the blood is the main research topic in this current study, we have not included fluid structure interaction simulation into this study. Similar approaches had been applied in other related studies to acquire accurate hemodynamics from simulation without FSI [14, 22,23,24].

Statistical analysis

The statistical analysis was performed by SPSS (SPSS 22, SPSS Inc., Chicago, USA). For continuous variables, data are expressed as the mean ± standard deviation (M ± SD). All tests were two-sided, and effects were considered significant at p < 0.05. Normality was tested using the Kolmogorov–Smirnov method, and variance homogeneity was then tested with the Levene’s test. Group differences were assessed by paired Student’s t test in normally distributed (Kolmogorov–Smirnov test) data. Otherwise, the paired Wilcoxon test was used.

To illustrate the accuracy of aorta reconstruction, morphometric parameters (diameters at three locations: aorta ascending, aorta descending, and stenosis) were analyzed by paired Student’s t-test. To demonstrate the accuracy of our simulation, the agreements between PSVCFD and PSVTTE, PSPCFD and PSPCC were characterized with the Bland–Altman plot and linear fitting.


Table 2 summarizes the reconstructed and measured geometric parameters of the aorta. Mean reconstructed diameters of the ascending aorta (inlet), descending aorta (outlet.4) and AI were 19. 2 ± 6.0 mm, 10.9 ± 2.6 mm, and 12.1 ± 3.5 mm respectively. The paired t test found no difference between measured diameters and simulated diameters.

Table 2 Calculated and measured geometric parameters

The data on PSPCC, PSVTTE, PSPCFD, and PSVCFD are given in Table 3. PSPCFD had an excellent correlation (Fig. 3a) with PSPCC (r = 0.918, p < 0.001). The mean PSPCC was 105.08 ± 15.38 mmHg, while mean PSPCFD was 106.48 ± 15 mmHg. The mean bias was 1.405 mmHg (Fig. 3b, 95% confidence interval − 7.237–10.04). PSVCFD was also excellently correlated with PSVTTE (Fig. 4a, r = 0.968, p < 0.001). The mean PSVTTE was 152.92 ± 64.36 cm/s, while mean PSVCFD was 145.24 ± 61.68 cm/s. The mean bias was − 7.68 cm/s (Fig. 4b, 95% confidence interval − 30.41 to 15.05).

Table 3 Measured and simulated PSV and PSP
Fig. 3

Validation of PSP derived from CFD. Demonstration of PSV and PSP validation. The image a was linear fitting of PSP, and image b was Bland–Altman plot of PSP. The reference line of Bland–Altman plots was mean difference ± 1.96 *SD

Fig. 4

Validation of PSV derived from CFD. Demonstration of PSV and PSP validation. The image a was linear fitting of PSV, and image b was Bland–Altman plot of PSV. The reference line of Bland–Altman plots was mean difference ± 1.96 *SD

During peak systole, on the one hand, lowest PSWSS (2.83 ± 1.23 Pa) was observed at right wall of ascending aorta (AscAo) and inner wall of the attachment of arterial ductus ligament adjacent to AI. On the other hand, highest PSWSS (16.26 ± 3.43 Pa) was localized at outer wall of aortic arch, ostium of BA (Fig. 5). Streamline was also visualized and demonstrated in Fig. 5, velocity skewed toward the inner wall of AscAo and the region between inner and right wall of DAo (Fig. 5). Helical flow was observed at arch (Fig. 5), and there was highest velocity observed at ostium of BA (Fig. 5).

Fig. 5

Distribution of PSWSS and streamline. Distribution of PSWSS was demonstrated in image a. Highest PSWSS was marked by red arrow, and lowest PSWSS was marked by blue one. Streamline at peak systolic was shown in image b, and the helical flow was marked by red arrow, and the projection of helical flow at cross section was demonstrated in image c


This study compared PSVCFD with PSVTTE and PSPCFD with PSPCC respectively. We demonstrated good agreements between reference standard method and CFD method.

Concerning the findings, there are additional factors that must be considered. Obviously, diameters of geometry reconstructed from MDCTA raw data will affect the CFD results, however, paired t-test found no difference between measured and reconstructed diameters, it means that geometry constructed from MDCTA raw data was identical to the real anatomy.

Given the results in the previous section, we observed that there are good agreements between PSVTTE versus PSVCFD and PSPCC versus PSPCFD. The result demonstrated that CFD method was accurate comparing to TTE and cardiac catheterization. The biases of PSV and PSP were small, and there are many factors contributed to the biases. Data acquisition with MDCTA and catheterization was performed sequentially and thus there could be slight day-to-day variations in stroke volumes and other hemodynamics parameters [6]. More importantly, PSP was measured in sedated patients, whereas MDCTA and TTE were measured when patients awake. But according to the results, the biases were acceptable.

During MDCTA processing, there is a level of uncertainty when reconstructing the 3D aorta geometry based on 2D MDCTA raw data. Precise reconstruction is critical for the outcome of CFD simulation. However, no significant differences between measured and reconstructed geometric parameters were observed in this study.

The step of simulation may also be a factor contributing to the biases. Lumped parameter model was widely used to implement the pressure at outlet boundary according to previous studies [9, 14, 18, 25,26,27]. However, modeling hemodynamics in aorta with LPM involves a challenging set of constraints. In addition, previous studies have found that no-slip wall boundary is suitable and it is wildly used in vessel modeling. However, accurate wall boundary definition is still an active topic of investigation [28, 29]. In this work, second order no-slip wall boundary was included in the simulation, with the rational to reduce the computational effort and complexity of simulation.

In most of cases, low PSWSS of aorta was localized at outer or right luminal surface of ascending aorta or inner wall of the attachment of arterial ductus ligament adjacent to AI, while high PSWSS of aorta was located at outer wall of aortic arch, ostium of BA. This result corresponds to the distribution of WSS reported in other studies [14, 24]. In the current study, low PSWSS was observed at outer or right wall of AscAo maybe explained by the skewing velocity profile towards the inner wall during peak systole, while the velocity imposed to outer wall of AscAo was low. On the contrary, the highest PSWSS observed on outer wall of aortic arch and ostium of BA. This maybe explained by the arc shape of aortic arch. In this study, PSWSS on inner wall was lower than PSWSS on AI, and the larger diameter of arterial ductus ligament attachment as compared to AI maybe a factor contributed to this phenomenon.

Study conducted by Chiu et al. proved that PSWSS played an important role for regulating the arrangement and function of endotheliocyte [30]. Previous study also suggests that region with low PSWSS was correlated with areas of atherosclerotic plaque [31], and other study also indicated that excessive PSWSS was related with aneurysm formation of endothelial cells [32]. In addition, PSWSS also played an important role in the formation of aortic dissection (AD). Thubrikar reported that elevated PSWSS was correlated with sites of intimal tears [33], Wen also reported that initial location of tears was coincident with the region of maximal WSS [34], and Nordon reported that low PSWSS was benefit for minimizing the propagation of the dissection [35]. These studies indicated that PSWSS was related to the formation and prognosis of AD. Hence, WSS derived from the MDCTA raw data may provide not only anatomic information on aorta, but also additional information on plaque and dissection development.

Uniform distribution and smooth laminar blood flow was observed at AscAo, while skewing and helical flow was observed at aorta during peak systole. The flow pattern in current study was consistent with other studies [36, 37]. The inertial force was higher than viscous force during peak systolic, which help the development of helical flow, a usual physiological phenomenon in our cardiovascular system [38]. Helical flow in aorta exhibits important hemodynamic effect in increasing flux of oxygen [39] and achieving normal level of PSWSS [40]. Moreover, helical blood flow hinders the development of atherosclerotic plaques [23]. Similarly, the luminal surface low-density lipoprotein concentration in the aortic arch can also be reduced by helical flow [41]. Flow pattern generated by CFD based on MDCTA raw data may provide a reliable and convenient method to evaluate flow information.

Study limitations

There were some limitations in the current study. First, the aortic valve morphology was not considered for its influence on the velocity profile. However, many studies have proven that the influence of valve morphology on velocity was small, and accurate simulation results had been achieved without considering aortic valvular morphology [14, 24, 42, 43]. Second, Constant LPM parameters were used to implement CFD in order to simplify the process of calculation despite unequal aortic morphology in each patient. However, their validated results suggested that the simulation error was negligible compared to other study [44, 45]. Third, this study had a limited number of subjects recruited in our study, and a follow up study with a larger number of subjects may provide a better valid assessment of aortic hemodynamic from MDCTA.


Our results showed that CFD scheme based on MDCTA raw data is an accurate and convenient method in obtaining flow velocity and pressure from aorta and displaying the distribution of WSS and flow pattern of aorta. The preliminary results from our study demonstrated the capability in combining clinical imaging data and novel CFD tools in infants with CHD and provide t provide another noninvasive approach for diagnose of CHD such as CoA, AS or PAH in future.



computational fluid dynamics


transthoracic echocardiography


peak systolic pressure


peak systolic velocity


PSP measured by CFD


PSV measured by CFD


PSP measured by CC


PSV measured by TTE


peak systolic pressure gradient


blood pressure


wall shear stress


peak systolic wall shear stress


lumped parameter model


congenital heart disease


multi-detector computed tomography angiography


cardiac magnetic resonance


pulmonary arterial hypertension


aortic valve stenosis


coarctation of aorta


ascending aorta


brachiocephalic artery


left common carotid artery


left subclavian artery


descending aorta


  1. 1.

    Molck MC, Simioni M, Vieira TP, Sgardioli IC, Monteiro FP, Souza J, Fett-Conte AC, Felix TM, Monlleo IL, Gil-da-Silva-Lopes VL. Genomic imbalances in syndromic congenital heart disease. Jornal de pediatria. 2017;93(5):497–507.

    Article  Google Scholar 

  2. 2.

    Kitabatake A, Inoue M, Asao M, Masuyama T, Tanouchi J, Morita T, Mishima M, Uematsu M, Shimazu T, Hori M, Abe H. Noninvasive evaluation of pulmonary hypertension by a pulsed Doppler technique. Circulation. 1983;68(2):302–9.

    Article  Google Scholar 

  3. 3.

    Suradi H, Hijazi ZM. Current management of coarctation of the aorta. Glob Cardiol Sci Pract. 2015;2015(4):44.

    Article  Google Scholar 

  4. 4.

    Otto CM, Burwash IG, Legget ME, Munt BI, Fujioka M, Healy NL, Kraft CD, Miyake-Hull CY, Schwaegler RG. Prospective study of asymptomatic valvular aortic stenosis. Clinical, echocardiographic, and exercise predictors of outcome. Circulation. 1997;95(9):2262–70.

    Article  Google Scholar 

  5. 5.

    Di Biase L, Burkhardt JD, Mohanty P, Sanchez J, Horton R, Gallinghouse GJ, Lakkireddy D, Verma A, Khaykin Y, Hongo R, et al. Periprocedural stroke and management of major bleeding complications in patients undergoing catheter ablation of atrial fibrillation: the impact of periprocedural therapeutic international normalized ratio. Circulation. 2010;121(23):2550–6.

    Article  Google Scholar 

  6. 6.

    Goubergrits L, Riesenkampff E, Yevtushenko P, Schaller J, Kertzscher U, Hennemuth A, Berger F, Schubert S, Kuehne T. MRI-based computational fluid dynamics for diagnosis and treatment prediction: clinical validation study in patients with coarctation of aorta. J Magn Reson Imaging JMRI. 2015;41(4):909–16.

    Article  Google Scholar 

  7. 7.

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

    Article  Google Scholar 

  8. 8.

    Koo BK, Erglis A, Doh JH, Daniels DV, Jegere S, Kim HS, Dunning A, DeFrance T, Lansky A, Leipsic J, Min JK. Diagnosis of ischemia-causing coronary stenoses by noninvasive fractional flow reserve computed from coronary computed tomographic angiograms. Results from the prospective multicenter DISCOVER-FLOW (diagnosis of ischemia-causing stenoses obtained via noninvasive fractional flow reserve) study. J Am Coll Cardiol. 2011;58(19):1989–97.

    Article  Google Scholar 

  9. 9.

    Liu X, Gao Z, Xiong H, Ghista D, Ren L, Zhang H, Wu W, Huang W, Hau WK. Three-dimensional hemodynamics analysis of the circle of Willis in the patient-specific nonintegral arterial structures. Biomech Model Mechanobiol. 2016;15(6):1439–56.

    Article  Google Scholar 

  10. 10.

    Qin Y, Wu J, Hu Q, Ghista DN, Wong KK. Computational evaluation of smoothed particle hydrodynamics for implementing blood flow modelling through CT reconstructed arteries. J X ray Sci Technol. 2017;25(2):213–32.

    Article  Google Scholar 

  11. 11.

    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(1):1350001.

    Article  Google Scholar 

  12. 12.

    Wong KK, 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.

    Article  Google Scholar 

  13. 13.

    Nie P, Wang X, Cheng Z, Ji X, Duan Y, Chen J. Accuracy, image quality and radiation dose comparison of high-pitch spiral and sequential acquisition on 128-slice dual-source CT angiography in children with congenital heart disease. Eur Radiol. 2012;22(10):2057–66.

    Article  Google Scholar 

  14. 14.

    Menon A, Wendell DC, Wang H, Eddinger TJ, Toth JM, Dholakia RJ, Larsen PM, Jensen ES, LaDisa JF. A coupled experimental and computational approach to quantify deleterious hemodynamics, vascular alterations, and mechanisms of long-term morbidity in response to aortic coarctation. J Pharmacol Toxicol Methods. 2012;65(1):18–28.

    Article  Google Scholar 

  15. 15.

    Sankaranarayanan M, Chua LP, Ghista DN, Tan YS. Computational model of blood flow in the aorto-coronary bypass graft. Biomed Eng Online. 2005;4:14.

    Article  Google Scholar 

  16. 16.

    Wendell DC, Samyn MM, Cava JR, Ellwein LM, Krolikowski MM, Gandy KL, Pelech AN, Shadden SC, LaDisa JF Jr. Including aortic valve morphology in computational fluid dynamics simulations: initial findings and application to aortic coarctation. Med Eng Phys. 2013;35(6):723–35.

    Article  Google Scholar 

  17. 17.

    Albert S, Balaban RS, Neufeld EB, Rossmann JS. Influence of the renal artery ostium flow diverter on hemodynamics and atherogenesis. J Biomech. 2014;47(7):1594–602.

    Article  Google Scholar 

  18. 18.

    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.

    Article  Google Scholar 

  19. 19.

    Wong KK, Thavornpattanapong P, Cheung SC, Tu J. Numerical stability of partitioned approach in fluid-structure interaction for a deformable thin-walled vessel. Comput Math Methods Med. 2013;2013:638519.

    MathSciNet  Article  MATH  Google Scholar 

  20. 20.

    Siguenza J, Pott D, Mendez S, Sonntag SJ, Kaufmann TAS, Steinseifer U, Nicoud F. Fluid-structure interaction of a pulsatile flow with an aortic valve model: a combined experimental and numerical study. Int J Numer Method Biomed Eng. 2017.

    Google Scholar 

  21. 21.

    Sharzehee M, Khalafvand SS, Han HC. Fluid-structure interaction modeling of aneurysmal arteries under steady-state and pulsatile blood flow: a stability analysis. Comput Methods Biomech Biomed Eng. 2018;21:219.

    Article  Google Scholar 

  22. 22.

    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.

    Article  Google Scholar 

  23. 23.

    Liu X, Pu F, Fan Y, Deng X, Li D, Li S. A numerical study on the flow of blood and the transport of LDL in the human aorta: the physiological significance of the helical flow in the aortic arch. Am J Physiol Heart Circ Physiol. 2009;297(1):H163–70.

    Article  Google Scholar 

  24. 24.

    LaDisa JF, Alberto Figueroa C, Vignon-Clementel IE, Jin Kim H, Xiao N, Ellwein LM, Chan FP, Feinstein JA, Taylor CA. Computational simulations for aortic coarctation: representative results from a sampling of patients. J Biomech Eng. 2011;133(9):091008.

    Article  Google Scholar 

  25. 25.

    Liu X, Zhang H, Ren L, Xiong H, Gao Z, Xu P, Huang W, Wu W. Functional assessment of the stenotic carotid artery by CFD-based pressure gradient evaluation. Am J Physiol Heart Circ Physiol. 2016;311(3):H645–53.

    Article  Google Scholar 

  26. 26.

    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.

    Article  Google Scholar 

  27. 27.

    Xu P, Liu X, Song Q, Chen G, Wang D, Zhang H, Yan L, Liu D, Huang W. Patient-specific structural effects on hemodynamics in the ischemic lower limb artery. Sci Rep. 2016;6:39225.

    Article  Google Scholar 

  28. 28.

    Nash RW, Carver HB, Bernabeu MO, Hetherington J, Groen D, Kruger T, Coveney PV. Choice of boundary condition for lattice-Boltzmann simulation of moderate-Reynolds-number flow in complex domains. Phys Rev E Stat Nonlin Soft Matter Phys. 2014;89(2):023303.

    Article  Google Scholar 

  29. 29.

    Boyd J, Buick JM, Cosgrove JA, Stansell P. Application of the lattice Boltzmann method to arterial flow simulation: investigation of boundary conditions for complex arterial geometries. Australas Phys Eng Sci Med. 2004;27(4):207–12.

    Article  Google Scholar 

  30. 30.

    Chiu JJ, Chien S. Effects of disturbed flow on vascular endothelium: pathophysiological basis and clinical perspectives. Physiol Rev. 2011;91(1):327–87.

    MathSciNet  Article  Google Scholar 

  31. 31.

    Wentzel JJ, Corti R, Fayad ZA, Wisdom P, Macaluso F, Winkelman MO, Fuster V, Badimon JJ. Does shear stress modulate both plaque progression and regression in the thoracic aorta? Human study using serial magnetic resonance imaging. J Am Coll Cardiol. 2005;45(6):846–54.

    Article  Google Scholar 

  32. 32.

    Dolan JM, Kolega J, Meng H. High wall shear stress and spatial gradients in vascular pathology: a review. Ann Biomed Eng. 2013;41(7):1411–27.

    Article  Google Scholar 

  33. 33.

    Thubrikar MJ, Agali P, Robicsek F. Wall stress as a possible mechanism for the development of transverse intimal tears in aortic dissections. J Med Eng Technol. 1999;23(4):127–34.

    Article  Google Scholar 

  34. 34.

    Wen CY, Yang AS, Tseng LY, Chai JW. Investigation of pulsatile flowfield in healthy thoracic aorta models. Ann Biomed Eng. 2010;38(2):391–402.

    Article  Google Scholar 

  35. 35.

    Nordon IM, Hinchliffe RJ, Loftus IM, Morgan RA, Thompson MM. Management of acute aortic syndrome and chronic aortic dissection. Cardiovasc Intervent Radiol. 2011;34(5):890–902.

    Article  Google Scholar 

  36. 36.

    LaDisa JF Jr, Dholakia RJ, Figueroa CA, Vignon-Clementel IE, Chan FP, Samyn MM, Cava JR, Taylor CA, Feinstein JA. Computational simulations demonstrate altered wall shear stress in aortic coarctation patients treated by resection with end-to-end anastomosis. Congenit Heart Dis. 2011;6(5):432–43.

    Article  Google Scholar 

  37. 37.

    Zhang Q, Gao B, Chang Y. The study on hemodynamic effect of series type LVAD on aortic blood flow pattern: a primary numerical study. Biomed Eng Online. 2016;15(Suppl 2):163.

    Article  Google Scholar 

  38. 38.

    Wong KK, Wang D, Ko JK, 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(1):35.

    Article  Google Scholar 

  39. 39.

    Liu X, Fan Y, Deng X. Effect of spiral flow on the transport of oxygen in the aorta: a numerical study. Ann Biomed Eng. 2010;38(3):917–26.

    Article  Google Scholar 

  40. 40.

    Caro CG, Cheshire NJ, Watkins N. Preliminary comparative study of small amplitude helical and conventional ePTFE arteriovenous shunts in pigs. J R Soc Interface. 2005;2(3):261–6.

    Article  Google Scholar 

  41. 41.

    Wang Z, Liu X, Kang H, Sun A, Fan Y, Deng X. Enhanced accumulation of LDLs within the venous graft wall induced by elevated filtration rate may account for its accelerated atherogenesis. Atherosclerosis. 2014;236(1):198–206.

    Article  Google Scholar 

  42. 42.

    Mirzaee H, Henn T, Krause MJ, Goubergrits L, Schumann C, Neugebauer M, Kuehne T, Preusser T, Hennemuth A. MRI-based computational hemodynamics in patients with aortic coarctation using the lattice Boltzmann methods: clinical validation study. J Magn Reson Imaging. 2017;45(1):139–46.

    Article  Google Scholar 

  43. 43.

    Kwon S, Feinstein JA, Dholakia RJ, LaDisa JF. Quantification of local hemodynamic alterations caused by virtual implantation of three commercially available stents for the treatment of aortic coarctation. Pediatr Cardiol. 2013;35(4):732–40.

    Article  Google Scholar 

  44. 44.

    Dwyer HA, Matthews PB, Azadani A, Jaussaud N, Ge L, Guy TS, Tseng EE. Computational fluid dynamics simulation of transcatheter aortic valve degeneration. Interact Cardiovasc Thorac Surg. 2009;9(2):301–8.

    Article  Google Scholar 

  45. 45.

    Kilner PJ, Yang GZ, Mohiaddin RH, Firmin DN, Longmore DB. Helical and retrograde secondary flow patterns in the aortic arch studied by three-directional magnetic resonance velocity mapping. Circulation. 1993;88(5 Pt 1):2235–47.

    Article  Google Scholar 

Download references

Authors’ contributions

YZ contributed in the research design and drafting of the manuscript. HL and ZY supervised and coordinated the study and also contributed in designing of research. YZ and RC acquired patient clinical information. YJ contributed the editing manuscript and tables. HL contributed in statistical analysis. YZ, RC, and JW contributed to the image acquisition, image reformation, interpretation, preparation of the figures and figure legends. All authors read and approved the final manuscript.



Competing interests

The authors declared that they have no competing interests.

Availability of data and materials

The datasets used or analyzed during the current study are available from the corresponding author on reasonable request.

Consent for publication

All participates had been informed before submission, and all authors agreed the publication.

Ethic approval and consent to participate

The study was approved by the institutional Research Ethics Committee following the ethical guidelines of the 1975 Declaration of Helsinki. Informed consent was obtained from all individual participants included in the study.


The current study was funded by national natural science foundation of China (No. 81771799).

Publisher’s Note

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

Author information



Corresponding authors

Correspondence to Zhuliang Yu or Hui Liu.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, 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 ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Zhu, Y., Chen, R., Juan, Y. et al. Clinical validation and assessment of aortic hemodynamics using computational fluid dynamics simulations from computed tomography angiography. BioMed Eng OnLine 17, 53 (2018).

Download citation


  • Multi-detector computed tomography angiography
  • Computational fluid dynamics
  • Aortic hemodynamics
  • Congenital heart disease