Bolus contaminant dispersion in oscillating flow in curved tubes.

The investigation of longitudinal dispersion of tracer substances in unsteady flows has biomechanical application in the study of heat and mass transport within the bronchial airways during normal, abnormal, and artificial pulmonary ventilation. To model the effects of airway curvature on intrapulmonary gas transport, we have measured local gas dispersion in axially uniform helical tubes of slight pitch during volume-cycled oscillatory flow. Following a small argon bolus injection into the flow field, the time-averaged effective diffusion coefficient (Deff/Dmol) for axial transport of the contaminant was evaluated from the time-dependent local argon concentration measured with a mass spectrometer. The value of (Deff/Dmol) is extracted from the curve of concentration versus time by two techniques yielding identical results. Experiments were conducted in two helical coiled tubes (delta = 0.031, lambda = 0.022 or delta = 0.085, lambda = 0.060) over a range of 2 < alpha < 15, 3 < A < 15, where delta is the ratio of tube radius to radius of curvature, lambda is the ratio of pitch height to radius of curvature, alpha is the Womersley parameter or dimensionless frequency, and A is the stroke amplitude or dimensionless tidal volume. Experimental results show that, when compared to transport in straight tubes, the effective diffusivity markedly increases in the presence of axial curvature. Results also compare favorably to mathematical predictions of bolus dispersion in a curved tube over the ranges of frequency and tidal volume studied.     

Blood flow in small curved tubes

Blood flow in small curved tubes is modeled by the two-fluid model where a relatively cell-free fluid layer envelops a fluid core of higher viscosity. The parameters in the model are successfully curve fitted to experimental data for straight tubes. The curved tube equations are then solved by perturbation theory. It was found that curvature in general lowers the tube resistance, but increases the shear stress near the inside wall.     

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.     

Steady flow through collapsible tubes: measurements of flow and geometry.

Compliant tubes attain a complex three-dimensional geometry when the external pressure exceeds the internal pressure and the tube is partially collapsed. A new technique for remote measurement of dynamic surfaces was applied to classical experiments with collapsible tubes. This work presents measurements of the three-dimensional structure of the tube as well as pressure and flow measurements during static loading and during steady-state fluid flow. Results are shown for two tubes of the same material and internal diameter but with different wall thicknesses. The measured tube laws compare well with previously published data and suggest the possible existence of a similarity tube law. The steady flow measurements did not compare well with the one-dimensional theoretical predictions.     

Numerical schemes for unsteady fluid flow through collapsible tubes

The study of fluid flow through compliant tubes is a fluid-structure type problem, in which a dynamic equilibrium is maintained between the fluid and the tube wall. The analogy between this flow and gas dynamics initiated the use of a number of numerical methods which were originally developed to solve compressible flow in rigid ducts. In this study we investigate the solutions obtained by applying the Lax-Wendroff and MacCormack schemes to one-dimensional incompressible flow through a straight collapsible tube. The time-evolving numerical results were compared with exact steady-state solutions. For boundary conditions which were held fixed after a prescribed rise time, the unsteady numerical solution converges to the exact steady-state solution with very good accuracy. The stability and accuracy of all the methods depend on the amount of viscous pressure loss dictated by wall friction. Flows with undamped oscillations cannot, however, be solved with these techniques.     

Blood flow in tapered tubes with biorheological applications

A steady laminar flow of blood in a uniform tapered tube has been examined. Blood rheology is assumed to be described by a polar fluid. The analytical expressions for velocities (both axial and radial), total angular velocity, wall shear and pressure drop have been obtained. In literature, the parameters N (coupling number) and L (length ratio) have been chosen independently. But, in the present analysis, it is found that they are interrelated. Variation of the flow variables with suspension concentration and tapered angle have been investigated. Some of the theoretical models for the flow through tapered tubes have been critically examined. The pressure-flow relationship has been studied numerically over the flow rate range 0.01-0.1 cc/sec and compared with experimental results. It has been shown that the existing experimental results are for the tapered tubes of larger diameter which correspond to the flow under Newtonian conditions. Finally, some biological implications and future developments of this theory have been indicated.     

Surfactant transport over airway liquid lining of nonuniform depth

Numerous effects (e.g., airway wall buckling, gravity, airway curvature, capillary instabilities) give rise to nonuniformities in the depth of the liquid lining of peripheral lung airways. The effects of such thickness variations on the unsteady spreading of a surfactant monolayer along an airway are explored theoretically here. Flow-induced film deformations are shown to have only a modest influence on spreading rates, motivating the use of a simplified model in which the liquid-lining depth is prescribed and the monolayer concentration satisfies a spatially inhomogeneous nonlinear diffusion equation. Two generic situations are considered: spreading along a continuous annular liquid lining of nonuniform depth, and spreading along a rivulet that wets the airway wall with zero contact angle. In both cases, transverse averaging at large times yields a one-dimensional approximation of axial spreading that is valid for the majority of the monolayer. However, a localized monolayer remains persistently two dimensional in a region at its leading edge having axial length scales comparable to the length scale of transverse depth variation. It is also shown how the transverse spreading of a monolayer may be arrested as it approaches a static contact line at the edge of a rivulet. Implications for Surfactant Replacement Therapy are discussed.     

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.     

