Oscillatory flow and gas transport through a symmetrical bifurcation

Axial gas transport due to the interaction between radial mixing and radially nonuniform axial velocities is responsible for gas transport in thick airways during High-frequency oscillatory ventilation (HFO). Because the airways can be characterized by a bifurcating tube network, the secondary flow in the curved portion of a bifurcating tube contributes to cross-stream mixing. In this study the oscillatory flow and concentration fields through a single symmetrical airway bifurcating tube model were numerically analyzed by solving three-dimensional Navier-Stokes and mass concentration equations with the SIMPLER algorithm. The simulation conditions were for a Womersley number, alpha = 9.1 and Reynolds numbers in the parent tube between 200 and 1000, corresponding to Dn2/alpha 4 in the curved portion between 2 and 80, where Dn is Dean number. For comparison with the results from the bifurcating tube, we calculated the velocity and concentration fields for fully developed oscillatory flow through a curved tube with a curvature rate of 1/10, which is identical to the curved portion of the bifurcating tube. For Dn2/alpha 4 < or = 10 in the curved portion of the bifurcating tube, the flow divider and area changes dominate the axial gas transport, because the effective diffusivity is greater than in either a straight or curved tube, in spite of low secondary velocities. However, for Dn2/alpha 4 > or = 20, the gas transport characteristics in a bifurcation are similar to a curved tube because of the significant effect of secondary flow.     

An experimental analysis of the flow field in a three-dimensional model of the human carotid artery bifurcation

Steady flow measurements were carried out in a rigid three-dimensional model of the human carotid artery bifurcation at a Reynolds number of 640 and a flow division ratio of 50/50. Both axial and secondary velocities were measured with a laser-Doppler anemometer. In the bulb opposite to the flow divider a zone with negative axial velocities was found with a maximal diameter of about 60% of the local diameter of the branch and a cross-sectional extent of about 25% of the local cross-sectional area. In the bulb the maximum axial velocity shifted towards the divider wall and at the end of the bulb an axial velocity plateau arose near the non-divider wall. Halfway through the bulb, secondary flow showed a vortex through which fluid flowed towards the divider wall near the bifurcation plane and back towards the non-divider wall near the upper walls.     

Development of wall surface tangent DPIV measurement techniques for arterial branch models.

Experimental techniques for measuring unsteady flow in a glass arterial bifurcation model have been developed to aid in quantifying three-dimensional wall shear fluctuations associated with arterial disease. The unique feature of the current technique is the use of a "curved" laser sheet, which was everywhere tangent to the inner wall of a daughter tube in an arterial bifurcation model. Surface tangent velocity vector field measurements were made to demonstrate the potential of this technique. Ensemble-averaged data showing weak secondary flows as well as statistical distributions of flow angles are presented. Measurements of this type may be used to estimate mean and instantaneous wall shear magnitude and direction, data that are necessary for understanding the importance of circumferential motions on arterial disease.     

Exploring high frequency temporal fluctuations in the terminal aneurysm of the basilar bifurcation

Cerebral aneurysms are a common cause of death and disability. Of all the cardiovascular diseases, aneurysms are perhaps the most strongly linked with the local fluid mechanic environment. Aside from early in vivo clinical work that hinted at the possibility of high-frequency intra-aneurysmal velocity oscillations, flow in cerebral aneurysms is most often assumed to be laminar. This work investigates, through the use of numerical simulations, the potential for disturbed flow to exist in the terminal aneurysm of the basilar bifurcation. The nature of the disturbed flow is explored using a series of four idealized basilar tip models, and the results supported by four patient specific terminal basilar tip aneurysms. All four idealized models demonstrated instability in the inflow jet through high frequency fluctuations in the velocity and the pressure at approximately 120?Hz. The instability arises through a breakdown of the inflow jet, which begins to oscillate upon entering the aneurysm. The wall shear stress undergoes similar high-frequency oscillations in both magnitude and direction. The neck and dome regions of the aneurysm present 180 deg changes in the direction of the wall shear stress, due to the formation of small recirculation zones near the shear layer of the jet (at the frequency of the inflow jet oscillation) and the oscillation of the impingement zone on the dome of the aneurysm, respectively. Similar results were observed in the patient-specific models, which showed high frequency fluctuations at approximately 112 Hz in two of the four models and oscillations in the magnitude and direction of the wall shear stress. These results demonstrate that there is potential for disturbed laminar unsteady flow in the terminal aneurysm of the basilar bifurcation. The instabilities appear similar to the first instability mode of a free round jet.     

A numerical analysis of steady flow in a three-dimensional model of the carotid artery bifurcation

