Next Article in Journal
M-YOLO: Traffic Sign Detection Algorithm Applicable to Complex Scenarios
Next Article in Special Issue
Magneto-Nanofluid Flow via Mixed Convection Inside E-Shaped Square Chamber
Previous Article in Journal
Fixed-Time Synchronization Analysis of Genetic Regulatory Network Model with Time-Delay
Previous Article in Special Issue
Significance of Chemical Reaction and Lorentz Force on Third-Grade Fluid Flow and Heat Transfer with Darcy–Forchheimer Law over an Inclined Exponentially Stretching Sheet Embedded in a Porous Medium
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analytical Solution for the MHD Flow of Non-Newtonian Fluids between Two Coaxial Cylinders

1
Shanghai Key Lab of Vehicle Aerodynamics and Vehicle Thermal Management Systems, Shanghai 201804, China
2
School of Automotive Studies, Tongji University, Shanghai 201804, China
3
Department of Mathematics, University of Baltistan Skardu, Skardu 16200, Pakistan
4
Directorate of Education Babylon, Ministry of Education, Babylon 51014, Iraq
5
School of Mechanical Engineering, Tongji University, Shanghai 201804, China
*
Author to whom correspondence should be addressed.
Symmetry 2022, 14(5), 953; https://doi.org/10.3390/sym14050953
Submission received: 17 March 2022 / Revised: 31 March 2022 / Accepted: 1 April 2022 / Published: 7 May 2022
(This article belongs to the Special Issue Symmetry in CFD: Convection, Diffusion and Dynamics)

Abstract

:
This paper deals with the MHD peristaltic flow of Williamson fluids through a porous medium between two joint cylinders. The fluid flow was considered to be that of a non-Newtonian fluid, i.e., a Williamson fluid. The inner tube was uniform, while the flexible outer tube had a Sine wave moving down its wall. The analytical solutions for velocity and temperature were obtained as functions (Bessell functions of the first and second types). The solution for velocity profile, temperature, and concentration distribution were obtained as functions of the physical parameters of the problem (Darcy number, magnetic parameter, Grasoff thermal number, Reynolds number, Prantl number, and Schmidt number) along with other physical parameters. The effect of the physical parameters was discussed graphically. A comparison with previously published graphical results was also carried out. The ambition of the present paper is to contribute to practical applications in geographical and physiological fluid dynamics, such as on sandstone, in the human lungs, on beach sand, on limestone, and in the bile duct. This study is based on theoretical research and can be helpful in the fields of fluid mechanics and mathematics.

1. Introduction

Research on non-Newtonian fluids is important for the development of fluids, lubricants, and plastic products. In fact, there are numerous fluids whose flow is considered non Newtonian, such as plasma, liquid metals, nuclear fuel slurries, Bingham fluids, blood, etc. The flow of these fluids does not follow the Newtonian law of viscosity; therefore, it is more appropriate to consider them non-Newtonian fluids. Various researchers have analyzed different types of non-Newtonian fluids. These include Ellis fluids, nanofluids, Jaffery fluids, and Williamson fluids [1,2,3,4,5]. Williamson fluid is non-Newtonian fluid with shear thinning property. Several investigations have been made on the Williamson fluid model. For instance, the peristaltic flow of a Williamson fluid in a curved channel under the influence of a magnetic field was studied by Rashid et al. [6]. They solved the coupled partial differential equation by using the analytical technique known as Homotopy Perturbation Method (HPM) under the assumption of long wavelength and low Reynolds number. After obtaining the solution, they graphically analyzed the effects of pertinent parameters on velocity, temperature profile, etc. In another study by Nadeem and Akram [7], the peristaltic flow of a Williamson fluid in an asymmetric channel was discussed and it was concluded that the pressure rise is not linear for large values of the Williamson parameter, whereas, for small values of this parameter, the pressure rise behaves like that of a Newtonian fluid. Abdulmajeed [8] discussed the Williamson fluid considering heat and mass transfer. I Their examination concluded that the Williamson parameter can play an important role in momentum transport. The second law of MHD slip flow for Williamson fluids was analyzed by Dadheech et al. [9]. The analysis discovered that rate of entropy generation increased in relation to the magnetic field parameter, while chemical reactions and opposite effects on the porosity parameter were observed. Asjad MI et al. [10] developed a mathematical model of Williamson fluid in the presence of bioconvection. The mathematical model was solved by using the Range–Kutta method with shooting, and it was found that the velocity increased with mix convection. Shashikumar NS et al. [11] explored the steady flow of a Williamson fluid in a micro-channel, in the presence of viscous dissipation, magnetic effects, and Joule dissipation.
In recent years, due to the countless importance of physiological and engineering applications, peristaltic transportation has attracted researchers’ interest. The term “peristaltic “ derives from the Greek word paristaltikos and indicates a sort of wave motion that appears within tubular structures and compels the movement of an object or medium. Some prominent examples of peristaltic motion are the movement of chyme in the gastrointestinal region, of fluid in the male reproductive system, the movement imparted by roller and finger pumps, and that of blood in small vessels. Researchers carried out several examinations on these types of motions. Particularly, Abbas et al. [12] discussed the peristaltic flow of blood using a nanofluid in the presance of magnetohydrodynamics (MHD) by using the perturbation method and numerical integration and determined pressure and friction force. They found that the magnetic parameter and the Brownian motion parameter have an important effect on the velocity profile and pressure distribution. Such findings are useful to control bleeding during surgery. In another study, the same authors [13] examined the entropy generation of the peristaltic flow of blood within compliant walls and found that the elastic parameters of the walls affected the temperature profile and can be modified to minimize entropy generation as blood flows. Abbas et al. [14] investigated the minimization of entropy generation of the peristaltic flow of a nanofluid with compliant walls. Sinha et al. [15] discussed the MHD of peristaltic flow and heat transfer in an asymmetrical channel. The peristaltic flow of a hybrid nanofluid along a curve was analyzed by Saleem Akbar et al. [16]. A. Riaz et al. [17] developed a mathematical model of peristaltic flow with an asymmetrical wavy motion with the influence of entropy generation. Abbas et al. [18] used an Artificial Neural Network to predict and minimize entropy generation by the peristaltic flow of a non-Newtonian nanofluid.
The flow of fluids through porous media is a significant area of research. It has great applications in geophysical fluid dynamics. Flow on sandstone, in the human lungs, in the gall bladder in the presence of stones, and that of blood in small vessels are well-known examples of porous media flow. The first studies on peristaltic porous media made by Afifi [19] and Bhatti et al. [20] analyzed the peristaltic flow of blood through a porous medium under the effect of slip and MHD. They discussed different cases in which the flow became that of a Newtonian fluid and found that the magnetic fields can be used for the magnetic targeting of drugs in cancer treatment.
Noreen et al. [21] explored heat transfer on electroosmotic flow through a porous medium. Pattnaik P. K. et al. [22] experimentally examined metallic nanoparticles (Cu, A l 2 O 3 , and SWCNT) with effects of MHD. More recently, other interesting research works related to porous media were carried. Harelip et al. [23] developed a continuous mathematical model with the time probability density function of a poly-dispersed fluid in a porous medium. M. Yang et al. [24] considered MHD for an Eyring-Powell fluid in a porous medium. They analyzed the effects of energy and viscosity on the flow. Nawaz M and M. Adil Sadiq [25] created a mathematical model for non-Newtonian hybrid nanoparticles and solved it by using the finite element method. Other more recent studies were conducted, as reported in the Reference section [26,27,28,29].
In all the above-mentioned studies, no attention was paid to the MHD peristaltic flow of Williamson fluids through a porous medium in two coaxial cylinders. Shaaban and Abou-zeid [30] considered the geometry of two coaxial cylinders for the MHD peristaltic flow of a non-Newtonian fluid in a porous medium with the effects of heat and mass transfer. Similarly, the MHD peristaltic flow of a micro polar bi-viscosity fluid thorough a porous medium and two coaxial tubes was investigated by Nabil and Mohamed [31]. Among the non-Newtonian fluids, the most used fluids are pseudoplastic fluids. The Williamson fluid is also a type of pseudo plastic fluid and includes polish, paints, whipped cream, blood etc. In the current study, the Williamson model was used in solving the governing equations. The coupled partial differential equation was solved with the Bessel function subject to the boundary condition of the flow through the coaxial cylinders. Temperature, concentration force, and velocity profile are discussed in relation to various parameters using graphs. The variation of streamline wave frames for some significance parameters is also displayed in graphs. A comparison analysis with previously published graphical results [30] with the same geometry was also carried out in the conclusion section. The findings of the present analysis can help understand the flow mechanisms of many fluids with physiological and geographical interest.