Physical Model of the Dynamic Instability in an Expanding Cell Culture

Collective cell migration is of great significance in many biological processes. The goal of this work is to give a physical model for the dynamics of cell migration during the wound healing response. Experiments demonstrate that an initially uniform cell-culture monolayer expands in a nonuniform manner, developing fingerlike shapes. These fingerlike shapes of the cell culture front are composed of columns of cells that move collectively. We propose a physical model to explain this phenomenon, based on the notion of dynamic instability. In this model, we treat the first layers of cells at the front of the moving cell culture as a continuous one-dimensional membrane (contour), with the usual elasticity of a membrane: curvature and surface-tension. This membrane is active, due to the forces of cellular motility of the cells, and we propose that this motility is related to the local curvature of the culture interface; larger convex curvature correlates with a stronger cellular motility force. This shape-force relation gives rise to a dynamic instability, which we then compare to the patterns observed in the wound healing experiments.     

Convection-diffusion interaction for oscillatory flow in a tapered tube

Transport of soluble material is analyzed for volume-cycle oscillatory flow in a tapered tube. The equations of motion are solved using a regular perturbation method for small taper angle and order unity amplitude over a range of the Womersley parameter. The transport equation is also solved by a regular perturbation method where uniform end concentrations and no wall flux are assumed. The time-averaged axial transport of solute is calculated for several tapered tubes. There is substantial modification of transport compared to the straight tube case and the results are interpreted with respect to pulmonary gas exchange.     

Blood flow in stented arteries: a parametric comparison of strut design patterns in three dimensions

BACKGROUND: Restenosis after stent implantation varies with stent design. Alterations in secondary flow patterns and wall shear stress (WSS) can modulate intimal hyperplasia via their effects on platelet and inflammatory cell transport toward the wall, as well as direct effects on the endothelium. METHOD OF APPROACH: Detailed flow characteristics were compared by estimating the WSS in the near-strut region of realistic stent designs using three-dimensional computational fluid dynamics (CFD), under pulsatile high and low flow conditions. The stent geometry employed was characterized by three geometric parameters (axial strut pitch, strut amplitude, and radius of curvature), and by the presence or lack of the longitudinal connector. RESULTS: Stagnation regions were localized around stent struts. The regions of low WSS are larger distal to the strut. Under low flow conditions, the percentage restoration of mean axial WSS between struts was lower than that for the high flow by 10-12%. The largest mean transverse shear stresses were 30-50% of the largest mean axial shear stresses. The percentage restoration in WSS in the models without the longitudinal connector was as much as 11% larger than with the connector The mean axial WSS restoration between the struts was larger for the stent model with larger interstrut spacing. CONCLUSION: The results indicate that stent design is crucial in determining the fluid mechanical environment in an artery. The sensitivity of flow characteristics to strut configuration could be partially responsible for the dependence of restenosis on stent design. From a fluid dynamics point of view, interstrut spacing should be larger in order to restore the disturbed flow; struts should be oriented to the flow direction in order to reduce the area of flow recirculation. Longitudinal connectors should be used only as necessary, and should be parallel to the axis. These results could guide future stent designs toward reducing restenosis.     

Persistent Symmetry Frustration in Pollen Tubes

Pollen tubes are extremely rapidly growing plant cells whose morphogenesis is determined by spatial gradients in the biochemical composition of the cell wall. We investigate the hypothesis (MP) that the distribution of the local mechanical properties of the wall, corresponding to the change of the radial symmetry along the axial direction, may lead to growth oscillations in pollen tubes. We claim that the experimentally observed oscillations originate from the symmetry change at the transition zone, where both intervening symmetries (cylindrical and spherical) meet. The characteristic oscillations between resonating symmetries at a given (constant) turgor pressure and a gradient of wall material constants may be identified with the observed growth-cycles in pollen tubes.     

Developing steady laminar flow through uniform straight tubes with varying wall cross curvature

Numerical calculations are used to determine not only the wall shear stress but also the entry length in a laminar steady flow of an incompressible Newtonian fluid. The fluid is conveyed through rigid straight tubes with axially uniform cross sections, which mimic collapsed vessels. For each tube configuration, the "Navier-Stokes" equations are solved using the finite element method. The numerical tests are performed with the same value of the volume flow-rate whatever the tube configuration for three "Reynolds numbers". The wall shear stress is computed and determined along the axis of the tube, then the entry length is estimated by introducing two indexes by using: (i) the axial fluid velocity, and (ii) the wall shear stress. The results are analysed in order to exhibit the mechanical environment of cultured endothelial cells in the flow chamber for which the test conditions will be well-defined. For example, in a tube configuration where the opposite walls are in contact for which the inner perimeter and the area of the cross section are respectively given by 45 mm and 37.02 mm(2), the computed entry lengths with the criteria defined by (i) and (ii) are equals to about 118 and 126 mm, respectively for R(e0) = 500.     

