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 mechanism for erythrocyte-mediated elevation of apparent viscosity by leukocytes in vivo without adhesion to the endothelium

In spite of the relatively small number of leukocytes in the circulation, they have a significant influence on the perfusion of such organs as skeletal muscle or kidney. However, the underlying mechanisms are incompletely understood. In the current study a combined in vivo and computational approach is presented in which the interaction of individual freely flowing leukocytes with erythrocytes and its effect on apparent blood viscosity are explored. The skeletal muscle microcirculation was perfused with different cell suspensions with and without leukocytes or erythrocytes. We examined a three-dimensional numerical model of low Reynolds number flow in a capillary with a train of erythrocytes (small spheres) in off-axis positions and single larger leukocytes in axisymmetric positions. The results indicate that in order to match the slower axial velocity of leukocytes in capillaries, erythrocytes need to position themselves into an off-axis position in the capillary. In such off-axis positions at constant mean capillary velocity, erythrocyte axial velocity matches on average the axial velocity of the leukocytes, but the apparent viscosity is elevated, in agreement with the whole organ perfusion observations. Thus, leukocytes influence the whole organ resistance in skeletal muscle to a significant degree only in the presence of erythrocytes.     

Velocity profiles of blood flow in microvessels measured by ten channels' dual-sensor method

Ten channels' dual-sensor method newly developed was applied to the measurement of velocity profiles in arterioles and venules in the rat mesentery. In some experiments, red blood cells (RBC), labelled in vitro with a fluorescein isothiocyanate (FITC), were injected to compare the velocity profiles obtained by the RBC visualization technique with those measured by the dual-sensor method. It was found that the velocity profile of the FITC labelled RBC in straight microvessels was blunt as compared to a parabola. The centerline velocity measured by the dual-sensor method was smaller than that of the FITC labelled RBC by about 20%. The velocity profiles were also measured at the curved arterioles and venules as well as at the bifurcation and the confluence. It was found that the velocities were higher along the inner wall at the curved portion and along the outer wall at the bifurcation of arterioles.     

Augmentation of axial dispersion by intermittent oscillatory flow.

The efficiency of axial gas dispersion during ventilation with high-frequency oscillation (HFO) is improved by manipulating the oscillatory flow waveform such that intermittent oscillatory flow occurs. We therefore measured the velocity profiles and effective axial gas diffusivity during intermittent oscillatory flow in a straight tube to verify the intermittency augmentation effect on axial gas transfer. The effective diffusivity was dependent on the flow patterns and significantly increased with an increase in the duration of the stationary phase. It was also found that the ratio of effective diffusivity to molecular diffusivity is two times greater than that in sinusoidal oscillatory flow. Moreover, turbulence during deceleration or at the beginning of the stationary phase further augments axial dispersion, with the effective diffusivity being over three times as large, thereby proving that the use of intermittent oscillatory flow effectively augments axial dispersion for ventilation with HFO.     

Pressure flow patterns in a cylinder with reabsorbing walls

The problem of fluid motion in renal tubules, in contrast to ordinary flow through cylinders with impermeable walls, is complicated by the existence of radial velocities generated by reabsorption processes. As a first approach to this problem, the Navier Stokes equations for axially symmetric, slow flow in an infinite cylinder whose walls reabsorb fluid are integrated. If the rate of reabsorption is constant, the solutions resemble the conventional Poiseuille flow, i.e., the longitudinal velocity profile is parabolic. In addition the drop in mean pressure is proportional to the mean axial flow, the length of tube between reference points, and inversely proportional to the fourth power of the radius. If the rate of reabsorption is a linear function of the distance from the origin, the presence of an additive term alters these relations. If, for example, the gradient in reabsorption is positive, the axial velocity profile tends to flatten and when the gradient is sufficiently large, the maximum velocity moves from the center of the stream toward the periphery, leaving a relative minimum at the center. In passing from the center of the tube to the walls, the radial velocity passes through a miximum, regardless of the reabsorption properties of the wall.     

Effects of curved inlet tubes on air flow and particle deposition in bifurcating lung models

