In this study, we investigated the effects of nonlinear thermal radiation and non-uniform heat source/sink in unsteady three-dimensional flow of Carreau and Casson fluids past a stretching surface in the presence of homogeneous–heterogeneous reactions. The transformed governing equations are solved numerically using Runge–Kutta based shooting technique. We obtained good accuracy of the present results by comparing with the already published literature. The influence of dimensionless governing parameters on velocity, temperature and concentration profiles along with the friction factors, local Nusselt and Sherwood numbers is discussed and presented graphically. We presented dual solutions for flow, heat and mass transfer in Carreau and Casson fluids. It is found that the heat and mass transfer rate in Casson fluid is significantly high while compared with the Carreau fluid.
MHD; Homogeneous–heterogeneous reaction; Nonlinear thermal radiation; Carreau fluid; Casson fluid; Non-uniform heat source/sink
- velocity components in directions respectively ()
- distance along the surface ()
- distance normal to the surface ()
- specific heat capacity at constant pressure ()
- dimensionless velocities
- Temperature of the fluid ()
- concentration of the fluid ()
- non uniform heat source/sink ()
- acceleration due to gravity ()
- thermal conductivity ()
- diffusion coefficient ()
- pressure ()
- similarity variable
- electrical conductivity ()
- Stefan–Boltzmann constant ()
- mean absorption coefficient ()
- volumetric thermal expansion ()
- concentration expansion coefficient ()
- dimensionless temperature
- density ()
- kinematic viscosity ()
- dynamic viscosity ()
- rate constants
- power-law index parameter
- non-uniform heat generation/absorption coefficients
- skin friction coefficient in x-direction
- skin friction coefficient in y-direction
- local Nusselt number
- local Sherwood number
- local Reynolds number
- Prandtl number
- Schmidt number
- unsteadiness parameter
- the ratio of temperatures
- Weissenberg number
- positive constant
- diffusion coefficients
- concentration of chemical species
- rate constants
- ratio of diffusion coefficients
- Casson fluid parameter
- magnetic induction parameter
- magnetic field parameter
- stretching ratio parameter
- strength of heterogeneous reaction parameter
- strength of homogeneous reaction parameter
- fluid
- condition at the wall
- condition at the free stream
Convection boundary layer flow over a stretching sheet has relevance in many engineering processes such as drawing of plastic films, tinning and annealing of copper wires and electrolyte paper production. Due to these applications Sakiadis [1] started an analysis on the flow past a stretching sheet. Thereafter, the researchers [2], [3], [4], [5] and [6] continued their research on the flow over a stretching sheet under various interesting aspects. On the other hand flow of non-Newtonian fluids encountered in several large-scale industrial applications including blood flows in micro-circulatory system, food and polymer processing, magma and ice flows. Due to flow diversity in the environment a single mathematical model does not overcome all the rheological fluid properties associated with non-Newtonian fluids. Thus various constitutive equations for such fluids are available in already existing literature [7] and [8]. Additionally, the power-law Carreau fluid is also one of the non-Newtonian fluids. Carreau fluid model is valid for viscous, high and low shear rates. Because of this advancement, it has benefitted in many technological and manufacturing flows. Keeping this into view, Zhu [9] discussed the mass transfer characteristics of Carreau fluid over a swarm of Newtonian drops. The time dependent Poiseuille flow of a Carreau fluid in the presence of slip effect was investigated by Georgiou [10] and concluded that the wavelength and amplitude of oscillations in radial direction are decreased with an increase in the slip effect. Abd El Naby et al. [11] considered the peristaltic flow characteristics of Carreau fluid in uniform tube and discussed the heat transfer characteristics of Carreau fluid.
Magnetohydrodynamic (MHD) is the mechanical property of fluids, which describes the motion of highly conducting fluid with existing magnetic field. The conducting fluids generate an electric current due to fluid flow and this force boosts up the mechanical properties of fluid. The peristaltic flow with induced magnetic field also has major applications in the psychological fluids: like peristaltic MHD compressor, blood pumping machines and the blood flows, and these were analytically studied by [12] and [13]. MHD convection flow of an unsteady EG-Nimonic nano fluid past a vertical plate was examined by Sandeep et al. [14]. Further, Raju et al. [15] extended this work for ferrofluids by considering the non-uniform heat source/sink and aligned magnetic field effects. The flow of peristaltic Carreau nanofluid past an asymmetric channel was numerically investigated by Akbar et al. [16] and found that increasing values of magnetic field parameter encourages the velocity profiles. Later on, the researchers [17], [18], [19] and [20] reported experimental and theoretical studies on non-Newtonian fluid flows with different flow configurations by considering different boundary conditions. Casson fluid is also a non-Newtonian fluid, which is a shear thinning liquid and exhibits the yield shear stresses. If yield stress is greater than the shear stress then it acts as a solid, whereas if yield stress lesser than the shear stress is applied then the fluid would start to move, for example honey, tomato sauce, fruit juices and human blood. It has various applications in fibrinogen, cancer homeo-therapy, protein and red blood cells form a chain type structure. Due to these applications many researchers are concentrating characteristics of Casson fluid, which are given in Refs. [21], [22] and [23].
The heat and mass transfer in the flows over a stretching sheet with homogeneous–heterogeneous reaction has a major role in metallurgy and chemical engineering industries, such as polymer production and food processing. Moreover, coupled heat and mass transfer problems in the presence of homogeneous–heterogeneous reaction are of importance in many processes, and therefore it is a considerable amount of attention in recent days. Therefore some of the possible applications can be found in the processes such as drying, damage of crops due to freezing, distribution of temperature and moisture over agricultural fields and groves of fruit trees, evaporation at the surface of a water body and energy transfer in a wet cooling tower. Hayat et al. [24] discussed Carreau fluid flow past a convectively heated stretching surface and concluded that the velocity profiles are improved with the material parameter. Convection flow of non-Newtonian MHD flow past a permeable exponentially stretching sheet was numerically investigated by Raju et al. [25]. The flow through stretching sheet has great attention due to its importance in various fields such as MHD accelerators, generators, pumps and flow meters, and design of cooling systems and these are given by Akbar et al. [26]. Jenny et al. [27] studied the Rayleigh–Benard flow for convection rolls in Carreau fluids and analyzed the momentum and heat transfer behavior of Carreau fluid. Jasmine Benazir et al. [31] examined an unsteady magnetohydrodynamic flow due to vertical cone with non-uniform heat generation/absorption. The variable conductivity effect on hydrodynamic convection flow due to cone in the presence of chemical reaction was investigated numerically by Rushikumar et al. [32]. Shehzad et al. [33] discussed the analytical solutions of slip effect on peristalsis flow over a curved channel in the presence of radial magnetic field. Dufour effect on radiative unsteady magnetohydrodynamic flow past a vertically accelerated wavy plate with variation of double diffusion effects was illustrated by Prakash et al. [34]. Recently, the researchers [35], [36], [37], [38] and [39] discussed the non-Newtonian fluid flow of stretching sheet with various effects and various boundary conditions. In these studies they found very interesting solutions as the non-uniform heat source/sink parameter has tendency to control the temperature profiles and also the non-Newtonian fluids are regulating the temperature profiles of the flow.
In most of the studies the radiation can be treated as a constant in the flow region. In real time it is very difficult to maintain the constant temperature entire flow region. But, in this study, we proposed a mathematical model for analyzing the effects of nonlinear thermal radiation on three-dimensional flow of Carreau and Casson fluids past a stretching surface with homogeneous–heterogeneous reactions and non-uniform heat source/sink. The nonlinear thermal radiation has importance in various industrial as well as science and technological applications. The transformed governing equations are solved numerically using Runge–Kutta based shooting technique. We presented dual solutions for the flow of Carreau and Casson fluids over a stretching sheet.
Consider an unsteady three dimensional flow of Carreau and Casson fluids past a stretching sheet in the presence of nonlinear thermal radiation and non-uniform heat source/sink. For improving mass transfer we also considered the homogeneous–heterogeneous reactions. The flow is restricted to z-direction as displayed in Fig. 1. In this study we skip the induced magnetic field and viscous dissipation effects. The flow is due to stretching surface. The rheological model for an isotropic flow of Casson fluid is [23] given as follows:
|
(1) |
|
(2) |
|
Figure 1. Schematic representation of physical model.
|
In the above equation and is the th component of the deformation rate, the product of the component of deformation rate with itself, is a critical value of this product based on the non-Newtonian model, is the plastic dynamic viscosity of non-Newtonian fluid, and is the yield stress of the fluid. The anonymous researchers have suggested the value of n = 1. However, in many applications this value is n > 1.
The flow is generated due to the linear stretching sheet. The extra stress tensor for Carreau fluid is given by [26]
|
(3) |
Here is the extra stress tensor, is the zero shear rate viscosity, is the time constant, is the power-law index and is defined as
|
(4) |
where Π is the second invariant strain tensor. According to above assumptions the flow analysis of Carreau and Casson fluids equations is given by
|
(5) |
|
(6) |
|
(7) |
where are the velocity components along the directions respectively. ν is the kinematic viscosity coefficient, is the Casson fluid parameter, is the time constant, is the density of the fluid and is the electric conductivity with the boundary conditions
|
(8) |
here are the stretching velocities near the surface. To convert the nonlinear partial differential equations for velocities, we are now introducing the similarity transformations as follows:
|
(9) |
Here in Eq. (9) automatically satisfy the continuity equation, and by using Eqs. (9), and (5), (6) and (7) are given by
|
(10) |
|
(11) |
The transformed boundary conditions are as follows:
|
(12) |
here A is the unsteadiness parameter, We is the Weissenberg number, M is the magnetic field parameter and is the stretching ratio parameter.
|
(13) |
The boundary layer thermal energy equation with nonlinear thermal radiation and non-uniform heat source/sink is given by
|
(14) |
with the boundary conditions
|
(15) |
The non-dimensional temperature parameters are given by
|
(16) |
where T is the fluid temperature, are near the fluid temperature and the far away from the fluid temperature, k is the thermal conductivity of the fluid, is the specific heat capacitance at constant pressure, is the concentration susceptibility and is the Stefan–Boltzmann constant.
The time dependent non-uniform heat source/sink is defined as follows:
|
(17) |
The positive values of of above equation correspond to heat generation and negative values corresponds to heat absorption.
Using Eqs. (16) and (17), (14) and (15) are reduced to
|
(18) |
With the transformed boundary conditions
|
(19) |
where Pr is the Prandtl number, A is the unsteadiness parameter, R is the thermal radiation parameter, and is the ratio of temperatures which are given by
|
The boundary layer equation for conservation of mass in the presence homogeneous–heterogeneous reactions is given by
|
(20) |
|
(21) |
The corresponding boundary conditions are
|
(22) |
The non-dimensional are given by
|
(23) |
where are the diffusion coefficient of the species , is the positive constant, are the concentrations of the chemical species, are the rate constants, and we assume that both the chemical reactions are isothermal. Using similarity transforms (23), Eqs. (20) and (21) reduce to
|
(24) |
|
(25) |
The corresponding boundary conditions are
|
(26) |
where Sc is the Schmidt number, is the strength of the heterogeneous-reaction parameter, δ is the ratio of diffusion coefficient, and K is the strength of the homogeneous reaction parameter, which are given by
|
(27) |
For the most of real time practical applications we assume that chemical diffusion coefficients are almost same size. So diffusion coefficients are equal. Then and in this case we get
|
(28) |
Then by substituting Eq. (25) in (21) and (22), we get
|
(29) |
Subject to the boundary conditions:
|
(30) |
For physical quantities of interest the friction factor coefficients along directions, local Nusselt and Sherwood numbers are given by
|
(31) |
|
(32) |
|
(33) |
where is the Reynolds number.
The set of nonlinear ordinary differential Eqs. (10), (11), (18) and (29) corresponding to the boundary conditions (12), (19) and (30) are solved numerically using Runge–Kutta based shooting technique (Sandeep and Sulochana [30]). Results display the influence of non-dimensional governing parameters on velocity, temperature and concentration profiles along with the friction factors, local Nusselt and Sherwood numbers. For numerical results we considered the non-dimensional parameter values as , These values are kept as common in entire study except the variations in respective figures and tables. In graphical results red color profiles indicate the flow of Carreau fluid and green color profiles indicate the flow of Casson fluid.
Figure 2, Figure 3, Figure 4 and Figure 5 depict the influence of Weissenberg number on velocity, temperature and concentration fields for both Carreau and Casson fluids. It is found that increasing values of the Weissenberg number enhances the thermal boundary layer and decreases the momentum concentration boundary layers. Physically, Weissenberg number is directly proportional to the time constant and inversely proportional to the viscosity. The time constant to viscosity ratio is higher for larger values of Weissenberg number. Hence, higher Weissenberg number causes to enhance the thermal boundary thickness.
|
Figure 2. Velocity filed for different values of Weissenberg number.
|
|
Figure 3. Velocity field for different values of Weissenberg number.
|
|
Figure 4. Temperature field for different values of the Weissenberg number.
|
|
Figure 5. Concentration field for different values of the Weissenberg number.
|
The dimensionless temperature distribution for different values of radiation parameter R is shown in Fig. 6 for both Carreau and Casson fluids. It reveals that the greater values of radiation parameter enhance the temperature boundary layer thickness. Generally, for higher values of radiation parameter produces more heat to working fluid that shows an enhancement in the temperature field. We have noticed an enhancement in the temperature profiles of both Carreau and Casson fluids due to increase in the radiation parameter. The ratio of temperature on temperature profiles is shown in Fig. 6. It is clear that increasing values of temperature ratio parameter improves the temperature profiles of the flow. This may happen due to increasing thermal conductivity of the flow (see Fig. 7).
|
Figure 6. Temperature field for different values of radiation parameter.
|
|
Figure 7. Temperature field for different values of ratio of temperature parameter.
|
The effects of the magnetic field on velocity, temperature and concentration fields are displayed in Figure 8, Figure 9, Figure 10 and Figure 11. We observed depreciation in the velocity, concentration fields and enhancement in the temperature field for increasing values of magnetic field parameter. This proves the general physical behavior of M that improved values of M depreciate the velocity fields. Physically, the drag force increases with an increase in the magnetic field and as a result depreciation occurs in the velocity field. The influence of unsteadiness parameter on velocity, temperature and concentration fields is exhibited in Figure 12, Figure 13, Figure 14 and Figure 15 for both Carreau and Casson fluids. We detect from the figure that the velocity and concentration fields are enhanced and declined the temperature field with increasing values of unsteadiness parameter. Physically, increasing values of unsteadiness parameter causes the less heat to transfer to the sheet. This may be the reason for decreasing sense in the temperature filed.
|
Figure 8. Velocity field for different values of magnetic field parameter.
|
|
Figure 9. Velocity field for different values of magnetic field parameter.
|
|
Figure 10. Temperature field for different values of magnetic field parameter.
|
|
Figure 11. Concentration field for different values of magnetic field parameter.
|
|
Figure 12. Velocity field for different values of an unsteadiness parameter.
|
|
Figure 13. Velocity field for different values of an unsteadiness parameter.
|
|
Figure 14. Temperature field for different values of an unsteadiness parameter.
|
|
Figure 15. Concentration field for different values of an unsteadiness parameter.
|
Figure 16, Figure 17, Figure 18 and Figure 19 present the effect of stretching ratio parameter on velocity, temperature and concentration fields for both Carreau and Casson fluid cases. The temperature field is suppressed and velocity and concentration fields are improved with increasing values of stretching ratio parameter. Generally, the stretching keeps more pressure on the sheet. Due to this reason we have seen a fall in temperature field and hike in velocity field. Figure 20 and Figure 21 demonstrate the effect of non-uniform heat source/sink parameter on temperature distribution of the flow for both Carreau and Casson fluid cases. It is clear that increasing values of space and temperature dependent heat source/sink parameters enhances the thermal boundary layer thickness of the flow over a stretching sheet for both Carreau fluid and Casson fluid cases. Physically, positive values of the non-uniform heat source/sink parameters act like a heat generator, which releases the heat energy to the flow and enhances the temperature profiles.
|
Figure 16. Velocity field for different values of stretching ratio parameter.
|
|
Figure 17. Velocity field for different values of stretching ratio parameter.
|
|
Figure 18. Temperature field for different values of stretching ratio parameter.
|
|
Figure 19. Concentration field for different values of stretching ratio parameter.
|
|
Figure 20. Temperature field for different values of non-uniform heat source/sink parameter.
|
|
Figure 21. Temperature field for different values of non-uniform heat source/sink parameter.
|
The strength of homogeneous and heterogeneous parameters on the concentration field is depicted in Figure 22 and Figure 23 for both Carreau and Casson fluids. It is clear that increasing values of homogeneous and heterogeneous parameters suppresses the concentration field for both cases. Physically, variation in the strength of homogeneous–heterogeneous reactions fluctuates the diffusivity of the flow. This causes to reduce the concentration profiles of the flow. Table 1 and Table 2 depict the validation of the present results by comparing with the existed literature under some special limited cases. We found a better agreement of the present results with the existed literature. This proves the validity of the present results along with the accuracy of the numerical technique we used in this study.
|
Figure 22. Concentration field for different values of strength of homogeneous reaction parameter.
|
|
Figure 23. Concentration field for different values of strength of heterogeneous reaction parameter.
|
β | M | ||||||
---|---|---|---|---|---|---|---|
Nadeem et al. [29] | Present results | Nadeem et al. [29] | Nadeem et al. [29] | Present results | Present results | ||
1 | 0 | −1.4142 | −1.4142 | −1.5459 | −0.6579 | −1.5459 | −0.6579 |
5 | −1.0954 | −1.0952 | −1.1974 | −0.5096 | −1.1980 | −0.5096 | |
−1.0049 | −1.0049 | −1.0932 | −0.4653 | −1.0932 | −0.4653 | ||
1 | 10 | −4.6904 | −4.6904 | −4.7263 | −2.3276 | −4.7263 | −2.3276 |
5 | −3.6331 | −3.6331 | −3.6610 | −1.8030 | −3.6610 | −1.8030 | |
−3.3165 | −3.3165 | −3.3420 | −1.6459 | −3.3420 | −1.6459 | ||
1 | 100 | −14.2127 | −14.212 | −14.2244 | −7.1004 | −14.224 | −7.1004 |
5 | −11.0091 | −11.009 | −11.0182 | −5.5000 | −11.017 | −5.4998 | |
−10.0490 | −10.049 | −10.058 | −5.0208 | −10.058 | −5.02079 |
Roy [28] | Raju et al.[20] | Present results | Roy [28] | Raju et al. [20] | Present results | |
---|---|---|---|---|---|---|
0 | 0.4305 | 0.4769 | 0.4305 | 0.557294 | 0.5180 | 0.5572 |
1 | 0.6127 | 0.6004 | 0.6004 | 0.721982 | 0.7005 | 0.7219 |
10 | 1.0175 | 1.0097 | 1.0172 | 1.170983 | 1.1494 | 1.1709 |
Table 3 and Table 4 display the variations in the friction factors, local Nusselt and Sherwood numbers for Carreau and Casson fluids for various values of non-dimensional governing parameters. It is noticed from the tables that the hike in the values of unsteadiness parameter enhances the friction factor coefficients and heat transfer rate of both Carreau and Casson fluids. We have seen exactly opposite results for increasing values of magnetic field parameter and Weissenberg number. Rise in the values of non-uniform heat source/sink parameter does not influence the friction factor and mass transfer rate but it reduces the Nusselt number for both Carreau and Casson fluid cases. We have observed a similar type of results for increasing values of thermal radiation parameter. Rise in the homogeneous–heterogeneous parameters shows a mixed response in mass transfer rate. This concludes that the homogeneous–heterogeneous reaction parameters help to control the concentration profiles of the flow.
Table 3.
Physical parameter values of for Casson fluid.
M | R | K | A | We | |||||||
---|---|---|---|---|---|---|---|---|---|---|---|
1 | −0.370542 | −0.218615 | 0.806053 | −0.018040 | |||||||
2 | −0.626189 | −0.324276 | 0.708147 | −0.014499 | |||||||
3 | −0.914227 | −0.427845 | 0.630938 | −0.012885 | |||||||
1 | −0.914227 | −0.427845 | 0.563670 | −0.012885 | |||||||
2 | −0.914227 | −0.427845 | 0.464686 | −0.012885 | |||||||
3 | −0.914227 | −0.427845 | 0.395295 | −0.012885 | |||||||
1 | −0.914227 | −0.427845 | 0.630938 | −0.012885 | |||||||
2 | −0.914227 | −0.427845 | 0.630933 | −0.004016 | |||||||
3 | −0.914227 | −0.427845 | 0.630933 | −0.001458 | |||||||
1 | −0.914227 | −0.427845 | 0.636191 | −0.079147 | |||||||
2 | −0.914227 | −0.427845 | 0.636191 | −0.094637 | |||||||
3 | −0.914227 | −0.427845 | 0.636191 | −0.101582 | |||||||
1 | −0.353005 | −0.268963 | 0.790376 | −0.017153 | |||||||
−0.022772 | 2 | −0.218155 | −0.194198 | 0.867563 | |||||||
3 | −0.172077 | −0.160872 | 0.901181 | −0.027919 | |||||||
1 | −5.929220 | −0.799438 | 0.401890 | −0.010760 | |||||||
3 | −8.094415 | −1.335176 | 0.363705 | −0.010528 | |||||||
5 | −10.63697 | −2.366287 | 0.328190 | −0.010329 | |||||||
0.2 | −1.356482 | −0.323332 | 0.638757 | −0.013330 | |||||||
0.4 | −8.613557 | −1.097922 | 1.079688 | −0.022118 | |||||||
0.6 | −32.08378 | −2.980574 | 1.402644 | −0.036072 | |||||||
1 | −0.368418 | −0.331556 | 0.647340 | −0.088065 | |||||||
2 | −0.368418 | −0.331556 | 0.546861 | −0.088065 | |||||||
3 | −0.368418 | −0.331556 | 0.446312 | −0.088065 |
M | R | K | A | We | |||||||
---|---|---|---|---|---|---|---|---|---|---|---|
1 | −0.234429 | −0.221416 | 0.848403 | −0.020765 | |||||||
2 | −0.354917 | −0.331612 | 0.767788 | −0.016204 | |||||||
3 | −0.478409 | −0.443038 | 0.701860 | −0.014134 | |||||||
1 | −0.478409 | −0.443038 | 0.630469 | −0.014134 | |||||||
2 | −0.478409 | −0.443038 | 0.525606 | −0.014134 | |||||||
3 | −0.478409 | −0.443038 | 0.452131 | −0.014134 | |||||||
1 | −0.478409 | −0.443038 | 0.701860 | −0.014134 | |||||||
2 | −0.478409 | −0.443038 | 0.701861 | −0.004388 | |||||||
3 | −0.478409 | −0.443038 | 0.701861 | −0.001587 | |||||||
1 | −0.499159 | −0.442978 | 0.699096 | −0.082839 | |||||||
2 | −0.499159 | −0.442978 | 0.699096 | −0.098872 | |||||||
3 | −0.499159 | −0.442978 | 0.699096 | −0.106066 | |||||||
1 | −0.281248 | −0.271720 | 0.812987 | −0.018332 | |||||||
2 | −0.198141 | −0.194860 | 0.875240 | −0.023661 | |||||||
3 | −0.162777 | −0.161132 | 0.904927 | −0.028765 | |||||||
1 | −0.880536 | −0.802796 | 0.580387 | −0.012285 | |||||||
3 | −1.473479 | −1.342562 | 0.536358 | −0.011855 | |||||||
5 | −2.595351 | −2.380762 | 0.494735 | −0.011495 | |||||||
0.2 | −0.354917 | −0.331612 | 0.767788 | −0.016204 | |||||||
0.4 | −1.409036 | −1.273276 | 1.216572 | −0.032431 | |||||||
0.6 | −5.205586 | −4.614127 | 1.535326 | −0.062058 | |||||||
1 | −1.356482 | −0.323332 | 0.545206 | −0.079966 | |||||||
2 | −1.356482 | −0.323332 | 0.494005 | −0.079966 | |||||||
3 | −1.356482 | −0.323332 | 0.442794 | −0.079966 |
In this study, we proposed a mathematical model for analyzing the effects of nonlinear thermal radiation on three-dimensional flow of Carreau and Casson fluids past a stretching surface with homogeneous–heterogeneous reactions and non-uniform heat source/sink. The nonlinear thermal radiation has importance in various industrial as well as science and technological applications. The transformed governing equations are solved numerically using Runge–Kutta based shooting technique. We presented dual solutions for the flow of Carreau and Casson fluids over a stretching sheet. The numerical findings are as follows:
Published on 12/04/17
Licence: Other
Are you one of the authors of this document?