Next Article in Journal
Effect of Localized Temperature Difference on Hydrogen Fermentation
Next Article in Special Issue
Estimation of the Wind Energy Potential in Various North Algerian Regions
Previous Article in Journal
Numerical Simulation of Combustion of Natural Gas Mixed with Hydrogen in Gas Boilers
Previous Article in Special Issue
Mechanism of Solute and Thermal Characteristics in a Casson Hybrid Nanofluid Based with Ethylene Glycol Influenced by Soret and Dufour Effects
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Using Non-Fourier’s Heat Flux and Non-Fick’s Mass Flux Theory in the Radiative and Chemically Reactive Flow of Powell–Eyring Fluid

1
Department of Science and Humanities, National University of Computer and Emerging Sciences, Lahore 44000, Pakistan
2
Section of Mathematics, International Telematic University Uninettuno, Corso Vittorio Emanuele II 39, 00186 Roma, Italy
3
Department of Statistics and Operations Research, College of Science, King Saud University, P.O. Box 2455, Riyadh 11451, Saudi Arabia
*
Authors to whom correspondence should be addressed.
Energies 2021, 14(21), 6882; https://doi.org/10.3390/en14216882
Submission received: 7 September 2021 / Revised: 6 October 2021 / Accepted: 8 October 2021 / Published: 20 October 2021
(This article belongs to the Special Issue Recent Advances in Solar Energy Collectors: Models and Applications)

Abstract

:
The behavior of convective boundary conditions is studied to delineate their role in heat and mass relegation in the presence of radiation, chemical reaction, and hydro-magnetic forces in three-dimensional Powell–Eyring nanofluids. Implications concerning non-Fourier’s heat flux and non-Fick’s mass flux with respect to temperature nanoparticle concentration were examined to discuss the graphical attributes of the principal parameters. An efficient optimal homotopy analysis method is used to solve the transformed partial differential equations. Tables and graphs are physically interpreted for significant parameters.

1. Introduction

Most of the fluid mechanics problems tacitly assume that the fluid is obeying the Fourier law [1] for heat transfer and Fick’s law for mass transfer. But it prognosticates the exhaustless speed of heat propagation. The frailty that heat disruption will be perceived immediately at other points of the medium is that attribute of Fourier constitutive law which ignores the principle of determinism in continuum mechanics. Cattaneo [2] proposed an answer to the unrealistic feature of the Fourier law by adding an intrinsic relaxation time. Although Cattaneo theory sustain the principle of determinism, it conflicts the Galilean postulate of frame-indifference. Therefore, Christov [3] remodel the heat transfer theory by replacing the Cattaneo theory on material time derivatives. The well known Catteneo-Christov (CC) model has attracted many researchers. Realizing the importance of CC theory, Straughan [4], [5] published his results with application to acoustic wave propagation and thermal convection. The flow of Maxwell fluid resulted by the inconstant stretching sheet incorporating CC model was analyzed by Hayat et al. [6]. The CC model has been investigated in the study of many nonlinear fluids [7,8,9,10,11].
Many biological and industrial operations involve complex fluids and their fragmentation. Non-Newtonian fluids have physical and worthy applications in coating flows, extrusion processes, mold filling and many others. The Eyring-Powell (EP) fluid model [12] fits to a notable class of non-linear fluid. The model is dominant to the other nonlinear models as it can depreciate to visco-fluids for immense and limited shear rates. It is viable in modeling pseudo-plastic attribute over different shear rates. Considering the importance of [12], Hayat et al. [13] analyzed the steady flow of [12]. The peristaltic flow of EP nanofluid was investigated by Akbar and Nadeem [14]. Boundary layer flow of [12] was examined numerically by Jalil et al. [15]. Malik et al. [16] also investigated the fluid flow characterized due to [12]. Mixed convected flow of [12] along a rotating cone was presented by Nadeem and Saleem [17]. The role of hybrid EP fluid [18] is observed for peristaltic transport. Hayat et al. [19] presented results for the axisymmetric radial flow of [12] over an impermeable stretching surface and the heat transfer process was analyzed through convective boundary conditions. Recently, Ibrahim [20] proposed the numerical solution for the rotating EP fluid flow in three dimensions with theory [3].
The dynamics of the flows that emerge from the conjoint interaction of fluid and magnetic fields is known as magnetohydrodynamics (MHD). The MHD study is highly motivated by its widespread application to the description of space and astrophysical plasmas, propelling liquid metals with magnets, magnetic resonance imaging for tumor diagnosis, the peristaltic activity of ureter scanned by giant magnetic resistive sensors and many more. Marsch [21] analyse the properties of MHD fluctuations in the interplanetary medium. An unsteady MHD boundary layer flow is examined numerically with viscous dissipation [22]. Maxwell nanofluid was discussed [23] with MHD effect.
Thermal radiation manifests a boundless demeanor of nature. Every material with a finite temperature emanates thermal radiation due to the motion of charged particles. Some promising applications of thermal radiation are in the field of physics and engineering disciplines such as boilers, radiative cooling, infrared sensing, and aeronautics. The heat transfer analysis with thermal radiation is significant in food processing and manufacturing industries as well. Realizing the importance of radiation Qasim et al. [24] carried out analysis of viscoelastic fluid with radiation effects. Mehmood and Fetecau [25] discussed the influence of thermal radiation on peristaltic transport of Sisko fluid. Siva and Govindarajan [26] documented the peristaltic transport under the influence of Soret effect and thermal radiation of hydromagnetic Newtonian fluid. A short time ago, Mallawi et al. [27] address the outcome of thermal radiation over a Riga plate with CC model.
Many of the modern propositions in technology are intended on making small devices. This can improve the efficiency scale and enhance the productivity. Similarly, advances are also occurring in fluid dynamics at a rapid pace known as micro fluidics and nanofluidic. Due to the enormous applications of nanofluids, many researchers have shown their interest in the studding effects of nanoparticles in non-Newtonian fluids with different physical aspects [28,29,30,31,32,33].
Above literature motivated us to target the analytical solutions of three-dimensional rotating EP fluid inclusive of MHD, radiation effects, and convective boundary conditions with [3]. Convergent series solution is devised by the optimal homotopy approach [33,34,35,36,37,38]. The impact of important parameters on velocity components, temperature, and concentration are interpreted graphically. At the same time, skin friction, Nusselt number, and Sherwood numbers are tabulated numerically.

