Next Article in Journal
Distributed Event-Triggered Control for Manipulator with Fixed-Time Disturbance Observer
Previous Article in Journal
POI Route Recommendation Model Based on Symmetrical Naive Bayes Classification Spatial Accessibility and Improved Cockroach Swarm Optimization Algorithm
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

MHD Time-Periodic Plane Poiseuille Flow of Generalized Burgers Fluids through a Porous Medium

by
Constantin Fetecau
1,*,
Shehraz Akhtar
2 and
Costică Moroşanu
3
1
Section of Mathematics, Academy of Romanian Scientists, 3 Ilfov, 050044 Bucharest, Romania
2
Department of Mathematics, The Islamia University of Bahawalpur, Rahim Yar Khan Campus, Rahim Yar Khan 64200, Pakistan
3
Department of Mathematics, “Alexandru Ioan Cuza” University, 700506 Iasi, Romania
*
Author to whom correspondence should be addressed.
Symmetry 2024, 16(4), 425; https://doi.org/10.3390/sym16040425
Submission received: 15 February 2024 / Revised: 5 March 2024 / Accepted: 1 April 2024 / Published: 3 April 2024
(This article belongs to the Section Mathematics)

Abstract

:
The time-periodic plane Poiseuille flow of electrically conducting incompressible generalized Burgers fluids through a porous medium is analytically and numerically investigated in the presence of a transverse uniform magnetic field. The main purpose is to provide analytical expressions for the dimensionless steady-state fluid velocity, non-trivial shear stress and Darcy’s resistance, which can be used to bring to light important characteristics concerning fluid behavior. Similar solutions corresponding to the Poiseuille flow of the same fluids induced by a constant pressure gradient are obtained as limiting cases of previous results. The present results reduce to known solutions from the literature when magnetic and porous effects are neglected, and their validation is graphically proved. The needed time to reach a steady state has been graphically determined. It was found that the steady state is later obtained in the absence of a magnetic field and porous medium. The impact of a magnetic field and porous medium on the fluid velocity, shear stress and flow resistance has been systematically examined and elucidated through graphical representations. The findings reveal that the presence of a magnetic field or porous medium results in a reduction in the fluid velocity, accompanied by an increase in the flow resistance.

1. Introduction

The channel flow, as it was named by Coleman et al. [1], is the flow of a fluid between two infinite parallel plates that are both at rest. Such a flow, induced by the sudden application of a constant pressure gradient or due to a time-periodic pressure gradient, was called the Poiseuille flow or time-periodic plane Poiseuille flow by Spurk and Axel [2] and Rajagopal [3], respectively. Analytical solutions for the velocity field corresponding to motions of the incompressible Newtonian fluids produced by constant and time-dependent pressure gradients have been established by Mendiburu et al. [4] and Erdogan [5]. The first exact steady-state or long-time solutions for the time-periodic plane Poiseuille flows of non-Newtonian fluids seem to be those of Rajagopal [3], corresponding to incompressible second-grade fluids. These solutions were later extended by Hayat et al. [6] to incompressible Oldroyd-B fluids and by Fetecau et al. [7] for incompressible generalized Burgers fluids (IGBFs). Exact analytical solutions for such motions between parallel plates have been determined by Danish et al. [8] for third-grade fluids. Interesting results regarding the unsteady Poiseuille flow of second-grade fluids with slip boundary conditions have been recently obtained by Baranovskii [9].
Magnetohydromagnetic (MHD) motions of fluids between parallel plates have significant applications to polymer technology, the petroleum industry, power generators and pumps, geothermal systems, electrostatic precipitation and many others. The interplay between a moving electrically conducive fluid and a magnetic field generates effects with different applications in chemistry, physics and engineering. The MHD effects on the formation of such a motion have been investigated early by Tao [10] and Katagiri [11]. Singh and Ram [12] studied the MHD flow of incompressible Newtonian fluids through a channel induced by a time-periodic pressure gradient using Laplace transform. Steady solutions for MHD Couette flow between infinite parallel plates were relatively recently obtained by Kiema et al. [13] using the Sumudu transform. At the same time, the dynamics of fluid motion through porous media finds varied applications in the fields of geophysics and astrophysical studies, petroleum industry, oil reservoir technology and agricultural engineering. Relatively recently, Safdar Rabia et al. [14] studied the MHD flow of Burgers fluids due to a pressure gradient through a porous medium. Other exact solutions for MHD motions of Burgers fluids through a porous medium have been obtained by Fetecau and Vieru [15]. Unfortunately, such solutions for MHD motions of IGBFs induced by a periodic or constant pressure gradient through a porous plane channel seem to be lacking in the literature.
The purpose of this work is to provide analytical expressions for the steady-state solutions corresponding to MHD motions of IGBFs induced by a periodic or constant pressure gradient in a porous plane channel. These solutions can be easily particularized to give similar solutions for incompressible Burgers, Oldroyd-B, Maxwell and Newtonian fluids performing the same motions. For validation purposes, the steady-state velocity fields are presented in two distinct forms, and their equivalence is established through graphical representations. The convergence of initial solutions (numerical solutions) to their steady-state components is graphically validated, and the requisite time to attain the steady state is determined. Additionally, the impact of a magnetic field and porous medium on fluid velocity, shear stress and flow resistance is underscored and comprehensively discussed. It is worth mentioning the fact that the profiles of these entities against z are symmetric with respect to the median plane.

2. Problem Presentation

