In vivo demonstration of flow recirculation and turbulence downstream of graded stenoses in canine arteries

The velocity field around arterial stenoses was investigated using a pulsed doppler velocimeter in vivo. Asymmetric zones of recirculation were identified by systolic flow reversal in the post-stenotic field in carotid and iliac arteries of anesthetised dogs. There was a close correlation between shear intensity and turbulence as estimated by the velocity difference between the jet and the recirculation zone and by maximum spectral width respectively. Under the conditions of these experiments, stenosis grade (% diameter reduction) dominated hemodynamic variables such as Reynolds number, oscillation and pulsatility in determining the intensity of turbulence. The method used does not appear to have sufficient resolution to distinguish between turbulence and discrete oscillating velocities (vorticity) nor to allow determination of wall shear stress though the pattern of change of the latter is similar to that found downstream of axisymmetric stenosis in models using steady flow.     

Numerical modeling of pulsatile turbulent flow in stenotic vessels

Pulsatile turbulent flow in stenotic vessels has been numerically modeled using the Reynolds-averaged Navier-Stokes equation approach. The commercially available computational fluid dynamics code (CFD), FLUENT, has been used for these studies. Two different experiments were modeled involving pulsatile flow through axisymmetric stenoses. Four different turbulence models were employed to study their influence on the results. It was found that the low Reynolds number k-omega turbulence model was in much better agreement with previous experimental measurements than both the low and high Reynolds number versions of the RNG (renormalization-group theory) k-epsilon turbulence model and the standard k-epsilon model, with regard to predicting the mean flow distal to the stenosis including aspects of the vortex shedding process and the turbulent flow field. All models predicted a wall shear stress peak at the throat of the stenosis with minimum values observed distal to the stenosis where flow separation occurred.     

Flow disturbance measurements through a constricted tube at moderate Reynolds numbers

Instantaneous velocities in the field distal to contoured axisymmetric stenoses were measured with a laser Doppler anemometer. Upstream flow conditions were steady and spanned a range of Reynolds numbers from 500 to 2000. Autocorrelation functions and spectra of the velocity were employed to describe the nature of fluid dynamic disturbances. Depending upon the degree of stenosis and the Reynolds number, the flow field contained disturbances of a discrete oscillation frequency, of a turbulent nature, or both. If turbulence was detected in a given experiment, it was always preceded upstream by velocity oscillations at discrete frequency arising from vortex shedding. For mild degrees of stenosis (50% area reduction or less) the intensity of flow disturbances was relatively low until the Reynolds number exceeded 1000, thus highlighting difficulties to be expected in employing flow disturbance detection as a diagnostic tool in the recognition of early atherosclerosis in major arteries. In view of the relatively high noise levels inherent in noninvasive Doppler ultrasound systems employed clinically, it seems unlikely that detection of stenoses of less than 50% area reduction is feasible unless the Reynolds numbers exceed 1000 or unless pulsatility introduces new unsteady flow features beyond those studied here.     

Comparison of LES of steady transitional flow in an idealized stenosed axisymmetric artery model with a RANS transitional model

