# Effects of Slip and Inclined Magnetic Field on the Flow of Immiscible Fluids (Couple Stress Fluid and Jeffrey Fluid) in a Porous Channel

Authors: Punnamchandar Bitla; Fekadu Yemataw Sitotaw
Journal: Journal of Applied Mathematics
 (2022)

---

## Abstract

In this paper, we study the flow of two immiscible fluids namely, couple stress fluid and Jeffrey fluid in a porous channel. Instead of the classical no-slip conditions on the boundaries, we used slip boundary conditions, which are more realistic and meaningful. In addition, we used inclined magnetic field effects on the fluid flow. The couple stress fluid and Jeffrey fluid are flowing adjacent to each other in the region I and in the region II, respectively, of the horizontal porous channel. The nondimensionalized governing equations are solved analytically by using slip conditions at the lower and upper boundaries and interface conditions at the fluid-fluid interface. The analytical expressions for the velocity components in both regions are obtained in closed form. The effects of slip parameter, Hartmann number, couple stress parameter, Jeffrey parameter, angle of inclination, and Darcy number on velocity components in both regions are investigated. In the absence of slip, couple stress parameter, and Jeffrey parameters, limiting cases are obtained and discussed.


---

## Body

## 1. Introduction

As the classical Newtonian fluid theory is inadequate to explain the fluids with additives and other nonlinear stress and strain relationships that occur in real fluids, many non-Newtonian fluid theories are proposed. Non-Newtonian fluids are categorized as visco-elastic fluids, time-dependent fluids, and time-independent fluids [1, 2].Stokes [3] in 1966 proposed the couple stress theory by introducing the nonsymmetry of the stress tensor that shows the effects of couple stresses. This fluid theory gained the attention of researchers due to its major applications in lubrication theory, liquid crystals, blood flows, etc. A review of couple stress fluid dynamics was reported by Stokes [4] in 1984. The presence of couple stress parameter decreases the velocity of the fluid and other parameters effect on the velocity of a single couple stress fluid studied by Devakar et al. [5], Ahmad [6], and Jangili et al. [7].Other non-Newtonian fluid which has been attracted much by the researchers in view of its simplicity is known as Jeffrey fluid [8]. The Jeffrey fluid model is able to describe the characteristics of relaxation and retardation time. Jeffrey’s fluid model is a significant generalization of the Newtonian fluid model as the latter one can be deduced as a special case of the former. Several scholars have studied Jeffrey fluid flow under different assumptions and conditions. Akbar et al. [9] have studied the effect of thermal and velocity slip for Jeffrey fluid symmetric channel. Akram and Nadeem [10] studied analytically the influence of induced magnetic field on the peristaltic motion of a Jeffrey fluid in an asymmetric channel. Shail [11] studied theoretically the possibility of using a two-phase system to obtain increased flow rates in an electromagnetic pump. Meyer and Garder [12] analyzed the mechanics of two immiscible fluids in porous media. Bhattacharya [13] studied the flow of immiscible fluids in a rigid channel with a time-dependent pressure gradient.The study of immiscible fluids like Newtonian versus Newtonian, Newtonian versus Non-Newtonian, and/or Non-Newtonian versus Non-Newtonian fluids occurs in industries, manufacturing process, groundwater hydrology, etc. Recently, several researchers studied the flows of immiscible fluids, to mention a few, such as Newtonian fluid-Newtonian fluid or Newtonian-Non-Newtonian fluid (Berg et al. [14], Padma Devi and Srinivas [15], Yadav and Sneha [16, 17], Deo and Ansari [18], Matias et al. [19], Malashetty and Umavathi [20], Ramana Murthy and Srinivas [21], and Punnamchandar and Fekadu [22]) in channel configurations with slip and/or no-slip boundary conditions.In general, the no-slip boundary condition is used when the fluid flows past rigid boundaries. However, it has been accepted now that a number of fluids such as polymerics and additive fluids slip/stick-slip on rigid boundaries. To describe the slip characteristics of a fluid on the solid surface, “Navier introduced a more general boundary condition – that the fluid velocity component tangential to the solid surface, relative to the solid surface, is proportional to the shear stress on the fluid-solid interface” (Thompson and Troia [23]). “Recent advances in experimental and mathematical investigation and results have provided evidence in support of the slip condition and try to show its effects on different fluid model” (Akram et al. [24], Punnamchandar and Iyengar [25], Ashmawy [26], and Srinivasacharya and Bindu [27]).To the extent, the authors have surveyed the combined effects of slip and inclined magnetic field on the flow of immiscible couple stress fluid and Jeffrey fluids in a porous medium channel have not been studied. Hence, in this paper, the flow of two immiscible fluids in a porous medium channel with an inclined magnetic field along with the slip at the boundaries is analyzed. The analytical expressions for the velocity components in both regions are obtained in closed form. The effects of various parameters on the two velocity profiles are shown graphically and briefly discussed.


## 2. Mathematical Formulation of the Problem

