Dynamics of Unsteady Flow of Chemically Reactive Upper-Convected Maxwell Fluid with Temperature-Dependent Viscosity: Keller Box Analysis
The heat and mass transfer characteristics of unsteady flow of incompressible chemically reactive upper-convected Maxwell fluid along a stretching surface in the presence of temperature-dependent viscosity has been studied. The theoretical analysis on heat and mass transfer over a stretching sheet is investigated for numerical analysis. With the use of stream function formulation, the governing boundary layer equations of momentum, energy, and concentration are reduced to a set of linked ordinary differential equations. The nonlinear ordinary differential equations are then solved numerically by using the Keller Box method. The physical behavior of governing parameters on velocity, temperature, and concentration profiles and the local skin friction coefficient and heat and mass transfer rates are graphed and tabulated. The physical impact of Maxwell parameter β, unsteadiness parameter M, Schmidt number Sc, Prandtl number Pr, reaction rate parameter γ, and variable viscosity parameter on the heat and mass transfer has been examined along the stretching surface numerically. The novelty of the present work is to examine the importance of destructive reaction and contractive reaction on the dynamics of upper-convected Maxwell fluid flow in the presence of temperature-dependent viscosity effects. It is observed an interesting behavior of temperature distribution and concentration profile is noted for lower value of viscosity parameter in the presence of chemical reaction. It is also found that skin friction and the rate of heat transfer are decreased by increasing the variable viscosity parameter .
Heat and mass transfer together with the momentum of laminar boundary flow on a stretching sheet is important both theoretically and practically due to its wide application in polymer technology. The flow characteristics of non-Newtonian fluids are quite different in comparison to Newtonian fluids. In order to obtain a clear idea of non-Newtonian fluids and their various applications, it is necessary to study their flow behavior. Heat is a sort of energy that is transferred between different temperature systems or objects. The mass transfer phenomenon is used in various scientific disciplines for different systems and mechanisms that involve molecular and convective transport of atoms and molecules. The dynamics of a material having the properties of elasticity and viscosity when undergoing deformation is a fundamental topic in fluid dynamics. This kind of material referred to as Maxwell fluid has attracted the attention of many researchers due to its wide industrial and technical applications. The upper-convected Maxwell model can be described as the generalization of the Maxwell material for the case of large deformation using the upper-convected time derivative in the coordinate system stretching with the fluid. The study of boundary layer flow over a stretching surface is always essential from a technological point due to its multiple applications in engineering for non-Newtonian fluid flow mechanisms such as printer’s ink, melts, slurries, and groceries. Siddheshwar et al.  worked on forced flow mechanism for Newtonian liquids along the various stretchable geometries numerically. Aman et al.  developed a fractional mathematical model of Maxwell fluid along the porous geometry with second-order slip effects numerically. A similarity analysis was presented by Abd El-Aziz  to investigate a special distribution of the stretching velocity and surface heat flux for unsteady boundary layers over a stretched sheet. Abel et al.  performed a nonisothermal analysis on the stretching shape with heat flux and variable viscosity effects numerically.
In all of the mentioned studies, fluid viscosity has been assumed to be a constant function of temperature within the boundary layer. However, it is known that the physical properties of the fluid may change significantly when exposed to internal generated temperature. Viscosity depends strongly on temperature. In liquids, it usually decreases with increasing temperature, whereas in gases, viscosity increases with increasing temperature. Hassanein  developed a mathematical analysis for continuous shapes to examine heat transfer with variable viscosity. Khan et al.  analyzed the viscosity change in the Maxwell fluid mechanism under the influence of solutal stratification effects. Khan et al.  discussed the technical note for thin film flow behavior over shrinking geometry with thermal properties. Mukhopadhyay and Bhattacharyya  discussed the effects of chemical reactions on unsteady mixed convection flow and heat transfer over a porous stretching sheet. Some authors [9, 10] studied the effects of ohmic dissipation on viscoelastic flow and heat transfer under different conditions. Cortell  examined the numerical phenomena on second-grade viscoelastic flow to predict the mass and momentum transfer behavior across the porous mechanism theoretically. He obtained valid solutions compared to the previous mentioned results. An analytical solution of chemically reacting species in viscous flow across the plane elastic geometry under homogeneous reactions has been illustrated by Andersson et al. . The technical problem about thermocapillary over a stretchable shape for thin film flow has been studied in [13, 14].
Elbashbeshy and Bazid  examined the velocity behavior of flow for unsteadiness along the stretching geometry numerically. Tsai et al.  studied the quiescent fluid medium problem with heat source effects across the unsteady stretching shape. The flow behavior of Maxwell fluid in a channel has been obtained theoretically in . Chamkha et al.  produced similar results in porous shape for mass and heat behavior with chemical reaction. The mass transfer characteristics along the stretching geometry under slip and suction/injection effects by using shooting technique have been obtained by Bhattacharyya et al. [19, 20]. Gupta and Gupta  illustrated the mass transfer phenomena with blowing or suction impacts along the stretching surface mathematically. Ishak et al.  investigated the heat mechanism on permeable geometry with stretching and wall temperature effects numerically. The combined effects of transpiration and suction/blowing on unsteady mixed convection flow via a stretching sheet were studied by Mukhopadhyay and Vajravelu .
The flow of fluids may be categorized into Newtonian fluid flow, viscoelastic fluid flow, a flow of non-Newtonian fluid with time-dependent viscosity, and the flow of non-Newtonian fluids with time-independent viscosity. However, there exists no fluid flow which accurately fits the definition of a Newtonian fluid flow (i.e., shear stress is proportional to the rate of deformation). By taking an idea from Mukhopadhyay and Bhattacharyya , the unsteady two-dimensional flow of Maxwell fluid in the presence of variable viscosity and chemical reaction has been studied. The governing partial differential equations are transformed into ordinary ones using appropriate transformations, and the reduced ordinary differential equations are solved numerically using the Keller Box scheme. With the use of graphical representations, the impacts of controlling factors on velocity, temperature, and concentration fields with their slopes are explored and analyzed.
2. Research Methodology
Over an unstable stretched sheet, we study laminar boundary-layer two-dimensional flow and mass transfer of an incompressible non-Newtonian Maxwell fluid. Let represent the concentration on the sheet surface and represent the concentration away from the sheet. Furthermore, the species’ reaction is a first-order homogeneous chemical reaction with a rate that varies with time. We assume that the fluid and mass fluxes are constant at time . At , the unsteady fluid and mass fluxes begin. The sheet emerges from a slit at the origin and moves with a nonuniform velocity , where are positive constants with dimensions (time)-1, is the initial stretching rate, and is the effective stretching rate that increases with time. The material properties of the extruded sheet may change over time in the case of polymer extrusion. The governing equations are given as follows :
In Figure 1, and are the components of velocity in and directions, are the fluid’s dynamic and kinematic viscosity, is the fluid’s concentration (of the species), and is the diffusing species’ diffusion coefficient. Furthermore, is the time-dependent reaction rate, where denotes destructive reaction and denotes constructive reaction, is a constant, is the period’s relaxation time, is a constant, denotes temperature, and denotes free-stream temperature.
The suitable boundary conditions for the problem are given as follows:
Here, is concentration at the sheet surface equal to , where is the concentration far away from the surface, is variable temperature, is constant, is the thermal conductivity of fluid, and is the specific heat.
2.1. Stream Function Formulation
The following physical dimensionless similarity transformations along with the stream function formulation were used by Mukhopadhyay [8, 23] to reduce the governing partial differential equations into the system of ordinary differential equations:where is the stream function and is the similarity variable.
With the help of the above relations, the governing partial differential equations are finally reduced to nonlinear ordinary differential equations given as follows:where Pr is the Prandtl number, is the similarity variable, is the dimensionless temperature, is the kinematic fluid viscosity, is the variable viscosity parameter, is the Schmidt number, is the unsteadiness parameter, is the Maxwell parameter, and is the reaction rate parameter. Moreover, shows the destructive reaction and shows reaction is constructive, but corresponds to no reaction.
The boundary conditions in (5) then becomes
The mathematical expression for skin friction is , the Nusselt number is , and the value of Schmidt number is defined as the ratio of momentum diffusivity and mass diffusivity. It is used to characterize fluid flows in which there are simultaneous momentum and mass diffusion convection process. The mathematical expression for the Sherwood number is . The values of , and are .
Thus, the values of skin friction, Nusselt number, and Sherwood number are given by .
2.2. Computational Technique
The coupled partial differential equations for the above model are reduced into a set of ordinary differential equations with the help of suitable stream function formulation. The reduced nonlinear ordinary differential equations (8)–(10) along with the boundary conditions in (11) are solved numerically by using the Keller Box scheme with an iterative method. We start by introducing new independent variables , and by using the following equation:
To overcome the difficulty of solving the equations, equations (8)–(11) become
The reduced boundary conditions are as follows:
We now consider the segment and with as the midpoint given as follows:
The central difference form and average form are given as follows:
By applying equations (18)–(22), governing equations (13)–(16) becomealong with boundary conditions
Now, by applying the iterative Newton–Raphson method for smooth algorithm, we get
As in the standard Newton–Raphson technique, we neglect all appearances of powers of greater than the first power and hence equations become
Again, by using equations (27)–(31) in equations (23)–(26), the reduced form of equations is given as follows:
We recall the boundary conditions that can be satisfied exactly with no iteration. Therefore, in order to keep up these correct values in all iterates, we take
Now, this is a very important step to arrange the above difference equations in matrix form. If it is carried out incorrectly, then either the method becomes very inefficient due to the absence of any discernible structure in the matrix or the matrix solution method breaks down due to having a singular matrix (determinant = 0) or submatrix. The matrix form is given as
3. Analysis and Discussion of Results
The numerical computations of Maxwell fluid are executed for several values of dimensionless parameters involved in the equations such as variable viscosity parameter , the Prandtl number (Pr), the Maxwell parameter , and the unsteadiness parameter M. To illustrate the computed results, some figures are plotted and physical explanations are given for velocity , concentration , and temperature with their slopes. Table 1 presents the validated numerical data for skin friction for accuracy of results by comparing them with existing results available in the study by Chamkha et al.  and Mukhopadhyay and Bhattacharyya . The validated results with good agreement are found by keeping unsteady flow of incompressible Newtonian flow fluid . The most favorable numerical results are obtained by adding the energy equation to the present phenomena.
Figures 2(a)–2(c)) present the physical behavior of , and profiles for numerous values of , and with some fixed parameters. From Figure 2(a), it is noted that the fluid velocity becomes maximum at a higher value of while the minimum value of velocity profile is examined at a lower value of . In Figure 2(b), an interesting behavior of the temperature profile is reported for various values of with chemical reaction effects. The concentration profile increases with lower value of but decreases as increases with prominent variations in Figure 2(c) in the presence of chemical reaction. The maximum value of concentration is noted at lower . The temperature profile is increased with suitable variations at but decreases as the parameter decreases significantly. In the presence of variable viscosity, the prominent behavior in concentration is illustrated for each value and temperature distribution is also found to be in good agreement.
Figures 3(a)–3(c) demonstrate fluid velocity profile, fluid temperature profile, and fluid concentration profile for the various values of Prandtl number Pr , and by keeping other parameters , and constant along the stretchable surface. The velocity profile increases with a lower value of Pr and gradually decreases with increasing Pr in Figure 3(a). The maximum increasing behavior in temperature is depicted at a lower value of Pr , but this decreases as the value of Pr increases with prominent variations in Figure 3(b) significantly. The similar trend in concentration is examined for each value of Pr graphically, but the concentration increases with decreasing Pr values by satisfying the given boundary conditions in Figure 3(c). In the presence of viscosity effects, a suitable change is examined in velocity and temperature but a similar trend is examined in concentration for . The suitable response in concentration is illustrated in the presence of . The fluid flow caused solely, and the concentration is maximum than free-stream concentration along the stretching sheet.
From Figures 4(a)–4(c), it is concluded that the fluid velocity is maximum along the surface of a given shape for a lower value of and the minimum value is obtained at the largest value of with good response. The velocity profile showed good variations in the presence of in Figure 4(a). In Figure 4(b), the temperature profile is represented graphically for different values of , and by keeping other parameters constant. In Figure 4(c), it is noted that the concentration distribution has a similar behavior at lower values of but increasing change is noted as increases. It is observed that temperature shows similar behavior for each value of and gradually increases with increasing . The boundary layer thickness is reduced by increasing the effect of Maxwell parameter .
In Figures 5(a)–5(c), the velocity, temperature, and concentration profiles are drafted against various values of unsteadiness parameter ,and with some fixed parameters , and Pr . The most favorable change in each profile is noted against the unsteadiness parameter in the presence of viscosity and chemical reaction effects. The velocity, temperature, and concentration increase at a minimum value of but decrease gradually with increasing with prominent variations in the presence of . The increasing effect of unsteadiness parameter exhibits the reduction in thickness of thermal boundary layer near the wall of stretching surface. In Figures 6(a)–6(c), the effects of viscosity parameter on skin friction , heat transfer , and mass transfer are clearly shown. The maximum value of and is noted at a lower value of viscosity parameter , and this gradually decreases with increasing in the presence of chemical reaction. The decreasing behavior in mass transfer is drafted by increasing but increases as decreases significantly.
The graphical and numerical outcomes of Maxwell fluid with variable viscosity and chemical reaction effects are executed for several values of dimensionless parameters involved in the equations such as the variable viscosity parameter , the Prandtl number (Pr), the Maxwell parameter , and the unsteadiness parameter M. To illustrate the computed results, some figures are plotted and physical explanations are given for velocity , concentration , and temperature with their slopes skin friction , mass transfer , and heat transfer , respectively. The main outcomes are given as follows:(i)The fluid velocity is depicted maximum at a higher value of , while the minimum value of velocity profile is examined at a lower value of (ii)The similar trend in concentration is examined for each value of Pr graphically but it increases with decreasing Pr by satisfying the given boundary conditions(iii)The increasing effect of the unsteadiness parameter exhibits the reduction in thickness of the thermal boundary layer near the wall of stretching surface(iv)The maximum increasing behavior in temperature is depicted at a lower value of Pr , but it decreases as the value of Pr increases with prominent variations.(v)The maximum value of and is noted at a lower value of viscosity parameter , and it gradually decreases with increasing in the presence of chemical reaction.
No data were used to support this study.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
P. Siddheshwar, G. Sekhar, and A. Chethan, “Flow and heat transfer in a Newtonian liquid with temperature dependent properties over an exponential stretching sheet,” Journal of Applied Fluid Mechanics, vol. 7, no. 2, pp. 367–374, 2014.View at: Google Scholar
S. Aman, Q. Al-Mdallal, and I. Khan, “Heat transfer and second order slip effect on mhd flow of fractional Maxwell fluid in a porous medium,” Journal of King Saud University Science, vol. 32, no. 1, pp. 450–458, 2020.View at: Publisher Site | Google Scholar
M. Abd El-Aziz, “Unsteady mixed convection heat transfer along a vertical stretching surface with variable viscosity and viscous dissipation,” Journal of the Egyptian Mathematical Society, vol. 22, no. 3, pp. 529–537, 2014.View at: Publisher Site | Google Scholar
M. S. Abel, S. K. Khan, and K. V. Prasad, “Study of visco-elastic fluid flow and heat transfer over a stretching sheet with variable viscosity,” International Journal of Non-linear Mechanics, vol. 37, no. 1, pp. 81–88, 2002.View at: Publisher Site | Google Scholar
I. A. Hassanien, “The effect of variable viscosity on flow and heat transfer on a continuous stretching surface,” ZAMM - Journal of Applied Mathematics and Mechanics/Zeitschrift für Angewandte Mathematik und Mechanik, vol. 77, no. 11, pp. 876–880, 1997.View at: Publisher Site | Google Scholar
M. Khan, M. Y. Malik, T. Salahuddin, S. Saleem, and A. Hussain, “Change in viscosity of Maxwell fluid flow due to thermal and solutal stratifications,” Journal of Molecular Liquids, vol. 288, Article ID 110970, 2019.View at: Publisher Site | Google Scholar
Y. Khan, Q. Wu, N. Faraz, and A. Yildirim, “The effects of variable viscosity and thermal conductivity on a thin film flow over a shrinking/stretching sheet,” Computers & Mathematics with Applications, vol. 61, no. 11, pp. 3391–3399, 2011.View at: Publisher Site | Google Scholar
S. Mukhopadhyay and K. Bhattacharyya, “Unsteady flow of a Maxwell fluid over a stretching surface in presence of chemical reaction,” Journal of the Egyptian Mathematical Society, vol. 20, no. 3, pp. 229–234, 2012.View at: Publisher Site | Google Scholar
K. R. Rajagopal, T. Y. Na, and A. S. Gupta, “Flow of a viscoelastic fluid over a stretching sheet,” Rheologica Acta, vol. 23, no. 2, pp. 213–215, 1984.View at: Publisher Site | Google Scholar
M. S. Abel, E. Sanjayanand, and M. M. Nandeppanavar, “Viscoelastic Mhd Flow and Heat Transfer over a Stretching Sheet with Viscous and Ohmic Dissipations,” Communications in Nonlinear Science and Numerical Simulation, vol. 13, no. 9, pp. 1808–1821, 2008.View at: Publisher Site | Google Scholar
R. Cortell, “Toward an understanding of the motion and mass transfer with chemically reactive species for two classes of viscoelastic fluid over a porous stretching sheet,” Chemical Engineering and Processing - Process Intensification, vol. 46, no. 10, pp. 982–989, 2007.View at: Publisher Site | Google Scholar
H. I. Andersson, O. R. Hansen, and B. Holmedal, “Diffusion of a chemically reactive species from a stretching sheet,” International Journal of Heat and Mass Transfer, vol. 37, no. 4, pp. 659–664, 1994.View at: Publisher Site | Google Scholar
B. S. Dandapat, B. Santra, and H. I. Andersson, “Thermocapillarity in a liquid film on an unsteady stretching surface,” International Journal of Heat and Mass Transfer, vol. 46, no. 16, pp. 3009–3015, 2003.View at: Publisher Site | Google Scholar
B. Dandapat, B. Santra, and K. Vajravelu, “The effects of variable fluid properties and thermo capillarity on the flow of a thin film on an unsteady stretching sheet,” International Journal of Heat and Mass Transfer, vol. 50, no. 5-6, pp. 991–996, 2007.View at: Publisher Site | Google Scholar
E. M. A. Elbashbeshy and M. A. A. Bazid, “Heat transfer over an unsteady stretching surface,” Heat and Mass Transfer, vol. 41, no. 1, pp. 1–4, 2004.View at: Publisher Site | Google Scholar
R. Tsai, K. H. Huang, and J. S. Huang, “Flow and heat transfer over an unsteady stretching surface with non-uniform heat source,” International Communications in Heat and Mass Transfer, vol. 35, no. 10, pp. 1340–1343, 2008.View at: Publisher Site | Google Scholar
T. Hayat and Z. Abbas, “Channel flow of a Maxwell fluid with chemical reaction,” Zeitschrift für Angewandte Mathematik und Physik, vol. 59, no. 1, pp. 124–144, 2008.View at: Publisher Site | Google Scholar
A. J. Chamkha, A. M. Aly, and M. A. Mansour, “Similarity solution for unsteady heat and mass transfer from a stretching surface embedded in a porous medium with suction/injection and chemical reaction effects,” Chemical Engineering Communications, vol. 197, no. 6, pp. 846–858, 2010.View at: Publisher Site | Google Scholar
K. Bhattacharyya, S. Mukhopadhyay, and G. C. Layek, “Slip effects on an unsteady boundary layer stagnation-point flow and heat transfer towards a stretching sheet,” Chinese Physics Letters, vol. 28, no. 9, Article ID 094702, 2011.View at: Publisher Site | Google Scholar
K. Bhattacharyya, “Effects of radiation and heat source/sink on unsteady mhd boundary layer flow and heat transfer over a shrinking sheet with suction/injection,” Frontiers of Chemical Science and Engineering, vol. 5, no. 3, pp. 376–384, 2011.View at: Publisher Site | Google Scholar
P. S. Gupta and A. S. Gupta, “Heat and mass transfer on a stretching sheet with suction or blowing,” Canadian Journal of Chemical Engineering, vol. 55, no. 6, pp. 744–746, 1977.View at: Publisher Site | Google Scholar
A. Ishak, R. Nazar, and I. Pop, “Heat transfer over an unsteady stretching permeable surface with prescribed wall temperature,” Nonlinear Analysis: Real World Applications, vol. 10, no. 5, pp. 2909–2913, 2009.View at: Publisher Site | Google Scholar
S. Mukhopadhyay and K. Vajravelu, “Effects of Transpiration and Internal Heat Generation/absorption on the Unsteady Flow of a Maxwell Fluid at a Stretching Surface,” Applied Mechanics, vol. 79, no. 4, Article ID 044508, 2012.View at: Publisher Site | Google Scholar