Next Article in Journal
Compiled Constructions towards Post-Quantum Group Key Exchange: A Design from Kyber
Previous Article in Journal
Factors Affecting the Behavioral Intention to Adopt Mobile Payment: An Empirical Study in Taiwan
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Hydrodynamical Study of Creeping Maxwell Fluid Flow through a Porous Slit with Uniform Reabsorption and Wall Slip

1
Department of Mathematics, Faculty of Science, Jiangsu University, Zhenjiang 212013, China
2
Department of Mathematics, COMSATS University Islamabad, Sahiwal 57000, Pakistan
3
Department of Mathematics, Pennsylvania State University, York Campus, 1031 Edgecomb Avenue, York, PA 17403, USA
4
Department of Mathematics and Statistics, International Islamic University, Islamabad 45320, Pakistan
*
Authors to whom correspondence should be addressed.
Mathematics 2020, 8(10), 1852; https://doi.org/10.3390/math8101852
Submission received: 19 September 2020 / Revised: 7 October 2020 / Accepted: 14 October 2020 / Published: 21 October 2020

Abstract

:
The present theoretical study investigates the influence of velocity slip characteristics on the plane steady two-dimensional incompressible creeping Maxwell fluid flow passing through a porous slit with uniform reabsorption. This two-dimensional flow phenomenon is governed by the mathematical model having nonlinear partial differential equations together with non-homogeneous boundary conditions. An analytical technique, namely the recursive approach, is used successfully to find the solutions of the problem. The explicit expressions for stream function, velocity components, pressure distribution, wall shear stress and normal stress difference have been derived. The axial flow rate, leakage flux and fractional reabsorption are also found out. The points of maximum velocity are identified. Non-dimensionalization is carried out and graphs are portrayed at different positions of the channel to show the impact of pertinent parameters: slip parameter, Maxwell fluid parameter and absorption parameter, on flow variables and found that the fluid velocity is affected significantly due to these parameters. This study provides a mathematical basis to understand the physical phenomenon for fluid flows through permeable boundaries which exists in different problems like gaseous diffusion, filtration and biological mechanisms.

1. Introduction

Flow-through permeable boundaries have enormous importance from many decades due to their tremendous applications in bio-sciences and engineering, such as processes like membrane filtration, desalination processes using reverse osmosis, transpiration cooling, blood flow, renal proximal tubule flow within a kidney and filtration of blood in hemodialysis of an artificial kidney are the key examples related to flows in permeable boundaries [1,2,3,4,5,6]. Berman [7] presented the complete description of the Newtonian fluid flow through a channel with permeable walls and the perturbation method is employed to elaborate the wall permeability impact on velocity and pressure distribution. Sellars [8] and Yuan [9] extended this work to discuss the channel flow problem for high suction Reynolds number and for small and large seepage rates, respectively. Afterward, many researchers [10,11,12] further investigated the flow-through channel with permeable walls.
Due to the practical applications of creeping flow, many researchers studied the slow flow of Newtonian as well as non-Newtonian fluids. Macey [4] considered blood as an incompressible Newtonian fluid through the renal tubule and assumed creeping flow. Marshall [6] investigated the slow flow through a porous channel and discussed the obtained approximate solutions for ultrafiltration rate in flat plate hemodialyzer. Siddiqui and Haroon [13,14,15,16] considered different absorption rates on walls of the permeable slit to analyze the creeping flow of a Newtonian fluid. Rajagopal [17] considered the second order fluid to discuss the creeping flow. Ullah et al. [18,19] discussed the slow flow by considering the slightly viscoelastic fluid with uniform reabsorption using a recursive approach. Khashan et al. [20,21,22] also considered the creeping flow of micro-polar and Jeffrey fluids through a porous channel. Motivated with this, we have considered the Maxwell fluid to investigate the slow flow through permeable slit having wall slip with uniform reabsorption and present the solutions using the recursive approach proposed by Langlois [23,24].
For many decades, authors have been interested with slip condition due to its practical implications in the polymer industry (polymer melts) [25,26], which shows a macroscopic wall slip and also it extends from technological applications to medical applications as it used in polishing artificial heart valves. Beavers and Joseph [27] experimentally confirmed that the slip condition exists at the porous boundary. Saffman [28] gave the mathematical justification to the study presented by Beavers [27] and later, Beavers and Kohler [29,30] revealed that the already presented analytical results and their experimental outcomes are in good agreement. Mikelic and Jager [31] presented the statistical and mathematical analysis to discuss the slip above the porous surface. Rao et al. [32] discussed the slippage effects on viscous fluid flow in the channel. Elshahed [33] considered the velocity slip on membrane surface for blood flow in microcirculation. Singh et al. [34] showed the importance of velocity slip on ultra-filtration. Makinde and Osalusi [35] has explained the slip effects on the magnetohydrodynamic Newtonian fluid flow in a channel with porous walls. Eldesoky et al. [36] investigated the magnetohydrodynamics and slip effects on pulsatile blood flow through a porous medium having permeable walls and Elshehawy et al. [37] also discussed the velocity slip impact on peristaltic flow of Maxwellian fluid. Ellahi et al. [38] considered the channel flow of non-Newtonian fluid along with slip condition. Hron et al. [39] also considered the Naviers slip conditions at the boundary and discussed its impact on non-Newtonian fluids. Hayat et al. [40] discussed the Oldroyd 6-constant fluid flowing between parallel plates under the influenced of slip condition.
The significant importance of non-Newtonian fluid arises in plastic manufacturing, lubrication theory, paint industry and in food processing. Most of the biologically important fluids are considered as non-Newtonian due to their higher molecular weight. There are three categories of non-Newtonian fluids, namely: differential, rate and integral type fluids. The models of rate type fluids [41] has its importance and exhibits the stress relaxation aspects of polymer solutions. The expressions for stresses in terms of velocity components can be easily found out due to the existence of their explicit expressions in constitutive equations of differential type fluids, whether in the case of rate type fluids, there exits an implicit relationship in constitutive equation and its difficult to express stresses in terms of velocity components. This fact is the key reason for the lack of literature on the flow of rate type fluids through the permeable slit in the two dimensional case. To contribute with this, we consider the slow flow of Maxwell fluid in the two dimensional case through a permeable slit. An effort has been put by different researchers to investigate the flow of an upper convected Maxwell fluid. Choi et al. [42] investigated the combined effects of inertia, suction and viscoelasticity in a porous channel. Sadeghy et al. [43] studied Sakiadis flow and showed that wall skin friction decreases when increasing the value of Deborah’s number. Abbas et al. [44] also discussed the boundary layer magnetohydrodynamic Maxwell fluid flow in a permeable channel.
We consider the creeping flow of Maxwell fluid through the permeable slit in accordance with slip condition at porous walls and uniform reabsorption. The modeled nonlinear PDEs in accordance with non-homogeneous boundary conditions are solved with the recursive approach proposed by Langlois [23,24], who successfully implemented this approach for one-dimensional slow non-Newtonian fluid flows. Ullah et al. [18,19] used this technique for the first time to study two-dimensional creeping flow of a slightly viscoelastic fluid through a permeable slit. Recently, Bhatti et al. [45] used this approach to investigate the slow flow of a second order fluid through a uniformly permeable circular tube and now we implement this approach to get solutions of creeping flow of a Maxwell fluid and discuss combined effects of absorption parameter, viscoelasticity and slip parameter on flow variables involved at different positions of the channel and also obtain the explicit expressions for axial and radial velocity components, total pressure difference, mean pressure drop, normal stresses difference, wall shear stress, leakage flux and fractional reabsorption.

