One-dimensional steady inviscid flow through a stenotic collapsible tube

A one-dimensional inviscid solution for flow through a compliant tube with a stenosis is presented. The model is used to represent an artery with an atherosclerotic plaque and to investigate a range of conditions for which arterial collapse may occur. The coupled equations for flow through collapsible tubes are solved using a Runge-Kutta finite difference scheme. Quantitative results are given for specific physiological parameters including inlet and outlet pressure, flow rate, stenosis size, length and stiffness. The results suggest that high-grade stenotic arteries may exhibit collapse with typical physiological pressures. Critical stenoses may cause choking of flow at the throat followed by a transition to supercritical flow with tube collapse downstream. Greater amounts of stenosis produced a linear reduction of flow rate and a shortening of the collapsed region. Changes in stenosis length created proportional changes in the length of collapse. Increasing the stiffness of the stenosis to a value greater than the nominal tube stiffness caused a greater amount of flow limitation and more negative pressures, compared to a stenosis with constant stiffness. These findings assist in understanding the clinical consequences of flow through atherosclerotic arteries.     

Numerical analysis for stability and self-excited oscillation in collapsible tube flow.

This paper describes numerical analysis of collapsible tube flow based on the one-dimensional distributed parameter model of Hayashi. In the present model the effect of flow separation at the collapsed part is replaced with simple viscous friction along the tube, so no ad-hoc modeling for flow separation in former studies is required. A stable semi-implicit numerical procedure based on the SIMPLE method is developed for the problem of flow and tube interaction. The numerical result for a characteristic self-excited oscillation agrees qualitatively with the experimental result. Nonlinear stability of the steady state dependent on the amplitude of the disturbance is numerically investigated and the result is compared with the linear stability analysis based on the former lumped parameter model. Finally, initiation of the self-excited oscillation is examined by applying the initial disturbance at the upstream end of the tube. The disturbance propagates in the downstream direction and is amplified to the self-excited oscillation.     

A mathematical model of flow through a collapsible tube--I. Model and steady flow results

A one-dimensional model is presented to describe the flow through a collapsible tube whose ends are tethered to rigid tubes and which is enclosed in a pressurized chamber. Results are presented for the special case of steady flow. Predicted pressure drop versus flow rate (delta P-Q) characteristics agree qualitatively with available experimental data. The significance of the model and of various physical parameters, in regard to the shape of these characteristics, is discussed.     

Flow of urine through the ureter: a collapsible, muscular tube undergoing peristalsis

In steady flow through nonuniform collapsible tubes a key concept is the compressive zone, at which flow limitation can occur at both high and low Reynolds numbers. Ureteral peristalsis can be considered as a series of compressive zones, corresponding to waves of active muscular contraction, that move at near-constant speed along the ureter towards the bladder. One-dimensional, lubrication-theory analysis shows that peristalsis can pump urine from kidney into the bladder only at relatively low mean rates of urine flow. Under these circumstances isolated boluses of urine are propelled steadily through the ureter (assumed uniform) by the contraction waves. At higher mean rates of flow the behavior depends on whether the frequency of peristalsis is higher or lower than a critical value. For frequencies above the critical value steady propagation of boluses that are in contact with contraction waves at both ends is possible. As the flow rate rises the urine begins to leak through the contraction waves and steady peristaltic flow breaks down. There is an upper limit to the mean flow rate that can be carried by steady peristalsis, which depends on the mechanical properties of the ureter. At high flow rates the peristaltic contractions do not pump but hinder the flow of urine through the ureter.     

Sensitivity and specificity of the computational model for maximal expiratory flow

The computational model for forced expiratory flow from human lungs of Lambert and associates (J. Appl. Physiol.: Respirat. Environ. Exercise Physiol. 52: 44-56, 1982) was used to investigate the sensitivity of maximal expiratory flow to lung properties. It was found that maximal flow is very sensitive to recoil pressure and airway areas but not very sensitive to lung volume, airway compliance, and airway length. Linear programming was used to show that a given air flow-pressure curves was compatible with a fairly wide range of airway properties. Additional data for maximal flow with a He-O2 mixture narrowed the range somewhat. It was shown that the flow-pressure curve contains more information about central than peripheral airways and that information about the latter is obtainable only from flows at recoils less than 2 cmH2O. Parameter ranges compatible with individual flow-pressure curves showed differences that demonstrated that such curves give some indication of individual central airway properties.     

A one-dimensional unsteady separable and reattachable flow model for collapsible tube-flow analysis

Taking into account both flow separation and reattachment observed in available experimental results on flows in a quasi-two-dimensional channel, we present a one-dimensional unsteady flow model, which is applicable to a flow in a collapsible tube. The flow model has been derived from the two-dimensional Navier-Stokes equations by introducing the concept of a dividing streamline, which divides a separated flow into a jet and a dead-water zone. We also present a criterion for the determination of a separation point. Numerical results show that the locations of the predicted separation points agree well with the experimental data. The predicted static pressure of the separated flow is almost constant downstream of the separation point and increases quickly just before the reattachment point as observed in the experiment. Finally, using the present flow model and the separation criterion, we examine the oscillatory behavior of an unsteady flow in a symmetric channel whose walls move sinusoidally.     

Microcontinuum model for pulsatile blood flow through a stenosed tube

The effects of polar nature of blood and pulsatility on flow through a stenosed tube have been analysed by assuming blood as a micropolar fluid. Linearized solutions of basic equations are obtained through consecutive applications of finite Hankel and Laplace transforms. The analytical expressions for axial and particle angular velocities, wall shear stress, resistance to flow and apparent viscosity have been obtained. The axial velocity profiles for Newtonian and micropolar fluids have been compared. The interesting observation of this analysis is velocity, in certain parts of cycle, for micropolar fluid is higher than Newtonain fluid. Variation of apparent viscosity eta a with tube radius shows both inverse Fahraeus-Lindqvist and Fahraeus-Lindqvist effects. Finally, the resistance to flow and wall shear stress for normal and diseased blood have been computed and compared.     

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.     

Influence of expiratory flow on closing capacity at low expiratory flow rates.

Single-breath oxygen (SBO2) tests at expiratory flow rates of 0.2, 0.5, and 1.01/s were performed by 10 normal subjects in a body plethysmograph. Closing capacity (CC)--the absolute lung volume at which phase IV began--increased significantly with increases in flow. Five subjects were restudied with a 200-ml bolus of 100% N2 inspired from residual volume after N2 washout by breathing 100% O2 and similar results were obtained. An additional five subjects performed SBO2 tests in the standing, supine, and prone positions; closing volume (CV)--the lung volume above residual volume at which phase IV began--also increased with increases of expiratory flow. The observed increase in CC with increasing flow did not appear to result from dependent lung regions reaching some critical "closing volume" at a higher overall lung volume. In normal subjects, the phase IV increase in NI concentration may be caused by the asynchronous onset of flow limitation occurring initially in dependent regions.     

A transient luminal chitinous matrix is required to model epithelial tube diameter in the Drosophila trachea

Epithelial tubes are found in many vital organs and require uniform and correct tube diameters for optimal function. Tube size depends on apical membrane growth and subapical cytoskeletal reorganization, but the cues that coordinate these events to ensure functional tube shape remain elusive. We find that epithelial tubes in the Drosophila trachea require luminal chitin polysaccharides to attain the correct diameter. Tracheal chitin forms a broad transient filament within the tubes during the restricted period of expansion. Loss of chitin causes tubular constrictions and cysts associated with irregular subapical cytoskeletal organization, without affecting epithelial integrity and polarity. Analysis of previously identified tube expansion mutants in genes encoding septate junction proteins further suggests that septate junction components may function in tubulogenesis through their role in luminal matrix assembly. We propose that the transient luminal protein/polysaccharide matrix is sensed by the epithelial cells and coordinates cytoskeletal organization to ensure uniform lumen diameter.     

Numerical simulation of pulsatile flow in a compliant curved tube model of a coronary artery

The endothelial cells (ECs) lining a blood vessel wall are exposed to both the wall shear stress (WSS) of blood flow and the circumferential strain (CS) of pulsing artery wall motion. These two forces and their interaction are believed to play a role in determining remodeling of the vessel wall and development of arterial disease (atherosclerosis). This study focused on the WSS and CS dynamic behavior in a compliant model of a coronary artery taking into account the curvature of the bending artery and physiological radial wall motion. A three-dimensional finite element model with transient flow and moving boundaries was set up to simulate pulsatile flow with physiological pressure and flow wave forms characteristic of the coronary arteries. The characteristic coronary artery curvature and flow conditions applied to the simulation were: aspect ratio (lambda) = 10, diameter variation (DV) = 6 percent, mean Reynolds number (Re) = 150, and unsteadiness parameter (alpha) = 3. The results show that mean WSS is about 50 percent lower on the inside wall than the outside wall while WSS oscillation is stronger on the inside wall. The stress phase angle (SPA) between CS and WSS, which characterizes the dynamics of the mechanical force pattern applied to the endothelial cell layer, shows that CS and WSS are more out of phase in the coronaries than in any other region of the circulation (-220 deg on the outside wall, -250 deg on the inside wall). This suggests that in addition to WSS, SPA may play a role in localization of coronary atherosclerosis.     

