# Partial Slip Impact on Double Diffusive Convection Flow of Magneto-Carreau Nanofluid through Inclined Peristaltic Asymmetric Channel

Authors: Safia Akram; Maria Athar; Khalid Saeed; Taseer Muhammad; Mir Yasir Umair
Journal: Mathematical Problems in Engineering
 (2021)

---

## Abstract

The significance of partial slip on double diffusive convection on magneto-Carreau nanofluid through inclined peristaltic asymmetric channel is examined in this paper. The two-dimensional and directional flow of a magneto-Carreau nanofluid is mathematically described in detail. Under the lubrication technique, the proposed model is simplified. The solutions of extremely nonlinear partial differential equations are calculated using a numerical technique. Graphical data are displayed using Mathematica software and Matlab to examine how temperature, pressure rise, concentration, pressure gradient, velocity profile, nanoparticle volume fraction, and stream functions behave on emerging parameters. It is noticed that as the velocity slip parameter is increased, the axial velocity at the channel’s center increases. Additionally, near the boundary, opposite behavior is observed. The temperature, concentration, and nanoparticle profile drops by increasing thermal slip, concentration slip, and nanoparticle slip parameter.


---

## Body

## 1. Introduction

Greek word “peristaltikos” is the origin of the word “peristalsis” that refers to being clasped and compressed. It is further employed to refer to a progressive contracting wave within a tube which may have a variable cross-sectional area. It is utilized by the bodies in physiology, to push and blend the guts inside a tube such as in uretra and other sorts of glandular vessels. The same principle is used by the roller and finger pumping devices. The theory of peristaltic transport has been deployed by multiple industries on different applications that may include sanitary ducts, blood pumping devices in the lungs. Ever since the most initial research on this topic by Latham [1], few other studies based on numerical models and experiments have also been reported under multiple conditions.In the recent times, the researchers in the domain of physiology have experienced a peristaltic kind of motion by the intrauterine duct flows. These are caused by myometrial contractions in both symmetric and asymmetric directions. Pozrikidis [2] developed a technique that was based on integral methods to further explore the peristaltic ducts in channels for the Stokes flow conditions to comprehend the process of fluid dynamics. He investigated the patterns and mean flowing rates caused by the amplitude and phase once the walls were deformed. More experiments were done by Eytan et al. [3] to characterize contractions in women who were not pregnant and found those complications as those were of different amplitudes and with a wide range of frequencies with variable wavelengths. There has been an identification that uterine cavity’s cross section breadth tends to increase to the fundus, and it does not occlude while contracting. Similarly, Eytan and Elad [4] managed to develop a model in 2D channels where wave trains tend to move with independence on the upper walls as well as the lower ones. The orientations have been deployed to further perform evaluations in nonpregnant women’s uterus. Calculations have been done to further calculate possible particle movements to better comprehend the process an embryo goes through. Mishra and Ramachandra Rao [5] further investigated the peristaltic transport for Newtonian fluids.The MHD flows in a channel holds much of importance with linkage to specific concerns related to the mobility of conductive fluids such as blood. This brings a need for research based on theory for the operational MHD compressor. Stud et al. [6] investigated the effects of a changing magnetic field on blood flowing channels. They managed to identify the fact that suitable field can accelerate the speed of the blood in the channel. Srivastava and Agarwal [7] studied other aspects such as after considering blood as conducting channel having a feature to suspend red blood cells in the plasma. Further studies [8–12] have also investigated other aspects related to MHD flows on peristalsis pumping.These days, nanotechnology is revolutionizing industrial era in the current century. For the same reason, many researchers have directed their work on modeling of such thermal conductivities and have investigated multiple nanofluids having different viscosities. Such fluids have suspended nanoparticles which are generally made of different metals, oxides, or carbon-based nanotubes. Most common ones include water, oil, and toluene. Choi [13] might have been the first researcher to coin the term nanofluid, where he mentioned that heat transfer mechanism is going to be most beneficial applications of the nanofluids. Further in [14], he explained that after adding a very meager number of nanoparticles to basic liquids, thermal conduction can be increased approximately double. Similarly, Masuda et al. [15], Lee et al. [16], Xuan and Li [17], and Xuan and Roetzel [18] showed the same enhancements by more than 20%. Therefore, for heat transfer mechanisms, nanofluids can be taken up as an interesting replacement [17]. In the recent times, studies have been done for nanofluids with peristalsis as the importance of this domain linked with biosciences and chemical and mechanical engineering cannot be undermined. Such mentions are referred as [19–23].In many industrial applications and engineering processes such as central receivers run on solar energy for creating potential through wind, cooling processes of electrical apparatus have collective forced and free convection in them. Transfer mechanism of heat and masses has effects on each other which also can yield cross-diffusion influences. The delta between the temperatures can generate mass transferring mechanism termed as Soret effect. On the contrary, when there is a heat transfer, and the difference causes an effect on the concentration gradient which is called the Dufour effect. Akram et al. [24] performed analysis in the existence of a tilted magnetic field as there are plenty of applications which are linked with peristaltic flows with double diffusion. Exact results can be acquired for a relatively new domain of nanoparticles, field of temperature, the concentration field, pressure rise, and the pressure gradient while considering the extent restrictions of longwave and a reduced Reynolds number. Authors in [25] have investigated the peristaltic flows that are induced by double diffusive convection naturally. This was to acquire analysis in an asymmetric porous channel. Current available literature studies on peristaltic flows with double diffused convection are available in references [26–29].The no-slip conditions are used in all of the abovementioned references. In some cases, there is surely a chance that a partial slip may happen between the fluids and the surface in motion. There are recorded instances when the fluids are particulate like there is emulsion or suspension [30]. In all these instances, boundary conditions are well explained by Navier condition [31], in which the proportionality is held between the relative slip and the local stress. Research on the effects on such slips on the peristaltic mobility has been stated in the literature [32–36].Considering the significance of the available and discussed literature and the applications linked with mixed convection in these peristaltic ducts, researchers have transformed their approaches and attention towards the theory related to double diffusion that have Carreau fluids as a fundamental fluid having slip boundaries. There is a strong chance that this investigation will head in the best manner for using the acquired data most efficiently. All the equations related to momentum, energies, and the concentration were required to be modelled with the suited physical settings such as lower wavenumbers and the Reynolds number. Equations have been solved numerically. All the results are depicted using graphs and further discussed.


