Numerical simulation of pressure-driven startup laminar flows through a planar T-junction channel

A start-up flow of a viscous incompressible fluid in a T-junction channel is studied numerically. The flow starting from rest is driven by a constant pressure drops suddenly applied between the entries and exits of a planar T-junction channel. The Navier-Stokes equations in primitive variables are solved numerically using finite-volume techniques. Predicted variations with time of the volume flow rates and the flow patterns are presented for several values of pressure drops. It has been shown that a start-up flow can pass through different regimes (or different flow direction) before asymptotically reaching steady state distribution.     

A micropolar fluid model of blood flow through a tapered artery with a stenosis

A nonlinear two‐dimensional micropolar fluid model for blood flow in a tapered artery with a single stenosis is considered. This model takes into account blood rheology in which blood consists of microelements suspended in plasma. The classical Navier–Stokes theory is inadequate to describe the microrotations or particles' spin of such suspension in a viscous medium. The governing equations involving unsteady nonlinear partial differential equations are solved using a finite difference scheme. A quantitative analysis performed through numerical computation shows that the axial velocity profile and the flow rate decrease and the wall shear stress increases once the artery is narrower in the presence of the polar effect. Furthermore, the taper angle certainly bears the potential to influence the velocity and the flow characteristics to considerable extent. Copyright © 2010 John Wiley & Sons, Ltd.     

Power law fluid model for blood flow through a tapered artery with a stenosis

In the present paper, blood flow through a tapered artery with a stenosis is analyzed, assuming the flow is steady and blood is treated as non-Newtonian power law fluid model. Exact solution has been evaluated for velocity, resistance impedance, wall shear stress and shearing stress at the stenosis throat. The graphical results of different types of tapered arteries (i.e. converging tapering, diverging tapering, non-tapered artery) have been examined for different parameters of interest. Some special cases of the problem are also presented.     

Convergence of gauge method for incompressible flow

A new formulation, a gauge formulation of the incompressible Navier-Stokes equations in terms of an auxiliary field and a gauge variable , , was proposed recently by E and Liu. This paper provides a theoretical analysis of their formulation and verifies the computational advantages. We discuss the implicit gauge method, which uses backward Euler or Crank-Nicolson in time discretization. However, the boundary conditions for the auxiliary field are implemented explicitly (vertical extrapolation). The resulting momentum equation is decoupled from the kinematic equation, and the computational cost is reduced to solving a standard heat and Poisson equation. Moreover, such explicit boundary conditions for the auxiliary field will be shown to be unconditionally stable for Stokes equations. For the full nonlinear Navier-Stokes equations the time stepping constraint is reduced to the standard CFL constraint . We also prove first order convergence of the gauge method when we use MAC grids as our spatial discretization. The optimal error estimate for the velocity field is also obtained.

     

A non-Newtonian fluid flow model for blood flow through a catheterized artery—Steady flow

In this paper the effects catheterization and non-Newtonian nature of blood in small arteries of diameter less than 100 μm, on velocity, flow resistance and wall shear stress are analyzed mathematically by modeling blood as a Herschel–Bulkley fluid with parameters n and θ and the artery and catheter by coaxial rigid circular cylinders. The influence of the catheter radius and the yield stress of the fluid on the yield plane locations, velocity distributions, flow rate, wall shear stress and frictional resistance are investigated assuming the flow to be steady. It is shown that the velocity decreases as the yield stress increases for given values of other parameters. The frictional resistance as well as the wall shear stress increases with increasing yield stress, whereas the frictional resistance increases and the wall shear stress decreases with increasing catheter radius ratio k (catheter radius to vessel radius). For the range of catheter radius ratio 0.3–0.6, in smaller arteries where blood is modeled by Herschel–Bulkley fluid with yield stress θ = 0.1, the resistance increases by a factor 3.98–21.12 for n = 0.95 and by a factor 4.35–25.09 for n = 1.05. When θ = 0.3, these factors are 7.47–124.6 when n = 0.95 and 8.97–247.76 when n = 1.05.     

Unsteady response of non-Newtonian blood flow through a stenosed artery in magnetic field

Current theoretical investigation of atherosclerotic arteries deals with mathematical models that represent non-Newtonian flow of blood through a stenosed artery in the presence of a transverse magnetic field. Here, the rheology of the flowing blood is characterised by a generalised Power law model. The distensibility of an arterial wall has been accounted for based on local fluid mechanics. A radial coordinate transformation is initiated to map cosine geometry of the stenosis into a rectangular grid. An appropriate finite difference scheme has been adopted to solve the unsteady non-Newtonian momentum equations in cylindrical coordinate system. Exploiting suitably prescribed conditions based on the assumption of an axial symmetry under laminar flow condition rendered the problem effectively to two dimensions. An extensive quantitative analysis has been performed based on numerical computations in order to estimate the effects of Hartmann number (M$M$), Power law index (n$n$), generalised Reynolds number (R_eG)$\left(R{e}_G\right)$, severity of the stenosis (δ)$\left(\delta \right)$ on various parameters such as flow velocity, flux and wall shear stress by means of their graphical representations so as to validate the applicability of the proposed mathematical model. The present results agree with some of the existing findings in the literature.     

