A two-fluid model for pulsatile flow in catheterized blood vessels

The pulsatile flow of blood through a catheterized artery is analyzed, assuming the blood as a two-fluid model with the suspension of all the erythrocytes in the core region as a Casson fluid and the peripheral region of plasma as a Newtonian fluid. The resulting non-linear implicit system of partial differential equations is solved using perturbation method. The expressions for shear stress, velocity, flow rate, wall shear stress and longitudinal impedance are obtained. The variations of these flow quantities with yield stress, catheter radius ratio, amplitude, pulsatile Reynolds number ratio and peripheral layer thickness are discussed. It is observed that the velocity distribution and flow rate decrease, while, the wall shear, width of the plug flow region and longitudinal impedance increase when the yield stress increases. It is also found that the velocity increases, but, the longitudinal impedance decreases when the thickness of the peripheral layer increases. The wall shear stress decreases non-linearly, while, the longitudinal impedance increases non-linearly when the catheter radius ratio increases. The estimates of the increase in the longitudinal impedance are considerably lower for the present two-fluid model than those of the single-fluid model.     

Pulsatile flow of Herschel-Bulkley fluid through stenosed arteries—A mathematical model

In this paper, the pulsatile flow of blood through stenosed artery is studied. The effects of pulsatility, stenosis and non-Newtonian behavior of blood, assuming the blood to be represented by Herschel-Bulkley fluid, are simultaneously considered. A perturbation method is used to analyze the flow assuming the thickness of plug core region to be non-uniform changing with axial distance. An expression for the variation of plug core radius with time and axial distance is obtained. The variation of pressure gradient with steady flow rate is given. Also the variation of wall shear stress distribution as well as resistance to flow with axial distance for different values of time and for different values of yield stress is given and the results analyzed.     

Two-fluid non-linear model for flow in catheterized blood vessels

Blood flow through a catheterized artery is analyzed, assuming the flow is steady and blood is treated as a two-fluid model with the suspension of all the erythrocytes in the core region as a Casson fluid and the plasma in the peripheral region as a Newtonian fluid. The expressions for velocity, flow rate, wall shear stress and frictional resistance are obtained. The variations of these flow quantities with yield stress, catheter radius ratio and peripheral layer thickness are discussed. It is noticed that the velocity and flow rate decrease while the wall shear stress and resistance to flow increase when the yield stress or the catheter radius ratio increases while all the other parameters were held fixed. It is found that the velocity and flow rate increase while the wall shear stress and frictional resistance decrease with the increase of the peripheral layer thickness. The estimates of the increase in the frictional resistance are significantly very small for the present two-fluid model than those of the single-fluid Casson model.     

Two-phase flow pattern identification using continuous hidden Markov model

This paper presents a new method for identifying two-phase flow regimes from the instantaneous local fluid phase signals using continuous hidden Markov model (CHMM). CHMM is known to be a very strong pattern identifier. Air–water two-phase flows were realized in a transparent vertical tube. The tube length was 2 m, and its inner diameter was 19 mm. The instantaneous local fluid phase signals were collected using a single step index multimode optical fiber probe located at the center and mid-length of the tube. Signal features required in CHMM implementation were extracted using an innovative method. Various aspects of hidden Markov modeling and their effects on the results were studied. The flow pattern results are in very good agreement with photographs of the flow captured during the experiments. In sum, the results show that hidden Markov model has a good potential in identifying two-phase flow patterns.     

An analysis model of pulsatile blood flow in arteries

IntroductionTheperiodicallypulsatilebloodflowinthearterycausesthecircumferentialandaxialmotionoftheelasticbloodvesselandinturntheoscillationofthevesselaffectsthatofthebloodflow .Womersley[1]resolvedsuccessfullythisfluid_solidcouplingproblembysolvingbothlinearNavier_Stokesequationsandthemotionequationsofthethin_walledelastictubeandgainedtheexpressionsofthebloodflowvelocitiesandthevasculardisplacements.Histheoryhasbeenthebasisforthequantitativeanalysisoftherelationshipofthearterialstructureandi…     

Numerical studies of blood flow in healthy,stenosed, and stented carotid arteries

Numerical analysis of pulsatile blood flow in healthy, stenosed, and stented carotid arteries is performed with the aim of identifying hemodynamic factors in the initiation, growth, and the potential of leading to severe occlusions of a diseased artery. The Immersed Finite Element Method is adopted for this study to conveniently incorporate various geometrical shapes of arteries without remeshing. Our computational results provide detailed quantitative analysis on the blood flow pattern, wall shear stress, particle residence time, and oscillatory shear index. The analysis of these parameters leads to a better understanding of blood clot formation and its localization in a stenosed and a stented carotid artery. A healthy artery is also studied to establish a baseline comparison. This analysis will assist in developing treatments for diseased arteries and novel stent designs to reduce restenosis. Copyright © 2008 John Wiley & Sons, Ltd.     

Macro-to-microchannel transition in two-phase flow: Part 1 - Two-phase flow patterns and film thickness measurements

The classification of macroscale, mesoscale and microscale channels with respect to two-phase processes is still an open question. The main objective of this study focuses on investigating the macro-to-microscale transition during flow boiling in small scale channels of three different sizes with three different refrigerants over a range of saturation conditions to investigate the effects of channel confinement on two-phase flow patterns and liquid film stratification in a single circular horizontal channel (Part 2 covers the flow boiling heat transfer and critical heat flux). This paper presents the experimental two-phase flow pattern transition data together with a top/bottom liquid film thickness comparison for refrigerants R134a, R236fa and R245fa during flow boiling in small channels of 1.03, 2.20 and 3.04 mm diameter. Based on this work, an improved flow pattern map has been proposed by determining the flow patterns transitions existing under different conditions including the transition to macroscale slug/plug flow at a confinement number of Co ≈ 0.3-0.4. From the top/bottom liquid film thickness comparison results, it was observed that the gravity forces are fully suppressed and overcome by the surface tension and shear forces when the confinement number approaches 1, Co ≈ 1. Thus, as a new approximate rule, the lower threshold of macroscale flow is Co = 0.3-0.4 while the upper threshold of symmetric microscale flow is Co ≈ 1 with a transition (or mesoscale) region in-between.     

Two-phase flow and pool boiling heat transfer in microgravity

Researches on two-phase flow and pool boiling heat transfer in microgravity, which included ground-based tests, flight experiments, and theoretical analyses, were conducted in the National Microgravity Laboratory/CAS. A semi-theoretical Weber number model was proposed to predict the slug-to-annular flow transition of two-phase gas–liquid flows in microgravity, while the influence of the initial bubble size on the bubble-to-slug flow transition was investigated numerically using the Monte Carlo method. Two-phase flow pattern maps in microgravity were obtained in the experiments both aboard the Russian space station Mir and aboard IL-76 reduced gravity airplane. Mini-scale modeling was also used to simulate the behavior of microgravity two-phase flow on the ground. Pressure drops of two-phase flow in microgravity were also measured experimentally and correlated successfully based on its characteristics. Two space experiments on pool boiling phenomena in microgravity were performed aboard the Chinese recoverable satellites. Steady pool boiling of R113 on a thin wire with a temperature-controlled heating method was studied aboard RS-22, while quasi-steady pool boiling of FC-72 on a plate was studied aboard SJ-8. Ground-based experiments were also performed both in normal gravity and in short-term microgravity in the drop tower Beijing. Only slight enhancement of heat transfer was observed in the wire case, while enhancement in low heat flux and deterioration in high heat flux were observed in the plate case. Lateral motions of vapor bubbles were observed before their departure in microgravity. The relationship between bubble behavior and heat transfer on plate was analyzed. A semi-theoretical model was also proposed for predicting the bubble departure diameter during pool boiling on wires. The results obtained here are intended to become a powerful aid for further investigation in the present discipline and development of two-phase systems for space applications.     

Numerical simulation of axisymmetric turbulent flow in combustors and diffusers

Numerical studies of turbulent flow in an axisymmetric 45° expansion combustor and bifurcated diffuser are presented. The Navier-Stokes equations incorporating a k–? model were solved in a non-orthogonal curvillinear co-ordinate system. A zonal grid method, wherein the flow field was divided into several subsections, was developed. This approach permitted different computational schemes to be used in the various zones. In addition, grid generation was made a more simple task. However, treatment of the zonal boundaries required special handling. Boundary overlap and interpolating techniques were used and an adjustment of the flow variables was required to assure conservation of mass flux. Three finite differencing methods—hybrid, quadratic upwind and skew upwind—were used to represent the convection terms. Results were compared with existing experimental data. In general, good agreement between predicted and measured values was obtained.     

Efficient implementation of the proper outlet flow conditions in blood flow simulations through asymmetric arterial bifurcations

We present an efficient implementation of the proper (in vivo) outlet boundary conditions in detailed, three‐dimensional (3D) and time‐periodic simulations of blood flow through arteries. This is achieved through the intermediate use of an approximate ‘simulant’ model of the outlet pressure/flow relationship corresponding to the full 3D and time‐dependent numerical simulation. This model allows us to efficiently couple the 3D outlet pressure/flow conditions to the equivalent relations due to the downstream arterial network, as obtained from a one‐dimensional approximate model in the form of Fourier frequency impedance coefficients. An adjustable time‐periodic function correction term in the simulant model requires input from the full 3D model that has to run iteratively until convergence. The advantage of the proposed numerical scheme is that it decouples the upstream detailed simulation from the downstream approximate network model offering exceptional versatility. This approach is demonstrated here in a series of detailed 3D simulations of blood flow, performed using the commercial software FLUENT?, through an asymmetric arterial bifurcation. Two cases are considered: first a healthy system patterned after the left main coronary arterial bifurcation, and second a diseased case where an occlusion has developed in one of the daughter vessels, resulting in strengthening the asymmetry of the bifurcation. Rapid convergence of the iterative process was achieved in both cases. Subtle changes occur in the shear patterns of the daughter vessels, whereas the flow distribution is quite different. In the presence of a stenosis additional regions of low shear develop due to inertial effects. Copyright © 2010 John Wiley & Sons, Ltd.     

Two-phase flow in heterogeneous porous media III: Laboratory experiments for flow parallel to a stratified system

Two-phase flow in stratified porous media is a problem of central importance in the study of oil recovery processes. In general, these flows are parallel to the stratifications, and it is this type of flow that we have investigated experimentally and theoretically in this study. The experiments were performed with a two-layer model of a stratified porous medium. The individual strata were composed of Aerolith-10, an artificial: sintered porous medium, and Berea sandstone, a natural porous medium reputed to be relatively homogeneous. Waterflooding experiments were performed in which the saturation field was measured by gamma-ray absorption. Data were obtained at 150 points distributed evenly over a flow domain of 0.1 × 0.6 m. The slabs of Aerolith-10 and Berea sandstone were of equal thickness, i.e. 5 centimeters thick. An intensive experimental study was carried out in order to accurately characterize the individual strata; however, this effort was hampered by both local heterogeneities and large-scale heterogeneities.The theoretical analysis of the waterflooding experiments was based on the method of large-scale averaging and the large-scale closure problem. The latter provides a precise method of discussing the crossflow phenomena, and it illustrates exactly how the crossflow influences the theoretical prediction of the large-scale permeability tensor. The theoretical analysis was restricted to the quasi-static theory of Quintard and Whitaker (1988), however, the dynamic effects described in Part I (Quintard and Whitaker 1990a) are discussed in terms of their influence on the crossflow.Roman Letters A interfacial area between the -region and the -region contained within V, m² - a vector that maps onto , m - b vector that maps onto , m - b vector that maps onto , m - B second order tensor that maps onto , m² - C second order tensor that maps onto , m² - E energy of the gamma emitter, keV - f fractional flow of the -phase - g gravitational vector, m/s² - h characteristic length of the large-scale averaging volume, m - H height of the stratified porous medium , m - i unit base vector in the x-direction - K local volume-averaged single-phase permeability, m² - K - {K}, large-scale spatial deviation permeability - { K} large-scale volume-averaged single-phase permeability, m² - K ^* large-scale single-phase permeability, m² - K ^** equivalent large-scale single-phase permeability, m² - K local volume-averaged -phase permeability in the -region, m² - K local volume-averaged -phase permeability in the -region, m² - K - {K } , large-scale spatial deviation for the -phase permeability, m² - K ^* large-scale permeability for the -phase, m² - l thickness of the porous medium, m - l characteristic length for the -region, m - l characteristic length for the -region, m - L length of the experimental porous medium, m - characteristic length for large-scale averaged quantities, m - n outward unit normal vector for the -region - n outward unit normal vector for the -region - n unit normal vector pointing from the -region toward the -region (n = - n ) - N number of photons - p pressure in the -phase, N/m² - p ₀ reference pressure in the -phase, N/m² - local volume-averaged intrinsic phase average pressure in the -phase, N/m² - large-scale volume-averaged pressure of the -phase, N/m² - large-scale intrinsic phase average pressure in the capillary region of the -phase, N/m² - - , large-scale spatial deviation for the -phase pressure, N/m² - p_c , capillary pressure, N/m² - p _c capillary pressure in the -region, N/m² - p capillary pressure in the -region, N/m² - {p _c } ^c large-scale capillary pressure, N/m² - q -phase velocity at the entrance of the porous medium, m/s - q -phase velocity at the entrance of the porous medium, m/s - S_wi irreducible water saturation - S /, local volume-averaged saturation for the -phase - S _i initial saturation for the -phase - S _r residual saturation for the -phase - S ^* { }^*/}^*, large-scale average saturation for the -phase - S saturation for the -phase in the -region - S saturation for the -phase in the -region - t time, s - v -phase velocity vector, m/s - v local volume-averaged phase average velocity for the -phase, m/s - {v } large-scale averaged velocity for the -phase, m/s - v local volume-averaged phase average velocity for the -phase in the -region, m/s - v local volume-averaged phase average velocity for the -phase in the -region, m/s - v -{v } , large-scale spatial deviation for the -phase velocity, m/s - v -{v } , large-scale spatial deviation for the -phase velocity in the -region, m/s - v -{v } , large-scale spatial deviation for the -phase velocity in the -region, m/s - V large-scale averaging volume, m³ - y position vector relative to the centroid of the large-scale averaging volume, m - {y}^c large-scale average of y over the capillary region, m Greek Letters local porosity - local porosity in the -region - local porosity in the -region - local volume fraction for the -phase - local volume fraction for the -phase in the -region - local volume fraction for the -phase in the -region - {}^* { }^*+{ }^*, large-scale spatial average volume fraction - { }^* large-scale spatial average volume fraction for the -phase - mass density of the -phase, kg/m³ - mass density of the -phase, kg/m³ - viscosity of the -phase, N s/m² - viscosity of the -phase, Ns/m² - V /V , volume fraction of the -region ( + =1) - V /V , volume fraction of the -region ( + =1) - attenuation coefficient to gamma-rays, m^-1 - -      

明渠气二相流的双流体模型

从解决平均二相流基本方程的封闭问题出发,分析了水利工程中气水二相流的特点,介绍了据此提出的可用于水利工程稀疏气泡流计算的双流体模型,并对明渠气水二相流进行了算例验证。     

Performance modeling and analysis of blood flow in elastic arteries

Nomenclature　　τ　　wallshearstressγshearrateτy yieldstressηc Cassonviscosityktheconsistencyindexnnon_Newtonianindexτp shearstressofthepthelementωangularvelocityRvessel’sradiusCwavespeedM　　magneticparameter (Hartmannnumber)u,w　velocitycomponentinther_andz_directions,respectivelyP　　pressureα　　unsteadinessparameter k , R meanparametersTp relaxationtimeofthepthelementρ densityIntroductionTheimportancetoatherogenesisofarterialflowphenomenasuchasflowseparation ,recirculationands…     

Well‐balanced discontinuous Galerkin method and finite volume WENO scheme based on hydrostatic reconstruction for blood flow model in arteries

The blood flow model in arteries admits the steady state solutions, for which the flux gradient is nonzero, and is exactly balanced by the source term. In this paper, by means of hydrostatic reconstruction, we construct a high order discontinuous Galerkin method, which exactly preserves the dead‐man steady state, which is characterized by a discharge equal to zero (analogue to hydrostatic equilibrium). Moreover, the method maintains genuine high order of accuracy. Subsequently, we apply the key idea to finite volume weighted essentially non‐oscillatory schemes and obtain a well‐balanced finite volume weighted essentially non‐oscillatory scheme. Extensive numerical experiments are performed to verify the well‐balanced property, high order accuracy, as well as good resolution for smooth and discontinuous solutions.     

An improved semi-empirical friction model for gas-liquid two-phase flow in horizontal and near horizontal pipes

Pressure drop and liquid hold-up are two very important fluid flow parameters in design and control of multiphase flow pipelines. Friction factors play an important role in the accurate calculation of pressure drop. Various empirical and semi-empirical closure relations exist in the literature to calculate the liquid-wall, gas-wall and interfacial friction in two-phase pipe flow.However most of them are empirical correlations found under special experimental conditions. In this paper by modification of a friction model available in the literature, an improved semiempirical model is proposed. The proposed model is incorporated in the two-fluid correlations under equilibrium conditions and solved. Pressure gradient and velocity profiles are validated against experimental data. Using the improved model, the pressure gradient deviation from experiments diminishes by about 3%; the no-slip condition at the interface is satisfied and the velocity profile is predicted in better agreement with the experimental data.     

Analytical and numerical solutions for axisymmetric flow of nanofluid due to non-linearly stretching sheet

This article reports the laminar axisymmetric flow of nanofluid over a non-linearly stretching sheet. The model used for nanofluid contains the simultaneous effects of Brownian motion and thermophoretic diffusion of nanoparticles. The recently proposed boundary condition is considered which requires the mass flux of nanoparticles at the wall to be zero. Analytic solutions of the arising boundary value problem are obtained by optimal homotopy analysis method. Moreover the numerical solutions are computed by Keller–Box method. Both the solutions are found in excellent agreement. The behavior of Brownian motion on the fluid temperature and wall heat transfer rate is insignificant. Further the nanoparticle volume fraction distribution is found to be negative near the vicinity of the stretching sheet.     

磁场和Hall电流对狭窄动脉中血液流动的影响

A micropolar model for blood simulating magnetohydrodynamic flow through a horizontally nonsymmetric but vertically symmetric artery with a mild stenosis is presented. To estimate the effect of the stenosis shape, a suitable geometry has been considered such that the horizontal shape of the stenosis can easily be changed just by varying a parameter referred to as the shape parameter. Flow parameters, such as velocity, the resistance to flow （the resistance impedance）, the wall shear stress distribution in the stenotic region, and its magnitude at the maximum height of the stenosis （stenosis throat）, have been computed for different shape parameters, the Hartmann number and the Hall parameter. This shows that the resistance to flow decreases with the increasing values of the parameter determining the stenosis shape and the Hail parameter, while it increases with the increasing Hartmann number. The wall shear stress and the shearing stress on the wall at the maximum height of the stenosis possess an inverse characteristic to the resistance to flow with respect to any given value of the Hartmann number and the Hall parameter. Finally, the effect of the Hartmann number and the Hall parameter on the horizontal velocity is examined.     

Numerical modeling of slug flow initiation in a horizontal channels using a two-fluid model

This paper presents a methodology for modeling slug initiation and growth in horizontal ducts. Transient two-fluid equations are solved numerically using a class of high-resolution shock capturing methods. The advantage of this method is that slug formation and growth in a stratified regime can be calculated directly from the solutions to the flow field differential equations. In addition, by using high-resolution shock capturing methods that do not contain numerical diffusion, the discontinuity generated by slugging in the flow field can be modeled with good accuracy. The two-fluid model is shown to be well-posed mathematically only under certain conditions. Under these circumstances, the two-fluid model is capable of correctly predicting and modeling the flow physics. When ill-posed, an unbounded instability occurs in the flow field solution, and the instability amplitude increases exponentially with decreasing mesh sizes. This work shows that there are three zones associated with slug formation. In addition, long wavelength slugs are shown to initiate from short wavelength waves. These short waves are generated at the interface of the two phases by the Kelvin-Helmholtz hydrodynamic instability. The results obtained through numerical modeling show good agreement with experimental results.     

The micropolar fluid model for blood flow through a tapered artery with a stenosis

A micropolar model for axisymmetric blood flow through an axially nonsymmetreic but radially symmetric mild stenosis tapered artery is presented. To estimate the effect of the stenosis shape, a suitable geometry has been considered such that the axial shape of the stenosis can be changed easily just by varying a parameter （referred to as the shape parameter）. The model is also used to study the effect of the taper angle Ф. Flow parameters such as the velocity, the resistance to flow （the resistance impedance）, the wall shear stress distribution in the stenotic region and its magnitude at the maximum height of the stenosis （stenosis throat） have been computed for different values of the shape parameter n, the taper angle Ф, the coupling number N and the micropolar parameter m. It is shown that the resistance to flow decreases with increasing the shape parameter n and the micropolar parameter m while it increases with increasing the coupling number N. So, the magnitude of the resistance impedance is higher for a micropolar fluid than that for a Newtonian fluid model. Finally, the velocity profile, the wall shear stress distribution in the stenotic region and its magnitude at the maximum height of the stenosis are discussed for different values of the parameters involved on the problem.