2. Formulation of the Problem

A two-dimensional isothermal, steady, incompressible creeping Maxwell fluid flow through a rectangular slit of length L, breadth W and width 2 H is considered, with the x-axis taken along the centerline of the rectangular slit and y-axis taken normal to it. The walls of the slit are porous with a constant velocity of ϵ V 0 along the slit length, where ϵ is a small dimensionless parameter. A slip condition is also considered at porous walls of the slit and we assumed a constant flow rate ϵ Q 0 , at x = 0 , (see Figure 1).
In the proposed problem, the two-dimensional flow through the slit has the following velocity profile,
V = ( u ( x , y ) , v ( x , y ) )
The conservation of mass and momentum equations for a creeping two-dimensional Maxwell fluid takes the form
u x + v y = 0 ,
p x = τ x x x + τ x y y ,
p y = τ x y x + τ y y y ,
τ x x + λ ( u τ x x x + v τ x x y 2 τ x x u x 2 τ x y u y ) = 2 μ u x ,
τ x y + λ ( u τ x y x + v τ x y y τ x x v x τ y y u y ) = μ ( u y + v x ) ,
τ y y + λ ( u τ y y x + v τ y y y 2 τ y y v y 2 τ x y v x ) = 2 μ v y ,
where, (u, v) are axial and radial components of velocity, the constant λ is the fluid relaxation time and μ is the viscosity of fluid. p and τ represent pressure and stresses of the Maxwell fluid, respectively.
The consequent boundary conditions are imposed as below by considering the upper half of the slit as due to symmetry along the centerline of the slit.
u = k γ u y , v = ϵ V 0 , at y = H ,
u y = 0 , v = 0 at y = 0
ϵ Q 0 = W H H ( u ( 0 , y ) ) d y .
where k represents the specific permeability of a porous medium and γ is a dimensionless constant. Here, the boundary condition (8) is the well-known slip condition given by Beavers and Joseph [27].

3. Solution of the Problem

The exact solution of the system of highly non-linear and coupled partial differential Equations (2)–(7) having six unknowns, namely, u ( x , y ) , v ( x , y ) , p ( x , y ) , τ x y ( x , y ) , τ x x ( x , y ) , τ y y ( x , y ) , together with non-homogeneous boundary conditions (8)–(10) is quite impossible and even analytical and numerical solutions of such equations are still a challenge for researchers. Here, we will use a recursive approach suggested by Langlois [23,24], which linearizes the above equations; so for this, we consider the flow field as a perturbation of a state of rest and the following is
V ( x , y ) = j = 1 j = ϵ j V ( j ) ( x , y ) ,
p ( x , y ) = C o n s t a n t + j = 1 j = ϵ j p ( j ) ( x , y ) ,
τ ( x , y ) = j = 1 j = ϵ j τ ( j ) ( x , y ) .
where ϵ is a small dimensionless number.
The above supposition set up a dynamical system of linear equations and boundary conditions for each [ u j , v j , p j , τ j ] , j = 1 , 2 , 3 and then [ u , v , p , τ ] can be found by Equations (11)–(13), which will represent the solution for a creeping Maxwell fluid. The equations obtained for the first order [ u 1 , v 1 , p 1 , τ 1 ] are identical to the Newtonian fluid flow. The expressions obtained at the second order [ u 2 , v 2 , p 2 , τ 2 ] are similar with the exception that there exists non homogeneous terms depending on [ u 1 , v 1 , p 1 , τ 1 ] . Likewise, [ u 3 , v 3 , p 3 , τ 3 ] are also similar but there also exits involvement of lower order solutions [ u 1 , v 1 , p 1 , τ 1 ] , [ u 2 , v 2 , p 2 , τ 2 ] and this process continues recursively. It is also observed that a linear dynamical system of equations are formed at each stage of analysis and explicitly used results from previous stages. There is not too much practical or conceptual difficulty arising in extending the theory for higher orders.
In result of using Equations (11)–(13) into Equations (2)–(10) and comparing the coefficients of like powers of ϵ , we get three sets of systems of PDEs at O ( ϵ ) , O ( ϵ 2 ) and O ( ϵ 3 ) .

3.1. System of Equations for the 1st Order

The first order system of equations is obtained by comparing coefficients of ϵ , as
u ( 1 ) x + v ( 1 ) y = 0 ,
p ( 1 ) x = τ x x ( 1 ) x + τ x y ( 1 ) y ,
p ( 1 ) y = τ x y ( 1 ) x + τ y y ( 1 ) y ,
τ x x ( 1 ) = 2 μ u ( 1 ) x , τ x y ( 1 ) = μ u ( 1 ) y + v ( 1 ) x , τ y y ( 1 ) = 2 μ v ( 1 ) y .
u ( 1 ) = k γ u ( 1 ) y , v ( 1 ) = V 0 , at y = H ,
u ( 1 ) y = 0 , v ( 1 ) = 0 , at y = 0 ,
0 H ( u ( 1 ) ( 0 , y ) ) d y = Q 0 2 W , at x = 0 .
Here, we consider ψ ( 1 ) ( x , y ) as
u ( 1 ) = ψ ( 1 ) ( x , y ) y , v ( 1 ) = ψ ( 1 ) ( x , y ) x .
We found that in view of Equation (21), Equation (14) is satisfied identically and after elimination of the pressure term, Equations (15), (16) and (18)–(20) in terms of ψ ( x , y ) become,
4 ψ ( 1 ) = 0 .
ψ ( 1 ) y = k γ 2 ψ ( 1 ) y 2 , ψ ( 1 ) x = v 0 , at y = H ,
2 ψ ( 1 ) y 2 = 0 , ψ ( 1 ) x = 0 , at y = 0 ,
ψ ( 1 ) ( 0 , H ) = Q 0 2 W , ψ ( 1 ) ( 0 , 0 ) = 0 .
The boundary conditions (23)–(25) suggest to consider ψ ( 1 ) ( x , y ) as
ψ ( 1 ) ( x , y ) = V 0 x M ( 1 ) ( y ) + N ( 1 ) ( y ) .
where M ( 1 ) ( y ) and N ( 1 ) ( y ) are unknowns, which can be obtained by using Equation (26) in Equations (22)–(25) and ultimately we have
M ( 1 ) ( y ) = 1 2 b 1 y 3 H 3 3 b 2 y H ,
N ( 1 ) ( y ) = b 1 Q 0 4 W 3 b 2 y H y 3 H 3 .
where ϕ = k H γ is the slip parameter and ( b 1 = 1 1 + 3 ϕ ) and b 2 = 1 + 2 ϕ are functions of ϕ .
Equation (26) in results of utilizing Equations (21), (27) and (28), we obtain
ψ ( 1 ) = b 1 2 ( Q 0 2 W x V 0 ) ( 3 b 2 y H y 3 H 3 ) ,
u ( 1 ) = 3 b 1 2 H ( Q 0 2 W x V 0 ) ( b 2 y 2 H 2 ) ,
v ( 1 ) = b 1 V 0 2 ( 3 b 2 y H y 3 H 3 ) ,
Here, we observed that the obtained first-order solution is strongly dependent on ϕ and no contribution of Maxwell fluid parameter δ is found. It is observed that the expressions for velocity components at first order with ( ϕ = 0 ) are similar to the results presented in [15].