Let us consider an electrically conducive, incompressible, generalized Burgers fluid (ECIGBF) at equilibrium, confined between two infinitely expansive horizontal parallel plates within a porous medium. The Cauchy stress tensor T and the extra-stress tensor T E corresponding to these rate-type fluids are mandated to satisfy to the following constitutive equations [7]:
T = p I + T E ,   ( 1 + a 1 D D t + a 2 D 2 D t 2 ) T E = 2 η ( 1 + a 3 D D t + a 4 D 2 D t 2 ) D   ,
where I is the identity tensor, D is the rate of deformation tensor [1], p is the hydrostatic pressure, η is the dynamic viscosity of the fluid, D/Dt denotes the time upper-convected derivative and a i   ( i = 1 , 2 , 3 , 4 ) are dimensional material constants. If a 4 = 0 ,   a 4 = a 2 = 0 or a 4 = a 3 = a 2 = 0 in the second part of Equation (1), the corresponding equations characterize other classes of rate type fluids, namely, Burgers, Oldroyd-B and Maxwell fluids, respectively. Burgers fluids, whose frame-indifferent three-dimensional constitutive equations have been developed by Krishnan and Rajagopal [16], correspond to a 4 = 0 . Newtonian fluids are obtained when all constants are zero.
Let us assume that at the moment t = 0 + , a time-periodic pressure gradient
p z = ρ P cos ( ω t )   or   p z = ρ P sin ( ω t ) ;   t > 0 ,
is applied along the y-direction of a suitable Cartesian coordinate system x, y and z whose z-axis is perpendicular to the plates (see Figure 1).
Here, ρ is the fluid density, ω is the frequency of oscillations and P is constant. At the same time, an external uniform magnetic field of magnitude B 0 begins to act vertically to the plates. We suppose that the fluid is finitely conducting, the magnetic Reynolds number is small enough (at least ionized fluids and metallic liquids satisfy this condition [17]) and there is no surplus electric charge distribution. Consequently, the induced magnetic field and Joule heating can be neglected. In addition, the Hall effects are neglected due to moderate values of the Hartman number. Since the plates are unbounded, the velocity vector corresponding to the resulting flow is given by the next relation [18,19,20]:
u = u ( z , t ) = u ( z , t ) e y ;   0 < z < d ,   t > 0 ,
where e y is the unit vector along the y-axis and d is the distance between the plates. For such motions, the continuity equation is identically verified.
We also assume that the extra-stress tensor T E , as well as the velocity vector u, is a function of z and t only. Introducing the velocity vector u ( z , t ) from Equation (3) in the second relation (1), and bearing in mind the fact that the fluid has been at rest up to the initial moment t = 0 , it is not difficult to show that the non-trivial shear stress τ ( z , t ) = T E y z ( z , t ) has to satisfy the partial differential equation
( 1 + a 1 t + a 2 2 t 2 ) τ ( z , t ) = η ( 1 + a 3 t + a 4 2 t 2 ) u ( z , t ) z ;   0 < z < d ,   t > 0 .
The balance of linear momentum reduces to the partial differential equation [19]
ρ u ( z , t ) t = τ ( z , t ) z p z σ B 0 2 u ( z , t ) + R ( z , t ) = 0 ;   0 < z < d ,   t > 0 .
Here, σ is the electrical conductivity of the fluid, while R ( z , t ) is the Darcy’s resistance, which has to satisfy the partial differential equation [19]
( 1 + a 1 t + a 2 2 t 2 ) R ( z , t ) = η φ k ( 1 + a 3 t + a 4 2 t 2 ) u ( z , t ) ;   0 < z < 1 ,   t > 0 ,
in which φ ( 0 < φ < 1 ) is porosity and k ( > 0 ) represents the permeability of the porous medium. The volume flux Q ( t ) across a plane normal to the flow direction per unit width of this plane can be determined by means of the relation
Q ( t ) = 0 d u ( z , t ) d z ;   t > 0 .
The adequate initial and boundary conditions are
u ( z , 0 ) = 0 ,   u ( z , t ) t | t = 0 = 0 ,   2 u ( z , t ) t 2 | t = 0 = 0 ;   0 < z < d ,
respectively, when
u ( 0 , t ) = u ( d , t ) = 0 ;   t > 0 .
Utilizing the following non-dimensional variables, functions and parameters
z = 1 d z ,   t = P d ν t ,   u = ν P d 2 u ,   τ = 1 ρ P d τ ,   R = 1 ρ P R , Q = ν P d 3 Q ,   p = 1 ρ P d p ,   ω = ν P d ω  
and eliminating the star notation, one finds the non-dimensional forms
( 1 + α t + β 2 t 2 ) τ ( z , t ) = ( 1 + γ t + δ 2 t 2 ) u ( z , t ) z ;   0 < z < 1 ,   t > 0 ,
Re u ( z , t ) t = τ ( z , t ) z p z Ha 2 u ( z , t ) + R ( z , t ) = 0 ;   0 < z < 1 ,   t > 0 ,
( 1 + α t + β 2 t 2 ) R ( z , t ) = K ( 1 + γ t + δ 2 t 2 ) u ( z , t ) ;   0 < z < 1 ,   t > 0 ,
of the motion equations.
In Equations (10), (12) and (13), ν = η / ρ is the kinematic viscosity of the fluid, while the Reynolds and Hartmann numbers Re and Ha, respectively, and the porous parameter K is given by the next relations:
Re = V d ν ,   Ha 2 = σ B 0 2 ρ   d 2 ν = d 2 η σ B 0 2 = M ,   K = φ k d 2 ,
where V = P d 2 / ν is the characteristic velocity. In addition, the dimensionless constants α ,   β ,   γ and δ are given by the relations
α = P d ν a 1 ,   β = ( P d ν ) 2 a 2 ,   γ = P d ν a 3 ,   δ = ( P d ν ) 2 a 4 .  
The non-dimensional volume flux Q ( t ) is given by the relation
Q ( t ) = 0 1 u ( z , t ) d z ;   t > 0 .
Eliminating τ ( z , t ) and R ( z , t ) between Equations (11)–(13), one finds the following partial differential equation for the dimensionless velocity field u ( z , t ) :
Re ( 1 + α t + β 2 t 2 ) u ( z , t ) t = ( 1 + γ t + δ 2 t 2 ) 2 u ( z , t ) z 2 ( 1 + α t + β 2 t 2 ) p z Ha 2 ( 1 + α t + β 2 t 2 ) u ( z , t ) K ( 1 + γ t + δ 2 t 2 ) u ( z , t ) ;   0 < z < 1 ,   t > 0 ,
The appropriate initial and boundary conditions are
u ( z , 0 ) = 0 ,   u ( z , t ) t | t = 0 = 0 ,   2 u ( z , t ) t 2 | t = 0 = 0 ;   0 < z < 1 ,
respectively, for
u ( 0 , t ) = u ( 1 , t ) = 0 ;   t > 0 ,
while the relations in (2) take the dimensionless forms
p z = cos ( ω t )   or   p z = sin ( ω t ) ;   t > 0 .
In the next sections, for their distinction, we denote by u c ( z , t ) , τ c ( z , t ) , R c ( z , t ) and u s ( z , t ) , τ s ( z , t ) , R s ( z , t ) the dimensionless start-up solutions corresponding to the two motion problems of ECIGBFs. These solutions characterize the fluid behavior some time after its initiation. After this time, the fluid behavior can be described by the dimensionless steady-state (long time) solutions u c p ( z , t ) , τ c p ( z , t ) , R c p ( z , t ) or u s p ( z , t ) , τ s p ( z , t ) , R s p ( z , t ) . This is the time needed to reach the steady state. The two notions “steady-state solution” and “long time solution” have been used in the literature, for instance, by Erdogan [21] and Prusa [22], respectively, for similar solutions describing fluid motion over an infinite plate that oscillates in its plane. They are different from the steady solutions corresponding to the Poiseuille flow due to a constant pressure gradient (which are obtained in Section 3.1). To ascertain the time evolution of a given motion, it is adequate to discern the dimensionless steady-state solutions. Accordingly, the ensuing section will be dedicated to the derivation of analytical expressions encapsulating these solutions. Subsequently, these expressions will be employed to elucidate the impact of both the magnetic field and porous medium on the necessary time for attaining a steady state, as well as on the fluid’s flow resistance.

