Nonquasi-steady character of pulsatile flow in human coronary arteries

This experiment was conducted to determine if the pulsatile flow through the proximal portion of the left coronary artery system in man exhibits quasi-steady characteristics. Steady and pulsatile flows were passed through an idealized model whose dimensions were based on a vascular cast. The mean Reynolds number was 180 and the unsteadiness number was 2.7. Velocity profiles were measured by laser Doppler anemometry at several locations along diameters in the parent and both daughter channels in the neighborhood of the "left main" bifurcation. Analysis of the results along one diameter in the "left main" channel shows that unsteady flow in the larger coronary arteries may not be simulated by a series of steady flow experiments.     

Hemodynamic impacts of various types of stenosis in the left coronary artery bifurcation: A patient-specific analysis

This study investigates the hemodynamic changes to various types of coronary stenosis in the left coronary artery bifurcation, based on a patient-specific analysis. Twenty two patients with left coronary artery disease were included in this study. All stenoses involving the left coronary artery bifurcation were classified into four types, according to their locations: A) left circumflex (LCx) and left anterior descending (LAD), B) LCx only, C) left main stem only, and D) LAD only. Computational fluid dynamics (CFD) was performed to analyze the flow and wall shear stress (WSS) changes in all reconstructed left coronary geometries. Our results showed that the flow velocity and WSS were significantly increased at stenotic locations. High WSS was found at >70% lumen stenosis, which ranged from 2.5 Pa to 3.5 Pa. This study demonstrates that in patients with more than 50% stenosis in the left coronary artery bifurcation, WSS plays an important role in providing information about the extent of coronary atherosclerosis in the left coronary artery branch.     

A numerical calculation of flow in a curved tube model of the left main coronary artery

The flow pattern in the left main coronary artery has been calculated using an idealized geometry and by numerically solving the full Navier-Stokes equations for a Newtonian fluid. Two different forms for the entrance velocity profile were used, one a time-varying, flat profile and the other a time-varying, less flat velocity profile. The results obtained demonstrate the presence of secondary motions for conditions simulating flow in the left main coronary artery, with maximum secondary flow velocities being on the order of three to four percent of the maximum axial velocity. This secondary flow phenomenon has an important influence on the wall shear stress distribution, in spite of the fact that there is virtually no alteration in the axial velocity profile. The maximum ratio of the outer wall shear stress to that on the inner wall is 1.4 at a Reynolds number of Re = 270, and it increases with increasing Reynolds number, reaching a value of 1.7 at Re = 810. Although there are significant differences in the results in the immediate vicinity of the inlet for the two different forms of the entrance velocity profile used, this difference does not persist far into the tube. Independent of the choice of the entrance velocity profile, it appears that there will be significant secondary flow effects on the wall shear stress.     

Effect of common carotid artery inlet length on normal carotid bifurcation hemodynamics

Controversy exists regarding the suitability of fully developed versus measured inlet velocity profiles for image-based computational fluid dynamics (CFD) studies of carotid bifurcation hemodynamics. Here, we attempt to resolve this by investigating the impact of the reconstructed common carotid artery (CCA) inlet length on computed metrics of "disturbed" flow. Twelve normal carotid bifurcation geometries were reconstructed from contrast-enhanced angiograms acquired as part of the Vascular Aging--The Link That Bridges Age to Atherosclerosis study (VALIDATE). The right carotid artery lumen geometry was reconstructed from its brachiocephalic origin to well above the bifurcation, and the CCA was truncated objectively at locations one, three, five, and seven diameters proximal to where it flares into the bifurcation. Relative to the simulations carried out using the full CCA, models truncated at one CCA diameter strongly overestimated the amount of disturbed flow. Substantial improvement was offered by using three CCA diameters, with only minor further improvement using five CCA diameters. With seven CCA diameters, the amounts of disturbed flow agreed unambiguously with those predicted by the corresponding full-length models. Based on these findings, we recommend that image-based CFD models of the carotid bifurcation should incorporate at least three diameters of CCA length if fully developed velocity profiles are to be imposed at the inlet. The need for imposing measured inlet velocity profiles would seem to be relevant only for those cases where the CCA is severely truncated.     

Effect of inlet velocity profiles on patient-specific computational fluid dynamics simulations of the carotid bifurcation

