The Influence of Drag on Nonlinear Oscillatory Flow through Concentric Annulus

B.Umadevi,C.V.Vinay and P.A.Dinesh

1Department of Mathematics, JSS Academy of Technical Education, Bangalore, India

2Department of Mathematics, Ramaiah Institute of Technology, Bangalore,India

ABSTRACT A mathematical model has been developed to study the effect of particle drag parameter and frequency parameter on velocity and pressure gradient in nonlinear oscillatory two phase flow.The main purpose is to apply the model to study the combined effect of introduction of the catheter and elastic properties of the arterial wall on the pulsatile nature of the blood flow.We model the artery as an isotropic thin walled elastic tube and the catheter as a coaxial flexible tube.Blood is modeled as an incompressible particulate viscous Newtonian fluid.Perturbation technique has been applied to find the approximations for velocity and pressure gradient up to second order.Numerical solutions are investigated with graphical presentations to understand the effects of drag parameter,frequency parameter and phase angle on velocity along radial direction and pressure gradient along axial directions.As the drag parameter increases,mean pressure gradient and mean velocity will be decreased.As frequency parameter increases mean velocity profile bends near the outer wall.Due to elastic nature of artery wall, a thin catheter experience small oscillations and a thick catheter remains stationary inside the artery.Finally,the effect of catheterization on various physiologically important flow rate characteristics—mean velocity,mean pressure gradient are studied for a range of different catheter sizes, particle drag parameter and frequency parameters.

KEYWORDS Particle drag; oscillatory flow;catheter

Nomenclature

(R,θ,Z):Cylindrical coordinate system

(U,0,W):Velocity vector of the fluid phase

(U′,0,W′):Velocity vector of particulate phase

P:Pressure

ρ:Density of the fluid

t:Time

N:Number density of the particles

ν:Kinematic viscosity of the fluid

K=3 μ πd:stokes drag term

m:Mass of the each drag particulate

(u,0,w):Non dimensional velocity of fluid phase

(u′,0,w′):Non dimensional velocity of particulate phase

d:Diameter of the particulate

(r,θ,z):Non-dimensional coordinates

P0:mean pressure

C0:characteristic wave speed

Ω:Angular velocity

D:Distensibility of the flexible tube

b:Non-dimensional boundary radius

R~max:Maximum radius

α:Womersley number

τ0:Phase lead angle

Rst:Steady streaming Reynolds number

ω:Drag parameter

k:Catheter Radius

w11:Mean axial velocity of first order

w20:Mean axial velocity of second order

p′20(z):Mean Pressure gradient of second order

φ:Phase difference angle

wc:Amplitude of catheter oscillation

s:Radial directional length

Z:Axial directional length

q(z,t):Non dimensional flow rate

1 Introduction

