Flow Chamber流动试验模型

从Navier-Stokes方程出发,运用无量纲化及虚宗量Bessel函数法求解,导得Flow Chamber系统中流体流动的速度场、切应力场等分布规律。并且,本文公式简化为定常状态时,无论是和Cao等人的理论结果,还是和他们的实验结果比较,都是相当一致的。     

Measurement and prediction of flow through a replica segment of a mildly atherosclerotic coronary artery of man

Pressure distributions were measured along a hollow vascular axisymmetric replica of a segment of the left circumflex coronary artery of man with mildly atherosclerotic diffuse disease. A large range of physiological Reynolds numbers from about 60 to 500, including hyperemic response, was spanned in the flow investigation using a fluid simulating blood kinematic viscosity. Predicted pressure distributions from the numerical solution of the Navier-Stokes equations were similar in trend and magnitude to the measurements. Large variations in the predicted velocity profiles occurred along the lumen. The influence of the smaller scale multiple flow obstacles along the wall (lesion variations) led to sharp spikes in the predicted wall shear stresses. Reynolds number similarity was discussed, and estimates of what time averaged in vivo pressure drop and shear stress might be were given for a vessel segment.     

Flow mixing enhancement from balloon pulsations in an intravenous oxygenator

Computational investigations of flow mixing and oxygen transfer characteristics in an intravenous membrane oxygenator (IMO) are performed by direct numerical simulations of the conservation of mass, momentum, and species equations. Three-dimensional computational models are developed to investigate flow-mixing and oxygen-transfer characteristics for stationary and pulsating balloons, using the spectral element method. For a stationary balloon, the effect of the fiber placement within the fiber bundle and the number of fiber rings is investigated. In a pulsating balloon, the flow mixing characteristics are determined and the oxygen transfer rate is evaluated. For a stationary balloon, numerical simulations show two well-defined flow patterns that depend on the region of the IMO device. Successive increases of the Reynolds number raise the longitudinal velocity without creating secondary flow. This characteristic is not affected by staggered or non-staggered fiber placement within the fiber bundle. For a pulsating balloon, the flow mixing is enhanced by generating a three-dimensional time-dependent flow characterized by oscillatory radial, pulsatile longitudinal, and both oscillatory and random tangential velocities. This three-dimensional flow increases the flow mixing due to an active time-dependent secondary flow, particularly around the fibers. Analytical models show the fiber bundle placement effect on the pressure gradient and flow pattern. The oxygen transport from the fiber surface to the mean flow is due to a dominant radial diffusion mechanism, for the stationary balloon. The oxygen transfer rate reaches an asymptotic behavior at relatively low Reynolds numbers. For a pulsating balloon, the time-dependent oxygen-concentration field resembles the oscillatory and wavy nature of the time-dependent flow. Sherwood number evaluations demonstrate that balloon pulsations enhance the oxygen transfer rate, even for smaller flow rates.     

Kinematic and Dynamic Characteristics of Pulsating Flow in 180^o Tube

Kinematic and dynamic characteristics of pulsating flow in a model of human aortic arch are obtained by a computational analysis. Three-dimensional flow processes are summarized by pressure distributions on the symmetric plane together with velocity and pressure contours on a few cross sections for systolic acceleration and deceleration. Without considering the effects of aortic tapering and the carotid arteries, the development of tubular boundary layer with centrifugal forces and pulsation are also analyzed for flow separation and backflow during systolic deceleration.     

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.     

The inverse Womersley problem for pulsatile flow in straight rigid tubes

In this study a numerical solution for the problem of pulsating flow in rigid tubes is described. The method applies to the case of known flow rate waveform, as opposed to Womersley solution where the pressure gradient was the known quantity. The solution provides the pressure gradient and wall shear stress waveforms as well as the instantaneous velocity profiles. Results show that the method can be used to study the blood flow characteristics in large arteries.     

Effects of couple stresses on the blood flow through an artery with mild stenosis

An analysis of the effects of couple stresses on the blood flow through thin artery in the presence of very mild stenosis has been carried out with the help of two nondimensional parameters, alpha (the length ratio parameter) and eta (the parameter characterizing the antisymmetric property of the couple stress tensor). It is shown that an increase in the couple stress (small value of alpha and eta), increases the resistance to the flow and the wall shear stress. These characteristics are further enhanced by the presence of the stenosis.     

Flow-induced wall shear stress in abdominal aortic aneurysms: Part I--steady flow hemodynamics

Numerical predictions of blood flow patterns and hemodynamic stresses in Abdominal Aortic Aneurysms (AAAs) are performed in a two-aneurysm, axisymmetric, rigid wall model using the spectral element method. Homogeneous, Newtonian blood flow is simulated under steady conditions for the range of Reynolds numbers 10 < or =Re < or =2265. Flow hemodynamics are quantified by calculating the distributions of wall pressure (p(w)), wall shear stress (tau(w)), Wall Shear Stress Gradient (WSSG). A correlation between maximum values of hemodynamic stresses and Reynolds number is established, and the spatial distribution of WSSG is considered as a hemodynamic force that may cause damage to the arterial wall at an intermediate stage of AAA growth. The temporal distribution of hemodynamic stresses in pulsatile flow and their physical implications in AAA rupture are discussed in Part II of this paper.     

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.     

An Analytical Model of the Counter-Current Heat Exchange Phenomena

An analytical model for the counter-current heat exchange mechanism in animals has been formulated and a solution has been obtained. The nondimensional parameters that govern the mechanism have been determined in terms of the properties of the animal. The normalized temperatures are functions of normalized distance and, in general, three nondimensional heat transfer conductances. Graphical results are presented for two representative physiological systems. These results allow a delineation of those situations in which counter-current heat transfer is important, and also a quantitative prediction of the heat transfer and temperature distributions. The theory is compared to the available experimental results.     

Enhancing Cell Viability with Pulsating Flow in a Hollow Fiber Bioartificial Liver

A pulsating flow of medium was used to alleviate diffusion and transport limitations in a hollow fiber bioreactor containing a human hepatoblastoma cell line. The strategy is easy to implement but effective. The pulsating flow is introduced by a solenoid pinch valve at the outlet of the bioreactor and regulated by a timing circuit. In a permeability test, the system with pulsating flow had much less membrane fouling as compared to the control, a conventional hollow fiber unit. In hepatocyte culture test runs, the pulsating-flow bioreactor demonstrated the ability to maintain a higher cell viability. Histological sections indicated significantly smaller necrotic regions in the pulsating-flow bioreactor as compared to the conventional unit.     

Analysis of flow in a cone-and-plate apparatus with respect to spatial and temporal effects on endothelial cells

Endothelial cells, covering the inner surface of vessels and the heart, are permanently exposed to fluid flow, which affects the endothelial structure and the function. The response of endothelial cells to fluid shear stress is frequently investigated in cone-plate systems. For this type of device, we performed an analytical and numerical analysis of the steady, laminar, three-dimensional flow of a Newtonian fluid at low Reynolds numbers. Unsteady oscillating and pulsating flow was studied numerically by taking the geometry of a corresponding experimental setup into account. Our investigation provides detailed information with regard to shear-stress distribution at the plate as well as secondary flow. We show that: (i) there is a region on the plate where shear stress is almost constant and an analytical approach can be applied with high accuracy; (ii) detailed information about the flow in a real cone-plate device can only be obtained by numerical simulations; (iii) the pulsating flow is quasi-stationary; and (iv) there is a time lag on the order of 10(-3) s between cone rotation and shear stress generated on the plate.     

An experimental study of pressure losses in pulsatile flows through rigid and pulsating stenosis

The time-dependent pressure curves of a pulsatile flow across rigid and pulsating stenoses were investigated experimentally in a laboratory simulator of the outflow tract of the heart right ventricle. The experiments were performed within the range of physiological conditions of frequency and flow rate. The experimental setup consisted of a closed flow system which was operated by a pulsatile pump, and a test chamber which enabled checking different modes of stenosis. Rigid constrictions were simulated by means of axisymmetric blunt-ended annular plugs with moderate-to-severe area reductions. The pulsating stenosis consisted of a short starling resistor device operated by a pulsating external pressure which was synchronized by the pulsatile flow. It was found that the shape of the time-dependent pressure curve upstream of the stenosis was different in the case of rigid stenosis than in the pulsating one. Potential clinical applications of the work may relate to diagnosis of the type of stenosis in the congenital heart disease known as Tetralogy of Fallot.     

Effect of variations in pressure and flow on the geometry of isolated canine carotid arteries

A study is described in which the effects of hemodynamics on arterial geometry are investigated in vitro. A novel perfusion apparatus is employed to deliver pulsatile flow through excised canine carotid arteries under carefully controlled conditions. Data of perfused vessel diameter and arterial wall thickness are derived from the radial displacement of the pulsating vessel as measured using a scanning laser micrometer whose accuracy is determined to be 0.0125 mm (0.0005 in). The results of 30 perfusion experiments suggest that the hemodynamic variables of transmural pressure, pulse pressure and flow rate influence vessel size and radial strain. The physiologic implications of these findings are discussed.     

Three-dimensional,unsteady simulation of alveolar respiration

A novel macroscopic gas transport model, derived from fundamental engineering principles, is used to simulate the three-dimensional, unsteady respiration process within the alveolar region of the lungs. The simulations, mimicking the single-breath technique for measuring the lung diffusing capacity for carbon-monoxide (CO), allow the prediction of the red blood cell (RBC) distribution effects on the lung diffusing capacity. Results, obtained through numerical simulations, unveil a strong relationship between the type of distribution and the lung diffusing capacity. Several RBC distributions are considered, namely: normal (random), uniform, center-cluster, and corner-cluster red cell distributions. A nondimensional correlation is obtained in terms of a geometric parameter characterizing the RBC distribution, and presented as a useful tool for predicting the RBC distribution effect on the lung diffusing capacity. The effect of red cell movement is not considered in the present study because CO does not equilibrate with capillary blood within the time spent by blood in the capillary. Hence, blood flow effect on CO diffusion is expected to be only marginal.     

Experiments on steady and oscillatory flows at moderate Reynolds numbers in a quasi-two-dimensional channel with a throat.

The study of steady and unsteady oscillatory static fluid pressures acting on the internal wall of a collapsible tube is essential for investigation of the complicated behavior observed when a flow is conveyed inside a tube. To examine the validity of two one-dimensional nonsteady theoretical flow models, this paper presents basic experimental observations of flow separation and reattachment and measured data on the static pressure distributions of the flow in a quasi-two-dimensional channel with a throat, together with information on the corresponding shape of the wall deflection and motion. For combinations of moderate Reynolds numbers and angles of the divergent segment of the channel, a smooth flow is separated from the wall downstream of the minimum cross section and reattached to the wall farther downstream. The measured data are compared with numerical results calculated by the two flow models.     

Flow and Oxygen Transfer Characteristics of an Intravenous Membrane Oxygenator: A Computational Study

Spectral element computational simulations of the conservation of mass, momentum and species equations are performed to investigate the flow and oxygen transfer characteristics of an Intravenous Membrane Oxygenator (IMO). The simulations consider a three-dimensional IMO computational model consisting of equally-spaced fibers, an elastic balloon with non-permeable walls positioned longitudinally within the vena cava, and a Newtonian and time-dependent incompressible flow. Flow characteristics and oxygen transfer parameters are determined for operating conditions of a stationary and a pulsating balloon. For the stationary balloon configuration the flow is two-dimensional, parallel, laminar and without secondary flows for the Reynolds number range of 5.7-455.2. Evaluations of the oxygen transfer characteristics for the stationary balloon indicate that the main transport mechanisms are diffusion and convection in the crosswise and streamwise directions, respectively. Additionally, evaluations of oxygen transfer rates and Sherwood numbers in this Reynolds number range indicate that the oxygen transfer rate reaches an asymptotic limit at relatively moderate Reynolds numbers. For the pulsating balloon, flow characteristic results demonstrate the existence of a strong secondary flow around the fiber, and between the balloon and the fiber. This secondary flow induces oscillatory crosswise and streamwise velocities and a seemingly random spanwise flow which enhances the flow mixing as well as the transport of oxygen from the fiber surface to the bulk flow.     