2. Mathematical Formulation

The mathematical design of the non-Newtonian fluid called Powell–Eyring fluid is investigated. For this purpose, the stress tensor of the fluid is taken from [12].
D = P I + μ W + 1 β ξ ˙ sinh ( 1 d ξ ˙ ) W ,
The dynamic viscosity and Powell-Eyring fluid parameters are represented by   μ , β   a n d   d . ξ ˙ = 1 2 t r W 2   and W = V ¯ + ( V ¯ ) T so that the second-order approximation of sin h 1 after using Taylor series expansion is
s i n h 1 ( 1 d ξ ˙ ) ( 1 d ξ ˙ ) 1 6 ( 1 d ξ ˙ ) 3 , | 1 2 ξ ˙ | 1
Hence Equation (1) takes the form
D = P I + ( μ + 1 β d ) W 1 6 β d 3 ( ξ ˙ ) 3 W
The chemical reaction effect of MHD steady and incompressible EP nano fluid in the essence of thermal radiation over two-sided stretching sheet with [3] is explored. The coordinate system fixed in the x ¯ y ¯ -plane such that the extending velocities at the layer in both x ¯ and y ¯ directions are u ¯ = a x ¯ ,   v ¯ = b y ¯ , where a > 0 and b > 0 are constants. Magnetic field B ¯ is functional perpendicular to the space z ¯ 0 , as the fluid is flowing in that direction. Furthermore, the fluid is rotating with uniform angular velocity ʊ about   z ¯ -axis. The sheet temperature is constant and presumed to be higher than the outside temperature T ¯ w > T ¯ . C ¯ w is the concentration of the nanoparticles at the surface and C ¯ is given as the ambient concentration as shown in Figure 1. Under these assumptions, the boundary layer flow is driven by the following conservation laws of mass, momentum, and energy.
u ¯ x ¯ + v ¯ y ¯ + w ¯ z ¯ = 0 ,
u ¯ u ¯ x ¯ + v ¯ u ¯ y ¯ + w ¯ u ¯ z ¯ 2 v ¯ = ( υ + 1 ρ β d ) ( 2 u ¯ z ¯ 2 ) 1 2 ρ β d 3 ( u ¯ z ¯ ) 2 2 u ¯ z ¯ 2 σ B ¯ 2 ρ u ¯ ,
u ¯ v ¯ x ¯ + v ¯ v ¯ y ¯ + w ¯ v ¯ z ¯ + 2 u ¯ = ( υ + 1 ρ β d ) ( 2 v ¯ z ¯ 2 ) 1 2 ρ β d 3 ( v ¯ z ¯ ) 2 2 v ¯ z ¯ 2 σ B ¯ 2 ρ v ¯ ,
u ¯ T ¯ x ¯ + v T ¯ y ¯ + w T ¯ z ¯ = α t ( 2 T ¯ z ¯ 2 ) χ q 1 ρ c p q ¯ r z ¯ + Λ { D ¯ B ( C ¯ z ¯ T ¯ z ¯ ) + D ¯ T T ¯ ( T ¯ z ¯ ) 2 } ,
u ¯ C ¯ x ¯ + v ¯ C ¯ y ¯ + w ¯ C ¯ z ¯ = D ¯ B ( 2 C ¯ z ¯ 2 ) + D ¯ T T ¯ ( 2 T ¯ z ¯ 2 ) Q ( C ¯ C ¯ ) ,
Cattaneo–Christov theory is incorporated in place of classical [3]. The equations involving heat flux χ q and mass flux χ J are given.
χ q + λ q ( χ q t ¯ + V ¯ . χ q χ q . V ¯ + ( . V ¯ ) χ q ) = K T ¯
χ J + λ J ( χ J t ¯ + V ¯ . χ J χ J . V ¯ + ( . V ¯ ) χ J ) = D ¯ B C ¯
where λ q is the thermal relaxation time and λ J is the concentration relaxation time. To proceed further, we use the Rosseland approximation for the radiative heat flux q r .
q ¯ r = 4 σ s 3 K e T ¯ 4 z ¯
where σ s   is the Stefan-Boltzmann constant and K e is the mean absorption coefficient.
Expanding T ¯ 4 into Taylor series about T ¯ and neglecting higher order terms
T ¯ 4 = 4 T ¯ 3 T ¯ 3 T ¯ 4 .
In accordance with the above supposition, the heat and mass transfer equation will reduce to
u ¯ T ¯ x ¯ + v ¯ T ¯ y ¯ + w ¯ T ¯ z ¯ = λ q [ u ¯ 2 2 T ¯ x ¯ 2 + v ¯ 2 2 T ¯ y ¯ 2 + w ¯ 2 2 T ¯ z ¯ 2 + 2 u ¯ v ¯ 2 T ¯ x ¯ y ¯ + 2 v ¯ w ¯ 2 T ¯ y ¯ z ¯ + 2 u ¯ w ¯ 2 T ¯ x ¯ z ¯ + ( u ¯ u ¯ x ¯ + v ¯ u ¯ y ¯ + w ¯ u ¯ z ¯ ) T ¯ x ¯ + ( u ¯ v ¯ x ¯ + v ¯ v ¯ y ¯ + w ¯ v ¯ z ¯ ) T ¯ y ¯ + ( u ¯ w ¯ x ¯ + v ¯ w ¯ y ¯ + w ¯ w ¯ z ¯ ) T ¯ z ¯ ] + α t 2 T ¯ z ¯ 2 16 σ s T ¯ 3 3 K e ρ c p 2 T ¯ z ¯ 2 + Λ [ D ¯ B ( C ¯ z ¯ T ¯ z ¯ ) + D ¯ T T ¯ ( T ¯ z ¯ ) 2 ] ,
u ¯ C ¯ x ¯ + v ¯ C ¯ y ¯ + w ¯ C ¯ z ¯ = λ J [ u ¯ 2 2 C ¯ x ¯ 2 + v ¯ 2 2 C ¯ y ¯ 2 + w ¯ 2 2 C ¯ z ¯ 2 + 2 u ¯ v ¯ 2 C ¯ x ¯ y ¯ + 2 v ¯ w ¯ 2 C ¯ y ¯ z ¯ + 2 u ¯ w ¯ 2 C ¯ x ¯ z ¯ + ( u ¯ u ¯ x ¯ + v ¯ u ¯ y ¯ + w ¯ u ¯ z ¯ ) C ¯ x ¯ + ( u ¯ v ¯ x ¯ + v ¯ v ¯ y ¯ + w ¯ v ¯ z ¯ ) C ¯ y ¯ + ( u ¯ w ¯ x ¯ + v ¯ w ¯ y ¯ + w ¯ w ¯ z ¯ ) C ¯ z ¯ ] + D ¯ B 2 C ¯ z ¯ 2 + D ¯ T T ¯ ( 2 T ¯ z ¯ 2 ) Q ( C ¯ C ¯ ) .
The boundary conditions associated with the study are
u ¯ = u ¯ ( x ¯ ) = a x ¯ , v ¯ = v ¯ ( y ¯ ) = b y ¯ , w ¯ = 0 , K T ¯ z ¯ = H 1 ( T ¯ w T ¯ ) , D ¯ B C ¯ z ¯ = H 2 ( C ¯ w C ¯ ) a t z ¯ = 0 u ¯ 0 , v ¯ 0 , T ¯ T ¯ , C ¯ C ¯ a s z ¯ .
Selecting the following similarity transformations
η = a υ z ¯ , u ¯ = a x ¯ f ¯ ( η ) , v ¯ = a y ¯ g ¯ ( η ) , w ¯ = a υ ( f ¯ ( η ) + g ¯ ( η ) ) , θ ¯ ( η ) = T ¯ T ¯ T ¯ w T ¯ , φ ¯ ( η ) = C ¯ C ¯ C ¯ w C ¯ .
The nondimensional equations in view of aforesaid conditions can be written as
( 1 + n ¯ ) f ¯ n ¯ e 1 ( f ¯ ) 2 f ¯ f ¯ 2 + ( f ¯ + g ¯ ) f ¯ + 2 γ Ω g ¯ M f ¯ = 0 ,
( 1 + n ¯ ) g ¯ n ¯ e 2 ( g ¯ ) 2 g ¯ g ¯ 2 + ( f ¯ + g ¯ ) g ¯ + 2 Ω γ f ¯ M g ¯ = 0 ,
( 1 + 4 3 N ¯ 1 ) θ ¯ + P r ( N ¯ b φ ¯ θ ¯ + N ¯ t θ ¯ 2 + ( f ¯ + g ¯ ) θ ¯ δ t ( f ¯ + g ¯ ) ( f ¯ + g ¯ ) θ ¯ + ( f ¯ + g ¯ ) 2 θ ¯ ) = 0 ,
φ ¯ + ( N ¯ t N ¯ b ) θ ¯ + S c ( ( f ¯ + g ¯ ) φ ¯ δ c ( f ¯ + g ¯ ) ( f ¯ + g ¯ ) φ ¯ + ( f ¯ + g ¯ ) 2 φ ¯ ε φ ¯ ) = 0 .
The corresponding boundary conditions are
f ¯ ( 0 ) = g ¯ ( 0 ) = 0 , f ¯ ( 0 ) = 1 , g ¯ ( 0 ) = A , θ ¯ ( 0 ) = α 1 ( 1 θ ¯ ( 0 ) ) , φ ¯ ( 0 ) = α 2 ( 1 φ ¯ ( 0 ) ) a t η = 0 . f ¯ ( ) 0 , g ¯ ( ) 0 , θ ¯ ( ) 0 , φ ¯ ( ) 0 a s η .
Here n ¯ = 1 μ β d ,   e 1 = a 3 x ¯ 2 2 υ d 2 ,   e 2 = b 3 y ¯ 2 2 υ d 2 are Eyring-Powell fluid parameters. A = b a is the stretching ratio parameter. = a is the rotation parameter and γ = y x ¯ ¯ . M 2 = σ B ¯ 2 a ρ is the magnetic parameter. ϵ = Q a is chemical reaction parameter. δ t = a λ q ,     δ c = a λ j denotes the nondimensional thermal and concentration relaxation parameter. N ¯ 1 = 4 σ s T ¯ 3 α t K e ρ c p represent thermal radiation parameter. P r = υ α t is Prandtl number. N ¯ b = Λ D ¯ B ( C ¯ w C ¯ ) υ is Brownian motion parameter, N ¯ t = Λ D ¯ T ( T ¯ w T ¯ ) υ T ¯ is thermophoresis parameter. α 1 = H 1 K υ a and α 2 = H 2 D ¯ B υ a are the heat transfer Biot number and mass transfer Biot number. The engineering components of skin friction coefficient C f , local Nusselt number N u x and the local Sherwood number S h x are represented as follows:
C f x = τ w x ρ u ¯ w x 2 , C f Y = τ w y ρ v ¯ w y 2 , N u x = x ¯ q ¯ w K ( T ¯ w T ¯ ) , S h x = x ¯ h ¯ w D ¯ B ( C ¯ w C ¯ ) ,
where q ¯ w surface heat flux, h ¯ w surface mass flux, wall shear stress along x ¯ -axis and y ¯ -axis are given by τ w x and τ w y .
τ w x = [ μ ( u ¯ z ¯ ) + 1 β d ( u ¯ z ¯ ) 1 6 β d 3 ( u ¯ z ¯ ) 3 ] z ¯ = 0 , τ w y = [ μ ( v ¯ z ¯ ) + 1 β d ( v ¯ z ¯ ) 1 6 β d 3 ( v ¯ z ¯ ) 3 ] z ¯ = 0 , h ¯ w = D ¯ B ( C ¯ z ¯ ) z ¯ = 0 , q ¯ w = K ( C ¯ z ¯ ) z ¯ = 0
By incorporating the above equations, we get
C f x R e x = ( 1 + n ¯ ) f ¯ ( 0 ) 1 3 n ¯ e 1 f ¯ 3 ( 0 ) , C f y R e y = ( 1 + n ¯ ) g ¯ ( 0 ) 1 3 n ¯ e 2 g ¯ 3 ( 0 ) , N u x R e x = θ ¯ ( 0 ) , S h x R e x = φ ¯ ( 0 ) .
R e x ,   R e y are local Reynolds numbers.

