Heat transfer of generalized second grade fluid with MHD, radiation and exponential heating using Caputo

Scientific Reports volume 13, Article number: 5220 (2023) Cite this article

427 Accesses

Metrics details

The aim of the present work is to apply the Caputo–Fabrizio fractional derivative to the heat transformation of unsteady incompressible second grade fluid. The effects of magneto hydro dynamic and radiation are analyzed. In governing equation of heat transfer nonlinear radiative heat is examined. Exponential heating phenomena is considered at boundary. Firstly, the dimensional governing equations with the initial & boundary conditions are converted into non-dimensional form. Exact analytical solutions are obtained for dimensionless fractional governing equations which consist of momentum and energy equations by using Laplace transform method. Special cases are investigated of the obtained solutions and it is noticed that some well-known results are achieved published in literature from these special cases. At the end, for graphical illustration the influences of different physical parameters like radiation, Prandtl, fractional parameter, Grashof numbers and Magneto hydro dynamic are checked graphically.

The theory of derivatives with fractional order has great importance in daily life. As integer order, the theory of non-integer order is also oldest. It is the branch of mathematics, a few years ago this concept was limited only in mathematics, but now a days the principles of fractional calculus have been often carried out to different fields such as fluid dynamics, bio engineering, electromagnetism, fluid mechanics, finance, electrochemistry, viscoelasticity, in biology the model of neurons, applied mathematics1. In fluid dynamics the non-integer derivative concept have been used to investigate viscoelastic process like polymers in the glassy state and glassy transition2. A few years ago, fractional order derivatives have been seen is an effective tool from which a suitable physical concepts generalization can be gained. There are so many other definitions of derivatives with non-integer order but the Caputo fractional and Riemann-Liouvilli fractional derivatives are used in different real world phenomena's3,4. Everyone knows that such methods show difficulties in application. Such as the derivative of a constant is non-zero in Riemann-Liouvilli fractional order derivative and also it has singular kernel. These difficulties are removed by Caputo and gave the concept in which the constant has zero derivative but have still singular kernel. After all these, Fabrizio & Caputo presented the idea of non-integer order derivative in which constant has derivative zero & without singular kernel. By Laplace technique Caputo-Febrizio fractional derivative is easy to find exact solution. Many existing fluid models have been examined and fractional order derivative has been developed. Some of well-known fluid models are presented here like, the Oldroyd-B, Maxwell, grade second, Burger and Jeffery fluid models etc. The Burger, Maxwell and Oldroyd models are rate type models, while the grade second are of differential type5. According to Tan et al.6 investigated the generalized unsteady flow of grade second non-Newtonian fluid between two parallel plates with the model of non- integer derivatives. Recently, Friedrich7, examined the fluid model of ordinary Maxwell fluid with fractional order derivative generalized the function of relaxation and retardation. In the earlier studies, Tan et al.8 analyzed a short note on non-integer Maxwell fluid with flow of unsteady viscoelastic fluid between two parallel plates. The non-integer viscoelastic Maxwell fluid model with one directional periodic fluid flow studied in9. The model of fractional Maxwell fluid of viscoelastic in pipe was examined by Yin et al.10. Brikman type fluid by Caputo fractional derivative is investigated in11. The effects of parameters in generalized second grade fluid is discussed in12. The Maxwell non-integer order derivative for the Stokes first problem is studied in13. Khan et al.14 studied the generalized modified Darcy's law with Oldroyd-B fluid to obtain exact solutions for Magnetohydrodynamic. Khan et al.15 studied Burgers fluid model of viscoelastic non-integer on accelerated flows. By using Caputo Fabrizio non-integer derivative studied heat transfer fluid of second grade over and oscillating perpendicular surface examined in16. Heat mass transfer investigated in the third grade fluid with chemical reaction upon a stretchable sheet fixed in a porous medium. Abbas et al.17 investigated the third grade fluid thermal diffusion with Darcy–Forchheimer relation upon a stretchable sheet. Analysis of heat transfer in Atangana–Baleanu derivative to Newtonian heating and convection flows of Caputo–Fabrizio with second grade fluids is investigated in18. Recently, by using derivative of Caputo Fabrizio non-integer and examined the exponential heating and Magnetohydrodynamic flow of second grade fluid in19. Saqib et al.20 studied Jeffery fluid flow by using Caputo–Fabrizio derivative and obtain exact solutions. Raptis et al.21 investigated MHD influence of thermal radiation over a stretchable sheet influence of heat radiation on MHD is studied in22. Purpose of this article is to discuss analysis of generalized second grade non-newtonian fluid on magneto hydrodynamic and heat radiation using Caputo–Fabrizio fractional derivative approach. On thermal aspect, exponential heating phenomena to be adopted.