A finite element approximation of steady flow in a rigid three-dimensional model of the carotid artery bifurcation is presented. A Reynolds number of 640 and a flow division ratio of about 50/50, simulating systolic flow, was used. To limit the CPU- and I/O-times needed for solving the systems of equations, a mesh-generator was developed, which gives full control over the number of elements into which the bifurcation is divided. A mini-supercomputer, based on parallel and vector processing techniques, was used to solve the system of equations. The numerical results of axial and secondary flow compare favorably with those obtained from previously performed laser-Doppler velocity measurements. Also, the influence of the Reynolds number, the flow division ratio, and the bifurcation angle on axial and secondary flow in the carotid sinus were studied in the three-dimensional model. The influence of the interventions is limited to a relatively small variation in the region with reversed axial flow, more or less pronounced C-shaped axial velocity contours, and increasing or decreasing axial velocity maxima.     

Some features of oscillatory flow in a model bifurcation

Oscillatory flow in the lung is studied using an order-of-magnitude analysis and flow visualization experiments in a single bifurcation with lung-like geometry. The results are used to obtain a classification scheme that identifies three major flow regimes, distinguished on the basis of whether the flow is dominated by unsteadiness, viscous effects, or the effects of convective acceleration. The unsteady regime is found to exist for values of a dimensionless stroke length (L/a, i.e., stroke volume/local cross-sectional area) less than or equal to 3 and for values of a dimensionless frequency (alpha 2 = alpha 2 omega/nu, where alpha is airway radius, omega the oscillatory frequency, and nu the kinematic viscosity) less than or equal to 10 in basic agreement with previous studies. The viscous regime is found when alpha 2(L/a)(a/R)1/2 less than 10 and alpha 2 less than 10 where R is the local radius of curvature in the bifurcation; the convective regime is found when alpha 2(L/a)(a/R)1/2 greater than 10 and L/a greater than 3. This same approach yields scaling laws for the magnitude of secondary flow velocities and shows that the ratio of secondary-to-axial velocity is small everywhere outside of the convective regime where it scales with (a/R)1/2. Comparison of these results to related simple flows shows that many of the features observed can be attributed to the effects of curvature, suggesting that the influence of the flow divider and of area change may be of lesser importance than previously thought.     

Comparison of axisymmetric and three-dimensional models for gas uptake in a single bifurcation during steady expiration

Reactive gas uptake is predicted and compared in a single bifurcation at steady expiratory flow in terms of Sherwood number using an axisymmetric single-path model (ASPM) and a three-dimensional computational fluid dynamics model (CFDM). ASPM is validated in a two-generation geometry by comparing the average gas-phase mass transfer coefficients with the experimental values. ASPM predicted mass transfer coefficients within 20% of the experimental values. The flow and concentration variables in the ASPM were solved using Galerkin finite element method and in the CFDM using commercial finite element software FIDAP. The simulations were performed for reactive gas flowing at Reynolds numbers ranging from 60 to 350 in both symmetric bifurcation for three bifurcation angles, 30 deg, 70 deg, and 90 deg, and in an asymmetric bifurcation. The numerical models compared with each other qualitatively but quantitatively they were within 0.4-8% due to nonfully developed flow in the parent branch predicted by the CFDM. The radially averaged concentration variation along the axial location matched qualitatively between the CFDM and ASPM but quantitatively they were within 32% due to differences in the flow field. ASPM predictions compared well with the CFDM predictions for an asymmetric bifurcation. These results validate the simplified ASPM and the complex CFDM. ASPM predicts higher Sherwood number with a flat velocity inlet profile compared to a parabolic inlet velocity profile. Sherwood number increases with the inlet average velocity, wall mass transfer coefficient, and bifurcation angle since the boundary layer grows slower in the parent and daughter branches.     

Numerical study on the effect of secondary flow in the human aorta on local shear stresses in abdominal aortic branches

Flow in the aortic arch is characterized primarily by the presence of a strong secondary flow superimposed over the axial flow, skewed axial velocity profiles and diastolic flow reversals. A significant amount of helical flow has also been observed in the descending aorta of humans and in models. In this study a computational model of the abdominal aorta complete with two sets of outflow arteries was adapted for three-dimensional steady flow simulations. The flow through the model was predicted using the Navier-Stokes equations to study the effect that a rotational component of flow has on the general flow dynamics in this vascular segment. The helical velocity profile introduced at the inlet was developed from magnetic resonance velocity mappings taken from a plane transaxial to the aortic arch. Results showed that flow division ratios increased in the first set of branches and decreased in the second set with the addition of rotational flow. Shear stress varied in magnitude with the addition of rotational flow, but the shear stress distribution did not change. No regions of flow separation were observed in the iliac arteries for either case. Helical flow may have a stabilizing effect on the flow patterns in branches in general, as evidenced by the decreased difference in shear stress between the inner and outer walls in the iliac arteries.     