3. Method of Solution

Optimal homotopy method [33,34,35,36,37] is adapted to obtain the solutions for the nonlinear Equations (16)–(19) together with boundary conditions. The initial guesses are given as
f ¯ 0 = 1 e η , g ¯ 0 = A ( 1 e η ) , θ ¯ 0 = ( α 1 1 + α 1 ) e η , φ ¯ 0 = ( α 2 1 + α 2 ) e η .
ψ f ¯ ( η ) = f ¯ f ¯ , ψ g ¯ ( η ) = g ¯ g ¯ , ψ θ ¯ ( η ) = θ ¯ θ ¯ , ψ φ ¯ ( η ) = φ ¯ φ ¯ .
The concept of minimizing the average square residual errors is utilized [33] to find the ideal values of nonzero auxiliary parameters h 0 f ,   h 0 g ,   h 0 θ , and h 0 φ which are actually responsible for defining the convergence region of homotopy series solutions.
E m f ¯ = 1 k + 1 j = 0 k { N f ¯ ( i = 0 m f ^ ( η ) , i = 0 m g ^ ( η ) , i = 0 m θ ^ ( η ) , i = 0 m φ ^ ( η ) ) η = j δ η } 2 d η ,
E m g ¯ = 1 k + 1 j = 0 k { N g ¯ ( i = 0 m f ^ ( η ) , i = 0 m g ^ ( η ) , i = 0 m θ ^ ( η ) , i = 0 m φ ^ ( η ) ) η = j δ η } 2 d η ,
E m θ ¯ = 1 k + 1 j = 0 k { N θ ¯ ( i = 0 m f ^ ( η ) , i = 0 m θ ^ ( η ) , i = 0 m φ ^ ( η ) ) η = j δ η } 2 d η .
E m φ ¯ = 1 k + 1 j = 0 k { N φ ¯ ( i = 0 m f ^ ( η ) , i = 0 m θ ^ ( η ) , i = 0 m φ ^ ( η ) ) η = j δ η } 2 d η .
E m t ¯ = E m f ¯ + E m g ¯ + E m θ ¯ + E m φ ¯
where E m t ¯ is the total of the square of residual error, δ η = 0.5 ,   k = 20 . Mathematica package BVPh2.0 has been used to reduce the average residual error. Table 1 is aligned to show the minimized values of the total residual error at various iterations.
In Table 2, the residual errors for f ¯ ,   g ¯ ,   θ ¯ ,   φ ¯   are given at three distinct iterations with the eighth-order optimal convergence control parameters. It is evident that the residual errors are reduced by raising iterations. Therefore, OHAM provides a procedure to select any set of local convergence control parameters to find convergent outcomes.