Consider a fully developed laminar incompressible, immiscible couple stress, and Jeffrey fluids flow through a porous medium between two parallel horizontal plates which are 2 h distance apart. Thex-axis is taken in the midway direction of the porous channel, and the y-axis is taken normal to the porous channel (Figure 1).Figure 1
A schematic representation of the flow of two immiscible fluids.Letβ1 and β2 be the slip coefficients of the fluids at the lower and the upper plates of the channel, respectively. The fluid flow is driven by a common pressure gradient. The couple stress fluid in the region I (−h≤y≤0) has density ρ1, shear viscosity μ1, and couple stress viscosity η1. Similarly, the Jeffrey fluid in the region II (0≤y≤h) has density ρ2, shear viscosity μ2 and λ1 is the material parameter of a Jeffrey fluid. The porous medium is rigid, isotropic, and homogeneous with permeabilitiesk1 and k2 in the two regions, respectively.A uniform magnetic fieldBois applied at an inclined angleθ with respect to the positive y-direction. The couple stress fluid and Jeffrey fluid are assumed to be slightly conducting. Hence, the induced magnetic field is negligible in comparison to the applied uniform inclined magnetic field. The Lorentz force is the only body force acting on the fluid, and there are no body couples.For region I−h≤y≤0,following Stokes [3] and Skelland [1], the flow of couple stress fluid in the region I has an extra stress tensor compare with Navier - Stokes equation is governed by the differential equations given by
(1)μ1d2u1dy2−η1d4u1dy4−σ1Bo2cos2θu1−μ1k1u1=−dpdx.For region II0≤y≤h,following Chabra and Richardson [2], Akbar et al. [9], and Akram and Nadeem [10], the flow of Jeffrey fluid in region II has an extra stress tensor S=μ2/1+λ1γ̇+γ¨, γ̇ is the shear rate. Using the above assumptions and this extra stress tensor, the governing differential equation reduced to
(2)μ21+λ1d2u2dy2−σ2B02cos2θu2−μ2k2u2=−dpdx,where dp/dx represents the driving force i.e., the common pressure gradient, μid2ui/dy2 represents the force due to shear stresses, η1d4u1/dy4 represents the force due to couple stresses, σiBo2ui represent the Lorentz force due to applied magnetic field, μi/kiui represents the retarding force due to porous media in both the fluid regions.To determine the velocity field componentsu1y and u2y in region I and region II described above, we adopt the following boundary and interface conditions:
(i)
aty=−h, the Navier slip relation says that the shear stress on the fluid at the boundary is proportional to the fluid velocity at the boundary [28]. A velocity slip boundary condition proposed by Navier in which slip velocity depends upon shear stress. Devakar et al. [5] and Pei-Ying et al. [29] presented the analytical solution for a couple stress fluids by taking velocity slip conditions. This slip conditions is given by(3)u1=β1du1dy−η1μ1d3u1dy3.Stokes [3] has proposed two types of boundary Conditions (A) and (B), respectively, and the vanishing of couple stresses on the boundary is referred to as Condition (A). This condition is adopted here as this is appropriate in the present context. The vanishing of couple stresses is given by
(4)d2u1dy2=0.(ii)
aty=0, at the fluid–fluid interface, Murthy and Srinivas [21] assumed that the velocity, shear stress, and couple stress components are continuous. This implies:Continuity of velocities,(5)u10−=u20+,u1=u2.Vanishing of couple stresses,(6)Mxy10−=Mxy20+,d2u1dy2=0.Continuity of shear stresses,(7)τxy10−=τxy20+,μ1du1dy−η1d3u1dy3=μ21+λ1du2dy.(iii)
aty=h, the slip boundary condition for Jeffrey fluid, is given by (Ramesh [30])(8)u2=−β21+λ1du2dy,By introducing the nondimensional quantities,x∗=x/h,y∗=y/y,ui∗=ui/U,p∗=p/ρ1U2 where U is the maximum velocity of the fluid in the channel and i=1,2, we have the governing equations, boundaries, and interface conditions (in the nondimensionalized form after dropping ∗‘s) asFor region I−1≤y≤0,
(9)1s1d4u1dy4−d2u1dy2+M2+1Dau1=ReP.Similarly for region II0≤y≤1,
(10)11+λ1d2u2dy2−1KDa+nσnμM2u2=−nρnμReP.(i)
aty=−1,(11)u1=1α1du1dy−1s1d3u1dy3,d2u1dy2=0(ii)
aty=0,(12)u1=u2nμ1+λ1du2dy=du1dy−1s1d3u1dy3,d2u1dy2=0(iii)
aty=1,(13)u2=−11+λ1α2du2dy,whereP=−dp/dx is the constant pressure gradient, Re=ρ1Uh/μ1 is the Reynolds number, αi=h/βii=1,2 are the slip parameters, s1=μ1h2/η1 is the couple stress parameter, Ha=B0hσ/μ1 is the Hartmann number, M=HaCosθ, Da=k/h2 is the Darcy number, nσ=σ2/σ1 is the electric conductivity ratio, nμ=μ2/μ1 is the viscosity ratio, and nρ=ρ2/ρ1 is the density ratio.


## 3. Solution of the Problem

