Impact of Partial Slip on Double Diffusion Convection of Sisko Nanofluids in Asymmetric Channel with Peristaltic Propulsion and Inclined Magnetic Field

The current article discusses the outcomes of the double diffusion convection of peristaltic transport in Sisko nanofluids along an asymmetric channel having an inclined magnetic field. Consideration is given to the Sisko fluid model, which can forecast both Newtonian and non-Newtonian fluid properties. Lubricating greases are the best examples of Sisko fluids. Experimental research shows that most realistic fluids, including human blood, paint, dirt, and other substances, correspond to Sisko’s proposed definition of viscosity. Mathematical modelling is considered to explain the flow behavior. The simpler non-linear PEDs are deduced by using an elongated wavelength and a minimal Reynolds number. The expression is also numerically calculated. The impacts of the physical variables on the quantities of flow are plotted graphically as well as numerically. The results reveal that there is a remarkable increase in the concentration, temperature, and nanoparticle fraction with the rise in the Dufour and thermophoresis variables.


Introduction
The study of non-Newtonian models is the pinnacle of contemporary research owing to the fact that it has tremendous practical utility in the manufacturing field. The common domains are electrical engineering or electricity generation, the energy industry, heatpower engineering, and the chemical and pharmaceutical industry. Moreover, many studies have focused on the application of non-Newtonian fluids in the biomedical industry, environmental management, and catalysis science and technology. The classification of non-Newtonian fluids is quite broad and one can find various categories of models because exclusive intrinsic relation to any category is insufficient to comprehensively elaborate all the properties of these fluids. The Sisko model is notable among all models because the constitutive equation of the model is exhaustive enough to reveal the fluid's behavior. It can elaborate the shear-thinning trend in the fluid as the shear-thinning rate is low and velocity reaches the constant value when the shear rate is high. The same phenomenon is used to study the blood circulation mechanism. Sisko, in 1958 [1], developed the model with the help of the power law model of fluid that is the combination of non-Newtonian and Newtonian fluids. Naturally occurring fluids fall in this category that has many practical applications such as the transport of greasy fluids.
Peristaltic techniques and their reactions have fundamental implications in human physiology, or one can say that they are the essential mechanisms on which all human systems function. The term peristalsis is coined from the Greek word peristaltikos meaning condensing and contracting. All the muscular movements in the human body are peristaltic transfer are some of the other significant fields of application. The layout of this paper is as follows: The motion equation along with the Sisko fluid equation is discussed in Section 2. Formulation of the problem is in Section 3, which also reflects on the governing problem in a linear fashion. Section 4 covers the numerical and graphical results and its explanation, while the last section is a comprehensive conclusion based on the main findings.

Mathematical Formulation
Let us consider peristaltic flux in non-Newtonian fluid in an asymmetric conduit with continual magnetic field. The assumptions are that the wave line on channel walls is moving with velocity c, and the angle of magnetic field is inclined at ω. A uniform magnetic field B = (B 0 sinω, B 0 cosω, 0) is imposed. When considering minute magnetic Reynolds number, the induced magnetic field is neglected. Moreover, the y = H 1 (top wall) and y = H 2 (lower wall) are held at temperature, concentration, and solute of nanoparticles with T = T 0 , C = C 0 , Θ = Θ 0 , and T = T 1 , C = C 1 , Θ = Θ 1 , respectively.
The structure of a wall is sketched in Figure 1, and mathematically it is defined as [19] Y = H 1 = a 1 + a 2 cos (X − ct) 2π λ , any normal cell, they are used to carry drug particles to the tumor-infected region through the blood stream. They have more power to absorb heat than any other microparticle; hence, they are used to change the magnetic current field within the human body that is required in cancer therapy. Contrast magnetic resonance imaging, hyperthermia, therapeutic procedures, and drug transfer are some of the other significant fields of application. The layout of this paper is as follows: The motion equation along with the Sisko fluid equation is discussed in Section 2. Formulation of the problem is in Section 3, which also reflects on the governing problem in a linear fashion. Section 4 covers the numerical and graphical results and its explanation, while the last section is a comprehensive conclusion based on the main findings.

