Dual Thermal Analysis of Fractional Convective Flow Through Aluminum Oxide and Titanium Dioxide Nanoparticles

The applications of free convection flow are perceived in many fields of engineering and science such as heat exchangers, solar energy, drying processes, electric components of communication lines, and thermal storage systems. Several of its applications are also found in conservation, ventilation systems, dehydration, concentration, etc. These forms of flow are frequently inspected with a vertical plate in various manufacturing developments i.e., petroleum industry, geothermal phenomena, thermal insulation, etc. The free convection MHD flow with permeable plate by the AB derivatives along with integral transform is studied in [1]. In [2], free convection with NF in the presence of a magnetic field is deliberated by utilizing a fractional approach. A mathematical study with NFs with radiation impact is studied in [3]. A natural convection fluid flow with a long and vertical cylinder is considered in [4]. They studied the impacts of energy as well as mass transfer with time. Different investigators applied various techniques to study convection fluid flow with diverse structures [5,6, 7, 8, 9, 10, 11].

Firstly, Choi presented NFs comprising nanoparticles in 1995. There are numerous uses of NFs in many fields for example fluid dynamics, in many engineering branches, and biomedical fields. NFs are considered the most excellent alternate ways to usual fluids [ 12, 13, 14]. The small particles in the base fluid which progress the competency of the characteristics of NFs for minimizing the system are known as nanoparticles and the key purpose of NFs is to realize the extreme possible heat conduction at a short nanoparticle concentration. The main conclusions have been distinguished because of the chemical configuration of the nanoparticles when are jumped in the base liquid i.e., reduced possibilities of erosion, thermal transmission, and solidity of the combination. These features play an exceptional part in increasing thermal transmission and energy proficiency in several fields i.e. biomedical instruments, microelectronics, and power generation [15]. Mud nanoparticles consume many practical uses in the penetrating of gasses and oil liquids because of their thermal conduction expressively. The growth of nanoparticles increases the thermal conduction and viscosity of NFs that oppose the intensifying temperature. Currently, because of their extensive uses in many branches of science and technology, NFs are now a fascinating and dominant field for research in fluid mechanics. The impacts of volume fraction and natural convection flow with NFs along with fractional calculus are studied in [16]. A viscous flow with NF (Titania-sodium cellulose) was discussed in [17]. To inspect and enhance the performance of solar accumulators exploited the NFs were studied by Farhana et al. [18]. Jamshed et al. [ 19, 20, 21, 22, 23, 24] applied different techniques to study Eyring NF, Casson NF, second-grade NF, Williamson hybrid NF and tetra hybrid binary NF in different channels. Many researchers applied various methods to study NF flow models with diverse constructions [ 25, 26, 27, 28, 29, 30, 31].

In 1695 [32], firstly, the concept of fractional derivative was given by Leibnitz and L’Hospital which is a proficient tool related to memory facts. Memory function narrates to the kernel of the time-fractional derivative that has not simulated a physical development. Fractional calculus deals with non-local differentiation and integration [33]. The numerical solution of a fractional Oldroyd B-fluid is achieved by the modified Bessel equation as well as the Laplace method [34]. They showed that shear stress is improved as dynamic viscosity is increased. Fractional derivatives are most suitable to the problems of physical nature i.e., earth quick vibrations, polymers, viscoelasticity, heat transfer problems, fluid flows, etc. Over time, various algorithms and definitions were determined by different mathematicians. To find the solution to various mathematical models, the researchers used different fractional derivatives i.e., Riemann-Liouville, Caputo, CF, and AB derivatives. Fractional models can define more proficiently the consequences of the real nature of world problems such as electromagnetic theory, diffusive transfer, electrical networks, fluid flows, rheology, and viscoelastic materials. Then due to a few complications and limitations, Caputo and Fabrizio proposed the latest non-integer order model named CF fractional derivative along with an exponential and non-singular kernel [ 35, 36, 37, 38, 39].

According to the author’s knowledge, there is no investigation on the study of mixed convection flow with Al₂O₃ and TiO₂ nanoparticles with water and blood-based NF in several situations which is a significant theoretical and practical study for the solution of important problems based on the fractional derivative. By getting motivation from these facts, our main purpose is to study a mixed convection flow with Al₂O₃ and TiO₂ nanoparticles with water and blood-based NF along with new definitions of fractional derivatives i.e., AB and CF fractional operators. A semi-analytical approach for AB and CF-based fractional models is applied by the Laplace transform technique along with Stehfest and Tzou’s numerical schemes. To improve the novelty of the recent work some particular cases of velocity profile are also deliberated whose physical importance is prominent in the literature. The graphical illustration for the under-discussed mathematical problem by changing diverse flow parameters is underlined.

2.

CHAPTER TITLE

We assume that a mixed convection fluid flow with Al₂O₃ and TiO₂ nanoparticles are flowing over an inclined plate with an inclination angle δ with the x-axis. Initially, when t = 0, the plate, as well as the fluid, is at rest and ambient medium temperature T_∞. When t = 0⁺, the plate moves by a constant value of velocity $\frac{g (t)}{μ}$ {{g\left( t \right)} \over \mu } where g(0) = 0, and temperature increases from T_∞ to T_w. By this motion of the plate, the fluid begins to move over the plate. Along with all these conditions, it is also supposed the slip impacts the boundaries of the plate. A magnetic field with an angle, θ is also utilized upon the plate as revealed in Fig 1.

We suppose the characteristics of the physical nanoparticles are from Table 1. The fluid velocity as well as the temperature depends on ξ and t. With Boussinesq’s approximation and in the absence of pressure gradient [7, 45], the governing equations are:

Momentum Equation: (1) $\begin{array}{l} ρ_{n f} \frac{\partial w (ξ, t)}{\partial t} = μ_{n f} (1 + α_{1} \frac{\partial}{\partial t}) \frac{\partial^{2} w (ξ, t)}{\partial ξ^{2}} + \\ g {(ρ β_{T})}_{n f} [T (ξ, t) - T_{\infty}] C o s δ - σ_{n f} B_{o}^{2} S i n θ w (ξ, t); ξ, t > 0. \end{array}$ \eqalign{ & {\rho _{nf}}{{\partial w\left( {\xi ,t} \right)} \over {\partial t}} = {\mu _{nf}}\left( {1 + {\alpha _1}{\partial \over {\partial t}}} \right){{{\partial ^2}w\left( {\xi ,t} \right)} \over {\partial {\xi ^2}}} + \cr & g{(\rho {\beta _T})_{nf}}\left[ {T\left( {\xi ,\;t} \right) - {T_\infty }} \right]Cos\delta - {\sigma _{nf}}B_o^2Sin\theta \,w\left( {\xi ,\;t} \right);\xi ,t > 0. \cr}

Tab. 1.

Thermophysical characteristics of base fluids (water and blood) and nanoparticles [ 6,38].

Material	H₂O	Blood	Al₂O₃	TiO₂
ρ(kgm⁻³)	997.1	1053	1600	4250
C_p(kg⁻¹k⁻¹)	0.4179	3594	796	686.2
K(Wm⁻¹k⁻¹)	0.613	0.492	3000	8.9528
B_T×10⁻⁵(k⁻¹⁾	21	0.18	44	0.90

The thermal balance Equation: (2) ${(ρ C_{p})}_{n f} \frac{\partial T (ξ, t)}{\partial t} = - \frac{\partial q_{1}}{\partial ξ}; ξ, t > 0.$ {\left( {\rho {C_p}} \right)_{nf}}{{\partial T\left( {\xi ,t} \right)} \over {\partial t}} = - {{\partial {q_1}} \over {\partial \xi }};\,\,\,\, \xi ,t > 0.

Fourier law [9]: (3) $q_{1} (ξ, t) = - k_{n f} \frac{\partial T (ξ, t)}{\partial ξ}$ {q_1}\left( {\xi ,\;t} \right) = - {k_{nf}}{{\partial T\left( {\xi ,t} \right)} \over {\partial \xi }} with the appropriate initial and boundary conditions (4) $w (ξ, 0) = 0, T (ξ, 0) = T_{\infty}; ξ > 0,$ w\left( {\xi ,\;0} \right) = 0,\,\,\,\, T\left( {\xi ,\;0} \right) = {T_\infty };\,\xi > 0, (5) ${w (0, t) - b \frac{\partial w (ξ, t)}{\partial ξ}|}_{ξ = 0} = \frac{g (t)}{μ}, T (0, t) = T_{w}; t > 0$ {\left. {w\left( {0,\;t} \right) - b{{\partial w\left( {\xi ,t} \right)} \over {\partial \xi }}} \right|_{\xi = 0}} = {{g\left( t \right)} \over \mu },\,\,\,\, T\left( {0,\;t} \right) = {T_w};t > 0 (6) $w (ξ, t) \to 0, T (ξ, t) \to T_{\infty}; ξ \to \infty, t > 0$ w\left( {\xi ,\;t} \right) \to 0,\,\,\,\, T\left( {\xi ,\;t} \right) \to {T_\infty };\,\,\,\, \xi \to \infty ,\,t > 0

The appropriate non-dimensional parameters are taken as (7) $\begin{array}{l} ξ^{*} = \frac{ξ v_{o}}{υ_{f}}, w^{*} = \frac{w}{U_{o}}, t^{*} = \frac{v_{o}^{2} t}{υ_{f}}, ϑ^{*} = \frac{T - T_{\infty}}{T_{w} - T_{\infty}}, \\ b^{*} = \frac{h}{k} b, q^{*} = \frac{q}{q_{o}}, q_{o} = \frac{k_{n f} (T_{w} - T_{\infty}) v_{o}}{υ_{f}}, g^{*} (t^{*}) = \\ \frac{1}{μ} \sqrt{\frac{t_{o}}{υ}} f (t_{o} t^{*}) . \end{array}$ \eqalign{ & {\xi ^*} = {{\xi {v_o}} \over {{\upsilon _f}}}, & {w^*} = {w \over {{U_o}}}, & {t^*} = {{v_o^2t} \over {{\upsilon _f}}}, & {\vartheta ^*} = {{T - {T_\infty }} \over {{T_w} - {T_\infty }}}, \cr & {b^*} = {h \over k}b, & {q^*} = {q \over {{q_o}}}, & {q_o} = {{{k_{nf}}\left( {{T_w} - {T_\infty }} \right){v_o}} \over {{\upsilon _f}}}, & {g^*}\left( {{t^*}} \right) = \cr & {1 \over \mu }\sqrt {{{{t_o}} \over \upsilon }} f({t_o}{t^*}). \cr}

By using the above non-dimensional parameters in Eq. (7), the governing Eqs. (1)–(3) and equivalent conditions (4)–(6) take the form as (8) $\begin{array}{l} \frac{\partial w (ξ, t)}{\partial t} = \frac{1}{Λ_{o} Λ_{1}} (1 + β_{1} \frac{\partial}{\partial t}) \frac{\partial^{2} w (ξ, t)}{\partial ξ^{2}} + \frac{Λ_{2}}{Λ_{o}} G r ϑ (ξ, t) C o s δ - \\ \frac{1}{Λ_{o}} M S i n θ w (ξ, t), \end{array}$ \eqalign{ & {{\partial w\left( {\xi ,t} \right)} \over {\partial t}} = {1 \over {{\Lambda _o}{\Lambda _1}}}\left( {1 + {{\rm{\beta }}_1}{\partial \over {\partial t}}} \right){{{\partial ^2}w\left( {\xi ,t} \right)} \over {\partial {\xi ^2}}} + {{{\Lambda _2}} \over {{\Lambda _o}}}Gr\,\vartheta \left( {\xi ,\;t} \right)Cos\delta - \cr & {1 \over {{\Lambda _o}}}M\,Sin\theta \,w\left( {\xi ,\;t} \right), \cr} (9) $Λ_{3} P r \frac{\partial ϑ (ξ, t)}{\partial t} = - \frac{\partial q (ξ, t)}{\partial ξ}; ξ, t > 0,$ {\Lambda _3}Pr{{\partial \vartheta \left( {\xi ,t} \right)} \over {\partial t}} = - {{\partial q\left( {\xi ,t} \right)} \over {\partial \xi }};\,\xi ,t > 0, (10) $q (ξ, t) = - Λ_{4} \frac{\partial ϑ (ξ, t)}{\partial ξ},$ q\left( {\xi ,\;t} \right) = - {\Lambda _4}{{\partial \vartheta \left( {\xi ,t} \right)} \over {\partial \xi }}, along with corresponding conditions (11) $w (ξ, 0) = 0, ϑ (ξ, 0) = 0; ξ > 0$ w\left( {\xi ,\;0} \right) = 0, \,\,\,\, \vartheta \left( {\xi ,\;0} \right) = 0;\,\xi > 0 (12) ${w (0, t) - b \frac{\partial w (ξ, t)}{\partial ξ}|}_{ξ = 0} = g (t), ϑ (0, t) = 1; t > 0,$ {\left. {w\left( {0,\;t} \right) - b{{\partial w\left( {\xi ,t} \right)} \over {\partial \xi }}} \right|_{\xi = 0}} = g\left( t \right), \,\,\,\, \,\,\,\, \vartheta \left( {0,\;t} \right) = 1;\,\,t > 0, (13) $w (ξ, t) \to 0, ϑ (ξ, t) \to 0; ξ \to \infty, t > 0$ w\left( {\xi ,\;t} \right) \to 0,\,\,\,\,\vartheta \left( {\xi ,\;t} \right) \to 0;\,\xi \to \infty ,\,t > 0 where: $\begin{array}{l} ρ_{n f} = (1 - φ) ρ_{f} + φ ρ_{S}, μ_{n f} = \frac{μ_{f}}{{(1 - φ)}^{2.5}}, \\ {(ρ β)}_{n f} = (1 - φ) {(ρ β)}_{f} + φ {(ρ β)}_{S}, {(ρ C_{p})}_{n f} = \\ (1 - φ) {(ρ C_{p})}_{f} + φ {(ρ C_{p})}_{s}, \\ \frac{k_{n f}}{k_{f}} = \frac{k_{s} + 2 k_{f} - 2 φ (k_{f} - k_{s})}{k_{s} + 2 k_{f} + 2 φ (k_{f} - k_{s})}, \frac{σ_{n f}}{σ_{f}} = 1 + {3 (\frac{σ_{s}}{σ_{f}} - \\ 1) φ} {(\frac{σ_{s}}{σ_{f}} + 2) - (\frac{σ_{s}}{σ_{f}} - 1) φ}^{- 1}, \\ Λ_{o} = (1 - φ) + φ \frac{ρ s}{ρ f}, Λ_{1} = \frac{1}{{(1 - φ)}^{2.5}}, Λ_{2} = \\ (1 - φ) + φ \frac{{(ρ β_{T})}_{S}}{{(ρ β_{T})}_{f}}, \\ Λ_{3} = (1 - φ) + φ \frac{{(ρ C_{p})}_{s}}{{(ρ C_{p})}_{f}}, Λ_{4} = \frac{k_{n f}}{k_{f}}, P r = \frac{{(μ C_{p})}_{f}}{k_{f}}, \\ G r = \frac{g {(β v)}_{f} (T_{w} - T_{\infty})}{U_{o}^{3}}, M = {(\frac{υ_{f}}{v_{o}})}^{2} \frac{σ_{n f} B_{o}^{2}}{ρ_{f} υ_{f}}, β_{1} = α_{1} υ_{f} {(\frac{υ_{f}}{v_{o}})}^{2} . \end{array}$ \eqalign{ & {\rho _{nf}} = \left( {1 - \varphi } \right){\rho _f} + \varphi {\rho _S},{\mu _{nf}} = {{{\mu _f}} \over {{{(1 - \varphi )}^{2.5}}}}, \cr &{(\rho \beta )_{nf}} = \left( {1 - \varphi } \right){(\rho \beta )_f} + \varphi {(\rho \beta )_S}, \,\,\,\, {(\rho {C_p})_{nf}} = \cr & (1 - \varphi ){(\rho {C_p})_f} + \varphi {(\rho {C_p})_s}, \cr & {{{k_{nf}}} \over {{k_f}}} = {{{k_s} + 2{k_f} - 2\varphi \left( {{k_f} - {k_s}} \right)} \over {{k_s} + 2{k_f} + 2\varphi \left( {{k_f} - {k_s}} \right)}},\,{{{\sigma _{nf}}} \over {{\sigma _f}}} = 1 + \left\{ {3\left( {{{{\sigma _s}} \over {{\sigma _f}}} - } \right.} \right. \cr & \left. {\left. 1 \right)\varphi } \right\}{\left\{ {\left( {{{{\sigma _s}} \over {{\sigma _f}}} + 2} \right) - \left( {{{{\sigma _s}} \over {{\sigma _f}}} - 1} \right)\varphi } \right\}^{ - 1}}, \cr & {\Lambda _o} = \left( {1 - \varphi } \right) + \varphi {{\rho s} \over {\rho f}}, \,\,\,\, {\Lambda _1} = {1 \over {{{(1 - \varphi )}^{2.5}}}},{\Lambda _2} = \cr & \left( {1 - \varphi } \right) + \varphi {{{{(\rho {\beta _T})}_S}} \over {{{(\rho {\beta _T})}_f}}}, \cr & {\Lambda _3} = \left( {1 - \varphi } \right) + \varphi {{{{(\rho {C_p})}_s}} \over {{{(\rho {C_p})}_f}}},\, \,\,\,\, {\Lambda _4} = {{{k_{nf}}} \over {{k_f}}}, \,\,\,\, Pr = {{{{(\mu {C_p})}_f}} \over {{k_f}}}, \cr & Gr = {{g{{(\beta v)}_f}\left( {{T_w} - {T_\infty }} \right)} \over {U_o^3}},M = {\left( {{{{\upsilon _f}} \over {{v_o}}}} \right)^2}{{{\sigma _{nf}}B_o^2} \over {{\rho _f}{\upsilon _f}}}, \,\,{\beta _1} = {\alpha _1}{\upsilon _f}{\left( {{{{\upsilon _f}} \over {{v_o}}}} \right)^2}. \cr}

2.1.

Formulation of governing equations by using non-singular kernels

To formulate the fractional model recent proposed definitions of fractional derivatives i.e., AB and CF derivatives. The AB derivative of order 0< β < 1 is defined as [41] (14) $^{A B} 𝔇_{t}^{β} h (t) = \frac{1}{Γ (1 - β)} \int_{0}^{t} E_{β} (\frac{β {(t - ε)}^{β}}{(t - ε)}) h^{'} (t) d ε; 0 < β < 1$ ^{AB}{\mathfrak D}_t^\beta {\rm h}(t) = {1 \over {\Gamma \left( {1 - \beta } \right)}}\mathop \smallint \nolimits_0^t {E_\beta }\left( {{{\beta {{(t - \varepsilon )}^\beta }} \over {(t - \varepsilon )}}} \right)h'(t)d\varepsilon ;0< \beta< 1 and E_β(z) is a Mittage-Leffler function defined by $E_{β} (z) = \sum_{r = 0}^{\infty} \frac{z^{β}}{Γ (γ β + 1)}; 0 < β < 1, z \in ℂ .$ E_\beta \left( z \right) = \mathop \sum \nolimits_{r = 0}^\infty \frac{{z^\beta }} {{\Gamma \left( {\gamma \beta + 1} \right)}};0 < \beta < 1, \,\,\,\,\,\, z \in {\Bbb C}.

The Laplace transform for the AB derivative is [42] (15) $L \{^{A B} 𝔇_{t}^{β} g (ξ, t)\} = \frac{q^{β} L [g (ξ, t)] - q^{β - 1} g (ξ, 0)}{(1 - β) q^{β} + β}$ {\cal L}\left\{ {^{AB}{\mathfrak D}_t^\beta g(\xi ,\;t)} \right\} = {{{q^\beta }{\cal L}[g(\xi ,t)] - {q^{\beta - 1}}g\left( {\xi ,0} \right)} \over {\left( {1 - \beta } \right){q^\beta } + \beta }} with ${\underset{β \to 1}{Lim}}^{A B} 𝔇_{t}^{β} g (ξ, t) = \frac{\partial g (ξ, t)}{\partial t} .$ {\mathop {{\rm{Lim}}}\limits_{\beta \to 1} ^{AB}}{\mathfrak D}_t^\beta g\left( {\xi ,\;t} \right) = {{\partial g\left( {\xi ,t} \right)} \over {\partial t}}.

The CF derivative of order 0< α < 1 is defined as [37,43] (16) $^{C F} 𝔇_{t}^{α} h (t) = \frac{1}{Γ (1 - α)} \int_{0}^{t} E x p (\frac{β {(t - ε)}^{β}}{(t - ε)}) h^{'} (t) d ε, 0 < α < 1,$ ^{CF}{\mathfrak D}_t^\alpha h(t) = {1 \over {\Gamma \left( {1 - \alpha } \right)}}\mathop \smallint \nolimits_0^t Exp\left( {{{\beta {{(t - {\rm{\varepsilon }})}^\beta }} \over {(t - {\rm{\varepsilon }})}}} \right)h'(t)d{\rm{\varepsilon }},\,0< \alpha< 1,

The Laplace transform for the CF derivative is [27,38] (17) $L {^{C F} 𝔇_{t}^{α} g (ξ, t)} = \frac{q L [g (ξ, t)] - g (ξ, 0)}{(1 - α) q + α}$ {\cal L}{\{ ^{CF}}{\mathfrak D}_t^\alpha g\left( {\xi ,\;t} \right)\} = {{q{\cal L}[g(\xi ,t)] - g\left( {\xi ,0} \right)} \over {\left( {1 - \alpha } \right)q + \alpha }} with $\underset{α \to 1}{Lim}^{C F} 𝔇_{t}^{α} g (ξ, t) = \frac{\partial g (ξ, t)}{\partial t} .$ \mathop {{\rm{Lim}}}\limits_{\alpha \to 1} {\,^{CF}}{\mathfrak D}_t^\alpha g\left( {\xi ,\;t} \right) = {{\partial g\left( {\xi ,t} \right)} \over {\partial t}}.

It is important to note that AB and CF fractional operators can also be extended significantly by letting β = 1 in Eq. (14) and α = 1 in Eq. (16) respectively.

The limitation of fractional parameters in derivatives, such as the AB and CF derivatives, to the interval (0,1) derives from their interpretation and physical significance of these values. This range permits for a smooth transition across integer-order derivatives, recording abnormalities as well as long-memory effects in processes, making it ideal for modelling phenomena using sub-diffusive behaviour as well as memory-dependent dynamics in fields such as time series analysis, signal processing, and anomalous diffusion.

3.

MODEL OF NANOFLUID WITH AB DERIVATIVE

The model to the problem with AB derivative can be expressed by substituting the ordinary derivative with AB derivative operator in Eqs. (8)–(10), we get (18) $\begin{array}{l} ^{A B} 𝔇_{t}^{β} w (ξ, t) = \frac{1}{Λ_{o} Λ_{1}} (1 + β_{1}^{A B} 𝔇_{t}^{β}) \frac{\partial^{2} w (ξ, t)}{\partial ξ^{2}} + \\ \frac{Λ_{2}}{Λ_{o}} G r ϑ (ξ, t) C o s δ - \frac{1}{Λ_{o}} M S i n θ w (ξ, t), \end{array}$ \eqalign{ & ^{AB}{\mathfrak D}_t^\beta w\left( {\xi ,\;t} \right) = {1 \over {{\Lambda _o}{\Lambda _1}}}\left( {1 + {{\rm{\beta }}_1}^{AB}{\mathfrak D}_t^\beta } \right){{{\partial ^2}w\left( {\xi ,t} \right)} \over {\partial {\xi ^2}}} + \cr & {{{\Lambda _2}} \over {{\Lambda _o}}}Gr\,\vartheta \left( {\xi ,t} \right)Cos\delta - {1 \over {{\Lambda _o}}}M\,Sin\theta \,w\left( {\xi ,\;t} \right), \cr} (19) $Λ_{3} P r^{A B} 𝔇_{t}^{β} ϑ (ξ, t) = - \frac{\partial q (ξ, t)}{\partial ξ}; ξ, t > 0,$ {\Lambda _3}P{r^{AB}}{\mathfrak D}_t^\beta \,\vartheta \left( {\xi ,\;t} \right) = - {{\partial q\left( {\xi ,t} \right)} \over {\partial \xi }};\,\xi ,\,t > 0, (20) $q (ξ, t) = - Λ_{4} \frac{\partial ϑ (ξ, t)}{\partial ξ} .$ q\left( {\xi ,\;t} \right) = - {\Lambda _4}{{\partial \vartheta (\xi ,t)} \over {\partial \xi }}.

3.1.

Temperature with ab derivative

By employing the Laplace transform on Eqs. (19) and (20), we get (21) $\frac{\partial^{2} \bar{ϑ} (ξ, q)}{\partial ξ^{2}} - \frac{Λ_{3} P r}{Λ_{4}} (\frac{q^{β}}{(1 - β) q^{β} + β}) \bar{ϑ} (ξ, q) = 0,$ {{{\partial ^2}\bar \vartheta \left( {\xi ,q} \right)} \over {\partial {\xi ^2}}} - {{{\Lambda _3}Pr} \over {{\Lambda _4}}}\left( {{{{q^\beta }} \over {(1 - \beta ){q^\beta } + \beta }}} \right)\bar \vartheta (\xi ,\;q) = 0, where ̄ϑ (ξ, q) is the Laplace transform for ϑ (ξ, t), and the transformed conditions after Laplace transform are as follows (22) $\begin{array}{l} \bar{ϑ} (ξ, q) = \frac{1}{q} \\ and \bar{ϑ} (ξ, q) \to 0 as ξ \to \infty . \end{array}$ \eqalign{ & \bar \vartheta (\xi ,\;q) = {1 \over q} \cr & {\rm{and}} & \bar \vartheta (\xi ,\;q) \to 0 & {\rm as} & \xi \to \infty . \cr}

With the above conditions of Eq. (22), we get the temperature as (23) $\bar{ϑ} (ξ, q) = \frac{1}{q} e^{- ξ} \sqrt{\frac{Λ_{3} \Pr}{Λ_{4}} (\frac{q β}{(1 - β) q^{β} + β}) .}$ \bar \vartheta (\xi ,\;q) = {1 \over q}{e^{ - \xi }}\,\sqrt {{{{\Lambda _3}\Pr } \over {{\Lambda _4}}}\left( {{{q\beta } \over {(1 - \beta ){q^\beta } + \beta }}} \right).} Eq (23) can be written as (24) $\bar{ϑ} (ξ, q) = \frac{1}{q} e^{- ξ} \sqrt{\frac{c_{1} q^{γ}}{q^{γ} + c_{2}}}$ \bar \vartheta \left( {\xi ,\;q} \right) = {1 \over q}{e^{ - \xi }}\sqrt {{{{c_1}{q^\gamma }} \over {{q^\gamma } + {c_2}}}} where: $c_{1} = \frac{Λ_{3} P r γ}{Λ_{4}}, c_{2} = β γ, γ = \frac{1}{1 - β} .$ {c_1} = {{{\Lambda _3}Pr\gamma } \over {{\Lambda _4}}},\,\,\,\,{c_2} = \beta \gamma ,\,\,\,\,\gamma = {1 \over {1 - \beta }}.

Eq (24) can also be written in summation form as (25) $\bar{ϑ} (ξ, q) = \frac{1}{q} + \sum_{\infty}^{a_{1} = 1} \sum_{\infty}^{a_{2} = 0} \frac{{(- ξ \sqrt{c_{1}})}^{a_{1}}}{a_{1}!} \frac{{(- c_{2})}^{a_{2}}}{q^{1}^{a_{2}}^{β}} \frac{Γ (\frac{a_{1}}{2} + a_{2})}{Γ (\frac{a_{1}}{2}) Γ (a_{2} + 1)}$ \bar \vartheta \left( {\xi ,\;q} \right) = {1 \over q} + \mathop \sum \nolimits_{{a_1} = 1}^\infty \mathop \sum \nolimits_{{a_2} = 0}^\infty {{( - \xi \sqrt {{c_1}} )^{a_1}} \over {{a_1}!}}{{( - {c_2})^{a_2}} \over {{q^1}{{^ + }^{{a_2}}}{\,^\beta }}}{{\Gamma \left( {{{{a_1}} \over 2} + {a_2}} \right)} \over {\Gamma \left( {{{{a_1}} \over 2}} \right)\Gamma ({a_2} + 1)}}

By taking the Laplace inverse of Eq. (25), we have (26) $ϑ (ξ, t) = 1 + \sum_{a_{1} = 1}^{\infty} \sum_{a_{2} = 0}^{\infty} \frac{{(- ξ \sqrt{c_{1}})}^{a_{1}}}{a_{1}!} \frac{Γ (\frac{a_{1}}{2} + a_{2})}{Γ (\frac{a_{1}}{2}) Γ (a_{2} + 1)} \frac{{(- c_{2})}^{α_{2} t α_{2} β}}{Γ (1 + a_{2} β)}$ \vartheta \left( {\xi ,t} \right) = 1 + \mathop \sum \nolimits_{a_1 = 1}^\infty \mathop \sum \nolimits_{a_2 = 0}^\infty \frac{{( - \xi \sqrt {c_1 } )^{a_1 } }}{{a_1 !}}\frac{{\Gamma \left( {\frac{{a_1 }}{2} + a_2 } \right)}}{{\Gamma \left( {\frac{{a_1 }}{2}} \right)\Gamma (a_2 + 1)}}\frac{{( - c_2 )^{\alpha _2 t\alpha _2 \beta } }}{{\Gamma \left( {1 + a_2 \beta } \right)}}

When β → 1, Eq. (26) becomes (27) $ϑ (ξ, t) = \frac{ξ (1 - e r f (\frac{|ξ| \sqrt{Λ_{3} P r}}{2 \sqrt{Λ_{4} t}}))}{|ξ|}; ξ, \sqrt{\frac{Λ_{3} P r}{Λ_{4}}} > 0.$ \vartheta \left( {\xi ,\;t} \right) = {{\xi \left( {1 - erf\left( {{{\left| \xi \right|\sqrt {{\Lambda _3}Pr} } \over {2\sqrt {{\Lambda _4}t} }}} \right)} \right)} \over {\left| \xi \right|}}; \,\,\,\, \xi ,\sqrt {{{{\Lambda _3}Pr} \over {{\Lambda _4}}}} > 0.

3.2.

Velocity with ab derivative

By using the Laplace transform on Eq. (18), we get (28) $\begin{array}{l} (\frac{q^{β}}{(1 - β) q^{β} + β}) \bar{w} (ξ, q) = \frac{1}{Λ_{o} Λ_{1}} (1 + β_{1} \frac{q^{β}}{(1 - β) q^{β} + β}) \frac{\partial^{2} \bar{w} (ξ, q)}{\partial ξ^{2}} + \\ \frac{Λ_{2}}{Λ_{o}} G r C o s δ \bar{ϑ} (ξ, q) - \frac{1}{Λ_{o}} M S i n θ \bar{w} (ξ, q) \end{array}$ \eqalign{ & \left( {{{{q^\beta }} \over {(1 - \beta ){q^\beta } + \beta }}} \right)\bar w\left( {\xi ,\;q} \right) = {1 \over {{\Lambda _o}{\Lambda _1}}}\left( {1 + {{\rm{\beta }}_1}{{{q^\beta }} \over {(1 - \beta ){q^\beta } + \beta }}} \right){{{\partial ^2}\bar w(\xi ,q)} \over {\partial {\xi ^2}}} + \cr & {{{\Lambda _2}} \over {{\Lambda _o}}}Gr\,Cos\delta \,\bar \vartheta \left( {\xi ,q} \right) - {1 \over {{\Lambda _o}}}M\,Sin\theta \,\bar w\left( {\xi ,\;q} \right) \cr} with the corresponding conditions ${\bar{w} (0, q) - b \frac{\partial \bar{w} (ξ, q)}{\partial ξ}|}_{ξ = 0} = G (q)$ {\left. {\bar w\left( {0,\;q} \right) - b{{\partial \bar w\left( {\xi ,q} \right)} \over {\partial \xi }}} \right|_{\xi = 0}} = G\left( q \right) and (29) $\bar{w} (ξ, q) \to 0 a s ξ \to \infty .$ \bar w(\xi ,\;q) \to 0\,\,\,\,as\,\,\,\,\xi \to \infty .

Eq (28) is solved by utilizing Eq. (29) and we get (30) $\begin{array}{l} \bar{w} (ξ, q) = \frac{1}{1 + b \sqrt{\frac{Λ_{O} Λ_{1} γ q^{β} + (q^{β} + β γ) Λ_{1} M S i n θ}{q^{β} + β γ + β_{1} q^{β} γ}}} \\ (\frac{Λ_{2} G r C o s δ}{Λ_{o} q} \frac{1}{\frac{Λ_{3} Pr}{Λ_{4}} (\frac{q^{β} γ}{q^{β} + β γ}) - \frac{Λ_{O} Λ_{1} γ q^{β} + (q^{β} + β γ) Λ_{1} M S i n θ}{q^{β} + β γ + β_{1} q^{β} γ}} (1 + \\ b \sqrt{\frac{Λ_{3} P r}{Λ_{4}} (\frac{q^{β} γ}{q^{β} + β γ})}) + G (q)) e^{- ξ \sqrt{\frac{Λ_{O} Λ_{1} γ q^{β} + (q^{β} + β γ) Λ_{1} M S i n θ}{q^{β} + β γ + β_{1} q^{β} γ}}} \\ - \frac{Λ_{2} G r C o s δ}{Λ_{o} q} \frac{1}{\frac{Λ_{3} P r}{Λ_{4}} (\frac{q^{β} γ}{q^{β} + β γ}) \frac{Λ_{O} Λ_{1} γ^{β} q + (q^{β} + β γ) Λ_{1} M S i n θ}{q^{β} + β γ + β_{1} q β γ}} e^{- ξ \sqrt{\frac{Λ_{3} P r}{Λ_{4}} (\frac{q^{β} γ}{q^{β} + β γ})}} . \end{array}$ \eqalign{ & \bar w\left( {\xi ,q} \right) = \frac{1}{{1 + b\sqrt {\frac{{\Lambda _O \Lambda _1 \gamma q^\beta + \left( {q^\beta + \beta \gamma } \right)\Lambda _1 MSin\theta }}{{q^\beta + \beta \gamma + \beta _1 q^\beta \gamma }}} }} \cr & \left( {\frac{{\Lambda _2 Gr\,Cos\delta }}{{\Lambda _o q}}} \right.\frac{1}{{\frac{{\Lambda _3 \Pr }}{{\Lambda _4 }}\left( {\frac{{q^\beta \gamma }}{{q^\beta + \beta \gamma }}} \right) - \frac{{\Lambda _O \Lambda _1 \gamma q^\beta + \left( {q^\beta + \beta \gamma } \right)\Lambda _1 M\,Sin\theta }}{{q^\beta + \beta \gamma + \beta _1 q^\beta \gamma }}}}\left( {1 + } \right. \cr & \left. {\left. {b\sqrt {\frac{{\Lambda _3 Pr}}{{\Lambda _4 }}\left( {\frac{{q^\beta \gamma }}{{q^\beta + \beta \gamma }}} \right)} } \right) + G(q)} \right)e^{ - \xi \sqrt {\frac{{\Lambda _O \Lambda _1 \gamma q^\beta + \left( {q^\beta + \beta \gamma } \right)\Lambda _1 M\,Sin\theta }}{{q^\beta + \beta \gamma + \beta _1 q^\beta \gamma }}} } \cr & - \frac{{\Lambda _2 Gr\,Cos\delta }}{{\Lambda _o q}}\frac{1}{{\frac{{\Lambda _3 Pr}}{{\Lambda _4 }}\left( {\frac{{q^\beta \gamma }}{{q^\beta + \beta \gamma }}} \right)\frac{{\Lambda _O \Lambda _1 \gamma ^\beta q + \left( {q^\beta + \beta \gamma } \right)\Lambda _1 M\,Sin\theta }}{{q^\beta + \beta \gamma + \beta _1 q\beta \gamma }}}}e^{ - \xi \sqrt {\frac{{\Lambda _3 Pr}}{{\Lambda _4 }}\left( {\frac{{q^\beta \gamma }}{{q^\beta + \beta \gamma }}} \right)} } . \cr}

When β → 1, Eq. (30) becomes (31) $\begin{array}{l} \bar{w} (ξ, q) = \\ \frac{1}{1 + b \sqrt{\frac{Λ_{O} Λ_{1} q + Λ_{1} M S i n θ}{1 + β_{1} q}}} (\frac{Λ_{2} G r C o s δ}{Λ_{o} q} \frac{1}{\frac{Λ_{3} P r}{Λ_{4}} q - \frac{Λ_{O} Λ_{1} q + Λ_{1} M S i n θ}{1 + β_{1} q}} (1 + \\ b \sqrt{\frac{Λ_{3} P r}{Λ_{4}} q}) + G (q)) e^{- ξ \sqrt{\frac{Λ_{O} Λ_{1} q + Λ_{1} M S i n θ}{1 + β_{1} q}}} - \\ \frac{Λ_{2} G r C o s δ}{Λ_{o} q} \frac{1}{\frac{Λ_{3} P r}{Λ_{4}} q - \frac{Λ_{O} Λ_{1} q + Λ_{1} M S i n θ}{1 + β_{1} q}} e^{- ξ \sqrt{\frac{Λ_{3} P r}{Λ_{4}} q}} \end{array}$ \eqalign{ & \bar w\left( {\xi ,q} \right) = \cr & {1 \over {1 + b\sqrt {{{{\Lambda _O}{\Lambda _1}q + {\Lambda _1}M\,Sin\theta } \over {1 + {\beta _1}q}}} }}\left( {{{{\Lambda _2}Gr\,Cos\delta } \over {{\Lambda _o}q}}} \right.{1 \over {{{{\Lambda _3}Pr} \over {{\Lambda _4}}}q - {{{\Lambda _O}{\Lambda _1}q + {\Lambda _1}M\,Sin\theta } \over {1 + {\beta _1}q}}}}\left( {1 + } \right. \cr & \left. {\left. {b\sqrt {{{{\Lambda _3}Pr} \over {{\Lambda _4}}}q} } \right) + G(q)} \right){e^{ - \xi \sqrt {{{{\Lambda _O}{\Lambda _1}q + {\Lambda _1}M\,Sin\theta } \over {1 + {\beta _1}q}}} }} - \cr & {{{\Lambda _2}Gr\,Cos\delta } \over {{\Lambda _o}q}}{1 \over {{{{\Lambda _3}Pr} \over {{\Lambda _4}}}q - {{{\Lambda _O}{\Lambda _1}q + {\Lambda _1}M\,Sin\theta } \over {1 + {\beta _1}q}}}}{e^{ - \xi \sqrt {{{{\Lambda _3}Pr} \over {{\Lambda _4}}}q} }} \cr}

The Laplace inverse of these solutions is determined numerically through Stehfest as well as Tzou’s approaches as in Tables 2–3.

Tab. 2.

A comparison of solutions with two diverse approaches

ξ	Temperature by Stehfest	Temperature by Tzou	Velocity by Stehfest	Velocity by Tzou
0.1	0.9471	0.9471	0.7001	0.6999
0.2	0.8971	0.8971	0.7856	0.7854
0.3	0.8496	0.8496	0.8532	0.8530
0.4	0.8046	0.8046	0.9053	0.9051
0.5	0.7619	0.7619	0.9438	0.9436
0.6	0.7215	0.7215	0.9707	0.9705
0.7	0.6831	0.6831	0.9875	0.9872
0.8	0.6468	0.6468	0.9956	0.9954
0.9	0.6123	0.6123	0.9964	0.9961

Tab. 3.

Numerical analysis of Nusselt number as well as skin friction for CF and AB derivatives

α, β	Nu by CF	Nu by AB	C_f by AB	C_f by CF
0.1	0.5352	0.5309	0.1836	0.1824
0.2	0.5276	0.5204	0.1751	0.1553
0.3	0.5151	0.5053	0.1611	0.1335
0.4	0.4972	0.4842	0.1423	0.1127
0.5	0.4730	0.4558	0.1203	0.0934
0.6	0.4411	0.4193	0.0965	0.0777
0.7	0.4000	0.3761	0.0728	0.0677
0.8	0.3502	0.3326	0.0513	0.0637
0.9	0.2965	0.2996	0.0359	0.0654

4.

MODEL OF NANOFLUID WITH CF DERIVATIVE

In the above section, the temperature, as well as velocity, is determined by using the AB-fractional derivative, now the modelling of governing equations by CF-fractional derivative may be expressed by replacing the ordinary derivative with CF-fractional derivative operative $^{C F} D_{t}^{α}$ ^{CF} {\mathfrak D}_t^\alpha , the governing equations for the CF-fractional derivative model are attained as (32) $\begin{array}{l} ^{C F} 𝔇_{t}^{α} w (ξ, t) = \frac{1}{Λ_{o} Λ_{1}} (1 + β_{1}^{C F} 𝔇_{t}^{α}) \frac{\partial^{2} w (ξ, t)}{\partial ξ^{2}} + \\ \frac{Λ_{2}}{Λ_{o}} G r ϑ (ξ, t) C o s δ - \frac{1}{Λ_{o}} M S i n θ w (ξ, t), \end{array}$ \eqalign{ & ^{CF}{\mathfrak D}_t^\alpha w\left( {\xi ,\;t} \right) = {1 \over {{\Lambda _o}{\Lambda _1}}}\left( {1 + {\rm{\beta }}_1 ^{CF}{\mathfrak D}_t^\alpha } \right){{{\partial ^2}w\left( {\xi ,t} \right)} \over {\partial {\xi ^2}}} + \cr & {{{\Lambda _2}} \over {{\Lambda _o}}}Gr\,\vartheta \left( {\xi ,\;t} \right)Cos\delta - {1 \over {{\Lambda _o}}}M\,Sin\theta \,w\left( {\xi ,\;t} \right), \cr} (33) $Λ_{3} P r^{C F} 𝔇_{t}^{α} ϑ (ξ, t) = - \frac{\partial q (ξ, t)}{\partial ξ}; ξ, t > 0,$ {\Lambda _3}P{r^{CF}}{\mathfrak D}_t^\alpha \vartheta \left( {\xi ,\;t} \right) = - {{\partial q\left( {\xi ,t} \right)} \over {\partial \xi }}; \,\,\, \xi ,t > 0, (34) $q (ξ, t) = - Λ_{4} \frac{\partial ϑ (ξ, t)}{\partial ξ} .$ q\left( {\xi ,\;t} \right) = - {\Lambda _4}{{\partial \vartheta \left( {\xi ,t} \right)} \over {\partial \xi }}.

4.1.

Temperature with cf derivative

By taking Laplace transform on Eqs. (33), and (34), we get (35) $\frac{\partial^{2} \bar{ϑ} (ξ, q)}{\partial ξ^{2}} - \frac{Λ_{3} \Pr}{Λ_{4}} (\frac{q}{(1 - α) q + α}) \bar{ϑ} (ξ, q) = 0.$ {{{\partial ^2}\bar \vartheta \left( {\xi ,q} \right)} \over {\partial {\xi ^2}}} - {{{\Lambda _3}Pr } \over {{\Lambda _4}}}\left( {{q \over {\left( {1 - \alpha } \right)q + \alpha }}} \right)\bar \vartheta \left( {\xi ,\;q} \right) = 0.

By using the corresponding conditions of Eq. (22), the solution of Eq. (35) is (36) $\bar{ϑ} (ξ, q) = \frac{1}{q} e^{- ξ} \sqrt{\frac{Λ_{3} P r}{Λ_{4}} (\frac{q}{(1 - α) q + α}) .}$ \bar \vartheta \left( {\xi ,\;q} \right) = {1 \over q}e{\;^{ - \xi }}\sqrt {{{{\Lambda _3}Pr} \over {{\Lambda _4}}}\left( {{q \over {(1 - \alpha )q + \alpha }}} \right).}

Eq (36) can be written as (37) $\bar{ϑ} (ξ, q) = \frac{1}{q} e^{- ξ \sqrt{\frac{d_{1} q}{q + d_{2}}}}$ \bar \vartheta \left( {\xi ,\;q} \right) = {1 \over q}{e^{ - \xi \sqrt {{{{d_1}q} \over {q + {d_2}}}} }}

Where $d_{1} = \frac{Λ_{3} P r γ}{Λ_{4}}, d_{2} = α γ, γ = \frac{1}{1 - α}$ {d_1} = {{{\Lambda _3}Pr\gamma } \over {{\Lambda _4}}},\,\,\,\,{d_2} = \alpha \gamma ,\,\,\,\,\gamma = {1 \over {1 - \alpha }}

Eq (37) may be written in summation form as (38) $\bar{ϑ} (ξ, q) = \frac{1}{q} + \sum_{a_{1} = 1}^{\infty} \sum_{a_{2} = 0}^{\infty} \frac{{(- ξ \sqrt{d_{1}})}^{a_{1}}}{a_{1}!} \frac{{(- d_{2})}^{a_{2}}}{q^{1 + a}^{_{2}}} \frac{Γ (\frac{a_{1}}{2} + a_{2})}{Γ (\frac{a_{1}}{2}) Γ (a_{2} + 1)}$ \bar \vartheta \left( {\xi ,\;q} \right) = {1 \over q} + \mathop \sum \nolimits_{{a_1} = 1}^\infty \mathop \sum \nolimits_{{a_2} = 0}^\infty {{{{( - \xi \sqrt {{d_1}} )}^{{a_1}}}} \over {{a_1}!}}{{{{( - {d_2})}^{{a_2}}}} \over {{q^{1 + a}}^{_2}}}{{\Gamma \left( {{{{a_1}} \over 2} + {a_2}} \right)} \over {\Gamma \left( {{{{a_1}} \over 2}} \right)\Gamma ({a_2} + 1)}}

By utilizing the Laplace inverse of Eq. (38), we have (39) $ϑ (ξ, t) = 1 + \sum_{a_{1} = 1}^{\infty} \sum_{a_{2} = 0}^{\infty} \frac{{(- ξ \sqrt{d_{1}})}^{a_{1}}}{a_{1}!} \frac{{(- d_{2})}^{a_{2}} t^{a_{2}}}{Γ (a_{2} + 1)} \frac{Γ (\frac{a_{1}}{2} + a_{2})}{Γ (\frac{a_{1}}{2}) Γ (a_{2} + 1)}$ \vartheta \left( {\xi ,\;t} \right) = 1 + \mathop \sum \nolimits_{{a_1} = 1}^\infty \mathop \sum \nolimits_{{a_2} = 0}^\infty {{{{( - \xi \sqrt {{d_1}} )}^{{a_1}}}} \over {{a_1}!}}{{{{( - {d_2})}^{{a_2}}}{t^{{a_2}}}} \over {\Gamma \left( {{a_2} + 1} \right)}}{{\Gamma \left( {{{{a_1}} \over 2} + {a_2}} \right)} \over {\Gamma \left( {{{{a_1}} \over 2}} \right)\Gamma \left( {{a_2} + 1} \right)}}

For the special case for ordinary solution replace α → 1 in Eq (39), and then the solution will be (40) $ϑ (ξ, t) = \frac{ξ (1 - e r f (\frac{|ξ| \sqrt{Λ_{3} P r}}{2 \sqrt{Λ_{4} t}}))}{|ξ|}; ξ, \sqrt{\frac{Λ_{3} P r}{Λ_{4}}} > 0.$ \vartheta \left( {\xi ,\;t} \right) = {{\xi \left( {1 - erf\left( {{{\left| \xi \right|\sqrt {{\Lambda _3}Pr} } \over {2\sqrt {{\Lambda _4}t} }}} \right)} \right)} \over {\left| \xi \right|}};\xi ,\sqrt {{{{\Lambda _3}Pr} \over {{\Lambda _4}}}} > 0.

4.2.

Velocity with cf derivative

By using the Laplace transform on Eqs. (32), we get (41) $\begin{array}{l} (\frac{s}{(1 - α) s + α}) \bar{w} (ξ, s) = \frac{1}{Λ_{o} Λ_{1}} (1 + β_{1} (\frac{s}{(1 - α) s + α})) \frac{\partial^{2} \bar{w} (ξ, s)}{\partial ξ^{2}} + \\ \frac{Λ_{2}}{Λ_{o}} G r \bar{ϑ} (ξ, s) C o s δ - \frac{1}{Λ_{o}} M S i n θ \bar{w} (ξ, s) . \end{array}$ \eqalign{ & \left( {{s \over {(1 - \alpha )s + \alpha }}} \right)\bar w\left( {\xi ,s} \right) = {1 \over {{\Lambda _o}{\Lambda _1}}}\left( {1 + {{\rm{\beta }}_1}\left( {{s \over {(1 - \alpha )s + \alpha }}} \right)} \right){{{\partial ^2}\bar w\left( {\xi ,s} \right)} \over {\partial {\xi ^2}}} + \cr & {{{\Lambda _2}} \over {{\Lambda _o}}}Gr\bar \vartheta \left( {\xi ,\;s} \right)Cos\delta - {1 \over {{\Lambda _o}}}M\,Sin\theta \bar w\left( {\xi ,\;s} \right). \cr}

By solving Eq (41) with the corresponding conditions in Eq. (29), we get (42) $\begin{array}{l} \bar{w} (ξ, q) = \frac{1}{1 + b \sqrt{\frac{Λ_{O} Λ_{1} γ q + (q + α γ) Λ_{1} M S i n θ}{q + α γ + β_{1} q γ}}} \\ (\frac{Λ_{2} G r C o s δ}{Λ_{o^{S}}} \frac{1}{\frac{Λ_{3} P r}{Λ_{4}} (\frac{q γ}{q + α γ}) - \frac{Λ_{O} Λ_{1} γ q + (q + α γ) Λ_{1} M S i n θ}{q + α γ + β_{1} q γ}} (1 + \\ b \sqrt{\frac{Λ_{3} P r}{Λ_{4}} (\frac{q γ}{q + α γ})} + G (q)) e^{- ξ \sqrt{\frac{Λ_{O} Λ_{1} γ q + (q + α γ) Λ_{1} M S i n θ}{q + α γ + β_{1} q γ}}} \\ - \frac{Λ_{2} G r C o s δ}{Λ_{o} q} \frac{1}{\frac{Λ_{3} P r}{Λ_{4}} (\frac{q γ}{q + α γ}) - \frac{Λ_{O} Λ_{1} γ q + (q + α γ) Λ_{1} M S i n θ}{q + α γ + β_{1} q γ}} e^{- ξ \sqrt{\frac{Λ_{3} P r}{Λ_{4}} (\frac{q γ}{q + α γ})}} . \end{array}$ \eqalign{ & \bar w\left( {\xi ,\;q} \right) = {1 \over {1 + b\sqrt {{{{\Lambda _O}{\Lambda _1}\gamma q + \left( {q + \alpha \gamma } \right){\Lambda _1}M\,Sin\theta } \over {q + \alpha \gamma + {\beta _1}q\gamma }}} }} \cr & \left( {{{{\Lambda _2}Gr\,Cos\delta } \over {{\Lambda _{{o^S}}}}}{1 \over {{{{\Lambda _3}Pr} \over {{\Lambda _4}}}\left( {{{q\gamma } \over {q + \alpha \gamma }}} \right) - {{{\Lambda _O}{\Lambda _1}\gamma q + \left( {q + \alpha \gamma } \right){\Lambda _1}M\,Sin\theta } \over {q + \alpha \gamma + {\beta _1}q\gamma }}}}} \right.\left( {1 + } \right. \cr & \left. {b\sqrt {{{{\Lambda _3}Pr} \over {{\Lambda _4}}}\left( {{{q\gamma } \over {q + \alpha \gamma }}} \right)} + G\left( q \right)} \right){e^{ - \xi \sqrt {{{{\Lambda _O}{\Lambda _1}\gamma q + \left( {q + \alpha \gamma } \right){\Lambda _1}M\,Sin\theta } \over {q + \alpha \gamma + {\beta _1}q\gamma }}} }} \cr & - {{{\Lambda _2}Gr\,Cos\delta } \over {{\Lambda _o}q}}{1 \over {{{{\Lambda _3}Pr} \over {{\Lambda _4}}}\left( {{{q\gamma } \over {q + \alpha \gamma }}} \right) - {{{\Lambda _O}{\Lambda _1}\gamma q + \left( {q + \alpha \gamma } \right){\Lambda _1}M\,Sin\theta } \over {q + \alpha \gamma + {\beta _1}q\gamma }}}}{e^{ - \xi \sqrt {{{{\Lambda _3}Pr} \over {{\Lambda _4}}}\left( {{{q\gamma } \over {q + \alpha \gamma }}} \right)} }}. \cr}

When α → 1, Eq. (42) becomes (43) $\begin{array}{l} \bar{w} (ξ, q) = \\ \frac{1}{1 + b \sqrt{\frac{Λ_{O} Λ_{1} q + Λ_{1} M S i n θ}{1 + β_{1} q}}} (\frac{Λ_{2} G r C o s δ}{Λ_{o} q} \frac{1}{\frac{Λ_{3} P r}{Λ_{4}} q - \frac{Λ_{O} Λ_{1} q + Λ_{1} M S i n θ}{1 + β_{1} q}} (1 + \\ b \sqrt{\frac{Λ_{3} P r}{Λ_{4}} q}) + G (q)) e^{- ξ \sqrt{\frac{Λ_{O} Λ_{1} q + Λ_{1} M S i n θ}{1 + β_{1} q}}} - \\ \frac{Λ_{2} G r C o s δ}{Λ_{o} q} \frac{1}{\frac{Λ_{3} P r}{Λ_{4}} q - \frac{Λ_{O} Λ_{1} q + Λ_{1} M S i n θ}{1 + β_{1} q}} e^{- ξ \sqrt{\frac{Λ_{3} P r}{Λ_{4}} q}} . \end{array}$ \eqalign{ & \bar w\left( {\xi ,q} \right) = \cr & {1 \over {1 + b\sqrt {{{{\Lambda _O}{\Lambda _1}q + {\Lambda _1}M\,Sin\theta } \over {1 + {\beta _1}q}}} }}\left( {{{{\Lambda _2}Gr\,Cos\delta } \over {{\Lambda _o}q}}{1 \over {{{{\Lambda _3}Pr} \over {{\Lambda _4}}}q - {{{\Lambda _O}{\Lambda _1}q + {\Lambda _1}M\,Sin\theta } \over {1 + {\beta _1}q}}}}\left( {1 + } \right.} \right. \cr & \left. {\left. {b\sqrt {{{{\Lambda _3}Pr} \over {{\Lambda _4}}}q} } \right) + G(q)} \right){e^{ - \xi \sqrt {{{{\Lambda _O}{\Lambda _1}q + {\Lambda _1}M\,Sin\theta } \over {1 + {\beta _1}q}}} }} - \cr & {{{\Lambda _2}Gr\,Cos\delta } \over {{\Lambda _o}q}}{1 \over {{{{\Lambda _3}Pr} \over {{\Lambda _4}}}q - {{{\Lambda _O}{\Lambda _1}q + {\Lambda _1}M\,Sin\theta } \over {1 + {\beta _1}q}}}}{e^{ - \xi \sqrt {{{{\Lambda _3}Pr} \over {{\Lambda _4}}}q} }}. \cr}

To find out the Laplace inverse, various researchers applied different numerical approaches to find the solution of diverse differential fractional models as in [44,45,46,47]. Consequently, here we will also utilize the Stehfest scheme to find the numerical solution of temperature as well as velocity numerically. Grave Stehfest scheme [48] can be expressed as (44) $w (ξ, t) = \frac{l n (2)}{t} \sum_{n = 1}^{M} v_{n} \bar{w} (ξ, n \frac{l n (2)}{t})$ w\left( {\xi ,\;t} \right) = {{ln\left( 2 \right)} \over t}\mathop \sum \nolimits_{n = 1}^M {v_n}\bar w\left( {\xi ,n{{ln\left( 2 \right)} \over t}} \right) where M be a non-negative integer, and (45) $v_{n} = {(- 1)}^{n + \frac{M}{2}} \sum_{p = [\frac{q + 1}{2}]}^{\min (q, \frac{M}{2})} \frac{p^{\frac{M}{2}} (2 p)!}{(\frac{M}{2} - p)! p! (p - 1)! (q - p)! (2 p - q)!}$ {v_n} = {( - 1)^{n + {M \over 2}}}\mathop \sum \nolimits_{p = \left[ {{{q + 1} \over 2}} \right]}^{\min \left( {q,{M \over 2}} \right)} {{{p^{{M \over 2}}}\left( {2p} \right)!} \over {\left( {{M \over 2} - p} \right)!p!\left( {p - 1} \right)!\left( {q - p} \right)!\left( {2p - q} \right)!}}

However, we also applied another estimation for temperature as well as velocity solutions, Tzou’s method for the comparison and validation of our numerical findings with the Stehfest [48] scheme. Tzou’s scheme [49] has the form as (46) $\begin{array}{l} w (ξ, t) = \frac{e^{4.7}}{t} [\frac{1}{2} \bar{w} (r, \frac{4.7}{t}) + \\ R e \{\sum_{j = 1}^{N} {(- 1)}^{k} \bar{w} (r, \frac{4.7 + k π i}{t})\}] \end{array}$ \eqalign{ & w\left( {\xi ,\;t} \right) = {{{e^{4.7}}} \over t}\left[ {{1 \over 2}\bar w\left( {r,{{4.7} \over t}} \right) + } \right. \cr & \left. {Re\left\{ {\mathop \sum \nolimits_{j = 1}^N {{( - 1)}^k}\bar w\left( {r,{{4.7 + k\pi i} \over t}} \right)} \right\}} \right] \cr} where i and Re(. ) are imaginary units and real portions and N > 1 is a natural number.

5.

PARTICULAR CASES

As the solution of the velocity field with AB and CF derivative in Eq. (30) and (42) correspondingly, is in a more general form. Consequently, to demonstrate some more physical perception of the problem, we will deliberate some particular cases for the function g(t) for the velocity whose physical explanation is prominent in the literature.

Case 1: g(t) = t

In this case, we take g(t) = t then the expressions of velocity with AB and CF derivative along with Eqs. (30) and (42) respectively will take the form as (47) $\begin{array}{l} \bar{w} (ξ, q) = \frac{1}{1 + b \sqrt{\frac{Λ_{O} Λ_{1} γ q^{β} + (q^{β} + β γ) Λ_{1} M S i n θ}{q^{β} + β γ + β_{1} q β γ}}}, \\ (\frac{Λ_{2} G r C o s δ}{Λ_{o} q} \frac{1}{\frac{Λ_{3} P r}{Λ_{4}} (\frac{q^{β} γ}{q^{β} + β γ}) - \frac{Λ_{O} Λ_{1} γ q^{β} + (q^{β} + β γ) Λ_{1} M S i n θ}{q^{β} + β γ + β_{1} q^{β} γ}} (1 + \\ b \sqrt{\frac{Λ_{3} P r}{Λ_{4}} (\frac{q^{β} γ}{q^{β} + β γ})}) + \frac{1}{q^{2}}) e^{- ξ \sqrt{\frac{Λ_{O} Λ_{1} γ q^{β} + (q^{β} + β γ) Λ_{1} M S i n θ}{q^{β} + β γ + β_{1} q^{β} γ}}} \\ - \frac{Λ_{2} G r C o s δ}{Λ_{o} q} \frac{1}{\frac{Λ_{3} P r}{Λ_{4}} (\frac{q^{β} γ}{q^{β} + β γ}) - \frac{Λ_{O} Λ_{1} γ q^{β} + (q^{β} + β γ) Λ_{1} M S i n θ}{q^{β} + β γ + β_{1} q^{β} γ}} e^{- ξ \sqrt{\frac{Λ_{3} P r}{Λ_{4}} (\frac{q^{β} γ}{q^{β} + β γ})}} . \end{array}$ \eqalign{ & \bar w\left( {\xi ,q} \right) = {1 \over {1 + b\sqrt {{{{\Lambda _O}{\Lambda _1}\gamma {q^\beta } + \left( {{q^\beta } + \beta \gamma } \right){\Lambda _1}M\,Sin\theta } \over {{q^\beta } + \beta \gamma + {\beta _1}q\beta \gamma }}} }}, \cr & \left( {{{{\Lambda _2}Gr\,Cos\delta } \over {{\Lambda _o}q}}{1 \over {{{{\Lambda _3}Pr} \over {{\Lambda _4}}}\left( {{{{q^\beta }\gamma } \over {{q^\beta } + \beta \gamma }}} \right) - {{{\Lambda _O}{\Lambda _1}\gamma {q^\beta } + \left( {{q^\beta } + \beta \gamma } \right){\Lambda _1}M\,Sin\theta } \over {{q^\beta } + \beta \gamma + {\beta _1}{q^\beta }\gamma }}}}\left( {1 + } \right.} \right. \cr & \left. {\left. {b\sqrt {{{{\Lambda _3}Pr} \over {{\Lambda _4}}}\left( {{{{q^\beta }\gamma } \over {{q^\beta } + \beta \gamma }}} \right)} } \right) + {1 \over {{q^2}}}} \right){e^{ - \xi \sqrt {{{{\Lambda _O}{\Lambda _1}\gamma {q^\beta } + \left( {{q^\beta } + \beta \gamma } \right){\Lambda _1}M\,Sin\theta } \over {{q^\beta } + \beta \gamma + {\beta _1}{q^\beta }\gamma }}} }} \cr & - {{{\Lambda _2}Gr\,Cos\delta } \over {{\Lambda _o}q}}{1 \over {{{{\Lambda _3}Pr} \over {{\Lambda _4}}}\left( {{{{q^\beta }\gamma } \over {{q^\beta } + \beta \gamma }}} \right) - {{{\Lambda _O}{\Lambda _1}\gamma {q^\beta } + \left( {{q^\beta } + \beta \gamma } \right){\Lambda _1}M\,Sin\theta } \over {{q^\beta } + \beta \gamma + {\beta _1}{q^\beta }\gamma }}}}{e^{ - \xi \sqrt {{{{\Lambda _3}Pr} \over {{\Lambda _4}}}\left( {{{{q^\beta }\gamma } \over {{q^\beta } + \beta \gamma }}} \right)} }}. \cr} and (48) $\begin{array}{l} \bar{w} (ξ, q) = \frac{1}{1 + b \sqrt{\frac{Λ_{O} Λ_{1} γ q + (q + α γ) Λ_{1} M S i n θ}{q + α γ + β_{1} q γ}}} \\ (\frac{Λ_{2} G r C o s δ}{Λ_{o} q} \frac{}{\frac{Λ_{3} P r}{Λ_{4}} (\frac{q γ}{q + α γ}) - \frac{Λ_{O} Λ_{1} γ q + (q + α γ) Λ_{1} M S i n θ}{q + α γ + β_{1} q γ}} (1 + \\ b \sqrt{\frac{Λ_{3} P r}{Λ_{4}} (\frac{q γ}{q + α γ})}) + \frac{1}{q^{2}}) e^{- ξ \sqrt{\frac{Λ_{O} Λ_{1} γ q + (q + α γ) Λ_{1} M S i n θ}{q + α γ + β_{1} q γ}}} \\ - \frac{Λ_{2} G r C o s δ}{Λ_{o} s} \frac{1}{\frac{Λ_{3} P r}{Λ_{4}} (\frac{q γ}{q + α γ}) - \frac{Λ_{O} Λ_{1} γ q + (q + α γ) Λ_{1} M S i n θ}{q + α γ + β_{1} q γ}} e^{- ξ \sqrt{\frac{Λ_{3} P r}{Λ_{4}} (\frac{q γ}{q + α γ})}} . \end{array}$ \eqalign{ & \bar w\left( {\xi ,q} \right) = {1 \over {1 + b\sqrt {{{{\Lambda _O}{\Lambda _1}\gamma q + \left( {q + \alpha \gamma } \right){\Lambda _1}M\,Sin\theta } \over {q + \alpha \gamma + {\beta _1}q\gamma }}} }} \cr & \left( {{{{\Lambda _2}Gr\,Cos\delta } \over {{\Lambda _o}q}}{{} \over {{{{\Lambda _3}Pr} \over {{\Lambda _4}}}\left( {{{q\gamma } \over {q + \alpha \gamma }}} \right) - {{{\Lambda _O}{\Lambda _1}\gamma q + \left( {q + \alpha \gamma } \right){\Lambda _1}M\,Sin\theta } \over {q + \alpha \gamma + {\beta _1}q\gamma }}}}\left( {1 + } \right.} \right. \cr & \left. {\left. {b\sqrt {{{{\Lambda _3}Pr} \over {{\Lambda _4}}}\left( {{{q\gamma } \over {q + \alpha \gamma }}} \right)} } \right) + {1 \over {{q^2}}}} \right){e^{ - \xi \sqrt {{{{\Lambda _O}{\Lambda _1}\gamma q + \left( {q + \alpha \gamma } \right){\Lambda _1}M\,Sin\theta } \over {q + \alpha \gamma + {\beta _1}q\gamma }}} }} \cr & - {{{\Lambda _2}Gr\,Cos\delta } \over {{\Lambda _o}s}}{1 \over {{{{\Lambda _3}Pr} \over {{\Lambda _4}}}\left( {{{q\gamma } \over {q + \alpha \gamma }}} \right) - {{{\Lambda _O}{\Lambda _1}\gamma q + \left( {q + \alpha \gamma } \right){\Lambda _1}M\,Sin\theta } \over {q + \alpha \gamma + {\beta _1}q\gamma }}}}{e^{ - \xi \sqrt {{{{\Lambda _3}Pr} \over {{\Lambda _4}}}\left( {{{q\gamma } \over {q + \alpha \gamma }}} \right)} }}. \cr & \cr}

Case 2: g(t) = Sin(ωt)

In this case, we take g(t) = Sin(ωt) where ω denotes the intensity of the shear stress, then the expressions for velocity with AB and CF derivative with Eqs. (30) and (42) correspondingly will take the form as (49) $\begin{array}{l} \bar{w} (ξ, q) = \frac{1}{1 + b \sqrt{\frac{Λ_{O} Λ_{1} γ q^{β} + (q^{β} + β γ) Λ_{1} M S i n θ}{q^{β} + β γ + β_{1} q^{β} γ}}} \\ (\frac{Λ_{2} G r C o s δ}{Λ_{o} q} \frac{1}{\frac{Λ_{3} P r}{Λ_{4}} (\frac{q^{β} γ}{q^{β} + β γ}) - \frac{Λ_{O} Λ_{1} γ q^{β} + (q^{β} + β γ) Λ_{1} M S i n θ}{q^{β} + β γ + β_{1} q^{β} γ}} (1 + \\ b \sqrt{\frac{Λ_{3} P r}{Λ_{4}} (\frac{q^{β} γ}{q^{β} + β γ})}) + \frac{ω}{ω^{2} + q^{2}}) e^{- ξ \sqrt{\frac{Λ_{O} Λ_{1} γ q^{β} + (q^{β} + β γ) Λ_{1} M S i n θ}{q^{β} + β γ + β_{1} q^{β} γ}}} \\ - \frac{Λ_{2} G r C o s δ}{Λ_{o} q} \frac{1}{\frac{Λ_{3} P r}{Λ_{4}} (\frac{q^{β} γ}{q^{β} + β γ}) - \frac{Λ_{O} Λ_{1} γ q^{β} + (q^{β} + β γ) Λ_{1} M S i n θ}{q^{β} + β γ + β_{1} q^{β} γ}} e^{- ξ \sqrt{\frac{Λ_{3} P r}{Λ_{4}} (\frac{q^{β} γ}{q^{β} + β γ})}} \end{array}$ \eqalign{ & \bar w\left( {\xi ,q} \right) = {1 \over {1 + b\sqrt {{{{\Lambda _O}{\Lambda _1}\gamma {q^\beta } + \left( {{q^\beta } + \beta \gamma } \right){\Lambda _1}M\,Sin\theta } \over {{q^\beta } + \beta \gamma + {\beta _1}{q^\beta }\gamma }}} }} \cr & \left( {{{{\Lambda _2}Gr\,Cos\delta } \over {{\Lambda _o}q}}{1 \over {{{{\Lambda _3}Pr} \over {{\Lambda _4}}}\left( {{{{q^\beta }\gamma } \over {{q^\beta } + \beta \gamma }}} \right) - {{{\Lambda _O}{\Lambda _1}\gamma {q^\beta } + \left( {{q^\beta } + \beta \gamma } \right){\Lambda _1}M\,Sin\theta } \over {{q^\beta } + \beta \gamma + {\beta _1}{q^\beta }\gamma }}}}\left( {1 + } \right.} \right. \cr & \left. {\left. {b\sqrt {{{{\Lambda _3}Pr} \over {{\Lambda _4}}}\left( {{{{q^\beta }\gamma } \over {{q^\beta } + \beta \gamma }}} \right)} } \right) + {\omega \over {{\omega ^2} + {q^2}}}} \right){e^{ - \xi \sqrt {{{{\Lambda _O}{\Lambda _1}\gamma {q^\beta } + \left( {{q^\beta } + \beta \gamma } \right){\Lambda _1}M\,Sin\theta } \over {{q^\beta } + \beta \gamma + {\beta _1}{q^\beta }\gamma }}} }} \cr & - {{{\Lambda _2}Gr\,Cos\delta } \over {{\Lambda _o}q}}{1 \over {{{{\Lambda _3}Pr} \over {{\Lambda _4}}}\left( {{{{q^\beta }\gamma } \over {{q^\beta } + \beta \gamma }}} \right) - {{{\Lambda _O}{\Lambda _1}\gamma {q^\beta } + \left( {{q^\beta } + \beta \gamma } \right){\Lambda _1}M\,Sin\theta } \over {{q^\beta } + \beta \gamma + {\beta _1}{q^\beta }\gamma }}}}{e^{ - \xi \sqrt {{{{\Lambda _3}Pr} \over {{\Lambda _4}}}\left( {{{{q^\beta }\gamma } \over {{q^\beta } + \beta \gamma }}} \right)} }} \cr} and (50) $\begin{array}{l} \bar{w} (ξ, q) = \frac{1}{1 + b \sqrt{\frac{Λ_{O} Λ_{1} γ q + (q + α γ) Λ_{1} M S i n θ}{q + α γ + β_{1} q γ}}} \\ (\frac{Λ_{2} G r C o s δ}{Λ_{o} q} \frac{1}{\frac{Λ_{3} P r}{Λ_{4}} (\frac{q γ}{q + α γ}) - \frac{Λ_{O} Λ_{1} γ q + (q + α γ) Λ_{1} M S i n θ}{q + α γ + β_{1} q γ}} (1 + \\ b \sqrt{\frac{Λ_{3} P r}{Λ_{4}} (\frac{q γ}{q + α γ})}) + \frac{ω}{ω^{2} + q^{2}}) e^{- ξ \sqrt{\frac{Λ_{O} Λ_{1} γ q + (q + α γ) Λ_{1} M S i n θ}{q + α γ + β_{1} q γ}}} \\ - \frac{Λ_{2} G C o s δ}{Λ_{o} q} \frac{1}{\frac{Λ_{3} P r}{Λ_{4}} (\frac{q γ}{q + α γ}) - \frac{Λ_{O} Λ_{1} γ q + (q + α γ) Λ_{1} M S i n θ}{q + α γ + β_{1} q γ}} e^{- ξ \sqrt{\frac{Λ_{3} P r}{Λ_{4}} (\frac{q γ}{q + α γ})}} \end{array}$ \eqalign{ & \bar w\left( {\xi ,\;q} \right) = {1 \over {1 + b\sqrt {{{{\Lambda _O}{\Lambda _1}\gamma q + \left( {q + \alpha \gamma } \right){\Lambda _1}M\,Sin\theta } \over {q + \alpha \gamma + {\beta _1}q\gamma }}} }} \cr & \left( {{{{\Lambda _2}Gr\,Cos\delta } \over {{\Lambda _o}q}}{1 \over {{{{\Lambda _3}Pr} \over {{\Lambda _4}}}\left( {{{q\gamma } \over {q + \alpha \gamma }}} \right) - {{{\Lambda _O}{\Lambda _1}\gamma q + \left( {q + \alpha \gamma } \right){\Lambda _1}M\,Sin\theta } \over {q + \alpha \gamma + {\beta _1}q\gamma }}}}\left( {1 + } \right.} \right. \cr & \left. {\left. {b\sqrt {{{{\Lambda _3}Pr} \over {{\Lambda _4}}}\left( {{{q\gamma } \over {q + \alpha \gamma }}} \right)} } \right) + {\omega \over {{\omega ^2} + {q^2}}}} \right){e^{ - \xi \sqrt {{{{\Lambda _O}{\Lambda _1}\gamma q + \left( {q + \alpha \gamma } \right){\Lambda _1}M\,Sin\theta } \over {q + \alpha \gamma + {\beta _1}q\gamma }}} }} \cr & - {{{\Lambda _2}G\,Cos\delta } \over {{\Lambda _o}q}}{1 \over {{{{\Lambda _3}Pr} \over {{\Lambda _4}}}\left( {{{q\gamma } \over {q + \alpha \gamma }}} \right) - {{{\Lambda _O}{\Lambda _1}\gamma q + \left( {q + \alpha \gamma } \right){\Lambda _1}M\,Sin\theta } \over {q + \alpha \gamma + {\beta _1}q\gamma }}}}{e^{ - \xi \sqrt {{{{\Lambda _3}Pr} \over {{\Lambda _4}}}\left( {{{q\gamma } \over {q + \alpha \gamma }}} \right)} }} \cr}

Case 3: g(t) = t Cos(t)

In this case, we take g(t) = t Cost with its Laplace $G (q) = \frac{q^{2} - 1}{{(q^{2} + 1)}^{2}}$ G\left( q \right) = {{{q^2} - 1} \over {{{({q^2} + 1)}^2}}} , then the expressions of velocity through AB and CF derivative along with Eqs. (30) and (42) correspondingly will take the form as (51) $\begin{array}{l} \bar{w} (ξ, q) = \frac{1}{1 + b \sqrt{\frac{Λ_{O} Λ_{1} γ q^{β} + (q^{β} + β γ) Λ_{1} M S i n θ}{q^{β} + β γ + β_{1} q^{β} γ}}} \\ (\frac{Λ_{2} G r C o s δ}{Λ_{o} q} \frac{1}{\frac{Λ_{3} P r}{Λ_{4}} (\frac{q^{β} γ}{q^{β} + β γ}) - \frac{Λ_{O} Λ_{1} γ q^{β} + (q^{β} + β γ) Λ_{1} M S i n θ}{q^{β} + β γ + β_{1} q^{β} γ}} (1 + \\ b \sqrt{\frac{Λ_{3} P r}{Λ_{4}} (\frac{q^{β} γ}{q^{β} + β γ})}) + \frac{q^{2} - 1}{{(1 + q^{2})}^{2}}) e^{- ξ \sqrt{\frac{Λ_{O} Λ_{1} γ q^{β} + (q^{β} + β γ) Λ_{1} M S i n θ}{q^{β} + β γ + β_{1} q^{β} γ}}} \\ - \frac{Λ_{2} G r C o s δ}{Λ_{o} q} \frac{1}{\frac{Λ_{3} P r}{Λ_{4}} (\frac{q^{β} γ}{q^{β} + β γ}) - \frac{Λ_{O} Λ_{1} γ q^{β} + (q^{β} + β γ) Λ_{1} M S i n θ}{q^{β} + β γ + β_{1} q^{β} γ}} e^{- ξ \sqrt{\frac{Λ_{3} P r}{Λ_{4}} (\frac{q^{β} γ}{q^{β} + β γ})}} \end{array}$ \eqalign{ & \bar w\left( {\xi ,q} \right) = {1 \over {1 + b\sqrt {{{\Lambda _O \Lambda _1 \gamma q^\beta + \left( {q^\beta + \beta \gamma } \right)\Lambda _1 M\,Sin\theta } \over {q^\beta + \beta \gamma + \beta _1 q^\beta \gamma }}} }} \cr & \left( {{{\Lambda _2 Gr\,Cos\delta } \over {\Lambda _o q}}{1 \over {{{\Lambda _3 Pr} \over {\Lambda _4 }}\left( {{{q^\beta \gamma } \over {q^\beta + \beta \gamma }}} \right) - {{\Lambda _O \Lambda _1 \gamma q^\beta + \left( {q^\beta + \beta \gamma } \right)\Lambda _1 M\,Sin\theta } \over {q^\beta + \beta \gamma + \beta _1 q^\beta \gamma }}}}\left( {1 + } \right.} \right. \cr & \left. {\left. {b\sqrt {{{\Lambda _3 Pr} \over {\Lambda _4 }}\left( {{{q^\beta \gamma } \over {q^\beta + \beta \gamma }}} \right)} } \right) + {{q^2 - 1} \over {(1 + q^2 )^2 }}} \right)e^{ - \xi \sqrt {{{\Lambda _O \Lambda _1 \gamma q^\beta + \left( {q^\beta + \beta \gamma } \right)\Lambda _1 M\,Sin\theta } \over {q^\beta + \beta \gamma + \beta _1 q^\beta \gamma }}} } \cr & - {{\Lambda _2 Gr\,Cos\delta } \over {\Lambda _o q}}{1 \over {{{\Lambda _3 Pr} \over {\Lambda _4 }}\left( {{{q^\beta \gamma } \over {q^\beta + \beta \gamma }}} \right) - {{\Lambda _O \Lambda _1 \gamma q^\beta + \left( {q^\beta + \beta \gamma } \right)\Lambda _1 M\,Sin\theta } \over {q^\beta + \beta \gamma + \beta _1 q^\beta \gamma }}}}e^{ - \xi \sqrt {{{\Lambda _3 Pr} \over {\Lambda _4 }}\left( {{{q^\beta \gamma } \over {q^\beta + \beta \gamma }}} \right)} } \cr} and (52) $\begin{array}{l} \bar{w} (ξ, q) = \frac{1}{1 + b \sqrt{\frac{Λ_{O} Λ_{1} γ q + (q + α γ) Λ_{1} M S i n θ}{q + α γ + β_{1} q γ}}} \\ (\frac{Λ_{2} G r C o s δ}{Λ_{o} q} \frac{1}{\frac{Λ_{3} P r}{Λ_{4}} (\frac{q γ}{q + α γ}) - \frac{Λ_{O} Λ_{1} γ q + (q + α γ) Λ_{1} M S i n θ}{q + α γ + β_{1} q γ}} (1 + \\ b \sqrt{\frac{Λ_{3} P r}{Λ_{4}} (\frac{q γ}{q + α γ})}) + \frac{q^{2} - 1}{{(1 + q^{2})}^{2}}) e^{- ξ \sqrt{\frac{Λ_{O} Λ_{1} γ q + (q + α γ) Λ_{1} M S i n θ}{q + α γ + β_{1} q γ}}} \\ - \frac{Λ_{2} G r C o s δ}{Λ_{o} q} \frac{1}{\frac{Λ_{3} P r}{Λ_{4}} (\frac{q γ}{q + α γ}) - \frac{Λ_{O} Λ_{1} γ q + (q + α γ) Λ_{1} M S i n θ}{q + α γ + β_{1} q γ}} e^{- ξ \sqrt{\frac{Λ_{3} P r}{Λ_{4}} (\frac{q γ}{q + α γ})}} . \end{array}$ \eqalign{ & \bar w\left( {\xi ,q} \right) = {1 \over {1 + b\sqrt {{{\Lambda _O \Lambda _1 \gamma q + \left( {q + \alpha \gamma } \right)\Lambda _1 M\,Sin\theta } \over {q + \alpha \gamma + \beta _1 q\gamma }}} }} \cr & \left( {{{\Lambda _2 Gr\,Cos\delta } \over {\Lambda _o q}}{1 \over {{{\Lambda _3 Pr} \over {\Lambda _4 }}\left( {{{q\gamma } \over {q + \alpha \gamma }}} \right) - {{\Lambda _O \Lambda _1 \gamma q + \left( {q + \alpha \gamma } \right)\Lambda _1 M\,Sin\theta } \over {q + \alpha \gamma + \beta _1 q\gamma }}}}\left( {1 + } \right.} \right. \cr & \left. {\left. {b\sqrt {{{\Lambda _3 Pr} \over {\Lambda _4 }}\left( {{{q\gamma } \over {q + \alpha \gamma }}} \right)} } \right) + {{q^2 - 1} \over {(1 + q^2 )^2 }}} \right)e^{ - \xi \sqrt {{{\Lambda _O \Lambda _1 \gamma q + \left( {q + \alpha \gamma } \right)\Lambda _1 M\,Sin\theta } \over {q + \alpha \gamma + \beta _1 q\gamma }}} } \cr & - {{\Lambda _2 Gr\,Cos\delta } \over {\Lambda _o q}}{1 \over {{{\Lambda _3 Pr} \over {\Lambda _4 }}\left( {{{q\gamma } \over {q + \alpha \gamma }}} \right) - {{\Lambda _O \Lambda _1 \gamma q + \left( {q + \alpha \gamma } \right)\Lambda _1 M\,Sin\theta } \over {q + \alpha \gamma + \beta _1 q\gamma }}}}e^{ - \xi \sqrt {{{\Lambda _3 Pr} \over {\Lambda _4 }}\left( {{{q\gamma } \over {q + \alpha \gamma }}} \right)} } . \cr}

Case 4: g(t)=t e^t

In the final case, we take g(t) = te^t with its Laplace $G (q) = \frac{1}{{(q - 1)}^{2}}$ G\left( q \right) = {1 \over {{{(q - 1)}^2}}} , then the expressions of velocity through AB and CF derivative along with Eqs. (30), and (42) correspondingly will take the form as (53) $\begin{array}{l} \bar{w} (ξ, q) = \frac{1}{1 + b \sqrt{\frac{Λ_{O} Λ_{1} γ q^{β} + (q^{β} + β γ) Λ_{1} M S i n θ}{q^{β} + β γ + β_{1} q^{β} γ}}} \\ (\frac{Λ_{2} G r C o s δ}{Λ_{o} q} \frac{1}{\frac{Λ_{3} P r}{Λ_{4}} (\frac{q^{β} γ}{q^{β} + β γ}) - \frac{Λ_{O} Λ_{1} γ q^{β} + (q^{β} + β γ) Λ_{1} M S i n θ}{q^{β} + β γ + β_{1} q^{β} γ}} (1 + \\ b \sqrt{\frac{Λ_{3} P r}{Λ_{4}} (\frac{q^{β} γ}{q^{β} + β γ})}) + \frac{1}{{(q - 1)}^{2}}) e^{- ξ \sqrt{\frac{Λ_{O} Λ_{1} γ q^{β} + (q^{β} + β γ) Λ_{1} M S i n θ}{q^{β} + β γ + β_{1} q^{β} γ}}} \\ - \frac{Λ_{2} G r C o s δ}{Λ_{o} q} \frac{1}{\frac{Λ_{3} P r}{Λ_{4}} (\frac{q^{β} γ}{q^{β} + β γ}) - \frac{Λ_{O} Λ_{1} γ q^{β} + (q^{β} + β γ) Λ_{1} M S i n θ}{q^{β} + β γ + β_{1} q^{β} γ}} e^{- ξ \sqrt{\frac{Λ_{3} P r}{Λ_{4}} (\frac{q^{β} γ}{q^{β} + β γ})}} \end{array}$ \eqalign{ & \bar w\left( {\xi ,q} \right) = {1 \over {1 + b\sqrt {{{{\Lambda _O}{\Lambda _1}\gamma {q^\beta } + \left( {{q^\beta } + \beta \gamma } \right){\Lambda _1}M\,Sin\theta } \over {{q^\beta } + \beta \gamma + {\beta _1}{q^\beta }\gamma }}} }} \cr & \left( {{{{\Lambda _2}Gr\,Cos\delta } \over {{\Lambda _o}q}}{1 \over {{{{\Lambda _3}Pr} \over {{\Lambda _4}}}\left( {{{{q^\beta }\gamma } \over {{q^\beta } + \beta \gamma }}} \right) - {{{\Lambda _O}{\Lambda _1}\gamma {q^\beta } + \left( {{q^\beta } + \beta \gamma } \right){\Lambda _1}M\,Sin\theta } \over {{q^\beta } + \beta \gamma + {\beta _1}{q^\beta }\gamma }}}}\left( {1 + } \right.} \right. \cr & \left. {\left. {b\sqrt {{{{\Lambda _3}Pr} \over {{\Lambda _4}}}\left( {{{{q^\beta }\gamma } \over {{q^\beta } + \beta \gamma }}} \right)} } \right) + {1 \over {{{(q - 1)}^2}}}} \right){e^{ - \xi \sqrt {{{{\Lambda _O}{\Lambda _1}\gamma {q^\beta } + \left( {{q^\beta } + \beta \gamma } \right){\Lambda _1}M\,Sin\theta } \over {{q^\beta } + \beta \gamma + {\beta _1}{q^\beta }\gamma }}} }} \cr & - {{{\Lambda _2}Gr\,Cos\delta } \over {{\Lambda _o}q}}{1 \over {{{{\Lambda _3}Pr} \over {{\Lambda _4}}}\left( {{{{q^\beta }\gamma } \over {{q^\beta } + \beta \gamma }}} \right) - {{{\Lambda _O}{\Lambda _1}\gamma {q^\beta } + \left( {{q^\beta } + \beta \gamma } \right){\Lambda _1}M\,Sin\theta } \over {{q^\beta } + \beta \gamma + {\beta _1}{q^\beta }\gamma }}}}{e^{ - \xi \sqrt {{{{\Lambda _3}Pr} \over {{\Lambda _4}}}\left( {{{{q^\beta }\gamma } \over {{q^\beta } + \beta \gamma }}} \right)} }} \cr} and (54) $\begin{array}{l} \bar{w} (ξ, q) = \frac{1}{1 + b \sqrt{\frac{Λ_{O} Λ_{1} γ q + (q + α γ) Λ_{1} M S i n θ}{q + α γ + β_{1} q γ}}} \\ (\frac{Λ_{2} G r C o s δ}{Λ_{o} q} \frac{1}{\frac{Λ_{3} P r}{Λ_{4}} (\frac{q γ}{q + α γ}) - \frac{Λ_{O} Λ_{1} γ q + (q + α γ) Λ_{1} M S i n θ}{q + α γ + β_{1} q γ}} (1 + \\ b \sqrt{\frac{Λ_{3} P r}{Λ_{4}} (\frac{q γ}{q + α γ})}) + \frac{1}{{(q - 1)}^{2}}) e^{- ξ \sqrt{\frac{Λ_{O} Λ_{1} γ q + (q + α γ) Λ_{1} M S i n θ}{q + α γ + β_{1} q γ}}} \\ - \frac{Λ_{2} G r C o s δ}{Λ_{o} q} \frac{1}{\frac{Λ_{3} P r}{Λ_{4}} (\frac{q γ}{q + α γ}) - \frac{Λ_{O} Λ_{1} γ q + (q + α γ) Λ_{1} M S i n θ}{q + α γ + β_{1} q γ}} e^{- ξ \sqrt{\frac{Λ_{3} P r}{Λ_{4}} (\frac{q γ}{q + α γ})}} . \end{array}$ \eqalign{ & \bar w\left( {\xi ,\;q} \right) = {1 \over {1 + b\sqrt {{{{\Lambda _O}{\Lambda _1}\gamma q + \left( {q + \alpha \gamma } \right){\Lambda _1}M\,Sin\theta } \over {q + \alpha \gamma + {\beta _1}q\gamma }}} }} \cr & \left( {{{{\Lambda _2}Gr\,Cos\delta } \over {{\Lambda _o}q}}{1 \over {{{{\Lambda _3}Pr} \over {{\Lambda _4}}}\left( {{{q\gamma } \over {q + \alpha \gamma }}} \right) - {{{\Lambda _O}{\Lambda _1}\gamma q + \left( {q + \alpha \gamma } \right){\Lambda _1}M\,Sin\theta } \over {q + \alpha \gamma + {\beta _1}q\gamma }}}}\left( {1 + } \right.} \right. \cr & \left. {\left. {b\sqrt {{{{\Lambda _3}Pr} \over {{\Lambda _4}}}\left( {{{q\gamma } \over {q + \alpha \gamma }}} \right)} } \right) + {1 \over {{{(q - 1)}^2}}}} \right){e^{ - \xi \sqrt {{{{\Lambda _O}{\Lambda _1}\gamma q + \left( {q + \alpha \gamma } \right){\Lambda _1}M\,Sin\theta } \over {q + \alpha \gamma + {\beta _1}q\gamma }}} }} \cr & - {{{\Lambda _2}Gr\,Cos\delta } \over {{\Lambda _o}q}}{1 \over {{{{\Lambda _3}Pr} \over {{\Lambda _4}}}\left( {{{q\gamma } \over {q + \alpha \gamma }}} \right) - {{{\Lambda _O}{\Lambda _1}\gamma q + \left( {q + \alpha \gamma } \right){\Lambda _1}M\,Sin\theta } \over {q + \alpha \gamma + {\beta _1}q\gamma }}}}{e^{ - \xi \sqrt {{{{\Lambda _3}Pr} \over {{\Lambda _4}}}\left( {{{q\gamma } \over {q + \alpha \gamma }}} \right)} }}. \cr}

6.

RESULTS AND DISCUSSION

The mixed convection fluid flow with Al₂O₃ and TiO₂ nanoparticles with water and blood as base fluids are investigated with a slip effect at the boundaries on an inclined plane under the magnetic field by using AB and CF derivative schemes. The solution for the considered problem is explored with the Laplace scheme and numerical approaches i.e., Stehfest and Tzou for the inversion phenomenon of the Laplace transform. To obtain some physical features of velocity attained with AB and CF derivatives, some particular cases are also discussed. For studying the impacts of diverse flow parameters i.e., fractional parameters (α, β), angle of inclination, magnetic parameter, volume fraction, Pr, and Grashof number, graphical diagrams are presented in Figs 2–9 through Mathematica.

The impact of fractional parameters (α, β), on temperature is shown in Figs. 2(a, b). By raising the estimations of α, β, the temperature illustrates decaying behaviour (at a small time) and an increasing trend at a large time, consequently, we see that the fractional parameters have dual behaviour (for small and large time) for temperature. This specifies the consequence of the CF and AB fractional operators that promise to illustrate the generalized memory and hereditary features. This is due to the different properties of CF (based on exponential function having no singularity) and AB (having non-singular and non-local kernel) fractional operators. From Fig. 3a, as improvement in the estimations of Pr shows that development in the viscosity of liquid declines the difference among thermal boundary layers of the liquid, so the temperature profile declines because of the rises in the estimations of Pr and in the same way, the comparison of two nanofluids is shown in Fig. 3b by considering other parameters constant and changing the fractional parameters α, β. We see that the temperature of the blood-based NF is smaller than the water-based nanofluid, which is due to the physical characteristics of certain nanoparticles.

From Figs. 4a and 4b, we see that the velocity of fluid also decelerates by enhancing the estimation of α, β for a small time but speeds up at a large time. This is also due to the diverse properties of CF (based on exponential function having no singularity) and AB (having non-singular and non-local kernel) fractional operators. The fluid velocity is increased as Gr grows as shown in 5a. The velocity increases because of enhancing the values Gr. The relative impact of the heat buoyant behaviour on the viscous force is investigated using Gr. Such effects happen because of the existence of buoyant forces. An enhancement in the Pr declines the fluid velocity due to development in the fluid viscosity as in 5b Pr is inversely related to thermal diffusivity, resulting in declining heat transmission. The effect of volume fraction (φ) on velocity is illustrated in Fig 6a with the variation in time. The increase in volume fraction enhances the viscous impact of fluid flow which slows down the velocity of the fluid.

In Figs. 7a and 7b, growing the estimation of the magnetic parameter slows down the fluid velocity. Larger estimates of M lead to decreased velocity. Physical applications for such observations are because of the Lorentz force which produces a resistance in the flow. Similarly, the behaviour of the inclination angle of the magnetic field is shown in Fig 7b. The growth in the inclination angle declines the influence of the magnetic field which conveys off the Lorentz force effect, so by growing the estimation of the inclination angle, the fluid velocity again decreases. For θ = π/2 (normal magnetic field), the velocity is maximum declined, such observations are because of the Lorentz force which produces resistance in the flow. The Lorentz force has the greatest influence on velocity, decreasing velocity.

The comparison of ordinary and fractional fluid velocity is illustrated in Fig. 8a and 8b for diverse values of fractional parameters. We observe that when the fractional parameter i.e., α, β → 1, the fractional fluid velocity almost overlaps with ordinary velocity, which represents the convergence of our obtained numerical solutions of the velocity profile. From the graphical illustration, we see that the results attained by the AB-fractional derivative show more growing behaviour than the CF-fractional derivative. This is also due to the different properties of CF- (based on exponential function having no singularity) and AB (having non-singular and non-local kernel) fractional operators.

The comparison of different nanofluids for velocity field is shown in Fig. 9a. We see that the enhancement in velocity, due to (H₂O-Al₂O₃) and (Blood-Al₂O₃ ) is more advanced, than (H₂O-TiO₂) and (Blood-TiO₂) based NFs but (H₂O-Al₂O₃) based NF has a higher velocity than (Blood-Al₂O₃ ) based NF, all this behaviour is due to the physical characteristics of certain nanoparticles. The comparison of numerical methods specifically Grave Stehfest as well as Tzou’s is considered in Fig. 9b. The curves of the Stehfests, as well as Tzou’s scheme, overlap each other in both cases, which also validates our present results. Furthermore, to make the validity of our attained solutions, the numerical comparison of temperature as well as velocity field through Stehfest and Tzou’s with Nusselt number as well as skin friction, are presented in Tables 2–3.

7.

CONCLUSIONS

We study mixed convection flows over an inclined plate with Al₂O₃ and TiO₂ nanoparticles with water and blood-based fluids along with new definitions of CF (based on exponential function having no singularity) and AB (having non-singular and non-local kernel) fractional operators in several circumstances which is a significant theoretical and practical study for the solution of important problems. A semi-analytical approach for AB and CF-based models is applied by the Laplace transform technique along with Stehfest and Tzou’s numerical schemes.

−

The temperature shows dual behaviour with different estimations of fractional parameters with diverse estimations (small and large) of the time.

−

The temperature displays decaying behaviour for large estimations of the Pr.

−

The velocity slows down by growing the estimation of M.

−

The velocity profile speeds up as increasing the estimations of Gr and declines for increasing values of volume fraction φ.

−

The enhancement in velocity, due to (H₂O-Al₂O₃) and (Blood-Al₂O₃ ) is more advanced, than (H₂O-TiO₂ ) and (Blood-TiO₂ ) based NFs.

−

Our obtained solutions through different numerical methods specifically Stehfest and Tzou’s are alike.

Consequently, our claimed results offer significant insights into industrial and engineering systems. These findings guide the development of thermal transfer technologies, assisting in the optimization of processes for better efficiency in applications such as cooling mechanisms and power generation. The research advances heat transfer processes, increasing the overall efficiency of industrial systems.

Langue:: Anglais

Périodicité:: 4 fois par an
Sujets de la revue:: Ingénierie, Ingénierie électrique, Électronique, Génie mécanique, Mécanique, Génie biologique et génie biomédicale, Biomécanique, Travaux publiques, Génie de l'environnement

RSS Feed de la revue

Dual Thermal Analysis of Fractional Convective Flow Through Aluminum Oxide and Titanium Dioxide Nanoparticles

Qasim Ali

Rajai S. Alassar

Irfan. A. Abro

Kashif. A. Abro

Publié en ligne: 31 mars 2025

Pages: 32 - 43

Reçu: 05 janv. 2024

Accepté: 11 mars 2024

DOI: https://doi.org/10.2478/ama-2025-0005

Mots clésNon-singularized derivative, Mixed convection flow with nanoparticles, Heat transfer of inclined plate, Integral transforms

© 2025 Qasim Ali et al., published by Sciendo

This work is licensed under the Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License.

Mots clés
Non-singularized derivative, Mixed convection flow with nanoparticles, Heat transfer of inclined plate, Integral transforms