2. Mathematical Formulation

The analysis considered the peristaltic flow of an incompressible Williamson electrically conducting fluid in coaxial cylinders. The cylindrical coordinate framework ( R , Z ) were utilized, where R is along the span of the cylinder, and Z corresponds to the center line of the cylinder, as shown in Figure 1. A uniform magnetic field B0 was applied and acted along the axis. The details of geometry with dimension is given in Table 1.
The equations used were:
H   ( Z ¯ ,   t ) = a + b S i n [ 2 π λ   ( Z ¯ c t ) ]
where a is the average radius of the undisturbed tube, b is the amplitude of the peristaltic wave, λ is the wavelength, c is the wave propagation speed, and t ¯ is time. In the fixed coordinates ( R ¯ , Z ¯ ) the flow between the two tubes is unsteady. It becomes steady in a wave frame (r, z) moving with the same speed as the wave in the Z-direction. The transformations between the two frames is:
r ¯ = R ¯   ,   z ¯ = Z ¯ c t ¯ ,                          
u ¯ = U ¯   ,   w ¯ = W ¯ c ,
where ( r ¯ ,   z ¯ ) and ( U ¯ , W ¯ ) are the velocity components in the moving and fixed frames, respectively.
The extra stress tensor S ¯ for a Williamson fluid is [32,33]:
S ¯ = μ 0 [ ( 1 + Γ γ ˙ ¯ ) ] A 1
where γ ˙ ¯ = 1 2   , and = t r ( A 1 ) 2 . Here, A 1 = V ¯ + ( V ¯ ) T represents the first Rivlin–Ericksen tensors, Γ is the time constant, and μ 0 is the zero-shear-rate viscosity. The stress component can be obtained by using Equation (2) as follows:
S ¯ z ¯ r ¯ = μ 0 [ ( 1 + Γ γ ˙ ¯ ) ] [ u ¯ z ¯ + w ¯ r ¯ ]
S ¯ r ¯ r ¯ = μ 0 [ ( 1 + Γ γ ˙ ¯ ) ] [ 2 u ¯ r ¯ ]
S ¯ z ¯ z ¯ = μ 0 [ ( 1 + Γ γ ˙ ¯ ) ] [ 2 w ¯ z ¯ ]
The governing equations of the flow can be written as [18]:
u ¯ r ¯ + v ¯ r ¯ + w ¯ r ¯ = 0 ,
ρ f ( u ¯   [ u ¯ r ¯ ] + w ¯ u ¯ z ¯ ) = p ¯ r ¯ + 1 r ¯ r ¯ ( r ¯ S ¯ r ¯ r ¯ ) + z ¯ ( S ¯ z ¯ r ¯ ) + S ¯ θ ¯ θ ¯ r ¯
ρ f ( u ¯ w ¯ r ¯ + w ¯ w ¯ z ¯ ) = p ¯ z ¯ + 1 r ¯ r ¯ ( r ¯ S ¯ r ¯ z ¯ ) + z ¯ ( S ¯ z ¯ z ¯ ) + ρ f g β T ( T T 0 ) + ρ f g β C ( C C 2 ) σ B 0 2 S i n 2 ( ξ ) w ¯ μ k w ¯ ,
u ¯ T r ¯ + w ¯ T z ¯ = K c p ρ f ( 2 T r ¯ 2 + 1 r ¯ T r ¯ + 2 T z ¯ 2 ) + 16 σ 0 T 2 E 3 k 0 c p ρ f 1 r ¯ ( r ¯ T r ¯ ) Q c p ρ f T
u ¯ C r ¯ + w ¯ C z ¯ = D m ( 2 C r ¯ 2 + 1 r ¯ C r ¯ + 2 C z ¯ 2 ) + D m k T T m 1 r ¯ ( 2 T r ¯ 2 + 1 r ¯ T r ¯ + 2 T z ¯ 2 )
The temperature equation is given by:
T t ¯ + u ¯ T r ¯ + w ¯ T z ¯ = K c p ρ ( 2 T r ¯ 2 + 1 r ¯ T r ¯ + 2 T z ¯ 2 ) 16 σ 0 T 2 E 3 k 0 c p ρ 1 r ¯ r ¯ ( r ¯ T r ¯ ) Q c p ρ T
The concentration equation is given by:
C t ¯ + u ¯ C R ¯ + w ¯ C Z ¯ = D m ( 2 C R ¯ 2 + 1 R ¯ C R ¯ + 2 C Z ¯ 2 ) + D m k T T m ( 2 T R ¯ 2 + 1 R ¯ T R ¯ + 2 T Z ¯ 2 )
where the pressure is denoted by p ¯ , the viscosity is μ , the velocity in the r ¯ and z ¯ directions is by u ¯ and w ¯ , respectively, the density id ρ . The appropriate boundary conditions are:
w ¯ = 1 ,   u ¯ = 0 ,     T = T 1 ,   C = C 1     a t     r = r 1 = ε w ¯ = 1 ,   u ¯ = 0 ,     T = T 0   ,   C = C 0   a t     r = r 2 = 1 +   S i n ( 2 π z ¯ )   }
We introduced the dimensionless parameters as follows:
r = r ¯ a 2   ,   z = z ¯ λ   ,   δ = a 2 λ ,   u = λ u ¯ a 2 c ,   w = w ¯ c ,   p = a 2 2 p μ 0 λ c   ,   r 1 = r 1 ¯ a 2 = ε < 1 , r 2 = r 2 ¯ a 2 = 1 + . sin ( 2 π z ) ,   = b a 2 ,   θ = T T 0 T 1 T 0 ,       D a = K a 2 2 ,   s = a 2 s ¯ μ 0 c   , W e = c Γ a 2 Φ = C C 0 C 1 C 0   , R e = c p a 2 μ , γ ˙ = γ ˙ ¯ D m , S c = μ c p K , M 1 2 = σ B 0     2 μ a 2 2 s i n 2 ( α ) , R n = ρ K 0   C p v 4 T 2 2 σ 0 G r = ρ g β T h 2 ( T T 0 ) μ U , S r = D K T ( T 1 T 0 ) U T m h ( C 1 C 2 ) , G c = ρ g β C h 2 ( T T 0 ) μ U ,   P r = μ 0 c p k }
where ( R n ) the is the radiation parameter, ( ) is the amplitude ratio, ( M 1 2 ) is the magnetic parameter, ( R e ) is the Reynolds number, ( D a ) is the Darcy number, ( S r ) is the Soret number, ( S c ) is the Brandt number, and ( δ ) is the dimensionless wave number. ( G r ) is the thermal Grashof number, and ( G c ) is by Solutol Grashof number.
Using Equation (14) in Equations (6)–(12), the dimensionless equations can be written as:
( u r + u r + w z ) = 0 ,
R e   δ 3 ( u   u r + w u z ) = p r + δ 2   ( S r z   ) z + δ r   ( r S r r   ) r δ S ϑ ϑ r ,
R e δ (   u   w r + w w z ) = p z + ( S r z ) r + ( S r z   ) r + δ   ( S z z   ) z ( M 1 2 + 1 D a ) w + G r ϑ + G c φ ,
δ ( u   ϑ r + w ϑ z ) = 1 p r ( 2 ϑ r 2 + 1 r ϑ r + δ 2 2 ϑ z 2 ) + 4 3 R n 1 r r ( r ϑ r ) Ω ϑ ,
δ ( u   φ r + w φ z ) = 1 S c ( 2 φ r 2 + 1 r φ r + δ 2 2 φ z 2 ) + S r ( 2 ϑ r 2 + 1 r ϑ r + δ 2 2 ϑ z 2 ) ,
By assuming a long wavelength (δ << 1) and a low Reynolds number (Re → 0), one obtains:
p r = 0 ,
1 r ( r S r z   ) r ( M 1 2 + 1 D a ) w + G r θ + G c Φ = p z
( 1 R e p r + 4 3 R n ) ( 2 θ r 2 + 1 r θ r ) Ω θ = 0 ,
1 S c ( 2 Φ r 2 + 1 r Φ r ) = S r ( 2 θ r 2 + 1 r θ r ) ,
where
S r r = S θ θ = S z z = 0 , and   S r z = [ 1 + W e ( w r ) ] ( w r )
By replacing ( S r z ) in Equation (21), we have:
r ( [ ( w r ) W e ( w r ) + ( w r ) ] ) + 1 r [ ( w r ) + W e ( w r ) ( w r ) ] ( M 1 2 + 1 D a ) w + G r θ + G c Φ = p z ,
The dimensionless boundary conditions are the following:
w = 1 ,   u = 0 ,     ϑ = 1     at     r = r 1 = ε w = 1 ,   u = 0 ,     ϑ = 0   at     r = r 2 = 1 +   S i n ( 2 π z )   }
Considering the stream function ψ :
u = 1 r ( ψ z )   and   w = 1 r ( ψ r ) ,