3. Analytical Expressions for the Dimensionless Steady-State Solutions

The dimensionless steady-state velocity fields u c p ( z , t ) and u s p ( z , t ) have to satisfy the governing Equation (17) and the boundary conditions (19). They are independent of the initial conditions (18), so that these conditions are ignored in what follows. In order to determine them and to be sure of their correctness, two different methods will be used. Firstly, let us use the steady-state complex velocity
u p ( z , t ) = u c p ( z , t ) + i u s p ( z , t ) ;   0 < z < 1 ,   t R ,
in which i is the imaginary unit. This dimensionless steady-state complex velocity has to fulfill the partial differential equation
Re ( 1 + α t + β 2 t 2 ) u p ( z , t ) t = ( 1 + γ t + δ 2 t 2 ) 2 u p ( z , t ) z 2 + ( 1 + α t + β 2 t 2 )   e i ω t Ha 2 ( 1 + α t + β 2 t 2 ) u p ( z , t ) K ( 1 + γ t + δ 2 t 2 ) u p ( z , t ) ;   0 < z < 1 ,   t R ,
with the boundary conditions
u p ( 0 , t ) = u p ( 1 , t ) = 0 ;   t > 0 .
Direct calculations show that u p ( z , t ) can be represented in the simple form
u p ( z , t ) = [ 1 cos h [ r ( z 1 / 2 ) ] cos h ( r / 2 ) ] q e i ω t r 2 ;   0 < z < 1 ,   t R ,
where the complex constants r and q are defined by the relations
r = ( Ha 2 + i ω Re ) ( 1 β ω 2 + i ω α ) + K ( 1 δ ω 2 + i ω γ ) 1 δ ω 2 + i ω γ ,   q = 1 β ω 2 + i ω α 1 δ ω 2 + i ω γ .
Consequently, the velocity fields u c p ( z , t ) and u s p ( z , t ) are given by the relations
u c p ( z , t ) =   e { [ 1 cos h [ r ( z 1 / 2 ) ] cos h ( r / 2 ) ] q e i ω t r 2 } ;   0 < z < 1 ,   t R ,
u s p ( z , t ) =   Im { [ 1 cos h [ r ( z 1 / 2 ) ] cos h ( r / 2 ) ] q e i ω t r 2 } ;   0 < z < 1 ,   t R ,
where e and Im denote the real and imaginary part, respectively, of that which follows.
Steady-state shear stresses τ c p ( z , t ) , τ s p ( z , t ) and Darcy’s resistances R c p ( z , t ) , R s p ( z , t ) corresponding to these motions, namely
τ c p ( z , t ) = e { sin h [ r ( z 1 / 2 ) ] cos h ( r / 2 ) e i ω t r } ;   0 < z < 1 ,   t R ,
τ s p ( z , t ) = Im { sin h [ r ( z 1 / 2 ) ] cos h ( r / 2 ) e i ω t r } ;   0 < z < 1 ,   t R ,
R c p ( z , t ) = K e { [ 1 cos h [ r ( z 1 / 2 ) ] cos h ( r / 2 ) ] q e i ω t r 2 } ;   0 < z < 1 ,   t R ,
R s p ( z , t ) = K Im { [ 1 cos h [ r ( z 1 / 2 ) ] cos h ( r / 2 ) ] q e i ω t r 2 } ;   0 < z < 1 ,   t R ,
have been obtained using Equations (11), (13), (26) and (27). Simple calculations show that u c p ( z , t ) , τ c p ( z , t ) , R c p ( z , t ) and u s p ( z , t ) , τ s p ( z , t ) , R s p ( z , t ) given by the above relations satisfy Equations (11)–(13) and the boundary conditions (19).
Equivalent expressions for the dimensionless steady-state velocity fields u c p ( z , t ) , u s p ( z , t ) can be obtained using the finite Fourier sine transform [23]. Multiplying Equation (17) by sin ( n π ) , integrating the obtained result between zero and one when p / z = cos ( ω   t ) and bearing in mind the identity (15) of Section 13 from the reference [23], one obtains the following ordinary differential equation for the finite Fourier sine transform u F n ( t ) of u ( z , t ) :
β Re d 3 u F n ( t ) d t 3 + b n d 2 u F n ( t ) d t 2 + c n d u F n ( t ) d t + ( λ n 2 + K e f f ) u F n ( t ) 1 + ( 1 ) n + 1 λ n [ ( 1 β ω 2 ) cos ( ω   t ) α ω sin ( ω   t ) ] = 0 ;   t > 0 ,   n = 1.2.3 ,
Here, λ n = n π with n = 1 , 2 , 3 , ,   K e f f = K + Ha 2 is the effective permeability and
b n = α Re   + β   Ha 2 + δ ( λ n 2 + K ) ,   c n = Re   + α Ha 2 + γ ( λ n 2 + K ) .
The steady-state component u F m p ( t ) of the start-up solution u F m ( t ) of the ordinary differential Equation (32) is given by the relation
u F m p ( t ) = 4 cos ( ω t ) ( λ m 2 + K e f f b m ω 2 ) ( 1 β ω 2 ) α ω 2 ( β ω 2 Re c m ) λ m [ ( λ m 2 + K e f f b m ω 2 ) 2 + ω 2 ( β ω 2 Re c m ) 2 ] 4 ω sin ( ω t ) α ( λ m 2 + K e f f b m ω 2 ) + ( 1 β ω 2 ) ( β ω 2 Re c m ) λ m [ ( λ m 2 + K e f f b m ω 2 ) 2 + ω 2 ( β ω 2 Re c m ) 2 ] ,
where m = 2 n 1 . Carrying out the inverse finite Fourier sine transform in the last equation, one finds the analytic expression of the steady state velocity u c p ( z , t ) , namely
u c p ( z , t ) = 4 cos ( ω t ) n = 1 ( λ m 2 + K e f f b m ω 2 ) ( 1 β ω 2 ) α ω 2 ( β ω 2 Re c m ) ( λ m 2 + K e f f b m ω 2 ) 2 + ω 2 ( β ω 2 Re c m ) 2 sin ( λ m z ) λ m 4 ω sin ( ω t ) n = 1 α ( λ m 2 + K e f f b m ω 2 ) + ( 1 β ω 2 ) ( β ω 2 Re c m ) ( λ m 2 + K e f f b m ω 2 ) 2 + ω 2 ( β ω 2 Re c m ) 2 sin ( λ m z ) λ m ;     0 < z < 1 ,     t R .
Analogous computations allow us to show that u s p ( z , t ) is given by the relation
u s p ( z , t ) = 4 ω cos ( ω t ) n = 1 α ( λ m 2 + K e f f b m ω 2 ) + ( 1 β ω 2 ) ( β ω 2 Re c m ) ( λ m 2 + K e f f b m ω 2 ) 2 + ω 2 ( β ω 2 Re c m ) 2   sin ( λ m z ) λ m + 4 sin ( ω t ) n = 1 ( λ m 2 + K e f f b m ω 2 ) ( 1 β ω 2 ) α ω 2 ( β ω 2 Re c m ) ( λ m 2 + K e f f b m ω 2 ) 2 + ω 2 ( β ω 2 Re c m ) 2   sin ( λ m z ) λ m ;   0 < z < 1 ,     t R .
Figure 2 clearly shows that the expressions of u c p ( z , t ) and u s p ( z , t ) given by Equations (26) and (35) and by Equations (27) and (36), respectively, are equivalent.
As was specified in Section 2, the steady-state solutions corresponding to incompressible Newtonian, Maxwell, Oldroyd-B and Burgers fluids can be obtained as limiting cases of the solutions to (26)–(31). Taking α = β = γ = δ = 0 in Equations (26)–(31), for instance, the dimensionless steady-state solutions
u N c p ( z , t ) = e { [ 1 cos h [ ( z 1 / 2 ) K e f f + i ω Re ] cos h ( K e f f + i ω Re / 2 ) ] e i ω t K e f f + i ω Re } ;   0 < z < 1 ,   t R ,
u N s p ( z , t ) = Im { [ 1 cos h [ ( z 1 / 2 ) K e f f + i ω Re ] cos h ( K e f f + i ω Re / 2 ) ] e i ω t K e f f + i ω Re } ;   0 < z < 1 ,   t R ,
τ N c p ( z , t ) = e { sin h [ ( z 1 / 2 ) K e f f + i ω Re ] cos h ( K e f f + i ω Re / 2 ) e i ω t K e f f + i ω Re } ;   0 < z < 1 ,   t R ,
τ N s p ( z , t ) = Im { sin h [ ( z 1 / 2 ) K e f f + i ω Re ] cos h ( K e f f + i ω Re / 2 ) e i ω t K e f f + i ω Re } ;   0 < z < 1 ,   t R ,
R N c p ( z , t ) = K e { [ 1 cos h [ ( z 1 / 2 ) K e f f + i ω Re ] cos h ( K e f f + i ω Re / 2 ) ] e i ω t K e f f + i ω Re } ;   0 < z < 1 ,   t R ,
R N s p ( z , t ) = K Im { [ 1 cos h [ ( z 1 / 2 ) K e f f + i ω Re ] cos h ( K e f f + i ω Re / 2 ) ] e i ω t K e f f + i ω Re } ;   0 < z < 1 ,   t R ,
corresponding to incompressible Newtonian fluids performing the same motions are obtained.
Equivalent forms for the velocity fields u N c p ( z , t ) and u N s p ( z , t ) , namely
u N c p ( z , t ) = 4 cos ( ω t ) n = 1 ( λ m 2 + K e f f ) sin ( λ m z ) λ m [ ( λ m 2 + K e f f ) 2 + ( ω Re ) 2 ]   + 4 ω Re sin ( ω t ) n = 1 sin ( λ m z ) λ m [ ( λ m 2 + K e f f ) 2 + ( ω Re ) 2 ] ;   0 < z < 1 ,   t R ,  
u N s p ( z , t ) = 4 ω Re cos ( ω t ) n = 1 sin ( λ m z ) λ m [ ( λ m 2 + K e f f ) 2 + ( ω Re ) 2 ]   + 4 sin ( ω t ) n = 1 ( λ m 2 + K e f f ) sin ( λ m z ) λ m [ ( λ m 2 + K e f f ) 2 + ( ω Re ) 2 ] ;   0 < z < 1 ,   t R ,  
have been obtained from Equations (35) and (36). To the best of our knowledge, the Newtonian solutions from Equations (35)–(44) are also new. From relations (37)–(40), (43) and (44), it results that the steady-state velocities u N c p ( z , t ) , u N s p ( z , t ) and the shear stresses τ N c p ( z , t ) , τ N s p ( z , t ) corresponding to these motions of incompressible Newtonian fluids do not depend on the Hartmann number Ha and the porosity parameter K independently but on the effective permeability K e f f = K + Ha 2 . Consequently, the investigation of these entities corresponding to MHD steady-state motions of ECIGBFs in the presence or absence of a porous medium is exactly the same task both analytically and computationally, and a two-parameter approach is superfluous.