Oscillatory flow is a widespread phenomenon and plays an important role in many fields,e.g.pneumatic propulsion, piston—driven flow, and acoustic oscillation are commonly used in mechanical engineering;pulsatile blood circulation, respiratory flow in lung, and capillary waves are of much interest in biomechanics; seasonal reversing wind, ocean circulations as well as tide flow are of high concern in meteorology, etc.More than mere oscillation or repetition, mass, momentum, and energy may be transferred via these reciprocating movements.Oscillatory flow of fluid enclosed in distensible vessels is of substantial significance due to its applications to various fields.The model can be used to understand(i) Blood flow in large arteries, (ii) wave propagation in fluid filled flexible tubes, which is essential for the study of acoustics or pulse propagation in arteries, (iii) flow in collapsible veins and (iv) peristaltic motion in the intestines.Many researchers for over 180 years are investigating the blood flow in arteries.Comparison of experimental measurements and mathematical modeling of these problems have been investigated.Recently, the study of flow with periodic variations has attracted much attention of researchers due to its various engineering and physiological applications.Oscillatory motion of a viscous liquid in a thin-walled elastic tube is investigated by Womersley [1].Further, Womersley [2] studied the elastic tube theory of pulse conduction and oscillatory flow in mammalian arteries.Saffman [3] studied the stability of laminar flow of dusty gas.Rubinow et al.[4] analyzed the flow of a viscous fluid in an elastic tube with application to blood flow.RamachandraRao et al.[5] studied the pulsatile flow in tubes of various cross section.Taylor et al.[6] presented a mathematical model to analyze the blood flow through arteries and expressed the different pressure radius relationships for elastic tube.RamachandraRao[7] investigated the oscillatory flow in an elastic tube of variable cross section.Analytical solution by the method of linear approximation to describe the velocity distribution for laminar periodic flow through porous walls is proposed by Chang et al.[8].Wang et al.[9] studied non-linear analysis of oscillatory folw with a non zero mean,in an elastic tube(artery).Srivastava[10]investigated the two phase model of blood folw through stenosed tubes in the presence of peripheral layer.Sarkar et al.[11]investigated the Correction to Flow Rate-Pressure Drop in Coronary Angioplasty: Steady Streaming Effect.Sarkar et al.[12] discussed the nonlinear analysis of oscillatory folw in the annulus of an elastic tube.In addition, they noticed that the velocity distribution in a small vessel depend significantly on geometry of the wall and its elastic nature.Jayaraman et al.[13] studied the Nonlinear analysis of arterial blood folw—steady streaming effect.Vajravelu et al.[14] considered the case of inserting a catheter into an elastic tube to observe the changes in blood folw pattern by taking Herschel–Bulkley fulid.Unsteady folw of a Jeffrey fluid in an elastic tube with a stenosis was considered by Sreedharamalle et al.[15].Sochi [16] projected the expression for the volumetric flow as a function of pressure in elastic tube using two pressure area constitutive relationships.Sankar et al.[17] analyzed the cooling of heat sources by natural convection heat transfer in a vertical annulus, Siddiquiet al.[18] investigated mathematical analysis on pulsatile flow through a catheterized stenosedartery.Khudayarov et al.[19] investigated the mathematical simulation of nonlinear oscillations of viscoelastic pipelines conveying fulid.Do et al.[20] studied the Navier’s slip condition on time dependent Darcy—Forchheimer nano fulid using Spectral relaxation method.Tsimpoukis et al.[21] investigated the nonlinear oscillatory fully-developed rarefied gas flow in plane geometry.Kiran et al.[22] studied the Computational analysis of conjugate buoyant convective transport in an annulus.Naveed et al.[23]investigated the mathematical analysis of novel coronavirus(2019-nCov)delay pandemic model.Azam et al.[24] studied the numerical modeling and theoretical analysis of a nonlinear advection-reaction epidemic system.Saqib et al.[25]investigated the Symmetric MHD channel folw of nonlocal fractional model of BTF containing hybrid nanoparticles.Bilal et al.[26] investigated Finite element method about heat transfer analysis of Newtonian material in triangular cavity with square cylinder.Rasool et al.[27] studied the MHD squeezed Darcy-Forchheimer nano fulid flow between two h-distance apart horizontal plates.Lu et al.[28]studied the oscillatory two-phase folw in microchannels.Bilal et al.[29] studied the Analytical treatment of radiative Casson fluid over an isothermal inclined Riga surface with aspects of chemically reactive species.

Blood is a concentrated suspension of several formed cellular elements, red blood cells (RBCs or erythrocytes), white blood cells (WBCs or leukocytes) and platelets (thrombocytes), in an aqueous polymeric and ionic solution, the plasma, composed of 93% water and 3% particles, namely electrolytes,organic molecules, numerous proteins (albumin, globulins and fbirinogen) and waste products.Plasma’s central physiological function is to transport these dissolved substances, nutrients, wastes and the formed cellular elements throughout the circulatory system.The primary function of erythrocytes is to transport oxygen and carbon dioxide.Leukocytes are roughly spherical and much larger than erythrocytes, but they exist in a smaller number in blood.Leukocytes are subdivided into granulocytes (65%), lymphocytes(30%), monocytes (5%) and natural killer cells.Granulocytes are further subdivided into neutrophils (95%),eosinophils(4%)and basophils(1%).The leukocytes play a vital role in fighting infection and thus are able to migrate out of the blood vessels and into the tissues.Thrombocytes are small discoid non-nucleated cell fragments, much smaller than erythrocytes and leukocytes.Thrombocytes are a vital component of the blood clotting mechanism.The total volume concentration of leukocytes and thrombocytes is only about 1%.Blood cells are continually produced by the bone marrow over a human’s life.

In the present investigation Blood is modeled as an incompressible particulate viscous Newtonian fluid,and the flow will be fully developed.A mathematical model has been developed to view the oscillatory nonlinear flow in the annulus formed by catheterized artery.The analysis of the problem has been given by perturbation method.The variation of velocity and pressure observed with the influence of the parameters such as elastic parameter,Womersley number and drag parameter.

2 Mathematical Formulation

The artery modeled as an isotropic thin walled elastic tube with mean radiusR0and the catheter as a coaxial flexible tube with radiuskR0(k＜1).Blood is considered as an incompressible particulate viscous Newtonian fluid.The flow is oscillatory and axial-symmetric.The oscillatory nature of the flow will have an influence on the instant position of the flexible catheter.The catheter movement will influence the flow there it is assumed that to have small constant amplitude and a small phase lead over the rate flow.This is necessary because,in general,in the case of unsteady boundary layer flows,the stress has a phase lead over the rate of velocity such that the amplitude is taken as very small.The oscillatory phenomena of the wall will generate longitudinal movement in the tube wall, on the other hand vascular tethering of the arteries has a damping effect on this longitudinal movement and it is assumed to be negligible, leading to simplification of the equations describing the wall motion.

Fig.1 shows a schematic diagram of the annular geometry and the cylindrical coordinate system(R,θ,Z) withZalong the axis of the tube.(U,0,W) be the velocity vector of the fluid phase and(U′,0,W′) be the velocity vector of particulate phase corresponding to the cylindrical coordinate system(R,θ,Z).The Navier–Stokes equations corresponding to this phenomenon from Saffman [3], Srivastava[10] and Sarkar et al.[12]are given by

where P is the pressure,ρ is the density the fluid,tis the time,Nis the number density of the particles,ν is the kinematic viscosity of the fluid,K= 3 μ πdis the stokes drag term,mis the mass of the each particulate anddis diameter of the particulate.

Figure 1: Physical configuration

whereWc(t)represents the pulsatile movement of the flexible catheter.

Non-dimensional transformation variables given as follows:where(u,0,w)be the non dimensional velocity of fluid phase,(u′,0,w′)be the non dimensional velocity of particulate phase,(r,θ,z) be the non dimensional coordinates,P0is the mean pressure,C0denotes a characteristic wave speed, Ω is the angular velocity, D is the distensibility of the flexible tube,bis the non-dimensional boundary radius is function of pressurep,qis the flow rate,is the maximum radius attained by the tube wall due to the influence of the oscillatory pressure on the elastic wall of the tube.We thus describe the parameters governing to the flow as Womersley number, α =and steady streaming Reynolds number asRst= ε2α2.IfRstis small, the steady streaming is the same to Stokes flow and when it is large the steady streaming may also have a tractable form.

In the Eq.(8)the axial viscous transport term is ignored because

The boundary conditions are

Introducing the non dimensional transformation to overcome the difficulty of moving boundary as

Thus the equations of motion and continuity in the transformed coordinate are

whereb1(s,z, τ) =(b(p) -k).

The consequent boundary conditions are given by

The relation between the cross sectional area of the tube and the transmutable pressure difference—has to be specified to understand the flow in inert tubes.That is in the present problem,specifyingb(p)will be given in tube law.

3 Method of Solution

Being nonlinear in nature,the Eq.(14)does not provide the closed form solution.By using perturbation method, the approximate solutions for velocity and pressure distribution up to second order have been analyzed.The diameter variation in arteries, resulting from the effect of pressure pulse, represented by ε is considered very small.This enables us to seek a solution by expanding for fluid velocity, particulate velocity in axial direction and radial direction, pressure and non-dimensional boundary radius i.e.,w,w′,p,b,u,u′in terms of ε ≤1 given below:

wi jrepresents thejth harmonic of theith order term in the perturbed expansion ofw.Similar definitions hold forw′i j,pi j,bi j,ui j&u′i j.ε ≤1 allows the frequency parameter α,to take moderate values withRst≈O(1 ).Again,considering purely oscillatory flow rate we can consider

where Re [ ] denotes the real part of a complex variable.Associated expressions can be written forw′11,u11,u′11,p11,b11.We take=wccos (τ - τ0) wherewc≤ 1 is the maximum amplitude of the flexible annular wall.τ0is the phase lead of this wall oscillation over the flow.In general thewcwill depend on the elastic properties of the catheter material and the tangential stress exerted by the fluid flow.However, in the present study we treatwc& τ0given.Implementing the Eqs.(19) and (20) in Eqs.(14)–(18) and equating the terms of varying orders of ε, we get a system of differential equations for various order terms.

3.1 Solution of O(ε)

The equations corresponding to first orderO(ε ) terms which are proportional toeiτare the following:

The corresponding boundary conditions are

After solving forU11&W11in terms of modified Bessel functions of complex arguments,solutions ofO(ε )for the Eqs.(14) and(19) are obtained as

The boundary conditions on radial velocity,i.e., Eq.(12),is

where

The ( `) denotes derivative with respect toz.F1(s),F2(s),F4(s),A4,C12are defined in the Appendix.f(z) be the unknown function which depends on the elastic nature of the tube and its response to the pressure gradient.

3.2 Steady Streaming Solution of O(ε2)for Arbitrary α and Small Rst

TheO(ε2)equations for the steady streaming components are given by

where

A double over bar denotes time averaged solution.

Boundary conditions:

The solution for pressure and velocity are obtained as follows:

where

H1(s),H2(s),defined in the Appendix.D1andD2are pure constants which are evaluated using the Eqs.(31)and(33).It is to be distinguished that the first term in Eq.(34)component due to the movement of the flexible inner tube.It is to be identified thatp′20is the correction to the mean pressure—an important result from the nonlinear steady streaming study.The effect ofp′20will depend on α,ε,wc, τ0,f′(z) and the functionF(z) =In addition, it was revealed thatF(z) is proportional to the wall movement and flow rate amplitude atO(ε ).

In the current study, the longitudinal oscillations of the catheter give rise to an additional term whose magnitude will depend uponwc, τ0,f′(z).Thus the interface of the amplitude of catheter oscillation(wc) and the amplitude of the wall movement which is proportional tof′(z) is first felt at theO(ε2)throughp′20andw20.

3.3 Dependence on Axial Position

f(z)is an unknown function and it is related to both pressure gradient and diameter variation, is assumed as

where the constants of integration ofA1,A2complex in nature and can be determined from the dynamic conditions of the flow.

The flow rate of wave form is determined by the pressure wave,movement of the wall,and in a small measure by the movement of the catheter.The non dimensional form of the flow rate can be written as

where the first term is due to the catheter movement and the second is due to the oscillatory pressure gradient.Q1will depend onwc, contributes to the modified amplitude,Q2is the amplitude in the absence of the oscillations of the flexible inner tube and φ(z) is the phase angle difference between the wall motion and flow rate in the absence of oscillation of the flexible inner tube atO(ε ).

A1,A2are estimated by setting up the flow rate amplitude when the inner tube is at rest and amplitude of wall motion.Q1has been fixed from the given values forwc.From experiments with and without the presence of the movement of the annulus and know the flow rate and diameter deviation at a particular point in the axial direction.

4 Results and Discussion

The present investigation reveals the influence of the particle drag parameter on nonlinear oscillatory flow through concentric annulus.The variation of velocity profile and pressure gradient have been analyzed with respect to the effects of drag parameter, frequency parameter and phase angle.In this section we discuss our results corresponding to values for the parameters in the model suitable to the physiological system.The frequency parameter α is taken from a small value of 5 to a moderate value up to 15,kis the ratio of the radii of the inner tube is varied from 0.2 to 0.5 andwcis given values from 0 to 0.4 and drag parameter ω is taken in the range 0.5–0.8.HereQ2(z)fixed as 0.5 and the amplitude of wall variationB1(z) as 0.05 atz= 0 for the purpose of analyzing the results of our mathematical model.The phase difference angle φ be the notion as an indicator of the impedance to the flow.Values for φ ranging from 0° to 90° are given to study the influence of wave reflections on the steady streaming induced mean pressure gradient and velocity profile.

4.1 Mean Pressure Gradient and Velocity Distribution

The induced mean pressure gradientp′20(z) along the axial direction for different values of ω drag parameter have been presented in the Figs.2–5.From these it reveals that as drag parameter increases mean pressure gradient decreases.It is seen from Fig.6 that the phase lead τ0of the wall of the oscillating catheter to flow, it does not have considerable influence on pressure gradient.The effect of catheter radius on mean pressure gradientp′20(z) can beseen in the Figs.7 and 8.From these, it reveals that as catheter radiuskincreases mean pressure gradient decreases.

Variations of mean axial velocityw11andw20along the radial direction for various values of drag parameter have been presented in the Figs.9–12, for various values of drag parameter and frequency parameter with fxied values of phase angle, and catheter radius.It discloses that the parabolic profile has been maintained that similar to annular flow in a rigid tube.It can be observed that as drag parameter increases then the velocity decreases, that is particulate drag resists the flow of the fluid.But as frequency parameter increases the parabolic profile shows kinks near the outer wall (artery wall) of the tube suggesting the influence of the elastic nature of the tube wall on the flow.The obtained results have been agreed with that of Wang et al.[9]and also Sarkar et al.[12].

Figure 2: Variation of mean pressure gradient along the axial length with drag parameter 0.5

Figure 3: Variation of mean pressure gradient along the axial length with drag parameter 0.6

Figure 4: Variation of mean pressure gradient along the axial length with drag parameter 0.7

Figure 5: Variation of mean pressure gradient along the axial length with drag parameter 0.8

Figure 6: Variation of mean pressure gradient along the axial length for different phase lead angles

Figure 7: Variation of mean pressure gradient along the axial length with Catheter radius 0.3

Figure 8: Variation of mean pressure gradient along the axial length with Catheter radius 0.4

Figure 9: Variation of mean axial velocity of first order along the radial direction for various values of drag parameter with frequency parameter value 5

Figure 10: Variation of mean axial velocity of first order along the radial-direction for various vales of drag parameter with frequency parameter value 10

Figure 11: Variation of mean axial velocity of second order along the radial-direction for various vales of drag parameter with frequency parameter value 5

Figure 12: Variation of mean axial velocity of second order along the radial-direction for various vales of drag parameter with frequency parameter value 10

5 Conclusion

To study the nonlinear oscillatory flow in the concentric annulus with the influence of particle drag and frequency parameter,a mathematical model has been developed.The model has been applied to investigate induced mean pressure gradient, mean velocity distribution—vary noticeably along axial, radial correspondingly.The computational results were found in good agreement with the bench mark results.As particle drag parameter increases then mean pressure gradient and mean velocity distribution decreases.Mean velocity proflie bends near the outer wall (artery wall) of the tube, as frequency parameter increases, suggesting the influence of the elastic nature of the tube wall on the flow.Also depending on catheter size, a slim catheter experience the minute oscillations due to the flow conditions is likely to influence in the same way as thicker catheter which remains quite stationary inside the artery.

This study gives an insight to enormous applications in the field of diagnosing, treating and certain surgical procedures related to the disorders/diseases which originate in the body relating to COVID-19 symptoms, cardiovascular, pulmonary, synovial systems etc., where as the different types of cardiovascular diseases include Aneurysms, Angina, Atherosclerosis, Stroke, different types of cerebrovascular disease,Heart Failure,Coronary Heart diseases and Myocardial infarction or Heart attacks.

Acknowledgement: The authors acknowledge the support of JSS Academy of Technical Education,Bangalore and Ramaiah Institute of Technology,Bangalore.

Funding Statement:The authors received no specific funding for this study.

Confilcts of Interest:The authors declare that they have no conflicts of interest to report regarding the present study.

Appendix

where