4. Effect of Parameters on Powell–Eyring Fluid Flow

The plots of velocities along x ¯ and y ¯ direction with coordinate η are executed for the important parameters involved in the flow.
Figure 2 and Figure 3 depict the decrease in velocity along both directions as the rotation in the flow increases. Velocity along y ¯ -direction is diminuting more compared to x ¯ -direction, this may be because of the retarding force which has more impact along y ¯ -direction. Figure 4 and Figure 5 narrates the behavior of Powell-Eyring fluid parameter n ¯   on the velocities. Velocity ascends by increasing non-Newtonian factors along x ¯ -direction, while along y ¯ -axis it increases near the surface and descends away from surface. Dimensionless velocities along x ¯ and y ¯ increase with the increasing ratio γ , as directed in Figure 6 and Figure 7. One of the most prime parameters to entertain for flow behavior is the magnetic parameter M . Velocity parallel to y ¯ -axis escalates by raising magnetic effects, but the reverse behavior is observed when parallel to x ¯ -axis as bent upon by Figure 8 and Figure 9. In Table 3, the skin friction coefficient on the surface is measured up for various values of n ¯ and e 1 . We noted that the local skin friction coefficient along x ¯ -axis is reduced for sufficiently large values of n ¯ , hence the smooth flow along x ¯ -direction and reverse is the behavior of skin friction co-efficient if e 1 is raised. Table 4 is given to analyze the skin friction co-efficient in absence of MHD, radiation, and chemical reaction. Skin friction coefficient along y ¯ -direction abates for leading values of e 2 . In accrual, the skin friction co-efficient exceeds for greater values of n ¯ ( e 2 0 )   and recedes for greater values of n ¯ ( e 2 > 0 ) as indicated in Table 5. Table 6 is constructed for the special behavior in the absence of MHD, radiation, and chemical reaction.

5. Effects of Parameters on Temperature and Concentration