## 2. Mathematical Formulation

We are in the process of investigating the incompressible trapping peristaltic flow of Carreau nanofluids in the 2D channel havingd˜1+d˜2 width. The coordinates are selected in such a way that the channel’s center is placed along the horizontal axis while vertical axis keeps the cross-sectional area. On the wall end, the speed of sinusoidal wave train stays persistent and does not change over time. Temperature, concentration of the solvents, and the nanoparticles of lower and the upper walls have also been considered as T1,T0, C1,C0, and Θ1,Θ0 respectively. The channel is inclined at an angle ζ. A persistent magnetic field is tested with angle β. In all this, the electrical field is taken as minimal as nil while Reynolds number is chosen as a low value that results in minimized or almost negligible induced magnetic field while comparing with the applied one. The mathematical expression for the surface of a wall is as follows [5]:(1)Y=H1=d˜1+a˜1cos2πλX−ct,Y=H2=−d˜2−b˜1cos2πλX−ct+φ,where λ, t, c, a˜1,b˜1, and d˜1,d˜2 are wavelength, time, speed of wave, wave amplitudes, and channel thickness, respectively. The range of phase difference φ is 0≤φ≤π, with φ=0 representing a symmetric channel without a phase wave and φ=π representing a channel with a phase wave. Furthermore, a˜1, d˜1, b˜1, d˜2, and φ meet the following condition: a˜12+b˜12+2a˜1b˜1cosφ≤d˜1+d˜22. The field of velocity for 2D channel flow is V=UX,Y,t,VX,Y,t,0.The expression of the Carreau nanofluid stress tensor is as follows [13]:(2)S=−μ∞+μ0−μ∞1+Γξ˙2n−1/2A˜,where Γ, μ∞, and μ0 denotes material constant, infinite viscosity, and zero shear rate viscosity, respectively. Also in the case where μ∞=0, the fundamental relation of S is(3)S=−μ01+Γξ˙2n−1/2A˜,The second invariant tensorξ˙ and first Rivlin–Ericksen tensor A˜ are defined as(4)A˜=∇V+∇V∗,ξ˙=12tracA˜2.ForΓ=0, the Carreau nanofluids are reduced to Newtonian nanofluids.In a fixed frame, the continuity, momentum, temperature, solute concentration, and nanoparticle fraction are defined [33] as(5)∂U∂X+∂V∂Y=0,(6)ρ∂∂t+U∂∂X+V∂∂YU=−∂p∂X−∂SXX∂X−∂SXY∂Y−σB02cosβUcosβ−Vsinβ+ρgsinζ+g1−Θ0ρf0βTT−T0+βCC−C0−ρp−ρf0Θ−Θ0,(7)ρ∂∂t+U∂∂X+V∂∂YV=−∂p∂Y−∂SYX∂X−∂SYY∂Y−ρgcosζ+σB02sinβUcosβ−Vsinβ,(8)ρcf∂T∂t+U∂T∂X+V∂T∂Y=k∂2T∂X2+∂2T∂Y2+ρcpDB∂Θ∂X∂T∂X+∂Θ∂Y∂T∂YDTT0∂T∂X2+∂T∂Y2+DTC∂2C∂X2+∂2C∂Y2,(9)∂C∂t+U∂C∂X+V∂C∂Y=Ds∂2C∂X2+∂2C∂Y2+DTC∂2T∂X2+∂2T∂Y2,(10)∂Θ∂t+U∂Θ∂X+V∂Θ∂Y=DB∂2Θ∂X2+∂2Θ∂Y2+DTT0∂2T∂X2+∂2T∂Y2,where the Carreau fluid stress tensor in component form is written as(11)SXX=−2μ01+n−12Γξ˙2∂U∂X,SXY=−μ01+n−12Γξ˙2∂U∂Y+∂V∂X,SYY=−2μ01+n−12Γξ˙2∂V∂Y,ξ˙2=2∂U∂X2+∂U∂Y+∂V∂X2+2∂V∂Y2.In the laboratory frameX,Y, the flow is unsteady, although the motion is steady in the coordinate system x,y. The transformation from one frame of reference to the other is described by(12)u=U−c,v=V,x=X−ct,y=Y,px,y=PX,Y,t.Define(13)x¯=xλ,y¯=yd˜1,v¯=vc,u¯=uc,δ=d˜1λ,d=d˜2d˜1,t¯=ctλ,h2=H2d˜2,h1=H1d˜1,b=b˜1d˜1,p¯=d˜12pμ0cλ,a=a˜1d˜1,Pr=ρcfυk,Re=ρfcd˜1μ0,υ=μ0ρf,Le=υDs,Ln=υDB,u=∂ψ∂y,v=−δ∂ψ∂x,θ=T−T0T1−T0,γ=C−C0C1−C0,Ω=Θ−Θ0Θ1−Θ0,We=cΓd˜1,S¯xx=Sxxλμ0c,S¯xy=Sxyd˜1μ0c,S¯yy=Syyd˜1μ0c,NCT=DCTT1−T0C1−C0Ds,NTC=DCTC1−C0kT1−T0,Grt=gd˜121−Θ0T1−T0ρfβTμ0c,Grc=g1−Θ0ρfβcC1−C0d˜12μ0c,GrF=gρp−ρfΘ1−Θ0μ0cd˜12,Nt=ρcpDTT1−T0T0k,Nb=ρcpDBΘ1−Θ0k,Fr=c2gd˜1,M=σμ0B0d˜1.Using equations (12) and (13), the equations (5)–(11) become(14)Reδψyψxy−ψxψyy=−∂p∂x−δ∂Sxx∂x−∂Sxy∂y−M2cosβψy+1cosβ+ψxδsinβ+ReFrsinζ+Grtθ+Grcγ−GrFΩ,(15)Reδ3ψxψxy−ψyψxx=−∂p∂y−δ2∂Sxy∂x−δ∂Syy∂y−δReFrcosζ+M2δsinβψy+1cosβ+ψxδsinβ,(16)Re Prδψyθx−ψxθy=θyy+δ2θxx+NTCδ2γxx+γyy+Nbδ2θxΩx+θyΩy+Ntδ2θx2+θy2,(17)ReδLeψyγx−ψxγy=δ2γxx+γyy+NCTδ2θxx+θyy,(18)ReδLnψyΩx−ψxΩy=δ2Ωxx+Ωyy+NtNbδ2θxx+θyy.The components of the stresses are represented in nondimensional form as(19)Sxx=−21+n−12We2ξ˙2∂2ψ∂x∂y,(20)Sxy=−1+n−12We2ξ˙2∂2ψ∂y2−δ2∂2ψ∂x2,(21)Syy=2δ1+n−12We2ξ˙2∂2ψ∂x∂y,(22)ξ˙2=2δ2∂2ψ∂x∂y2+∂2ψ∂y2−δ2∂2ψ∂x22+2δ2∂2ψ∂x∂y2.where ψ, Pr, Re, γ, p, θ, δ, Ω, We, Le, Nb, M, Nt, GrF, Grt, Grc, NCT, NTC, and Ln denotes stream function, Prandtl number, Reynolds number, solutal (species) concentration in dimensionless form, pressure, dimensionless temperature, wave number, nanoparticle fraction, Weissenberg number, Lewis number, Brownian motion, Hartmann number, thermophoresis parameter, Grashof number of nanoparticle, thermal Grashof number, solutal Grashof number, Soret parameter, Dufour parameter, and nanofluid Lewis number, respectively.Using estimations of long wavelengthδ≪1 and low Reynolds number, equations (14)–(22) are reduced to the form as(23)∂p∂x=∂∂y1+n−12We2∂2ψ∂y22∂2ψ∂y2−M2cos2βΨy+1+ReFrsinζ+Grtθ+Grcγ−GrFΩ,(24)−∂p∂y=0,(25)∂2θ∂y2+NTC∂2γ∂y2+Nb∂θ∂y∂Ω∂y+Nt∂θ∂y2=0,(26)∂2γ∂y2+NCT∂2θ∂y2=0,(27)∂2Ω∂y2+NtNb∂2θ∂y2=0.After removing pressure from equations (23) and (24), the governing equation for stream function is(28)∂2∂y21+n−12We2∂2ψ∂y22∂2ψ∂y2−M2cos2β∂2ψ∂y2+Grt∂θ∂y+Grc∂γ∂y−GrF∂Ω∂y=0.The mean flowQ is calculated in dimensionless form as follows:(29)Q=1+F+d,where(30)F=∫h2xh1x∂ψ∂ydy=ψh1x−ψh2x,here(31)h1x=1+acos2πx,h2x=−d−bcos2πx+φ.The foregoing system of PDE’s (25)–(28) will be solved under the following nondimensional boundary conditions:(32)ψ=F2,∂Ψ∂y=−η1Sxy−1ony=h1x,ψ=−F2,∂Ψ∂y=η1Sxy−1ony=h2x,θ+η2∂θ∂y=0,ony=h1x,θ−η2∂θ∂y=1,ony=h2x,γ+η3∂γ∂y=0,ony=h1x,γ−η3∂γ∂y=1,ony=h2x,Ω+η4∂Ω∂y=0,ony=h1x,Ω−η4∂Ω∂y=1,ony=h2x.There are no-slip conditions ifη1, η2, η3, η4=0 in the above conditions.