Mathematical Formulation
Let us consider peristaltic flux in non-Newtonian fluid in an asymmetric conduit with continual magnetic field. The assumptions are that the wave line on channel walls is moving with velocity , and the angle of magnetic field is inclined at . A uniform magnetic field = ( , , 0) is imposed. When considering minute magnetic Reynolds number, the induced magnetic field is neglected. Moreover, the = (top wall) and = (lower wall) are held at temperature, concentration, and solute of nanoparticles with = , = , = , and = , = , = , respectively.
The structure of a wall is sketched in Figure 1, and mathematically it is defined as [19] = =̄ +̄ ( − ) 2 , Here stands for wave propagation direction, represents velocity propagation, (̄ , ̄ ) act as wave amplitudes, ̄ +̄ is channel width, denotes wavelength, Here X stands for wave propagation direction, c represents velocity propagation, (a 2 , a 4 ) act as wave amplitudes, a 1 + a 3 is channel width, λ denotes wavelength, t represents time, The phase difference φ range is 0 ≤ φ ≤ π, φ = 0 implies a symmetric channel, without waves phase, and φ = π indicates a channel with in waves phase. In addition, a 1 , a 2 , a 3 , a 4 , and φ satisfy the constraint a 2 2 + a 2 4 + 2a 2 a 4 cosφ ≤ (a 1 + a 3 ) 2 . The Sisko fluid stress tensor is defined by [1]  The velocities U and V along X and Y are taken for the current flow, so the continuity equation, momentum, temperature, nanoparticles fraction, and the solute concentration of an incompressible fluid at a fixed frame, is given as ∂U ∂X where component form of stresses is defined as where g, k, T, d dt , ρ f , ρ p , ρ f 0 , p, (ρc) p , (ρc) f , C, β C , Θ, β T , D s , D TC , D B , D CT , D T represent acceleration, thermal conductivity, temperature, material derivative, fluid density, nanoparticle mass density, fluid density at T 0 , pressure, nanoparticle heat capacity, fluid heat capacity, solutal concentration, volumetric coefficient of solutal expansion, nanoparticle fraction, volumetric coefficient of thermal expansion, solutal diffusively, Dufour diffusively, Brownian diffusion, Soret diffusively, thermophoretic diffusion, respectively. Now, using the well-known Galilean transformation leads to Dimensionless quantities are defined as Here Ω, Le, Pr, γ, G rt , δ, N CT , G rF , N t , θ, N b , Ln, N TC , G rc denote nanoparticle fraction, Lewis number, Prandtl number, solutal (species) concentration, thermal Grashof number, wave number, Soret parameter, nanoparticle Grashof number, Reynolds number, thermophoresis parameters, temperature, Brownian motion, nanofluid Lewis number, parameter of Dufour, and solutal Grashof number, respectively.
Using Equations (10) and (11), Equation (3) is satisfied automatically and Equations (4)-(9) after bars dropping in wave frame become ReδPr ReLeδ where Now, using approximation of low Reynolds number and long wavelength, the Equations (12)- (17) are reduced to the form as where By eliminating the pressure from Equations (18) and (19), we obtained the following equations The existing problem boundary conditions are expressed in wave frame as follows The case of no slip boundary conditions exists when we consider Γ 1 , Γ 2 , Γ 3 , Γ 4 = 0 in Equations (27)- (30).
The flow rate Q in dimensionless form is represented as [19] here

Numerical Solution and Graphical Outcomes
Due to the non-linear system and coupled behavior, exact solutions to the systems of PDEs (18), (20)- (22), and (25) are difficult to find. As a result, we can use Mathematica's NDSolve command to find the numerical solutions to Equations (18), (20)- (22), and (25) as well as the boundary conditions (26)- (29). NDSolve is a Wolfram language function that solves numerical differential equations in a broad sense. In addition to various ordinary differential equations, it can handle a variety of partial differential equations. This command uses interpolating function objects to iteratively find the solutions. As a result of the numerical solution, the impact of the evolving parameters of fluid quantities is examined via graphical illustration.

Effects of Hartmann Number (M)
The plots in Figure 2a-c are drawn to see how the Hartmann number (M) affects the velocity, pressure gradient, and pressure rise. In the interval when y ∈ [−0.8, −0.1], the velocity magnitude increases as M becomes larger, while in the interval when y ∈ [−0.1, 0.1], the reverse trend is observed (see Figure 2a). This occurs because as the magnetic number grows, the Lorentz force, which acts as a retarding force, increases, resulting in the fluid motion decelerating. The characteristic of pressure rises on M is depicted in Figure 2b. It is obvious that escalating M enhances the pressure rise of the retrograde zone (∆p > 0, Q < 0), peristaltic zone (∆p > 0, Q > 0), and free (∆p = 0) pumping segment, while increasing M lessens the pressure rise in the augmented (Q > 0, ∆p < 0) region. As the Hartmann number rises, the pressure gradient magnitude goes up as well (see Figure 2c). . It is obvious that escalating enhances the pressure rise of the retrograde zone ( > 0, < 0), peristaltic zone ( > 0, > 0), and free ( = 0) pumping segment, while increasing lessens the pressure rise in the augmented ( > 0, < 0) region. As the Hartmann number rises, the pressure gradient magnitude goes up as well (see Figure 2c).

Effects of Slip Parameters ( , , , )
Figure 3a-j have been created to analyze the effect of slip factors on the velocity, temperature, solvent concentration, nanoparticle fraction, pressure rise, and pressure gradient. Figure 3a,b explain the effects of velocity on (slip parameter of velocity) and (slip parameter of nanoparticles). It is illustrated in Figure 3a that rising values increase the velocity magnitude at the channel's center, but the opposite behavior occurs near the channel walls. In Figure Figure 3f that raising the velocity slip parameter reduces the pressure rise in the peristaltic zone ( > 0, > 0), retrograde zone ( > 0, < 0), and free ( = 0) pumping segments, while increasing tends to increase the pressure rise in the augmented ( > 0, < 0) area. It is significant to note that when enhancing the slip parameters of concentration , the pressure rise in all peristaltic pumping zones (retrograde, peristaltic, free, and augmented pumping areas) tends to decrease (see Figure 3g). The effects of the slip factors of velocity, temperature, and nanoparticles ( , and ) on the pressure gradient are shown in Figure 3h,i. The pressure gradient falls at the channel's center when the slip factors of velocity are increased, but the opposite behavior is observed near the channel walls, as shown in Figure 3h. As shown in Figure 3i, raising the temperature slip factor lowers the pressure gradient. Contrary to this, the slip factor of nanoparticles has the opposite effect (see Figure 3j).   Figure 3f that raising the velocity slip parameter Γ 1 reduces the pressure rise in the peristaltic zone (∆p > 0, Q > 0), retrograde zone (∆p > 0, Q < 0), and free (∆p = 0) pumping segments, while increasing Γ 1 tends to increase the pressure rise in the augmented (Q > 0, ∆p < 0) area. It is significant to note that when enhancing the slip parameters of concentration Γ 3 , the pressure rise in all peristaltic pumping zones (retrograde, peristaltic, free, and augmented pumping areas) tends to decrease (see Figure 3g). The effects of the slip factors of velocity, temperature, and nanoparticles (Γ 1 , Γ 2 and Γ 4 ) on the pressure gradient are shown in Figure 3h,i. The pressure gradient falls at the channel's center when the slip factors of velocity are increased, but the opposite behavior is observed near the channel walls, as shown in Figure 3h. As shown in Figure 3i, raising the temperature slip factor Γ 2 lowers the pressure gradient. Contrary to this, the slip factor of nanoparticles Γ 4 has the opposite effect (see Figure 3j).

Effects of Brownian Motion ( )
The Brownian motion's, , impacts on the axial velocity, temperature, nanoparticle fraction, pressure rise, solvent concentration, and pressure gradient are assessed in Figure  4a-f. Figure 4a demonstrates the Brownian motion's impact on the velocity profile. As noticed in Figure 4a, the magnitude value of the velocity declines in the region where ∈ Figure 3. Slip parameters' (Γ 1 , Γ 2 , Γ 3 , Γ 4 ) impact on velocity, temperature, solutal concentration, nanoparticle fraction, pressure rise, and pressure gradient. Other parameters' values are as follows:

Effects of Brownian Motion (N b )
The Brownian motion's, N b , impacts on the axial velocity, temperature, nanoparticle fraction, pressure rise, solvent concentration, and pressure gradient are assessed in Figure 4a-f. Figure 4a demonstrates the Brownian motion's impact on the velocity profile. As noticed in Figure 4a, the magnitude value of the velocity declines in the region where y ∈ [−0.8, −0.2], i.e., opposes the backflow. As a result, the Brownian motion actually slows down the flow, whereas the situation is reversed when y ∈ [−0.2, 0.3]. The velocity magnitude rises in the given area as the Brownian motion parameter rises. It is worth noting in Figure 4b-d that increasing the Brownian motion parameter values lowers the temperature and nanoparticle fraction profiles (see Figure 4b,c), whereas increasing the Brownian motion parameter values raises the solvent concentration profile (see Figure 4d). Figure 4e depicts the role of Brownian motion on the pressure rise. The pressure rise drops in all peristaltic pumping zones as the values of the Brownian motion increase (see Figure 4e). The consequences of the pressure gradient on the Brownian motion are seen in Figure 4f. It should be noticed that as N b increases, the pressure gradient decreases, with the highest pressure gradient occurring at x = 0.5.  Figure 4b-d that increasing the Brownian motion parameter values lowers the temperature and nanoparticle fraction profiles (see Figure 4b,c), whereas increasing the Brownian motion parameter values raises the solvent concentration profile (see Figure 4d). Figure  4e depicts the role of Brownian motion on the pressure rise. The pressure rise drops in all peristaltic pumping zones as the values of the Brownian motion increase (see Figure 4e). The consequences of the pressure gradient on the Brownian motion are seen in Figure 4f.
It should be noticed that as increases, the pressure gradient decreases, with the highest pressure gradient occurring at = 0.5.

Effects of Thermophoresis Parameter (N t )
Figure 5a-f are shown to illustrate the effects of the thermophoresis (N t ) parameter on the axial velocity, temperature, concentration, pressure rise, nanoparticle fraction, and pressure gradient. Figure 5a depicts the velocity curve for N t values. The velocity magnitude grows when the fluid moves in the region y ∈ [−0.8, −0.2] due to the growing N t parameter, whereas the opposite effect is noted when the fluid moves in the region y ∈ [−0.2, 0.3]. Here, with the rise in the thermophoresis (N t ) parameter, the magnitude of the velocity drops. Thermophoresis is a valuable source, which is caused by a temperature difference between the hot gas and the cold objects. It also controls particle migration to the cold wall. It must be noted that the temperature distribution varies as the thermophoresis parameter varies. The plots in Figure 5b show that temperature drops as the thermophoresis parameter N t rises. However, in the case of the concentration and nanoparticle fraction, the opposite effect has been demonstrated. In this scenario the temperature rises as the thermophoresis parameter N t rises (see Figure 5c,d). Figure 5e reveals the visual aspect of the pressure rise on the thermophoresis (N t ) parameter. Increasing the thermophoresis parameter enhances the pressure rise in all peristaltic pumping zones, as seen in Figure 5e. The impact of the pressure gradient on thermophoresis (N t ) is depicted graphically in Figure 5f. As seen in Figure 5f, the pressure gradient widens as the thermophoresis parameter rises.
Here, with the rise in the thermophoresis ( ) parameter, the magnitude of the velocity drops. Thermophoresis is a valuable source, which is caused by a temperature difference between the hot gas and the cold objects. It also controls particle migration to the cold wall. It must be noted that the temperature distribution varies as the thermophoresis parameter varies. The plots in Figure 5b show that temperature drops as the thermophoresis parameter rises. However, in the case of the concentration and nanoparticle fraction, the opposite effect has been demonstrated. In this scenario the temperature rises as the thermophoresis parameter rises (see Figure 5c,d). Figure 5e reveals the visual aspect of the pressure rise on the thermophoresis ( ) parameter. Increasing the thermophoresis parameter enhances the pressure rise in all peristaltic pumping zones, as seen in Figure 5e. The impact of the pressure gradient on thermophoresis ( ) is depicted graphically in Figure 5f. As seen in Figure 5f, the pressure gradient widens as the thermophoresis parameter rises.
(a)  Figure 6a,b show the effects of the pressure rise on the nanoparticle and thermal Grashof numbers. Figure 6a,b show that increasing the number of nanoparticles Grashof ( ) and thermal Grashof ( ) increases the pressure rise in all peristaltic pumping zones. Figure 6c,d depict the roles of the nanoparticles and thermal Grashof numbers on the pressure rise. These graphs show that raising and leads to an increase in the pressure gradient. In Figure 6e, the role of the nanoparticle Grashof ( ) is investigated. Figure 6e illustrates that as increases, the velocity magnitude tends to rise in the region ∈ [−0.8, −0.1], whereas velocity shows the opposite tendency in the region ∈ [−0.1,0.3]. The influence of the thermal Grashof number ( ) on velocity is seen in Figure 6f. The marginal effect of the thermal buoyancy force and viscous hydrodynamic force is represented by this parameter. The viscous forces dominate the peristaltic regime for < 1, while the viscous forces dominate the peristaltic regime for

Effects of Nanoparticle Grashof (G rF ) and Thermal Grashof (G rt )
Figure 6a-f show the effects of the nanoparticle and thermal Grashof numbers on flow quantities. Figure 6a,b show the effects of the pressure rise on the nanoparticle and thermal Grashof numbers. Figure 6a,b show that increasing the number of nanoparticles Grashof (G rF ) and thermal Grashof (G rt ) increases the pressure rise in all peristaltic pumping zones. Figure 6c,d depict the roles of the nanoparticles and thermal Grashof numbers on the pressure rise. These graphs show that raising G rF and G rt leads to an increase in the pressure gradient. In Figure 6e, the role of the nanoparticle Grashof (G rF ) is investigated. Figure 6e illustrates that as G rF increases, the velocity magnitude tends to rise in the region y ∈ [−0.8, −0.1], whereas velocity shows the opposite tendency in the region y ∈ [−0.1, 0.3]. The influence of the thermal Grashof number (G rt ) on velocity is seen in Figure 6f. The marginal effect of the thermal buoyancy force and viscous hydrodynamic force is represented by this parameter. The viscous forces dominate the peristaltic regime for G rt < 1, while the viscous forces dominate the peristaltic regime for G rt > 1. As the thermal Grashof number rises, the velocity magnitude drops in the region y ∈ [−0.8, −0.1] and vice versa for y ∈ [−0.1, 0.3]. In most cases, thermal buoyancy tends to slow down the flow throughout the regime.

Trapping Effects
Trapping seems to be a remarkable phenomenon of propelled peristaltic flows. Trapping occurs when the inner fluid moves, which creates a mass along with a peristaltic wave streamline. Streamlines grab a fluid mass bolus to drive it forward with peristaltic waves at high flow rates and substantial occlusions. To investigate the phenomenon of Figure 6. Thermal Grashof (G rt ), nanoparticle Grashof (G rF ) impacts on pressure rise, pressure gradient, and velocity. Other parameters' values are as follows: N CT = 0.9, N TC = 1.4,

Trapping Effects
Trapping seems to be a remarkable phenomenon of propelled peristaltic flows. Trapping occurs when the inner fluid moves, which creates a mass along with a peristaltic wave streamline. Streamlines grab a fluid mass bolus to drive it forward with peristaltic waves at high flow rates and substantial occlusions. To investigate the phenomenon of streamlines, Figures 7-10 are drawn. Figure 7 shows a decrease in the trapped bolus volume of the lower part of the channel, whereas the amount and volume of the trapped bolus decreases in the upper half due to growing values of the slip parameter of velocity Γ 1 . In Figure 8 it is noted that due to growing values of the slip parameter of temperature Γ 2 , the quantity and size of the trapped bolus increases in the lower half of the channel, whereas the trapped bolus size decreases in the upper half. From Figure 9 it is illustrated that in the upper section of the channel, raising the solutal Grashof number G rc increases the number of trapped boluses, but the opposite impact is shown in the lower part. As seen in Figure 10, the number of boluses in the upper area of the channel decreases as the Dufour parameter is increased; however, the trapped bolus volume grows in the lower half.  Figure 7 shows a decrease in the trapped bolus volume of the lower part of the channel, whereas the amount and volume of the trapped bolus decreases in the upper half due to growing values of the slip parameter of velocity . In Figure 8 it is noted that due to growing values of the slip parameter of temperature , the quantity and size of the trapped bolus increases in the lower half of the channel, whereas the trapped bolus size decreases in the upper half. From Figure 9 it is illustrated that in the upper section of the channel, raising the solutal Grashof number increases the number of trapped boluses, but the opposite impact is shown in the lower part. As seen in Figure 10, the number of boluses in the upper area of the channel decreases as the Dufour parameter is increased; however, the trapped bolus volume grows in the lower half.   Figure 7 shows a decrease in the trapped bolus volume of the lower part of the channel, whereas the amount and volume of the trapped bolus decreases in the upper half due to growing values of the slip parameter of velocity . In Figure 8 it is noted that due to growing values of the slip parameter of temperature , the quantity and size of the trapped bolus increases in the lower half of the channel, whereas the trapped bolus size decreases in the upper half. From Figure 9 it is illustrated that in the upper section of the channel, raising the solutal Grashof number increases the number of trapped boluses, but the opposite impact is shown in the lower part. As seen in Figure 10, the number of boluses in the upper area of the channel decreases as the Dufour parameter is increased; however, the trapped bolus volume grows in the lower half.

Concluding Remarks
This section covers the concluding remarks on the ongoing problem. A theoretical evaluation is shown to explore the partial slip impact on the double diffusion convection of peristaltic transport in Sisko nanofluids along an asymmetric channel by taking an inclined magnetic field. To understand the flow dynamics of the current problem, mathematical modeling is considered. Numerical solutions are proposed for the problem under analysis. The key findings are as follows:


The temperature and nanoparticle profiles drop as the Brownian motion is increased, while the concentration profile rises.  The velocity magnitude at the channel's center grows as the slip factor of velocity ( ) increases, but the opposite behavior occurs at the channel walls.  As the slip parameters of concentration and nanoparticles are enhanced, the solvent concentration and nanoparticle fraction increase.  The temperature drops as the thermophoresis parameter rises.  The amount of bolus in the upper channel decreases as the Dufour parameter is increased; however, the volume of the trapped bolus grows in the bottom half.  The trapped bolus grows in terms of size and number in the bottom half of the channel as the slip parameter of temperature increases, while the size of the trapped bolus decreases in the top half.

Concluding Remarks
This section covers the concluding remarks on the ongoing problem. A theoretical evaluation is shown to explore the partial slip impact on the double diffusion convection of peristaltic transport in Sisko nanofluids along an asymmetric channel by taking an inclined magnetic field. To understand the flow dynamics of the current problem, mathematical modeling is considered. Numerical solutions are proposed for the problem under analysis. The key findings are as follows:


The temperature and nanoparticle profiles drop as the Brownian motion is increased, while the concentration profile rises.  The velocity magnitude at the channel's center grows as the slip factor of velocity ( ) increases, but the opposite behavior occurs at the channel walls.  As the slip parameters of concentration and nanoparticles are enhanced, the solvent concentration and nanoparticle fraction increase.  The temperature drops as the thermophoresis parameter rises.  The amount of bolus in the upper channel decreases as the Dufour parameter is increased; however, the volume of the trapped bolus grows in the bottom half.  The trapped bolus grows in terms of size and number in the bottom half of the channel as the slip parameter of temperature increases, while the size of the trapped bolus decreases in the top half.

Concluding Remarks
This section covers the concluding remarks on the ongoing problem. A theoretical evaluation is shown to explore the partial slip impact on the double diffusion convection of peristaltic transport in Sisko nanofluids along an asymmetric channel by taking an inclined magnetic field. To understand the flow dynamics of the current problem, mathematical modeling is considered. Numerical solutions are proposed for the problem under analysis.
The key findings are as follows:

•
The temperature and nanoparticle profiles drop as the Brownian motion is increased, while the concentration profile rises.

•
The velocity magnitude at the channel's center grows as the slip factor of velocity (Γ 1 ) increases, but the opposite behavior occurs at the channel walls.
• As the slip parameters of concentration and nanoparticles are enhanced, the solvent concentration and nanoparticle fraction increase.

•
The temperature drops as the thermophoresis parameter rises.

•
The amount of bolus in the upper channel decreases as the Dufour parameter is increased; however, the volume of the trapped bolus grows in the bottom half.