3. Solution of the Problem

The dimensionless partial differential equations were solved using the Bessel functions of the first and second types, with the help of the software Mathematica (version-11). The Bessel functions can be applied on the following type of equation to obtain an exact solution
x 2   d 2 y d x 2 + x d y d x ( x 2 + p 2 ) y = 0
Comparing Equation (27) with Equations (22)–(24), we obtained temperature, concentration, and velocity profile, as shown here.
The solution of the temperature Equation (22) is:
θ = c 1 I 0 [ A r ] + c 2 B 0 [ A r ]
The constant are
  c 1 = B 0 [ h A ] I 0 ϵ [ [ A ] B 0 [ h A ] ( h A ) B 0 ( A ) ]
c 2 = I 0 [ h A ] I 0 ϵ [ ( h A ) B 0 ( A ) ( A ) B 0 ( h A ) ]
From the concentration Equation (23), we obtain
Φ = S c S r + c 3 ln ( r ) + c 4
With c 3 = 1 + S c S r L n   ( r 1 ( 1 / r 2 ) ) , and   c 4 = c 3 ln ( r 2 ) .
Equation (24) shows that p depends on z only. Equation (24) can be written as
2 w r 2 + 1 r w r + 1 r ( w r ) 2 ( M 1 2 + 1 D a ) = p z G c θ G r ϕ
With H = ( M 1 2 + 1 D a ) , and Y = ( p z G r ϑ G c φ + ( M 1 2 + 1 D a ) ) ,
The general solution of Equation (25) is
w = H Y + k 1 I 0 [ F r ] + k 2 B 0 [ F r ]
Considering the boundary conditions provided in Equation (25), we have
k 1 = ( H + Y ) ( B 0 [ h Y ] B 0 [ Y ϵ ] ) H ( I 0 [ Y ϵ ] B 0 [ h Y ] I 0 [ h Y ] B 0 [ Y ϵ ] )
k 2 = ( H + Y ) ( I 0 [ h Y ] + I 0 [ Y ϵ ] ) Y ( I 0 [ Z ϵ ] B 0 [ h Z ] I 0 [ h Y ] B 0 [ Y ϵ ] )
The modified Bessel functions of the first and second kind of zero order are I 0 , B 0 . The solution was obtained by using the boundary condition, which was provided in Equation (25); with the help of (MATHEMATICA-11) software, we obtained the constants c 1 ,   c 2 , c 3 ,   c 4 ,   k 1 and k 2 .

