Influence of Cattaneo-Christov model on Darcy-Forchheimer flow of Micropolar Ferrofluid over a stretching/shrinking sheet

https://doi.org/10.1016/j.icheatmasstransfer.2019.104385Get rights and content

Abstract

This work deals with the study of Darcy-Forchheimer flow of micropolar ferrofluid on a porous and dynamic (stretching/shrinking) sheet under the influence of thermal radiations subjected to both suction and injection. The effects of the external electric and magnetic fields are considered as well. Water is used as a base fluid and Fe3O4 (iron oxide) as electro-magnetite nanoparticles. The mathematical equations developed in this study are based on the Cattaneo-Christov model consisting of coupled nonlinear partial differential equations. These equations are transformed into a set of coupled ordinary differential equations (ODEs) by using similarity transformations. These ODEs are solved by applying the standard mathematical technique of homotopy analysis (HAM). The effects produced by different parameters on the velocity, micro-rotational velocity and temperature profiles are shown graphically for positive as well as negative mass transfer flow and for both stretching and shrinking cases. It has been observed that the velocity profile increases (decreases) with the increasing electric field strength and microrotation parameter during the stretching (shrinking) of the surface in both suction (S > 0) and injection (S < 0) cases. Furthermore, similar results has been observed for the velocity profile with the increasing inertial coefficient, porosity, magnetic and boundary parameters during flow over the stretching (shrinking) surface for both S > 0 and S < 0. The micro-rotational velocity increases with higher values of microrotation parameter for stretching, while decreases for the shrinking of the surface. The temperature profile displays an increasing trend with the increasing values of heat energy source and sink terms and thermal radiation parameter for stretching as well as shrinking of the sheet for S > 0 as well as S < 0. The temperature profile also changes with the variation in thermal relaxation parameter and Prandlt number.

Introduction

There are numerous scientific challenges for effective heat transfer in various fields like medicine, food, electronics, chemical sensors, solar cells, fuel cells, batteries, etc. These challenges can be handled by applying the recently developed field of nanotechnology. One of the major achievements of nanotechnology is the introduction of nanofluids. Nanofluid is a wonderful addition of Modern science and Technology. Due to nanometer dimensions, the nanofluid can flow through micro-size channels. Because of the enhanced convection between the base liquid surfaces and the nanoparticles, the nano-suspensions exhibit large thermal conduction. The suspension of the nanoparticles in the base fluid increases the heat conduction capacity, thermal conductivity, effective surface area, interparticles collisions, and the various interactions between the constituent particles. The nanofluid utilization for the first time was introduced by Choi in 1995 [1] for enhancing the thermal conductivity of the classical fluid. After its introduction, the nanofluid research has attracted many investigators throughout the world due its fascinating thermal properties and possible applications in different fields like medicine, electronics, petrochemical, refining, heating, air-conditioning and transportation. The experimental investigation of two different kinds of nanoparticles γ-alumina (Al2O3) and titaniumdioxide (TiO2) on the turbulent heat transfer flow was given by Pak and Cho [2]. It has been observed that the addition of the suspended nanoparticles caused an increase in the convective heat transfer coefficient of water. Eastman et al. [3] has claimed that the heat transfer coefficient in Cu-water suspended nanoparticles is higher than as compared with for pure water by analyzing some initial suspended nanoparticles experiments. Qiang and Yimin [4] have found better results than Eastman et al. [3] by experimentally studying the convective heat transfer characteristics of Cu- water based nanofluid flow. The coefficient of convective heat transfer of carbon nanotubes nanofluid in the presence of constant heat flux was experimentally studied by Rashidi and Nezamabad [5]. They found a considerable increase in the convective heat transfer coefficient of the carbon nanotubes nanofluid and also found it as function of axial distance from the inlet. Mahanta and Abramson [6] investigated the thermal conductivities of graphene and graphene oxide and obtained that the higher thermal conductivity of multilayer graphene platelets can be associated with the interlayer coupling due to covalent interactions provided by oxygen atoms. Sun et al. [7] have examined the convective heat transfer properties of water based on ferrofluid nanofluid inside copper tubes and showed that the addition of dispersants resulted a significant enhancement of the stability of the nanofluids. Walvekar et al. [8] have experimentally investigated the heat transfer enhancement during the turbulent flow of carbon nanotube (CNT) nanofluids by varying the CNT concentration and found that heat transfer enhances as a function of CNT concentration and temperature. Dawar et al. [9] have recently examined the radiative heat transfer and the MHD carbon nanotubes Casson nanofluid flow between two rotating plates and investigated the effects produced by various parameters of interest on the velocity and temperature profiles.