Region I (−1≤y≤0), from Equation (9), the velocity component differential equation in region I is given by
(14)u1y=C11coshλ11y+C12sinhλ11y+C13coshλ12y+C14sinhλ12y+P1where λ1i=−1i−1s1/2±1/2s12−4m1,i=1,2,
(15)m1=M2+1Das1;P1=RePs1m1.Region II (0≤y≤1), from Equation (10), the velocity component differential equation in region II is given by
(16)u2y=C21coshλ21y+C22sinhλ21y+P2,whereλ2i=−1i−1m2fori=1,2,
(17)m2=ησημM2+1KDa1+λ1,P2=nρnμm21+λ1ReP.The velocity componentsu1y and u2y involve six constants C11,C12,C13,C14,C21,andC22. The expressions of these are given in appendix A.Evaluation of volumetric flow rate : the expression of volumetric flow rate through the channel in dimensionless form can be evaluated using the following relation [31]:
(18)Qy=∫−10u1ydy+∫01u2ydy,Qy=sinhλ11λ11C11+1−coshλ11λ11C12+sinhλ12λ12C13+1−coshλ12λ12C14+P1+sinhλ21λ21C21+coshλ21−1λ21C22+P2.Evaluation of shearing stress : the expression for nondimensional shearing stress experienced by fluids in the upper and lower region of the horizontal porous channel [31], respectively, is given by:
(19)τxy1=du1dy−1s1d3u1dy3aty=−1,τxy2=−11+λ1du2dyaty=1.τxy1=C12coshλ11λ11−coshλ11λ113s1,C11−sinhλ11λ11+sinhλ11λ113s1,+C14coshλ12λ12−coshλ12λ123s1,+C13−sinhλ12λ12+sinhλ12λ123s1,τxy2=−11+λ1sinhλ21C21λ21,−11+λ1sinhλ21C22λ21.


## 4. Particular Solution for the Effects of Slip and Inclined Magnetic Field on the Immiscible Fluid Flow in a Horizontal Porous Channel

### 4.1. Immiscible Fluids (CSF and Newtonian Fluid) Flow in the Assumed Channel in the Existence of a Uniform Inclined Magnetic Field forλ1=0

Region I (−1≤y≤0), the governing differentialequation is given by
(20)1s1d4u1dy4−d2u1dy2+M2+1Dau1=ReP.Therefore, the general solution is given by(21)u1y=C11coshλ11y+C12sinhλ11y+C13coshλ12y+C14sinhλ12y+P1,whereλ1i=−1i−1s1/2±1/2s12−4m1fori=1,2.Region II (0≤y≤1), the governing differential equation is given by
(22)d2u2dy2−1KDa+nσnμM2u2=−nρnμReP.Therefore, the general solution is given by(23)u2y=C21coshλ21y+C22sinhλ21y+P2,whereλ21=m2;m1=M2+1/Das1,P1=RePs1/m1,(24)m2=ησημM2+1KDa,P2=nρnμm2ReP.


### 4.2. Immiscible Newtonian and Jeffery Fluids Flow in the Assumed Channel in the Existence of a Uniform Inclined Magnetic Field forη1=0 or s1⟶∞

Region I (-1≤y≤0), the governing differential equation is given by
(25)−d2u1dy2+Ha2+1Dau1=ReP.Therefore, the general solution is given by(26)u1y=C11coshλ11y+C12sinhλ11y+P1.Region II (0≤y≤1), the governing differential equation is given by
(27)11+λ1d2u2dy2−1KDa+nσnμM2u2=−nρnμReP.Therefore, the general solution is given by(28)u2y=C21coshλ21y+C22sinhλ21y+P2whereλ11=m1;λ2i=−1i−1m2fori=1,2,m1=M2+1/Da,m2=ησ/ημM2+1/KDa1+λ1,P1=ReP/m_1, and P2=n_ρ/n_μm_21+λ_1ReP.


### 4.3. Immiscible Newtonian and Newtonian Fluid Flow in the Channel in the Existence of Uniform Inclined Magnetic Field forη1=0 or s1⟶∞ and λ1=0

Region I (-1≤y≤0), the governing differential equation is given by
(29)−d2u1dy2+M2+1Dau1=ReP.Therefore, the general solution is given by(30)u1y=C11coshλ11y+C12sinhλ11y+P1.Region II (0≤y≤1), the governing differential equation is given by
(31)d2u2dy2−1KDa+nσnμM2u2=−nρnμReP.Therefore, the general solution is given by(32)u2y=C21coshλ21y+C22sinhλ21y+P2,whereλ11=m1,λ21=m2,m1=M2+1/Da,P1=ReP/m1,m2=ησ/ημM2+1/KDa,andP2=nρ/nμm2ReP.


### 4.4. Comparative Study with Kumar and Agarwal [32]

Assuming the two immiscible fluids are Newtonian fluids in such a way that the fluid in region I is nonconducting and the fluid in region II is conducting fluids. From our model equation given in Equations (1) and (2), let the couple stress viscosity is zero (η1=0 or s1⟶∞), and the Jeffrey fluid parameter is also zero (λ1=0), and the angle of inclination θ=0 (which implies M=Ha) then we have the following.Region I (-1≤y≤0), the differential equation is given by
(33)−d2u1dy2+Ha2+1Dau1=ReP.Therefore, the general solution is given by(34)u1y=C11coshλ11y+C12sinhλ11y+P1.Region II (0≤y≤1), the differential equation is given by
(35)d2u2dy2−1Dau2=−nρnμReP.Therefore, the general solution is given by(36)u2y=C21coshλ21y+C22sinhλ21y+P2,whereλ11=m1;λ21=m2;m1=Ha2+1/Da,
(37)P1=RePm1;m2=1Da;P2=nρnμm2ReP.Now by sustitutingthe following variables in to the equations(34) and (36) (P=Ps,Ha=M,Da=K,andnμ=1/β;(38)nρ=1α,m1=A1,m2=A2,Re=R1,C11=C1,C12=C2,C21=C3,andC22=C4), after simplification
(39)u1y=C1coshA1y+C2sinhA1y+Re1PsA1,(40)u2y=C3coshA2y+C4sinhA2y+αR1PsA2.These equations (39) and (40) are exactly the same with the work of Kumar and Agarwal [32] which is the solution of steady flow.