In vivo bifurcating airways are complex and the airway segments leading to the bifurcations are not always straight, but curved to various degrees. How do such curved inlet tubes influence the motion as well as local deposition and hence the biological responses of inhaled particulate matter in lung airways? In this paper steady laminar dilute suspension flows of micron-particles are simulated in realistic double bifurcations with curved inlet tubes, i.e., 0 degrees < or =theta< or =90 degrees, using a commercial finite-volume code with user-enhanced programs. The resulting air-flow patterns as well as particle transport and wall depositions were analyzed for different flow inlet conditions, i.e., uniform and parabolic velocity profiles, and geometric configurations. The curved inlet segments have quite pronounced effects on air-flow, particle motion and wall deposition in the downstream bifurcating airways. In contrast to straight double bifurcations, those with bent parent tubes also exhibit irregular variations in particle deposition efficiencies as a function of Stokes number and Reynolds number. There are fewer particles deposited at mildly curved inlet segments, but the particle deposition efficiencies at the downstream sequential bifurcations vary much when compared to those with straight inlets. Under certain flow conditions in sharply curved lung airways, relatively high, localized particle depositions may take place. The findings provide necessary information for toxicologic or therapeutic impact assessments and for global lung dosimetry models of inhaled particulate matter.     

Evaluation of magnetic resonance velocimetry for steady flow

Whole body magnetic resonance (MR) imaging has recently become an important diagnostic tool for cardiovascular diseases. The technique of magnetic resonance phase velocity encoding allows the quantitative measurement of velocity for an arbitrary component direction. A study was initiated to determine the ability and accuracy of MR velocimetry to measure a wide range of flow conditions including flow separation, three-dimensional secondary flow, high velocity gradients, and turbulence. A steady flow system pumped water doped with manganese chloride through a variety of test sections. Images were produced using gradient echo sequences on test sections including a straight tube, a curved tube, a smoothly converging-diverging nozzle, and an orifice. Magnetic resonance measurements of laminar and turbulent flows were depicted as cross-sectional velocity profiles. MR velocity measurements revealed such flow behavior as spatially varying velocity, recirculation and secondary flows over a wide range of conditions. Comparisons made with published experimental laser Doppler anemometry measurements and theoretical calculations for similar flow conditions revealed excellent accuracy and precision levels. The successful measurement of velocity profiles for a variety of flow conditions and geometries indicate that magnetic resonance imaging is an accurate, non-contacting velocimeter.     

Pulsatile blood flow effects on temperature distribution and heat transfer in rigid vessels

The effect of blood velocity pulsations on bioheat transfer is studied. A simple model of a straight rigid blood vessel with unsteady periodic flow is considered. A numerical solution that considers the fully coupled Navier-Stokes and energy equations is used for the simulations. The influence of the pulsation rate on the temperature distribution and energy transport is studied for four typical vessel sizes: aorta, large arteries, terminal arterial branches, and arterioles. The results show that: the pulsating axial velocity produces a pulsating temperature distribution; reversal of flow occurs in the aorta and in large vessels, which produces significant time variation in the temperature profile. Change of the pulsation rate yields a change of the energy transport between the vessel wall and fluid for the large vessels. For the thermally important terminal arteries (0.04-1 mm), velocity pulsations have a small influence on temperature distribution and on the energy transport out of the vessels (8 percent for the Womersley number corresponding to a normal heart rate). Given that there is a small difference between the time-averaged unsteady heat flux due to a pulsating blood velocity and an assumed nonpulsating blood velocity, it is reasonable to assume a nonpulsating blood velocity for the purposes of estimating bioheat transfer.     

Estimation of viscous dissipation inside an erythrocyte during aspirational entry into a micropipette.

Viscous dissipation inside the erythrocyte during its aspirational entry into a micropipette is analyzed. The motion of the intracellular fluid is approximated by a flow into the micropipette orifice from a half space (the portion of the erythrocyte outside the micropipette). The stream function and intracellular pressure (p) in the half space are obtained as a function of radial and axial positions near the orifice. Solution of the boundary value problem for a uniform stream entering a circular hole gives p = 2 eta HQ/pi R3p, where eta H is the intracellular viscosity, Q is the total discharge, and Rp is the pipette radius. The results indicate that the moving erythrocyte membrane helps to drive the intracellular fluid into the orifice. For normal erythrocytes, p is only approximately 0.5% of the total aspiration pressure (delta P). The contribution of p to delta P, however, may become significant when there is a large increase in eta H due to a markedly elevated intracellular hemoglobin concentration or an alteration of the physical state of hemoglobin.     

Bioartificial kidney. I. Theoretical analysis of convective flow in hollow fiber modules: application to a bioartificial hemofilter

Analytical expressions describing convective flow in a continuous arteriovenous hollow fiber hemofilter were developed. In the lumen of the hollow fiber membrane, existing analytical expressions were applied to describe velocity profiles and pressure. For flow in the shell (the extracapillary space separating the fibers), analytical expressions for the radial and axial velocity profiles and pressure distribution were derived by first finding the stream function. The expressions are based on a similarity solution. Previous analyses of ultrafiltration have either ignored osmotic pressure or assumed constant shell pressure. In this paper, the axial variation in lumen pressure, shell pressure, and osmotic pressure were accounted for. The predicted filtration rates agree well with the experimental results. This flow model is general enough to describe flow in hollow fiber membrane systems employed as bioreactors (e.g., for cell cultures and as bioartificial organs) and as separators (e.g., ultrafiltration and microfiltration) operating in the open-shell mode. The results were applied to determine the design of an optimally functioning bioartificial hemofilter for use ex vivo or in vivo.     

