Steady Flow Visualization in a Rigid Model of the Aortic Bifurcation: Application to Atherosclerosis

Hemodynamics have long been implicated in atherogenesis. The studiesreported here seek to explain the mechanisms for the formation ofatherosclerotic plaque in an aortic bifurcation. Flow studies were made ina model constructed from plexiglass to represent an aortic bifurcation. Under steady flow conditions at inflow Reynolds numbers of 80–1250,the streamline flow patterns and the boundary layer separation zones wereinvestigated in relation to the location of atherosclerotic plaques clinicallyfound at regions in the human aortic bifurcation. The streamline flowswere visualized by a slow injection of dye over the cross section of the tubeentrance and along the tube walls. The studies revealed a complex flowfield where secondary flows, induced by the centrifugal and viscous forces,cause the fluid to move towards the inner walls of the aortic bifurcation. The effect was more clearly seen with increasing Reynolds number. Boundary layer separation zones were observed to occur at the outercorners of the branching. The nature of the separation zone formed wasfound to be dependent on Reynolds number. The residence time of fluidparticles within such a separation zone was estimated by measuring thewashout time of a bolus of dye injected at strategic locations along the tubewalls. The residence time was found to decrease exponentially withincreasing Reynolds number. These observations provide strong support forthe role of flow separation in the accumulation of LDL and plateletaggregation within the aortic bifurcation.     

Low Reynolds number equi-bifurcation flow in a two-dimensional channel

The fluid flow in some physiological vessels such as the blood flow in blood vessels and the air flow through bronchi and bronchioles in the lungs undergoes a large number of bifurcations. The understanding of the bifurcation flow is of importance for a better comprehension of its effect in the blood and the air circulatory systems of the living body. The Reynolds number of flow in large blood vessels and bronchi is high and fluid inertia plays a dominant role in the bifurcation flow in such vessels. In small caliber blood vessels such as arterioles and capillaries, and bronchioles, the Reynolds number of flow is quite low and the effect of fluid inertia is negligible compared to the pressure and shear forces. In order to have a quantitative understanding of the bifurcation flow at low Reynolds numbers, the low Reynolds number equi-bifurcation flow in a two-dimensional channel at zero bifurcation angle is studied based on the Stokes approximation. The solution of the problem is posed as an infinite series, where the truncated version is used in numerical calculations. The results of this analysis is discussed in connection with the bifurcation flow of blood in small caliber blood vessels and that of the air in bronchioles in the lung.     

Branch angle and flow into a symmetric bifurcation

Arterial branches are found to be a major site for formation of arterial plaque. In this study, we investigate the role of the bifurcation angle on the flow into a symmetric bifurcation. Specially, how the changes in the bifurcation angle influences the distribution of axial wall shear in the bifurcation model. The flow in a range of branch opening half-angle of pi/25< or =theta< or =pi/4 are numerically simulated. The flow in the above models is calculated for the inlet flow Reynolds numbers of 250, 500, 1000, and 2000. It is found that at higher values of the opening angle of the bifurcation, the possibility and severity of flow separation at the appropriate wall location increases.     

LDA measurements and numerical prediction of pulsatile laminar flow in a plane 90-degree bifurcation

An experimental and numerical study of pulsatile laminar flow in a plane 90-degree bifurcation is presented. Detailed LDA velocity measurements of the oscillatory flow field have been carried out. The numerical predictions, which are based on an iterative, finite-difference numerical procedure using primitive dependent variables, are in good agreement with the measurements. The results show that one separation zone is established near the bottom wall of the main duct and another near the upstream wall of the branch. The location and size of the separation zones vary within the cycle and are influenced by the Reynolds number, the flow rate ratio, and the Stokes number.     

Finite element analysis of two dimensional steady flow in model arterial bifurcation

The purpose of the investigation reported in this paper is to determine theoretically the fluid dynamic field in models of common iliac arterial bifurcation and to identify the flow features which might influence the predominant occurrence of atherosclerotic lesions at such sites. This has been accomplished by numerically simulating fluid flow through 90 degrees symmetric bifurcations with branch-to-trunk area ratios of 0.8-1.414 and for Reynolds numbers ranging from 100 to 400. The analysis predicts flow reversal along the outer wall in models with area ratios over unity for high Reynolds number range, while no flow reversal occurred in models with area ratio below unity; a low shear zone along the outer wall and high shear stresses at the divider lip. Adverse pressure gradients are observed along the outer wall downstream of the corner point, the magnitudes increased with Reynolds number for a given branch to area ratio. Biological implication of the results is discussed with specific reference to the sites of atherosclerotic lesions found in man for these geometries.     

Computational Analysis of Micron-particle Deposition in a Human Triple Bifurcation Airway Model