The effect of unsteadiness on the flow through stenoses and bifurcations

This paper is concerned with the influence of a stenosis or a bifurcation on the flow through a tube. In particular the effect of unsteadiness is investigated using simple pulsatile and physiological type flows (Fig. 1). The experimental investigations reported herein are concerned with velocity measurements and flow visualizations. (see formula in text) These measurements, performed in a 60 degrees bifurcation, have permitted the reconstruction of the three-dimensional velocity profiles. The importance of the secondary flow in the branching is analyzed for various values of the flow parameters. Results of tests show a strong influence of unsteadiness on flow characteristics and then on hemodynamic factors. One conclusion is the following: if hemodynamic factors play an important role in the problems of atherosclerosis, then, for macrocirculation studies, it is necessary to take into account unsteadiness and, in particular, the actual shape of the flow-time forcing function.     

Blood flow multiscale phenomena

The cardiovascular disease is one of most frequent cause deaths in modern society. The objective of this work is analyse the effect of dynamic vascular geometry (curvature, torsion, bifurcation) and pulsatile blood nature on secondary flow, wall shear stress and platelet deposition. The problem was examined as multi-scale physical phenomena using perturbation analysis and numerical modelling. The secondary flow determined as influence pulsatile pressure, vascular tube time-dependent bending and torsion on the main axial flow. Bifurcation and branching phenomena are analysed experimentally through, blood-like fluid pulsatile flow across elastic rubber-like Y-model model. The problem complex geometry near branching in platelet deposit modelling is resolved numerically as Falker-Skan flow.     

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.     

An LDA and flow visualization study of pulsatile flow in an aortic bifurcation model

The structure of pulsatile flow in a rigid aortic bifurcation model was studied by means of a flow visualization technique and three-dimensional laser-Doppler anemometry. The model was made of glass, having the same shape as that of the average human aortic bifurcation. It was installed into a mock circulatory loop that generated physiological pulsatile flow. Flow separation was observed during accelerated and decelerated flow periods. Double helical flow existed inside the flow separation in the early accelerated flow period. In the decelerated flow period, disturbed flow appeared behind the separation zone. Flow was strongly disturbed during the back flow period, and then was gradually stabilized in the forward flow period. The flow separation and the disturbances released from the flow separation zone greatly influenced near-wall velocities along the lateral wall. The wave form of the near-wall velocity in the flow separation zone was much different from that observed in the aortic portion and behind the separation zone; for example, the magnitude of the negative peak velocity in the direction of the tube axis was larger than that of the positive one, and mean velocity in a cycle was very low. This abnormal phasic change of the near-wall velocity may be associated with atherogenesis. The three-dimensional velocity measurement is very useful for the detailed analysis of near-wall velocity patterns.     

Bubble splitting in bifurcating tubes: a model study of cardiovascular gas emboli transport.

The transport of long gas bubbles, suspended in liquid, through symmetric bifurcations, is investigated experimentally and theoretically as a model of cardiovascular gas bubble transport in air embolism and gas embolotherapy. The relevant dimensionless parameters in the models match the corresponding values for arteries and arterioles. The effects of roll angle (the angle the plane of the bifurcation makes with the horizontal), capillary number (a dimensionless indicator of flow), and bubble volume (or length) on the splitting of bubbles as they pass through the bifurcation are examined. Splitting is observed to be more homogenous at higher capillary numbers and lower roll angles. It is shown that, at nonzero roll angles, there is a critical value of the capillary number below which the bubbles do not split and are transported entirely into the upper branch. The value of the critical capillary number increases with roll angle and parent tube diameter. A unique bubble motion is observed at the critical capillary number and for slightly slower flows: the bubble begins to split, the meniscus in the lower branch then moves backward, and finally the entire bubble enters the upper branch. These findings suggest that, in large vessels, emboli tend to be transported upward unless flow is unusually strong but that a more homogeneous distribution of emboli occurs in smaller vessels. This corresponds to previous observations that air emboli tend to lodge in the upper regions of the lungs and suggests that relatively uniform infarction of tumors by gas embolotherapy may be possible.     

Secondary velocity fields in the conducting airways of the human lung