The mathematical models which researchers employ for studying different behaviours of nanofluids are Buongiorno and Tiwari-Das (TD) models [10,11]. Inertia, Magnus effect, diffusiophoresis, fluid drainage, gravity, thermophoresis, and Brownian diffusion, are the seven different mechanisms in which slip take place among the solid and fluid phases in the Buongiorno model. Buongiorno studied all these seven mechanisms and obtained that when turbulent effect is absent, the theromphoresis and the Brownian diffusion are the significant and leading mechanisms in the nanofluid dynamics. The Buongiorno model has been used in various studies [[12], [13], [14], [15], [16], [17]]. Later on, the Buongiorno model has been further extended by different researchers to investigate variety of problems related to the nanofluid flow on stretching as well as shrinking surfaces. For example, Zaimi et al. [18] studied numerically the boundary layer and heat transfer flow of a nanofluid over a permeable stretching (shrinking) surface. By analyzing the effects of various parameters, they found that the suction widens the range of stretching (shrinking) parameter for which the solution exists. Hayat et al. [19] has analytically studied the three dimensional MHD nanofluid flow on a stretching surface in the presence of heat generation (absorption) and convective heat transfer by considering the effects of Brownian motion and thermophoresis. They found that the temperature distribution varies directly with Biot number and magnetic parameter. Muhammad et al. [20] has investigated the MHD boundary layer flow of Maxwell nanofluid on a porous stretching surface based on Darcy-Forchheimer model and studied the temperature and concentrations profiles by varying the different relevant parameters. Muhammad et al. [21] studied the three dimensional Darcy-Forchheimer nanofluid flow on a bidirectional stretching surface by including the thermophoretic and Brownian diffusion effects in the energy and volume fraction of nanoparticles expressions. They explained the results of their work through different graphs and tables. Hayat et al. [22] examined the MHD flow of electrically conducting nanofluid over a porous stretching surface by taking into account the Brownian motion and thermophoretic diffusion of nanoparticles and studied graphically the impacts produced by different physical parameters on the velocity, temperature and nanoparticle concentration. Hayat et al. [23] undertook the three dimensional MHD flow of a couple stress electrically conducting nanofluid on a bidirectional stretching surface in the presence of thermophoresis and Brownian motion effects. They have graphically shown the effects produced by different emerging parameters on the temperature and nanoparticles concentration and concluded that the thickness of the thermal boundary layer is an increasing function of radiative effect. Further relevant study can be found in Hayat et al. [24,25]. In contrast to the Buongiorno model which mainly concentrates on the the thermophoresis effects and Brownian motion, the TD model considers the nanoparticles volume fraction. Under different physical conditions, many researchers use TD model to study the different heat transfer characteristics of the nanofluids. The convective heat transfer of Ag and Cu-water base nanofluids on a stretching (shrinking) surface during steady flow of a boundary layer has been studied by Yacob et al. [26]. They found that the heat energy transfer rate decreases with increasing convective parameter while increases with higher values of nanoparticle volume fraction. It was also observed that heat transfer rate at the Cu-water nanofluid surface is larger than at the surface of Ag-water nanofluid, although Ag is a good thermal conductor as compared to Cu. Furthermore, Hussanan et al. [27] has examined the magnetite ferrofluid flow on a stretching (shrinking) surface subjected to suction and injection in order to investigate the micro-structure and inertial properties of the substructure of iron oxide (Fe3O4) nanoparticles mixed with water (base fluid). During this study, they observed that micro-rotation velocity changes remarkably with positive and negative values of mass transfer parameter and also concluded that micropolar ferrofluid has higher velocity as compared to classical nanofluid. The TD model has further been applied to other relevant research and can be found in the references [[28], [29], [30], [31], [32], [33], [34]].