Steady laminar axisymmetric inhalation flow and wall deposition of micron-size particles in representative triple bifurcation airways have been simulated using a commercial finite-volume code with user-enhanced programs. Assuming spherical non-interacting particles (3 μm≤ d p ≤7 μm), various inlet Reynolds numbers (Re=500-2000) and Stokes numbers (St=0.02-0.23) were considered. The resulting particle deposition patterns were analyzed and then summarized in terms of deposition efficiencies, i.e. DE=DE(Re,St) Surprisingly high DE-values occur at relatively low Reynolds numbers (e.g., Re=500 ) in the third bifurcation. The quantitative results are of interest to researchers either conducting health risk assessment studies for inhaled particulate pollutants or analyzing drug aerosol inhalation and deposition at desired lung target sites.     

Flow in t-bifurcations: Effect of the sharpness of the flow divider

The local geometry of a bifurcation has been hypothesized to be a potential geometrical risk factor for the development of atherosclerosis. While flow division and branch area ratios clearly affect the flow field, the importance of the flow divider shape is not as clear. A fast spectral element computational fluid mechanics (CFD) solver was used to simulate flow through 90 ° T-bifurcations with three different flow divider shapes. Other factors, such as flow partition, area ratio, and bifurcation angle, were kept constant. A Reynolds number range of 15 to 350 was studied to bracket experimental results in the literature. The variation in the sharpness of the corners had a dramatic effect on both the flow field and wall shear stress distribution in the side branch, but little effect on the flow in the main tube. The magnitude of reverse velocities and wall shear stress in the side branch increased linearly over a physiological range of Reynolds number and corner shape. This paper verifies the accuracy and usefulness of spectral element CFD in studying three-dimensional hemodynamics.     

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.     

The effect of asymmetry in abdominal aortic aneurysms under physiologically realistic pulsatile flow conditions

In the abdominal segment of the human aorta under a patient's average resting conditions, pulsatile blood flow exhibits complex laminar patterns with secondary flows induced by adjacent branches and irregular vessel geometries. The flow dynamics becomes more complex when there is a pathological condition that causes changes in the normal structural composition of the vessel wall, for example, in the presence of an aneurysm. This work examines the hemodynamics of pulsatile blood flow in hypothetical three-dimensional models of abdominal aortic aneurysms (AAAs). Numerical predictions of blood flow patterns and hemodynamic stresses in AAAs are performed in single-aneurysm, asymmetric, rigid wall models using the finite element method. We characterize pulsatile flow dynamics in AAAs for average resting conditions by means of identifying regions of disturbed flow and quantifying the disturbance by evaluating flow-induced stresses at the aneurysm wall, specifically wall pressure and wall shear stress. Physiologically realistic abdominal aortic blood flow is simulated under pulsatile conditions for the range of time-average Reynolds numbers 50 < or = Rem < or = 300, corresponding to a range of peak Reynolds numbers 262.5 < or = Repeak < or = 1575. The vortex dynamics induced by pulsatile flow in AAAs is depicted by a sequence of four different flow phases in one period of the cardiac pulse. Peak wall shear stress and peak wall pressure are reported as a function of the time-average Reynolds number and aneurysm asymmetry. The effect of asymmetry in hypothetically shaped AAAs is to increase the maximum wall shear stress at peak flow and to induce the appearance of secondary flows in late diastole.     

Numerical investigation of the non-Newtonian blood flow in a bifurcation model with a non-planar branch

The non-Newtonian fluid flow in a bifurcation model with a non-planar daughter branch is investigated by using finite element method to solve the three-dimensional Navier–Stokes equations coupled with a non-Newtonian constitutive model, in which the shear thinning behavior of the blood fluid is incorporated by the Carreau–Yasuda model. The objective of this study is to investigate the influence of the non-Newtonian property of fluid as well as of curvature and out-of-plane geometry in the non-planar daughter vessel on wall shear stress (WSS) and flow phenomena. In the non-planar daughter vessel, the flows are typified by the skewing of the velocity profile towards the outer wall, creating a relatively low WSS at the inner wall. In the downstream of the bifurcation, the velocity profiles are shifted towards the flow divider. The low WSS is found at the inner walls of the curvature and the lateral walls of the bifurcation. Secondary flow patterns that swirl fluid from the inner wall of curvature to the outer wall in the middle of the vessel are also well documented for the curved and bifurcating vessels. The numerical results for the non-Newtonian fluid and the Newtonian fluid with original Reynolds number and the corresponding rescaled Reynolds number are presented. Significant difference between the non-Newtonian flow and the Newtonian flow is revealed; however, reasonable agreement between the non-Newtonian flow and the rescaled Newtonian flow is found. Results of this study support the view that the non-planarity of blood vessels and the non-Newtonian properties of blood are an important factor in hemodynamics and may play a significant role in vascular biology and pathophysiology.     