Viscous energy dissipation in a model of the human bronchial tree

The viscous energy dissipation in a two generation model of the human bronchial tree is determined from inspiratory velocity and static pressure data obtained for large Reynolds numbers (10⁴ < Re < 10⁵). This dissipation is found to be an increasing function of both Re and distance downstream from the inlet of the model. The ratio of the dissipation in the model to the energy dissipation in an equivalent straight pipe system is determined. This ratio, Z^*, for the model is compared to values in the literature for lower (laminar) Re. There is more dissipation in the branched model than in a straight pipe (Z^* > 1) and turbulence keeps Z^* at roughly a fixed value for large Reynolds numbers (10⁴ < Re < 10⁵). Z^* values for curved pipes are also compared to the branching system values. It is found that the energy dissipation for the branched model behaves similarly to that in curved pipes.     

The effect of blood viscoelasticity on pulsatile flow in stationary and axially moving tubes

An analytical solution for pulsatile flow of a generalized Maxwell fluid in straight rigid tubes, with and without axial vessel motion, has been used to calculate the effect of blood viscoelasticity on velocity profiles and shear stress in flows representative of those in the large arteries. Measured bulk flow rate Q waveforms were used as starting points in the calculations for the aorta and femoral arteries, from which axial pressure gradient ▿P waves were derived that would reproduce the starting Q waves for viscoelastic flow. The ▿P waves were then used to calculate velocity profiles for both viscoelastic and purely viscous flow. For the coronary artery, published ▿P and axial vessel acceleration waveforms were used in a similar procedure to determine the separate and combined influences of viscoelasticity and vessel motion.Differences in local velocities, comparing viscous flow to viscoelastic flow, were in all cases less than about 2% of the peak local velocity. Differences in peak wall shear stress were less than about 3%.In the coronary artery, wall shear stress differences between viscous and viscoelastic flow were small, regardless of whether axial vessel motion was included. The shape of the wall shear stress waveform and its difference, however, changed dramatically between the stationary and moving vessel cases. The peaks in wall shear stress difference corresponded with large temporal gradients in the combined driving force for the flow.     

A model investigation of the velocity and pressure spectra in vascular murmurs

A physical model consisting of an axisymmetrical jet in a rigid plexiglass pipe was used to study the flow and pressure fluctuations downstream from an aortic stenosis. The fluctuating velocity components, u and v, at several locations in the steady liquid jet were measured using a laser Doppler anemometer system. Simultaneous wall pressure fluctuations were monitored by an array of nine miniature pressure transducers wall mounted in the axial direction. This paper presents the detailed measurements of mean velocity profiles, turbulent intensity distributions and RMS pressure fluctuations. The energy spectra obtained for the pressure fluctuations and the u and v velocity components are compared. Contrary to earlier works, we found that the differences between peak frequencies of the pressure spectra and the characteristic frequencies of the velocity spectra vary with positions downstream from the nozzle. These differences are discussed in light of pseudosound generation by the eddy structures in the stenotic flow field.     

Analysis of blood flow in an out-of-plane CABG model

Coronary artery bypass graft (CABG) is a routine surgical treatment for ischemic and infarcted myocardium. A large number of CABG fail postoperatively because of intimal hyperplasia within months or years. The cause of this failure is thought to be partly related to the flow patterns and shear stresses acting on the endothelial cells. An accurate representation of the flow field and associated wall shear stress (WSS) requires a detailed three-dimensional (3D) model of the CABG. The purpose of this study is to present a detailed analysis of blood flow in a 3D aorto/left CABG, bypassing the occluded left anterior descending coronary (LAD) artery. The analysis takes into account the influence of the out-of-plane geometry of the graft. The finite volume technique was employed to model the 3D blood flow pattern to determine the velocity and WSS distributions. This study presents the flow field distributions of the velocity and WSS at four instances of the cardiac cycle, two in systole and two in diastole. Our results reveal that the CABG geometry has a significant effect on the velocity distribution. The axial velocity profiles at different instances of the cardiac cycle exhibit strong skewing; significant secondary flow and vortex structures are seen in the in-plane velocity patterns. The maximum WSS on the bed of the occluded LAD artery opposite to the graft junction is 14 Pa in middiastole, whereas there is a significantly lower and more uniform distribution of WSS on the bed of the anastomosis. The present results indicate that nonplanarity of the blood vessel along with the inflow conditions has a substantial effect on the fluid mechanics of CABG that contribute to the patency of graft.     