## 3. Results and Discussion

This section depicts the mathematical and realistic consequences of the present topic under study. The systems of PDE’s (23), (25)–(27), and (28) are coupled and nonlinear in nature. So, finding the exact solution to these equations is extremely difficult. As a result, the built-in command MATHEMATICA is used to calculate the numeric solution (NDSolve). To check the accuracy of the proposed numerical results and to compare our findings to the available literature, a slew of graphs was created.


### 3.1. Flow Characteristics

The profile of velocity under the stimulation of Hartmann numberM, Brownian motion Nb, Weissenberg number We, and velocity slip parameter η1 is shown in Figures 1(a)–1(d). The behaviour of velocity field for increasing values of Hartmann number is plotted in Figure 1(a). We can notice that increasing M tends to resist the flow, resulting in a decrease in fluid velocity. This occurs because when a magnetic field is produced, a resistive force (or “Lorentz force”) is created, which operates in the opposite direction of the flow, causing the velocity field to decrease (see Figure 1(a)). It is noted in Figure 1(b)) that the magnitude of velocity distribution decreases at the left side of the channel, whereas it decreases at the right side of the channel by raising the Brownian motion parameter Nb. It is depicted in Figure 1(c) that when Weissenberg number grows, velocity will increase. However, near the channel walls behaviour is opposite. This occurs because We has a direct relationship with the time constant and an inverse relationship with the channel width. The axial velocity increases near the channel’s center as the velocity slip parameter η1 is raised, as seen in Figure 1(d). Furthermore, opposite behavior is found near the boundaries.Figure 1
Impact of velocity onM, Nb, We, and η1.
(a)(b)(c)(d)


### 3.2. Pressure Rise and Pressure Gradient

Pumping is a hallmark aspect of peristalsis. The effects of several emergent factors on pressure rise per wavelength and pressure gradient are described in this section. Figures2(a)–2(d) show the change in pressure rise against flow rate Q for various values of the parameters n (power law index), We (Weissenberg number),M (Hartmann number), and η1 (velocity slip parameter). From these figures, it is noted that in retrograde Δp>0,Q<0, peristaltic Δp>0,Q>0, and free Δp=0 pumping regions, the pressure rise increases due to the increasing values of n, We, M, and η1 and whereas in copumping Δp<0,Q>0 region pressure rise drops due to the increasing values of n, We, M, and η1. The fluctuation of axial pressure gradient for different values of the η1 (velocity slip parameter), Nt (thermophoresis parameter), ζ (inclination angle), and We (Weissenberg number) is shown in Figures 3(a)–3(d). We observe that the pressure gradient is quite modest in the larger region of the channel at x∈0,0.3 and x∈0.5,1.0, indicating that the flow can readily pass without imposing a high-pressure gradient, where a substantially bigger pressure gradient is necessary to maintain the same flux in a small region of the channel at x∈0.3,0.5, particularly for the widest region at x=0.5. This is in line with the current physical situation. We also notice that as the velocity slip parameter η1 is increased, the axial pressure gradient falls (see Figure 3(a)). When Nt and ζ is increased, the pressure gradient increases (see Figures 3(b) and 3(c)). Furthermore, as seen in Figure 3(d), the magnitude value of axial pressure gradient grows as the Weissenberg number (We) increases.Figure 2
Impact of pressure rise onn, We, M, andη1.
(a)(b)(c)(d)Figure 3
Impact of pressure gradient onη1, Nt, ζ, and We.
(a)(b)(c)(d)


### 3.3. Temperature, Concentration, and Nanoparticle Volume Fraction

This section analyses the impact of many salient parameters on temperature, concentration, and nanoparticle volume fraction. Figures4(a)–4(c) shows the temperature profile for various values of thermal slip η2 parameter, Brownian motion Nb parameter, and Soret parameter NCT. It is noted in Figure 4(a) that when y∈−0.1,0.2, the temperature profile drops by increasing thermal slip parameter η2 whereas when y∈0.2,1 the temperature profile increases. The temperature increases by increasing Brownian motion parameter and Soret parameter (see Figures 4(b) and 4(c)). It is because Soret parameters have a direct relationship with temperature. Furthermore, nanoparticles travel efficiently from the wall to the fluid via enhancing Brownian motion. As a result of this impact, temperature rises as Nb and NCT rise. The concentration profile is shown in Figures 5(a)–5(c) for various values of the concentration slip parameter η3, Soret parameter NCT, and thermophoresis parameter Nt. It is noted in Figures 5(a) that increasing concentration slip parameter η3 lowers the concentration profile in the region y∈−0.1,0.6, while raising it in the region y∈0.6,1. The concentration decreases by increasing Soret parameter NCT and thermophoresis parameter Nt (see Figures 5(b) and 5(c)). It is due to the influence of random motion with micromixing and erratic collision activity of solid nanoparticles, which disperses the solid nanoparticles and so lowers the solute concentration. Figures 6(a)–6(c) show the nanoparticle fraction for various values of the nanoparticle slip parameter η4, Dufour parameter NTC, and Brownian motion parameter Nb. It is noticed in Figure 6(a) that increasing nanoparticle slip parameter η4 lowers the nanoparticle profile in the region where y∈−0.1,0.6, while raising it in the region where y∈0.6,1. When the Dufour parameter NTC is increased, the nanoparticle fraction falls (see Figure 6(b)). It is observed in Figure 6(c) that the nanoparticle volume fraction profile grows as Brownian motion parameter increases.Figure 4
Impact of temperature onη2, Nb, andNCT.
(a)(b)(c)Figure 5
Impact of concentration onη2, NCT, and Nt.
(a)(b)(c)Figure 6
Impact of nanoparticle onη4, NTC, and Nb.
(a)(b)(c)


### 3.4. Trapping

Another remarkable characteristic of peristaltic movement is trapping. It is the phenomenon of closed streamlines generating a bolus of fluid that circulates inside. Streamlines grab the fluid mass bolus and drive it forward with peristaltic waves at high flow rates and substantial occlusions. In Figures7–10, streamlines are plotted to show the phenomenon of trapping. The behaviour of streamlines for various Weissenberg numbers is depicted in Figure 7. It is observed from these streamlines that, by raising Weissenberg numbers, the number of trapped bolus grows in the upper half of the channel while the size of trapped bolus increases in the lower half. Figure 8 depicts the streamlines when the velocity slip parameter β1 is increased. As illustrated in Figure 8, the size of trapped bolus reduces in the upper half of the channel, whereas the number and size of trapped bolus decrease in the lower half by increasing β1 values. As the thermal Grashof number rises, the size of the trapped bolus grows larger (see Figure 9), and as the Hartmann number rises, the size of the trapped bolus shrinks (see Figure 10). The current study’s findings are compared to the existing literature in Table 1.Figure 7
Impact of streamlines onWe. (a) We=0.02 and (b) We=0.04.Figure 8
Impact of streamlines onβ1. (a) β1=0.01 and (b) β1=0.05.Figure 9
Impact of streamlines onGrt. (a) Grt=1.5 and (b) Grt=2.0.Figure 10
Impact of streamlines onM. (a) M = 1.5 and (b) M = 2.0.Table 1
The comparison of velocity profile with existing literature.
yPresent workMishra and Rao [5]0.919093−1−10.796856−0.5996340.5302070.674618−0.243748−0.188060.552380.04937340.04344740.4301420.2616440.1785430.3079050.349980.2251190.1856670.2651350.1852510.0634290.0540730.0553249−0.0588087−0.239407−0.174226−0.181046−0.596941−0.519578−0.303284−1−1


## 3.1. Flow Characteristics

The profile of velocity under the stimulation of Hartmann numberM, Brownian motion Nb, Weissenberg number We, and velocity slip parameter η1 is shown in Figures 1(a)–1(d). The behaviour of velocity field for increasing values of Hartmann number is plotted in Figure 1(a). We can notice that increasing M tends to resist the flow, resulting in a decrease in fluid velocity. This occurs because when a magnetic field is produced, a resistive force (or “Lorentz force”) is created, which operates in the opposite direction of the flow, causing the velocity field to decrease (see Figure 1(a)). It is noted in Figure 1(b)) that the magnitude of velocity distribution decreases at the left side of the channel, whereas it decreases at the right side of the channel by raising the Brownian motion parameter Nb. It is depicted in Figure 1(c) that when Weissenberg number grows, velocity will increase. However, near the channel walls behaviour is opposite. This occurs because We has a direct relationship with the time constant and an inverse relationship with the channel width. The axial velocity increases near the channel’s center as the velocity slip parameter η1 is raised, as seen in Figure 1(d). Furthermore, opposite behavior is found near the boundaries.Figure 1
Impact of velocity onM, Nb, We, and η1.
(a)(b)(c)(d)


## 3.2. Pressure Rise and Pressure Gradient

Pumping is a hallmark aspect of peristalsis. The effects of several emergent factors on pressure rise per wavelength and pressure gradient are described in this section. Figures2(a)–2(d) show the change in pressure rise against flow rate Q for various values of the parameters n (power law index), We (Weissenberg number),M (Hartmann number), and η1 (velocity slip parameter). From these figures, it is noted that in retrograde Δp>0,Q<0, peristaltic Δp>0,Q>0, and free Δp=0 pumping regions, the pressure rise increases due to the increasing values of n, We, M, and η1 and whereas in copumping Δp<0,Q>0 region pressure rise drops due to the increasing values of n, We, M, and η1. The fluctuation of axial pressure gradient for different values of the η1 (velocity slip parameter), Nt (thermophoresis parameter), ζ (inclination angle), and We (Weissenberg number) is shown in Figures 3(a)–3(d). We observe that the pressure gradient is quite modest in the larger region of the channel at x∈0,0.3 and x∈0.5,1.0, indicating that the flow can readily pass without imposing a high-pressure gradient, where a substantially bigger pressure gradient is necessary to maintain the same flux in a small region of the channel at x∈0.3,0.5, particularly for the widest region at x=0.5. This is in line with the current physical situation. We also notice that as the velocity slip parameter η1 is increased, the axial pressure gradient falls (see Figure 3(a)). When Nt and ζ is increased, the pressure gradient increases (see Figures 3(b) and 3(c)). Furthermore, as seen in Figure 3(d), the magnitude value of axial pressure gradient grows as the Weissenberg number (We) increases.Figure 2
Impact of pressure rise onn, We, M, andη1.
(a)(b)(c)(d)Figure 3
Impact of pressure gradient onη1, Nt, ζ, and We.
(a)(b)(c)(d)


## 3.3. Temperature, Concentration, and Nanoparticle Volume Fraction

This section analyses the impact of many salient parameters on temperature, concentration, and nanoparticle volume fraction. Figures4(a)–4(c) shows the temperature profile for various values of thermal slip η2 parameter, Brownian motion Nb parameter, and Soret parameter NCT. It is noted in Figure 4(a) that when y∈−0.1,0.2, the temperature profile drops by increasing thermal slip parameter η2 whereas when y∈0.2,1 the temperature profile increases. The temperature increases by increasing Brownian motion parameter and Soret parameter (see Figures 4(b) and 4(c)). It is because Soret parameters have a direct relationship with temperature. Furthermore, nanoparticles travel efficiently from the wall to the fluid via enhancing Brownian motion. As a result of this impact, temperature rises as Nb and NCT rise. The concentration profile is shown in Figures 5(a)–5(c) for various values of the concentration slip parameter η3, Soret parameter NCT, and thermophoresis parameter Nt. It is noted in Figures 5(a) that increasing concentration slip parameter η3 lowers the concentration profile in the region y∈−0.1,0.6, while raising it in the region y∈0.6,1. The concentration decreases by increasing Soret parameter NCT and thermophoresis parameter Nt (see Figures 5(b) and 5(c)). It is due to the influence of random motion with micromixing and erratic collision activity of solid nanoparticles, which disperses the solid nanoparticles and so lowers the solute concentration. Figures 6(a)–6(c) show the nanoparticle fraction for various values of the nanoparticle slip parameter η4, Dufour parameter NTC, and Brownian motion parameter Nb. It is noticed in Figure 6(a) that increasing nanoparticle slip parameter η4 lowers the nanoparticle profile in the region where y∈−0.1,0.6, while raising it in the region where y∈0.6,1. When the Dufour parameter NTC is increased, the nanoparticle fraction falls (see Figure 6(b)). It is observed in Figure 6(c) that the nanoparticle volume fraction profile grows as Brownian motion parameter increases.Figure 4
Impact of temperature onη2, Nb, andNCT.
(a)(b)(c)Figure 5
Impact of concentration onη2, NCT, and Nt.
(a)(b)(c)Figure 6
Impact of nanoparticle onη4, NTC, and Nb.
(a)(b)(c)


## 3.4. Trapping

Another remarkable characteristic of peristaltic movement is trapping. It is the phenomenon of closed streamlines generating a bolus of fluid that circulates inside. Streamlines grab the fluid mass bolus and drive it forward with peristaltic waves at high flow rates and substantial occlusions. In Figures7–10, streamlines are plotted to show the phenomenon of trapping. The behaviour of streamlines for various Weissenberg numbers is depicted in Figure 7. It is observed from these streamlines that, by raising Weissenberg numbers, the number of trapped bolus grows in the upper half of the channel while the size of trapped bolus increases in the lower half. Figure 8 depicts the streamlines when the velocity slip parameter β1 is increased. As illustrated in Figure 8, the size of trapped bolus reduces in the upper half of the channel, whereas the number and size of trapped bolus decrease in the lower half by increasing β1 values. As the thermal Grashof number rises, the size of the trapped bolus grows larger (see Figure 9), and as the Hartmann number rises, the size of the trapped bolus shrinks (see Figure 10). The current study’s findings are compared to the existing literature in Table 1.Figure 7
Impact of streamlines onWe. (a) We=0.02 and (b) We=0.04.Figure 8
Impact of streamlines onβ1. (a) β1=0.01 and (b) β1=0.05.Figure 9
Impact of streamlines onGrt. (a) Grt=1.5 and (b) Grt=2.0.Figure 10
Impact of streamlines onM. (a) M = 1.5 and (b) M = 2.0.Table 1
The comparison of velocity profile with existing literature.
yPresent workMishra and Rao [5]0.919093−1−10.796856−0.5996340.5302070.674618−0.243748−0.188060.552380.04937340.04344740.4301420.2616440.1785430.3079050.349980.2251190.1856670.2651350.1852510.0634290.0540730.0553249−0.0588087−0.239407−0.174226−0.181046−0.596941−0.519578−0.303284−1−1


## 4. Conclusions

The importance of partial slip on double diffusive convection on magneto-Carreau nanofluid in inclined peristaltic asymmetric channel is discussed. The problem is numerically solved. The following is a summary of the main findings from the above analysis:The velocity drops as the Hartmann number rises, whereas the velocity rises as the Weissenberg number risesAs the velocity slip parameterη1 is increased, the axial velocity at the center of channel increasesThe pressure gradient tends to grow as the thermophoresis parameter, inclination angle, and Weissenberg number increaseThe temperature rises as the Brownian motion parameter and the Soret parameter are increasedBy increasing the Soret and thermophoresis parameters, the concentration dropsThe nanoparticle fraction grows when the Brownian motion parameter is raised, whereas the nanoparticle fraction decreases as the Dufour parameter is increasedAs the thermal Grashof number rises, the size of the trapped bolus grows larger, and as the Hartmann number rises, the size of the trapped bolus shrinks



---

*Source: 2475846-2021-11-24.xml*