Abstract

The principle aim of the current communication is to scrutinize the impact of distinguished effects like variable thermal conductivity and variable molecular diffusivity on non-Newtonian Reiner–Philippoff fluid moving over a stretchable surface. The process of heat transfer is carried out in the presence of nonlinear thermal radiation, variable thermal conductivity, and heat generation/absorption. Furthermore, the study of mass transfer phenomena is carried out in the existence of variable molecular diffusivity. The PDEs regarding our model are renovated into ODEs by utilizing similarity transformation. Furthermore, the dimensionless model is tackled with the help of the RK4 method in conjunction with the shooting technique. The effects of different physical parameters that emerged during the numerical simulation on mass transfer rate, heat transfer rate, and velocity field are portrayed in the form of tables and graphs. It is noteworthy that an elevation in the heat source/sink parameters causes a reduction in the temperature profile. Moreover, a positive variation in the species diffusivity parameter augments the mass fraction field. A variation in the fluid parameter is found to be significantly affecting the shear thinning and shear thickening behaviour of the fluid. Reliability of the numerical outcomes is judged by comparing the obtained outcomes with the already available literature. The article is unique in its sense that the heat and mass transfer analysis of Reiner–Philippoff fluid under the aforementioned effects has not been investigated yet.

1. Introduction

Non-Newtonian fluids have attracted the attention of scientists and engineers in the past years due to their enormous applications in the field of energy and technology. Daily life examples of non-Newtonian fluids are cosmetics, paper production, fiber technology, ketchup, toothpaste, paint, shampoo and blood, wall paint, greases, lubricants, plastic, drilling, mud, etc. All non-Newtonian fluids on the basis of their behaviour in shear are not predicted by a single relation. Several mathematical models have been developed to understand the behaviour of shear stress and strain phenomenon in non-Newtonian fluids. Among all, the Sisko model, Carreau viscosity model, Powell–Eyring model, viscoelastic model, Ellis model, cross viscosity model, and Reiner–Philippoff model are most important models to understand the nature of such fluids. The Reiner–Philippoff model belongs to a class of pseudoplastic/shear thinning fluid. A few number of research studies are available in the literature regarding the boundary layer nature of Reiner–Philippoff fluid. Kapur and Gupta [1] adopted the Karman–Pohlhausen method to find the numerical solution of 2D boundary layer flow of Reiner–Philippoff fluid passing through an inlet length of a straight channel. Ghoshal [2] pondered the impact of non-Newtonian Reiner–Philippoff fluid passing through a circular tube. Na [3] utilized the finite difference method to analyze the boundary layer flow of Reiner–Philippoff fluid over different body shapes like stretching wedge and semi-infinite plate. Yam et al. [4] presented a numerical analysis of boundary layer flow of Reiner–Philippoff fluid flow across a stretching wedge. Ahmed [5] examined the impact of nanoparticles on Reiner–Philippoff nanofluid past a stretching sheet. Ahmed et al. [6] adopted the nonlinear shooting method to find the numerical solution of Reiner–Philippoff fluid past a nonlinear stretching sheet with variable thickness. Reddy et al. [7] successfully applied the shooting method to Reiner–Philippoff fluid flow across a stretching sheet accompanied with transverse applied magnetic field effect. Kumar et al. [8] examined the conduct of Cattaneo–Christov heat diffusion and magnetic field on Reiner–Philippoff fluid over a stretchable surface. Reddy et al. [9] studied the Darcy–Forchheimer flow of Reiner–Philippoff fluid past a linear stretching sheet accompanied with viscous dissipation and thermal radiation.