Longitudinal tension variation in collapsible channels: a new mechanism for the breakdown of steady flow.

There are several mechanisms potentially involved in the breakdown of steady fluid flow in a collapsible tube under external pressure. Here we investigate one that has received little attention in the past: the fact that the longitudinal tension in the tube wall, T, decreases with distance downstream as a consequence of the viscous shear stress exerted by the fluid. If the tube is long enough, or the initial tension small enough, T may fall to zero before the end of the collapsible tube, and unsteady motion would presumably then ensue; this is what we mean by "breakdown." We study the phenomenon theoretically, when the flow Reynolds number is of order one, using lubrication theory in a symmetric two-dimensional channel in which the collapsible tube is replaced by membranes occupying a segment of each wall. The resulting nonlinear ordinary differential equations are solved numerically for values of the dimensionless parameters that cover all the qualitatively different types of solution (e.g., in which the channel is distended over all its length, collapsed over all its length, or distended in the upstream part and collapsed downstream). Reducing the longitudinal tension has a marked effect on the shape of the collapsible segment, causing it to become much more deformed for the same flow rate and external pressure. Indeed, the wall slope is predicted to become very large when the downstream tension is very small, so the model is not self-consistent then. Nevertheless, the parameter values for which T becomes zero are mapped out and are expected to be qualitatively useful.(ABSTRACT TRUNCATED AT 250 WORDS)     

Laser-Doppler measurements of velocities just downstream of a collapsible tube during flow-induced oscillations

The flow field less than one diameter downstream of the end of a collapsible tube executing self excited oscillations was examined using a two-component fiber-optic laser-Doppler anemometer. The time-averaged Reynolds number of the flow was 11,000. With the tube oscillating periodically, results obtained during many cycles of oscillation were combined to yield surface plots of the axial component over the cross section at 16 phases of the cycle. By combining measurements obtained with the laser probe in two different orientations, secondary flow vectors over the cross section were likewise constructed for 16 phases. The measurements showed strongly phasic turbulence intensity, with the phase of high intensity coinciding with the time of maximal tube collapse. Reverse flow occurred during much of the cycle, at places in the cross section that agree with our previous observations of laminar and turbulent steady flow through a rigid simulated collapsed tube.     

A numerical calculation of flow in a curved tube model of the left main coronary artery

The flow pattern in the left main coronary artery has been calculated using an idealized geometry and by numerically solving the full Navier-Stokes equations for a Newtonian fluid. Two different forms for the entrance velocity profile were used, one a time-varying, flat profile and the other a time-varying, less flat velocity profile. The results obtained demonstrate the presence of secondary motions for conditions simulating flow in the left main coronary artery, with maximum secondary flow velocities being on the order of three to four percent of the maximum axial velocity. This secondary flow phenomenon has an important influence on the wall shear stress distribution, in spite of the fact that there is virtually no alteration in the axial velocity profile. The maximum ratio of the outer wall shear stress to that on the inner wall is 1.4 at a Reynolds number of Re = 270, and it increases with increasing Reynolds number, reaching a value of 1.7 at Re = 810. Although there are significant differences in the results in the immediate vicinity of the inlet for the two different forms of the entrance velocity profile used, this difference does not persist far into the tube. Independent of the choice of the entrance velocity profile, it appears that there will be significant secondary flow effects on the wall shear stress.     

Steady flow through a double converging-diverging tube model for mild coronary stenoses

Velocity profiles and the pressure drop across two mild (62 percent) coronary stenoses in series have been investigated numerically and experimentally in a perspex-tube model. The mean flow rate was varied to correspond to a Reynolds number range of 50-400. The pressure drop across two identical (62 percent) stenoses show that for low Reynolds numbers the total effect of two stenoses equals that of two single stenoses. A reduction of 10 percent is found for the higher Reynolds numbers investigated. Numerical and experimental results obtained for the velocity profiles agree very well. The effect of varying the converging angle of a single mild (62 percent) coronary stenosis on the fluid flow has been determined numerically using a finite element method. Pressure-flow relation, especially with respect to relative short stenoses, is discussed.