The most important model for the transfer of heat due to conduction is the Fourier's law [35]. Cattaneo [36] in 1948 recommended a thermal relaxation time factor to the Fourier's model to improve the heat transfer. But the main difficulty is the different thermal relaxation times of the different materials. Keeping this fact in mind, a very useful time derivative model was developed for the effective heat transfer mechanism by Christov [37], which is called the Cattaneo-Christov heat flux model. Meraj Mustafa [38] applied it to study the viscoelastic rotating fluid flow bounded by a stretching surface. He found an inverse relation between velocity and viscoelastic fluid parameter and also observed that fluid temperature varies inversely with the relaxation time and with the Prandtl number. The heat flux model of Cattaneo-Christov was also applied by Hayat et al. [39] in order to study the boundary layer flow on a dynamic surface of variable thickness and to explore the various aspects of heat transformation with varying thermal conductivity. They analyzed and explain the dependence of temperature and velocity distributions on different pertinent parameters. Ali and Sandeep [40] have used Cattaneo-Christov model and numerically studied the heat transfer due to radiations of the MHD Casson-ferrofluid. The different studies where Cattaneo-Christov model has been used in [[41], [42], [43]].

In this work, we examine the Cattaneo-Christov ferrofluid model on a porous and stretching (shrinking) surface by taking into account the effects of thermal radiations subjected to suction and injection. The effects of electric and magnetic fields are considered as well. In this brief study, water is taken as base fluid mixed with Fe3O4 (iron oxide) electro-magnetite nanoparticles. The mathematical equations are developed by using the Cattaneo-Christov model of nanofluid. We solve these model equations through the standard homtopy analysis procedure by using the similarity variables. In Section 2, the mathematical model consisting of different equations is outlined. Section 3 is devoted to the solution of these models equations by homotopy analysis. The results and discussion are presented in Section 4. We have concluded our work in Section 5.

Section snippets

Problem statement

We consider the boundary layer two-dimensional flow of a micropolar ferrofluid on a linear porous and stretching (shrinking) sheet. The x-axis is taken parallel while the y-axis is taken perpendicular to the sheet. Both magnetic and electric fields are assumed along the y-axis. The stretching (shrinking) velocity of the sheet is taken as Uw(x) = αx, where α > 0 (<0) for the stretching (shrinking) of the sheet. Due to the assumed geometry and conditions, the system of equations which govern the

Homotopy analysis method

In order to solve Eqs. (11), (12), (13) by applying the boundary (14), (15), we use the standard mathematical technique (HAM) [49,50]. With the auxiliary parameter ℏ, we adjust and control the solution convergence.

The initial assumptions are taken in the following form [51]:f0ξ=αeξ+S+α,g0ξ=αδeξ,θ0ξ=1eξ

The operators Lf, Lg and Lθ, which are linear in nature are defined as:Lff=d3fdξ3df,Lgg=d2gdξ2g,Lθθ=d2θdξ2θhaving the following propertiesLfε1+ε2eξ+ε3eξ=0,Lgε4+ε5eξ=0,Lθε6eξ+ε7eξ=0where ε

Results and discussions

We present graphically here the impact of the various parameters on f'(ξ) (velocity profile), g(ξ) (micro-rotational velocity profile), and θ(ξ) (temperature profile). To see the variations produced by different embedded parameters on these three profiles for the stretching and shrinking (α > 0, α < 0) cases, Fig. 3, Fig. 4, Fig. 5, Fig. 6, Fig. 7, Fig. 8, Fig. 9, Fig. 10, Fig. 11, Fig. 12, Fig. 13, Fig. 14, Fig. 15, Fig. 16, Fig. 17, Fig. 18, Fig. 19, Fig. 20, Fig. 21, Fig. 22, Fig. 23, Fig. 24

Conclusions