Radiation is actually the emission or transmission of energy in the form of waves or particles through space or any material medium. Objects emit radiation when electrons from high energy orbit falls down to lower energy orbits. The energy is emitted in the form of electromagnetic radiation. Radiation has immense applications in agriculture, space exploration, law enforcement, geology, glass generation, gas turbines, polymer preparation, spacecraft working at high temperature, furnace design, electricity generation, etc. The linear thermal radiation can be achieved by linearizing the Rosseland approximation [10] with the help of Taylor’s series. Linearized Rosseland approximation has unique Prandtl number [11], whereas in the case of nonlinear Rosseland approximation, three key factors are involved in the problem, termed as Prandtl number, radiation parameter, and temperature ratio parameter. Nonlinear thermal radiation has various applications in industry like polymer production, nuclear reactors, and thermal furnaces. Rosseland approximation [11] for the radiative heat flux has been helpful in computing the heat transfer rate associated with an optically thick fluid. Pantokratoras scrutinized the impact of natural convection along a vertical isothermal plate with linear and nonlinear Rosseland approximation [10]. Linearized version of Rosseland approximation for thermal radiation is used where temperature difference is small. Kho et al. [12] studied the impact of thermal radiation on Williamson nanofluid past a stretching sheet with constant wall temperature. Hayat et al. [13] investigated two-dimensional stagnation point flow of Maxwell nanofluid past a permeable stretchable surface accompanied with thermal radiation. The impact of thermal radiation and chemical reaction on third grade fluid across an exponential stretching sheet was examined by Hayat et al. [14]. Waqas et al. [15] scrutinized the Jeffrey nanofluid flow over stretching sheet embedded with thermal radiation, stratification, and convective boundary conditions. Waqas et al. [16] modelled viscoelastic nanofluid including buoyancy forces, thermal radiation, and convective boundary conditions. Nonlinear thermal radiation is used for both small as well as high temperature difference. Some new literature regarding nonlinear thermal radiation may be found in Refs. [1726].

The ability of material to conduct heat is called thermal conductivity. Thermal conductivity is a property of material that varies with temperature. There are two mechanisms behind thermal conductivity in fluids. First, when the collision of molecules increases, the exchange of energy increases which ultimately helps to transport heat energy through the medium. Second, the random movement of molecules increases in the presence of thermal conductivity. As random movement of molecules increases, it transports heat energy in a particular direction. It is observed that materials having high thermal conductivity are used in heat sink while materials with low thermal conductivity are used as thermal insulation. Thermal conductivity has important applications like electrolytes, steam generators, concrete heating, laminating, catalysis, and molding blow. Hayat et al. [27] scrutinized the conduct of variable thermal conductivity and Cattaneo–Christov heat flux on stratified Oldroyd-B fluid flow towards linearly stretched surface. Khan et al. [28] explored the impact of variable thermal conductivity and nonlinear thermal radiation on nanofluid flow across a moving thin needle. Reddy et al. [29] considered the Williamson nanofluid moving over a stretching surface along with MHD and variable thermal conductivity. Hayat et al. [30] treated the stratified Powell–Eyring fluid moving towards a nonlinear stretching surface along with variable thermal conductivity and Cattaneo–Christov heat flux. Kumar and Varma [31] studied the effect of variable thermal conductivity and variable molecular diffusivity on a nanofluid moving over a variable porous stretching sheet. Recently, some researchers have discussed variable thermal conductivity with different geometries [3236].

