computational fluid dynamics ( cfd ) is an interdisciplinary of computer technology and fluid physics . fluid dynamics parameters within coronary arteries , including blood pressure distribution , wall shear stress ( wss ) distribution and velocity distribution , can be derived by combining coronary artery three - dimensional data , which were obtained by computed tomography ( ct ) scans , with cfd . adding functional features to morphological features offers a new method for non - invasive assessment of myocardial perfusion . the main purpose of this study was to explore technical routes of assessing the left coronary artery pressure gradient , wss distribution and blood flow velocity distribution , combining three - dimensional coronary model which was based on high - resolution dual - source ct with cfd simulation . this was a retrospective study enrolled with three patients who underwent dual - source coronary computed tomography angiography ( ccta ) between march and april in 2014 . case 1 was a 35-year - old male patient with no obvious stenosis in the left anterior descending ( lad ) . case 2 was a 65-year - old male patient with mild stenosis in the lad ( a stenosis of 20% in the proximal - middle segments of lad diagnosed by ccta ) . case 3 was a 60-year - old male patient with severe stenosis in the lad ( a stenosis of 90% in the proximal segments of lad diagnosed by ccta , fig . ccta was performed by using retrospective electrocardiographic gating and dual - source ct ( siemens somatom definition , germany ) . the scan was acquired during injection of 70 ml of lomeprol ( bracco ) at 4.0 ml / s , followed by 40 ml of saline at the same rate . scans were automatically triggered by using bolus tracking when the contrast attenuation reached 100 hu in the ascending aorta . scan parameters were as follows : tube current , 300 - 500 ma depending on body mass index and tube voltage both fixed at 120 kv . oral b - blockers were administered by targeting a heart rate of less than 75 beats per minute , and sublingual nitrates were administered immediately before the scan to optimize coronary vasodilatation . two methods were used to build three - dimensional luminal model of the left coronary arteries . the first method was to input complete original axial images into software mimics 10.01 ( materialise , belgium ) and extract soft tissues to build complete coronaries artery model by identifying gray scales . the second one was to rebuild axial images as follows : field of view ( fov ) 200  200 mm , the left coronary arteries only without chambers or myocardium , and then generate the model in software mimics directly ( fig . three - dimensional luminal model of the left coronary arteries was generated in software mimics . the three - dimensional model was input into software icemcfd 12.0 ( ansys , usa ) for grid generation . five kinds of densities of mesh were undertaken to test the relationship between overall numbers of grid and computation which was described in a graph : overall numbers as abscissa and pressure descending of the lad as ordinate ( fig . the total number of grid ranged from about 240,000 to 4,660,000 , and part of the diagram is shown in figure 3 ( case 2 as the example ) . so choosing the minimum numbers of grids which would not change computation could raise efficiency and guarantee veracity . part of the diagram of five kinds of densities of mesh of case 2 ( a ) . the relationship between overall numbers of grid and computation which was described in a graph ( b ) : overall numbers as abscissa and pressure descending of the lad as ordinate . pressure gradient change was in line with the trend of overall numbers of grid and differences were getting smaller . according to the graph of coronary arteries flow changes during cardiac cycle which was reported by chaichana et al , quantity of flow of the entry of the left coronary artery ( inlet ) was determined ( fig . the exports were set as free flow boundary ( out-1 to 11 ) , and the flow allocation was determined by diameter of each export depending on murray law . the following assumptions were made as such [ 4 - 6 ] : 1 ) the vessel walls were smooth , rigid and non - porous ; 2 ) the blood was considered incompressible , newtonian fluid obeying navier - stokes equations ; 3 ) the blood viscosity was 0.0035 pas , and blood density equaled 1,060 kg / m ; and 4 ) the flow was unsteady laminar . after grid division and setting boundary conditions , the left coronary arteries model was imported into software fluent 12.0 ( ansys , usa ) to simulate flow pressure gradient ( the entry was one ) , wss distribution and blood flow velocity distribution of the lad . the graph of coronary arteries flow changes during cardiac cycle was reported by chaichana et al : time as abscissa ( s ) and quantity of flow as ordinate ( ml / min ) . flow of the entry of the left coronary artery ( inlet ) at different moment could be acquired . the three segments of lad ( proximal , middle and distal segment ) were divided into three parts averagely and severally , and then mean values of each segment were computed according to simulation and ct values , excluding luminal wall , plaques or calcifications . flow pressure gradient and ct values of both ends of stenosis of cases 2 and 3 were recorded . coronary enhancement ratio ( cer ) was ct values above mentioned / ct values of aortic root . ccta was performed by using retrospective electrocardiographic gating and dual - source ct ( siemens somatom definition , germany ) . the scan was acquired during injection of 70 ml of lomeprol ( bracco ) at 4.0 ml / s , followed by 40 ml of saline at the same rate . scans were automatically triggered by using bolus tracking when the contrast attenuation reached 100 hu in the ascending aorta . scan parameters were as follows : tube current , 300 - 500 ma depending on body mass index and tube voltage both fixed at 120 kv . oral b - blockers were administered by targeting a heart rate of less than 75 beats per minute , and sublingual nitrates were administered immediately before the scan to optimize coronary vasodilatation . two methods were used to build three - dimensional luminal model of the left coronary arteries . the first method was to input complete original axial images into software mimics 10.01 ( materialise , belgium ) and extract soft tissues to build complete coronaries artery model by identifying gray scales . the second one was to rebuild axial images as follows : field of view ( fov ) 200  200 mm , the left coronary arteries only without chambers or myocardium , and then generate the model in software mimics directly ( fig . three - dimensional luminal model of the left coronary arteries was generated in software mimics . the three - dimensional model was input into software icemcfd 12.0 ( ansys , usa ) for grid generation . five kinds of densities of mesh were undertaken to test the relationship between overall numbers of grid and computation which was described in a graph : overall numbers as abscissa and pressure descending of the lad as ordinate ( fig . the total number of grid ranged from about 240,000 to 4,660,000 , and part of the diagram is shown in figure 3 ( case 2 as the example ) . so choosing the minimum numbers of grids which would not change computation could raise efficiency and guarantee veracity . part of the diagram of five kinds of densities of mesh of case 2 ( a ) . the relationship between overall numbers of grid and computation which was described in a graph ( b ) : overall numbers as abscissa and pressure descending of the lad as ordinate . pressure gradient change was in line with the trend of overall numbers of grid and differences were getting smaller .   according to the graph of coronary arteries flow changes during cardiac cycle which was reported by chaichana et al , quantity of flow of the entry of the left coronary artery ( inlet ) was determined ( fig . the exports were set as free flow boundary ( out-1 to 11 ) , and the flow allocation was determined by diameter of each export depending on murray law . the following assumptions were made as such [ 4 - 6 ] : 1 ) the vessel walls were smooth , rigid and non - porous ; 2 ) the blood was considered incompressible , newtonian fluid obeying navier - stokes equations ; 3 ) the blood viscosity was 0.0035 pas , and blood density equaled 1,060 kg / m ; and 4 ) the flow was unsteady laminar . after grid division and setting boundary conditions , the left coronary arteries model was imported into software fluent 12.0 ( ansys , usa ) to simulate flow pressure gradient ( the entry was one ) , wss distribution and blood flow velocity distribution of the lad . the graph of coronary arteries flow changes during cardiac cycle was reported by chaichana et al : time as abscissa ( s ) and quantity of flow as ordinate ( ml / min ) . flow of the entry of the left coronary artery ( inlet ) at different moment could be acquired . the three segments of lad ( proximal , middle and distal segment ) were divided into three parts averagely and severally , and then mean values of each segment were computed according to simulation and ct values , excluding luminal wall , plaques or calcifications . flow pressure gradient and ct values of both ends of stenosis of cases 2 and 3 were recorded . coronary enhancement ratio ( cer ) was ct values above mentioned / ct values of aortic root . it took 30 h with first method and 10 h with second method in reconstructing three - dimensional luminal model in software mimics . pressure gradient changes were in line with the trend of overall numbers of grid and differences were getting smaller . in case 2 cfd simulations of the lad included pressure gradient which was the coronary artery pressure relative to the entry ( fig . there was little descending in the pressure gradient of cases 1 and 2 . but compared to the normal left circumflex artery ( lca ) , the descending in the pressure gradient of case 3 between both ends of stenosis of lad was large . wss values in lumen without plaques or stenosis in three cases were low and uniform . at regions opposite to the flow divider , dominant low wss values occurred . in case 3 , wss values of severe stenosis were higher . blood flow velocity distributed uniformly in the lad of cases 1 and 2 , and center axial flow was faster than border flow . but there was obvious variation in lad of case 3 : acceleration of flow velocity on both ends of stenosis , the fastest flow velocity in the center of stenosis and obvious turbulence after end of stenosis . pressure distribution simulated diagram of cases 2 and 3 ( the entry was one ) ( a1-a2 ) . but the descending in the color difference which was pressure gradient of case 3 between proximal and distal of stenosis of lad was large . noteworthy reduction of pressure was observed in the distal segment of left septal branch of case 2 . it showed the wss distribution simulated diagram of cases 2 and 3 ( pa ) ( b1-b2 ) . wss values in lumen without plaques or stenosis were low and uniform . at regions opposite to the flow divider , dominant low wss values occurred . in case 3 , wss values of severe stenosis were higher . it showed the flow velocity distribution simulated diagram of cases 2 and 3 ( m / s ) ( c1-c2 ) . blood flow velocity distributed uniformly in the lad , and center axial flow was faster than border flow . but there was obvious variation in lad of case 3 : acceleration of flow velocity on both ends of stenosis , the fastest flow velocity in the center of stenosis , obvious turbulence after end of stenosis . table 1 lists the mean values of flow pressure gradient and cer in the lad of three cases . differences of pressure gradient and cer values between both ends of severe stenosis of case 3 were larger than mild stenosis of case 2 . the difference of time between two methods in building three - dimensional luminal model was large , indicating that it shortened much time with partial axial images and simplified the process . choosing appropriate numbers of grid which guarantee quality of computation would shorten process too . compared to ffrct which was based on super computers in dicover - flow and defacto study [ 7 , 8 ] , this study applied workstation with relevant softwares and cost lower . the trends of flow pressure gradient and cer values of lad of three cases were roughly same . it detected pressure difference between both ends of stenosis of case 3 in cfd simulation sensitively , which was in line with cer values of corresponding segments . it was reported that cer values were significantly different between both ends of stenosis of more than 50% by nagao et al [ 8 , 9 ] . according to steigner et al and wong et al , ct values gradient in coronary arteries could predict myocardial ischemia related to stenosis [ 10 , 11 ] . thus , cfd simulation might predict myocardial ischemia too and show it to clinicians and patients with pseudo - color images intuitively . noteworthy reduction of pressure was observed in the distal segment of left septal branch of case 2 because of direct oppression from myocardium and high resistance . besides , colors of left main coronary artery and lad were inconformity in the two cases shown in figure 5a1-a2 . it would not affect evaluation of pressure gradient but it did remind researchers to avoid similar bias in follow - up work . high curvature and severe stenosis affected the velocity distribution at the flow divider , giving rise to high wss values , which led to hemorrhage and rupture in atherosclerotic lesions [ 6 , 12 , 13 ] . the assumption made about the nature of the flow was that it was incompressible newtonian fluid according to studies domestic and overseas at present . some researchers underwent studies to compare newtonian fluid models and non - newtonian fluid models which reported similar wss distributions and different velocity distributions . as non - newtonian fluid model , axial velocity at curve regions tended to be more gentle , limiting diffusion of contrast agent [ 4 , 14 - 17 ] . even so , the study offered technical assistance to cfd simulation of non - newtonian fluid model and flexible vessel walls . ccta images based on dual - source ct were acquired during several cardiac cycles . in follow - up studies , 320-detector row ct will be used to get time - uniformity coronary arteries images . but its temporal resolution is low so it requires lower heart rates during examination . with a heart rate of less than 70 beats per minute , otherwise it needs to scan two or more cycles and multi - sectors reconstructions to improve resolution , with more time of exposure and radiation dose . parameters were assumptions according to theories and short of actual measurement including initial conditions , flow distribution ratio and outflow boundary conditions . three - dimensional luminal model of the left coronary arteries were not reconstructed with multi - phases limiting to temporal resolution . however , chow et al   reported that ct values gradient might be more sensitive and precise under stressing state , which will be solved in follow - up studies . parameters were assumptions according to theories and short of actual measurement including initial conditions , flow distribution ratio and outflow boundary conditions . three - dimensional luminal model of the left coronary arteries were not reconstructed with multi - phases limiting to temporal resolution . however , chow et al   reported that ct values gradient might be more sensitive and precise under stressing state , which will be solved in follow - up studies . the study was supported by national natural science foundation of china ( no . 81271570 ) .
backgroundadding functional features to morphological features offers a new method for non - invasive assessment of myocardial perfusion . this study aimed to explore technical routes of assessing the left coronary artery pressure gradient , wall shear stress distribution and blood flow velocity distribution , combining three - dimensional coronary model which was based on high resolution dual - source computed tomography ( ct ) with computational fluid dynamics ( cfd ) simulation.methodsthree cases of no obvious 
 stenosis , mild stenosis and severe stenosis in left anterior descending ( lad ) were enrolled . 
 images acquired on dual - source ct were input into software mimics , icemcfd and fluent to simulate pressure gradient , wall shear stress distribution and blood flow velocity distribution . measuring 
 coronary enhancement ratio of coronary artery was to compare with pressure gradient.resultsresults conformed to theoretical values and showed difference between normal and abnormal samples.conclusionsthe study verified essential parameters and basic techniques in blood flow numerical simulation preliminarily . 
 it was proved feasible .