In the current study, we have investigated the Cattaneo-Christov model of electromagnetite ferrofluid on a porous stretching (shrinking) sheet under the influence of the thermal radiation subjected to both suction and injection. The developed mathematical equations are tackled with the help of similarity variables through HAM. We have investigated the effects produced by different parameters on the velocity, micro-rotation velcoity and temperature functions. The most important findings of the

Declaration of Competing Interest

The author declares that they have no competing interests.

Acknowledgments

This project was funded by the Deanship of Scientific Research (DSR), King Abdulaziz University, Jeddah, under grant No. (DF-866-130-1441). The authors, therefore, gratefully acknowledge DSR for technical and financial support.

References (53)

  • T. Hayat et al.

    On darcy-forchheimer flow of carbon nanotubes due to a rotating disk

    Int. J. Heat Mass Transf.

    (2017)
  • M. Sheikholeslami et al.

    MHD forced convection flow of nanofluid in a porous cavity with hot elliptic obstacle by means of lattice boltzmann method

    Int. J. Mech. Sci.

    (2018)
  • A. Hussanan et al.

    Microstructure and inertial characteristics of a magnetite ferrofluid over a stretching/shrinking sheet using effective thermal conductivity model

    J. Mol. Liq.

    (2018)
  • K. Vajravelu et al.

    Convective heat transfer in the flow of viscous ag–water and cu–water nanofluids over a stretching surface

    Int. J. Therm. Sci.

    (2011)
  • M. Hamad

    Analytical solution of natural convection flow of a nanofluid over a linearly stretching sheet in the presence of magnetic field

    Int. Commun. Heat Mass Tran.

    (2011)
  • M. Hamad et al.

    Magnetic field effects on free convection flow of a nanofluid past a vertical semi-infinite flat plate

    Nonlinear Anal.

    (2011)
  • M. Sheikholeslami et al.

    Analytical investigation of NHD nanofluid flow in a semi-porous channel

    Powder Technol.

    (2013)
  • M. Sheikholeslami et al.

    Heat transfer of Cu-water nanofluid flow between parallel plates

    Powder Technol.

    (2013)
  • R. Kandasamy et al.

    Impact of chemical reaction on cu, al2o3 and swcnts–nanofluid flow under slip conditions

    Eng. Sci. Technol. Int. J.

    (2016)
  • C. Christov

    On frame indifferent formulation of the Maxwell–cattaneo model of finite-speed heat conduction

    Mech. Res. Commun.

    (2009)
  • M. Ali et al.

    Cattaneo-christov model for radiative heat transfer of magnetohydrodynamic casson-ferrofluid: a numerical study

    Results Phys.

    (2017)
  • M. Waqas et al.

    Cattaneo-christov heat flux model for flow of variable thermal conductivity generalized burgers fluid

    J. Mol. Liq.

    (2016)
  • T. Hayat et al.

    Three-dimensional flow of nanofluid with cattaneo–christov double diffusion

    Results Phys.

    (2016)
  • M. Turkyilmazoglu

    Flow of a micropolar fluid due to a porous stretching sheet and heat transfer

    Int. J. Non-Linear Mech.

    (2016)
  • S.U. Choi et al.

    Enhancing Thermal Conductivity of Fluids With Nanoparticles

    (1995)
  • B.C. Pak et al.

    Hydrodynamic and heat transfer study of dispersed fluids with submicron metallic oxide particles

    Exp. Heat Transf. Int. J.

    (1998)
  • Cited by (66)

    • Dynamics of magneto-electric hybrid nanoparticles with chemically reacting and radiated moving plate: Entropy analysis

      2022, International Communications in Heat and Mass Transfer
      Citation Excerpt :

      They emphasized that the understanding the basics of how energy moves through nanofluids is important for making nanofluids that are very energy efficient for a variety of heat transfer applications. Subsequently, other researchers [2–6] considered stretching surface as a geometry and studied various nanofluids under varying boundary conditions. Maleki et al. [7–9] discussed different nanofluid flows over porous surface in the presence/absence of slip parameters.

    View all citing articles on Scopus
    View full text