Internal heat generation/absorption generates heat energy inside a body by the nuclear process or chemical reaction. There are two models available in literature to scrutinize the effect of heat source/sink on the body. In the first model, internal heat source/sink is uniformly temperature dependent, i.e., . The second model is temperature- and space-dependent heat source/sink model, also called nonuniform heat source/sink model mentioned in energy equation (5). Possible heat generation effects may change the temperature distribution and furthermore the particle deposition. In the present study, the effect of nonuniform heat source/sink has been debated in detail. Heat source/sink in moving fluids has various applications in industry like reactor safety analysis, spent nuclear fuel, combustion analysis, and metal waste. Hayat et al. [37] investigated the impact of heat generation and radiative radiation on fluid flow over a rotating disk along with velocity and thermal slip conditions. The impact of heat generation and magnetic dipole on ferromagnetic Williamson fluid flow over a stretchable surface was examined by Waqas [38]. Salem and Abd El-Aziz [39] scrutinized the effect of internal heat generation/absorption and Hall current on hydromagnetic flow over a stretching surface. Hayat et al. [40] pondered on the behaviour of internal heat generation/absorption, thermal radiation, and double stratification on Maxwell nanofluid fluid moving over a stretchable surface. The behaviour of a two-dimensional laminar boundary layer flow over a stretching surface accompanied with heat source/sink and linear thermal radiation was contemplated by Devi et al. [41]. Muthtamilselvan and Ramya analyzed the impact of heat source/sink on unsteady laminar flow of an incompressible electrically fluid flow towards a porous stretching sheet [42]. The impact of mixed convection, heat generation/absorption, and thermal stratification on an Oldroyd-B fluid past a stretchable surface was examined by Hayat et al. [43]. Few recent attempts in this direction may be represented by [4450].

In the light of aforementioned literature, the purpose of current study is to analyze the behaviour of Reiner–Philippoff fluid past a stretching sheet accompanied with nonlinear thermal radiation, heat source/sink, variable thermal conductivity, and variable molecular diffusivity. The present research fills the gap not adequately addressed in existing literature with the points enumerated underneath.(i)Heat and mass transfer analysis under the effects of variable thermal conductivity, nonlinear thermal radiation, nonuniform heat source/sink, and variable molecular diffusivity on Reiner–Philippoff has not been done yet.(ii)Prior to this article, no comparison analysis in the case of Reiner–Philippoff has been available in the literature to check the reliability of numerical results obtained.(iii)In the available literature, no work has been reported other than [5] to discuss the mass transfer analysis in the case of Reiner–Philippoff fluid.(iv)The skin friction coefficient and especially local Sherwood number in the form of figures and tables on Reiner–Philippoff have not been deeply investigated in the available literature.

2. Mathematical Formulation

One of the traditional portrayals of the stress-deformation conduct of Reiner–Philippoff fluid [59] is

In the above equation, the expression in right side named as the flow function. The distinguished symbols such as , , , and indicate the shear stress, reference shear stress, zero shear viscosity, and limiting viscosity. Reiner–Philippoff fluid belongs to a class of non-Newtonian fluids which exhibits all three, dilatant, Newtonian, and pseudoplastic type, behaviours. The depicted flow function [5] is portrayed bywhere and . The behaviour of the fluid varies with change in . Fluid behaves like Newtonian in case of , dilatant for , and pseudoplastic for .

Figure 1 reflects the geometrical conduct of two-dimensional non-Newtonian Reiner–Philippoff fluid moving over a stretched sheet with stretching velocity acting along . At the surface of the sheet, temperature and concentration are represented by and . The symbols and indicate the ambient temperature and ambient concentration, respectively. The heat generation/absorption along thermal boundary layer is manifested by . The effects like nonlinear thermal radiation and variable thermal are expected to affect the heat transfer process in the energy equation while the variable molecular diffusivity is assumed to affect the mass transfer phenomenon. Equations (3) and (4) depict the continuity and fluid flow phenomenon, whereas heat and mass transfer analysis in the present Reiner–Philippoff model is carried out by (5) and (6). The equations regarding continuity, momentum, energy, and concentration [59] are enumerated underneath:

The boundary conditions associated with the above model are

The variable thermal conductivity and variable molecular diffusivity [21] are given by

The Rosseland radiative heat flux [10, 11] used in (5) is premeditated bywhere the symbols and stand for Stefan–Boltzmann constant and the mean absorption coefficient, respectively.

A similarity transformation [5, 6] under the effect of the stretching velocity is given below:

Using the above variables, fd1(1) and (4)fd4–(6) have been transformed into the ODEs bestowed underneath:

Furthermore, the boundary conditions in the dimensionless form are

The distinguished parameters in (11)–(14) are given by

The surface drag coefficient is given by

The dimensionless form of surface drag coefficient is given by