Computational analysis of micron-particle deposition in a human triple bifurcation airway model

Steady laminar axisymmetric inhalation flow and wall deposition of micron-size particles in representative triple bifurcation airways have been simulated using a commercial finite-volume code with user-enhanced programs. Assuming spherical non-interacting particles (3 microm相似文献   

Three-dimensional steady flow through a bifurcation

Steady flow of an incompressible, Newtonian fluid through a symmetric bifurcated rigid channel was numerically analyzed by solving the three-dimensional Navier-Stokes equations. The upstream Reynolds number ranged from 100 to 1500. The bifurcation was symmetrical with a branch angle of 60 deg and the area ratio of the daughter to the mother vessel was 2.0. The numerical procedure utilized a coordinate transformation and a control volume approach to discretize the equations to finite difference form and incorporated the SIMPLE algorithm in performing the calculation. The predicted velocity pattern was in qualitative agreement with experimental measurements available in the literature. The results also showed the effect of secondary flow which can not be predicted using previous two-dimensional simulations. A region of reversed flow was observed near the outer wall of the branch except for the case of the lowest Reynolds number. Particle trajectory was examined and it was found that no fluid particles remained within the recirculation zone. The shear stress was calculated on both the inner and the outer wall of the branch. The largest wall shear stress, located in the vicinity of the apex of the branch, was of the same order of magnitude as the level that can cause damage to the vessel wall as reported in a recent study.     

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.     

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.     

A numerical study of the shape of the surface separating flow into branches in microvascular bifurcations.

The shape of the separating surface formed by the streamlines entering the branches of microvascular bifurcations plays a major role in determining the distribution of red blood cells and other blood constituents downstream from the bifurcation. Using the finite element method, we determined the shape of the surface through numerical solution of three dimensional Navier-Stokes equations for fluid flow at low Reynolds numbers in a T-type bifurcation of circular tubes. Calculations were done for a wide range of daughter branch to parent vessel diameter ratios and flow ratios. The effect of Reynolds number was also studied. Our numerical results are in good agreement with previously reported experimental data of Rong and Carr (Microvascular Research, Vol. 39, pp. 186-202, 1990). The numerical results of this study will be used to predict the concentration of blood constituents downstream from microvascular bifurcations providing that the inlet concentration profile is known.     

Development of a general method for designing microvascular networks using distribution of wall shear stress

In the present study, theoretical formulations for calculation of optimal bifurcation angle and relationship between the diameters of mother and daughter vessels using the power law model for non-Newtonian fluids are developed. The method is based on the distribution of wall shear stress in the mother and daughter vessels. Also, the effect of distribution of wall shear stress on the minimization of energy loss and flow resistance is considered. It is shown that constant wall shear stress in the mother and daughter vessels provides the minimum flow resistance and energy loss of biological flows. Moreover, the effects of different wall shear stresses in the mother and daughter branches, different lengths of daughter branches in the asymmetric bifurcations and non-Newtonian effect of biological fluid flows on the bifurcation angle and the relationship between the diameters of mother and daughter branches are considered. Using numerical simulations for non-Newtonian models such as power law and Carreau models, the effects of optimal bifurcation angle on the pressure drop and flow resistance of blood flow in the symmetric bifurcation are investigated. Numerical simulations show that optimal bifurcation angle decreases the pressure drop and flow resistance especially for bifurcations at large Reynolds number.     

Validating CFD predictions of respiratory aerosol deposition: effects of upstream transition and turbulence

A number of computational fluid dynamics (CFD) studies have explored local deposition patterns of inhaled aerosols in the respiratory tract. These studies have highlighted the effects of multiple physiologic, geometric, and particle characteristics on deposition. However, very few studies have reported local or sub-branch quantitative comparisons to in vitro particle deposition data. The objective of this study is to numerically investigate the effects of transition and turbulence on highly localized particle deposition in a respiratory double bifurcation model in order to quantitatively validate CFD results. To perform the validations, local comparisons have been made to a specific in vitro case study of 10 microm particles depositing in a model of respiratory generations G3-G5. To achieve this objective, two geometric cases have been considered. The first case includes only the double bifurcation model. The second case includes a portion of the experimental particle delivery geometry, where transitional flow is expected. To evaluate the effectiveness of two-equation turbulence models in this system, the flow field solutions have been computed using laminar, standard k-omega, and low Reynolds number (LRN) k-omega approximations. Results indicate that even though the Reynolds number remained below the critical limit required for full turbulence, transition and turbulence have a significant impact on the flow field and local particle deposition patterns. For the experimental case considered, turbulence impacted the local deposition of 10 microm particles primarily by influencing the initial velocity and particle profiles. As such, both the laminar and LRN k-omega flow models provided good local quantitative matches to the in vitro deposition data, provided that the correct initial particle profile was specified. Implications of this study include the need for local quantitative validations of particle deposition results, the importance of correct inlet conditions, and the need to consider upstream effects in experimental and computational studies of the respiratory tract. Applications of these results to realistic respiratory geometries will require consideration on upstream flow conditions in the lung, transient flow, and intermittent turbulent structures.     

Pulsatile blood flow and oxygen transport past a circular cylinder

The fundamental study of blood flow past a circular cylinder filled with an oxygen source is investigated as a building block for an artificial lung. The Casson constitutive equation is used to describe the shear-thinning and yield stress properties of blood. The presence of hemoglobin is also considered. Far from the cylinder, a pulsatile blood flow in the x direction is prescribed, represented by a time periodic (sinusoidal) component superimposed on a steady velocity. The dimensionless parameters of interest for the characterization of the flow and transport are the steady Reynolds number (Re), Womersley parameter (alpha), pulsation amplitude (A), and the Schmidt number (Sc). The Hill equation is used to describe the saturation curve of hemoglobin with oxygen. Two different feed-gas mixtures were considered: pure O(2) and air. The flow and concentration fields were computed for Re=5, 10, and 40, 0< or =A< or =0.75, alpha=0.25, 0.4, and Schmidt number, Sc=1000. The Casson fluid properties result in reduced recirculations (when present) downstream of the cylinder as compared to a Newtonian fluid. These vortices oscillate in size and strength as A and alpha are varied. Hemoglobin enhances mass transport and is especially important for an air feed which is dominated by oxyhemoglobin dispersion near the cylinder. For a pure O(2) feed, oxygen transport in the plasma dominates near the cylinder. Maximum oxygen transport is achieved by operating near steady flow (small A) for both feed-gas mixtures. The time averaged Sherwood number, Sh, is found to be largely influenced by the steady Reynolds number, increasing as Re increases and decreasing with A. Little change is observed with varying alpha for the ranges investigated. The effect of pulsatility on Sh is greater at larger Re. Increasing Re aids transport, but yields a higher cylinder drag force and shear stresses on the cylinder surface which are potentially undesirable.     

Numerical simulations of the pulsating flow of cerebrospinal fluid flow in the cervical spinal canal of a Chiari patient

The flow of cerebrospinal fluid (CSF) in a patient-specific model of the subarachnoid space in a Chiari I patient was investigated using numerical simulations. The pulsating CSF flow was modeled using a time-varying velocity pulse based on peak velocity measurements (diastole and systole) derived from a selection of patients with Chiari I malformation. The present study introduces the general definition of the Reynolds number to provide a measure of CSF flow instability to give an estimate of the possibility of turbulence occurring in CSF flow. This was motivated by the fact that the combination of pulsating flow and the geometric complexity of the spinal canal may result in local Reynolds numbers that are significantly higher than the commonly used global measure such that flow instabilities may develop into turbulent flow in these regions. The local Reynolds number was used in combination with derived statistics to characterize the flow. The results revealed the existence of both local unstable regions and local regions with velocity fluctuations similar in magnitude to what is observed in fully turbulent flows. The results also indicated that the fluctuations were not self-sustained turbulence, but rather flow instabilities that may develop into turbulence. The case considered was therefore believed to represent a CSF flow close to transition.     

Computational analysis of flow structure and particle deposition in a single asthmatic human airway bifurcation

This paper aims to improve current understanding of flow structure and particle deposition in asthmatic human airways. A single, symmetric airway bifurcation, corresponding to generations 10–11 of Weibel’s model, is investigated through validated numerical simulations. The parent airway segment is modelled as a smooth circular tube. The child segments are considered asthmatic and their cross-section is modelled as a constricted tube with sinusoidal folds uniformly distributed along the circumference. The flow structure and particle deposition pattern for normal (i.e., healthy) and asthmatic airway bifurcations are compared and discussed. The numerical results reveal that the secondary flow in the asthmatic airway bifurcation is much stronger than in the healthy one, resulting in higher particle deposition. The effects of size of the lumen area and number of folds on particle deposition and pressure drop are also investigated. It is found that particle deposition efficiency is significantly affected by lumen area of the asthmatic segment (the smaller the lumen area, the higher the particle deposition efficiency). The effect of number of folds is small. Particle deposition efficiency also increases with Reynolds number. The pressure drop in the asthmatic airway bifurcation depends mainly on size of the lumen area. The effect of number of folds becomes important for strongly collapsed airways.