## 4.1. Immiscible Fluids (CSF and Newtonian Fluid) Flow in the Assumed Channel in the Existence of a Uniform Inclined Magnetic Field forλ1=0

Region I (−1≤y≤0), the governing differentialequation is given by
(20)1s1d4u1dy4−d2u1dy2+M2+1Dau1=ReP.Therefore, the general solution is given by(21)u1y=C11coshλ11y+C12sinhλ11y+C13coshλ12y+C14sinhλ12y+P1,whereλ1i=−1i−1s1/2±1/2s12−4m1fori=1,2.Region II (0≤y≤1), the governing differential equation is given by
(22)d2u2dy2−1KDa+nσnμM2u2=−nρnμReP.Therefore, the general solution is given by(23)u2y=C21coshλ21y+C22sinhλ21y+P2,whereλ21=m2;m1=M2+1/Das1,P1=RePs1/m1,(24)m2=ησημM2+1KDa,P2=nρnμm2ReP.


## 4.2. Immiscible Newtonian and Jeffery Fluids Flow in the Assumed Channel in the Existence of a Uniform Inclined Magnetic Field forη1=0 or s1⟶∞

Region I (-1≤y≤0), the governing differential equation is given by
(25)−d2u1dy2+Ha2+1Dau1=ReP.Therefore, the general solution is given by(26)u1y=C11coshλ11y+C12sinhλ11y+P1.Region II (0≤y≤1), the governing differential equation is given by
(27)11+λ1d2u2dy2−1KDa+nσnμM2u2=−nρnμReP.Therefore, the general solution is given by(28)u2y=C21coshλ21y+C22sinhλ21y+P2whereλ11=m1;λ2i=−1i−1m2fori=1,2,m1=M2+1/Da,m2=ησ/ημM2+1/KDa1+λ1,P1=ReP/m_1, and P2=n_ρ/n_μm_21+λ_1ReP.


## 4.3. Immiscible Newtonian and Newtonian Fluid Flow in the Channel in the Existence of Uniform Inclined Magnetic Field forη1=0 or s1⟶∞ and λ1=0

Region I (-1≤y≤0), the governing differential equation is given by
(29)−d2u1dy2+M2+1Dau1=ReP.Therefore, the general solution is given by(30)u1y=C11coshλ11y+C12sinhλ11y+P1.Region II (0≤y≤1), the governing differential equation is given by
(31)d2u2dy2−1KDa+nσnμM2u2=−nρnμReP.Therefore, the general solution is given by(32)u2y=C21coshλ21y+C22sinhλ21y+P2,whereλ11=m1,λ21=m2,m1=M2+1/Da,P1=ReP/m1,m2=ησ/ημM2+1/KDa,andP2=nρ/nμm2ReP.


## 4.4. Comparative Study with Kumar and Agarwal [32]

Assuming the two immiscible fluids are Newtonian fluids in such a way that the fluid in region I is nonconducting and the fluid in region II is conducting fluids. From our model equation given in Equations (1) and (2), let the couple stress viscosity is zero (η1=0 or s1⟶∞), and the Jeffrey fluid parameter is also zero (λ1=0), and the angle of inclination θ=0 (which implies M=Ha) then we have the following.Region I (-1≤y≤0), the differential equation is given by
(33)−d2u1dy2+Ha2+1Dau1=ReP.Therefore, the general solution is given by(34)u1y=C11coshλ11y+C12sinhλ11y+P1.Region II (0≤y≤1), the differential equation is given by
(35)d2u2dy2−1Dau2=−nρnμReP.Therefore, the general solution is given by(36)u2y=C21coshλ21y+C22sinhλ21y+P2,whereλ11=m1;λ21=m2;m1=Ha2+1/Da,
(37)P1=RePm1;m2=1Da;P2=nρnμm2ReP.Now by sustitutingthe following variables in to the equations(34) and (36) (P=Ps,Ha=M,Da=K,andnμ=1/β;(38)nρ=1α,m1=A1,m2=A2,Re=R1,C11=C1,C12=C2,C21=C3,andC22=C4), after simplification
(39)u1y=C1coshA1y+C2sinhA1y+Re1PsA1,(40)u2y=C3coshA2y+C4sinhA2y+αR1PsA2.These equations (39) and (40) are exactly the same with the work of Kumar and Agarwal [32] which is the solution of steady flow.


## 5. Results and Discussion

### 5.1. Velocity Profile