Consider the incompressible non-Newtonian second grade fluid. Initially for time t = 0 the temperature T∞ and velocity is zero. As time start for t = 0+, the fluid velocity becomes \(fH(t)e^{i\omega t}\), here H(t) is the unit step function and temperature reaches \(T_{\infty } + T_{\omega } (1 - ae^{ - bt} )\). According to all these assumptions temperature and velocity both are the function of space variable "y" and time "t" only. Now by usual Boussinesq's approximation16, the unsteady flow is governed by the following set of partial differential equations. The schematic diagram used in fluid flow problem is represented geometrically by Fig. 1.

Geometry of the problem.

For radiation approximation of Rosseland is used23, we have

By neglecting the higher terms with the help of Taylor series we express T4 as a linear function,

Initial and boundary condition:

Dimensionless Variable:

After dimensionless Eqs. (1)–(5), we get

Now using the Caputo–Fabrizio time derivative in Eqs. (7) and (8) we get the system given below:

To obtain the solutions of the governing equations we use the Laplace transform technique. First we find the solution of energy equations because momentum equation depends on it. Now taking the Laplace transform of Eq. (12) with initial and boundary conditions Eqs. (9) and (10), we have the following equations:

Solving Eq. (15) with help of Eq. (16) we a transformed solution which is given

Now to find the exact analytical solution of energy equation we take the inverse Laplace transform of Eq. (17) by using the Appendices A1 and A2, we obtain the solution which is presented in Eq. (18)

To find the velocity equation we take the Laplace Transform of Eq. (12) with initial boundary conditions Eqs. (9) and (10), now the transformed equation with transformed initial boundary conditions are given in Eqs. (19) and (20):

Solving Eq. (19) with help of Eq. (20) we get the transformed which is presented in Eq. (21)

where; \(a_{1} = \frac{M + \gamma }{{1 + \alpha_{2} \gamma }},a_{2} = \frac{M\alpha \gamma }{{1 + \alpha_{2} \gamma }},a_{3} = \frac{\alpha \gamma }{{1 + \alpha_{2} \gamma }},h_{1} = - \frac{Gr\xi }{{\Pr \gamma + \Pr \gamma^{2} \alpha_{2} - M\xi - \gamma \xi }},\)

To obtain the exact analytical solution of momentum equation we take the inverse Laplace transform of Eq. (21) we get the solution which is given in Eq. (22) by using the Appendices A1, A2 and A4,

(i) In the absence of radiation effect \(N = 0\) and neglecting the exponential heating of the plate.

In Eq. (8), when we put \(N = 0\), then we obtain the solution in the form as given under:

where \(\varphi \left( {y, \, t, \, p_{r} \gamma , \, \alpha \gamma } \right) \) is defined by the Appendix (A1).

The result is uniform to that in the published literature achieved by Shah and Khan16.

Also neglecting the radiation effect in Eq. (7) we get the solution for velocity equation which is given as under:

where \(a_{1} = \frac{\gamma }{{1 + \alpha_{2} \gamma }}, \, a_{2} = \, \alpha a_{1}\) and \(\varphi \left( {y, \tau ,a_{1} ,a_{2} } \right)\) is defined by the Appendix (A1).

The result is uniform to the one in published literature achieved by Shah and Khan16.