3.1. Limiting Case ω 0 (Poiseuille Flow Due to a Constant Pressure Gradient)

Making ω 0 in Equations (37), (39) and (41), one obtains the steady components
u C p ( z ) = 1 K e f f [ 1 cos h [ ( z 1 / 2 ) K e f f ] cos h ( K e f f / 2 ) ] ;   0 < z < 1 ,   t > 0 ,
τ C p ( z ) = sin h [ ( z 1 / 2 ) K e f f ] K e f f cos h ( K e f f / 2 ) ;   0 < z < 1 ,   t > 0 ,
R C p ( z ) = K K e f f { 1 cos h [ ( z 1 / 2 ) K e f f ] cos h ( K e f f / 2 ) } ;   0 < z < 1 ,   t > 0 ,
of the dimensionless start-up solutions u C ( z , t ) ,   τ C ( z , t ) and R C ( z , t ) corresponding to the MHD Poiseuille flow of ECIGBFs through a porous medium induced by a constant pressure gradient ρ P . As expected, these solutions are identical to those corresponding to electrically conducting incompressible Newtonian fluids performing the same motion. This is possible because the governing equations corresponding to the steady motions of incompressible Newtonian and non-Newtonian fluids are identical.
Substituting the steady dimensionless velocity u C p ( z ) from Equation (45) in Equation (16), one obtains the corresponding expression of the non-dimensional volume flux, namely
Q C p = K e f f cos h ( K e f f / 2 ) 2 sin h ( K e f f / 2 ) K e f f K e f f cos h ( K e f f / 2 ) .
It results that the volume flux Q C p across a plane normal to the flow direction as well as the fluid velocity u C p ( z ) do not depend on Ha and K independently.