The plots for temperature with coordinate η are executed for the important parameters involved in the heat/mass transfer.
Figure 10 is plotted for the analysis of parameter N ¯ 1 on θ ¯ and φ ¯ fields. It can be seen that the temperature of the fluid positively increases as the radiation parameters vary. Further, Figure 11 shows that the volume fraction boundary layer thickens near the sheet, while the reverse effect is noticed away from the sheet with the rising values of radiation parameter. Temperature profile drops with a raise in the magnetic field strength and the concentration field also decline with booming magnetic parameter as illustrated in Figure 12 and Figure 13. Chemical reactions restrained the thermal and concentration fields (Figure 14 and Figure 15). Figure 16 demonstrates typical θ ¯ and φ ¯ profiles for diverse heat transfer Biot number α 1 . The bold curves show the temperature profile, while the dotted curves mark the concentration profile. We can see that an increase in Biot number depicts an increase in the heat transfer coefficient which is responsible for an increase in temperature and thermal boundary layer thickness. Low heat transfer rate is seen for smaller values of heat transfer coefficient, when it gradually increases, the fluid becomes hotter and has a high rate of heat transfer. This concrete behavior can be seen from Figure 17, that its the concentration of nano fluid accumulated due to an increase in the mass transfer Biot number α 2 . Thermophoresis parameter increases the temperature and concentration profiles while Brownian motion parameter enhances the temperature profile, while the concentration profile decreases with an increase in Brownian motion parameter as explained in Figure 18 and Figure 19. Scrutinizing Table 7, we can see that the heat transfer and mass transfer are enhanced significantly by strengthening the magnetic effect. The same outcomes can be observed by mounting the radiation parameter, heat and mass transfer Biot numbers.

6. Conclusions

For this study of chemically reactive flow of Powell–Eyring nanofluid using non-Fourier’s heat flux and non-Fick’s mass flux theory with radiation, the following developments can be compiled as:
The amplitude of velocity along the x ¯ -direction decreases, while the velocity along the y ¯ -direction increases with the escalating magnetic effect. Temperature rises with the consideration of nanoparticles in the fluid, but the opposite is true for the concentration profile.
Rotational parameter decreases the velocity boundary layer thickness in the directions.
The Powell–Eyring fluid parameter n ¯ has the property of reducing the Nusselt number in both directions.
Heat transfer rate—i.e., the Nusselt number increases as rotation parameter —began to increase. Same response can be seen for increasing magnetic parameter M .
Mass transfer rate—i.e., the Sherwood number S h —remains constant for different values of magnetic parameter M , while it rises for chemical reaction parameter ϵ .

Author Contributions

Conceptualization, H.F.; Data collection, H.F. and H.A.; Formal analysis, H.F. and S.A.; Investigation, H.A.; Methodology, H.F.; Software, S.T.S. and H.F.; Writing—original draft, H.F. and S.A.; Review, S.T.S.; Funding, H.A.; Writing—review and editing, H.A. and S.T.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research has been supported by Research Supporting Project number (RSP-2021/167), King Saud University, Riyadh, Saudi Arabia.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All the data available in this manuscript.

Acknowledgments

Research Supporting Project number (RSP-2021/167), King Saud University, Riyadh, Saudi Arabia

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

D Extra   stress   tensor
u ¯ ,   v ¯ Velocities   along   dimensional   coordinate   axis   ( x ¯   and   y ¯ -   direction )
T ¯ Fluid   temperature   inside   the   boundary   layer
η Dimensionless   similarity   variable   ( distance )
θ ¯ Temperature   function   without   dimensions  
Rotation   parameter  
D ¯ B , D ¯ T Brownian   and   thermophoresis   diffusion   coefficient
S c Schmidt   number
δ t Parameter   responsible   for   thermal   relaxation   time
δ c Parameter   responsible   for   concentration   relaxation   time