An understanding of flow and dispersion in the human respiratory airways is necessary to assess the toxicological impact of inhaled particulate matter as well as to optimize the design of inhalable pharmaceutical aerosols and their delivery systems. Secondary flows affect dispersion in the lung by mixing solute in the lumen cross section. The goal of this study is to measure and interpret these secondary velocity fields using in vitro lung models. Particle image velocimetry experiments were conducted in a three-generational, anatomically accurate model of the conducting region of the lung. Inspiration and expiration flows were examined under steady and oscillatory flow conditions. Results illustrate secondary flow fields as a function of flow direction, Reynolds number, axial location with respect to the bifurcation junction, generation, branch, phase in the oscillatory cycle, and Womersley number. Critical Dean number for the formation of secondary vortices in the airways, as well as the strength and development length of secondary flow, is characterized. The normalized secondary velocity magnitude was similar on inspiration and expiration and did not vary appreciably with generation or branch. Oscillatory flow fields were not significantly different from corresponding steady flow fields up to a Womersley number of 1 and no instabilities related to shear were detected on flow reversal. These observations were qualitatively interpreted with respect to the simple streaming, augmented dispersion, and steady streaming convective dispersion mechanisms.     

Frequency dependence of dynamic curvature effects on flow through coronary arteries

The flow through a curved tube model of a coronary artery was investigated computationally to determine the importance of time-varying curvature on flow patterns that have been associated with the development of atherosclerosis. The entry to the tube was fixed while the radius of curvature varied sinusoidally in time at a frequency of 1 or 5 Hz. Angiographic data from other studies suggest that the radius of curvature waveform contains significant spectral content up to 6 Hz. The overall flow patterns were similar to those observed in stationary curved tubes; velocity profile skewed toward the outer wall, secondary flow patterns, etc. The effects of time-varying curvature on the changes in wall shear rate were expressed by normalizing the wall shear rate amplitude with the shear rate calculated at the static mean radius of curvature. It was found that the wall shear rate varied as much as 94 percent of the mean wall shear rate at the mid wall of curvature for a mean curvature ratio of 0.08 and a 50 percent change in radius of curvature. The effects of 5 Hz deformation were not well predicted by a quasi-static approach. The maximum values of the normalized inner wall shear rate amplitude were found to scale well with a dimensionless parameter equivalent to the product of the mean curvature ratio (delta), normalized change in radius of curvature (epsilon), and a Womersley parameter (alpha). This parameter was less successful at predicting the amplitudes elsewhere in the tube, thus additional studies are necessary. The mean wall shear rate was well predicted with a static geometry. These results indicate that dynamic curvature plays an important role in determining the inner wall shear rates in coronary arteries that are subjected to deformation levels of epsilon delta alpha > 0.05. The effects were not always predictable with a quasi-static approach. These results provide guidelines for constructing more realistic models of coronary artery flow for atherogenesis research.     

Hemodynamics in the carotid artery bifurcation: a comparison between numerical simulations and in vitro MRI measurements

The presence of atherosclerotic plaques has been shown to be closely related to the vessel geometry. Studies on postmortem human arteries and on the experimental animal show positive correlation between the presence of plaque thickness and low shear stress, departure of unidirectional flow and regions of flow separation and recirculation. Numerical simulations of arterial blood flow and direct blood flow velocity measurements by magnetic resonance imaging (MRI) are two approaches for the assessment of arterial blood flow patterns. In order to verify that both approaches give equivalent results magnetic resonance velocity data measured in a compliant anatomical carotid bifurcation model were compared to the results of numerical simulations performed for a corresponding computational vessel model. Cross sectional axial velocity profiles were calculated and measured for the midsinus and endsinus internal carotid artery. At both locations a skewed velocity profile with slow velocities at the outer vessel wall, medium velocities at the side walls and high velocities at the flow divider (inner) wall were observed. Qualitative comparison of the axial velocity patterns revealed no significant differences between simulations and in vitro measurements. Even quantitative differences such as for axial peak flow velocities were less than 10%. Secondary flow patterns revealed some minor differences concerning the form of the vortices but maximum circumferential velocities were in the same range for both methods.     

Effects of gas properties and waveform asymmetry on gas transport in a branching tube network

Local gas transport coefficients, quantifying longitudinal dispersion through a symmetrical constant-diameter tube network, have been measured during oscillation with both symmetrical and nonsymmetrical waveforms. Experiments were carried out over a range of conditions that would prevail in the central to lower airways during high-frequency ventilation at moderate frequency (5 Hz) and tidal volume (15-80 ml). Gas transport coefficients resulting from oscillation of three different resident-trace gas pairs were measured using a new analytic technique. This technique allowed rapid determination of the transport coefficient distribution along the entire network. Results demonstrate a small but significant influence attributable to changes in gas properties that is similar to that found in a straight tube and indicate that augmented dispersion is an important mechanism of axial transport. Gas transport coefficients were found to be unaffected by changes in flow waveform symmetry, suggesting that previously reported improvements in gas exchange associated with decreasing inspiratory to expiratory time ratios are not due to a change in local conditions such as asymmetry in the velocity profile.