By using Mathcad software different physical parameters are sketch to analyze the effects of fluid velocity and temperature. The parameter alpha α in Fig. 2, Prandtl number Pr in Fig. 3 and heat radiation in Fig. 4 are drawn for temperature field, while for velocity field alpha α in Fig. 5, Prandtl number Pr in Fig. 6, Magneto hydro dynamic MHD in Fig. 7 and Grashof number Gr in Fig. 8 are presented.

Graph of temperature for different value of alpha α.

Graph of temperature for different value of Prandtl number Pr.

Graph of temperature for different value of radiation N.

Graph of velocity for different value of alpha α in case of cosine and sine oscillation.

Graph of velocity for different value of Prandtl number Pr in case of cosine and sine oscillation.

Graph of velocity for different value of M in case of cosine and sine oscillation.

Graph of velocity for different value of Grashof number Gr in case of cosine and sine oscillation.

Figure 2 is sketch to check the effects of temperature and alpha α in which we saw this temperature is rises by growing the value of α, the boundary layer thermal thickness is rises with parameter alpha α and time t. Figure 3 is sketch to check the influence of temperature and Prandtl Pr in which we observed that the temperature decreases by rising value of Prandtl Pr, the thickness of thermal boundary layer is decreasing with the parameter Prandtl number Pr and time t and diffusivity of temperature is large. Figure 4 is sketch to check the effects of temperature and heat radiation N it has been investigated, by growing the small value of heat radiation N the temperature is also rises. The graph is plotted temperature versus y. Figure 5 is drawn to check the influence of alpha α, both cases of sine and cosine oscillation are discussed in which we investigate this fluid velocity decreases by increasing the value of alpha α. This graph shows both the effects of cosine and sine oscillation for fluid. Effects of sine oscillation are greater than cosine oscillation by increasing the time t. Figure 6 is drawn to examine Prandtl Pr effects over fluid velocity, individually the cases of sine & cosine oscillation were considered in which we emphasized this by growing the small values of Prandtl Pr number, velocity decreases. Figure 7 is drawn to study the behaviour of Magnetohydrodynamic M, both the cases of sine and cosine oscillation are considered in which we noticed by small values of MHD increasing the velocity has decreases. Effects of sine oscillation are greater than cosine oscillation by increasing the time t. Figure 8 is drawn to observe the influence of Grashof number, both cases of sine and cosine oscillation are considered in which we noticed that the fluid velocity is increases by growing the value of Gr. The effects of sine oscillation are greater than cosine oscillation by increasing the time t. In Fig. 9 we compared the obtained solutions as limiting cases with those obtained by Shah and khan16.

Velocity and temperature profile with comparison of published literature published by Shah & Khan.

Numerical results of skin-friction and Nusselt number at the plate \(\left( {y = 0} \right)\) are presented in Tables 1 and 2 for different values of \(\left( t \right)\), \(\left( \alpha \right)\), \(\left( N \right)\), \(\left( M \right)\),\(\left( {\Pr } \right)\) and \(\left( {Gr} \right)\). It is observed from Table 1 that skin friction \(\left( \tau \right)\) increases with an increase in \(\left( t \right)\), \(\left( N \right)\) and \(\left( {Gr} \right)\) while the result is reversed with increase in \(\left( \alpha \right)\), \(\left( M \right)\) and \(\left( {\Pr } \right)\) in Table 1. Numerical results of Nusselt number \(\left( {Nu} \right)\) at the plate \(\left( {y = 0} \right)\) are expressed in Tables 2 for different values of \(\left( t \right)\), \(\left( \alpha \right)\), \(\left( N \right)\) and \(\left( {\Pr } \right)\). Table 2 shows that the Nusselt number Nu which determines the rate of heat transfer at the plate increases as \(\left( \alpha \right)\) and \(\left( {\Pr } \right)\) progresses while the result is reversed with increase in \(\left( t \right)\) and \(\left( N \right)\).

Unsteady free convection flow of generalized second grade fluid over an infinite vertical plate is studied. The flow is analyzed under the effect of magneto hydrodynamic and radiation together with heat transfer. Moreover On the thermal aspects of the infinite vertical plate, we are taking into account the exponential heating phenomena. The Caputo–Fabrizio derivative has been applied to the set of dimensionless governing equation. Exact solution of the problem is obtained through Laplace transorm technique.The profiles (temperature and velocity) are analyzed graphically for both sine and cosine oscillations of the plate for distinct physical parameters.

It is observed that.

By increasing the fractional parameter α and radiation N, temperature is also increases.

With the increases of Prandtl number the temperature can be decreased.

Velocity is decreases by increasing α parameter and hence velocity and temperature have opposite behavior for α parameter.

With large value of Prandtl number the fluid velocity tends to decreasing.

The motion of the fluid is decreasing for growing value of MHD.

The velocity is growing by large value of Gr.

The datasets analyzed during the current study available from the corresponding author on reasonable request.

Specific heat

Gravity (acceleration)

Parameter of radiation

Grashof number

Wall temperature

Viscosity (kinematic)

Fluid density

Electrical conductivity

Fluid temperature

Prandtl number

Coefficient of mean absorption

Constant (Stefan–Boltzmann)

Grade second parameter

Uniform magnetic field

Thermal expansion volumetric coefficient

Time

Laplace transform parameter

Thermal conductivity

Kulish, V. V. & Lage, J. L. Application of fractional calculus to fluid mechanics. J. Fluids Eng. 124(3), 803–806 (2002).

Article Google Scholar

Debnath, L. Recent applications of fractional calculus to science and engineering. Int. J. Math. Math. Sci. 2003(54), 3413–3442 (2003).

Article MathSciNet MATH Google Scholar

Hilfer, R. Threefold introduction to fractional derivatives. Anomal. Transport Found. Appl. 23, 17–73 (2008).

Article Google Scholar

Gorenflo, R., Mainardi, F., Moretti, D. & Paradisi, P. Time fractional diffusion: A discrete random walk approach. Nonlin. Dyn. 29(1), 129–143 (2002).

Article MathSciNet MATH Google Scholar

Caputo, M. & Fabrizio, M. A new definition of fractional derivative without singular kernel. Progr. Fract. Differ. Appl. 1(2), 1–13 (2015).

Google Scholar

Wenchang, T. & Mingyu, X. Unsteady flows of a generalized second grade fluid with the fractional derivative model between two parallel plates. Acta Mech. Sin. 20(5), 471–476 (2004).

Article MathSciNet Google Scholar

Friedrich, C. H. R. Relaxation and retardation functions of the Maxwell model with fractional derivatives. Rheol. Acta 30(2), 151–158 (1991).

Article CAS Google Scholar

Wenchang, T., Wenxiao, P. & Mingyu, X. A note on unsteady flows of a viscoelastic fluid with the fractional Maxwell model between two parallel plates. Int. J. Non-Linear Mech. 38(5), 645–650 (2003).

Article ADS MATH Google Scholar

Hayat, T., Nadeem, S. & Asghar, S. Periodic unidirectional flows of a viscoelastic fluid with the fractional Maxwell model. Appl. Math. Comput. 151(1), 153–161 (2004).

MathSciNet MATH Google Scholar

Yin, Y. & Zhu, K. Q. Oscillating flow of a viscoelastic fluid in a pipe with the fractional Maxwell model. Appl. Math. Comput. 173(1), 231–242 (2006).

MathSciNet MATH Google Scholar

Sene, N. Analytical investigations of the fractional free convection flow of Brinkman type fluid described by the Caputo fractional derivative. Res. Phys. 37, 105555 (2022).

Google Scholar

Yavuz, M., Sene, N. & Yıldız, M. Analysis of the influences of parameters in the fractional second-grade fluid dynamics. J. Math. 10(7), 1125 (2022).

Article Google Scholar

Jamil, M., Rauf, A., Zafar, A. A. & Khan, N. A. New exact analytical solutions for Stokes’ first problem of Maxwell fluid with fractional derivative approach. Comput. Math. Appl. 62(3), 1013–1023 (2011).

Article MathSciNet MATH Google Scholar

Khan, M., Hayat, T. & Asghar, S. Exact solution for MHD flow of a generalized Oldroyd-B fluid with modified Darcy's law. Int. J. Eng. Sci. 44(5–6), 333–339 (2006).

Article MathSciNet MATH Google Scholar

Khan, M., Ali, S. H. & Qi, H. On accelerated flows of a viscoelastic fluid with the fractional Burgers’ model. Nonlinear Anal. Real World Appl. 10(4), 2286–2296 (2009).

Article MathSciNet MATH Google Scholar

Shah, N. A. & Khan, I. Heat transfer analysis in a second grade fluid over and oscillating vertical plate using fractional Caputo–Fabrizio derivatives. Eur. Phys. J. C 76(7), 1–11 (2016).

Article Google Scholar

Abbas, A., Shafqat, R., Jeelani, M. B. & Alharthi, N. H. Convective heat and mass transfer in third-grade fluid with Darcy-Forchheimer relation in the presence of thermal-diffusion and diffusion-thermo effects over an exponentially inclined stretching sheet surrounded by a porous medium. Adv. CFD Convect. Heat Transf. 10(4), 776–791 (2022).

CAS Google Scholar

Siddique, I., Tlili, I., Bukhari, S. M. & Mahsud, Y. Heat transfer analysis in convective flows of fractional second grade fluids with Caputo–Fabrizio and Atangana-Baleanu derivative subject to Newtonion heating. Mech. Time-Dependent Mater. 25(3), 291–311 (2021).

Article ADS Google Scholar

Haq, S. U., Shah, S. I. A., Jan, S. U. & Khan, I. MHD flow of generalized second grade fluid with modified Darcy's law and exponential heating using fractional Caputo–Fabrizio derivatives. Alex. Eng. J. 60(4), 3845–3854 (2021).

Article Google Scholar

Saqib, M., Ali, F., Khan, I., Sheikh, N. A. & Jan, S. A. A. Exact solutions for free convection flow of generalized Jeffrey fluid: A Caputo–Fabrizio fractional model. Alex. Eng. J. 57(3), 1849–1858 (2018).

Article Google Scholar

Raptis, A., Perdikis, C. & Takhar, H. S. Effect of thermal radiation on MHD flow. Appl. Math. Comput. 153(3), 645–649 (2004).

MathSciNet MATH Google Scholar

Siddheshwar, P. G. & Mahabaleswar, U. S. Effects of radiation and heat source on MHD flow of a viscoelastic liquid and heat transfer over a stretching sheet. Int. J. Non-Linear Mech. 40(6), 807–820 (2005).

Article ADS MATH Google Scholar

Hayat, T. & Qasim, M. Radiation and magnetic field effects on the unsteady mixed convection flow of a second grade fluid over a vertical stretching sheet. Int. J. Numer. Methods Fluids 66(7), 820–832 (2011).

Article ADS MathSciNet MATH Google Scholar

Download references

Shaheed Benazir Bhutto Women University Peshawar, Peshawar, 25000, Khyber Pakhtunkhwa, Pakistan

Sehra Sehra & Afshan Noor

Department of Mathematics, Islamia College Peshawar, Peshawar, 25000, Khyber Pakhtunkhwa, Pakistan

Sami Ul Haq & Saeed Ullah Jan

Department of Mathematics, College of Science Al-Zulfi, Majmaah University, Al-Majmaah, 11952, Saudi Arabia

Ilyas Khan

Research Centre, Future University in Egypt, New Cairo, 11835, Egypt

Abdullah Mohamed

You can also search for this author in PubMed Google Scholar

S.S. designed the study; A.N. conducted the experiments with technical assistance from S.H., S.U.J and I.K. analyzed the data and wrote the paper. A.M. computed special case results with discussion and revised manuscript.

Correspondence to Ilyas Khan.

The authors declare no competing interests.

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and Permissions

Sehra, S., Noor, A., Haq, S.U. et al. Heat transfer of generalized second grade fluid with MHD, radiation and exponential heating using Caputo–Fabrizio fractional derivatives approach. Sci Rep 13, 5220 (2023). https://doi.org/10.1038/s41598-022-22665-4

Download citation

Received: 15 May 2022

Accepted: 18 October 2022

Published: 30 March 2023

DOI: https://doi.org/10.1038/s41598-022-22665-4

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.