^{1}

^{*}

^{2}

^{1}

^{1}

In this research, we modeled MHD third grade blood flow in a stenosed artery. The blood viscosity and the density have been modeled into the shear thinning/thickening parameters, the most important rheological properties of blood. We used regular perturbation method and obtained the flow characteristics such as the flow velocity, the volume flow rate, the shear stress and the resistance to the flow considering a single layered stenosed artery. The results however showed that there is significant increase in volume flow rate and the velocity with increase in the magnetic field intensity H and the shear thinning Λ and reduces with increase in the shear thickening Ω .

The term viscosity (the ability of fluids to resist gradual deformation by shear or tensile stresses) is commonly used in fluid mechanics to characterize the shear properties of a fluid. It cannot be enough to describe non-Newtonian fluids like blood. Blood is a body fluid in humans and other animals, which delivers necessary substances such as nutrients and oxygen to the cells through some body canals called arteries and transports metabolic waste products away from those same cells through veins. So far, blood is best studied through several other rheological properties that relate stress and strain rate tensors under many different flow conditions such as oscillatory shear or extensional flow which is measured using different devices or rheometers. These properties are better studied using tensor-valued constitutive equations, which are common in the field of continuum mechanics. According to [

Arteries transport oxygenated blood to different tissues and organs of the body. This is essential for metabolic activities and sustaining life functions. Any narrowing of an artery will hamper the flow of blood and depending on the extent of the narrowing, known as stenosis, the artery may be partially or completely occluded (blocked) [

A non-Newtonian fluid is one whose flow curve (shear stress versus shear rate) is nonlinear or does not pass through the origin, i.e. where the apparent viscosity, shear stress divided by shear rate, is not constant at a given temperature and pressure but is dependent on flow conditions such as flow geometry, shear rate, etc. and sometimes even on the kinematic history of the fluid element under consideration.

The study of non-Newtonian fluid flow is important due to the nonlinear fluid rheology which is of special interest and has practical applications in industry and engineering. Several fluid models have been suggested to disclose the nature of non-Newtonian fluids. Among them is the differential type. The differential type fluid model has a simplest subclass known as the second-grade fluid, which describes the normal stress differences but cannot predict the shear thinning/thickening phenomena. However, the third-grade fluid model is capable of predicting both the normal stresses and the shear thinning/thickening phenomena, the most important rheological properties of blood. In [

The use of electromagnetic fields in the healing arts dates back as far as the fifteenth century, although the magnetic effects of “lodestone” were first described by the shepherd Magnes circa 1000 BCE in the region known today as Turkey [^{th} century, Mesmer began treating hysteria and other disorders (today recognized to be psychosomatic in origin) with lodestones [

In the 1930s, it was found that haemoglobin has magnetic properties that are different depending on whether it is carrying oxygen or not. When the haemoglobin is not carrying oxygen it is more sensitive or paramagnetic than oxygenated blood [

The effects of applying a magnetic field to liquid water have been intensively studied since 1980. it has been shown that the water vaporization rate, an essential process for all biological processes, is significantly affected by the application of a static magnetic field in both air and oxygen environments [

Although, the constitutive equations for third grade fluids have various complexities, many researchers have examined the flows of third-grade fluid under various aspects. We will take the study under a pulsatile condition and therefore begin by considering the flow of electrically conducting fluid (blood) in a stenosed artery.

The continuity and momentum equations are given as follows

d i v u ¯ = 0 (1)

ρ ¯ T ∂ ∂ t ¯ u ¯ T = − ∂ ∂ z ¯ P ¯ − 1 r ∂ ∂ r ¯ ( r ¯ τ ¯ T ) + μ ¯ 0 M ∂ ∂ z ¯ H ¯ (2)

Generally, for both Newtonian and non-Newtonian fluids, the Cauchy stress tensor is obtained by

τ = − P I + ∑ j = 1 n S j (3)

S j , j = 1 , 2 , 3 are called the stress tensors, P is the pressure force due to fluid flow. The pressure force is directed against the flow direction, which is the axial direction to induce the fluid flow. For the third grade fluid we have n = 3 and the first three tensors S j are given by

S 1 = μ ¯ A 1 (4)

S 2 = α 1 A 2 + α 2 A 1 2 (5)

S 3 = β 1 A 3 + β 2 ( A 2 A 1 + A 1 A 2 ) + β 3 ( t r A 1 2 ) A 1 (6)

where μ ¯ is the coefficient of sheer viscosity and α i , ( i = 1 , 2 ) , β i , ( i = 1 , 2 , 3 ) are material constants. A n are called Rivlin Ericksen tensors [

A n = D D t A n − 1 + A n − 1 ( ∇ u ¯ ) + ( ∇ u ¯ ) T A n − 1 , n > 1 (7)

A 1 = ( ∇ u ¯ ) + ( ∇ u ¯ ) T (8)

When β j = 0 , ( j = 1 , 2 , 3 ) , then, the above model reduces to second grade fluid model and if α i = 0 , ( i = 1 , 2 ) and β j = 0 , ( j = 1 , 2 , 3 ) , the model reduces to classical Navier stokes viscous fluid model [

Let the velocity field for the fluid flow be given in vector field. We assume the flow to be in the positive z-direction. This implies that the pressure gradient must be negative and the following must hold.

1) The velocity field is independent of the coordinate z and q. That is, u ¯ r r , u ¯ θ θ , u ¯ z z , u ¯ r θ , u ¯ r z , u ¯ θ r , u ¯ θ z , u ¯ z θ = 0.

2) The extra stress being the result of the velocity field, are also independent of the z and q. That is, τ r r , τ θ θ , τ z z , τ r θ , τ θ r = τ θ z = 0.

3) The sheer stress τ θ z and τ z θ on planes through the axis of the pipe are zero due to symmetry. This implies that the velocity field and the shear stresses are functions of r alone.

Thus, the constitutive equation of motion for a third grade fluid flow is

τ = μ ¯ ∂ u ¯ ∂ r ¯ + α 1 ∂ 2 u ¯ ∂ t ¯ ∂ r ¯ + 2 β 3 ( ∂ u ¯ ∂ r ¯ ) 3 (9)

Since the pressure gradient is a function of z ¯ and t ¯ we assume

− ∂ ∂ z ¯ p ¯ = q ¯ 0 + A 1 cos ω ¯ t ¯ , t ¯ ≥ 0 (10)

The boundary conditions are given as follows.

τ ¯ T is finite at r ¯ = 0 (11)

u ¯ T = 0 at r ¯ = R ¯ ( z ¯ ) (12)

The geometry of the stenosis is given in

R ¯ ( z ¯ ) = { R ¯ 0 − ( δ ¯ E 2 ) { 1 + cos [ 2 π l ¯ 0 ( z ¯ − d ¯ − l ¯ 0 2 ) ] } in d ¯ ≤ z ¯ ≤ d ¯ + l ¯ 0 R ¯ 0 inthenormalarteryregion (13)

where δ ¯ E is the height of the stenosis.

The key to solving modern problems is mathematical modeling. This process involves keeping certain elements, neglecting some, and approximating yet others.

We present the following non-dimensional variables

z = z ¯ R ¯ 0 , r = r ¯ R ¯ 0 , R ( z ) = R ¯ ( z ¯ ) R ¯ 0 , q ( z ) = q ¯ ( z ¯ ) R ¯ 0 , δ E = δ E ¯ R ¯ 0 , t = ω ¯ t ¯ , α T 2 = R 0 2 ¯ ω ¯ ρ ¯ T μ ¯ T , l = A ¯ 1 q ¯ 0 , u T = μ ¯ T q ¯ 0 R ¯ 0 2 / 4 μ ¯ T , R p = R ¯ p R ¯ 0 , τ T = τ ¯ T q ¯ 0 R ¯ 0 / 2 , τ N = τ ¯ N q ¯ 0 R ¯ 0 / 2 , H = H ¯ q ¯ 0 } (14)

Thus, Substituting (14) into (2), (9), (10) and (13) gives

α T 2 4 ∂ ∂ t u T = ( 1 + l cos ( t ) ) − 1 2 r ∂ ∂ r ( r τ T ) + F ∂ ∂ z H (15)

τ T = 1 2 ∂ u T ∂ r + α 2 Ω 1 ∂ 2 u T ∂ t ∂ r + α 2 Λ 1 ( ∂ u T ∂ r ) 3 ,

if 0 ≤ r ≤ R (16)

and the boundary conditions becomes

t is finite and ∂ u ∂ r = 0 (17)

u T = 0 at r = R (18)

R ( z ) = { 1 − ( δ E 2 ) { 1 + cos [ 2 π l 0 ( z − d − l 0 2 ) ] } in d ≤ z ≤ d + l 0 1 inthenormalarteryregion . (19)

Let us expand the flow characteristics in (15) and (16) in the perturbation series below taking α T 2 as the perturbation parameter.

τ T = τ T 0 + α T 2 τ T 1 + ⋯ u T = u T 0 + α T 2 u T 1 + ⋯

If we equate the coefficients of α 2 to have;

α T 0 : 2 r ∂ ∂ r [ r ( τ T 0 ) ] = 4 ( 1 + l cos ( t ) ) + F ∂ ∂ z H (20)

τ T 0 = 1 2 ∂ ∂ r ( u T 0 ) (21)

α T 2 : − 2 r ∂ ∂ r [ r ( τ T 1 ) ] = ∂ ∂ t ( u T 0 ) (22)

1 2 ∂ ∂ r ( u T 1 ) + Ω 1 ∂ 2 ∂ t ∂ r ( u T 0 ) + Λ 1 ( ∂ ∂ r ( u T 0 ) ) 3 = τ T 1 (23)

Thus, the solution set is given by

τ T = τ T 0 + α T 2 τ T 1 + ⋯ u T = u T 0 + α T 2 u T 1 + ⋯ } (24)

τ T = ( 1 + l cos ( t ) ) ( r ) + 1 4 F ∂ ∂ z H ( r ) + 1 8 α T 2 l sin ( t ) ( r 3 − 2 R 2 r ) (25)

u T = ( 1 + l cos ( t ) ) ( r 2 − R 2 ) + 1 16 Λ 1 α T 2 ( F ∂ ∂ z H ) 3 ( R 4 − r 4 ) + 3 4 Λ 1 α T 2 ( 1 + l cos ( t ) ) ( F ∂ ∂ z H ) 2 ( R 4 − r 4 ) + 2 Ω 1 α T 2 l sin ( t ) ( r 2 − R 2 ) + 1 4 F ∂ ∂ z H ( r 2 − R 2 ) + 4 Λ 1 α T 2 ( 1 + l cos ( t ) ) 3 ( R 4 − r 4 ) + 3 Λ 1 α T 2 ( 1 + l cos ( t ) ) 2 ( F ∂ ∂ z H ) ( R 4 − r 4 ) + 1 16 α T 2 l sin ( t ) ( 4 R 2 r 2 − 3 R 4 − r 4 ) (26)

The volume flow rate is given by

Q = 2 π ∫ 0 R ( z ) u T ( r , t ) r d r

Q = 2 π ( − 1 4 ( 1 + l cos ( t ) ) ( R 4 ) − 1 16 F ∂ ∂ z H ( R 4 ) + 1 32 Λ 1 α T 2 ( F ∂ ∂ z H ) 3 ( R 6 ) − 1 2 Ω 1 α T 2 l sin ( t ) ( R 4 ) + 1 4 Λ 1 α T 2 ( 1 + l cos ( t ) ) ( F ∂ ∂ z H ) 2 ( R 6 ) + 4 3 Λ 1 α T 2 ( 1 + l cos ( t ) ) 3 ( R 6 ) − 1 24 α T 2 l sin ( t ) ( R 6 ) + 2 3 Λ 1 α T 2 ( 1 + l cos ( t ) ) 2 ( F ∂ ∂ z H ) ( R 6 ) ) (27)

The resistance to the flow is given by

Λ = − ∂ ∂ z ( P ) Q

Λ = ( 1 + l cos ( t ) ) 2 π ( − 1 4 ( 1 + l cos ( t ) ) ( R 4 ) − 1 16 F ∂ ∂ z H ( R 4 ) + 1 32 Λ 1 α T 2 ( F ∂ ∂ z H ) 3 ( R 6 ) − 1 2 Ω 1 α T 2 l sin ( t ) ( R 4 ) + 1 4 Λ 1 α T 2 ( 1 + l cos ( t ) ) ( F ∂ ∂ z H ) 2 ( R 6 ) + 4 3 Λ 1 α T 2 ( 1 + l cos ( t ) ) 3 ( R 6 ) − 1 24 α T 2 l sin ( t ) ( R 6 ) + 2 3 Λ 1 α T 2 ( 1 + l cos ( t ) ) 2 ( F ∂ ∂ z H ) ( R 6 ) ) (28)

The wall shear stress is given by

τ w = τ T | r = R (29)

τ w = ( 1 + l cos ( t ) ) ( R ) + 1 4 F ∂ ∂ z H ( R ) − 1 8 α T 2 l sin ( t ) ( R 4 ) (30)

We present this paper with the view to analyze the effect of magnetic field on pulsatile third grade blood flow in a stenosed artery. We have used the following parameter values; 0.2 is used for the dimensionless amplitude l. The Womersley numbers α T is taken as 1.2. The Womersley number denotes the ratio of unsteady inertial forces to viscous forces in the flow. It ranges from as large as about 20 in the aorta, significantly greater than 1 in all large arteries, to as small as 10^{−3} in the capillaries [

and thickening respectively. Figures 8-10 show the variation of the volume flow rate with time at different values of the shear thinning, shear thickening and magnetic field intensity H. The graphs shift away from the origin when the shear thinning and the magnetic field intensity increases and shift closer to the origin as the shear thickening increase. This shows that the volume flow rate increases with increase in shear thinning and magnetic field but reduces over time, the volume flow rate hits the zero bar after some time, implying a possible case of occlusion. This result is in agreement with the laboratory results of the work of [

Figures 11-13 show variation of velocity profile of the third grade fluid along the radial distance for different values of shear thickening, shear thinning and magnetic field H,with d = 2 , l = 0.2 , l = 0.2 , l = 1 , δ T = 0.15 . The velocity curve shift away from the origin when the shear thinning & magnetic field increase (

Figures 14-16 show the variation of the total resistance to the volumetric flow rate in the artery geometry along the axial direction for different values of shear thickening, shear thinning and magnetic field intensity H, with d = 2 , l = 0.2 , l = 0.2 , l = 1 , δ T = 0.15 . We observed that increasing shear thickening offers high resistance to the volume flow rate (

source to drive the flow of blood thereby weakening the muscular walls of the heart resulting to myocardial infarction (in which part of the heart muscle dies) or cardiac arrest (in which blood flow stops altogether). In [

opines that this is the cause of the sudden break down of the entire circulatory system leading to death in recent time. Increasing shear thinning reduces the resistance to flow (

field intensity increases (

observed that the wall shear stress behave similar to the volume flow rate. The graphs shift away from the origin as the magnetic field increases with a periodic

relation in

Any fluid that does not obey Newton’s Law of Viscosity is a non-Newtonian fluid. The viscosity of such fluids is mostly dependent on their shear rate or the shear rate history. In fluid mechanics, the fluid property called viscosity is used

to characterize the shear properties of a fluid, which is inadequate to describe non-Newtonian fluids. Non-Newtonian fluids are best studied through several other rheological properties that relate stress and strain rate tensors under many different flow conditions. An example of a non-Newtonian fluid with the shear thinning and thickening properties is blood which we considered as third grade fluid. This application is highly valued for its importance in understanding flows in defective arteries within the body, as it allows us to relate the magnetic field and blood to further understand the flow characteristics of blood through a stenosed artery.

Ikpakyegh, L.N., Okedayo, G.T., Aboiyar, T. and Onah, E.S. (2018) Analysis of Pulsatile Magnetohydrodynamic (MHD) Third Grade Blood Flow in a Stenosed Artery. American Journal of Computational Mathematics, 8, 78-95. https://doi.org/10.4236/ajcm.2018.81007

R ¯ ( z ¯ ) - Radius of the stenosed artery

δ ¯ E - Maximum height of the stenosis

l 0 ¯ - The length of stenosis.

d ¯ - Location of the stenosis.

z ¯ - The axial distance.

r ¯ - The radial distance.

τ ¯ T - The shear stress for the Third grade fluid.

ρ ¯ T - The density of the Third grade fluid.

p ¯ - The pressure force.

μ ¯ T - Viscosity of the Third grade fluid.

U ¯ T - Velosity of the Third grade fluid.

M - Magnetisation of the fluid.

μ 0 - The magnetic permeability.

H ¯ - The magnetic field intensity.

R ¯ p - The plug core radius.