3.2. System of Equations for the 2nd Order

On equating terms involving ϵ 2 , we have the following system of equations:
u ( 2 ) x + v ( 2 ) y = 0 ,
p ( 2 ) x = τ x x ( 2 ) x + τ x y ( 2 ) y ,
p ( 2 ) y = τ x y ( 2 ) x + τ y y ( 2 ) y ,
τ x x ( 2 ) + λ ( u ( 1 ) τ x x ( 1 ) x + v ( 1 ) τ x x ( 1 ) y 2 τ x x ( 1 ) u ( 1 ) x 2 τ x y ( 1 ) u ( 1 ) y ) = 2 μ u ( 2 ) x ,
τ x y ( 2 ) + λ ( u ( 1 ) τ x y ( 1 ) x + v ( 1 ) τ x y ( 1 ) y τ x x ( 1 ) v ( 1 ) x τ y y ( 1 ) u ( 1 ) y ) = μ ( v ( 2 ) x + u ( 2 ) y ) ,
τ y y ( 2 ) + λ ( u ( 1 ) τ y y ( 1 ) x + v ( 1 ) τ y y ( 1 ) y 2 τ y y ( 1 ) v ( 1 ) y 2 τ x y ( 1 ) v ( 1 ) x ) = 2 μ v ( 2 ) y ,
u ( 2 ) = k γ u ( 2 ) y , v ( 2 ) = 0 , at y = H ,
u ( 2 ) y = 0 , v ( 2 ) = 0 , at y = 0 ,
0 H ( u ( 2 ) ( 0 , y ) ) d y = 0 , at x = 0 .
We consider ψ ( 2 ) ( x , y ) as
u ( 2 ) = ψ ( 2 ) ( x , y ) y , v ( 2 ) = ψ ( 2 ) ( x , y ) x .
Equation (32) is satisfied identically. Using the first order solution (29)–(31) in Equations (35)–(37) and then eliminating pressure from Equations (33) and (34), the resulting equation in terms of stream function becomes
4 ψ ( 2 ) = 0 ,
The boundary conditions at this order take the following form
ψ ( 2 ) y = k γ 2 ψ ( 2 ) y 2 , ψ ( 2 ) x = 0 , at y = H ,
2 ψ ( 2 ) y 2 = 0 , ψ ( 2 ) x = 0 , at y = 0 ,
ψ ( 2 ) ( 0 , H ) = 0 , ψ ( 2 ) ( 0 , 0 ) = 0 .
The solution of the system of Equations (42)–(45) becomes zero due to the homogeneous boundary conditions for any supposed form of ψ ( 2 ) ( x , y ) . Therefore,
ψ ( 2 ) = 0 ,
u ( 2 ) = 0 ,
v ( 2 ) = 0 ,
Thus, at this order, we do not see any contribution of Newtonian or viscoelastic effects towards the velocity field.

3.3. System of Equations for the 3rd Order