4. Results and Discussion

4.1. Concentration Profile

Figure 2 shows the effects of the parameters Sc and Sr on the concentration profile Φ . It was observed that Φ decreases with the increase on Sr with r in the range from 0 to 1.37622, beyond which, it increased with the increase of Sr. The curves of Φ were obtained for S r = 0.1, 0.2, and 0.3 when the other parameters were constant. Figure 2b shows the effect of S c on the concentration distribution, and it is clear that Φ decreased as S c increased, and also S c affected the relation between Φ and r. This relation appeared approximately inverse for small values of Sc, but, for large values of Sc, Φ increased with increasing, till a finite value of r (minimum value), after which it decreased. These phenomena are due to the fact that the solute diffusion in fluids is always proportional to the diffusion coefficient. Therefore, a decrease in the concentration is due to a decrease of the diffusion coefficient. The same behavior was shown for S c in Figure 2b. It can be explained by Equation (22), since both parameter are inversely proportional to the concentration. As a consequence for it recuses for higher values of Sr and Sc.

4.2. Temperature Profile

Figure 3 shows the effects of the parameters P r , R e , and R n on the temperature profile. The graphs show that θ increased with the increase in any one of these parameters, while the others were fixed. It should be noted that for higher values of the Prandtl number, thermal conductivity was lower, which caused a reduction in conduction and boundary layer thickness. Further, the temperature increased when the radiation parameter increased. This was due to fact that the surface heat flux became large under the influence of thermal radiation, which resulted in higher temperature inside the boundary layer region. In Figure 3, the effect of R e , R n , and P r on the temperature distribution T is displayed, and this figure shows that the temperature increased with the increasing of these parameters. The relation between T and r was approximately linear for large values of P r ,   R e ,   and   R n > 9 , but for large values of Pr (>0.7), T increased with increasing r till a finite value of r (maximum value), after which T decreased with increasing r.

4.3. Velocity Profile

The effect of Grashof number ( G c ) is illustrated in Figure 4. It was observed that the fluid velocity ( W ) decreased with the increase of ( G c ) , and started to increase when r   was   2.27 . It is seen that with the increase of ( G r ) , the velocity ( W ) declined, but W became larger for r > 2.1726 . Figure 5 shows the variation of the velocity ( W ) with ( r ) for various values of the Weissenberg number ( W e ) . This figure shows that the fluid velocity increased with the increase of the Weissenberg number in the range of [0.1–0.32], whereas the opposite was observed in the range of [0.32–0.60]. To examine a viscoelastic flow, the Weissenberg number was used. We inspected the effect of forces from elastic to viscous. For high values of the Weissenberg number, the flow of fluid particles encountered resistance, and as a consequence, the velocity decreased. In Figure 5 we see that ( W ) decreased with the increase of ( ε ) . Figure 6 illustrates the velocity profile for various values of the magnetic and Darcy parameters. We concluded that the velocity decreased for increasing values of the magnetic field parameter ( M ) . The reason for this phenomenon is that, due to the increase in magnetic the parameters, the Lorentz force was generated and created a resistance force in the flow; as a result, the velocity profile decreased. Further, this figure shows that the velocity decreased by increasing M in the range of 0.10 < r < 0.27 ; ad after that, the velocity increased by increasing M. It is also clear that, in all curves in this figure, the maximum value of ( W ) occurred at 0.23 < r < 0.25 . In Figure 6, the effect of D a on the velocity distribution is displayed. It is clear that ( W ) increased by increasing D a in the range of 0.10 < r < 0.22 , and ( W ) decreased by increasing D a in the range of 0.22 < r < 0.30 . In this figure, we can also see that, although there was a big difference between the values of D a , we observed small differences between the three curves of ( W ) , which were taken at D a = 0.1 ,   0.5 ,   0.9 , when the other parameters were constant.
The effects of Grashof number ( G c ) , Solutol Grashof number ( G r ) , Weissenberg number ( W e ) , Darcy number ( D a ) , magnetic field parameter ( M ) , and ( ε ) on trapping can be seen in Figure 7, Figure 8, Figure 9, Figure 10 and Figure 11. Trapping is basically the development of an interior circulating bolus of fluid by closed streamlines that is pushed ahead along with the peristaltic wave.
Figure 7 shows that the size of the trapped bolus increased with the increase in G r , and near the upper wall, it disappeared when G r = 4.1 . Figure 8 shows that a new region of the trapped bolus vanished at G c   = 3.92 nearby the flat wall of the channel, while the trapped bolus size increased for a large value of G c .
It is observed in Figure 9 that the size of the trapped bolus declined with the enhancement of the magnetic parameter M , while at M   = 4.82 , a new trapped bolus appeared close to the flat wall of the channel and increased in dimension with the increase of M . Figure 10 reveals that close to the upper wall at D a = 3.99 , the trapped bolus vanished into a wave, while the size of the trapped bolus increased with the increase of D a .
It appears in Figure 11 that the size of the trapped bolus became larger for larger values of the parameter W e , and at W e = 3.9 , the trapped bolus disappeared into wave near the region of the upper wall.

