Viscous flow through slowly expanding or contracting porous walls with low seepage Reynolds number: a model for transport of biological fluids through vessels

In this article, the problem of laminar, isothermal, incompressible and viscous flow in a rectangular domain bounded by two moving porous walls, which enable the fluid to enter or exit during successive expansions or contractions, is investigated. The governing non-linear equations and their associated boundary conditions are transformed into a highly non-linear ordinary differential equation. The series solution of the problem is obtained by utilising the homotopy perturbation method. Graphical results are presented to investigate the influence of the non-dimensional wall dilation rate and seepage Reynolds number (Re) on the velocity, normal pressure distribution and wall shear stress. Since the transport of biological fluids through contracting or expanding vessels is characterised by low seepage Res, the current study focuses on the viscous flow driven by small wall contractions and expansions of two weakly permeable walls.     

Viscous flow through slowly expanding or contracting porous walls with low seepage Reynolds number: a model for transport of biological fluids through vessels

In this article, the problem of laminar, isothermal, incompressible and viscous flow in a rectangular domain bounded by two moving porous walls, which enable the fluid to enter or exit during successive expansions or contractions, is investigated. The governing non-linear equations and their associated boundary conditions are transformed into a highly non-linear ordinary differential equation. The series solution of the problem is obtained by utilising the homotopy perturbation method. Graphical results are presented to investigate the influence of the non-dimensional wall dilation rate and seepage Reynolds number (Re) on the velocity, normal pressure distribution and wall shear stress. Since the transport of biological fluids through contracting or expanding vessels is characterised by low seepage Res, the current study focuses on the viscous flow driven by small wall contractions and expansions of two weakly permeable walls.     

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.     

Fluid Dynamics of Ventricular Filling in the Embryonic Heart

The vertebrate embryonic heart first forms as a valveless tube that pumps blood using waves of contraction. As the heart develops, the atrium and ventricle bulge out from the heart tube, and valves begin to form through the expansion of the endocardial cushions. As a result of changes in geometry, conduction velocities, and material properties of the heart wall, the fluid dynamics and resulting spatial patterns of shear stress and transmural pressure change dramatically. Recent work suggests that these transitions are significant because fluid forces acting on the cardiac walls, as well as the activity of myocardial cells that drive the flow, are necessary for correct chamber and valve morphogenesis. In this article, computational fluid dynamics was used to explore how spatial distributions of the normal forces acting on the heart wall change as the endocardial cushions grow and as the cardiac wall increases in stiffness. The immersed boundary method was used to simulate the fluid-moving boundary problem of the cardiac wall driving the motion of the blood in a simplified model of a two-dimensional heart. The normal forces acting on the heart walls increased during the period of one atrial contraction because inertial forces are negligible and the ventricular walls must be stretched during filling. Furthermore, the force required to fill the ventricle increased as the stiffness of the ventricular wall was increased. Increased endocardial cushion height also drastically increased the force necessary to contract the ventricle. Finally, flow in the moving boundary model was compared to flow through immobile rigid chambers, and the forces acting normal to the walls were substantially different.     

Velocity field of pulsatile flow in a porous tube

This paper describes velocity fields for fully developed periodic laminar flow in a rigid tube with a porous wall. We obtained an analytical solution of the flow by the linear approximation of the Navier-Stokes equation. Unlike the previous works with a constant seepage rate along the axis, we used a wall velocity which contained hydraulic permeation constant Lp. The axial velocity profile shows a local maximum velocity near the wall at a large Womersley number alpha. This suggests that concentration polarization in porous tubular membrane may be reduced at high frequencies if a membrane device is operated under pulsatile flow conditions. The magnitude of wall permeation velocity decreases linearly along the tube axis because the damping of the pressure difference between the inside and the outside of the tube is very small.     

Peristaltic Transport of Carreau-Yasuda Fluid in a Curved Channel with Slip Effects

The wide occurrence of peristaltic pumping should not be surprising at all since it results physiologically from neuro-muscular properties of any tubular smooth muscle. Of special concern here is to predict the rheological effects on the peristaltic motion in a curved channel. Attention is focused to develop and simulate a nonlinear mathematical model for Carreau-Yasuda fluid. The progressive wave front of peristaltic flow is taken sinusoidal (expansion/contraction type). The governing problem is challenge since it has nonlinear differential equation and nonlinear boundary conditions even in the long wavelength and low Reynolds number regime. Numerical solutions for various flow quantities of interest are presented. Comparison for different flow situations is also made. Results of physical quantities are interpreted with particular emphasis to rheological characteristics.     

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.     

Analytic Approximate Solutions to the Boundary Layer Flow Equation over a Stretching Wall with Partial Slip at the Boundary

Analytic approximate solutions using Optimal Homotopy Perturbation Method (OHPM) are given for steady boundary layer flow over a nonlinearly stretching wall in presence of partial slip at the boundary. The governing equations are reduced to nonlinear ordinary differential equation by means of similarity transformations. Some examples are considered and the effects of different parameters are shown. OHPM is a very efficient procedure, ensuring a very rapid convergence of the solutions after only two iterations.     

Interaction of peristaltic flow with pulsatile flow in a circular cylindrical tube

The effect of pulsatile flow on peristaltic transport in a circular cylindrical tube is analysed. The flow of a Newtonian viscous incompressible fluid in a flexible circular cylindrical tube on which an axisymmetric travelling sinusoidal wave is imposed, is considered. The initial flow in the tube is induced by an arbitrary periodic pressure gradient. A perturbation solution with amplitude ratio (wave amplitude/tube radius) as a parameter is obtained when the frequency of the travelling wave and that of the imposed pressure gradient are equal. The interaction effects of periodic wall induced flow and periodic pressure imposed flow are visualized through the presence of substantially different components of steady and higher harmonic oscillating flow in the first order flow solution. Numerical results show a strong variation of steady state velocity profiles with boundary wave number and Reynolds number and a strong phase shift behaviour of the flow in the radial direction.     

A one-dimensional viscous-inviscid strong interaction model for flow in indented channels with separation and reattachment

A new one-dimensional model is presented for the calculation of steady and unsteady flow through an indented two-dimensional channel with separation and reattachment. It is based on an interactive boundary layer approach, where the equations for the boundary layer flow near the channel walls and for an inviscid core flow are solved simultaneously. This approach requires no semi-empirical inputs, such as the location of separation and reattachment, which is an advantage over other existing one-dimensional models. Because of the need of an inviscid core alongside the boundary layers, the type of inflow as well as the length of the channel and the value of the Reynolds number poses some limitations on the use of the new model. Results have been obtained for steady flow through the indented channel of Ikeda and Matsuzaki. In further perspective, it is discussed how the present model, in contrast to other one-dimensional flow models, can be extended to calculate the flow in nonsymmetrical channels, by considering different boundary layers on each of the walls.     

On Three-Dimensional Flow and Heat Transfer over a Non-Linearly Stretching Sheet: Analytical and Numerical Solutions

This article studies the viscous flow and heat transfer over a plane horizontal surface stretched non-linearly in two lateral directions. Appropriate wall conditions characterizing the non-linear variation in the velocity and temperature of the sheet are employed for the first time. A new set of similarity variables is introduced to reduce the boundary layer equations into self-similar forms. The velocity and temperature distributions are determined by two methods, namely (i) optimal homotopy analysis method (OHAM) and (ii) fourth-fifth-order Runge-Kutta integration based shooting technique. The analytic and numerical solutions are compared and these are found in excellent agreement. Influences of embedded parameters on momentum and thermal boundary layers are sketched and discussed.     

Peristaltic flow of a couple stress fluid in an asymmetric channel

This paper presents an analysis of the peristaltic flow of a couple stress fluid in an asymmetric channel. The asymmetric nature of the flow is introduced through the peristaltic waves of different amplitudes and phases on the channel walls. Mathematical modelling corresponding to a two-dimensional flow has been carried out. The flow analysis is presented under long wavelength and low Reynolds number approximations. Closed form solutions for the axial velocity, stream function and the axial pressure gradient are given. Numerical computations have been carried out for the pressure rise per wavelength, friction forces and trapping. It is noted that there is a decrease in the pressure when the couple stress fluid parameter increases. The variation of the couple stress fluid parameter with the size of the trapped bolus is also similar to that of pressure. Furthermore, the friction force on the lower channel wall is greater than that on the upper channel wall.     

Effect of Joule Heating and Thermal Radiation in Flow of Third Grade Fluid over Radiative Surface

This article addresses the boundary layer flow and heat transfer in third grade fluid over an unsteady permeable stretching sheet. The transverse magnetic and electric fields in the momentum equations are considered. Thermal boundary layer equation includes both viscous and Ohmic dissipations. The related nonlinear partial differential system is reduced first into ordinary differential system and then solved for the series solutions. The dependence of velocity and temperature profiles on the various parameters are shown and discussed by sketching graphs. Expressions of skin friction coefficient and local Nusselt number are calculated and analyzed. Numerical values of skin friction coefficient and Nusselt number are tabulated and examined. It is observed that both velocity and temperature increases in presence of electric field. Further the temperature is increased due to the radiation parameter. Thermal boundary layer thickness increases by increasing Eckert number.     

Salt-induced Contraction of Bacterial Cell Walls

Intact Bacillus megaterium cells were found to contract as much as 26% in terms of dextran-impermeable volume when transferred from water to unbuffered, non-plasmolyzing NaCl solutions. This shrinkage appeared to be primarily due to electrostatic wall contraction rather than to any osmotic response of the cells. A variety of salts (but not sucrose) added to water suspensions of isolated cell walls caused protons to be released from the walls with resultant lowering of suspension pH and contraction of the structures. In effect, B. megaterium walls behaved as flexible, amphoteric polyelectrolytes, and their compactness in aqueous suspensions was affected by changes in environmental ionic strength and pH. Isolated walls were most compact in low ionic strength media with a pH of about 4, a value close to the apparent isoelectric pH of wall peptidoglycan. Electrostatic attractions appeared to play a major role in determining the compactness of highly contracted walls, and the walls responded to increased environmental ionic strength by expanding. In contrast, electrostatic repulsions were dominant in highly expanded walls, and increased environmental ionic strength induced wall contraction. Walls of whole bacteria also shrank when the cells were plasmolyzed. This second type of contraction seemed to result from relief of wall tension during plasmolysis, and it could be induced with nonionic solutes. Thus, cell wall tone in B. megaterium appeared to be set both by mechanical tension and by electrostatic interactions among wall ions.     

On the influence of wall properties and Poiseuille flow in peristalsis

The present study extends the two-dimensional analysis of peristaltic motion by Fung and Yih to include an elastic or viscoelastic wall, and a Poiseuille flow. This fluid-solid interaction problem is investigated by considering equations of motion of both the fluid and the deformable boundaries. The wall characteristics appear in their equations of motion, which are solved to represent boundary conditions of fluid motion. The influence of Poiseuille flow on pure peristalisis is also investigated.

The phenomenon of the ‘mean flow reversal’ is found to exist both at the center and at the boundaries of the channel. When the walls of the channel are elastic, pure peristalsis involves flow reversal only at the center. This position may shift drastically to the boundaries, if viscous damping forces are considered.     

Resistance to Water Flow in Xylem Vessels

Experimental data on flow resistances in xylem vessels withdifferent lumen wall surface sculptures are presented. The techniqueinvolved using determinable forces at menisci to pull waterthrough isolated undamaged metaxylem and protoxylem vesselswhich were empty but had water-saturated walls. In the horizontalorientation, the surface tension forces moved the water at velocitiesthat the resisting viscous forces at the vessel walls wouldallow since inertial forces were found negligible. A high speedcamera was used to determine the meniscus translation rates.Vessel diameters as well as average dimensions of the microscopicinternal surface irregularities were measured with respect toaxial position from the inlet. From these, flow resistanceswere determined in terms of dimensionless friction factor, f,as functions of Reynolds number, Re. It was found that, at certain helical ring thicknesses and spacing,resistance to flow was lowest. Deviations from these parametervalues cause dramatic increases in resistance to flow. Resultsare applicable to normal flow in plants, i.e. without meniscipresent.     

Three-dimensional convective alveolar flow induced by rhythmic breathing motion of the pulmonary acinus

Low Reynolds number flows (Re<1) in the human pulmonary acinus are often difficult to assess due to the submillimeter dimensions and accessibility of the region. In the present computational study, we simulated three-dimensional alveolar flows in an alveolated duct at each generation of the pulmonary acinar tree using recent morphometric data. Rhythmic lung expansion and contraction motion was modeled using moving wall boundary conditions to simulate realistic sedentary tidal breathing. The resulting alveolar flow patterns are largely time independent and governed by the ratio of the alveolar to ductal flow rates, Qa/Qd. This ratio depends uniquely on geometrical configuration such that alveolar flow patterns may be entirely determined by the location of the alveoli along the acinar tree. Although flows within alveoli travel very slowly relative to those in acinar ducts, 0.021%相似文献   

Some features of oscillatory flow in a model bifurcation

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

Transport of macromolecules in arterial wallin vivo: A mathematical model and analytical solutions

A mathematical model has been developed to simulatein vivo transmural accumulation of an intravenously injected tracer in the aortic wall of experimental animals. Parameters have been included to represent the following processes that affect tracer distribution: permeation of the blood-tissue interface, diffusion through the layers of the artery wall,convective solute drag through the same, and degradation. Of particular interest for thein vivo situation situation is the inclusion of boundary conditions that account for the variation in the plasma concentration of injected tracer as a function of time. Two analytical solutions are presented. The first describes a system in which two boundaries must be delineated; it pertains if the tracer is allowed to circulate until it enters the avascular media of the artery wall both across its luminal boundary and from the capillaries in its outer layer. The second applies to shorter duration experiments in which entry across only the luminal boundary is considered. A limiting case of the solution for short circulation times is presented, compared with a previously published solution, and examined for its potential utility in parameter estimation. Because of its treatment of time-dependent boundary conditions, the model has unique application toin vivo experiments related to macromolecular transport in atherosclerosis that may otherwise elude proper interpretation. This work was supported by National Institutes of Health Grants HL-29582 and HL-07242.     

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.