Equating terms involving ϵ 3 and the obtained system of equations along with boundary conditions are
u ( 3 ) x + v ( 3 ) y = 0 ,
p ( 3 ) x = τ x x ( 3 ) x + τ x y ( 3 ) y ,
p ( 3 ) y = τ x y ( 3 ) x + τ y y ( 3 ) y ,
τ x x ( 3 ) + λ ( u ( 1 ) τ x x ( 2 ) x + v ( 1 ) τ x x ( 2 ) y 2 τ x x ( 2 ) u ( 1 ) x 2 τ x y ( 2 ) u ( 1 ) y + u ( 2 ) τ x x ( 1 ) x + v ( 2 ) τ x x ( 1 ) y 2 τ x x ( 1 ) u ( 2 ) x 2 τ x y ( 1 ) u ( 2 ) y ) = 2 μ u ( 3 ) x ,
τ x y ( 3 ) + λ ( u ( 1 ) τ x y ( 2 ) x + v ( 1 ) τ x y ( 2 ) y τ x x ( 2 ) v ( 1 ) x τ y y ( 2 ) u ( 1 ) y + u ( 2 ) τ x y ( 1 ) x + v ( 2 ) τ x y ( 1 ) y τ x x ( 1 ) v ( 2 ) x τ y y ( 1 ) u ( 2 ) y ) = μ ( u ( 3 ) y + v ( 3 ) x ) ,
τ y y ( 3 ) + λ ( u ( 1 ) τ y y ( 2 ) x + v ( 1 ) τ y y ( 2 ) y 2 τ y y ( 2 ) v ( 1 ) y 2 τ x y ( 2 ) v ( 1 ) x + τ y y ( 1 ) x ( u ( 2 ) ) + τ y y ( 1 ) y ( v ( 2 ) ) 2 τ y y ( 1 ) v ( 2 ) y 2 τ x y ( 1 ) v ( 2 ) x ) = 2 μ v ( 3 ) y ,
u ( 3 ) = k γ u ( 3 ) y , v ( 3 ) = 0 , at y = H ,
u ( 3 ) y = 0 , v ( 3 ) = 0 , at y = 0 ,
0 H ( u ( 3 ) ( 0 , y ) ) d y = 0 , at x = 0 .
The ψ ( 3 ) ( x , y ) in the third order is define as
u ( 3 ) = ψ ( 3 ) ( x , y ) y , v ( 3 ) = ψ ( 3 ) ( x , y ) x .
Equation (49) is satisfied identically. On substituting the first order solution (29)–(31) and the second order solution (46)–(48) along with Equation (58) into the third order system (53)–(55) and then eliminating the pressure by cross differentiation of Equations (50) and (51), we arrive at Equation (59).
4 ψ ( 3 ) = 72 b 1 3 λ 2 V 0 2 H 2 Q 0 2 W x V 0 y 3 H 7 ,
Similarly,
4 ψ ( 3 ) = 72 b 1 3 δ 2 Q 0 2 W x V 0 y 3 H 7 ,
where δ = λ V 0 H is the Maxwell fluid parameter. The boundary conditions at this order are
ψ ( 3 ) y = k γ 2 ψ ( 3 ) y 2 , ψ ( 3 ) x = 0 , at y = H ,
2 ψ ( 3 ) y 2 = 0 , ψ ( 3 ) x = 0 , at y = 0 ,
ψ ( 3 ) ( 0 , H ) = 0 , ψ ( 3 ) ( 0 , 0 ) = 0 .
For the solution of Equations (60)–(63), we supposed the ψ ( 3 ) as
ψ ( 3 ) = 72 b 1 3 δ 2 ( Q 0 2 W x V 0 ) M ( 3 ) ( y ) + N ( 3 ) ( y ) ,
where M ( 3 ) ( y ) and N ( 3 ) ( y ) are unknowns and obtained after considerable calculations in utilizing Equation (64) in Equations (60)–(63) as
M ( 3 ) i v ( y ) = y 3 H 7 ,
M ( 3 ) ( 0 ) = M ( 3 ) ( H ) = M ( 3 ) ( 0 ) = 0 , M ( 3 ) ( H ) = k γ M ( 3 ) ( H ) ,
N ( 3 ) i v ( y ) = 0 .
N ( 3 ) ( 0 ) = N ( 3 ) ( H ) = N ( 3 ) ( 0 ) = 0 , N ( 3 ) ( H ) = k γ N ( 3 ) ( H ) .
Solving Equations (65)–(68), we obtain
M ( 3 ) ( y ) = 1 840 y H 7 3 b 1 b 3 840 y H 3 + 2 b 1 b 4 840 y H ,
N ( 3 ) ( y ) = 0 .
where b 3 = 1 + 7 ϕ , b 4 = 1 + 9 ϕ .
Using Equations (69) and (70) in Equation (64), ultimately, we have
ψ ( 3 ) = 3 35 b 1 3 δ 2 ( Q 0 2 W x V 0 ) y H 7 3 b 1 b 3 y H 3 + 2 b 1 b 4 y H ,
u ( 3 ) = 3 35 H b 1 3 δ 2 ( Q 0 2 W x V 0 ) 7 y H 6 9 b 1 b 3 y H 2 + 2 b 1 b 4 ,
v ( 3 ) = 3 35 b 1 3 δ 2 V 0 y H 7 3 b 1 b 3 y H 3 + 2 b 1 b 4 y H ,
The contribution of δ is observed in the solution obtained in the third order for velocity field, which is the key parameter of investigation of the present study.
Summing up all the obtained solutions, we have
ψ = 1 35 ( Q o 2 W x V 0 ) ( 3 b 1 3 δ 2 y H 7 + ( 35 b 1 2 + 9 b 1 4 b 3 δ 2 ) ( y H ) 3 + ( 105 b 1 b 2 2 6 b 1 4 b 4 δ 2 ) ( y H ) ) ,
u = 1 35 H ( Q o 2 W x V 0 ) ( 21 b 1 3 δ 2 y H 6 + 3 ( 35 b 1 2 + 9 b 1 4 b 3 δ 2 ) ( y H ) 2 + ( 105 b 1 b 2 2 6 b 1 4 b 4 δ 2 ) ) ,
v = V 0 35 ( 3 b 1 3 δ 2 y H 7 + ( 35 b 1 2 + 9 b 1 4 b 3 δ 2 ) ( y H ) 3 + ( 105 b 1 b 2 2 6 b 1 4 b 4 δ 2 ) ( y H ) ) ,
Equations (75)–(77) give the expression of velocity profile strongly influenced by slip parameter ϕ and Maxwell parameter δ for creeping Maxwell fluid flow through a permeable slit.
Here, we present the expressions for remaining flow properties.
(a)
The axial velocity is maximum along the center line of the slit as
u m a x = 1 35 H Q 0 2 W x V 0 105 b 1 b 2 2 6 b 1 4 b 4 δ 2 .
(b)
The maximum radial velocity occurs at the slit walls, i.e.,
v m a x = V 0 b 1 70 ( 35 + 105 b 2 + 6 b 1 2 ( 1 + 3 b 1 b 3 2 b 1 b 4 ) δ 2 ) .
(c)
The axial flow rate is obtained using the summarized solution as
Q ( x ) = 2 W 0 H ( u ( x , y ) ) d y , = b 1 ( Q 0 2 W V 0 x ) 70 35 + 105 b 2 + 6 b 1 2 ( 1 + 3 b 1 b 3 2 b 1 b 4 ) δ 2 .
Here, the dependence of Q ( x ) on δ is only due to the presence of ϕ , if ϕ = 0 then
Q ( x ) = Q 0 2 W V 0 x
(d)
The leakage flux q ( x ) is obtained as
q ( x ) = d Q d x , = b 1 W V 0 35 ( 35 + 105 b 2 + 6 b 1 2 ( 1 + 3 b 1 b 3 2 b 1 b 4 ) δ 2 ) ) .
(e)
The fractional reabsorption F a is obtained as
F a = 1 Q ( 0 ) Q ( 0 ) Q ( L ) , = 2 L W V 0 Q 0 .
The contribution of δ in leakage flux q ( x ) is only due to ϕ and F a is only influenced by absorption parameter.
(f)
Pressure Distribution
Here, we will find out the pressure for each order by utilizing Equations (30) and (31) into Equations (15)–(17) to get the pressure for the first order
p ( 1 ) x = 3 b 1 μ Q 0 2 W x V 0 2 H 3 W ,
p ( 1 ) y = 3 b 1 μ y V 0 H 3 .
Equation (82) on integration with respect to x yields
p ( 1 ) = 3 b 1 μ x Q 0 W x 2 V 0 2 H 3 W + B ( y ) ,
where the unknown function B ( y ) is to be determined as follows:
Differentiating Equation (84) with respect to y, and on comparing with Equation (83), gives
B ( y ) = 3 b 1 μ y V 0 H 3 ,
Then, Equation (84) becomes
p ( 1 ) p 0 ( 1 ) = 3 b 1 μ x Q 0 W x 2 V 0 2 H 3 W 3 b 1 μ y 2 V 0 2 H 3 ,
where p 0 ( 1 ) = p ( 0 , 0 ) is the pressure at ( x , y ) = ( 0 , 0 ) , i.e., the entrance region of the slit.
The mean pressure is obtained as
p ¯ ( 1 ) ( x ) = 1 H 0 H ( p ( 1 ) p 0 ( 1 ) ) d y , = b 1 μ 2 H 3 x W x V 0 Q 0 H 2 W V 0 ,
The pressure drop over slit length L is obtained as
Δ p ¯ 1 ( L ) = p ¯ 1 ( 0 ) p ¯ 1 ( L ) , = 3 b 1 μ L Q 0 L W V 0 2 H 3 W .
Similarly, we can obtain
p ( 2 ) p 0 ( 2 ) = 3 b 1 2 δ μ 2 H 3 b 2 x Q 0 + W x v 0 H 2 W 3 b 2 v 0 y H 2 + 2 v 0 y H 4 ,
p ¯ ( 2 ) ( x ) = 3 b 1 2 δ μ 10 H 15 b 2 x Q 0 + W x v 0 H 2 W + ( 2 5 b 2 ) v 0 ,
Δ p ¯ 2 ( L ) = 9 b 1 2 b 2 L δ μ Q 0 L W v 0 2 H 3 W
p ( 3 ) p 0 ( 3 ) = 9 b 1 3 δ 2 μ 70 H ( 28 V 0 ( y H ) 6 + ( 3 35 b 2 2 + 2 b 1 b 3 V 0 ) y H 2 3 x Q 0 W x V 0 H 2 W ( 35 b 2 2 2 b 1 b 3 ) ) ,
p ¯ ( 3 ) ( x ) = 9 b 1 3 δ 2 μ 70 H ( 3 x Q 0 W x V 0 H 2 W ( 35 b 2 2 + 2 b 1 b 3 ) + ( 4 35 b 2 2 + 2 b 1 b 3 ) V 0 )
Δ p ¯ 3 ( L ) = 27 b 1 3 L δ 2 μ Q 0 L W v 0 70 H 3 W 35 b 2 2 2 b 1 b 3
Here, we note that the pressure field has considerable contribution from viscoelastic effect for the second order even though the velocity field had zero contribution at that order.
The summarized forms of total pressure difference, mean pressure and pressure drop are given as
p ( x , y ) p ( 0 , 0 ) = 3 b 1 μ v 0 70 H ( 84 b 1 2 δ 2 y H 6 + 70 b 1 δ y H 4 + ( 35 105 b 1 b 2 δ + 315 b 1 2 b 2 2 + 18 b 1 3 b 3 δ 2 ) ) y H 2 ) + 3 b 1 x μ Q 0 W x v 0 70 H 3 W ( 35 105 b 1 b 2 δ + 315 b 1 2 b 2 2 + 18 b 1 3 b 3 δ 2 )
p ¯ ( x ) = b 1 μ ( 3 x Q 0 W x v 0 70 H 3 W ( 35 105 b 1 b 2 δ + ( 315 b 1 2 b 2 2 + 18 b 1 3 b 3 ) δ 2 ) + v 0 70 H ( 35 + ( 42 b 1 105 b 1 b 2 ) δ + 9 b 1 2 ( 4 35 b 2 2 + 2 b 1 b 3 ) δ 2 ) )
Δ p ¯ ( L ) = 3 b 1 L μ Q 0 + L W v 0 70 H 3 W ( 35 105 b 1 b 2 δ + ( 315 b 1 2 b 2 2 + 18 b 1 3 b 3 ) δ 2 )
We noticed that p ( x , y ) , Δ p ¯ ( L ) and p ¯ ( x ) are varying with slip coefficient ϕ , and Maxwell parameter δ . In the limiting case, when ϕ = 0 , δ = 0 , we recovered the exact expressions of p ( x , y ) , Δ p ¯ ( L ) and p ¯ ( x ) of a Newtonian fluid [15].
(g)
The wall shear stress is obtained as
τ w = τ x y | y = H = 3 b 1 μ Q 0 2 W x V 0 70 H 2 W ( 35 + 35 b 1 ( 4 + 3 b 2 ) δ + 9 b 1 2 ( 28 + 35 b 2 2 2 b 1 b 3 ) δ 2 ) .
(h)
The expressions for normal stress differences are given as
τ n = τ x x τ y y = 3 b 1 μ v 0 70 H ( 1596 b 1 2 δ 2 ( y H ) 6 35 b 1 ( 4 δ + 144 b 1 b 2 δ 2 + 9 b 1 δ 2 Q 0 2 W x v 0 2 H 2 W 2 v 0 2 ) ( y H ) 4 + ( 140 420 b 1 b 2 δ 72 b 1 2 ( 70 b 2 2 + b 1 b 3 ) δ 2 + 105 b 1 δ Q 0 2 W x v 0 2 H 2 W 2 v 0 2 ) ( y H ) 2 140 b 2 4 b 1 2 ( 315 b 2 3 4 b 1 b 4 ) δ 2 )
τ n | y = H = 3 b 1 μ v 0 70 H ( 140 ( 1 + b 2 ) + 140 b 1 ( 1 3 b 2 ) δ + 4 b 1 2 ( 399 315 ( 2 + b 2 ) 2 b 2 18 b 1 b 3 + 4 b 1 b 4 ) δ 2 105 b 1 δ ( 1 + 3 b 1 δ ) Q 0 2 W x v 0 2 H 2 W 2 v 0 2 ) ,
τ n | y = 0 = 3 b 1 μ v 0 70 H ( 140 b 2 + 4 b 1 2 315 b 2 3 + 4 b 1 b 4 δ 2 ) .