Validation of the Model

The presented mathematical model was developed for non-Newtonian Williamson fluids and can be reduced to non-Newtonian fluids as W e   0 . In that case, we used the same model and mathematical equation previously described [30,31] for non-Newtonian fluids. For this purpose, the velocity profile for different values of magnetic parameter and Darcy number is presented in Figure 12. It can be seen from this comparison graph that the velocity decreased as M increased in the range of 0.5 < r < 1 . The effect of D a on velocity distribution was similar. It is clear that w decreased as D a increased in the range of 0.5 < r < 1 , and w decreased as D a increased. Both results are similar to those obtained by Shabaan [30].

5. Conclusions

The magnetohydrodynamic peristaltic flow of a Williamson fluid for varying temperature, velocity, and concentration values through a porous medium was considered in this study. It is important to mention that if the current Williamson model is used in rectangular coordinates, we obtain the same results as those previously reported [23,24]. The following graphical results that we obtained with the described model are the same those as obtained by Shaban et al. [21].
  • The behavior of the velocity profile was the same as the values of magnetic parameter and Da parameter increased.
  • The temperature distribution was same for higher values of the Prandtl and M parameters.
  • The dimensionless concentration was the same for different values of Sc and Sr
The essential features of this analysis are the following:
  • For higher values of boundary layer thickness and thermal conductivity, the temperature θ increased; the relation between θ and r was approximately linear for large values of P r ,   R e ,   and   R n > 9 .
  • The velocity decreased by increasing M in the range of 0.10 < r < 0.27 ; after that, the velocity increased by increasing M.
  • The velocity decreased for higher values of the parameters G r and G c , but at r = 0.20 , it started to increase.
  • The concentration force showed a dual behavior for various values of the parameters Sc and S r because solute diffusion in fluids is always proportional to the diffusion coefficient. Therefore, a decrease in concentration field was due to a decrease in the diffusion coefficient.
  • Similar effects were observed in trapping for the parameter G r , and the size of the trapped bolus decreased with the increase of G r . It is further concluded that a trapped bolus developed in a new region near the flat wall of the channel and increased in size with the increase of G r .
  • The performance of the parameters D a ,   M ,   G c , and W e was similar on the trap and it was revealed that near the trap, the bolus increased with the increase of W .

Author Contributions

Conceptualization, L.C. and W.S.K.; methodology, M.A.A.; software, M.A.A.; validation, B.S.; formal analysis, L.C.; investigation, W.S.K.; resources, B.S.; data curation, B.S.; writing—original draft preparation, M.A.A.; writing—review and editing, M.A.A.; visualization, B.S.; supervision, L.C.; project administration W.S.K.; funding acquisition, L.C. and B.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by [Shanghai Major Science Popularization] grant number [20DZ2306500] And Special Project for Capacity Improvement of Shanghai Professional Technical Service Platform, grant number [19DZ2290400].

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

P r Prandtl number
ρ Density
p Pressure
T ˜ Temperature
σ Electrical conductivity
B 0 Magnetic field
C p Specific heat
KThermal conductivity
GcSolutol Grashof number
U, WExpression of velocity in r and z directions, Respectively
G r Thermal Grashof number
θ Dimensionless temperature of the model
S r Soret number
δ Dimensionless wave number
Da Darcy Number
μ o Zero share stress viscosity
ϕ Concentration force
M Dimensionless magnetic parameter
μ Fluid viscosity
Ω Stagnation speed
D a Darcy number
λ 1 The ratio of relaxation to retardation time
tTime
φ Amplitude rate
RnRadiation parameter
ω Frequency of oscillation
Γ Time Constant Parameter
A 1 Rivlin Ericksen Tensor
WeWeissenberg number
ScSuction/Injection parameter
ϵ Variation of viscosity with temperature
ReReynolds number