Dynamic response of pulsatile flow of blood in a stenosed tapered artery

The aim of this paper is to throw some light on the rheological study of pulsatile blood flow in a stenosed tapered arterial segment. Arterial wall is considered to be rigid and flexible separately for improving the similarity to the in vivo situation. The streaming blood is considered to be Newtonian. The governing nonlinear equations of motion are sought using the well‐known stream function‐vorticity method and are solved numerically by finite difference technique. Important rheological parameters, such as axial velocity component, wall shear stress, and flow separation region are estimated in the neighborhood of the stenosis. Effects of stenosis height, vessel tapering, and wall flexibility on the blood flow are investigated properly and are explained in detail through their graphical representations.     

Mathematical modelling of unsteady flow of a Sisko fluid through an anisotropically tapered elastic arteries with time-variant overlapping stenosis

In the present investigation, we have studied the influence of heat and chemical reactions on blood flow through anisotropically tapered elastic artery with time-variant overlapping stenosis. The nature of blood in small arteries are analyzed mathematically by considering it as a Sisko fluid. The analysis is carried out for an artery with a mild stenosis. Analytical expressions for the axial velocity, the stream function, the temperature distribution, the concentration of fluid, the pressure gradient, the resistance impedance and the wall shear stress distribution have been computed numerically and the results were studied for various values of the physical parameters, such as the Sisko parameter b^∗, the power index n, the taper angle ?  , the maximum height of stenosis δ^∗${\delta }^{\ast }$, the Soret number S_r, the Brickmann number B_r, the total mass of the vessel and the surrounding tissues M and the longitudinal contributions of the viscous and elastic constraints to the total tethering C and K respectively. The obtained results show that the magnitude of the axial velocity is higher for a Newtonian fluid than that for a Sisko fluid and it decreases by increasing of the power index n also the transmission of axial velocity curves through a tethered tube is substantially higher than that through the free tube. The wall shear stress distribution and resistance impedance profiles with the time have an oscillation form through the tapered overlapping stenosed arteries and this oscillation decaying as the time increases. The temperature profile increase by increasing the Sisko parameter b^∗ and the power index n   but the concentration profile has an opposite behavior as compared to the temperature profile. For a fixed flux, the magnitude of the pressure drop for a shear-thinning fluid (n<1)$(n<1)$ is much larger than that through a shear-thickening (n>1)$(n>1)$. The stream lines separate and the trapping bolus appear by increasing the maximum height of the stenosis δ^∗${\delta }^{\ast }$. The trapping bolus increase in size toward the line center of the tube as the power index n increases and it appear gradually by increasing the Sisko parameter b^∗. Finally the size of trapped bolus for the stream lines in the free isotropic tube is smaller than those in the tethered tube.     

Local recovery of a solenoidal vector field by an extension of the thin-plate spline technique

We show how one may interpolate a vector-valued function in two or three dimensions, whose value is (wholly or partly) known at a sufficient (but not large) number of points disposed in almost any configuration, under the condition that the interpolating function has zero divergence. The technique is based on the theory of thin-plate splines. One may use a similar scheme in the case where the data consist of flux integrals (or other linear functionals) of the unknown function.     

Numerical simulation of a flexible fiber deformation in a viscous flow by the immersed boundary-lattice Boltzmann method

In this paper, deformation of a mass-less elastic fiber with a fixed end, immersed in a two-dimensional viscous channel flow, is simulated numerically. The lattice-Boltzmann method (LBM) is used to solve the Newtonian flow field and the immersed-boundary method (IBM) is employed to simulate the deformation of the flexible fiber interacting with the flow. The results of this unsteady simulation including fiber deformation, fluid velocity field, and variations of the fiber length are depicted in different time-steps through the simulation time. Similar trends are observed in plots representing length change of fibers with different values of stretching constant. Also, the numerical solution reaches a steady state equivalent to the fluid channel flow over a flat plate.     

旋风分离器内流动的三维分析

本文详细阐述了旋风分离器内流动在球坐标系中的数学表述和结果,应用质量守恒定律和定常流动的运动定律,在轴对称的考虑下,用流函数方法详尽推导了流动的三个速度分量.此讨论是从三维的整体观点来全面分析流动状况的.此外,对文[1]中的一些结果作了必要的修正.     

On two‐dimensional unsteady incompressible fluid flow in an infinite strip

In this paper, we study two‐dimensional incompressible fluid flow in an infinite strip. The stream function form of Navier–Stokes equation is considered, which keeps the physical boundary condition and avoids some difficulties in numerical simulations. The existence and uniqueness of global solution are proved. Some results on the regularity of solution are obtained. Copyright © 2000 John Wiley & Sons, Ltd.     

A uniformly stable family of mixed hp-finite elements with continuous pressures for incompressible flow

A new family of mixed hp-finite elements is presented for thediscretization of planar Stokes flow on meshes of curvilinear,quadrilateral elements. The elements involve continuous pressuresand are shown to be stable with an inf–sup constant boundedbelow independently of the mesh-size h and the spectral orderp. The spaces have balanced approximation properties—theorders of approximation in h and p are equal for both the velocityand the pressure. This is the first example of a uniformly stablemethod with continuous pressures for spectral element discretizationof Stokes equations, valid for geometrically refined meshesand curvilinear elements.     

Analysis of viscoelastic fluid flow with temperature dependent properties in plane Couette flow and thin annuli

The steady state flow in very thin annuli has been studied analytically for the case where the annular gap is much smaller than the radius of the inner cylinder and for the outer cylinder rotating at constant angular speed and the inner cylinder at rest. The cylinders were subjected to two different thermal boundary conditions. The exponential effect of temperature on the relaxation time and the viscosity coefficient was accounted into the governing differential equations using Nahme’s law. Effects of viscous dissipation as well as εDe² (viscoelastic index for SPTT constitutive equation) on the dimensionless velocity and temperature profiles have been investigated. Results show that while the properties of the fluid depend on temperature, the velocity and temperature profiles are different compared to those obtained with constant physical properties. The Nahme–Griffith number increases whereas εDe² as a viscoelastic index decreases when temperature dependent physical properties are considered. In addition, the results indicate that the viscous dissipation has a sensible effect on heat transfer and the Nusselt number decreases with an increase in the Nahme–Griffith number.     

Numerical simulation of laminar flow past a circular cylinder

The present paper focuses on the analysis of two- and three-dimensional flow past a circular cylinder in different laminar flow regimes. In this simulation, an implicit pressure-based finite volume method is used for time-accurate computation of incompressible flow using second order accurate convective flux discretisation schemes. The computation results are validated against measurement data for mean surface pressure, skin friction coefficients, the size and strength of the recirculating wake for the steady flow regime and also for the Strouhal frequency of vortex shedding and the mean and RMS amplitude of the fluctuating aerodynamic coefficients for the unsteady periodic flow regime. The complex three dimensional flow structure of the cylinder wake is also reasonably captured by the present prediction procedure.     

Global classical solution for a three-dimensional viscous liquid-gas two-fluid flow model with vacuum

This is a continuation of the paper （J. Math. Phys., 52（2011）, 093102）. We consider the Cauchy problem to the three-dimensional viscous liquid-gas two-fluid flow model. The global existence of classical solution is proved, where the initial vacuum is allowed.     

Series solution for unsteady axisymmetric flow and heat transfer over a radially stretching sheet

This paper deals with the unsteady axisymmetric flow and heat transfer of a viscous fluid over a radially stretching sheet. The heat is prescribed at the surface. The modelled non-linear partial differential equations are solved using an analytic approach namely the homotopy analysis method. Unlike perturbation technique, this approach gives accurate analytic approximation uniformly valid for all dimensionless time. The explicit expressions for velocity, temperature and skin friction coefficient are developed. The influence of time on the velocity, temperature and skin friction coefficient is discussed.     

A one dimensional model of blood flow through a curvilinear artery

We present a one-dimensional model describing the blood flow through a moderately curved and elastic blood vessel. We use an existing two dimensional model of the vessel wall along with Navier−Stokes equations to model the flow through the channel while taking factors, namely, surrounding muscle tissue and presence of external forces other than gravity into account. Our model is obtained via a dimension reduction procedure based on the assumption of thinness of the vessel relative to its length. Results of numerical simulations are presented to highlight the influence of different factors on the blood flow.     

Three-dimensional Couette flow of a dusty fluid with heat transfer

This work is focused on the mathematical modeling of three-dimensional Couette flow and heat transfer of a dusty fluid between two infinite horizontal parallel porous flat plates. The problem is formulated using a continuum two-phase model and the resulting equations are solved analytically. The lower plate is stationary while the upper plate is undergoing uniform motion in its plane. These plates are, respectively, subjected to transverse exponential injection and its corresponding removal by constant suction. Due to this type of injection velocity, the flow becomes three dimensional. The closed-form expressions for velocity and temperature fields of both the fluid and dust phases are obtained by solving the governing partial differential equations using the perturbation method. A selective set of graphical results is presented and discussed to show interesting features of the problem.     

Viscous potential flow analysis of interfacial stability with mass transfer through porous media

The interfacial stability with mass transfer, surface tension, and porous media between two rigid planes will be investigated in the view of viscous potential flow analysis. A general dispersion relation is obtained. For Kelvin-Helmholtz instability, it is found that the stability criterion is given by a critical value of the relative velocity. On the other hand, in the absence of gravity the problem reduces to Brinkman model of the stability of two fluid layers between two rigid planes. Vanishing of the critical value of the relative velocity gives rise to a new dispersion relation for Rayleigh-Taylor instability. Formulas for the growth rates and neutral stability curve are also given and applied to air-water flows. The effects of viscosity, porous media, surface tension, and heat transfer are also discussed in relation to whether the system is potentially stable or unstable. The Darcian term, permeability’s and porosity effects are also concluded for Kelvin-Helmholtz and Rayleigh-Taylor instabilities. The relation between porosity and dimensionless relative velocity is also investigated.