References

  1. Fourier, J.B.J. Théorie Analytique De La Chaleur; F. Didot: Paris, France, 1822. [Google Scholar]
  2. Cattaneo, C. Sulla conduzione del calore. Atti Sem. Mat. Fis. Univ. Modena 1948, 3, 83–101. [Google Scholar]
  3. Christov, C. On frame indifferent formulation of the Maxwell–Cattaneo model of finite-speed heat conduction. Mech. Res Commun. 2009, 36, 481–486. [Google Scholar] [CrossRef]
  4. Straughan, B. Acoustic waves in a Cattaneo-Christov gas. Phys. Lett. A. 2010, 26, 2667–2669. [Google Scholar] [CrossRef]
  5. Straughan, B. Thermal convection with the Cattaneo–Christov model. Int. J. Heat Mass Transf. 2010, 53, 95–98. [Google Scholar] [CrossRef]
  6. Hayat, T.; Farooq, M.; Alsaedi, A.; Al-Solamy, F. Impact of Cattaneo-Christov heat flux in the flow over a stretching sheet with variable thickness. AIP Adv. 2015, 5, 087159. [Google Scholar] [CrossRef]
  7. Salahuddin, T.; Malik, M.Y.; Hussain, A.; Bilal, S.; Awais, M. Mhd flow of Cattanneo–Christov heat flux model for williamson fluid over a stretching sheet with variable thickness:using numerical approach. J. Magn. Magn. Mater. 2016, 401, 991–997. [Google Scholar] [CrossRef]
  8. Hayat, T.; Nadeem, S. Flow of 3D Eyring-Powell fluid by utilizing Cattaneo-Christov heat flux model and chemical processes over an exponentially stretching surface. Results Phys. 2018, 8, 397–403. [Google Scholar] [CrossRef]
  9. Abbasi, F.M. On Cattaneo-Christov heat flux model for Carreau fluid flow over a slendering sheet. Results Phys. 2017, 7, 310–319. [Google Scholar] [CrossRef]
  10. Abbasi, F.M.; Shehzad, S.A.; Hayat, T.; Alsaedi, A.; Hegazy, A. Influence of Cattaneo-Christov heat flux in flow of an Oldroyd-B fluid with variable thermal conductivity. Int. J. Numer Methods Heat Fluid Flow. 2016, 26, 2271–2282. [Google Scholar] [CrossRef]
  11. Hayat, T.; Muhammad, T.; Alsaedi, A.; Mustafa, M. A comparative study for flow of viscoelastic fluids with Cattaneo-Christov heat flux. PLoS ONE 2016, 11, e0155185. [Google Scholar] [CrossRef] [Green Version]
  12. Powell, R.E.; Eyring, E. Mechanism for the Relaxation Theory of Viscosity. Nature 1944, 154, 427–428. [Google Scholar] [CrossRef]
  13. Hayat, T.; Iqbal, Z.; Qasim, M.; Obaidat, S. Steady flow of an Eyring Powell fluid over a moving surface with convective boundary conditions. Int. J. Heat Mass Transf. 2012, 55, 1817–1822. [Google Scholar] [CrossRef]
  14. Akbar, N.S.; Nadeem, S. Characteristics of heating scheme and mass transfer on the peristaltic flow for an Eyring–Powell fluid in an endoscope. Int. J. Heat Mass Transf. 2012, 55, 375–383. [Google Scholar] [CrossRef]
  15. Jalil, M.; Asghar, S. Flow and heat transfer of Powell–Eyring Fluid over a stretching surface: A Lie Group Analysis. J. Fluids Eng. 2013, 135, 121201. [Google Scholar] [CrossRef]
  16. Malik, M.Y.; Hussain, A.; Nadeem, S. Boundary layer flow of an Eyring-Powell model fluid due to a stretching cylinder with variable viscosity. Sci. Iran. 2013, 20, 313–321. [Google Scholar]
  17. Nadeem, S.; Saleem, S. Mixed convection flow of Eyring–Powell fluid along a rotating cone. Results Phys. 2014, 4, 54–62. [Google Scholar] [CrossRef] [Green Version]
  18. Riaz, A.; Ellahi, R.; Sait, M.S. Role of hybrid nanoparticles in thermal performance of peristaltic flow of Eyring-Powell fluid model. J. Therm. Anal. Calorim. 2020. [Google Scholar] [CrossRef]
  19. Hayat, T.; Makhdoom, S.; Awais, M.; Saleem, S.; Rashidi, M.M. Axisymmetric Powell-Eyring fluid flow with convective boundary condition: Optimal analysis. Appl. Math Mech. 2016, 37, 919–928. [Google Scholar] [CrossRef]
  20. Ibrahim, W. Three dimensional rotating flow of Powell-Eyring nanofluid with non-Fourier’s heat flux and non-Fick’s mass flux theory. Results Phys. 2018, 8, 569–577. [Google Scholar] [CrossRef]
  21. Marsch, E.; Tu, C.Y. Intermittency, non-Gaussian statistics and fractal scaling of MHD fluctuations in the solar wind. Nonlin. Process. Geophys. 1997, 4, 101–124. [Google Scholar] [CrossRef] [Green Version]
  22. Muhammad, Q.; Afridi, M.I.; Wakif, A.; Thoi, T.N.; Hussanan, A. Second Law Analysis of Unsteady MHD Viscous Flow over a Horizontal Stretching Sheet Heated Non-Uniformly in the Presence of Ohmic Heating: Utilization of Gear-Generalized Differential Quadrature Method. Entropy 2019, 21, 240. [Google Scholar] [CrossRef] [Green Version]
  23. Nadeem, S.; Akhtar, S.; Abbas, N. Heat transfer of Maxwell base fluid flow of nanomaterial with MHD over a vertical moving surface. Alex. Eng. J. 2020, 59, 1847–1856. [Google Scholar] [CrossRef]
  24. Qasim, M.; Hayat, T.; Obaidat, S. Radiation effect on the mixed convection flow of a viscoelastic fluid along an inclined stretching sheet. Z. Naturforsch. A 2012, 67, 195–202. [Google Scholar] [CrossRef]
  25. Mehmood, O.U.; Fetecau, C. A note on radiative heat transfer to peristaltic flow of Sisko fluid. Appl. Bionics and Biomech. 2015, 2015. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Siva, E.P.; Govindarajan, A. Thermal radiation and Soret effect on MHD peristaltic transport through a tapered asymmetric channel with convective boundary conditions. GJPAM 2016, 12, 213–221. [Google Scholar]
  27. Mallawi, F.O.M.; Bhuvaneswari, M.; Sivasankaran, S.; Eswaramoorthi, S. Impact of double-stratification on convective flow of a non-Newtonian liquid in a Riga plate with Cattaneo-Christov double-flux and thermal radiation. Ain Shams Eng. J. 2021, 12, 969–981. [Google Scholar] [CrossRef]
  28. Buongiorno, J. Convective transport in nanofluids. ASME J. Heat Transf. 2005, 128, 240–250. [Google Scholar] [CrossRef]
  29. Kuznetsov, A.V.; Nield, D.A. Natural convective boundary-layer flow of a nanofluid past a vertical plate. Int. J. Therm. Sci. 2010, 49, 243–247. [Google Scholar] [CrossRef]
  30. Khan, W.A.; Pop, I. Boundary-layer flow of a nanofluid past a stretching sheet. Int. J. Heat Mass Transf. 2010, 53, 2477–2483. [Google Scholar] [CrossRef]
  31. Nadeem, S.; Haq, R.; Khan, Z.H. Numerical solution of non-newtonian nanofluid flow over a stretching sheet. Appl. Nano Sci. 2014, 4, 625–631. [Google Scholar] [CrossRef] [Green Version]
  32. Ahmed, Z.; Nadeem, S.; Saleem, S.; Ellahi, R. Numerical study of unsteady flow and heat transfer CNT-based MHD nanofluid with variable viscosity over a permeable shrinking surface. Int. J. Numer. Method H. 2019, 29, 4607–4623. [Google Scholar] [CrossRef]
  33. Alvi, N.; Latif, T.; Qamar, H.; Asghar, S. Peristalsis of nanofluid with temperature dependent viscosity. J. Comp. Theor. Nano Sci. 2017, 14, 1417–1430. [Google Scholar] [CrossRef]
  34. Liao, S.J. An optimal homotopy-analysis approach for strongly nonlinear differential equations. Commun. Nonlin. Sci. Num. Sim. 2010, 15, 2003–2016. [Google Scholar] [CrossRef]
  35. Nadeem, S.; Mehmood, R.; Akbar, N.S. Optimized analytical solution for oblique flow of a Casson-nano fluid with convective boundary conditions. Int. J. Therm. Sci. 2014, 78, 90–100. [Google Scholar] [CrossRef]
  36. Saleem, S.; Firdous, H.; Nadeem, S.; Khan, A.U. Convective heat and mass transfer in magneto Walter’s B Nanofluid flow induced by a rotating cone. Arab. J. Sci. Eng. 2019, 44, 1515–1523. [Google Scholar] [CrossRef]
  37. Saleem, S.; Awais, M.; Nadeem, S.; Mustafa, T.; Sandeep, S. Theoretical analysis of upper-convected Maxwell fluid flow with Cattaneo-Christov heat flux model. Chin. J. Phys. 2017, 55, 1615–1625. [Google Scholar] [CrossRef]
  38. Awais, M.; Saleem, S.; Hayat, T.; Irum, S. Hydromagnetic couple-stress nanofluid flow over a moving convective wall: OHAM analysis. Acta Astronaut. 2016, 129, 271–276. [Google Scholar] [CrossRef]