In this study, two different turbulence methodologies are investigated to predict transitional flow in a 75% stenosed axisymmetric experimental arterial model and in a slightly modified version of the model with an eccentric stenosis. Large eddy simulation (LES) and Reynolds-averaged Navier-Stokes (RANS) methods were applied; in the LES simulations eddy viscosity subgrid-scale models were employed (basic and dynamic Smagorinsky) while the RANS method involved the correlation-based transitional version of the hybrid k-ε/k-ω flow model. The RANS simulations used 410,000 and 820,000 element meshes for the axisymmetric and eccentric stenoses, respectively, with y(+) less than 2 viscous wall units for the boundary elements, while the LES used 1,200,000 elements with y(+) less than 1. Implicit filtering was used for LES, giving an overlap between the resolved and modeled eddies, ensuring accurate treatment of near wall turbulence structures. Flow analysis was carried out in terms of vorticity and eddy viscosity magnitudes, velocity, and turbulence intensity profiles and the results were compared both with established experimental data and with available direct numerical simulations (DNSs) from the literature. The simulation results demonstrated that the dynamic Smagorinsky LES and RANS transitional model predicted fairly comparable velocity and turbulence intensity profiles with the experimental data, although the dynamic Smagorinsky model gave the best overall agreement. The present study demonstrated the power of LES methods, although they were computationally more costly, and added further evidence of the promise of the RANS transition model used here, previously tested in pulsatile flow on a similar model. Both dynamic Smagorinsky LES and the RANS model captured the complex transition phenomena under physiological Reynolds numbers in steady flow, including separation and reattachment. In this respect, LES with dynamic Smagorinsky appeared more successful than DNS in replicating the axisymmetric experimental results, although inflow conditions, which are subject to caveats, may have differed. For the eccentric stenosis, LES with Smagorinsky coefficient of 0.13 gave the closest agreement with DNS despite the known shortcomings of fixed coefficients. The relaminarization as the flow escaped the influence of the stenosis was amply demonstrated in the simulations, graphically so in the case of LES.     

Velocity field studies at surgically imposed arterial stenoses on the abdominal aorta in pigs

In order to describe velocity profiles and the size of deterministic and non-deterministic velocity disturbances at arterial stenoses, symmetrical and asymmetrical stenoses with intended area reductions of 50% (‘moderate’) and 85% (‘severe’) were applied on the abdominal aorta in six pigs. Blood velocities were registered by hot-film anemometry in 21 measuring points distributed across the vessel cross-sectional area in one pre-stenotic and three post-stenotic positions. Signal analysis included ensemble averaging, the high-pass filtering technique, and three-dimensional visualization. None of the stenoses affected the pre-stenotic velocity field. Downstream moderate stenoses flow separation and vortex formation were present. Moderate asymmetric stenoses induced turbulence in the post-stenotic velocity field. Immediately downstream of severe stenoses a prominent post-stenotic jet was present. Farther downstream, a multitude of coherent vortices and turbulence dominated the flow field. The transverse distribution of turbulence intensity parallelled with the peak systolic velocity profile, whereas transverse profiles of the relative turbulence intensity (turbulence intensity/mean velocity) revealed peak values in flow field locations with high velocity gradients. Velocity parameters for symmetric and asymmetric severe stenoses were highly comparable. However, the exact degree of stenosis was significantly higher for symmetrical (85%) than for asymmetrical (76%) stenoses. Therefore, recalling that stenosis severity strongly influences the development of velocity disturbances, this indicates that asymmetry of a stenosis is a predictor for blood velocity disturbances.     

Low Reynolds number turbulence modeling of blood flow in arterial stenoses.

Moderate and severe arterial stenoses can produce highly disturbed flow regions with transitional and or turbulent flow characteristics. Neither laminar flow modeling nor standard two-equation models such as the kappa-epsilon turbulence ones are suitable for this kind of blood flow. In order to analyze the transitional or turbulent flow distal to an arterial stenosis, authors of this study have used the Wilcox low-Re turbulence model. Flow simulations were carried out on stenoses with 50, 75 and 86% reductions in cross-sectional area over a range of physiologically relevant Reynolds numbers. The results obtained with this low-Re turbulence model were compared with experimental measurements and with the results obtained by the standard kappa-epsilon model in terms of velocity profile, vortex length, wall shear stress, wall static pressure, and turbulence intensity. The comparisons show that results predicted by the low-Re model are in good agreement with the experimental measurements. This model accurately predicts the critical Reynolds number at which blood flow becomes transitional or turbulent distal an arterial stenosis. Most interestingly, over the Re range of laminar flow, the vortex length calculated with the low-Re model also closely matches the vortex length predicted by laminar flow modeling. In conclusion, the study strongly suggests that the proposed model is suitable for blood flow studies in certain areas of the arterial tree where both laminar and transitional/turbulent flows coexist.     

Direct numerical simulation of transitional flow in a stenosed carotid bifurcation

The blood flow dynamics of a stenosed, subject-specific, carotid bifurcation were numerically simulated using the spectral element method. Pulsatile inlet conditions were based on in vivo color Doppler ultrasound measurements of blood velocity. The results demonstrated the transitional or weakly turbulent state of the blood flow, which featured rapid velocity and pressure fluctuations in the post-stenotic region of the internal carotid artery (ICA) during systole and laminar flow during diastole. High-frequency vortex shedding was greatest downstream of the stenosis during the deceleration phase of systole. Velocity fluctuations had a frequency within the audible range of 100-300Hz. Instantaneous wall shear stress (WSS) within the stenosis was relatively high during systole ( approximately 25-45Pa) compared to that in a healthy carotid. In addition, high spatial gradients of WSS were present due to flow separation on the inner wall. Oscillatory flow reversal and low pressure were observed distal to the stenosis in the ICA. This study predicts the complex flow field, the turbulence levels and the distribution of the biomechanical stresses present in vivo within a stenosed carotid artery.     

Experimental Study of the Vortex-Induced Vibration of Drilling Risers under the Shear Flow with the Same Shear Parameter at the Different Reynolds Numbers

A considerable number of studies for VIV under the uniform flow have been performed. However, research on VIV under shear flow is scarce. An experiment for VIV under the shear flow with the same shear parameter at the two different Reynolds numbers was conducted in a deep-water offshore basin. Various measurements were obtained by the fiber bragg grating strain sensors. Experimental data were analyzed by modal analysis method. Results show several valuable features. First, the corresponding maximum order mode of the natural frequency for shedding frequency is the maximum dominant vibration mode and multi-modal phenomenon is appeared in VIV under the shear flow, and multi-modal phenomenon is more apparent at the same shear parameter with an increasing Reynolds number under the shear flow effect. Secondly, the riser vibrates at the natural frequency and the dominant vibration frequency increases for the effect of the real-time tension amplitude under the shear flow and the IL vibration frequency is the similar with the CF vibration frequency at the Reynolds number of 1105 in our experimental condition and the IL dominant frequency is twice the CF dominant frequency with an increasing Reynolds number. In addition, the displacement trajectories at the different locations of the riser appear the same shape and the shape is changed at the same shear parameter with an increasing Reynolds number under the shear flow. The diagonal displacement trajectories are observed at the low Reynolds number and the crescent-shaped displacement trajectories appear with an increasing Reynolds number under shear flow in the experiment.     

Attenuation of flow disturbances in tapered arterial grafts

Flow disturbances in tapered arterial grafts of angles of taper between 0.5 and 1.0 deg were measured in vitro using a pulsed ultrasound Doppler velocimeter. The increase in transition Reynolds numbers with angle of taper and axial distance was determined for steady flow. The instantaneous centerline velocities were measured distal to a 50 percent area stenosis (as a model of a proximal anastomosis), in steady and pulsatile flow, from which the disturbance intensities were calculated. A significant reduction in post-stenotic disturbance intensity was recorded in the tapered grafts, relative to a conventional cylindrical graft. In pulsatile flow with a large backflow component, however, there was an increase in disturbance intensity due to diverging flow during flow reversal. This was observed only in the 1.0 deg tapered graft. These findings indicate that taper is an important consideration in the design of vascular prostheses.     

A numerical and experimental investigation of transitional pulsatile flow in a stenosed channel

In the present paper, a closely coupled numerical and experimental investigation of pulsatile flow in a prototypical stenotic site is presented. Detailed laser Doppler velocimetry measurements upstream of the stenosis are used to guide the specification of velocity boundary conditions at the inflow plane in a series of direct numerical simulations (DNSs). Comparisons of the velocity statistics between the experiments and DNS in the post-stenotic area demonstrate the great importance of accurate inflow conditions, and the sensitivity of the post-stenotic flow to the disturbance environment upstream. In general, the results highlight a borderline turbulent flow that sequentially undergoes transition to turbulence and relaminarization. Before the peak mass flow rate, the strong confined jet that forms just downstream of the stenosis becomes unstable, forcing a role-up and subsequent breakdown of the shear layer. In addition, the large-scale structures originating from the shear layer are observed to perturb the near wall flow, creating packets of near wall hairpin vortices.     

Distribution of mass transfer rate and wall shear stress behind simple rectangular stenosis in pulsating flow

The distributions of mass transfer rate and wall shear stress in sinusoidal laminar pulsating flow through a two-dimensional asymmetric stenosed channel have been studied experimentally and numerically. The distributions are measured by the electrochemical method. The measurement is conducted at a Reynolds number of about 150, a Schmidt number of about 1000, a nondimensional pulsating frequency of 3.40, and a nondimensional flow amplitude of 0.3. It is suggested that the deterioration of an arterial wall distal to stenosis may be greatly enhanced by fluid dynamic effects.     

Disorder distal to modeled stenoses in steady and pulsatile flow

Measurements of the velocity and energy spectra were made in the distal region of modeled stenoses in a rigid tube with both steady and pulsatile water flows. Reynolds numbers of 318–2540 and a pulsatile flow frequency parameter of 15 were employed. The effects of the degree of stenosis, the stenosis geometry and the presence or absence of the downstream confining wall on the development of flow disturbances were investigated. Visualization of the distal flow patterns in stenotic and free jets illustrated the existence of complex fields which included vortex shedding, highly turbulent regions, and recirculation zones. Significant flow disorder was created by a mild stenosis in pulsatile, but not in steady, flow. Nondimensionalization employing the stenosis diameter and flow velocity in the throat of the constriction correlates the vortex shedding frequency and energy spectra within a limited postestenotic region.     

Mathematical modelling of flow through an irregular arterial stenosis.

A mathematical model of flow through an irregular arterial stenosis is developed. The model is two-dimensional and axi-symmetric with the stenosis outline obtained from a three-dimensional casting of a mildly stenosed artery. Agreement between modelled and experimental pressure drops (obtained from an axi-symmetric machined stenosis with the same profile) is excellent. Results are also obtained for a smooth stenosis model, similar to that used for most mathematical modelling studies. This model overestimates the pressure drop across the stenosis, as well as the wall shear stress and separation Reynolds number. Also, the smooth model predicts one instead of three recirculation zones present in the irregular model. The original stenosis is modified to increase the severity from 48 and 87% areal occlusion, while maintaining the same general shape. This has the effect of increasing the pressure drop by an order of magnitude and decreasing the number of recirculation zones to one, with a lower separation Reynolds number.     

Effects of vibration on steady flow downstream of a stenosis

A physical model was used to test the effects of vibrations on the position of transition to turbulence (zt) downstream of a constriction. Constrictions were inserted in a length of clear plastic (Tygon) tubing and vibrated. Water was the fluid medium and flow was visualized with India ink. The frequency and velocity of vibration were monitored. For a given constriction and flow velocity, there was a band of frequencies which caused zt to move upstream. This band corresponded to frequencies of flow disturbances as measured with a hot-film anemometer without vibration. Both flow-visualization and hot-film frequencies were correlated via a Strouhal number to the Reynolds number and contraction ratio of the flow. Values of zt decreased with increasing vibration amplitude. The critical Reynolds number for turbulence was also decreased by vibration. These results are of importance in the diagnosis of vascular disease and the design of physical models of stenotic flows.     

Correction to: Characterization of anisotropic turbulence behavior in pulsatile blood flow

Turbulent-like hemodynamics with prominent cycle-to-cycle flow variations have received increased attention as a potential stimulus for cardiovascular diseases. These turbulent conditions are typically evaluated in a statistical sense from single scalars extracted from ensemble-averaged tensors (such as the Reynolds stress tensor), limiting the amount of information that can be used for physical interpretations and quality assessments of numerical models. In this study, barycentric anisotropy invariant mapping was used to demonstrate an efficient and comprehensive approach to characterize turbulence-related tensor fields in patient-specific cardiovascular flows, obtained from scale-resolving large eddy simulations. These techniques were also used to analyze some common modeling compromises as well as MRI turbulence measurements through an idealized constriction. The proposed method found explicit sites of elevated turbulence anisotropy, including a broad but time-varying spectrum of characteristics over the flow deceleration phase, which was different for both the steady inflow and Reynolds-averaged Navier–Stokes modeling assumptions. Qualitatively, the MRI results showed overall expected post-stenotic turbulence characteristics, however, also with apparent regions of unrealizable or conceivably physically unrealistic conditions, including the highest turbulence intensity ranges. These findings suggest that more detailed studies of MRI-measured turbulence fields are needed, which hopefully can be assisted by more comprehensive evaluation tools such as the once described herein.

     