4. Discussion

In the present article, we have successfully implemented the recursive approach to linearize the full momentum equations describing the incompressible slow flow of a Maxwell fluid between two porous walls and analytical solutions for stream function, velocity components, total pressure difference, explicit expressions for normal and shear stresses, volumetric flow rate, fractional absorption and leakage flux are calculated. To predict the impact of pertinent parameters like absorption parameter K, slip parameter ϕ and Maxwell fluid parameter δ on flow variables, we first used the following dimensionless parameters and then we established different graphs after removing star along the three different places along slit length, i.e., entrance region (x = 0.1), middle region (x = 0.5) and exit region (x = 0.9).
x * = x H , ζ = y H , u * = u Q 0 / W H , v * = v Q 0 / W H , p * = p μ V 0 / H , τ * i j = τ i j μ V 0 / H K = 2 H W V 0 Q 0 , Q * ( x ) = Q ( x ) Q 0 .
Figure 2, Figure 3 and Figure 4 are drawn to see the impact of porosity K on u ( ζ ) for the creeping flow of Newtonian and non-Newtonian (Maxwell) fluids. It is justified that when we increase the value of porosity K, the wall permeability increases and more fluid passed through the slit walls which causes an increase in magnitude of the seepage velocity. Keeping this fact in mind, we see that these figures depicted a decrease in axial velocity u ( ζ ) on increasing porosity K and the magnitude of axial velocity also reduces downstream throughout the channel. For large values of K, a reversal flow may occur near the exit region of the slit. A standard parabolic velocity profile for Poiseuille flow in the absence of porosity throughout the channel is also observed in Figure 2. An increase is also observed in magnitude of axial velocity profile due to the presence of ϕ in Figure 4. A clear difference can be observed in velocity profiles for slow flow of Newtonian and Maxwell fluids in Figure 3 and Figure 4.
Figure 5 and Figure 6 show the impact of Maxwell fluid parameter δ on the axial velocity u ( ζ ) of creeping Maxwell fluid with and without slip parameter, keeping K = 0.2 . It is viewed that on increasing value of δ , the shear thickening property is observed along the centerline of the slit as axial velocity u ( ζ ) decreases along the centerline of the slit and due to reabsorption and slip effects on the slit wall, the shear thinning property of the fluid near the wall is observed as u ( ζ ) increases near the slit wall. An increase in the magnitude of u ( ζ ) can be seen due to parameter ϕ in Figure 6.
Figure 7 and Figure 8 reveal the impact of parameter ϕ on u ( ζ ) in Newtonian and non-Newtonian cases. It is viewed that the velocity of fluid at the boundary is not the same velocity with fluid particles closet to the wall. In both cases with increasing the slip parameter ϕ , u ( ζ ) increases. The variation in the radial velocity component v ( ζ ) can be seen in Figure 9 and we have observed that v ( ζ ) increases as we increased the parameters K and ϕ , whereas it is decreasing due to δ .
Figure 10, Figure 11 and Figure 12 are drawn to show the impact of pertinent parameters on pressure p ( x , y ) , shear stress τ w and normal stress difference τ n of a creeping Maxwell fluid. It is noted that in Figure 10, the pressure difference is increasing due to K and decreasing downstream on increasing δ and ϕ , whereas Figure 11 illustrate that the wall shear stress τ w has opposite effects as compared to pressure difference due to these parameters.
It is also noted in Figure 12a that porosity K does not affect the normal stress difference τ n away from boundaries but τ n decreases near the walls with K and also in Figure 12b,c τ n decreases on increasing δ and ϕ . Figure 13 shows the variation of Q ( x ) and it is depicted that Q ( x ) decreases as we increased porosity K whereas it is increasing as we increase parameters δ and ϕ .