Figure 1. Three-dimensional bi-directional stretching sheet of Powell–Eyring fluid flow.
Figure 1. Three-dimensional bi-directional stretching sheet of Powell–Eyring fluid flow.
Energies 14 06882 g001
Figure 2. Velocity along x ¯ -direction for variant Ω .
Figure 2. Velocity along x ¯ -direction for variant Ω .
Energies 14 06882 g002
Figure 3. Velocity along x ¯ -direction for variant Ω .
Figure 3. Velocity along x ¯ -direction for variant Ω .
Energies 14 06882 g003
Figure 4. Velocity along x ¯ -direction for variant n ¯ .
Figure 4. Velocity along x ¯ -direction for variant n ¯ .
Energies 14 06882 g004
Figure 5. Velocity along y ¯ -direction for variant n ¯ .
Figure 5. Velocity along y ¯ -direction for variant n ¯ .
Energies 14 06882 g005
Figure 6. Velocity along x ¯ -direction for variant γ .
Figure 6. Velocity along x ¯ -direction for variant γ .
Energies 14 06882 g006
Figure 7. Velocity along y ¯ -direction for variant γ .
Figure 7. Velocity along y ¯ -direction for variant γ .
Energies 14 06882 g007
Figure 8. Velocity along x ¯ -direction for variant M .
Figure 8. Velocity along x ¯ -direction for variant M .
Energies 14 06882 g008
Figure 9. Velocity along x ¯ -direction for variant M .
Figure 9. Velocity along x ¯ -direction for variant M .
Energies 14 06882 g009
Figure 10. Temperature profile for variant N ¯ 1 .
Figure 10. Temperature profile for variant N ¯ 1 .
Energies 14 06882 g010
Figure 11. Concentration profile for variant N ¯ 1 .
Figure 11. Concentration profile for variant N ¯ 1 .
Energies 14 06882 g011
Figure 12. Temperature profile for variant M .
Figure 12. Temperature profile for variant M .
Energies 14 06882 g012
Figure 13. Mass Concentration profile for variant M .
Figure 13. Mass Concentration profile for variant M .
Energies 14 06882 g013
Figure 14. Temperature profile for variant ϵ .
Figure 14. Temperature profile for variant ϵ .
Energies 14 06882 g014
Figure 15. Mass concentration profile for variant ϵ .
Figure 15. Mass concentration profile for variant ϵ .
Energies 14 06882 g015
Figure 16. Temperature and concentration profiles for variant   α 1 .
Figure 16. Temperature and concentration profiles for variant   α 1 .
Energies 14 06882 g016
Figure 17. Temperature and concentration profiles for variant   α 2 .
Figure 17. Temperature and concentration profiles for variant   α 2 .
Energies 14 06882 g017
Figure 18. Temperature and concentration profiles for variant N ¯ t .
Figure 18. Temperature and concentration profiles for variant N ¯ t .
Energies 14 06882 g018
Figure 19. Temperature and concentration profiles for variant N ¯ b . .
Figure 19. Temperature and concentration profiles for variant N ¯ b . .
Energies 14 06882 g019
Table 1. Optimal convergence control parameter values and total averaged squared residual errors
Table 1. Optimal convergence control parameter values and total averaged squared residual errors
Order m h 0 f ¯ h 0 g ¯ h 0 θ ¯ h 0 φ ¯ E m t ¯ CPU Time(sec)
2−0.58−0.87−1.90−1.871.89   × 10   4 14.47
4−0.53−1.01−1.85−1.842.73   × 10   5 312.73
6−0.47−1.08−1.85−1.831.10   × 10   5 4452.7
8−0.43−1.12−1.85−1.846.63   × 10   6 33756.6
Table 2. Individual averaged squared residual errors with optimized values at m = 8 from Table 1
Table 2. Individual averaged squared residual errors with optimized values at m = 8 from Table 1
Order m E m f ¯ E m g ¯ E m θ ¯ E m φ ¯ CPU Time(sec)
41.21   × 10   5 1.67   × 10   5 7.46   × 10   6 6.13   × 10   6 171.38
88.23   × 10   8 4.53   × 10   7 3.77   × 10   6 2.32   × 10   6 314.38
123.91   × 10   9 4.10   × 10   8 2.32   × 10   6 1.07   × 10   8 506.22
Table 3. Values for C f x R e x , when e 2 = 0.1 ,   γ = 0.1 ,   δ = 0.2 ,   α 1 ,   α 2 = 0.1 ,   M = 0.1 ,   Ω = 0.5 ,   P r = 0.7 ,   S c = 0.2 .
Table 3. Values for C f x R e x , when e 2 = 0.1 ,   γ = 0.1 ,   δ = 0.2 ,   α 1 ,   α 2 = 0.1 ,   M = 0.1 ,   Ω = 0.5 ,   P r = 0.7 ,   S c = 0.2 .
n ¯ 0.10.20.30.40.5
e 1
00.3790−0.2230−1.0345−2.0887−3.4222
0.10.4052−0.1505−0.8864−1.8077−2.8865
0.20.4309−0.0803−0.7474−1.5610−2.4658
0.30.4561−0.0123−0.6158−1.3391−2.1205
0.40.48070.0533−0.4907−1.1355−1.8237
0.50.50480.1168−0.3714−0.9456−1.5578
Table 4. Values for C f x R e x , when   e 2 = 0.1 ,   γ = 0.1 ,   δ = 0.2 ,   α 1 ,   α 2 = 0.1 ,   M = 0 ,   N ¯   1 = 0 ,   Ω = 0.5 ,   P r = 0.7 ,   S c = 0.2 .
Table 4. Values for C f x R e x , when   e 2 = 0.1 ,   γ = 0.1 ,   δ = 0.2 ,   α 1 ,   α 2 = 0.1 ,   M = 0 ,   N ¯   1 = 0 ,   Ω = 0.5 ,   P r = 0.7 ,   S c = 0.2 .
n ¯ 0.10.20.30.40.5
e 1
0−0.2609−0.9519−1.8591−3.0163−4.4604
0.1−0.2304−0.8671−1.6750−2.6393−3.7018
0.2−0.2004−0.7862−1.5090−2.3293−3.1479
0.3−0.1709−0.7086−1.3572−2.0674−2.7315
0.4−0.1419−0.6339−1.2168−1.8405−2.4049
0.5−0.1134−0.5620−1.0854−1.6381−2.1353
Table 5. Values for C f y R e y , when e 1 = 0.1 ,   γ = 0.1 ,   δ = 0.2 ,   α 1 ,   α 2 = 0.1 ,   M = 0.1 ,   Ω = 0.5   ,   P r = 0.7 ,   S c = 0.2 .
Table 5. Values for C f y R e y , when e 1 = 0.1 ,   γ = 0.1 ,   δ = 0.2 ,   α 1 ,   α 2 = 0.1 ,   M = 0.1 ,   Ω = 0.5   ,   P r = 0.7 ,   S c = 0.2 .
n ¯ 0.10.20.30.40.5
e 2
055.029762.85671.314680.447690.2985
0.1−360.339−886.582−1558.52−2409.23−3478.23
0.2−771.59−1818.44−3146.07−4818.55−6912.48
0.3−1178.74−2732.87−4691.9−7148.79−10214.9
0.4−1581.82−3630.03−6196.52−9401.23−13388.0
0.5−1980.84−4510.1−7660.48−11577.1−16434.2
Table 6. Values for C f y R e y , when   e 1 = 0.1 ,   γ = 0.1 ,   δ = 0.2 ,   α 1 ,   α 2 = 0.1 ,   M = 0 ,   N ¯ 1 = 0 ,   Ω = 0.5 ,   P r = 0.7 ,   S c = 0.2 .
Table 6. Values for C f y R e y , when   e 1 = 0.1 ,   γ = 0.1 ,   δ = 0.2 ,   α 1 ,   α 2 = 0.1 ,   M = 0 ,   N ¯ 1 = 0 ,   Ω = 0.5 ,   P r = 0.7 ,   S c = 0.2 .
n ¯ 0.10.20.30.40.5
e 2
052.623160.182168.364177.21186.766
0.1−310.516−772.798−1366.53−2122.21−3076.35
0.2−669.87−1589.57−2762.36−4247.19−6114.77
0.3−1025.46−2390.28−4119.65−6298.98−9030.92
0.4−1377.29−3175.1−5438.91−8278.8−11827.2
0.5−1725.4−3944.16−6720.67−10187.9−14506.0
Table 7. Computed values of heat transfer rate and mass transfer rate for different values of magnetic parameter   M , radiation parameter N ¯ 1 , chemical reaction parameter ϵ , heat transfer biot number α 1 ,   mass transfer biot number   α 2 .
Table 7. Computed values of heat transfer rate and mass transfer rate for different values of magnetic parameter   M , radiation parameter N ¯ 1 , chemical reaction parameter ϵ , heat transfer biot number α 1 ,   mass transfer biot number   α 2 .
M N ¯ 1 ϵ α 1 α 2 θ ¯ ( 0 ) φ ¯ ( 0 )
0.10.100.10.1 0.239077
0.5 0.249323
1.0 0.260205
1.5 0.271679
2.0 0.283674
0 0.1 0.1377120.14419
0.5 0.1442310.145819
1.0 0.1516940.147705
1.5 0.1601440.149863
2.0 0.1696280.152304
0.1 00.500.722521
1.001.44597
2.002.48147
1.00.52.464061.59123
1.02.529212.86485
2.02.602054.53856
2.00.54.629062.03578
1.04.776893.58565
2.04.940915.5765
0 0.10.10.3654090.234357
0.5 0.4273310.277738
1.0 0.587680.354092
1.5 0.9457520.481151
2.0 1.662870.68243
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Firdous, H.; Saeed, S.T.; Ahmad, H.; Askar, S. Using Non-Fourier’s Heat Flux and Non-Fick’s Mass Flux Theory in the Radiative and Chemically Reactive Flow of Powell–Eyring Fluid. Energies 2021, 14, 6882. https://doi.org/10.3390/en14216882

AMA Style

Firdous H, Saeed ST, Ahmad H, Askar S. Using Non-Fourier’s Heat Flux and Non-Fick’s Mass Flux Theory in the Radiative and Chemically Reactive Flow of Powell–Eyring Fluid. Energies. 2021; 14(21):6882. https://doi.org/10.3390/en14216882

Chicago/Turabian Style

Firdous, Hina, Syed Tauseef Saeed, Hijaz Ahmad, and Sameh Askar. 2021. "Using Non-Fourier’s Heat Flux and Non-Fick’s Mass Flux Theory in the Radiative and Chemically Reactive Flow of Powell–Eyring Fluid" Energies 14, no. 21: 6882. https://doi.org/10.3390/en14216882

APA Style

Firdous, H., Saeed, S. T., Ahmad, H., & Askar, S. (2021). Using Non-Fourier’s Heat Flux and Non-Fick’s Mass Flux Theory in the Radiative and Chemically Reactive Flow of Powell–Eyring Fluid. Energies, 14(21), 6882. https://doi.org/10.3390/en14216882

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop