Computation of Stokes Flow in a Channel with a Collapsible Segment

The numerical solution of Stokes flow in two-dimensional channel in which a segment of one wall is formed by an elastic membrane under longitudinal tension and the remaining channel boundary is rigid is considered. This model problem is being used to gain an understanding of the complex interactions that occurs between the fluid flow and the wall mechanics when fluid flows through a collapsible tube, examples of which are widespread in physiology. Previous work by Pedley considered a similar system using lubrication theory in which the wall slopes are assumed small. The results showed that as the longitudinal wall tension is reduce, the downstream end of the collapsible segment becomes ever steeper, thus violating the assumptions. Here, lubrication theory is abandoned and a numerical solution of the full governing equations, including the complete expression for wall curvature, is sought using an iterative scheme. The effect of the variation in wall tension due to the fluid shear stresses at the compliant boundary is also included.Results are presented for a range of transmural (internal minus external) pressures and wall tensions. It is found, however, that as the wall tension is reduced, the iterative scheme considered fails to converge. This similar behaviour to that seen by Silliman & Scriven in viscous free-surface flows. Possible reasons for this breakdown together with alternative solution strategies are discussed.     

Effective equations describing the flow of a viscous incompressible fluid through a long elastic tubeEquations efficaces décrivant l'écoulement d'un fluide visqueux incompressible dans un tuyau long et élastique

We study the flow of a viscous incompressible fluid through a long and narrow elastic tube whose walls are modeled by the Navier equations for a curved, linearly elastic membrane. The flow is governed by a given small time dependent pressure drop between the inlet and the outlet boundary, giving rise to creeping flow modeled by the Stokes equations. By employing asymptotic analysis in thin, elastic, domains we obtain the reduced equations which correspond to a Biot type viscoelastic equation for the effective pressure and the effective displacement. The approximation is rigorously justified by obtaining the error estimates for the velocity, pressure and displacement. Applications of the model problem include blood flow in small arteries. We recover the well-known Law of Laplace and provide a new, improved model when shear modulus of the vessel wall is not negligible. To cite this article: S. ?ani?, A. Mikeli?, C. R. Mecanique 330 (2002) 661–666.     

Modelling the flow of power law fluids in a packed bed using a volume-averaged equation of motion

The development of a theoretical model for the prediction of velocity and pressure drop for the flow of a viscous power law fluid through a bed packed with uniform spherical particles is presented. The model is developed by volume averaging the equation of motion. A porous microstructure model based on a cell model is used. Numerical solution of the resulting equation is effected using a penalty Galerkin finite element method. Experimental pressure drop values for dilute solutions of carboxymethylcellulose flowing in narrow tubes packed with uniformly sized spherical particles are compared to theoretical predictions over a range of operating conditions. Overall agreement between experimental and theoretical values is within 15%. The extra pressure drop due to the presence of the wall is incorporated directly into the model through the application of the no-slip boundary condition at the container wall. The extra pressure drop reaches a maximum of about 10% of the bed pressure drop without wall effect. The wall effect increases as the ratio of tube diameter to particle diameter decreases, as the Reynolds number decreases and as the power law index increases.     

Effect of elastic wall motion on oscillatory flow of micropolar fluid in an annular tube

Summary Oscillatory flow of a micropolar fluid in an annular tube is investigated. The outer wall of the tube is taken to be elastic and the variation in the diameter of the elastic wall due to pulsatile nature of pressure gradient is assumed to be small. The wall motion is governed by a tube law. The nonlinear equations governing the fluid flow and the tube law are solved using perturbation analysis. The steady-streaming phenomenon due to the interaction of convected inertia with viscous effects is studied. The analysis, is carried out for zero mean flow rate. It presents the effects of the elastic nature of the wall combined with micropolar fluid parameters on the mean pressure gradient and wall shear stress for different catheter sizes and frequency parameters. It is found that the effect of micropolarity is of considerable importance for small steady-streaming Reynolds number. Also, it is observed that the relationship between mean pressure gradient and the flow rate depends on the amplitude of the diameter variation, flow rate waveforms and the phase difference between them.     

Numerical simulation of a turbulent flow with droplets injection in annular heated air tube using the Reynolds stress model

This work presents computational fluid dynamics (CFD) simulations of single-phase and two-phase flow. The droplets are injected in annular heated air tube. The numerical simulation is performed by using a commercial CFD code witch uses the finite-volume method to discretize the equations of fluid flow. The Reynolds-averaged Navier–Stokes equations with Reynolds stress model were used in the computation. The governing equations are solved by using a SIMPLE algorithm to treat the pressure terms in the momentum equations. The results of prediction are compared with the experimental data.     

Peristaltic transport of rheological fluid: model for movement of food bolus through esophagus

Fluid mechanical peristaltic transport through esophagus is studied in the paper.A mathematical model has been developed to study the peristaltic transport of a rheological fluid for arbitrary wave shapes and tube lengths.The Ostwald-de Waele power law of a viscous fluid is considered here to depict the non-Newtonian behaviour of the fluid.The model is formulated and analyzed specifically to explore some important information concerning the movement of food bolus through esophagus.The analysis is carried out by using the lubrication theory.The study is particularly suitable for the cases where the Reynolds number is small.The esophagus is treated as a circular tube through which the transport of food bolus takes place by periodic contraction of the esophageal wall.Variation of different variables concerned with the transport phenomena such as pressure,flow velocities,particle trajectory,and reflux is investigated for a single wave as well as a train of periodic peristaltic waves.The locally variable pressure is seen to be highly sensitive to the flow index "n".The study clearly shows that continuous fluid transport for Newtonian/rheological fluids by wave train propagation is more effective than widely spaced single wave propagation in the case of peristaltic movement of food bolus in the esophagus.     

膨胀性非牛顿流体多孔介质流动的模拟分析

在表征体元尺度采用格子Boltzmann方法分析膨胀性非牛顿流体在多孔介质中的流动,基于二阶矩模型在演化方程中引入表征介质阻力的作用力项,求解描述渗流模型的广义Navier-Stokes方程．采用局部法计算形变速率张量,通过循环迭代得到非牛顿粘度和松弛时间．对多孔介质的Poiseuille流动进行分析,通过比较发现结果与孔隙尺度的解析解十分吻合,并且收敛较快,表明方法合理有效．分析了渗透率和幂律指数对速度和压力降的影响,研究结果表明,膨胀性流体的多孔介质流动不符合达西规律,压力降的增加幅度小于渗透率的减小幅度．当无量纲渗透率Da小于10-5时,流道中的速度呈现均匀分布,并且速度分布随着幂律指数的减小趋于平滑．压力降随着幂律指数的增加而增加,Da越大幂律指数对压力降的影响越明显．     

A unified approach to simulation and upscaling of single‐phase flow through vuggy carbonates

Numerical modeling of flow through vuggy porous media, mainly vuggy carbonates, is a challenging endeavor. Firstly, because the presence of vugs can significantly alter the effective porosity and permeability of the medium. Secondly, because of the co‐existence of porous and free flow regions within the medium and the uncertainties in defining the exact boundary between them. Traditionally, such heterogeneous systems are modeled by the coupled Darcy–Stokes equations. However, numerical modeling of flow through vuggy porous media using coupled Darcy–Stokes equations poses several numerical challenges particularly with respect to specification of correct interface condition between the porous and free‐flow regions. Hence, an alternative method, a more unified approach for modeling flows through vuggy porous media, the Stokes–Brinkman model, is proposed here. It is a single equation model with variable coefficients, which can be used for both porous and free‐flow regions. This also makes the requirement for interface condition redundant. Thus, there is an obvious benefit of using the Stokes–Brinkman equation, which can be reduced to Stokes or Darcy equation by the appropriate choice of parameters. At the same time, the Stokes–Brinkman equation provides a smooth transition between porous and free‐flow region, thereby taking care of the associated uncertainties. A numerical treatment for upscaling Stokes–Brinkman model is presented as an approach to use Stokes–Brinkman model for multi‐phase flow. Numerical upscaling methodology is validated by analyzing the error norm for numerical pressure convergence. Stokes–Brinkman single equation model is tested on a series of realistic data sets, and the results are compared with traditional coupled Darcy–Stokes model. Copyright © 2011 John Wiley & Sons, Ltd.     