Numerical analysis of flow through a severely stenotic carotid artery bifurcation

The results of computational simulations may supplement MR and other in vivo diagnostic techniques to provide an accurate picture of the hemodynamics in particular vessels, which may help demonstrate the risks of embolism or plaque rupture posed by particular plaque deposits. In this study, a model based on an endarterectomy specimen of the plaque in a carotid bifurcation was examined. The flow conditions include steady flow at Reynolds numbers of 300, 600, and 900 as well as unsteady pulsatile flow. Both dynamic pressure and wall shear stress are very high, with shear values up to 70 N/m2, proximal to the stenosis throat in the internal carotid artery, and both vary significantly through the flow cycle. The wall shear stress gradient is also strong along the throat. Vortex shedding is observed downstream of the most severe occlusion. Two turbulence models, the Chien and Goldberg varieties of k-epsilon, are tested and evaluated for their relevance in this geometry. The Chien model better captures phenomena such as vortex shedding. The flow distal to stenosis is likely transitional, so a model that captures both laminar and turbulent behavior is needed.     

Velocity profiles and streamlines of a revolution post-stenotic flow.

Velocity fields have been measured in models simulating arterial stenoses for continuous and revolution flows. A pulsed Doppler velocimeter allows for velocity readings in the entire tube and in the wall area. Streamlines are determined by numerical solving of the system of equations defining the current function. Velocity profiles and streamlines are presented and discussed either for steady or for unsteady flows, with different Reynolds numbers and variable degrees of stenosis. There is, in the wall area, a recirculating zone made of a well-defined rouleau. Its length varies increasingly according to the increasing severity of the stenosis. The stability of axial flow depends on the input profile, the degree of stenosis and the Reynolds number. Plotting streamlines allows to describe accurately the flow; its quantitative aspect offers advantages with respect to conventional visualization mode.     

Velocity measurements in steady flow through axisymmetric stenoses at moderate Reynolds numbers

The velocity field in the neighborhood of axisymmetric constrictions in rigid tubes was investigated using laser Doppler anemometry and flow visualization. Upstream flow conditions were steady; and Reynolds numbers were in the range 500-2000, values which are representative of the larger arteries in humans. Stenoses of 25, 50 and 75% area reduction were studied. Velocity profiles are presented in sufficient detail to allow comparison with computational biofluid dynamics models. Wall shear stresses were estimated from the near wall velocity gradient, and the nature of observed poststenotic flow disturbances is discussed. Results indicate that flow disturbances of discrete oscillation frequency may be more valuable than turbulence as an indicator of early stages of stenosis development. Additionally, despite the fact that poststenotic turbulence exists for the higher degrees of stenosis and Reynolds numbers, the resulting wall shear stresses are only three to four times greater than the Poiseuille value and are considerably less than the wall shear stress within the stenosis itself.     

To what extent does a minimal atherosclerotic plaque alter the arterial wall shear stress distribution? A model study by an electrochemical method

An electrochemical surface shear stress measurement was applied to a model of very thin unilateral arterial stenosis (height of 1/8 of the model pipe diameter with very smooth surface). Three dimensional wall shear stress distribution was measured under steady flow field from a relatively low Reynolds number, Re = 270, to a high Reynolds number, Re = 1200. There was a characteristic high and low wall shear distribution pattern around the stenosis. There were also remarkable high shear stress areas on the opposite wall and both side walls of the stenosis. It was clearly shown that three dimensional structure of the flow field, hence, the wall shear stress distribution, is affected by a minimal change on the arterial wall.