Dynamics of vesicles in a wall-bounded shear flow

We report a detailed study of the behavior (shapes, experienced forces, velocities) of giant lipid vesicles subjected to a shear flow close to a wall. Vesicle buoyancy, size, and reduced volume were separately varied. We show that vesicles are deformed by the flow and exhibit a tank-treading motion with steady orientation. Their shapes are characterized by two nondimensional parameters: the reduced volume and the ratio of the shear stress with the hydrostatic pressure. We confirm the existence of a force, able to lift away nonspherical buoyant vesicles from the substrate. We give the functional variation and the value of this lift force (up to 150 pN in our experimental conditions) as a function of the relevant physical parameters: vesicle-substrate distance, wall shear rate, viscosity of the solution, vesicle size, and reduced volume. Circulating deformable cells disclosing a nonspherical shape also experience this force of viscous origin, which contributes to take them away from the endothelium and should be taken into account in studies on cell adhesion in flow chambers, where cells membrane and the adhesive substrate are in relative motion. Finally, the kinematics of vesicles along the flow direction can be described in a first approximation with a model of rigid spheres.     

Flow and Oxygen Transfer Characteristics of an Intravenous Membrane Oxygenator: A Computational Study

Abstract

Spectral element computational simulations of the conservation of mass, momentum and species equations are performed to investigate the flow and oxygen transfer characteristics of an Intravenous Membrane Oxygenator (IMO). The simulations consider a three-dimensional IMO computational model consisting of equally-spaced fibers, an elastic balloon with non-permeable walls positioned longitudinally within the vena cava, and a Newtonian and time-dependent incompressible flow. Flow characteristics and oxygen transfer parameters are determined for operating conditions of a stationary and a pulsating balloon. For the stationary balloon configuration the flow is two-dimensional, parallel, laminar and without secondary flows for the Reynolds number range of 5.7-455.2. Evaluations of the oxygen transfer characteristics for the stationary balloon indicate that the main transport mechanisms are diffusion and convection in the crosswise and streamwise directions, respectively. Additionally, evaluations of oxygen transfer rates and Sherwood numbers in this Reynolds number range indicate that the oxygen transfer rate reaches an asymptotic limit at relatively moderate Reynolds numbers. For the pulsating balloon, flow characteristic results demonstrate the existence of a strong secondary flow around the fiber, and between the balloon and the fiber. This secondary flow induces oscillatory crosswise and streamwise velocities and a seemingly random spanwise flow which enhances the flow mixing as well as the transport of oxygen from the fiber surface to the bulk flow.     

Effects of elastic property of the wall on flow characteristics through arterial stenoses

Hemodynamic characteristics of blood flow through arterial stenoses are numerically investigated in this work. The blood is assumed as a Newtonian fluid and the pulsatile nature of flow is modeled by using measured values of the flowrate and pressure for the canine femoral artery. An isotropic elastic and incompressible material is assumed for the wall at each axial section, but a non-uniform distribution of the shear modulus in axial direction is used to model the high stiffness of the wall at the stenosis location. Full Navier equations for a thick wall are used as the governing equations for the wall displacements. A continuous grid extending over the flow field and the wall is considered and governing equations are transformed for use in the computational domain. Discretized forms of the transformed wall and flow equations, which are coupled through the boundary conditions at their interface, are obtained by control volume method and simultaneously solved using the well-known SIMPLER algorithm. To study the effects of wall deformability, solutions are obtained for both rigid and elastic walls. The results indicate that deformability of the wall causes an increase in the time average of pressure drop, but a decrease in the maximum wall shear stress. Displacement and stress distributions in the wall are presented.