5. Conclusions

In this article, we have presented the slow flow of a Maxwell fluid through the permeable slit under the influence of a slip condition and the recursive approach is consider to get an analytical solution. The key features of the current investigation were the impact of parameters K, ϕ and δ on flow variables involved. Moreover, the recursive approach is used for the first time to deal with a rate type Maxwell fluid creeping flow and the obtained results are very similar to existing results presented by [15] in a limiting case with certain restricted conditions. The following outcomes are drawn from this study:
  • If the Maxwell fluid parameter δ = 0 and slip parameter ϕ = 0 , then results obtained by Haroon [15] are recovered.
  • The axial velocity u ( ζ ) of creeping Maxwell fluid decreases downstream along the slit length on increasing porosity K and also for increasing values of K, backward flow can be seen near the exit region of the slit.
  • The axial velocity profile has an increasing behavior near the slit walls and its decreasing trend is observed along the centerline of the slit with increasing δ .
  • The shear thickening and thinning behavior of the Maxwell fluid is observed along centerline and near the walls of the slit, respectively.
  • Along the slit length, the magnitude of u ( ζ ) decreases as the fluid moves from the entrance to exit region of the slit.
  • An increase in axial velocity of a creeping Maxwell fluid is observed due to the increasing value of ϕ .
  • The slip parameter ϕ significantly influenced the magnitude of axial and radial velocities in comparison to other parameters.
  • A decreasing trend in pressure profile for increasing values of ϕ and δ , whereas pressure is increasing with increasing porosity parameter K.
  • The wall shear stress τ w is increasing significantly by increasing δ and ϕ , but with K and τ w decreasing.
  • The contribution of δ in axial flow rate and leakage flux is only due to the presence of a slip parameter ϕ .

Author Contributions

Conceptualization, H.U. and A.M.S.; methodology, H.U., D.L. and A.M.S.; software, H.U.; T.H.; validation, H.U., T.H. and K.M.; writing—original draft preparation, H.U.; writing—review and editing, H.U., A.M.S., K.M.; supervision, D.L.; project administration, D.L.; funding acquisition, D.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by China Postdoctoral Science Foundation grant number 2019M651715.

Acknowledgments

The authors are thankful to Jiangsu University for providing excellent research facilities.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
u, vComponents of velocity field
x, yCartesian coordinates
LLength of the slit
HWidth of slit
WBreadth of slit
V 0 Uniform velocity
QAxial flow rate at any point x
μ Coefficient of viscosity
λ fluid relaxation time
δ Maxwell fluid parameter (Deborah number)
KPorosity parameter
ϕ Slip parameter
ψ Stream function
p ( x , y ) Pressure in the slit
τ w Wall shear stress
τ n Normal stresses difference
q ( x ) Leakage flux
F a Fractional reabsorption

References

  1. Nikolay, V. Desalination Engineering: Planning and Design; McGraw-Hill Professional: New York, NY, USA, 2013. [Google Scholar]
  2. Espedal, M.S.; Mikelic, A. Filtration in Porous Media and Industrial Application: Lectures Given at the 4th Session of the Centro Internazionale Matematico Estivo (CIME) Held in Cetraro, Italy, 24–29 August 1998; Springer: New York, NY, USA, 2007. [Google Scholar]
  3. Macey, R.I. Pressure flow patterns in a cylinder with reabsorbing walls. Bull. Math. Biophys. 1963, 25, 1–9. [Google Scholar] [CrossRef]
  4. Macey, R.I. Hydrodynamics in the renal tubule. Bull. Math. Biophys. 1965, 27, 117. [Google Scholar] [CrossRef] [PubMed]
  5. Marshall, E.; Trowbridge, E. Flow of a Newtonian fluid through a permeable tube: The application to the proximal renal tubule. Bull. Math. Biol. 1974, 36, 457–476. [Google Scholar] [CrossRef]
  6. Marshall, E.; Trowbridge, E.; Aplin, A. Flow of a Newtonian fluid between parallel flat permeable plates, The application to a flat plate hemodialyzer. Math. Biosci. 1975, 27, 119–139. [Google Scholar] [CrossRef]
  7. Berman, A.S. Laminar flow in channels with porous walls. J. Appl. Phys. 1953, 24, 1232–1235. [Google Scholar] [CrossRef]
  8. Sellars, J.R. Laminar flow in channels with porous walls at high suction Reynolds numbers. J. Appl. Phys. 1955, 26, 489–490. [Google Scholar] [CrossRef] [Green Version]
  9. Yuan, S. Further investigation of laminar flow in channels with porous walls. J. Appl. Phys. 1956, 27, 267–269. [Google Scholar] [CrossRef]
  10. Wah, T. Laminar flow in a uniformly porous channel. Aeronaut. Q. 1964, 15, 299–310. [Google Scholar] [CrossRef]
  11. Terrill, R. Laminar flow in a uniformly porous channel with large injection. Aeronaut. Q. 1965, 16, 323–332. [Google Scholar] [CrossRef]
  12. Karode, S.K. Laminar flow in channels with porous walls revisited. J. Membr. Sci. 2001, 191, 237–241. [Google Scholar] [CrossRef]
  13. Siddiqui, A.M.; Haroon, T.; Shahzad, A. Hydrodynamics of viscous fluid through porous slit with linear absorption. Appl. Math. Mech. 2016, 37, 361–378. [Google Scholar] [CrossRef]
  14. Haroon, T.; Siddiqui, A.M.; Shahzad, A. Stokes flow through a slit with periodic reabsorption: An application to renal tubule. Alex. Eng. J. 2016, 55, 1799–1810. [Google Scholar] [CrossRef] [Green Version]
  15. Haroon, T.; Siddiqui, A.; Shahzad, A. Creeping flow of viscous fluid through a proximal tubule with uniform reabsorption: A mathematical study. Appl. Math. Sci. 2016, 10, 795–807. [Google Scholar] [CrossRef]
  16. Haroon, T.; Siddiqui, A.; Shahzad, A.; Smeltzer, J. Steady creeping slip flow of viscous fluid through a permeable slit with exponential reabsorption1. Appl. Math. Sci. 2017, 11, 2477–2504. [Google Scholar]
  17. Rajagopal, K. On the creeping flow of the second-order fluid. J. Non-Newton. Fluid Mech. 1984, 15, 239–246. [Google Scholar] [CrossRef]
  18. Ullah, H.; Sun, H.; Siddiqui, A.M.; Haroon, T. Creeping flow analysis of slightly non-Newtonian fluid in a uniformly porous slit. J. Appl. Anal. Comput. 2019, 9, 140–158. [Google Scholar] [CrossRef]
  19. Ullah, H.; Siddiqui, A.M.; Sun, H.; Haroon, T. Slip effects on creeping flow of slightly non-Newtonian fluid in a uniformly porous slit. J. Braz. Soc. Mech. Sci. Eng. 2019, 41, 412. [Google Scholar] [CrossRef]
  20. Kahshan, M.; Siddiqui, A.; Haroon, T. A micropolar fluid model for hydrodynamics in the renal tubule. Eur. Phys. J. Plus 2018, 133, 546. [Google Scholar] [CrossRef]
  21. Kahshan, M.; Lu, D.; Siddiqui, A. A Jeffrey fluid model for a porous-walled channel: Application to flat plate dialyzer. Sci. Rep. 2019, 9, 1–18. [Google Scholar] [CrossRef] [Green Version]
  22. Lu, D.; Kahshan, M.; Siddiqui, A. Hydrodynamical study of micropolar fluid in a porous-walled channel: Application to flat plate dialyzer. Symmetry 2019, 11, 541. [Google Scholar] [CrossRef] [Green Version]
  23. Langlois, W. A Recursive Approach to the Theory of Slow, Steady-State Viscoelastic Flow. Trans. Soc. Rheol. 1963, 7, 75–99. [Google Scholar] [CrossRef]
  24. Langlois, W. The recursive theory of slow viscoelastic flow applied to three basic problems of hydrodynamics. Trans. Soc. Rheol. 1964, 8, 33–60. [Google Scholar] [CrossRef]
  25. Léger, L.; Hervet, H.; Massey, G.; Durliat, E. Wall slip in polymer melts. J. Phys. Condens. Matter 1997, 9, 7719. [Google Scholar] [CrossRef]
  26. Atwood, B.; Schowalter, W. Measurements of slip at the wall during flow of high-density polyethylene through a rectangular conduit. Rheol. Acta 1989, 28, 134–146. [Google Scholar] [CrossRef]
  27. Beavers, G.S.; Joseph, D.D. Boundary conditions at a naturally permeable wall. J. Fluid Mech. 1967, 30, 197–207. [Google Scholar] [CrossRef]
  28. Saffman, P.G. On the boundary condition at the surface of a porous medium. Stud. Appl. Math. 1971, 50, 93–101. [Google Scholar] [CrossRef]
  29. Beavers, G.S.; Sparrow, E.M.; Magnuson, R.A. Experiments on coupled parallel flows in a channel and a bounding porous medium. J. Basic Eng. 1970, 92, 843–848. [Google Scholar] [CrossRef]
  30. Kohler, J.T. An Investigation of Laminar Flow through a Porous Walled Channel: I. A Perturbation Solution Assuming Slip at the Permeable Wall. II. An Experimental Measurement of the Velocity Distribution Utilizing a Dye Tracer Technique. Ph.D. Thesis, University of Massachusetts Amherst, Amherst, MA, USA, 1974. [Google Scholar]
  31. Mikelic, A.; Jäger, W. On the interface boundary condition of Beavers, Joseph, and Saffman. SIAM J. Appl. Math. 2000, 60, 1111–1127. [Google Scholar] [CrossRef]
  32. Rao, I.; Rajagopal, K. The effect of the slip boundary condition on the flow of fluids in a channel. Acta Mech. 1999, 135, 113–126. [Google Scholar] [CrossRef]
  33. Elshahed, M. Blood flow in capillary under starling hypothesis. Appl. Math. Comput. 2004, 149, 431–439. [Google Scholar] [CrossRef]
  34. Singh, R.; Laurence, R.L. Influence of slip velocity at a membrane surface on ultrafiltration performance I. Channel flow system. Int. J. Heat Mass Transf. 1979, 22, 721–729. [Google Scholar] [CrossRef]
  35. Makinde, O.; Osalusi, E. MHD steady flow in a channel with slip at the permeable boundaries. Rom. J. Phys. 2006, 51, 319. [Google Scholar]
  36. Eldesoky, I.M. Unsteady MHD pulsatile blood flow through porous medium in a stenotic channel with slip at the permeable walls subjected to time dependent velocity (injection/suction). In Proceedings of the International Conference on Mathematics and Engineering Physics, Kobry Elkobbah, Egypt, 29–31 May 2014; Volume 7, pp. 1–25. [Google Scholar]
  37. El-Shehawy, E.; El-Dabe, N.; El-Desoky, I. Slip effects on the peristaltic flow of a non-Newtonian Maxwellian fluid. Acta Mech. 2006, 186, 141–159. [Google Scholar] [CrossRef]
  38. Ellahi, R. Effects of the slip boundary condition on non-Newtonian flows in a channel. Commun. Nonlinear Sci. Numer. Simul. 2009, 14, 1377–1384. [Google Scholar] [CrossRef]
  39. Hron, J.; Le Roux, C.; Málek, J.; Rajagopal, K. Flows of incompressible fluids subject to Naviers slip on the boundary. Comput. Math. Appl. 2008, 56, 2128–2143. [Google Scholar] [CrossRef] [Green Version]
  40. Hayat, T.; Khan, M.; Ayub, M. The effect of the slip condition on flows of an Oldroyd 6-constant fluid. J. Comput. Appl. Math. 2007, 202, 402–413. [Google Scholar] [CrossRef] [Green Version]
  41. Rajagopal, K.R.; Srinivasa, A.R. A thermodynamic frame work for rate type fluid models. J. Non-Newton. Fluid Mech. 2000, 88, 207–227. [Google Scholar] [CrossRef]
  42. Choi, J.; Rusak, Z.; Tichy, J. Maxwell fluid suction flow in a channel. J. Non-Newton. Fluid Mech. 1999, 85, 165–187. [Google Scholar] [CrossRef]
  43. Sadeghy, K.; Najafi, A.H.; Saffaripour, M. Sakiadis flow of an upper-convected Maxwell fluid. Int. J. Non-Linear Mech. 2005, 40, 1220–1228. [Google Scholar] [CrossRef]
  44. Abbas, Z.; Sajid, M.; Hayat, T. MHD boundary-layer flow of an upper-convected Maxwell fluid in a porous channel. Theor. Comput. Fluid Dyn. 2006, 20, 229–238. [Google Scholar] [CrossRef]
  45. Bhatti, K.; Siddiqui, A.M.; Bano, Z. Application of Recursive Theory of Slow Viscoelastic Flow to the Hydrodynamics of Second-Order Fluid Flowing through a Uniformly Porous Circular Tube. Mathematics 2020, 8, 1170. [Google Scholar] [CrossRef]
Figure 1. Geometry of the problem.
Figure 1. Geometry of the problem.
Mathematics 08 01852 g001
Figure 2. Behavior of u ( ζ ) with K at (a) x = 0.1 , (b) x = 0.5 , (c) x = 0.9 , when ϕ = δ = 0.0 .
Figure 2. Behavior of u ( ζ ) with K at (a) x = 0.1 , (b) x = 0.5 , (c) x = 0.9 , when ϕ = δ = 0.0 .
Mathematics 08 01852 g002
Figure 3. Behavior of u ( ζ ) with K at (a) x = 0.1 , (b) x = 0.5 , (c) x = 0.9 , when ϕ = 0.2 , δ = 0.0 .
Figure 3. Behavior of u ( ζ ) with K at (a) x = 0.1 , (b) x = 0.5 , (c) x = 0.9 , when ϕ = 0.2 , δ = 0.0 .
Mathematics 08 01852 g003
Figure 4. Behavior of u ( ζ ) with K at (a) x = 0.1 , (b) x = 0.5 , (c) x = 0.9 , when ϕ = δ = 0.2 .
Figure 4. Behavior of u ( ζ ) with K at (a) x = 0.1 , (b) x = 0.5 , (c) x = 0.9 , when ϕ = δ = 0.2 .
Mathematics 08 01852 g004
Figure 5. Behavior of u ( ζ ) with δ at (a) x = 0.1 , (b) x = 0.5 , (c) x = 0.9 , when K = 0.2 , ϕ = 0.0 .
Figure 5. Behavior of u ( ζ ) with δ at (a) x = 0.1 , (b) x = 0.5 , (c) x = 0.9 , when K = 0.2 , ϕ = 0.0 .
Mathematics 08 01852 g005
Figure 6. Behavior of u ( ζ ) with δ at (a) x = 0.1 , (b) x = 0.5 , (c) x = 0.9 , when K = 0.2 , ϕ = 0.1 .
Figure 6. Behavior of u ( ζ ) with δ at (a) x = 0.1 , (b) x = 0.5 , (c) x = 0.9 , when K = 0.2 , ϕ = 0.1 .
Mathematics 08 01852 g006
Figure 7. Behavior of u ( ζ ) with ϕ at (a) x = 0.1 , (b) x = 0.5 , (c) x = 0.9 , when K = 0.2 , δ = 0.0 .
Figure 7. Behavior of u ( ζ ) with ϕ at (a) x = 0.1 , (b) x = 0.5 , (c) x = 0.9 , when K = 0.2 , δ = 0.0 .
Mathematics 08 01852 g007
Figure 8. Behavior of u ( ζ ) with ϕ at (a) x = 0.1 , (b) x = 0.5 , (c) x = 0.9 , when K = 0.2 , δ = 0.2 .
Figure 8. Behavior of u ( ζ ) with ϕ at (a) x = 0.1 , (b) x = 0.5 , (c) x = 0.9 , when K = 0.2 , δ = 0.2 .
Mathematics 08 01852 g008
Figure 9. Behavior of v ( ζ ) with (a) K when ϕ = δ = 0.2 , (b) δ when K = ϕ = 0.2 , (c) δ when K = ϕ = 0.2 .
Figure 9. Behavior of v ( ζ ) with (a) K when ϕ = δ = 0.2 , (b) δ when K = ϕ = 0.2 , (c) δ when K = ϕ = 0.2 .
Mathematics 08 01852 g009
Figure 10. Behavior of p ( x , 0 ) p ( 0 , 0 ) due to (a) K when ϕ = δ = 0.2 , (b) δ when K = δ = 0.2 , (c) ϕ when K = δ = 0.2 .
Figure 10. Behavior of p ( x , 0 ) p ( 0 , 0 ) due to (a) K when ϕ = δ = 0.2 , (b) δ when K = δ = 0.2 , (c) ϕ when K = δ = 0.2 .
Mathematics 08 01852 g010
Figure 11. Behavior of τ w with (a) K when ϕ = δ = 0.2 , (b) δ when K = δ = 0.2 , (c) ϕ when K = δ = 0.2 .
Figure 11. Behavior of τ w with (a) K when ϕ = δ = 0.2 , (b) δ when K = δ = 0.2 , (c) ϕ when K = δ = 0.2 .
Mathematics 08 01852 g011
Figure 12. Behavior of τ n due to (a) K when ϕ = δ = 0.2 , (b) δ when K = δ = 0.2 , (c) ϕ when K = δ = 0.2 .
Figure 12. Behavior of τ n due to (a) K when ϕ = δ = 0.2 , (b) δ when K = δ = 0.2 , (c) ϕ when K = δ = 0.2 .
Mathematics 08 01852 g012
Figure 13. Behavior of Q ( x ) due to (a) K when ϕ = δ = 0.2 , (b) δ when K = δ = 0.2 , (c) ϕ when K = δ = 0.2 .
Figure 13. Behavior of Q ( x ) due to (a) K when ϕ = δ = 0.2 , (b) δ when K = δ = 0.2 , (c) ϕ when K = δ = 0.2 .
Mathematics 08 01852 g013
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ullah, H.; Lu, D.; Siddiqui, A.M.; Haroon, T.; Maqbool, K. Hydrodynamical Study of Creeping Maxwell Fluid Flow through a Porous Slit with Uniform Reabsorption and Wall Slip. Mathematics 2020, 8, 1852. https://doi.org/10.3390/math8101852

AMA Style

Ullah H, Lu D, Siddiqui AM, Haroon T, Maqbool K. Hydrodynamical Study of Creeping Maxwell Fluid Flow through a Porous Slit with Uniform Reabsorption and Wall Slip. Mathematics. 2020; 8(10):1852. https://doi.org/10.3390/math8101852

Chicago/Turabian Style

Ullah, Hameed, Dianchen Lu, Abdul Majeed Siddiqui, Tahira Haroon, and Khadija Maqbool. 2020. "Hydrodynamical Study of Creeping Maxwell Fluid Flow through a Porous Slit with Uniform Reabsorption and Wall Slip" Mathematics 8, no. 10: 1852. https://doi.org/10.3390/math8101852

APA Style

Ullah, H., Lu, D., Siddiqui, A. M., Haroon, T., & Maqbool, K. (2020). Hydrodynamical Study of Creeping Maxwell Fluid Flow through a Porous Slit with Uniform Reabsorption and Wall Slip. Mathematics, 8(10), 1852. https://doi.org/10.3390/math8101852

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