Steady flow through a slender tapered vessel with a partially permeable wall and closed end

We study the flow of a viscous fluid through a slender tapered tube whose radius may reduce to zero. The vessel is closed at the end, so that the flow is made possible owing to the fact that a portion of the tube wall is permeable. The smallness of the tube aspect ratio is exploited using an upscaling technique leading to a degenerate differential equation for pressure. Solutions are found either in explicit form or as power series expansions. This class of flows may represent, though in a largely approximated way, the blood flow though a coronary artery.     

Motion through spherical droplet with non-homogenous porous layer in spherical container

The problem of the creeping flow through a spherical droplet with a nonhomogenous porous layer in a spherical container has been studied analytically. Darcy's model for the flow inside the porous annular region and the Stokes equation for the flow inside the spherical cavity and container are used to analyze the flow. The drag force is exerted on the porous spherical particles enclosing a cavity, and the hydrodynamic permeability of the spherical droplet with a non-homogeneous porous layer is calculated. Emphasis is placed on the spatially varying permeability of a porous medium, which is not covered in all the previous works related to spherical containers. The variation of hydrodynamic permeability and the wall effect with respect to various flow parameters are presented and discussed graphically. The streamlines are presented to discuss the kinematics of the flow. Some previous results for hydrodynamic permeability and drag forces have been verified as special limiting cases.     

New Correlative Models to Improve Prediction of Fracture Permeability and Inertial Resistance Coefficient

Presence of fracture roughness and occurrence of nonlinear flow complicate fluid flow through rock fractures. This paper presents a qualitative and quantitative study on the effects of fracture wall surface roughness on flow behavior using direct flow simulation on artificial fractures. Previous studies have highlighted the importance of roughness on linear and nonlinear flow through rock fractures. Therefore, considering fracture roughness to propose models for the linear and nonlinear flow parameters seems to be necessary. In the current report, lattice Boltzmann method is used to numerically simulate fluid flow through different fracture realizations. Flow simulations are conducted over a wide range of pressure gradients through each fracture. It is observed that creeping flow at lower pressure gradients can be described using Darcy’s law, while transition to inertial flow occurs at higher pressure gradients. By detecting the onset of inertial flow and regression analysis on the simulation results with Forchheimer equation, inertial resistance coefficients are determined for each fracture. Fracture permeability values are also determined from Darcy flow as well. According to simulation results through different fractures, two parametric expressions are proposed for permeability and inertial resistance coefficient. The proposed models are validated using 3D numerical simulations and experimental results. The results obtained from these two proposed models are further compared with those obtained from the conventional models. The calculated average absolute relative errors and correlation coefficients indicate that the proposed models, despite their simplicity, present acceptable outcomes; the models are also more accurate compared to the available methods in the literature.     

Effect of yield stress on the flow of a Casson fluid in a homogeneous porous medium bounded by a circular tube

The effect of yield stress on the flow characteristics of a Casson fluid in a homogeneous porous medium bounded by a circular tube is investigated by employing the Brinkman model to account for the Darcy resistance offered by the porous medium. The non-linear coupled implicit system of differential equations governing the flow is first transformed into suitable integral equations and are solved numerically. Analytical solution is obtained for a Newtonian fluid in the case of constant permeability, and the numerical solution is verified with that of the analytic solution. The effect of yield stress of the fluid and permeability of the porous medium on shear stress and velocity distributions, plug flow radius and flow rate are examined. The minimum pressure gradient required to start the flow is found to be independent of the permeability of the porous medium and is equal to the yield stress of the fluid.     

Flow of a viscous fluid past a heterogeneous porous sphere at low Reynolds numbers

A flow past a heterogeneous porous sphere is investigated by using the perturbation theory. The flow through the sphere is divided into two zones, which are fully saturated with the viscous fluid, and the flow in these zones is governed by the Brinkman equation. The space outside the sphere, where a clear fluid flows, is also divided into two zones: the Navier–Stokes zone and the Oseen flow zone. The solutions on the interface inside the sphere are matched with the condition proposed by Merrikh and Mohammad. The stream function in the Navier–Stokes zone is matched with that on the sphere surface by the condition proposed by Ochoa-Tapia and Whitaker. It is found that the drag on the spherical shell decreases as the permeability toward the sphere boundary increases.     

Effective Permeability of a Porous Medium with Spherical and Spheroidal Vug and Fracture Inclusions

Vugs and fractures are common features of carbonate formations. The presence of vugs and fractures in porous media can significantly affect pressure and flow behavior of a fluid. A vug is a cavity (usually a void space, occasionally filled with sediments), and its pore volume is much larger than the intergranular pore volume. Fractures occur in almost all geological formations to some extent. The fluid flow in vugs and fractures at the microscopic level does not obey Darcy’s law; rather, it is governed by Stokes flow (sometimes is also called Stokes’ law). In this paper, analytical solutions are derived for the fluid flow in porous media with spherical- and spheroidal-shaped vug and/or fracture inclusions. The coupling of Stokes flow and Darcy’s law is implemented through a no-jump condition on normal velocities, a jump condition on pressures, and generalized Beavers–Joseph–Saffman condition on the interface of the matrix and vug or fracture. The spheroidal geometry is used because of its flexibility to represent many different geometrical shapes. A spheroid reduces to a sphere when the focal length of the spheroid approaches zero. A prolate spheroid degenerates to a long rod to represent the connected vug geometry (a tunnel geometry) when the focal length of the spheroid approaches infinity. An oblate spheroid degenerates to a flat spheroidal disk to represent the fracture geometry. Once the pressure field in a single vug or fracture and in the matrix domains is obtained, the equivalent permeability of the vug with the matrix or the fracture with matrix can be determined. Using the effective medium theory, the effective permeability of the vug–matrix or fracture–matrix ensemble domain can be determined. The effect of the volume fraction and geometrical properties of vugs, such as the aspect ratio and spatial distribution, in the matrix is also investigated. It is shown that the higher volume fraction of the vugs or fractures enhances the effective permeability of the system. For a fixed-volume fraction, highly elongated vugs or fractures significantly increase the effective permeability compared with shorter vugs or fractures. A set of disconnected vugs or fractures yields lower effective permeability compared with a single vug or fracture of the same volume fraction.     

Shear flow over a porous particle

Linear axisymmetric Stokes flow over a porous spherical particle is investigated. An exact analytic solution for the fluid velocity components and the pressure inside and outside the porous particle is obtained. The solution is generalized to include the cases of arbitrary three-dimensional linear shear flow as well as translational-shear Stokes flow. As the permeability of the particle tends to zero, the solutions obtained go over into the corresponding solutions for an impermeable particle. The problem of translational Stokes flow around a spherical drop (in the limit a gas bubble or an impermeable sphere) was considered, for example, in [1,2]. A solution of the problem of translational Stokes flow over a porous spherical particle was given in [3]. Linear shear-strain Stokes flow over a spherical drop was investigated in [2].Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 3, pp. 113–120, May–June, 1995.     

NUMERICAL SIMULATION OF STEADY FLOW IN A COMPLIANT TUBE OR CHANNEL WITH TAPERED WALL THICKNESS

Flow through compliant tubes with linear taper in wall thickness is numerically simulated by finite element analysis. Two models are examined: a compliant channel and an axisymmetric tube. For verification of the numerical method, flow through a compliant stenotic vessel is simulated and compared to existing experimental data. Steady two-dimensional flow in a collapsible channel with initial tension is also simulated and the results are compared with numerical solutions from the literature. Computational results for an axisymmetric tube show that as cross-sectional area falls with a reduction in downstream pressure, flow rate increases and reaches a maximum when the speed index (mean velocity divided by wave speed) is near unity at the point of minimum cross-sectional area, indicative of wave-speed flow limitation or “choking” (flow speed equals wave speed) in previous one-dimensional studies. For further reductions in downstream pressure, the flow rate decreases. Cross-sectional narrowing is significant but localized. For the particular wall and fluid properties used in these simulations, the area throat is located near the downstream end when the ratio of downstream-to-upstream wall thickness is 2; as wall taper is increased to 3, the constriction moves to the upstream end of the tube. In the planar two-dimensional channel, area reduction and flow limitation are also observed when outlet pressure is decreased. In contrast to the axisymmetric case, however, the elastic wall in the two-dimensional channel forms a smooth concave surface with the area throat located near the mid-point of the elastic wall. Though flow rate reaches a maximum and then falls, the flow does not appear to be choked.     

Peristaltic flow of MHD Jeffrey fluid through finite length cylindrical tube

The present investigation studies the peristaltic flow of the Jeffrey fluid through a tube of finite length. The fluid is electrically conducting in the presence of an applied magnetic field. Analysis is carried out under the assumption of long wavelength and low Reynolds number approximations. Expressions of the pressure gradient, volume flow rate, average volume flow rate, and local wall shear stress are obtained. The effects of relaxation time, retardation time, Hartman number on pressure, local wall she...     

High-Velocity Laminar and Turbulent Flow in Porous Media

We model high-velocity flow in porous media with the multiple scale homogenization technique and basic fluid mechanics. Momentum and mechanical energy theorems are derived. In idealized porous media inviscid irrotational flow in the pores and wall boundary layers give a pressure loss with a power of 3/2 in average velocity. This model has support from flow in simple model media. In complex media the flow separates from the solid surface. Pressure loss effects of flow separation, wall and free shear layers, pressure drag, flow tube velocity and developing flow are discussed by using phenomenological arguments. We propose that the square pressure loss in the laminar Forchheimer equation is caused by development of strong localized dissipation zones around flow separation, that is, in the viscous boundary layer in triple decks. For turbulent flow, the resulting pressure loss due to average dissipation is a power 2 term in velocity.     

Numerical modelling of shear-dependent mass transfer in large arteries

A numerical scheme for the simulation of blood flow and transport processes in large arteries is presented. Blood flow is described by the unsteady 3D incompressible Navier–Stokes equations for Newtonian fluids; solute transport is modelled by the advection–diffusion equation. The resistance of the arterial wall to transmural transport is described by a shear-dependent wall permeability model. The finite element formulation of the Navier–Stokes equations is based on an operator-splitting method and implicit time discretization. The streamline upwind/Petrov–Galerkin (SUPG) method is applied for stabilization of the advective terms in the transport equation and in the flow equations. A numerical simulation is carried out for pulsatile mass transport in a 3D arterial bend to demonstrate the influence of arterial flow patterns on wall permeability characteristics and transmural mass transfer. The main result is a substantial wall flux reduction at the inner side of the curved region. © 1997 John Wiley & Sons, Ltd.     

A Two-Layer Mathematical Model of Blood Flow in Porous Constricted Blood Vessels

In this paper, we discussed a mathematical model for two-layered non-Newtonian blood flow through porous constricted blood vessels. The core region of blood flow contains the suspension of erythrocytes as non-Newtonian Casson fluid and the peripheral region contains the plasma flow as Newtonian fluid. The wall of porous constricted blood vessel configured as thin transition Brinkman layer over layered by Darcy region. The boundary of fluid layer is defined as stress jump condition of Ocha-Tapiya and Beavers–Joseph. In this paper, we obtained an analytic expression for velocity, flow rate, wall shear stress. The effect of permeability, plasma layer thickness, yield stress and shape of the constriction on velocity in core & peripheral region, wall shear stress and flow rate is discussed graphically. This is found throughout the discussion that permeability and plasma layer thickness have accountable effect on various flow parameters which gives an important observation for diseased blood vessels.