On the propagation of a disturbance through a viscous liquid flowing in a distensible tube of appreciable mass

The velocity of propagation of a disturbance wave in a liquid flowing in a distensible tube is computed. The mathematical model is more general than those used in previous analyses: the tube wall properties are realistic; the convective part of the axial inertia forces is taken into account; radial inertia forces of both the fluid and tube wall are present; viscous stresses are present. Four parameters influencing the velocity of propagation are obtained and discussed. Curves are plotted illustrating the effects of the parameters. Contrary to the results of previous analyses, viscous effects are shown to be appreciable in blood flow. It is also shown that radial inertia effects can be important in laboratory set-ups. The material presented in this paper was adapted from the Ph.D. thesis written by the author at Harvard University.     

Numerical simulation of steady flow in a model of the aortic bifurcation.

The governing equations of steady flow of an incompressible viscous fluid through a 3-D model of the aortic bifurcation are solved with the finite element method. The effect of Reynolds number on the flow was studied for a range including the physiological values (200 < or = Re < or = 1600). The symmetrical bifurcation, with a branch angle of 70 degrees and an area ratio of 0.8, includes a tapered transition zone. Secondary flows induced by the tube curvature are observed in the daughter tubes. Transverse currents in the transition zone are generated by the combined effect of diverging and converging walls. Flow separation depends on both the Reynolds number and the inlet wall shear.     

A numerical study of flow in curved tubes simulating coronary arteries

Numerical simulations of pulsatile flow in coronary arteries which take into account the curvature associated with the bending of arteries over the surface of the heart are presented for resting, excited and drug induced states. The study was motivated by reported observations of atherosclerotic plaque localization on the inner curvature of coronary arteries. The simulated flow field appears quasi-steady under resting conditions with wall shear stress always highest on the outside wall and only a single secondary flow vortex in the half tube. However, reversal of wall shear stress direction at the inside wall does occur under resting flow conditions and this is not a quasi-steady characteristic. The flow field is markedly unsteady under excited conditions with wall shear stress sometimes peaking on the inside wall and an increase in the magnitude of wall shear stress reversal on the inside wall. However, only a single secondary flow vortex in the half tube is observed. Implications of the simulations for the role of fluid mechanics in coronary artery atherosclerosis are also discussed.     

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

The pulsatile flow of non-Newtonian fluid in a bifurcation model with a non-planar daughter branch is investigated numerically by using the Carreau-Yasuda model to take into account the shear thinning behavior of the analog blood fluid. The objective of this study is to deal with the influence of the non-Newtonian property of fluid and of out-of-plane curvature in the non-planar daughter vessel on wall shear stress (WSS), oscillatory shear index (OSI), and flow phenomena during the pulse cycle. The non-Newtonian property in the daughter vessels induces a flattened axial velocity profile due to its shear thinning behavior. The non-planarity deflects flow from the inner wall of the vessel to the outer wall and changes the distribution of WSS along the vessel, in particular in systole phase. Downstream of the bifurcation, the velocity profiles are shifted toward the flow divider, and low WSS and high shear stress temporal oscillations characterized by OSI occur on the outer wall region of the daughter vessels close to the bifurcation. Secondary motions become stronger with the addition of the out-of-plane curvature induced by the bending of the vessel, and the secondary flow patterns swirl along the non-planar daughter vessel. A significant difference between the non-Newtonian and the Newtonian pulsatile flow is revealed during the pulse cycle; however, reasonable agreement between the non-Newtonian and the rescaled Newtonian flow is found. Calculated results for the pulsatile flow 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.     

Effects of wall motion and compliance on flow patterns in the ascending aorta

Helical flows have been observed in the ascending aorta in vivo, and geometric curvature has been suggested to be a major contributing factor. We employed magnetic resonance imaging (MRI) and velocity mapping to develop a computational model to examine the effects of curvature and also wall compliance and movement upon flow patterns. In the computational model, MRI-derived geometry and velocities were imposed as boundary conditions, which included both radial expansion-contraction and translational motion of the wall. The computed results were in agreement with the MR data only when full wall motion was included in the model, suggesting that the flow patterns observed in the ascending aorta arise not only from geometric curvature of the arch but also from the motion of the aorta resulting from its attachment to the beating heart.     

The effects of wall thickness, axial strain and end proximity on the pressure-area relation of collapsible tubes

Segments of silicone rubber tube were suspended between rigid pipes and subjected to slowly varying transmural pressure covering a range from slight distension to collapse with osculation. The local inside cross-sectional area at a chosen axial site was simultaneously measured via catheter by an electrical impedance method. Pressure-area relations were recorded thus at various axial sites, under varying conditions of axial tube wall tension, in tubes of two different wall thickness (0.3 and 0.4 of mean radius). Unsupported tube segment length was also varied by means of an insert device. The relations were used to calculate the variation of wave velocity with area according to Young's equation. First opposite wall contact during collapse was shown to occur at a smaller fraction of undistended circular cross-sectional area than in the thin-walled tubes investigated previously by others.