References

  1. Bhatti, M.M.; Abbas, M.A.; Rashidi, M.M. Entropy generation in blood flow with heat and mass transfer for the Ellis fluid model. Heat Transf. Res. 2018, 49, 747–760. [Google Scholar] [CrossRef]
  2. Bhatti, M.M.; Abbas, M.A.; Rashidi, M.M. Entropy generation for peristaltic blood flow with casson model and consideration of magnetohydrodynamics effects. Walailak J. Sci. Technol. (WJST). 2017, 14, 451–461. [Google Scholar]
  3. Raza, R.; Mabood, F.; Naz, R.; Abdelsalam, S.I. Thermal transport of radiative Williamson fluid over stretchable curved surface. Therm. Sci. Eng. Prog. 2021, 23, 100887. [Google Scholar] [CrossRef]
  4. Jia, Z.; Xuan, Y. Analysis on interaction among solar light and suspended nanoparticles in nanofluids. J. Quant. Spectrosc. Radiat. Transf. 2021, 269, 107692. [Google Scholar]
  5. Sun, J.; Guo, L.; Jing, J.; Tang, C.; Lu, Y.; Fu, J.; Ullmann, A.; Brauner, N. Investigation on laminar pipe flow of a non-Newtonian Carreau-Extended fluid. J. Pet. Sci. Eng. 2021, 205, 108915. [Google Scholar] [CrossRef]
  6. Rashid, M.; Ansar, K.; Nadeem, S. Effects of induced magnetic field for peristaltic flow of Williamson fluid in a curved channel. Phys. A Stat. Mech. Appl. 2020, 553, 123979. [Google Scholar] [CrossRef]
  7. Nadeem, S.; Akram, S. Peristaltic flow of a Williamson fluid in an asymmetric channel. Communications in Nonlinear Sci. Num. Simul. 2010, 15, 1705–1716. [Google Scholar] [CrossRef]
  8. Almaneea, A. Numerical study on heat and mass transport enhancement in MHD Williamson fluid via hybrid nanoparticles. Alex. Eng. J. 2022, 61, 8343–8354. [Google Scholar] [CrossRef]
  9. Dadheech, A.; Parmar, A.; Agrawal, K.; Al-Mdallal, Q.; Sharma, S. Second law analysis for MHD slip flow for Williamson fluid over a vertical plate with Cattaneo-Christov heat flux. Case Stud. Therm. Eng. 2022, 33, 101931. [Google Scholar] [CrossRef]
  10. Asjad, M.I.; Zahid, M.; Inc, M.; Baleanu, D.; Almohsen, B. Impact of activation energy and MHD on Williamson fluid flow in the presence of bioconvection. Alex. Eng. J. 2022, 61, 8715–8727. [Google Scholar] [CrossRef]
  11. Shashikumar, N.; Madhu, M.; Sindhu, S.; Gireesha, B.; Kishan, N. Thermal analysis of MHD Williamson fluid flow through a microchannel. Int. Commun. Heat Mass Transf. 2021, 127, 105582. [Google Scholar] [CrossRef]
  12. Bhatti, M.M.; Rashidi, M.M.; Abbas, M.A. Analytic Study of Drug Delivery in Peristaltically Induced Motion of Non-Newtonian Nanofluid. J. Nanofluids 2016, 5, 920–927. [Google Scholar] [CrossRef]
  13. Abbas, M.A.; Bai, Y.; Rashidi, M.M.; Bhatti, M.M. Analysis of Entropy Generation in the Flow of Peristaltic Nanofluids in Channels with Compliant Walls. Entropy 2016, 18, 90. [Google Scholar] [CrossRef]
  14. Sinha, A.; Shit, G.C.; Ranjit, N.K. Peristaltic transport of MHD flow and heat transfer in an asymmetric channel: Effects of variable viscosity, velocity-slip and temperature jump. Alex. Eng. J. 2015, 54, 691–704. [Google Scholar] [CrossRef] [Green Version]
  15. Saleem, A.; Akhtar, S.; Alharbi, F.M.; Nadeem, S.; Ghalambaz, M.; Issakhov, A. Physical aspects of peristaltic flow of hybrid nano fluid inside a curved tube having ciliated wall. Results Phys. 2020, 19, 103431. [Google Scholar] [CrossRef]
  16. Riaz, A.; Bhatti, M.M.; Ellahi, R.; Zeeshan, A.; Sait, S.M. Mathematical Analysis on an Asymmetrical Wavy Motion of Blood under the Influence Entropy Generation with Convective Boundary Conditions. Symmetry 2020, 12, 102. [Google Scholar] [CrossRef] [Green Version]
  17. Abbas, M.A.; Bég, O.A.; Zeeshan, A.; Hobiny, A.; Bhatti, M.M. Parametric analysis and minimization of entropy generation in bioinspired magnetized non-Newtonian nanofluid pumping using artificial neural networks and particle swarm optimization. Therm. Sci. Eng. Prog. 2021, 24, 100930. [Google Scholar] [CrossRef]
  18. Afifi, N.A.S. Study of Peristaltic Flow for Different Cases. Ph.D. Thesis, Ain Shams University, Cairo, Egypt, 1998. [Google Scholar]
  19. Bhatti, M.M.; Abbas, M.A. Simultaneous effects of slip and MHD on peristaltic blood flow of Jeffrey fluid model through a porous medium. Alex. Eng. J. 2016, 55, 1017–1023. [Google Scholar] [CrossRef] [Green Version]
  20. Noreen, S.; Tripathi, D. Heat transfer analysis on electroosmotic flow via peristaltic pumping in non-Darcy porous medium. Therm. Sci. Eng. Prog. 2019, 11, 254–262. [Google Scholar] [CrossRef]
  21. Pattnaik, P.K.; Parida, S.K.; Mishra, S.R.; Abbas, M.A.; Bhatti, M.M. Analysis of Metallic Nanoparticles (Cu, Al2O3, and SWCNTs) on Magnetohydrodynamics Water-Based Nanofluid through a Porous Medium. J. Math. 2022, 2022, 3237815. [Google Scholar] [CrossRef]
  22. Hareli, S.; Nave, O.; Gol’dshtein, V. The Evolutions in Time of Probability Density Functions of Poly dispersed Fuel Spray—The Continuous Mathematical Model. Appl. Sci. 2021, 11, 9739. [Google Scholar] [CrossRef]
  23. Yang, M.; Abbas, M.A.; Khudair, W.S. Energy and Temperature-Dependent Viscosity Analysis on Magnetized Eyring-Powell Fluid Oscillatory Flow in a Porous Channel. Energies 2021, 14, 7829. [Google Scholar] [CrossRef]
  24. Nawaz, M.; Sadiq, M.A. Unsteady heat transfer enhancement in Williamson fluid in Darcy-Forchheimer porous medium under non-Fourier condition of heat flux. Case Stud. Therm. Eng. 2021, 28, 101647. [Google Scholar] [CrossRef]
  25. Hayat, T.; Saleem, A.; Tanveer, A.; Alsaadi, F. Numerical study for MHD peristaltic flow of Williamson nanofluid in an endoscope with partial slip and wall properties. Int. J. Heat Mass Transf. 2017, 114, 1181–1187. [Google Scholar] [CrossRef]
  26. Bhatti, M.M.; Arain, M.B.; Zeeshan, A.; Ellahi, R.; Doranehgard, M.H. Swimming of Gyrotactic Microorganism in MHD Williamson nanofluid flow between rotating circular plates embedded in porous medium: Application of thermal energy storage. J. Energy Storage 2022, 45, 103511. [Google Scholar] [CrossRef]
  27. Li, Y.-X.; Alshbool, M.H.; Lv, Y.-P.; Khan, I.; Khan, M.R.; Issakhov, A. Heat and mass transfer in MHD Williamson nanofluid flow over an exponentially porous stretching surface. Case Stud. Therm. Eng. 2021, 26, 100975. [Google Scholar] [CrossRef]
  28. Ahmed, K.; Akbar, T.; Muhammad, T.; Alghamdi, M. Heat transfer characteristics of MHD flow of Williamson nanofluid over an exponential permeable stretching curved surface with variable thermal conductivity. Case Stud. Therm. Eng. 2021, 28, 101544. [Google Scholar] [CrossRef]
  29. Salmi, A.; Madkhali, H.A.; Nawaz, M.; Alharbi, S.O.; Alqahtani, A. Numerical study on non-Fourier heat and mass transfer in partially ionized MHD Williamson hybrid nanofluid. Int. Commun. Heat Mass Transf. 2022, 133, 105967. [Google Scholar] [CrossRef]
  30. Shaaban, A.A.; Abou-Zeid, M.Y. Effects of Heat and Mass Transfer on MHD Peristaltic Flow of a Non-Newtonian Fluid through a Porous Medium between Two Coaxial Cylinders. Math. Probl. Eng. 2013, 2013, 819683. [Google Scholar] [CrossRef]
  31. Eldabe, N.T.; Mohamed, Y.A. Magnetohydrodynamic peristaltic flow with heat and mass transfer of mi-cropolar biviscosity fluid through a porous medium between two co-axial tubes. Arab. J. Sci. Eng. 2014, 39, 5045–5062. [Google Scholar] [CrossRef]
  32. Khudair, W.S.; Dheia Gaze, S.A. Influence of Heat Transfer on MHD Oscillatory Flow for Williamson Fluid with Variable Viscosity through a Porous Medium. Int. J. Fluid Mech. Therm. Sci. 2018, 4, 11–17. [Google Scholar] [CrossRef] [Green Version]
  33. Wissam, S.K.; Al-Khafajy, D.G. Influence of heat transfer on Magneto hydrodynamics oscillatory flow for Williamson fluid through a porous medium. Iraqi J. Sci. 2018, 59, 389–397. [Google Scholar]
