Abstract

The double-diffusive convection is a significant physical phenomenon that arises in fluid mechanics. It is primarily associated with a convection process in which two dissimilar density gradients with varying diffusion rates are considered. The primary goal of this study is to investigate the effects of double-diffusivity convection and partial slip with an inclined magnetic field on peristaltic propulsion in an asymmetric channel for Oldroyd-4 constants nanofluids. The flow of an Oldroyd-4 constant nanofluid is mathematically modeled in the presence of double-diffusivity convection and a tilted magnetic field. Lubrication methodology is applied to simplify the highly nonlinear system of partial differential equations (PDEs). The numerical scheme is used to calculate the solution of coupled nonlinear PDEs. Furthermore, the effect of changing the parameters associated with slip, thermophoresis, Brownian motion, Grashof number of nanoparticles, Hartmann number, pumping, and trapping are investigated in this article. It is noticed that the temperature rises as the Brownian motion and thermophoresis constraints increases. This is because the growth in the Brownian motion parameter indicates the increase in the kinetic energy of nanoparticles which results in warming up the nanofluid. Also, concentration falls as the Brownian motion and thermophoresis constraints increases.

1. Introduction

The peristaltic transportation of biological fluids has huge number of applications in biomedical industry and hence received huge attention from the researchers in last two decades. Transport phenomenon of biological fluids is based on peristaltic propulsion. This phenomenon can be internally observed in the gastrointestinal tract, urine flow, blood flow, male reproductive tract, esophageal swallowing, and ureter and externally as worm's movement. This natural process is caused by periodic relaxation and contraction of muscles which produce a sinusoidal wave along the channel or tube walls. Some of the worthwhile theoretical and experimental studies considering the peristaltic propulsion for numerous basic fluid models are cited in references [16]. These studies have further been generalized for many non-Newtonian fluid models like the Carreau model [7], Maxwell model [8], Johnson–Segalman model [9], Williamson model [10], Casson model [11], six-constant Jeffreys model [12], Walter's B model [13], hyperbolic tangent model [14], Herschel–Bulkley model [15], Jeffrey model [16], Sisko fluid [17], Oldroyd-4 model [18], Phan–Thien–Tanner model [19], and Burgers' model [20]. The authors of these research studies extended the work to new findings which may widen the horizon of their usage to biomedical engineering and technology.

In the field of biomagnetic fluid dynamics (BFD), the study of peristaltic flow is based on magnetohydrodynamic (MHD) effects. The study of these biological fluids (examples of biofluids include blood, urine, and chyme) flows is important in bioengineering and medical sciences. These fluids are broadly found in living organisms and influence of magnetic field effects the flow greatly. Also, in peristaltic MHD compressor and blood pump machines, MHD effects on conductive physiological fluids become vital [2125].

Nanotechnology is considered important in improving as well as revolutionizing information technology, industry sectors, homeland security, energy, food safety, medicine, transportation, and environmental discipline. In this century, many researchers and mathematicians are working on the development of mathematics and physics of nanofluid mechanics. This is due to the fact that nanotechnology is used in industry to get the optimum outputs in constrained environment. It is, therefore, an attractive research area in the modern fluid mechanics. Nanofluids are composed by dissemination of nanosized materials/particles in the base fluids (both viscous and inviscid liquids). It is known that the term nanofluid was introduced first ever by J. C. Maxwell, a Scottish scientist, in the late nineteenth century. Mostly, modern work is based on the analysis of nanofluids made by the Choi [26]. From this productive study, a vast number of applications of nanotechnology can be found in microchannel cooling and reduction/enhancement in heat transfer. The bond of peristalsis with nanofluids has many utilizations in biomedical science (e.g., radiotherapy for cancer cure and drug delivery), chemical, and mechanical engineering (pumps and transportation of chemicals). Recently, such studies of nanoparticles along with different flow geometries have been examined extensively in physiological flows, such as references [2735].

The double-diffusive convection is an important physical phenomenon arising in fluid mechanics. It is mainly associated with such a convection process where two dissimilar density gradients, having a diverse diffusion rate, are considered. The literature review tells us that none of the analysis has been considered for the double diffusion convection with assumptions of creeping phenomena and low Reynolds number. The double diffusion convection and peristaltic pumping have many applications of the innate mechanism in industrial and chemical engineering. Because of these applications, some authors have contributed to this area with various fluid models. Few are cited in references [3642].

Incorporation of partial slips in fluid flows is essential in the study of polymers and polishing the artificial heart valve. To the best of our knowledge, the slip effects were first used in the peristalsis by Chu and Fang [43]. Later, Akbar et al. [44] further added these effects to examine the influence on peristaltic flow of nanofluid. Also, these effects on hydromagnetic driven peristaltic flow were investigated by Abbasi et al. in reference [45]. More recent works on these effects can be viewed in references [4650].

From the abovementioned discussion, the heat convection impact and magnetic flux on double diffusion convection cannot be neglected. The study of the Oldroyd-4 constant nanofluid is studied in the literature but effects of double diffusion and partial slip with inclined MHD are not studied yet. So, the rationale of our current research is to show how magnetic field and slip boundaries affect peristaltic flow and heat transfer with double diffusion convection.

2. Flow Equations

The equations that describe flow in an incompressible fluid are as follows [36, 37]:

In the abovementioned equations, temperature is represented by , concentration, by , nanoparticle volume fraction, by , represents velocity, illustrates material time derivative, stands for body force, describes stress tensor, is acceleration, represents base fluid density, particles density is represented , describes fluid density at , denotes the thermal capacity of nanoparticles, refers fluid heat capacity, describes fluid volumetrical thermal expansion coefficient, stands for fluid volumetrical solutal expansion coefficient, denotes Dufour diffusively, Soret diffusively is denoted by , is coefficient of Brownian diffusion, coefficient of thermophoretic diffusion, stands for solutal diffusively, and thermal conductivity is represented by

The Oldroyd- constant fluid stress tensor in [18] is defined by the following equations:Where represents relaxation times, is retardation time, is denoted by viscosity, is used for transpose, and is Rivilin–Ericksen tensors.

3. Mathematical Formulation

We assume that an incompressible peristaltic flow and electrically conducting of Oldroyd- constant nanofluid in dimensional conduit having width equal to , in the Cartesian coordinates system. The channel's center is supposed to be along horizontal line, and cross-sectional area is assumed to be beside the vertical line. Conduit boundary is assumed to be moving with constant speed and shape like sinusoidal wave train. Temperatures, solvent concentrations, and nanoparticle concentrations of the lower and upper walls are ( ), ( ), and ( ), respectively. Fixed magnetic field at an angle is applied on the flow. It is assumed that electric field is zero and the Reynolds number is low so that it produces insignificant induced magnetic field as compared to applied magnetic field.

The geometric shape of wall is specified in Figure 1 and mathematical expression is defined as follows[3]:where illustrates wave amplitudes, is channel width, illustrates wavelength, indicates time, and is wave velocity. The range of phase difference ( is . When , the channel is symmetric without a phase wave and at , the channel with a phase wave. Furthermore, the constraints , , , , and satisfy the condition The velocity field in directional and dimensional flow is.

The motion equations in dimensional incompressible flow comprising Oldroyd- constant nanofluid are as follows:

Using the Galilean transformation between fixed and wave frames as follows:

Defining the dimensionless quantities as follows:Where , and stands for thermal Grashof number, Reynolds number, nanoparticle Grashof number, Soret parameter, solutal Grashof number, Prandtl number, Dufour parameter, solutal (species) concentration, Hartmann number, thermophoresis parameter, nanoparticle fraction, Lewis number, temperature, Brownian motion, Lewis number of nanofluid, and wave number, respectively.

Now, using (10) and (11), (4) is identically satisfied and equations (12)–(16) (after omitting bars) in the wave frame becomes as follows:

Now, imposing constraints of (low Reynolds number) and (long wavelength), the equations. (19)-(23) are now reduced as follows:

Now, taking pressure out of (17) and (18) yields the following expression as follows:where the dimensionless equation for is obtained from equation (7)and expressed as follows:in which and . Now, if the model of Oldroyd- constant fluid reduces to viscous fluid.

The expression for (mean flow) is computed in the dimensionless form as follows:where

For the problem under investigation, the slip boundary conditions in dimensionless forms are defined as follows:

If in the parameters, then no slip conditions exist.

3.1. Special Cases

In the absence of slip conditions ( ), , and the findings of reference [4] can also be recovered as a limited case of existing problem.

4. Numerical Solution and Graphical Outcomes

Numerical simulations become essential as analytical solutions can only be found for a limited number of cases. Numerical simulations provide us an alternative mean to understand the problem and its solution without wasting the real resources that are involved in the study. We can develop a comprehensive understanding of the flow situations using modern tools available in the form of software’s like: MATLAB, MATHEMATICA, and ANSYS. The main objective of the present work is to evaluate the consequence of double-diffusion convection and partial slip on magneto-Oldroyd- constants nanofluids with peristaltic propulsion in an asymmetric channel. The exact solutions of the equations (24)-(29) are not easy to evaluate due to coupled and highly nonlinear characteristics. Therefore, the regressive equations are numerically solved in MATHEMATICA by using built-in command ND-Solve. Solutions are utilized in obtaining pressure gradient, pressure rise, and streamlines for different flow parameters. The graphical results are also established employing numerical solutions to verify that numerical data are accurate and to examine the impact of a variety of flow parameters.

4.1. Effects of the Hartmann Number

The Hartmann number is described as ratio of electromagnetic force to viscous force. It is a regular occurrence in fluid flows passing via magnetic fields. Figures 2(a) to 2(c) explain the impact of velocity profile, pressure rise, and pressure gradient on the Hartmann number . As Hartmann number enhances, the magnitude of the velocity profile significantly reduces in channel's center and tends to increase near peristaltic walls (see Figure 2(a)). Furthermore, the profile of velocity has a parabolic shape. In fact, enhancing magnetic number causes the Lorentz force to increase, which tends as a retarding force and slowing fluid motion. Figure 2(b) depicts response to pressure rise. It is illustrated in Figure 2(b) that increasing causes pressure rise to increase in peristaltic (), retrograde (), and free () pumping zones. Moreover, it reduces in the augmented () pumping zone due to higher values of . The pressure gradient continues to decrease as the Hartmann number rises (see Figure 2(c)).

4.2. Effects of Nanoparticle , Solutal , and Thermal Grashof Numbers

In fluid dynamics and heat transfer, the Grashof numbers are described by ratio of buoyancy to viscous force that acts on a fluid. The consequence of nanoparticle, solutal, and thermal Grashof numbers on velocity profile is illustrated in Figures 3(a) to 3(c). The magnitudes of the flow velocity rises when by enhancing , while opposite phenomenon occurs when . Here magnitude of velocity field drops (see Fig. 3(a)). This happens because nanoparticles viscosity drops, causing the velocity to decrease. In Figure 3(b) and 3(c), it is noted that magnitudes of fluid velocity drops when but it increases when by enhancing and . In most instances, thermal buoyancy serves to slow down the flow in the regime. The behavior of pressure rise for , , and is demonstrated in Figures 3(d) to 3(f). It is noted in Figure 3(d) that in all peristaltic regions (augmented (), free ( ), augmented (), and retrograde ()) the pressure rise drops by rising values. On the other hand, and exhibit the opposite trend. Here, pressure rises in all peristaltic zones are increased by boosting and values (see Figure 3(e) and 3(f)). The roll of pressure gradient for and are explained in Figure 3(g) and 3(h). It is indicated in Figure 3(g) that the pressure gradient significantly increases when values of increases. The pressure gradient tends to reduce when increases (see Figure 3(h)).

4.3. Effects of Soret and Dufour Parameters

The outcomes of Soret and Dufour constraints are shown in Figures 4(a) to 4(h). It is shown in Figures 4(a) and 4(b) that and have a similar behavior on the velocity profile, as already explained in Figure 3(a). Figures 4(c) and 4(d) depict the impact of and on temperature profile. The temperature increases by increasing and . It is only because temperature has a direct connection with the constraints of Soret and Dufour. The concentration and nanoparticle fraction profiles decrease by enhancing and values (see Figures 4(e) to 4(h)). Its because random motion reacts with micromixing and random collision tendency of solid nanoparticle, spreading the solid nanoparticles and lessening solute concentration.

4.4. Effects of the Brownian Motion () Parameter

Figures 5(a) to 5(d) look at the impact of the Brownian motion on and . It is clear from Figure 5(a) that by increasing values pressure gradient drops. The temperature and nanoparticle fraction increase by increasing values (see Figures 5(b) and 5(c)). When the Brownian motion gets more vigorous, the nanoparticles are effectively transferred from wall to fluid. The temperature rises as rises as an outcome of this interaction. Furthermore, nanoparticle volume fraction has a direct relationship with . The adverse trends are noted for the case of concentration. Here, by rising values, the concentration drops (see Figure 5(d)). In nature, the nanofluid is just a two-phase fluid, and stochastic mobility of isolated nanoparticles increases energy exchange rates while dropping concentrations in the fluid flow.

4.5. Effects of the Thermophoresis () Parameter

The roll of velocity, concentration, temperature, nanoparticle fraction, and pressure gradient on are indicated in Figures 6(a)–6(e). The outcomes of altering the thermophoresis coefficient on the flow velocity are seen in Figure 6(a). When value of thermophoresis parameters rises the magnitude of flow velocity, it increases in the zone . Furthermore, it tends to fall when . The velocity of fluid is maximum near the channel's center. The role of on temperature is shown in Figure 6(b). It shows that and have an identical behavior of fluid temperature, as already shown in Figure 5(b). The slightly different effects are noted for the case of concentration and nanoparticle fraction (see Figures 6(c) and 6(d)). Here, due to the rising tendency of , the concentration and nanoparticle fraction decreases. The occurrence of pressure gradient is elaborated in Figure 6(e). It is noted in Figure 6(e) that pressure gradient falls as thermophoresis parameter rises.

4.6. Effects of Slip Parameters ()

The outcomes of parameters of slip-on velocity, concentration, temperature, nanoparticle fraction pressure gradient, and pressure rise are shown in Figures 7(a) to 7(g). It is noted in Figure 7(a) that nature of velocity curve is parabolic. Moreover, by increasing parameter of velocity slip, the magnitude of flow velocity tends to fall when and but reverse effects are noted when , here magnitude of flow velocity increases. The impact of temperature slip is shown in Figure 7(b). It is indicated in Figure 7(b) that by enhancing temperature slip temperature drops in the region but it tends to rise in region . The kinetic energy of particles of fluid increases due to slip which rises the fluid temperature. Figures 7(c) and 7(d) depict the influence of parameters of slip concentration and slip nanoparticle fraction . There is a fall in concentration and nanoparticle fraction in the region due to the rising values of slip parameter of concentration and slip factor of nanoparticle fraction. Furthermore, opposite effects are noted in the region (see Figure 7(c) and 7(d)). The particle of fluid is interrupted less by the walls of channel, so concentration drops as the value of increases. Hence, the rate of mass transfer of nanoparticles is slowed. The role of pressure rise on velocity slip is drawn in Figure 7(e). It is indicated in Figure 7(e) that pressure rise reduces in augmented and retrograde pumping zones but increases in augmented region by increasing velocity slip constraints. The pressure gradient increases due to the increasing behaviour of velocity slip constraints and nanoparticle slip factor (see Figures 7(f) and 7(g)).

4.7. Effects of Non-Newtonian Parameters ( and )

To discuss the roll of non-Newtonian parameters and on pressure rise, pressure gradient, and velocity, Figures 8(a) to 8(e) are displayed. It is exhibited in Figures 8(a) and 8(b) that and shows similar behavior on pressure rise. It is illustrated in Figures 8(a) and 8(b) that by increasing and , pressure rise increases in peristaltic and free and retrograde pumping areas but reduces in the augmented region. The pressure gradient is maximum at channel's center, but near the channel wall's pressure gradient drops due to the increasing values of non-Newtonian parameters and (see (Figures 8(c) and 8(d)). It is shown in Figure 8(e) that and have an identical behavior on fluid velocity, as already explain in Figure 2(a).

4.8. Trapping Phenomenon

Trapping is an unusual occurrence in peristaltic propelling flows. It is begun by the development of a fluid mass that internally moves and is enclosed by streamlines of peristaltic wave. Streamlines capture the mass bolus of fluid and move it forward using waves of peristaltic at high flow rates and significant occlusions. It is indicated in Figure 9 that due to the increasing behavior of velocity slip factor size of trapped bolus increases and amount of trapped bolus reduces in both upper and lower portion of channel. By increasing non-Newtonian parameters, streamlines show that the trapped bolus size enhances (see Figure 10). In Figure 11, it is illustrated that the trapped bolus volume grows as the thermal Grashof number rises. The reverse behavior is noted for the Hartmann number case. Here, the size reduces by rising values of (see Figure 12).

Table 1 shows the comparison with the existing literature.

5. Conclusion

The main objectives of this paper are to examine the effect of double-diffusion convection and partial slip with a tilted magnetic field on peristaltic movement in an asymmetric channel for the Oldroyd-4 constants nanofluids. A numerical methodology is used to solve nonlinear system of PDEs. The impact of several physiological parameters on flow quantities is visually depicted. Based on our analysis, it was found that the nanoparticle fraction and concentration profile decrease by boosting slip factors of concentration and nanoparticles. Also, the velocity profiles can be controlled by adjusting the parameters under observation (like slip parameter, non-Newtonian parameter, and diffusivity parameters). Throughout the velocity profiles, it is observed that velocity tend to increase, reach its maximum value, and then decreases to meet the boundary condition while moving from one end to the other. It was also found that the profile of temperature rises as the Brownian motion and thermophoresis constraints increases. This is because growth in the Brownian motion parameter indicates the increase in the kinetic energy of nanoparticles which results in warming up the nanofluid. Physically, thermophoresis demonstrates the temperature difference of nanofluid in the channel walls, and as a result, the temperature rises as the values of the associated parameters are increased. Also, concentration drops as the Brownian motion and thermophoresis constraints increases. Another main finding of the study is that shape of trapped bolus decreases as the Hartmann number increases. This means that by making electromagnetic forces dominant as compared to viscous forces can be effective in drug delivery. Lastly, the confined bolus size grows as the thermal Grashof number rises since buoyancy force becomes dominant as compared to viscous forces.

Nomenclature

:Solutal concentration
:Hartmann number
:Temperature
Pr:Prandtl number
:Grashof number of nanoparticles
:Lewis number
:Nanoparticle volume fraction
:Brownian motion parameter
:Brownian diffusion coefficient
:Thermophoretic diffusion coefficient
Re:Reynolds number
:Thermal Grashof number
:Solutal diffusively
:Dufour parameter
:Thermophoresis parameter
:Nanofluid Lewis number
:Heat capacity of fluid
:Heat capacity of nanoparticle
:Soret parameter
:Dufour diffusively
:Solutal Grashof number
:Soret diffusively

Small alphabets

:Axial velocity
:Transverse velocity
:Acceleration due to gravity
:Wave amplitude
:Channel width
:Pressure
:Thermal conductivity
:Time
:Wave amplitudes
:Propagation of velocity

Greek symbols

:Wavelength
:Stream function
:Nanoparticle mass density
:Nanoparticle heat capacity
:Solutal concentration
:Magnetic field inclination angle
:Volumetric coefficient of thermal expansion
:Volumetric coefficient of solutal expansion
:Concentration slip parameter
:Nanoparticles slip parameter
:Velocity slip parameter
:Temperature slip parameter
:Wave number
:Temperature
:Fluid density
:Fluid density at .

Data Availability

No data were used to support this study.

Conflicts of Interest

All authors declare that there are no conflicts of interest for this manuscript.

Acknowledgments

The authors extend their appreciation to the Deanship of Scientific Research, University of Hafr Al Batin, for funding this work through the research group project no. (0033-1443-S).