Pressure drop and flow rate measurements in a human aortic bifurcation cast for steady and pulsatile flow

Pressure drop and flow rate measurements in a rigid cast of a human aortic bifurcation under both steady and physiological pulsatile flow conditions are reported. Integral momentum and mechanical energy balances are used to calculate impedance, spatially averaged wall shear stress and viscous dissipation rate from the data. In the daughter branches, steady flow impedance is within 30% of the Poiseuille flow prediction, while pulsatile flow impedance is within a factor of 2 of fully developed, oscillatory, straight tube flow theory (Womersley theory). Estimates of wall shear stress are in accord with measurements obtained from velocity profiles. Mean pressure drop and viscous dissipation rate are elevated in pulsatile flow relative to steady flow at the mean flow rate, and the exponents of their Reynolds number dependence are in accord with available theory.     

Nonisothermal theory of the positive column of an electric discharge in the axial magnetic field

A nonisothermal model of the positive column allowing for electron energy balance is analyzed. The influence of the axial magnetic field on the characteristics of the cylindrical positive column of a low-pressure discharge is investigated in the hydrodynamic approximation. It is shown that the magnetic field affects the plasma density distribution, plasma velocity, and electron energies. The radial dependences of the plasma density, electron energy, and plasma velocity, as well as the azimuthal velocities of electrons and ions, are calculated for helium at different values of the magnetic field strength. It is established that inertia should be taken into account in the equations for the azimuthal motion of electrons and ions. The results obtained in the hydrodynamic approximation differ significantly from those obtained in the framework of the common diffusion model of the positive column in the axial magnetic field. It is shown that the distributions of the plasma density and radial plasma velocity in the greater part of the positive column tend to those obtained in the diffusion approximation at higher values of the axial magnetic field and gas density, although substantial differences remain in the near-wall region.     

Dynamics of blood flow: modeling of Fåhraeus and Fåhraeus–Lindqvist effects using a shear-induced red blood cell migration model

Blood flow in micro capillaries of diameter approximately 15–500 μm is accompanied with a lower tube hematocrit level and lower apparent viscosity as the diameter decreases. These effects are termed the Fåhraeus and Fåhraeus–Lindqvist effects, respectively. Both effects are linked to axial accumulation of red blood cells. In the present investigation, we extend previous works using a shear-induced model for the migration of red blood cells and adopt a model for blood viscosity that accounts for the suspending medium viscosity and local hematocrit level. For fully developed hematocrit profiles (i.e., independent of axial location), the diffusion fluxes due to particle collision frequency and viscosity gradients are of equal magnitude and opposite directions. The ratio of the diffusion coefficients for the two fluxes affects both the Fåhraeus and Fåhraeus–Lindqvist effects and is found related to the capillary diameter and discharge hematocrit using a well-known data-fit correlation for apparent blood viscosity. The velocity and hematocrit profiles were determined numerically as functions of radial coordinate, tube diameter, and discharge hematocrit. The velocity profile determined numerically is consistent with the derived analytical expression and the results are in good agreement with published numerical results and experimental data for hematocrit ratio and hematocrit and velocity profiles.     

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.     

Co-firing of eucalyptus bark and rubberwood sawdust in a swirling fluidized-bed combustor using an axial flow swirler

Co-combustion of eucalyptus bark (as shredded fuel) and rubberwood sawdust was conducted in a swirling fluidized-bed combustor using a 22-vane axial flow swirler. During the co-firing tests, the fuel blend feed rate was maintained at about 60 kg/h, while the mass/energy fraction of the blended fuels was variable. Excess air supply ranged from 20% to 80% for each fuel option via variation of primary air, while secondary air was injected tangentially into the bed splash zone at a constant flowrate. For comparison, tests with pure rubberwood sawdust at similar operating conditions were performed as well. Temperature and concentrations of O₂, CO and NO were measured along radial and axial directions in the combustor as well as at stack. For all fuel firing options, the radial and axial temperature profiles in the reactor were found to be weakly dependent on operating conditions. However, the gas concentration profiles exhibited apparent effects of fuel properties, excess air and secondary air injection, which resulted in variable emission characteristics of the combustor. For the sawdust energy fraction in the fuel blend of about 0.85, CO and NO emissions can be controlled at acceptable levels (both complying with the national emission limits) by maintaining excess air between 50% and 55%. Under such conditions, the co-combustion of high-moisture eucalyptus bark and rubberwood sawdust in the proposed combustor occurs in a stable regime with high, 99.6%, combustion efficiency.