The heat transfer rate iswhere the term indicates the heat flux which is manifested by

The heat transfer rate in the dimensionless form is taken as

The local Sherwood number is given by the following formula:

The expression represents the mass flux which is delineated by

The dimensionless form of Sherwood number is

3. Solution Methodology

The above nonlinear system of equations (11)–(14) along with the boundary conditions (15) can be tackled by the utilization of the numerical technique termed as shooting method [51] for distinguished parameters arising during numerical simulation of the problem. Figure 2 is designed to describe the flowchart procedure of the shooting method.

For numerical solution, the unbounded domain has been replaced by where is a real number chosen in such a way that the solution does not show any significant variations for . It is noteworthy that assures the expected level of convergence for all the numerical outcomes delineated in this article. The momentum equations (11) and (12) will be tackled collectively by the shooting method and then the temperature and concentration equations will be tackled by using as a known function. Denoting by , by , by , and the missing initial condition by , the momentum equations (11) and (12) are converted into the following system of first-order ODEs:

The above system (25) has been handled numerically with the assistance of the RK4. Furthermore, the missing initial conditions are updated with the help of Newton’s scheme until the criteria stated below are met.where the symbol is a positive number having value and .

In order to solve the temperature equation (13), it is converted into the following system comprising of the first-order differential expressions (25) signifying by and by and using as a known function. The following system of ODEs together with the initial conditions is achieved.

The system of equation (27) is treated the same way as (25) to obtain and .

The concentration equation (14) is transformed into the first-order ODEs by denoting by and by and taking as a known function. The following resulting system of equations is achieved.

Applying the same procedure on (28), the numerical solution of (14) can be obtained in the form of and .

Table 1 presents the comparison analysis of current results with those reported by Reddy et al. [7] by fixing distinguished parameters , , , and , whereas Table 2 depicts the comparison of local Nusselt number for present outcomes with Reddy et al. [9] by fixing , , , , and .

4. Results and Discussion

The current section revealed the impact of various parameters (raised during numerical simulation of the ODEs) on velocity, temperature, concentration profiles, skin friction coefficient, Nusselt number, and Sherwood number; moreover, such effects are portrayed in the form of tables. Table 3 reflects the influence of distinguished physical parameters that were raised during the numerical simulation of the model on the surface drag coefficient, heat transfer, and mass transfer. It is noted that a magnification in the Bingham number creates an augmentation in skin friction coefficient but situation is quite opposite in the case of the Reiner–Philippoff fluid parameter . In the case of heat transfer coefficient, a positive variation in the Reiner–Philippoff fluid parameter , thermal radiation , temperature ratio parameter , and Prandtl number produces an enlargement in the heat transfer rate but a reverse behaviour is monitored for the case of the remaining parameters such as the Bingham number , thermal conductivity , species diffusivity , heat generation/absorption coefficients and , and the Schmidt number . The major factors responsible for an improvement in the mass transfer rate are and . Furthermore, a conflicting effect on the mass transfer rate is being monitored in the case of and . Table 4 shows the fluctuation of distinguished parameters on the Nusselt number under the absence and presence of thermal radiation parameter by fixing and . From the table, it is quite clear that production of heat is more in nonlinear thermal radiation as compared to absence of thermal radiation.

Figure 3 highlights the impact of Bingham number on the velocity field. Bingham number is actually the ratio of the yield stress to the viscous stress. It is observed that the apparent viscosity increases by increasing shear rate. As a result, velocity profile depreciates. Figure 4 shows the impact of Reiner–Philippoff fluid parameter on the velocity profile. This parameter depicts the ratio of the zero shear viscosity to the upper Newton limiting viscosity. It is revealed from the figure that mounting the parameter guides to an abatement in the velocity field. It happens because the increasing values of cause a decrement in the viscosity which allows the fluid to move more freely. Figure 5 is portrayed to establish the impact of the Bingham number on the shear stress field. It is noted that an escalation in the parameter leads to an augmentation in the shear stress field. Physically, viscosity of the fluid decreases in the case of small Bingham number which consequently drives to an enhancement in the fluid shear stress. As a result, shear stress profile increases. As a result, shear stress field increases. The portrayal for nondimensional shear stress field against the Reiner–Philippoff fluid parameter is illustrated in Figure 6. It can be remarked that the shear stress field depreciates owing to an improvement in the parameter . It is quite clear that when the viscosity of fluid heightens, the shear stress rate of the fluid abates which guides to a reduction in the shear stress field. Figures 7 and 8 are plotted to illustrate the impact of heat generation/absorption parameters on the temperature field. It is revealed that enhancing the heat generation and absorption parameters and , more heat is generated inside the fluid which upsurges the thermal boundary layer thickness and ultimately guides to an augmentation in the temperature field. The variation of the temperature field versus the Prandtl number is plotted in Figure 9. Prandtl number is defined as the ratio of momentum diffusivity to thermal diffusivity. Physically, an improvement in the Prandtl number weakens thermal diffusivity which leads to an abatement in the thermal boundary layer thickness, which lessens the temperature within the boundary layer which causes a decrement in the temperature profile. Figure 10 is displayed to analyze conduct of temperature profile against the thermal conductivity parameter . It is found that when the collision of molecules increases, the exchange of energy increases which transports heat energy through the fluid and guides to an improvement in the variable thermal conductivity and temperature profile. Figure 11 is prepared for the analysis of the temperature field due to the diverse values of the thermal radiation parameter. Physically, nonlinear thermal radiation comprises of thermal radiation, temperature ratio parameter, and Prandtl number. Nonlinear thermal radiation is used where high temperature difference is required like glass, polymers, and nuclear reactors. It can be noted that an uplift in the radiation parameter delivers more heat to the fluid which leads to an embellishment in the fluid temperature. The impact of the temperature ratio parameter on the temperature field is portrayed in Figure 12. It is quite understood that by enhancing temperature ratio parameter, more heat is generated which furthermore improves nonlinear thermal radiation and temperature field. Figure 13 is displayed to study the effect of the species diffusivity on the mass fraction field. It is noticed that the growing value of elevates the boundary layer thickness. It is quite clear that species diffusivity is proportional to concentration. A positive variation in the species diffusivity parameter heightens the mass fraction field. Figure 14 is plotted to observe the behaviour of the Schmidt number on the mass fraction field. Schmidt number is defined as the ratio of momentum diffusivity to molecular diffusivity. So, an augmentation in the Schmidt number causes a decrement in the concentration field because molecular diffusivity is inversely related to . Figure 15 is sketched to demonstrate the conduct of Bingham number versus Reiner–Philippoff fluid parameter on the skin friction coefficient. Bingham number is defined as the ratio of yield stress to viscous stress. Physically, when Bingham number increases, the viscosity of fluid decreases (shear thinning) which ultimately guides to an abatement in the Reiner–Philippoff fluid parameter (shear thinning) and fluid moves easily near the surface and an augmentation in the skin friction coefficient takes place. Figure 16 provides the explanation regarding the impact of Reiner–Philippoff fluid parameter in three different cases dilatant , Newtonian , and pseudoplastic versus Bingham number on the skin friction coefficient. It is quite understood that if viscosity of fluid decreases, its surrounding temperature drops gradually. At this stage, the drag coefficient will increase with increase in certain degree of viscosity but opposite behaviour is observed where viscosity increases gradually. In the case of shear thinning (), the viscosity of fluid decreases (shear thinning) which guides to an abatement in the Bingham number (shear thinning) and surface drag coefficient but opposite behaviour is observed in the case of (shear thickening). Figure 17 specifically explains the impact of Prandtl number versus radiation parameter on the heat transfer rate. It is observed that a positive variation in Prandtl number and radiation parameter guides to an improvement in the Nusselt number. Both thermal radiation and Prandtl number are prominent factors of nonlinear thermal radiation which is used where high temperature difference is required. In the presence of nonlinear thermal radiation, heat transfer rate increases. Figure 18 is sketched to interrogate the impact of Prandtl number and temperature ratio parameter on the Nusselt number. Both Prandtl number and temperature ratio parameter are important constituents of nonlinear thermal radiation. The temperature of the fluid increases in the presence of temperature ratio parameter and Prandtl number. As a result, heat transfer rate escalates. Figures 19 and 20 are portrayed to interrogate the impact of space- and temperature-dependent heat source/sink parameter versus Prandtl number on the Nusselt number. It is canvassed that more heat is generated inside the fluid in the case of and . Thermal diffusivity booms because of heat generated inside the fluid by heat source/sink parameter. As a result, Prandtl number decreases which is ratio of momentum diffusivity to thermal diffusivity which leads to a reduction in the heat transfer rate. Figure 21 is designed to discuss the impact of Reiner–Philippoff parameter, i.e., , , and versus Bingham number on the Nusselt number. In the case of shear thinning fluid (), the viscosity of fluid decreases which eventually guides to an increment in temperature and heat transfer rate booms as well but situation is opposite in the thickening fluid where an augmentation in viscosity drives to an abatement in temperature and heat transfer rate as well. Figure 22 is portrayed to investigate the impact of Schmidt number versus variable molecular diffusivity on the Sherwood number. It is observed that the gradual enhancement in tends to weak molecular diffusivity and boundary layer thickness. As a result of this, an augmentation in drives to an improvement in the Sherwood number. Figure 23 is sketched to scrutinize the conduct of variable molecular diffusivity versus Schmidt number on Sherwood number. It is observed that species diffusivity is directly proportional to concentration. It is perceived that a positive variation in species diffusivity guides to an augmentation in concentration. As a result, Sherwood number enhances.

5. Final Remarks

The research examines the impact of nonlinear thermal radiation, heat source/sink, variable thermal conductivity, and variable molecular diffusivity on Reiner–Philippoff fluid past a stretching sheet. The concluding remarks of the present research are enumerated underneath.(i)A mitigation in the shear stress field occurs on account of an escalation in the Bingham number but a transverse behaviour is monitored in the case of fluid parameter .(ii)Surface drag coefficient abates in the case of dilatant fluid but improves in the case of pseudoplastic fluid.(iii)An uplift in the thermal conductivity parameter produces a random movement among the molecules owing to an improvement in the temperature field.(iv)Temperature ratio parameter , radiation parameter , and Prandtl number are three prominent factors of nonlinear thermal radiation responsible for an improvement in the temperature field and local Nusselt number.(v)Heat transfer rate augments in the case of shear thinning fluid but depreciates in the case of shear thickening fluid.(vi)The upshots revealed that a positive variation in the heat generation/absorption parameters and generates more heat which produces an enlargement in the temperature field.(vii)A positive variation in the species diffusivity parameter causes an increment in the mass fraction field.(viii)Augmentation in Schmidt number and variable molecular diffusivity decreases the local Sherwood number.

Nomenclature

:Variable molecular diffusivity parameter
:Concentration of fluid
:Wall concentration
:Ambient concentration
:Fluid parameter
:Mass flux
:Heat generation/absorption
:Schmidt number
:Nusselt number
:Prandtl number
:Temperature ratio parameter
:Shear stress
:Shear viscosity
:Stefan–Boltzmann constant
:Heat capacity of nanoparticles
:Variable thermal conductivity parameter
:Ambient Brownian diffusion coefficient
:Ambient temperature
:Stretching velocity
:Fluid parameter
:Ambient fluid thermal conductivity
:Heat generation/absorption
:Surface heat flux
:Sherwood number
:Radiation parameter
:Radiative heat flux
:Reference shear stress
:Limiting viscosity
:Absorption coefficient.

Data Availability

The data supporting the findings of this study are available within the article.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.