The flow of immiscible fluids (couple stress fluid and Jeffery fluid) in a horizontal porous channel in the presence of a uniform inclined magnetic field is considered. The simplified model equations are solved analytically, and the expressions for velocities in both are obtained in closed form. The effects of pertinent parameters on the velocityu1y in region I and on the velocity u2y in region II are studied and presented through Figures 2–12. The set of fixed values of all parameters (Ha=3,Re=2,P=1.1,nσ=1.0,nμ=0.9,nρ=0.9,a1=10,a2=10,θ=π/4,Da=1,s1=2,K=1,andλ1=0.3) is carefully chosen in accordance with previous literatures (see [30]) when a sparticular parameter is varied to see the variation. Further, if any, variation in the parameters is indicated in that particular figure.Figure 2
Effect of couple stress parameterS1 on velocity components.Figure 3
Variation of velocity profile with Jeffrey parameterλ1.Figure 4
Effect of Hartmann numberHa on velocity components.Figure 5
Variation of velocities with the angle of inclinationθ.Figure 6
Effect of Darcy numberDa=k/h2 on velocity components.Figure 7
Effect of viscosity ratioημ=μ2/μ1 on velocity components.Figure 8
Effect of slip parameterα=α1=α2 on velocity components.Figure 9
Variation of velocity with electrical conductivity ratio.Figure 10
Variation of velocity of Newtonian and Jeffrey fluids with Jeffrey fluid parameterλ1.Figure 11
Velocity variation of couple stress and Newtonian fluid with couple stress fluid parameters1.Figure 12
Velocity variation of two immiscible Newtonian fluids with Hartmann numberHa.The effect of the couple stress parameters1 on the velocity components is shown in Figure 2. It is seen that as s1=μ1h2/η1 is increasing (i.e., as η is decreasing); the velocity components are increasing. We conclude from Figure 2 that the velocity in the case of a couple stress fluids is smaller than that in the Newtonian fluid case. Thus, the presence of couple stresses in the fluid decreases the velocity. This may be due to the fact that the couple stresses spend some energy to rotate the particles, thereby decreasing the particle velocity. We observe that the flow velocities are strongly dependent on the non-Newtonian material parameter s1. As s1⟶∞ (asη⟶0); the velocity component u2y in region I corresponds to the Newtonian fluid. This particular velocity profile represents the flow of immiscible fluids (Newtonian fluid and Jeffrey fluid) in a porous medium with slip boundary conditions.The effects of Jeffrey parameterλ1 on the velocity profile is shown in Figure 3. It is observed that on increasing Jeffrey parameter λ1, the velocity of the fluid in the region II increases. Further, the maximum values of fluid velocity occur at the interface for small Jeffrey parameter λ1, and it shifts to region II for large Jeffrey parameter λ1. This is because of the force due to shear stresses decreases when the Jeffrey parameter λ1 increases. We observe that the fluid velocity in the region II is strongly dependent on the non-Newtonian material parameter λ1. Further, it is seen that, as λ1 approaches zero, the flow corresponds to that of immiscible couple stress fluid and Newtonian fluid in a porous channel.The effect of Hartmann number on the distribution of velocity components is as shown in Figure4. The velocity profile attains the maximum for Ha=0; this indicates the fluid flow in the absence of a magnetic field attains the maximum velocity. As Ha is increasing, it is seen that the velocity is decreasing. This shows that the imposed magnetic field has a retarding effect on the flow. Therefore, increasing the Hartmann number increases the effect of resistive force, which further decreases the flow velocity.In Figure5, the effect of the angle of inclination of the magnetic field on the velocity profile is shown. We observe that the velocity profile increases with the increasing of the angle inclination θ relative to y-axis of the applied magnetic field. This is due to the fact that the magnetic field has maximum effect in the transverse direction of the flow to retard the motion of the fluids.From Figure6, it is seen that as the permeability parameter Da=k/h2 is increasing (i.e., as the permeability k is increasing), the velocity is increasing, and the increase of the permeability of the porous medium reduces the drag force and hence causes the flow velocity to increase. As Da⟶∞, the velocity corresponds to the immiscible fluids (couple stress fluid and Jeffrey fluid) flow between two parallel plates in the absence of porous media.The effect of the shear viscosity ratioημ=μ2/μ1 is depicted in Figure 7. As ημ is increasing, the velocity is decreasing. For ημ<1.00 the velocity in region II is more than in region I. It confirms that the fluid with less viscosity has more velocity. For ημ≥1.00, the velocity components are symmetric about the interface. The maximum velocity of the fluid occurs at the interface of the regions when ημ≥1.00, while the velocity of the fluid in region II first increases and gains maximum value and then decreases toward the wall of the plate of region II.Figure8 shows the effect of slip parameter α=h/β1=h/β2=α1=α2 on the velocity profile. As α is increasing, the velocity is decreasing in both regions at the same rate, and hence the velocity profiles become flatten due to the decreasing shearing force from the slip boundary. Further, when θ=0 andα⟶∞β⟶0 represents the flow with no-slip boundary conditions which is similar result with Devakar [5] and Pei-Ying et al. [29].The influence of electrical conductivity rationσ=σ2/σ1 on the velocity profile is shown in Figure 9. It is seen that the velocity of fluids in both regions decreases with the increase in the electrical conductivity ratio. The increasing in conductivity ratio means increasing the resistive force to the fluid flow. It is important to note that the velocity of the fluid in region II is greater than the velocity of the fluid in region I, and the maximum velocity occurs in region II for nσ=σ2/σ1≤1.25.


### 5.2. Limiting Cases

(i)
if the couple stress parameters1⟶∞ or η1=0, then the problem corresponds to the flow of two immiscible Newtonian fluid and Jeffrey fluid flow in the porous channel with aunifrom inclined magneticfield(ii)
if Jeffrey fluid parameterλ1=0, then the problem corresponds to the flow of two immiscible couple stress and Newtonian fluids through porous channel with a unifrom inclined magnetic field(iii)
if the couple stress parameters1⟶∞ or η1=0 and Jeffrey fluid parameter λ1=0, then the problem corresponds to the flow of two immiscible Newtonian fluids flow in the porous channel with a unifrom inclined magnetic fieldWe present some results in the absence of couple stress and/or Jeffrey fluid parameters through Figures10–12. In Figure 10, as the Jeffrey fluid parameter λ1increases, the Jeffrey fluid velocity is more than the Newtonian fluids and λ1=0; the two fluids are Newtonian fluids, and hence the two velcities are symmetric with respect the fluid-fluid interface. In Figure 11 we observe the couple stress fluid velocity lesser than Newtonian fluid velocity.. The effect of Hartmann number Ha on the distributions of velocity profiles of the two immiscible Newtonian fluids is shown in Figure 12. As Ha increases, the velocity decreases. The velocity profiles attain maximum value for Ha=0, this indicates the fluid flow in the absence of a magnetic field attains the maximum velocity. Furthermore, the profile is symmetric with respect to the fluid-fluid interface.


### 5.3. Flow Rate and Shear Stresses

The variations of volumetric flow rate and shear stresses on the upper and lower wall of the channel with different physical parameters are shown in Tables1–5.Table 1
Variation of volumetric flow rate with conductivity ratio (nσ), slip parameter (α), Hartmann number (Ha) and viscosity ratio (nμ).
nσQα1=α2=αQ0.251.0024750.88490.750.89079100.845441.250.80541150.829981.750.73793250.8166HaQnμQ01.292610.251.176040.51.187360.51.0041711.0981510.8158820.95511.50.70028Table 2
Variation of volumetric flow rate with Darcy number (Da), angle of inclination (θ), and Jeffery fluid parameter (λ1).
DaQθQλ1Q10.84544θ=00.6300900.8097520.94787θ=π/60.7217610.8999830.98787θ=π/31.0215720.9463241.00917θ=π/21.2926130.97628Table 3
Variation of shear wall stress with Hartmann number (Ha), viscosity ratio (nμ), slip parameter (α) and conductivity ratio (nσ) at the lower plate τxy1 and at the upper plate τxy2.
Haτxy1τxy2nμτxy1τxy201.664370.143760.251.273120.20650.51.560110.133410.51.253620.1412111.471240.124610.751.229810.111251.51.394520.1170511.207230.09331ατxy1τxy2nστxy1τxy251.22119-0.07220.251.2661-0.0702101.21602-0.0120.751.23056-0.1082151.21399-0.00721.251.203120.00366201.21291-0.0051.751.181210.07319Table 4
Variation of shear stress with Jeffery fluid parameter (λ1) and angle of inclination (θ) at the lower plate τxy1 and at the upper plate τxy2.
λ1τxy1τxy2θτxy1τxy201.20383-0.065800.991510.0184611.23084-0.0431π/61.08820.0120921.23899-0.0821π/31.39452-0.014531.24184-0.1331π/21.66437-0.0779Table 5
Variation of shear wall stress Darcy number (Da) and density ratio (nρ) at the lower plate τxy1 and at the upper plate τxy2.
Daτxy1τxy2nρτxy1τxy211.21602-0.03230.51.157130.0743921.32484-0.01911.01.23074-0.00731.36709-0.03371.51.30434-0.177441.38956-0.03122.01.37795-0.3365From Table1, the volumetric flow rate in the horizontal porous channel decreases with the increase of conductivity ratio (nσ), slip parameter (α), Hartmann number (Ha), and viscosity ratio (nμ) in immiscible couple stress and Jeffery fluid flow. This is obvious, for example, as the magnetic parameter increases, the velocity decreases which in turn lead to the decrease in the volumetric flow rate. The effects of the other pertinent parameters on the volumetric flow rate are in tune with the velocity.From Table2, the volumetric flow rate in the porous channel increases with the increase of Darcy number (Da), angle of inclination (θ), and Jeffery fluid parameter (λ1). This is obvious, for example, as Darcy number increases, the velocity increases which lead to the increase in the volumetric flow rate. The effects of the other pertinent parameters on the volumetric flow rate are in tune with the velocity.From Table3, the shear wall stress of immiscible couple stress and Jeffery fluid flow in the horizontal porous channel decreases with the increase magnetic field (M), viscosity ratio (nμ), slip parameter (α=α1=α2), and conductivity ratio (nσ). It is also noticed that the shear stress on the wall of the lower region is greater than the shear stress on the wall of the upper region of the horizontal porous channel.From Table4, with the increasing of Jeffery parameter (λ1) and angle of inclination (θ) the shear stress is increasing at the lower wall and is decreasing at the upper wall. It is also noticed that the shear stress on the wall of the lower region is greater than the shear stress on the wall of the upper region of the horizontal porous channel.From Table5, the increasing of Darcy number (Da) and density ratio (nρ) makes the shear stress on the wall of the region increase and decrease in the upper region. It is also noticed that the shear stress on the wall of the lower region is greater than the shear wall stress on the wall of the upper region of the horizontal porous channel.


## 5.1. Velocity Profile

The flow of immiscible fluids (couple stress fluid and Jeffery fluid) in a horizontal porous channel in the presence of a uniform inclined magnetic field is considered. The simplified model equations are solved analytically, and the expressions for velocities in both are obtained in closed form. The effects of pertinent parameters on the velocityu1y in region I and on the velocity u2y in region II are studied and presented through Figures 2–12. The set of fixed values of all parameters (Ha=3,Re=2,P=1.1,nσ=1.0,nμ=0.9,nρ=0.9,a1=10,a2=10,θ=π/4,Da=1,s1=2,K=1,andλ1=0.3) is carefully chosen in accordance with previous literatures (see [30]) when a sparticular parameter is varied to see the variation. Further, if any, variation in the parameters is indicated in that particular figure.Figure 2
Effect of couple stress parameterS1 on velocity components.Figure 3
Variation of velocity profile with Jeffrey parameterλ1.Figure 4
Effect of Hartmann numberHa on velocity components.Figure 5
Variation of velocities with the angle of inclinationθ.Figure 6
Effect of Darcy numberDa=k/h2 on velocity components.Figure 7
Effect of viscosity ratioημ=μ2/μ1 on velocity components.Figure 8
Effect of slip parameterα=α1=α2 on velocity components.Figure 9
Variation of velocity with electrical conductivity ratio.Figure 10
Variation of velocity of Newtonian and Jeffrey fluids with Jeffrey fluid parameterλ1.Figure 11
Velocity variation of couple stress and Newtonian fluid with couple stress fluid parameters1.Figure 12
Velocity variation of two immiscible Newtonian fluids with Hartmann numberHa.The effect of the couple stress parameters1 on the velocity components is shown in Figure 2. It is seen that as s1=μ1h2/η1 is increasing (i.e., as η is decreasing); the velocity components are increasing. We conclude from Figure 2 that the velocity in the case of a couple stress fluids is smaller than that in the Newtonian fluid case. Thus, the presence of couple stresses in the fluid decreases the velocity. This may be due to the fact that the couple stresses spend some energy to rotate the particles, thereby decreasing the particle velocity. We observe that the flow velocities are strongly dependent on the non-Newtonian material parameter s1. As s1⟶∞ (asη⟶0); the velocity component u2y in region I corresponds to the Newtonian fluid. This particular velocity profile represents the flow of immiscible fluids (Newtonian fluid and Jeffrey fluid) in a porous medium with slip boundary conditions.The effects of Jeffrey parameterλ1 on the velocity profile is shown in Figure 3. It is observed that on increasing Jeffrey parameter λ1, the velocity of the fluid in the region II increases. Further, the maximum values of fluid velocity occur at the interface for small Jeffrey parameter λ1, and it shifts to region II for large Jeffrey parameter λ1. This is because of the force due to shear stresses decreases when the Jeffrey parameter λ1 increases. We observe that the fluid velocity in the region II is strongly dependent on the non-Newtonian material parameter λ1. Further, it is seen that, as λ1 approaches zero, the flow corresponds to that of immiscible couple stress fluid and Newtonian fluid in a porous channel.The effect of Hartmann number on the distribution of velocity components is as shown in Figure4. The velocity profile attains the maximum for Ha=0; this indicates the fluid flow in the absence of a magnetic field attains the maximum velocity. As Ha is increasing, it is seen that the velocity is decreasing. This shows that the imposed magnetic field has a retarding effect on the flow. Therefore, increasing the Hartmann number increases the effect of resistive force, which further decreases the flow velocity.In Figure5, the effect of the angle of inclination of the magnetic field on the velocity profile is shown. We observe that the velocity profile increases with the increasing of the angle inclination θ relative to y-axis of the applied magnetic field. This is due to the fact that the magnetic field has maximum effect in the transverse direction of the flow to retard the motion of the fluids.From Figure6, it is seen that as the permeability parameter Da=k/h2 is increasing (i.e., as the permeability k is increasing), the velocity is increasing, and the increase of the permeability of the porous medium reduces the drag force and hence causes the flow velocity to increase. As Da⟶∞, the velocity corresponds to the immiscible fluids (couple stress fluid and Jeffrey fluid) flow between two parallel plates in the absence of porous media.The effect of the shear viscosity ratioημ=μ2/μ1 is depicted in Figure 7. As ημ is increasing, the velocity is decreasing. For ημ<1.00 the velocity in region II is more than in region I. It confirms that the fluid with less viscosity has more velocity. For ημ≥1.00, the velocity components are symmetric about the interface. The maximum velocity of the fluid occurs at the interface of the regions when ημ≥1.00, while the velocity of the fluid in region II first increases and gains maximum value and then decreases toward the wall of the plate of region II.Figure8 shows the effect of slip parameter α=h/β1=h/β2=α1=α2 on the velocity profile. As α is increasing, the velocity is decreasing in both regions at the same rate, and hence the velocity profiles become flatten due to the decreasing shearing force from the slip boundary. Further, when θ=0 andα⟶∞β⟶0 represents the flow with no-slip boundary conditions which is similar result with Devakar [5] and Pei-Ying et al. [29].The influence of electrical conductivity rationσ=σ2/σ1 on the velocity profile is shown in Figure 9. It is seen that the velocity of fluids in both regions decreases with the increase in the electrical conductivity ratio. The increasing in conductivity ratio means increasing the resistive force to the fluid flow. It is important to note that the velocity of the fluid in region II is greater than the velocity of the fluid in region I, and the maximum velocity occurs in region II for nσ=σ2/σ1≤1.25.


## 5.2. Limiting Cases

(i)
if the couple stress parameters1⟶∞ or η1=0, then the problem corresponds to the flow of two immiscible Newtonian fluid and Jeffrey fluid flow in the porous channel with aunifrom inclined magneticfield(ii)
if Jeffrey fluid parameterλ1=0, then the problem corresponds to the flow of two immiscible couple stress and Newtonian fluids through porous channel with a unifrom inclined magnetic field(iii)
if the couple stress parameters1⟶∞ or η1=0 and Jeffrey fluid parameter λ1=0, then the problem corresponds to the flow of two immiscible Newtonian fluids flow in the porous channel with a unifrom inclined magnetic fieldWe present some results in the absence of couple stress and/or Jeffrey fluid parameters through Figures10–12. In Figure 10, as the Jeffrey fluid parameter λ1increases, the Jeffrey fluid velocity is more than the Newtonian fluids and λ1=0; the two fluids are Newtonian fluids, and hence the two velcities are symmetric with respect the fluid-fluid interface. In Figure 11 we observe the couple stress fluid velocity lesser than Newtonian fluid velocity.. The effect of Hartmann number Ha on the distributions of velocity profiles of the two immiscible Newtonian fluids is shown in Figure 12. As Ha increases, the velocity decreases. The velocity profiles attain maximum value for Ha=0, this indicates the fluid flow in the absence of a magnetic field attains the maximum velocity. Furthermore, the profile is symmetric with respect to the fluid-fluid interface.


## 5.3. Flow Rate and Shear Stresses

The variations of volumetric flow rate and shear stresses on the upper and lower wall of the channel with different physical parameters are shown in Tables1–5.Table 1
Variation of volumetric flow rate with conductivity ratio (nσ), slip parameter (α), Hartmann number (Ha) and viscosity ratio (nμ).
nσQα1=α2=αQ0.251.0024750.88490.750.89079100.845441.250.80541150.829981.750.73793250.8166HaQnμQ01.292610.251.176040.51.187360.51.0041711.0981510.8158820.95511.50.70028Table 2
Variation of volumetric flow rate with Darcy number (Da), angle of inclination (θ), and Jeffery fluid parameter (λ1).
DaQθQλ1Q10.84544θ=00.6300900.8097520.94787θ=π/60.7217610.8999830.98787θ=π/31.0215720.9463241.00917θ=π/21.2926130.97628Table 3
Variation of shear wall stress with Hartmann number (Ha), viscosity ratio (nμ), slip parameter (α) and conductivity ratio (nσ) at the lower plate τxy1 and at the upper plate τxy2.
Haτxy1τxy2nμτxy1τxy201.664370.143760.251.273120.20650.51.560110.133410.51.253620.1412111.471240.124610.751.229810.111251.51.394520.1170511.207230.09331ατxy1τxy2nστxy1τxy251.22119-0.07220.251.2661-0.0702101.21602-0.0120.751.23056-0.1082151.21399-0.00721.251.203120.00366201.21291-0.0051.751.181210.07319Table 4
Variation of shear stress with Jeffery fluid parameter (λ1) and angle of inclination (θ) at the lower plate τxy1 and at the upper plate τxy2.
λ1τxy1τxy2θτxy1τxy201.20383-0.065800.991510.0184611.23084-0.0431π/61.08820.0120921.23899-0.0821π/31.39452-0.014531.24184-0.1331π/21.66437-0.0779Table 5
Variation of shear wall stress Darcy number (Da) and density ratio (nρ) at the lower plate τxy1 and at the upper plate τxy2.
Daτxy1τxy2nρτxy1τxy211.21602-0.03230.51.157130.0743921.32484-0.01911.01.23074-0.00731.36709-0.03371.51.30434-0.177441.38956-0.03122.01.37795-0.3365From Table1, the volumetric flow rate in the horizontal porous channel decreases with the increase of conductivity ratio (nσ), slip parameter (α), Hartmann number (Ha), and viscosity ratio (nμ) in immiscible couple stress and Jeffery fluid flow. This is obvious, for example, as the magnetic parameter increases, the velocity decreases which in turn lead to the decrease in the volumetric flow rate. The effects of the other pertinent parameters on the volumetric flow rate are in tune with the velocity.From Table2, the volumetric flow rate in the porous channel increases with the increase of Darcy number (Da), angle of inclination (θ), and Jeffery fluid parameter (λ1). This is obvious, for example, as Darcy number increases, the velocity increases which lead to the increase in the volumetric flow rate. The effects of the other pertinent parameters on the volumetric flow rate are in tune with the velocity.From Table3, the shear wall stress of immiscible couple stress and Jeffery fluid flow in the horizontal porous channel decreases with the increase magnetic field (M), viscosity ratio (nμ), slip parameter (α=α1=α2), and conductivity ratio (nσ). It is also noticed that the shear stress on the wall of the lower region is greater than the shear stress on the wall of the upper region of the horizontal porous channel.From Table4, with the increasing of Jeffery parameter (λ1) and angle of inclination (θ) the shear stress is increasing at the lower wall and is decreasing at the upper wall. It is also noticed that the shear stress on the wall of the lower region is greater than the shear stress on the wall of the upper region of the horizontal porous channel.From Table5, the increasing of Darcy number (Da) and density ratio (nρ) makes the shear stress on the wall of the region increase and decrease in the upper region. It is also noticed that the shear stress on the wall of the lower region is greater than the shear wall stress on the wall of the upper region of the horizontal porous channel.


## 6. Conclusion

Effects of slip and inclined uniform magnetic field on immiscible fluids (couple stress fluid and Jeffrey fluid) flow in a porous channel are studied. We observe that, as the Hartmann number, slip parameter and viscosity ratios are increasing, the velocity components are decreasing in both regions. The velocity increases in region II as the Jeffery parameter is increasing, as the velocity in region I is constant. As the couple stress parameter is increasing, the velocity is increasing. As the Darcy numberDa and the Reynolds number Re are increasing; the velocity is increasing. The shear stress on the wall of the lower region is greater than the shear stress on the wall of the upper region of the porous channel.



---

*Source: 2799773-2022-09-22.xml*