Patient-specific computational fluid dynamics (CFD) is a powerful tool for researching the role of blood flow in disease processes. Modern clinical imaging technology such as MRI and CT can provide high resolution information about vessel geometry, but in many situations, patient-specific inlet velocity information is not available. In these situations, a simplified velocity profile must be selected. We studied how idealized inlet velocity profiles (blunt, parabolic, and Womersley flow) affect patient-specific CFD results when compared to simulations employing a "reference standard" of the patient's own measured velocity profile in the carotid bifurcation. To place the magnitude of these effects in context, we also investigated the effect of geometry and the use of subject-specific flow waveform on the CFD results. We quantified these differences by examining the pointwise percent error of the mean wall shear stress (WSS) and the oscillatory shear index (OSI) and by computing the intra-class correlation coefficient (ICC) between axial profiles of the mean WSS and OSI in the internal carotid artery bulb. The parabolic inlet velocity profile produced the most similar mean WSS and OSI to simulations employing the real patient-specific inlet velocity profile. However, anatomic variation in vessel geometry and the use of a nonpatient-specific flow waveform both affected the WSS and OSI results more than did the choice of inlet velocity profile. Although careful selection of boundary conditions is essential for all CFD analysis, accurate patient-specific geometry reconstruction and measurement of vessel flow rate waveform are more important than the choice of velocity profile. A parabolic velocity profile provided results most similar to the patient-specific velocity profile.     

In-vitro pulsatile flow visualization studies in a pulmonary artery model

In-vitro pulsatile flow visualization studies were conducted in an adult-sized pulmonary artery model to observe the effects of valvular pulmonic stenosis on the flow fields of the main, left and right pulmonary arteries. The flow patterns revealed that as the degree of stenosis increased, the jet-type flow created by the valve became narrower, and it impinged on the far (distal) wall of the left pulmonary artery further downstream from the junction of the bifurcation. This in turn led to larger regions of disturbed turbulent flow, as well as helical-type secondary flow motions in the left pulmonary artery, compared to the right pulmonary artery. The flow field in the main pulmonary artery also became more disturbed and turbulent, especially during peak systole and the deceleration phase. The flow visualization observations have been valuable in helping to conduct further quantitative studies such as pressure and velocity field mapping. Such studies are important to understanding the fluid mechanics characteristics of the main pulmonary artery and its two major branches.     

Blood velocity patterns in poststenotic regions and velocity waveforms for myocardial inflow associated with coronary artery stenosis in dogs.

Velocity profiles across a vessel were investigated in poststenotic regions of the canine left coronary artery by our 80-channel 20 MHz ultrasound velocimeter. The velocity waveform in a small artery just before its penetration into myocardium was measured by our laser Doppler method. The poststenotic velocity configuration was characterized by a narrow region of high velocity with diastolic reverse flow near the wall which may dissipate energy. The velocity waveform in the distal small arteries exhibited increased systolic reverse flow with decreased diastolic forward flow, resulting in a remarkable reduction of coronary inflow into the myocardium.     

Numerical Analysis of Blood Flow through COVID-19 Infected Arteries

Computational Fluid Dynamics has become relevant in the study of hemodynamics, where clinical results are challenging to obtain. This paper discusses a 2-Dimensional transient blood flow analysis through an arterial bifurcation for patients infected with the Coronavirus. The geometry considered is an arterial bifurcation with main stem diameter 3 mm and two outlets. The left outlet (smaller) has a diameter of 1.5 mm and the right outlet (larger), 2 mm. The length of the main stem, left branch and right branch are fixed at 35 mm, 20 mm and 25 mm respectively. Viscosity change that occurs in the blood leads to different parametrical changes in blood flow. The blood flow towards the smaller branch is significantly affected by the changed blood viscosity. Extended regions of high pressure and increased velocity towards the larger outlet are obtained. The Time Averaged Wall Shear Stress (TAWSS) for the corona affected artery is found to be 10.4114 Pa at a 90° angle of bifurcation as compared to 2.45002 Pa of the normal artery. For varying angles of bifurcation, an angle of 75° was found to have a maximum Time Averaged Wall Shear Stress of 2.46076 Pa and 10.42542 Pa for normal and corona affected artery, respectively.     

Velocity distribution along an elastic model of human arterial tree

An experimental investigation of an elastic model of the human arterial tree, has been performed for physiological type flow by pulsed Doppler ultrasonic velocimetry. The arterial tree model, fabricated in clear polyurethane, includes the aortic arch, with a Starr-Edwards ball valve mounted in the root of the aorta, the descending aorta and the iliac bifurcation. Our study showed that the velocity profile, a few centimeters beyond the valve, is skewed, with higher velocities towards the top and the inner wall (anatomically the posterior and left lateral wall). An inward shift of the maximum velocity and reverse flow are denoted along the inner wall of the aortic arch. The velocity profiles in the descending aorta are blunted. Downstream from the vertex of the iliac bifurcation, there is vorticity creation, but the branching effect is quickly damped by the pulsatility of the flow and the elasticity of the wall.     

Which diameter and angle rule provides optimal flow patterns in a coronary bifurcation?

The branching angle and diameter ratio in epicardial coronary artery bifurcations are two important determinants of atherogenesis. Murray's cubed diameter law and bifurcation angle have been assumed to yield optimal flows through a bifurcation. In contrast, we have recently shown a 7/3 diameter law (HK diameter model), based on minimum energy hypothesis in an entire tree structure. Here, we derive a bifurcation angle rule corresponding to the HK diameter model and critically evaluate the streamline flow through HK and Murray-type bifurcations. The bifurcations from coronary casts were found to obey the HK diameter model and angle rule much more than Murray's model. A finite element model was used to investigate flow patterns for coronary artery bifurcations of various types. The inlet velocity and pressure boundary conditions were measured by ComboWire. Y-bifurcation of Murray type decreased wall shear stress-WSS (10%-40%) and created an increased oscillatory shear index-OSI in atherosclerosis-prone regions as compared with HK-type bifurcations. The HK-type bifurcations were found to have more optimal flow patterns (i.e., higher WSS and lower OSI) than Murray-type bifurcations which have been traditionally believed to be optimized. This study has implications for changes in bifurcation angles and diameters in percutaneous coronary intervention.     

Quantification of absolute coronary flow reserve and relative fractional flow reserve in a swine animal model using angiographic image data

Coronary flow reserve (CFR) and fractional flow reserve (FFR) are important physiological indexes for coronary disease. The purpose of this study was to validate the CFR and FFR measurement techniques using only angiographic image data. Fifteen swine were instrumented with an ultrasound flow probe on the left anterior descending artery (LAD). Microspheres were gradually injected into the LAD to create microvascular disruption. An occluder was used to produce stenosis. Contrast material injections were made into the left coronary artery during image acquisition. Volumetric blood flow from the flow probe (Q(q)) was continuously recorded. Angiography-based blood flow (Q(a)) was calculated by using a time-density curve based on the first-pass analysis technique. Flow probe-based CFR (CFR(q)) and angiography-based CFR (CFR(a)) were calculated as the ratio of hyperemic to baseline flow using Q(q) and Q(a), respectively. Relative angiographic FFR (relative FFR(a)) was calculated as the ratio of the normalized Q(a) in LAD to the left circumflex artery (LC(X)) during hyperemia. Flow probe-based FFR (FFR(q)) was measured from the ratio of hyperemic flow with and without disease. CFR(a) showed a strong correlation with the gold standard CFR(q) (CFR(a) = 0.91 CFR(q) + 0.30; r = 0.90; P < 0.0001). Relative FFR(a) correlated linearly with FFR(q) (relative FFR(a) = 0.86 FFR(q) + 0.05; r = 0.90; P < 0.0001). The quantification of CFR and relative FFR(a) using angiographic image data was validated in a swine model. This angiographic technique can potentially be used for coronary physiological assessment during routine cardiac catheterization.     

Flow measurements in a human femoral artery model with reverse lumen curvature

Flow visualization and wall pressure measurements were made in a smooth reverse curvature model that conformed to the gentle "s" shape of a left femoral artery angiogram of a patient in a clinical trial. Observed lesion localization at the inner (lesser) curvatures appeared to be associated with secondary flows in the wall vicinity directed toward the inner curvatures that tended to reverse direction in the flow entering the reverse curvature region. Moderate flow resistance increases of about 20 percent above the Poiseuille flow relation were found at the higher physiological Reynolds numbers Re above about 600-700 and thus Dean numbers for steady flow. For pulsatile flow simulation, flow resistances did not increase up to the largest Re of 470 tested. Apparently, the large variations in velocity during the cardiac cycle disrupted the stronger secondary flow patterns observed at the higher Reynolds numbers for steady flow.     

Primary percutaneous coronary intervention in the left main stem of a monocoronary artery

In a 71-year-old female with evolving anterior wall myocardial infarction, coronary angiography revealed a monocoronary artery which arose from the right sinus of Valsalva. Originating from a short common trunk, the left main stem showed a thrombotic lesion that occluded the left anterior descending coronary artery while the circumflex artery was obstructed. Intracoronary administration of abciximab, followed by stenting of the transition between the left anterior descending coronary artery and the main stem, and final kissing balloon inflation of the bifurcation resulted in an excellent angiographic result and favourable clinical outcome. (Neth Heart J 2009;17:274–6.)     

Steady flow velocity measurements in a pulmonary artery model with varying degrees of pulmonic stenosis

Velocity and flow visualization studies were conducted in an adult size pulmonary artery model with varying degrees of valvular stenosis, using a two dimensional laser Doppler anemometer system. Velocity measurements in the main, left and right branches of the pulmonary artery revealed that as the degree of pulmonic stenosis increased, the jet type flow created by the valve hit the distal wall of the LPA farther downstream from the junction of the bifurcation. This in turn led to higher levels of turbulent and disturbed flow, and larger secondary flow motion in the LPA compared to the RPA. The high levels of turbulence measured in the main and left pulmonary arteries with the stenotic valves, could lead to the clinically observed phenomenon of post stenotic dilatation in the MPA extending into the LPA.     

Comparison of intima-media thickness and ophthalmic artery resistance index for assessing subclinical atherosclerosis in HIV-1-infected patients

Background

The determination of coronary flow reserve (CFR) is an essential concept at the moment of decision-making in ischemic heart disease. There are several direct and indirect tests to evaluate this parameter. In this sense, dobutamine stress echocardiography is one of the pharmacological method most commonly used worldwide. It has been previously demonstrated that CFR can be determined by this technique. Despite our wide experience with dobutamine stress echocardiography, we ignored the necessary heart rate to consider sufficient the test for the analysis of CFR. For this reason, our main goal was to determine the velocity of coronary flow in each stage of dobutamine stress echocardiography and the heart rate value necessary to double the baseline values of coronary flow velocity in the territory of the left anterior descending (LAD) coronary artery.

Methods

A total of 33 consecutive patients were analyzed. The patients included had low risk for coronary artery disease. All the participants underwent dobutamine stress echocardiography and coronary artery flow velocity was evaluated in the distal segment of LAD coronary artery using transthoracic color-Doppler echocardiography.

Results

The feasibility of determining CFR in the territory of the LAD during dobutamine stress echocardiography was high: 31/33 patients (94%). Mean CFR was 2.67 at de end of dobutamine test. There was an excellent concordance between delta HR (difference between baseline HR and maximum HR) and the increase in the CFR (correlation coefficient 0.84). In this sense, we found that when HR increased by 50 beats, CFR was ≥ 2 (CI 93-99.2%). In addition, 96.4% of patients reached a CFR ≥ 2 (IC 91.1 - 99%) at 75% of their predicted maximum heart rate.

Conclusions

We found that the feasibility of dobutamine stress echocardiography to determine CFR in the territory of the LAD coronary artery was high. In this study, it was necessary to achieve a difference of 50 bpm from baseline HR or at least 75% of the maximum predicted heart rate to consider sufficient the test for the analysis of CFR.     

The effects of non-Newtonian viscoelasticity and wall elasticity on flow at a 90 degrees bifurcation

To study the fundamentals of hemodynamics in arteries, the flow parameters: pulsatility, elasticity and non-Newtonian viscoelasticity were considered in detail in a 90 degrees-T-bifurcation of a rigid and elastic model. The velocity distribution 2.5 mm behind the bifurcation in the straight tube was measured with a laser-Doppler-anemometer. The fluid used was an aqueous glycerine solution and a viscoelastic Separan mixture. Flow visualization studies were done with a sheet of laser light in the plane of the bifurcation. The velocity distribution was measured for both steady and pulsatile flows with a laser-Doppler-anemometer in a backward scattered way. From the velocity measurements the shear gradients were calculated. Substantial differences were found in the flow behavior of Newtonian and non-Newtonian fluids, especially behind the bifurcation in the main tube, where secondary flows and flow separation started. Also, differences due to the elastic and rigid wall could be seen. Very high shear gradients were found in the flow between main flow and the separation zone which can lead to a damage of the blood cells.     

Computation of hemodynamics in the left coronary artery with variable angulations

The purpose of this study was to investigate the hemodynamic effect of variations in the angulations of the left coronary artery, based on simulated and realistic coronary artery models. Twelve models consisting of four realistic and eight simulated coronary artery geometries were generated with the inclusion of left main stem, left anterior descending and left circumflex branches. The simulated models included various coronary artery angulations, namely, 15°, 30°, 45°, 60°, 75°, 90°, 105° and 120°. The realistic coronary angulations were based on selected patient's data with angles ranging from narrow angles of 58° and 73° to wide angles of 110° and 120°. Computational fluid dynamics analysis was performed to simulate realistic physiological conditions that reflect the in vivo cardiac hemodynamics. The wall shear stress, wall shear stress gradient, velocity flow patterns and wall pressure were measured in simulated and realistic models during the cardiac cycle. Our results showed that a disturbed flow pattern was observed in models with wider angulations, and wall pressure was found to reduce when the flow changed from the left main stem to the bifurcated regions, based on simulated and realistic models. A low wall shear stress gradient was demonstrated at left bifurcations with wide angles. There is a direct correlation between coronary angulations and subsequent hemodynamic changes, based on realistic and simulated models. Further studies based on patients with different severities of coronary artery disease are required to verify our results.     

Assessment of atrial regional and global electromechanical function by tissue velocity echocardiography: a feasibility study on healthy individuals

Background

Myocardial contrast echocardiography and coronary flow velocity pattern with a rapid diastolic deceleration time after percutaneous coronary intervention has been reported to be useful in assessing microvascular damage in patients with acute myocardial infarction.

Aim

To evaluate myocardial contrast echocardiography with harmonic power Doppler imaging, coronary flow velocity reserve and coronary artery flow pattern in predicting functional recovery by using transthoracic echocardiography.

Methods

Thirty patients with anterior acute myocardial infarction underwent myocardial contrast echocardiography at rest and during hyperemia and were quantitatively analyzed by the peak color pixel intensity ratio of the risk area to the control area (PIR). Coronary flow pattern was measured using transthoracic echocardiography in the distal portion of left anterior descending artery within 24 hours after recanalization and we assessed deceleration time of diastolic flow velocity. Coronary flow velocity reserve was calculated two weeks after acute myocardial infarction. Left ventricular end-diastolic volumes and ejection fraction by angiography were computed.

Results

Pts were divided into 2 groups according to the deceleration time of coronary artery flow pattern (Group A; 20 pts with deceleration time ≧ 600 msec, Group B; 10 pts with deceleration time < 600 msec). In acute phase, there were no significant differences in left ventricular end-diastolic volume and ejection fraction (Left ventricular end-diastolic volume 112 ± 33 vs. 146 ± 38 ml, ejection fraction 50 ± 7 vs. 45 ± 9 %; group A vs. B). However, left ventricular end-diastolic volume in Group B was significantly larger than that in Group A (192 ± 39 vs. 114 ± 30 ml, p < 0.01), and ejection fraction in Group B was significantly lower than that in Group A (39 ± 9 vs. 52 ± 7%, p < 0.01) at 6 months. PIR and coronary flow velocity reserve of Group A were higher than Group B (PIR, at rest: 0.668 ± 0.178 vs. 0.248 ± 0.015, p < 0.0001: during hyperemia 0.725 ± 0.194 vs. 0.295 ± 0.107, p < 0.0001; coronary flow velocity reserve, 2.60 ± 0.80 vs. 1.31 ± 0.29, p = 0.0002, respectively).

Conclusion

The preserved microvasculature detecting by myocardial contrast echocardiography and coronary flow velocity reserve is related to functional recovery after acute myocardial infarction.