3.2. Reductions to Known Results from the Existing Literature

The steady-state solutions for same motions of ECIGBFs when the magnetic field or porous medium is absent can be obtained by making Ha = 0 or K = 0 , respectively, in the general results. In the absence of magnetic and porous effects, the proper solutions are obtained by taking Ha = K = 0 in the respective relations. In this case, the velocity fields u c p ( z , t ) and u s p ( z , t ) given by the relations (26) and (27) take the simplified forms
u c p ( z , t ) = 1 ω Re sin ( ω t ) + 1 ω Re e { cos h [ ( z 1 / 2 ) i ω   q Re ] cos h ( i ω   q Re / 2 ) i e i ω t } ;   0 < z < 1 ,   t R ,
u s p ( z , t ) = 1 ω Re cos ( ω t ) + 1 ω Re Im { cos h [ ( z 1 / 2 ) i ω   q Re ] cos h ( i ω   q Re / 2 ) i e i ω t } ;   0 < z < 1 ,   t R .
The dimensional forms u c p d i m ( z , t ) and u s p d i m ( z , t ) of the non-dimensional velocities u c p ( z , t ) and u s p ( z , t ) from Equations (49) and (50), namely
u c p d i m ( z , t ) = P ω sin ( ω t ) + P ω e { cos h [ m ( z d / 2 ) ] cos h ( m d / 2 ) i e i ω t } ;   0 < z < 1 ,   t R ,
u s p d i m ( z , t ) = P ω cos ( ω t ) + P ω Im { cos h [ m ( z d / 2 ) ] cos h ( m d / 2 ) i e i ω t } ;   0 < z < 1 ,   t R ,
are identical to those obtained by Fetecau et al. [7] in Equations (38) and (39). In the above relations, the imaginary constant m is given by the relation
m = ω 2   ω     [ a 3 a 1 + ( a 1 a 4 a 2 a 3 ) ω 2 ] + i [ 1 ( a 2 + a 4 ) ω 2 + ( a 1 a 3 + a 2 a 4 ω 2 ) ω 2 ( 1 a 4 ω 2 ) 2 + ( a 3 ω ) 2 .
Furthermore, taking the limit of the relations (45), (46) and (48) when the effective permeability K e f f 0 , we find the steady solutions
u C p ( z ) = ( z 2 z ) / 2 ,   τ C p ( z ) = z + 1 / 2 ,   Q C p = 1 / 12 ;   0 < z < 1 ,  
corresponding to the steady Poiseuille flow of ECIGBFs or incompressible Newtonian fluids generated by a constant pressure gradient p / z = ρ P in the absence of magnetic and porous effects. The dimensional forms of these last solutions are given by the next relations:
u C p d i m ( z ) = 1 2 μ ( z 2 d z ) p z ,   τ C p d i m ( z ) = ( z 2 d / 2 ) p z ,   Q C p d i m = d 3 12 μ   p z ;   0 < z < 1 .
The first of the last three relations, as expected, is identical to that in Equation (11) obtained by Mendiburu et al. [4]. This result, as expected, is also in accordance with the relation (6.19) of Spurk and Aksel [2] when both walls are at rest. Taking d = 2 b in the last Equation (55), we recover the result obtained by Erdogan [5] (Equation (20)), namely
Q C p d i m = 2 b 3 3 η p z .
The dimensional average velocity w C p a d i m = Q C p d i m / d and the maximum velocity w C p max d i m , as they result from Equations (55)3 and (55)1, respectively, are given by the relations
w C p a d i m = d 2 12 μ p z ,   w C p max d i m = d 2 8 μ p z = 3 2 w C p a d i m .
These results are also in accordance with those of Spurk and Axel ((Section 6.1) in ref. [2]).

4. Numerical Results and Discussing

In this work, we provide exact analytical expressions for the dimensionless steady-state velocity, non-trivial shear stress and Darcy’s resistance corresponding to the time-periodic plane Poiseuille flow of ECIGBFs through a porous medium in the presence of a uniform transversal magnetic field. These expressions have been used to find similar solutions for the Poiseuille flow of same fluids due to a constant pressure gradient. The results that have been found here can be particularized to give steady-state solutions for incompressible Burgers, Oldroyd-B, Maxwell and Newtonian fluids performing the same motions. Moreover, most of them reduce to known results from the existing literature.
To gain a deeper understanding of the acquired results, Figure 3, Figure 4, Figure 5, Figure 6, Figure 7, Figure 8, Figure 9 and Figure 10 have been prepared for increasing values of time t, distinct values for K and Ha, fixed values of the frequency ω and of the material constants α , β , γ and δ : namely, ω = π / 6 , α = 0.7 ,   β = 0.6 , γ = 0.5 and δ = 0.4 Figure 3, Figure 4, Figure 5 and Figure 6, which show the convergence of starting solutions u c ( z , t ) , u s ( z , t ) to their steady-state components u c p ( z , t ) ,   u s p ( z , t ) , give the time needed to reach a steady state. This is the time after which the diagrams of the starting solutions overlap with those of their steady-state components. From these figures, it also results that the fluid velocity diminishes for increasing values of the parameter M or K. Consequently, the fluid flows slower in the presence of a magnetic field or porous medium. A similar result has been recently obtained by Kulkarni [24] for the MHD flow of a water-based hybrid nanofluid about a permeable rotating disk.
The variations in time of the steady-state velocities u c p ( z , t ) and u s p ( z , t ) given by Equations (26) and (27) at the middle of the channel are together depicted in Figure 7 and Figure 8 at distinct values of the Hartmann number Ha and porous parameter K. In each instance, the phase difference and the oscillatory behavior of the two motions are easily visualized. Moreover, as expected, the amplitudes of oscillations are identical for motions due to cosine or sine oscillations of the pressure gradient, and they decline for increasing values of Ha or K. Consequently, the fluid velocity diminishes in the presence of a magnetic field or porous medium. It means that, as we previously mentioned, the fluid moves faster in the absence of a magnetic field or porous medium.
Profiles of the steady-state velocity u C p ( z ) and of the shear stress τ C p ( z ) are presented in Figure 9 at three values of the effective permeability K e f f . The fluid velocity, whose maximum values are touched at the middle of the channel in each case, is a decreasing function with regard to this parameter.
It results that the fluid flows slower in the presence of a magnetic field or porous medium. The shear stress is positive in the first half of the channel and negative in the other part. However, in absolute value it is a decreasing function with respect to K e f f on the whole flow domain and takes a zero value at middle of the channel.
In the last Figure 10 are depicted the profiles of Darcy’s resistance R C p ( z ) for increasing values of Ha and K, respectively. Their symmetry with respect to the median plane is clearly visualized. In all situations, as expected, the Darcy’s resistance has the minimum value in z = 0.5 , has maximum values equal to zero on the two plates and is an increasing function with respect to the Hartman number Ha and porous parameter K. It means that the fluid flows faster if the magnetic field or porous medium is absent. Of course, this outcome is in accord with those coming from Figure 7, Figure 8 and Figure 9.

5. Conclusions

Isothermal Poiseuille flows of ECIGBFs through an infinite plane channel have been investigated when magnetic and porous effects are taken into consideration. Closed-form expressions, which are new in the literature, have been established for steady-state solutions corresponding to Poiseuille flows induced by oscillatory and constant pressure gradients. Their correctness was graphically proved by means of Figure 1 or comparison with known results from the literature. More exactly, as it results from Section 3.2., known solutions from the existing literature were obtained as limiting cases of present results. In addition, the obtained results have been graphically used to bring to light the influence of a magnetic field and porous medium on the steady state and flow resistance of fluid. The main outcomes that have been obtained here are as follows:
-
The requisite time to attain the steady state diminishes with increasing values of the Hartmann number Ha or the porosity parameter K. Consequently, the steady state is achieved sooner in presence of a magnetic field or porous medium.
-
Fluid velocity is a decreasing function with respect to the effective permeability. It means that the flow velocity increases in the absence of a magnetic field or porous medium.
-
As expected, spatial variations in the fluid velocity, non-trivial shear stress and Darcy’s resistance are symmetric with respect to the median plane.

Author Contributions

Conceptualization, C.F. and S.A.; Methodology, C.F., S.A. and C.M.; Software, C.F. and S.A.; Validation, C.F., S.A. and C.M.; Writing—review and editing, C.F. and S.A. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Data are contained within the article.

Acknowledgments

The authors wish to express their sincere appreciation to reviewers and to the Academic Editor for their meticulous evaluation, valuable insights and constructive recommendations pertaining to the initial version of the manuscript.

Conflicts of Interest

The authors declare no conflicts of interest.

Nomenclature

TCauchy stress tensor
T E Extra-stress tensor
IIdentity tensor
DRate of deformation tensor
a 1 ,   a 2 ,   a 3 ,   a 4 Dimensional material constants
p Hydrostatic pressure
x , y , z Cartesian coordinates
u ( z , t ) Fluid velocity
uVelocity vector
B 0 Magnitude of the applied magnetic field
R ( z , t ) Darcy’s resistance
Q ( t ) Volume flux
kPermeability of porous medium
ReReynolds number
HaHartmann number
KPorous parameter
K e f f Effective permeability
MMagnetic parameter
τ ( z , t ) Non-trivial shear stress
φ Porosity
α , β ,   γ ,   δ Dimensionless constants
η Dynamic viscosity
ν Kinematic viscosity
ρ Fluid density
ω Frequency of oscillations
σ Electrical conductivity

References

  1. Coleman, B.D.; Markovitz, H.; Noll, W. Viscometric Flows of Non-Newtonian Fluids. Theory and Experiment; Springer: Berlin/Heidelberg, Germany, 1966. [Google Scholar]
  2. Spurk, J.H.; Aksel, N. Fluid Mechanics, 2nd ed.; Springer: Berlin/Heidelberg, Germany, 2008. [Google Scholar]
  3. Rajagopal, K.R. A note on unsteady unidirectional flows of a non-Newtonian fluid. Int. J. Non-Linear Mech. 1982, 17, 369–373. [Google Scholar] [CrossRef]
  4. Mendiburu, A.A.; Carrocci, L.R.; Carvalho, J.A. Analytical solution for transient one-dimensional Couette flow considering constant and time-dependent pressure gradients. Eng. Thermica (Therm. Eng.) 2009, 8, 92–98. [Google Scholar]
  5. Erdogan, M.E. On the unsteady unidirectional flows generated by impulsive motion of a boundary or sudden application of a pressure gradient. Int. J. Non-Linear Mech. 2002, 37, 1091–1106. [Google Scholar] [CrossRef]
  6. Hayat, T.; Siddiqui, A.M.; Asghar, S. Some simple flows of an Oldroyd-B fluid. Int. J. Eng. Sci. 2001, 39, 135–147. [Google Scholar] [CrossRef]
  7. Fetecau, C.; Hayat, T.; Fetecau, C. Steady-state solutions for some simple flows of generalized Burgers fluids. Int. J. Non-Linear Mech. 2006, 41, 880–887. [Google Scholar] [CrossRef]
  8. Danish, M.; Kumar, S.; Kumar, S. Exact analytical solutions for the Poiseuille and Couette-Poiseuille flow of third grade fluid between parallel plates. Commun. Nonlinear Sci. Numer. Simulat. 2012, 17, 1089–1097. [Google Scholar] [CrossRef]
  9. Baranovskii, E. Analytical solutions to the unsteady Poiseuille flow of a second grade fluid with slip boundary conditions. Polymers 2024, 16, 179. [Google Scholar] [CrossRef] [PubMed]
  10. Tao, L.N. Magneto hydrodynamic effects on the formation of Couette flow. J. Aerospace Sci. 1960, 27, 334–338. [Google Scholar] [CrossRef]
  11. Katagiri, M. Flow formation in Couette motion in magneto hydrodynamics. Phys. Soc. Jpn. 1962, 17, 393–396. [Google Scholar] [CrossRef]
  12. Singh, C.B.; Ram, P.C. Unsteady magneto-hydrodynamic fluid flow through a channel. J. Sci. Res. Banaras Hindu Univ. 1978, 28, 1977–1978. [Google Scholar]
  13. Kiema, D.W.; Manyonge, W.A.; Bitok, J.K.; Adenyah, R.K.; Barasa, J.S. On the steady MHD Couette flow between two infinite parallel plates in an uniform transverse magnetic field. J. Appl. Math. Bioinform. 2015, 5, 87–99. [Google Scholar]
  14. Safdar, R.; Imran, M.; Tahir, M.; Sadiq, N.; Imran, M.A. MHD flow of Burgers’ fluid under the effect of pressure gradient through a porous material pipe. Punjab Univ. J. Math. 2018, 50, 73–90. [Google Scholar]
  15. Fetecau, C.; Vieru, D. Symmetric and non-symmetric flows of Burgers’ fluids through porous media between parallel plates. Symmetry 2021, 13, 1109. [Google Scholar] [CrossRef]
  16. Krishnan, J.M.; Rajagopal, K.R. A thermodynamic framework for the constitutive modeling of asphalt concrete: Theory and application. J. Mater. Civil Eng. 2004, 16, 155–166. [Google Scholar] [CrossRef]
  17. Cramer, K.R.; Pai, S.I. Magnetofluid Dynamics for Engineers and Applied Physicists; McGraw-Hill: New York, NY, USA, 1973. [Google Scholar]
  18. Sultan, Q.; Nazar, M.; Imran, M.; Alghamdi, U.A. Flow of generalized Burgers fluid between parallel walls induced by rectified sine pulses stress. Bound. Value Probl. 2014, 2014, 152. [Google Scholar] [CrossRef]
  19. Khan, M.; Malik, R.; Anjum, A. Exact solutions of MHD second Stokes flow of generalized Burgers fluid. Appl. Math. Mech.–Engl. 2015, 36, 211–224. [Google Scholar] [CrossRef]
  20. Fetecau, C.; Akhtar, S.; Morosanu, C. Permanent solutions for MHD motions of generalized Burgers fluids adjacent to an unbounded plate subject to oscillatory shear stresses. Symmetry 2023, 15, 1683. [Google Scholar] [CrossRef]
  21. Erdogan, M.E. A note on an unsteady flow of a viscous fluid due to an oscillating plane wall. Int. J. Non-Linear Mech. 2000, 35, 1–6. [Google Scholar] [CrossRef]
  22. Prusa, V. Revisiting Stokes first and second problems for fluids with pressure-dependent viscosity. Int. J. Eng. Sci. 2010, 48, 2054–2065. [Google Scholar] [CrossRef]
  23. Sneddon, I.N. Fourier Transforms; McGraw-Hill: New York, NY, USA, 1951. [Google Scholar]
  24. Kulkarni, M. Mixed convective magnetized GO-MoS2/H2O hybrid nanofluid flow about a permeable rotating disk. Asia-Pac. J. Chem. Eng. 2023, 18, e2923. [Google Scholar] [CrossRef]
Figure 1. Flow geometry.
Figure 1. Flow geometry.
Symmetry 16 00425 g001
Figure 2. Equivalence of expressions of u c p ( z , t ) and u s p ( z , t ) given by Equations (26), (35) and (27), (36) when f ( ) ,   g ( ) are H ( ) , 4 H ( ) , respectively, for α = 0.8 ,   M = 0.7 ,   K = 0.5 ,   Re = 100   .
Figure 2. Equivalence of expressions of u c p ( z , t ) and u s p ( z , t ) given by Equations (26), (35) and (27), (36) when f ( ) ,   g ( ) are H ( ) , 4 H ( ) , respectively, for α = 0.8 ,   M = 0.7 ,   K = 0.5 ,   Re = 100   .
Symmetry 16 00425 g002
Figure 3. Convergence of starting solution u c ( z , t ) (numerical solution) to u c p ( z , t ) given by Equation (26) for α = 0.7 ,   β = 0.6 ,   γ = 0.5 ,   δ = 0.4 ,   ω = π / 6 ,   Re = 10 ,   K = 0.5 for two values of M.
Figure 3. Convergence of starting solution u c ( z , t ) (numerical solution) to u c p ( z , t ) given by Equation (26) for α = 0.7 ,   β = 0.6 ,   γ = 0.5 ,   δ = 0.4 ,   ω = π / 6 ,   Re = 10 ,   K = 0.5 for two values of M.
Symmetry 16 00425 g003
Figure 4. Convergence of starting solution u c ( z , t ) (numerical solution) to u c p ( z , t ) given by Equation (26) for α = 0.7 ,   β = 0.6 ,   γ = 0.5 ,   δ = 0.4 ,   ω = π / 6 ,   Re = 10 ,   M = 0.8 for two values of K.
Figure 4. Convergence of starting solution u c ( z , t ) (numerical solution) to u c p ( z , t ) given by Equation (26) for α = 0.7 ,   β = 0.6 ,   γ = 0.5 ,   δ = 0.4 ,   ω = π / 6 ,   Re = 10 ,   M = 0.8 for two values of K.
Symmetry 16 00425 g004
Figure 5. Convergence of starting solution u s ( z , t ) (numerical solution) to u s p ( z , t ) given by Equation (27) for α = 0.7 ,   β = 0.6 ,   γ = 0.5 ,   δ = 0.4 ,   ω = π / 6 ,   Re = 10 ,   K = 0.5 for two values of M.
Figure 5. Convergence of starting solution u s ( z , t ) (numerical solution) to u s p ( z , t ) given by Equation (27) for α = 0.7 ,   β = 0.6 ,   γ = 0.5 ,   δ = 0.4 ,   ω = π / 6 ,   Re = 10 ,   K = 0.5 for two values of M.
Symmetry 16 00425 g005
Figure 6. Convergence of starting solution u s ( z , t ) (numerical solution) to u s p ( z , t ) given by Equation (27) for α = 0.7 ,   β = 0.6 ,   γ = 0.5 ,   δ = 0.4 ,   ω = π / 6 ,   Re = 10 ,   M = 0.8 for two values of K.
Figure 6. Convergence of starting solution u s ( z , t ) (numerical solution) to u s p ( z , t ) given by Equation (27) for α = 0.7 ,   β = 0.6 ,   γ = 0.5 ,   δ = 0.4 ,   ω = π / 6 ,   Re = 10 ,   M = 0.8 for two values of K.
Symmetry 16 00425 g006
Figure 7. Variations in time of velocities u c p ( z , t ) and u s p ( z , t ) given by Equations (26) and (27) at the middle of channel when α = 0.7 ,   β = 0.6 ,   γ = 0.5 ,   δ = 0.4 ,   ω = π / 6 ,   Re = 100 ,   K = 0.5   for three values of Ha.
Figure 7. Variations in time of velocities u c p ( z , t ) and u s p ( z , t ) given by Equations (26) and (27) at the middle of channel when α = 0.7 ,   β = 0.6 ,   γ = 0.5 ,   δ = 0.4 ,   ω = π / 6 ,   Re = 100 ,   K = 0.5   for three values of Ha.
Symmetry 16 00425 g007
Figure 8. Variations in time of velocities u c p ( z , t ) and u s p ( z , t ) given by Equation (26) and (27) at the middle of channel when α = 0.7 ,   β = 0.6 ,   γ = 0.5 ,   δ = 0.4 ,   ω = π / 6 ,   Re = 100 ,   Ha = 1 for three values K.
Figure 8. Variations in time of velocities u c p ( z , t ) and u s p ( z , t ) given by Equation (26) and (27) at the middle of channel when α = 0.7 ,   β = 0.6 ,   γ = 0.5 ,   δ = 0.4 ,   ω = π / 6 ,   Re = 100 ,   Ha = 1 for three values K.
Symmetry 16 00425 g008
Figure 9. Profiles of u C p ( z ) and τ C p ( z ) given by Equations (45) and (46) at three values of K e f f .
Figure 9. Profiles of u C p ( z ) and τ C p ( z ) given by Equations (45) and (46) at three values of K e f f .
Symmetry 16 00425 g009
Figure 10. Profiles of the Darcy’s resistance R C p ( z ) given by Equation (47) for K = 0.5 and three values of Ha and for Ha = 1 and three values of K.
Figure 10. Profiles of the Darcy’s resistance R C p ( z ) given by Equation (47) for K = 0.5 and three values of Ha and for Ha = 1 and three values of K.
Symmetry 16 00425 g010
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Fetecau, C.; Akhtar, S.; Moroşanu, C. MHD Time-Periodic Plane Poiseuille Flow of Generalized Burgers Fluids through a Porous Medium. Symmetry 2024, 16, 425. https://doi.org/10.3390/sym16040425

AMA Style

Fetecau C, Akhtar S, Moroşanu C. MHD Time-Periodic Plane Poiseuille Flow of Generalized Burgers Fluids through a Porous Medium. Symmetry. 2024; 16(4):425. https://doi.org/10.3390/sym16040425

Chicago/Turabian Style

Fetecau, Constantin, Shehraz Akhtar, and Costică Moroşanu. 2024. "MHD Time-Periodic Plane Poiseuille Flow of Generalized Burgers Fluids through a Porous Medium" Symmetry 16, no. 4: 425. https://doi.org/10.3390/sym16040425

APA Style

Fetecau, C., Akhtar, S., & Moroşanu, C. (2024). MHD Time-Periodic Plane Poiseuille Flow of Generalized Burgers Fluids through a Porous Medium. Symmetry, 16(4), 425. https://doi.org/10.3390/sym16040425

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