Figure 1. Geometry of the flow.
Figure 1. Geometry of the flow.
Symmetry 14 00953 g001
Figure 2. The variation of the concentration Φ vs. r at Ω = 1 , R e = 3 , P r = 2 , R n = 0.5 , = 0.3 , ϵ = 0.3 , i = 1 , z = 0.1 for: (a) S c = 0.1 (b) S r = 0.07 .
Figure 2. The variation of the concentration Φ vs. r at Ω = 1 , R e = 3 , P r = 2 , R n = 0.5 , = 0.3 , ϵ = 0.3 , i = 1 , z = 0.1 for: (a) S c = 0.1 (b) S r = 0.07 .
Symmetry 14 00953 g002
Figure 3. The variation of the temperature θ vs. r, at Ω = 1 , P r = 2 ,   ϵ = 0.3 , = 0.3 , z = 0.1 . (a) Pr = 2, Rn = 2 (b) Re = 0.1, Pr = 2, (c) Re = 0.1, Rn = 2.
Figure 3. The variation of the temperature θ vs. r, at Ω = 1 , P r = 2 ,   ϵ = 0.3 , = 0.3 , z = 0.1 . (a) Pr = 2, Rn = 2 (b) Re = 0.1, Pr = 2, (c) Re = 0.1, Rn = 2.
Symmetry 14 00953 g003
Figure 4. Velocity profile for G c and G r with Ω = 1 , ϵ = 0.2 , ϕ = 0.5 , W e = 1 , R e = 2 , P r = 2 , R n = 2 , p z = 1.5 , S c = 0.5 ,   S r = 0.07 , M = 1 , D a = 0.8 , ξ = π 4 .
Figure 4. Velocity profile for G c and G r with Ω = 1 , ϵ = 0.2 , ϕ = 0.5 , W e = 1 , R e = 2 , P r = 2 , R n = 2 , p z = 1.5 , S c = 0.5 ,   S r = 0.07 , M = 1 , D a = 0.8 , ξ = π 4 .
Symmetry 14 00953 g004
Figure 5. Velocity profile for W e and ϵ with Ω = 1 , G c = 1 , ϕ = 0.5 , G r = 1 , R e = 2 , P r = 2 , R n = 2 , p z = 1.5 , S c = 0.5 ,   S r = 0.07 , M = 1 , D a = 0.8 , ξ = π 4 .
Figure 5. Velocity profile for W e and ϵ with Ω = 1 , G c = 1 , ϕ = 0.5 , G r = 1 , R e = 2 , P r = 2 , R n = 2 , p z = 1.5 , S c = 0.5 ,   S r = 0.07 , M = 1 , D a = 0.8 , ξ = π 4 .
Symmetry 14 00953 g005
Figure 6. Velocity profile for M and D a with Ω = 1 , ϵ = 0.2 , ϕ = 0.5 , W e = 1 , R e = 2 , P r = 2 , R n = 2 , p z = 1.5 , S c = 0.5 , S r = 0.07 , G c = 1 , G r = 1 , ξ = π 4 .
Figure 6. Velocity profile for M and D a with Ω = 1 , ϵ = 0.2 , ϕ = 0.5 , W e = 1 , R e = 2 , P r = 2 , R n = 2 , p z = 1.5 , S c = 0.5 , S r = 0.07 , G c = 1 , G r = 1 , ξ = π 4 .
Symmetry 14 00953 g006
Figure 7. Streamlines in the wave frame for various values of G r when Ω = 0.9 , ϵ = 0.3 , ϕ = 0.1 , W e = 0.1 , R e = 1 , R n = 2 , P r = 2 , S c = 0.5 , S r = 0.1 , G c = 0.1 , M = 1.1 , D a = 0.9 , ξ = π 4 ; ( a )   G r = 0.1 ,   ( b )   G r = 0.2 .
Figure 7. Streamlines in the wave frame for various values of G r when Ω = 0.9 , ϵ = 0.3 , ϕ = 0.1 , W e = 0.1 , R e = 1 , R n = 2 , P r = 2 , S c = 0.5 , S r = 0.1 , G c = 0.1 , M = 1.1 , D a = 0.9 , ξ = π 4 ; ( a )   G r = 0.1 ,   ( b )   G r = 0.2 .
Symmetry 14 00953 g007
Figure 8. Streamlines in the wave frame for various values of G c when ϕ = 0.1 , ϵ = 0.3 , W e = 0.1 ,   Ω = 0.9 , R e = 1 , R n = 2 , P r = 2 , S c = 0.5 , S r = 0.1 , G r = 2 , M = 1.1 , D a = 0.9 , ξ = π 4 ; ( a )   G c = 3 ,   ( b )   G c = 3.5 .
Figure 8. Streamlines in the wave frame for various values of G c when ϕ = 0.1 , ϵ = 0.3 , W e = 0.1 ,   Ω = 0.9 , R e = 1 , R n = 2 , P r = 2 , S c = 0.5 , S r = 0.1 , G r = 2 , M = 1.1 , D a = 0.9 , ξ = π 4 ; ( a )   G c = 3 ,   ( b )   G c = 3.5 .
Symmetry 14 00953 g008
Figure 9. Streamlines in the wave frame for various values of M when ϕ = 0.1 , ϵ = 0.3 , W e = 0.1 , R e = 1 , R n = 2 , P r = 2 , S c = 0.5 , S r = 0.1 , G c = 0.1 , G r = 2 , D a = 0.9 , ξ = π 4 ; ( a )   M = 0.1 ,   ( b )   M = 0.2 .
Figure 9. Streamlines in the wave frame for various values of M when ϕ = 0.1 , ϵ = 0.3 , W e = 0.1 , R e = 1 , R n = 2 , P r = 2 , S c = 0.5 , S r = 0.1 , G c = 0.1 , G r = 2 , D a = 0.9 , ξ = π 4 ; ( a )   M = 0.1 ,   ( b )   M = 0.2 .
Symmetry 14 00953 g009
Figure 10. Streamlines in the wave frame for various values of W e when ϵ = 0.3 , M = 1.1 , R e = 1 , R n = 2 , Ω = 0.9 , P r = 2 , ϕ = 0.1 , S c = 0.5 , S r = 0.1 , D a = 0.9 , G r = 2 , G c = 0.1 , ξ = π 4 ; ( a )   W e = 0.3 ,   ( b )   W e = 0.4 .
Figure 10. Streamlines in the wave frame for various values of W e when ϵ = 0.3 , M = 1.1 , R e = 1 , R n = 2 , Ω = 0.9 , P r = 2 , ϕ = 0.1 , S c = 0.5 , S r = 0.1 , D a = 0.9 , G r = 2 , G c = 0.1 , ξ = π 4 ; ( a )   W e = 0.3 ,   ( b )   W e = 0.4 .
Symmetry 14 00953 g010
Figure 11. Streamlines in the wave frame for various values of ϵ when D a = 0.9 , R n = 2 , ϕ = 0.1 , W e = 0.1 , Ω = 0.9 , R e = 1 , P r = 2 , S c = 0.5 , S r = 0.1 ,   G r = 2 , G r = 2 , M = 1 , ξ = π 4 ; ( a )   ϵ = 0.2 ,   ( b )   ϵ = 0.3 .
Figure 11. Streamlines in the wave frame for various values of ϵ when D a = 0.9 , R n = 2 , ϕ = 0.1 , W e = 0.1 , Ω = 0.9 , R e = 1 , P r = 2 , S c = 0.5 , S r = 0.1 ,   G r = 2 , G r = 2 , M = 1 , ξ = π 4 ; ( a )   ϵ = 0.2 ,   ( b )   ϵ = 0.3 .
Symmetry 14 00953 g011
Figure 12. Graphical results of the Velocity Profile against Magnetic parameter and Darcy number for a non-Newtonian fluid ( W e   0 ).
Figure 12. Graphical results of the Velocity Profile against Magnetic parameter and Darcy number for a non-Newtonian fluid ( W e   0 ).
Symmetry 14 00953 g012
Table 1. Details of the flow geometry.
Table 1. Details of the flow geometry.
ParametersDimensionsValues
r 1 Radius of the inner cylinder a 2 2
r 2 Radius of the outer cylinder H   ( Z ¯ ,   t ) = a + b S i n [ 2 π λ   ( Z ¯ c t ) ]
c = 2 ε 0   λ ln   ( r 2 r 1 ) 2   2 ( 8.85 × 10 2 ) ( 10 × 10 2 ) l n 4 m m 2 m m Capacitance of the Cylinder
B 0 Applied magnetic field
a r 1 + r 2 2 Average radius
R, Z(Cylindrical coordinates)R is the radical direction, and Z lies along the center line of
the inner and outer tubes
λ Length 10 × 10 2 m
bAmplitude of the peristaltic waveAssumed as a constant
C Wave propagation
t ¯ Time
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Chen, L.; Abbas, M.A.; Khudair, W.S.; Sun, B. Analytical Solution for the MHD Flow of Non-Newtonian Fluids between Two Coaxial Cylinders. Symmetry 2022, 14, 953. https://doi.org/10.3390/sym14050953

AMA Style

Chen L, Abbas MA, Khudair WS, Sun B. Analytical Solution for the MHD Flow of Non-Newtonian Fluids between Two Coaxial Cylinders. Symmetry. 2022; 14(5):953. https://doi.org/10.3390/sym14050953

Chicago/Turabian Style

Chen, Li, Munawwar Ali Abbas, Wissam Sadiq Khudair, and Bo Sun. 2022. "Analytical Solution for the MHD Flow of Non-Newtonian Fluids between Two Coaxial Cylinders" Symmetry 14, no. 5: 953. https://doi.org/10.3390/sym14050953

APA Style

Chen, L., Abbas, M. A., Khudair, W. S., & Sun, B. (2022). Analytical Solution for the MHD Flow of Non-Newtonian Fluids between Two Coaxial Cylinders. Symmetry, 14(5), 953. https://doi.org/10.3390/sym14050953

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop