Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Darcy-Forchheimer nanofluidic flow manifested with Cattaneo-Christov theory of heat and mass flux over non-linearly stretching surface

Abstract

This research article aims to disclose the features of nanofluidic flow manifested with Cattaneo-Christov model of heat and mass flux over non-linearly stretching surface. An incompressible visco-elastic nanofluid saturates the given porous medium through Darcy-Forchheimer relation. A non-uniformly induced magnetic effect is considered to accentuate the electro-magnetic and thermal conductivity of the base fluid. The model is restricted to small magnetic Reynolds. Boundary layer assumptions are incorporated for the given flow model. Governing equations are remodeled into non-linear ordinary differential equations through transformations. So formulated nonlinear system is solved through homotopy analysis method (HAM) to achieve series solutions for velocity field, concentration of nanoparticles and temperature distribution. It is noticed that the temperature distribution and corresponding thermal boundary layer pattern shows declination for Cattaneo-Christov model of heat and mass flux as compared to classical Fourier’s law of heat flux/conduction. Furthermore, the intensive resistance offered by the addition of porosity factor in the flow model results in rise of temperature profile, however, opposite behavior is noticed in concentration of nanoparticles. The wall-drag intensity, the heat flux and the mass flux are discussed on the premise of numerical information obtained upon simulation of the problem.

1 Introduction

Due to considerable importance of nanofluids in industrial, chemical, environmental, geological and may other setups, the concept of suspension of metallic nano-size particles in typical base fluids presented by Choi [1], termed as nanofluids, has become a subject of immense interest for numerous researchers especially physicists for the past many years. The repute of nano-science in present era enforces us to interpret the behavior of fluids and nanofluids over various surfaces including stretching sheets, cylinders and plates etc. to analyze the drag force variation, heat and mass flux mechanism and other important industrial aspects (see for example [2] [3] [4] [5] [6] [7] [8] [9] [10] [11] [12] [13] [14] [15]) and cross references cited therein where authors have focused on the heat and mass flux attributes of nanofluids. Industrial and environmental setups that are likely to involve the heat and mass flux such as geological setups, chemical actor-reactors and other heating systems are the main focus in convection through porous medium. The porosity factor involved in such systems is merely to induce the incremental inertial effect that results in variation of skin friction (wall-drag). Such formulations are helpful in building heat insulation related materials, energy storage setups, fossil fuel beds, disposal reactors for nuclear wastes and many others. The concept was initially introduced by Darcy and the law, known as Darcy’s law, is valid for a restricted range of velocity. In the case of porous media, it becomes almost impossible to ignore the high effects of inertial force and resistance offered by the porosity factor. Thus, a squared velocity term within the context of Darcian velocity is involved in the momentum equation to analyze the inertial effects. Numerous articles have been reported on this subject explaining various influential parameters. For example, Hayat et al. [16] reported Darcy model flow using Cattaneo-Christov theory of heat and mass flux over a linearly stretching sheet. Reportedly, the momentum of the fluid and the associated boundary layer resulted a declination for augmented values of porosity factor. In another article, Bakar et al. [17] disclosed the fluid flow analysis at stagnation point. The fluid momentum and the thermal boundary layer resulted in reducing behavior for augmented porosity factor. Taseer et al. [18] has introduced a revised model for Maxwell nanofluid using Darcy Forchheimer relation. Similar results have been reported in this study as compared to [17]. Alshomrani et al. [19] involved a convectively heated stretching surface manifested with Darcy flow model to analyze the impact of porosity and enhanced frictional force on fluid flow. Mandal and Seth [20] reported a Casson type nanofluid flow with Darcy-Forchheimer relation. The magnetic effect resulted in obvious declination of the fluid momentum. Hayat et al. [21] implemented Darcy model for curvy stretching sheet witnessing an increase in thermal layer for incremental nature of porosity factor. Recently, Sheikholeslami et al. [22] reported the impact of MHD (Lorentz forces) on Fe3O4-H2O Ferro-fluid flow using permeable semi-annulus. The porosity factor resulted in higher resistance offered to the fluid motion following a low velocity field. In another article, Sheikholeslami et al. [23] analyzed the behavior of water based nanofluid in a porous media/enclosure with quite similar results. Dogonchi et al. [24] reported a CuH2O based nanofluid flow through porous media between a hot and cold cylindrical framework. The impact of thermal radiation and porosity factor in CuOH2O based nanofluid has been analyzed by Dogonchi et al. [25]. A novel correlation for the averaged Nusselt number has been presented in this study. Quite similar but more effective results have been reported in the study presented by Dogonchi et al. [26] involving an annulus subject to thermal radiation.

Recently, the engagement of Cattaneo-Christov model in nanofluids to analyze the heat and mass flux mechanism is trending as numerous researchers are working on such formulations. Heat flux is a natural process that occurs where there is a temperature difference within a system or within the systems. Usually, the law of thermodynamics and heat conduction reported by Fourier [27] has been extensively used for many years to analyze the heat and mass characteristics however, it restricts the energy and concentration equations to parabolic type equations that means an initial disturbance would lead to instant experience by the system, called a paradox of heat and mass flux. This restriction was removed by Cattaneo [28] by enforcing a modification with relaxation time factor. This term, therefore, covers/overcomes the heat and mass flux paradox. Later on, Christov [29] further improved the version of model presented by Cattaneo by replacing time derivative with Oldroyd upper-convective derivative. The theory, thus, termed as Cattaneo-Christov heat and mass flux theory (CC-model). Related to this discussion, Haddad [30] analyzed thermal instability using CC-model in Brinkman porous medium. Li et al. [31] utilized CC-model in heat and mass flux analysis of visco-elastic fluid flow over stretching surface subject to slip conditions. A decreasing velocity profile is noticed for enhanced magnetic effect however, the result is opposite for the thermal profile. Sui et al. [32] disclosed the influence of CC-model on Maxwell nanofluidic flow through stretching surface subject to slip conditions. Relatively small impact of relaxation parameters is noticed on velocity profile as compared to thermal and salute layers. Ganji and Dogonchi [33] reported a squeezing flow of nanofluids confined in parallel plates involving CC-model and thermal radiation. The results indicated that thermal distribution is lesser for CC-model as compared to the classical Fourier’s law. Upadhay et al. [34] discussed the impact of CC-model on heat and mass transfer attributes in Powell-Eyring type nanofluidic flow. The thermal relaxation time parameter resulted in an improvement in heat flux rate. Acharya et al. [35] reported impact of CC-model in a magnetized upper convective Maxwell type nanofluidic flow past an inclined stretching sheet. The study was conducted in the context of generalized Fourier and Ficks perspectives.

The theory of fluid flow, heat and mass transfer analysis over a stretching sheet/surface and porous media is an important aspect. It’s numerous applications in engineering procedures like paper production, aluminous plate cooling procedures, plastic sheet extrusions etc. are remarkably best-known to the readers. This immense interest in the involvement of stretching sheet in fluid flow resulted in numerous research studies reported in the last two decades. For example, Sajid et al. [36] mentioned viscous flow of steady incompressible fluid instigated by curvy extended stretching sheet. It is reported that the wall-drag force (skin-friction) on a curvy surface is lesser than a flat surface. The studies mentioned above apprehend the argument that more focus has been emphasized on linear stretching surfaces and non-linearity in stretching rate has been neglected upto a decent level though, the stretching with nonlinear pattern is reasonably important in practical applications of fluids/nanofluids. The concept is new maneuver and trending recently. Therefore, keeping in view the context of non-linearity in stretching rates, recently, Rasool et al. [37] presented MHD nanofluidic flow over non-linearly stretching surface/sheet reporting some interesting jumps in the heat and mass flux rates. An innovative study on MHD nanofluidic flow over vertically stretching surface/sheet has been reported by Alarifi et al. [38]. The model involved heat source to balance the base temperature of the surface. A report on second grade nanofluid flow past a porous media has be presented by Khan et al. [39]. The study incorporated stretching surface/sheet, heat and mass transfer attributes. In a similar study, Khan et al. [40] reported an analysis on heat-mass flux attributes as well as thermophoresis and Brownian diffusion in MHD nanofluid flow over stretching surface. Palwasha et al. [41] reported a non-Newtonian nanofluid flow through porous media using microorganisms magneto-tactics. Some recent studies can be seen in ([42] [43] [44] [45] [46]) and cross references cited therein.

The MHD (Magnetohydrodynamics) is another very important aspect in engineering, chemical and environmental setups involving fluid flow analysis. Instruments such as bearings, generators, chemical actor/reactors, pumps and many others are tormented by MHD. The applications of MHD are also found in industrial systems where a high speed machine is working having tiny size. The temperature range is witnessed between zero to 350-degrees on the scale of Celcius. Over the years, the concept of MHD has been involved in fluid flow analysis to help improve the thermal and electro-magnetic conductivity of the subjected fluid/nanofluid. Some related studies have been reported in this text. For example, Singh et al. [47] presented an MHD flow using variable thermal conductivity. The analysis involved stagnation point formulation. MHD convective heat-mass transfer under Soret-Dufour effects has been reported by Chatterjee et al. [48] using Power-Law model and porosity. Lund et al. [49] reported a study on the dual solution for MHD Williamson fluid flow with slippage. Khan et al. [50] reported MHD mixed convection in second grade nanofluidic flow considering Browniian diffusion and thermophoresis together with Hall effects over a stretching sheet/surface. Khan et al. [51] involved MHD in nanoliquid thin film flow in a cylinder. Zuhra et al. [52] discussed the heat-mass flux attributes in second grade MHD nanofluid flow saturated with gyrotactic micro-organisms and nanoparticles. The influence of inclined MHD on heat and mass transfer attributes in Carreau nanoliquid flow has been reported by Khan et al. [53]. MHD CuH2O based natural convection in complex shaped enclosure has been reported by Dogonchi et al. [54]. A similar study on copper-water nanofluidic flow through a horizontal semi-cylinder has been reported by Dogonchi et al. [55].

In all the studies mentioned above, the major emphasis has been given to linearity in stretching rates. Here in this research, the non-linearity in stretching rate of the surface has been targeted to interpret various aspects of fluid flow analysis. Furthermore, no such study is found in literature involving non-linear stretching surface, Darcy-Forchheimer relation and Cattaneo-Christov model of heat-mass flux all together. The article is organized as follows: Firstly, a Darcy-Forchheimer relation along with Cattaneo-Christov model of heat and mass flux is implemented on steady, incompressible and viscoelastic nanofluid flow bounded by a flat non-linearly stretching sheet/surface. Boundary layer assumptions, Brownian diffusion and thermophoresis are attended. Secondly, the governing equations are remodeled into non-linear ordinary differential equations through transformations. Thirdly, HAM [56] is used to get the final series solutions. Fourthly, graphs are plotted to investigate the variation in velocity field, temperature distribution and concentration of the nanoparticles. Finally, a correlation is given for elevated values of different parameters to help the audience in understanding the variation in skin-friction (wall-drag), heat and mass flux rates.

2 Problem formulation

Consider a viscoelastic incompressible nanofluidic flow manifested with Cattaneo-Christov model of heat and mass flux over non-linearly stretching surface. The nanofluid saturates the given porous medium through Darcy-Forchheimer relation. A non-uniformly induced magnetic effect is involved to accentuate the electro-magnetic and thermal conductivity of the base fluid. The model is restricted to small Reynolds to dismiss the influence of induced magnetic effect. Cattaneo-Christov model is employed to involve the effect of modified Fourier’s law. The stretching velocity is u = Uw = axn where n = 1 implies linear rate of stretching and n > 1 implies non-linearity in stretching of the sheet. The stretching rate is assumed to be non-linear for the given flow model. The sheet extends along x–direction while y–direction is taken surface normal to it. There is no fluid movement along the y–axis and the steady flow is assumed along x–axis only. Schematic can be seen in Fig 1.

The governing equations are as follows, (1) (2) (3) (4)

According to Cattaneo-Charistov heat flux theory, see for example [35], we have, (5) and (6) (5) & (6) are the equations for heat and mass flux with relaxation time λ1 for thermal and λ2 for solute levels. Classical Fourier’s law can be deduced by setting λ1 = λ2 = 0 in the above mentioned equations. Considering the natural incompressibility i.e. and assumption of steady flow i.e. in (5-6) yields, (7) and, (8) The respective equations of energy and concentration of nanoparticles are, therefore, (9) (10) where, (11) and (12) with boundary conditions, (13) (14) Defining, (15)

Using (15) in (1)–(2) & (9)–(10), we have (16) (17) (18) (19) Here, γi for i = 1, 2 is the relaxation time parameter for temperature and concentration of nanoparticles. Fr is local inertial force, M is magnetic effect, λ is known for porosity, Pr is known as Prandtl factor, Le is known as Lewis factor, k1 is known as viscoelastic factor, Nb is called Brownian diffusion factor and Nt is called thermophoretic force factor. Mathematically, (20) Quantities of physical and industrial interest such as wall drag (Cf), heat flux (local Nusselt) (Nu) and mass flux (local Sherwood) (Sh) numbers are given by: (21) where qw, τw and qm are the surface heat flux, shear stress and surface mass flux, respectively. Using the definitions of τw, qw and qm and simplifying, (22) Therefore, the non-dimensional forms are given below: (23) where Rex = axn+1/ν is the local Reynolds number.

3 Methodology

Homotopy analysis method (HAM) is an efficient solver for nonlinear systems where small and large physical parameters are involved for analysis of variation in the respective flow profiles. Thus, it is more suitable and economical as compared to that of perturbation method for solving non-linear systems of equations. It is applied in most of non-linear systems developed in engineering, industrial, science and finance problems. Define, (24) (25) such that, (26) where Ki, i = 1, 2, ⋯, 7, are constants. Subsequently, the zeroth order deformation problems: (27) (28) (29) correspond to, (30) with boundary conditions, (31) where and are auxiliary-parameters and z ∈ [0, 1] is known as embedding factor. and are non-linear operators. The Taylor’s series expansion results in, (32) where for or . The speedy and smooth convergence of the series solutions is strictly dependent on choice . For z = 0, 1, (33) Consequently, the ith deformations: (34) (35) (36) correspond to, (37) with boundary conditions, (38) where Φi = 1 for i > 1, otherwise 0. Finally,

The general solutions are, (39) where and are special solutions.

4 Convergence analysis

The auxiliary parameters involved in series solutions for the velocity field (f′), temperature distribution (θ) and concentration of the nanoparticles (ϕ) for the problem under consideration are termed as convergence control parameters. These parameters are critical in choosing appropriate values to speed-up the convergence. The intervals of interest for f, θ and ϕ are presented in Fig 2. One can see that the intervals of convergence are [-1.30, -0.10], [-1.50, -0.10] and [-1.50, -0.10], respectively. Data upto 40th order approximations is listed in Table 1. 15th order of approximations are sufficient to achieve convergence in velocity field whereas temperature and concentration requires 10th order approximations.

5 Results & discussion

The given non-linear system of Eqs (16)–(18) with boundary conditions (19) is solved by HAM for series solutions. Numerical data for wall-drag (skin-friction) coefficient (Cf), local Nusselt (heat flux) (Nu) and local Sherwood (mass flux) (Sh) is compiled in Tables 2 & 3, respectively. Tables 4 & 5 are the results of comparison/validation of skin friction, local Nusselt and local Sherwood with Rasool et al. [37], respectively. Table 6 is the correlation of the skin-friction while Table 7 is the correlation of local Nusselt and local Sherwood for pertinent fluid parameters. In distinction, a reducing behavior in skin friction (Cf) and mass flux (Sh) is noted for visco-elastic nanofluid parameter (k1). Furthermore, the Brownian diffusion rate (Nb) and thermophoresis effect (Nt) appeared as reducing factors for heat flux (Nu). Influence of various fluid parameters on the flow profiles i.e. velocity field, temperature distribution and concentration of the nanoparticles in base conventional fluid, is shown in Figs 313. The influences of porosity and inertia on the non-dimensional velocity f′(η) and corresponding variations are plotted in Figs 3 and 4, respectively. The influence of porosity parameter λ presents a decreasing behavior in the respective profile plotted in Fig 3. Physically, the existence of porous factor results in increment of resistance offered by medium to the fluid motion that causes declination in the fluid momentum and connected boundary layer thickness reduces. Fig 4 shows the change in velocity profile for various incremental values of inertial coefficient. Similar to the porosity factor, reduction in flow profile is witnessed with incremental values of Fr. A similar trend is noticed in velocity profile for various values of visco-elastic nano-fluid parameter k1 plotted in Fig 5. Physically, the incremental values of k1 result in increasing the fluid viscosity that reduces the flow motion. Figs 611. The visco-elastic nanofluid parameter results in increasing thermal field and the associated boundary layer shows increasing thickness. The results are plotted in Fig 6. The relaxation time parameter for the Temperature γ1 reduces the temperature profile and the associated thickness of thermal boundary layer as well. Physically, it confirms that incremental nature of thermal relaxation time parameter requires more time to shift the heat from intensively packed fluid particles to the low energetic fluid particles. One can call it a demonstration of characteristics of non-conducting fluid material. Thus, a decay in temperature profile is noted. The results are plotted in Fig 7. Fig 8 shows the variation noted in thermal distribution and the associated boundary layer for incremental values of the induced non-uniform magnetic effect. An increasing trend is noticed in temperature profile for stronger magnetic effect. Physically, the retardation offered to the fluid motion by sudden jumps created by Magnetic field increase the particles’ collision which is responsible for increasing trend in temperature boundary layer. Figs 9 and 10 retrieve the increasing trend of Brownian diffusion parameter and thermophoretic force for associated thermal profile. Physically, the erratic motion of fluid packets appeared due to increasing trend in Brownian motion, thus an enhancing nature of temperature profile is noted. Further, the increase in thermophoretic force produces more intensive and vigorous thermophoretic influence causing the nano-particles to move away from the stretching sheet. This development induces boost in the temperature profile. Figs 1113 are the plots of variation in concentration of nanoparticles against the concentration relaxation time parameters, Thermophoretic force and Prandtl number. Fig 11 is specifically plotted for variation in concentration profile against the incremental values of relaxation parameter for concentration of the nanoparticles. A mixed trend is noted as plotted in the respective figure. Physically, the relaxation parameter allows sufficient time to the nanoparticles to dilute in the base fluid that results in an increasing trend with the passage of time. Fig 12 shows the behavior of concentration profile for incremental values of thermophoretic parameter. An increasing trend is noted for higher values of thermophoresis parameter. Physically, during the thermophoresis, nanoparticles are forced to move from hotter region to the colder region, therefore, the hot particles saturated near the sheet, start moving away from the sheet surface. This movement results in an augmentation in the concentration distribution. The un-deniable fact that the Prandtl number induces a decreasing impact on concentration of nanoparticles is seen in Fig 13. Physically, the incremental values of Prandtl number correspond to a weaker thermal diffusivity. Thus, the concentration of the nanoparticles reduces for higher Prandtl.

thumbnail
Fig 3. Variation in velocity field for incremental values of λ.

https://doi.org/10.1371/journal.pone.0221302.g003

thumbnail
Fig 4. Variation in velocity field for incremental values of Fr.

https://doi.org/10.1371/journal.pone.0221302.g004

thumbnail
Fig 5. Variation in velocity field for incremental values of k1.

https://doi.org/10.1371/journal.pone.0221302.g005

thumbnail
Fig 6. Variation in temperature field for incremental values of k1.

https://doi.org/10.1371/journal.pone.0221302.g006

thumbnail
Fig 7. Variation in velocity field for incremental values of γ1.

https://doi.org/10.1371/journal.pone.0221302.g007

thumbnail
Fig 8. Variation in temperature field for incremental values of M.

https://doi.org/10.1371/journal.pone.0221302.g008

thumbnail
Fig 9. Variation in temperature field for incremental values of Nt.

https://doi.org/10.1371/journal.pone.0221302.g009

thumbnail
Fig 10. Variation in temperature field for incremental values of Nb.

https://doi.org/10.1371/journal.pone.0221302.g010

thumbnail
Fig 11. Variation in concentration of nanoparticles for incremental values of γ2.

https://doi.org/10.1371/journal.pone.0221302.g011

thumbnail
Fig 12. Variation in concentration of nanoparticles for incremental values of Nt.

https://doi.org/10.1371/journal.pone.0221302.g012

thumbnail
Fig 13. Variation in concentration of nanoparticles for incremental values of Pr.

https://doi.org/10.1371/journal.pone.0221302.g013

thumbnail
Table 2. Numerical results/data of skin-friction for both linear and non-linear cases.

https://doi.org/10.1371/journal.pone.0221302.t002

thumbnail
Table 3. Numerical results/data of local Nusselt and local Sherwood for both non-linear (n = 1.2) and linear (n = 1) cases at Pr = 1, k1 = 0.2.

https://doi.org/10.1371/journal.pone.0221302.t003

thumbnail
Table 5. Comparison/validation of results with Rasool et al. [37] setting γ1 = γ2 = 0 = k1, n = 1.2.

https://doi.org/10.1371/journal.pone.0221302.t005

Conclusion

A locally similar analysis on Darcy Forchheimer visco-elastic nanofluid flow bounded by a non-linearly stretching sheet/surface manifested with Cattaneo-Christov theory of heat—mass flux has been carried out in this analytic research article. The key findings are itemized below:

  • Increasing values of interial coefficient and porosity result in declination of the velocity field and the associated momentum boundary layer.
  • The visco-elastic nanofluid parameter shows reduction in the velocity field however, an increment is noted in the thermal layer for augments in the aforementioned parameter.
  • The intensive resistance offered by the addition of porosity factor in the flow model results in rise of temperature profile, however, opposite behavior is noticed in concentration of nanoparticles.
  • Larger Prandtl number shows reduction in concentration of the nanoparticles.
  • Thermal relaxation time parameter allows more heat to be absorbed by the system. Thus, a decay is noted.
  • Solute relaxation time parameter shows mixed behavior in concentration of nanoparticles.
  • Skin-friction force is intensive for higher thermal relaxation time parameter.
  • Local Nusslt shows increasing behavior for higher values of thermal relaxation time parameter, however, the values are higher for non-linear case as compared with linear case.
  • Local Sherwood is increasing function of concentration relaxation time parameter, however, the values are lower for linear case as compared to the non-linear case.

Nomenclature:

References

  1. 1. Choi SUS. Enhancing thermal conductivity of fluids with nanoparticles. D. App. Non-Newtonian Flows. 1995; 231(66): 99–105.
  2. 2. Dogonchi AS, Chamkha AJ, Ganji DD. A numerical investigation of magneto-hydrodynamic naturalconvection of Cu–water nanofluid in a wavy cavity using CVFEM. Journal of Thermal Analysis and Calorimetry. 2019; 135: 2599–2611.
  3. 3. Dogonchi AS, Ismael MA, Chamkha AJ, Ganji DD. Numerical analysis of natural convection of Cu-water nanofluid filling triangular cavity with semi-circular bottom wall. Journal of Thermal Analysis and Calorimetry. 2019; 135: 3485–3497.
  4. 4. Dogonchi AS, Hashim . Heat transfer by natural convection of Fe3O4-water nanofluid in an annulus between a wavy circular cylinder and a rhombus. International Journal of Heat and Mass Transfer. 2019; 130: 320–332.
  5. 5. Zuhra S, Khan NS, Shah Z, Islam S, Bonyah E. Simulation of bioconvection in the suspension of second grade nanofluid containing nanoparticles and gyrotactic microorganisms. AIP Advances, 2018; 8: 105–210.
  6. 6. Khan NS, Zuhra S, Shah Z, Bonyah E, Khan W, Islam S. Slip flow of Eyring-Powell nanoliquid film containing graphene nanoparticles. AIP Advances. 2018; 8: 115–302.
  7. 7. Khan NS. Bioconvection in second grade nanofluid flow containing nanoparticles and gyrotactic microorganisms. Braz J Phys. 2018; 48: 227.
  8. 8. Mehryan SAM, Kashkooli FM, Soltani M, Raahemifar K. Fluid flow and heat transfer analysis of a nanofluid containing motile gyrotactic micro-organisms passing a nonlinear stretching vertical sheet in the presence of a non-uniform magnetic field; Numerical approach. PLoS ONE. 2016; 11(6): e0157598.
  9. 9. Hayat T, Aziz A, Muhammad T, Alsaedi A. A revised model for Jeffrey nanofluid subject to convective condition and heat generation/absorption. PLoS ONE. 2017; 12(2): e0172518. pmid:28231298
  10. 10. Ramzan M, Bilal M, Chung JD. Radiative flow of Powell-Eyring magneto-nanofluid over a stretching cylinder with chemical reaction and double stratification near a stagnation point. PLoS ONE. 2017; 12(1): e0170790. pmid:28129356
  11. 11. Rasool G, Zhang T, Shafiq A, Durur H. Influence of chemical reaction on Marangoni convective flow of nanoliquid in the presence of Lorentz forces and thermal radiation: A numerical investigation. J. Adv. Nanotechnology. 2019; 1(1): 32–49.
  12. 12. Rasool G, Zhang T, Shafiq A. Marangoni effect in second grade forced convective flow of water based nanofluid. J. Adv. Nanotechnology. 2019; 1(1): 50–61.
  13. 13. Ali I, Rasool G, Alrashed S. Numerical simulations of reaction–diffusion equations modeling prey–predator interaction with delay. Int. J. Biomath. 2018; 11(4): 1850054.
  14. 14. Rasool G, Zhang T, Shafiq A, Second grade nanofluidic flow past a convectively heated vertical Riga plate. Physica Scripta. 2019; Accepted.
  15. 15. Mustafa M, Mushtaq A, Hayat T, Alsaedi A. Rotating flow of magnetite-water nanofluid over a stretching surface inspired by non-linear thermal radiation. PLoS ONE. 2016; 11(2): e0149304. pmid:26894690
  16. 16. Hayat T, Muhammad T, Al-Mezal S, Liao SJ. Darcy-Forchheimer flow with variable thermal conductivity and Cattaneo-Christov heat flux. Int. J. Numer. Methods Heat Fluid Flow. 2016; 26: 2355–2369.
  17. 17. Bakar SA, Arifin NM, Nazar R, Ali FM, Pop I. Forced convection boundary layer stagnation-point flow in Darcy-Forchheimer porous medium past a shrinking sheet. Front. Heat Mass Transf. 2016; 7: 38.
  18. 18. Muhammad T, Alsaedi A, Shehzad SA, Hayat T. A revised model for Darcy-Forchheimer flow of Maxwell nanofluid subject to convective boundary condition. Chin. J. Phys. 2017; 55: 963–976.
  19. 19. Alshomrani AS, Ullah MZ. Effects of homogeneous-heterogeneous reactions and convective condition in Darcy-Forchheimer flow of carbon nanotubes. J. Heat Transf. 2019; 141: 012405.
  20. 20. Seth GS, Mandal PK. Hydromagnetic rotating flow of Casson fluid in Darcy-Forchheimer porous medium. In MATEC Web of Conferences; EDP Sciences: Les Ulis, France. 2018; 192: 02059.
  21. 21. Hayat T, Saif RS, Ellahi R, Muhammad T, Ahmad B. Numerical study for Darcy-Forchheimer flow due to a curved stretching surface with Cattaneo-Christov heat flux and homogeneous-heterogeneous reactions. Results Phys. 2017; 7: 2886–2892.
  22. 22. Sheikholeslami M, Arabkoohsar A, Khan I, Shafee A, Li Z. Impact of Lorentz forces on Fe3O4-water Ferro-fluid entropy and energy treatment within a permeable semi annulus. Journal of Cleaner Production. 2019; 221: 885–898.
  23. 23. Sheikholeslami M, Shah Z, Shafee A, Khan I, Tlili I. Uniform magnetic force impact on water based nanofuid thermal behavior in a porous enclosure with ellipse shaped obstacle. Scientific Reports. 2019; 9: 1196.
  24. 24. Dogonchi AS, Sheremet MA, Ganji DD, Pop I. Free convection of copper–water nanofluid in a porous gap between hot rectangular cylinder and cold circular cylinder under the effect of inclined magnetic field. Journal of Thermal Analysis and Calorimetry. 2019; 135(2): 1171–1184.
  25. 25. Dogonchi AS, Waqas M, Seyyedi SM. Numerical simulation for thermal radiation and porous medium characteristics in flow of CuOH2O nanofluid. J Braz. Soc. Mech. Sci. Eng. 2019; 41: 249.
  26. 26. Dogonchi AS, Waqas M, Seyyedi SM, Hashemi-Tilehnoee M, Ganji DD. CVFEM analysis for Fe3O4H2O nanofluid in an annulus subject to thermal radiation. Int. J. Heat and Mass Transfer. 2019; 132: 473–483.
  27. 27. Fourier JBJ, TheÂorie Analytique De La Chaleur. Paris 1822.
  28. 28. Cattaneo C. Sulla conduzione del calore. Atti Semin Mat Fis Univ Modena Reggio Emilia. 1948; 3: 83–101.
  29. 29. Christov CI. On frame indifferent formulation of the Maxwell-Cattaneo model of finite-speed heat conduction. Mech Res Commun. 2009; 36: 481–486.
  30. 30. Haddad . Thermal instability in Brinkman porous media with Cattaneo-Christov heat flux. Int J Heat Mass Transfer. 2014; 68: 659–668.
  31. 31. Li J, Zheng L, Liu L. MHD viscoelastic flow and heat transfer over a vertical stretching sheet with Cattaneo-Christov heat flux effects. J Mol Liq. 2016; 221: 19–25.
  32. 32. Sui J, Zheng L, Zhang X. Boundary layer heat and mass transfer with Cattaneo-Christov double-diffusion in upper-convected Maxwell nanofluid past a stretching sheet with slip velocity. Int J Ther Sci. 2016; 104: 461–8.
  33. 33. Dogonchi AS, Ganji DD. Impact of Cattaneo-Christov heat flux on MHD nanofluid flow and heat transfer between parallel plates considering thermal radiation effect. J Taiwan Inst Chem Eng. 2017; 80: 52–63.
  34. 34. Upadhay MS, Mahesha , Raju C. Cattaneo-Christov on heat and mass transfer of unsteady Eyring-Powell dusty nanofluid over sheet with heat and mass flux conditions. Inf Med Unlocked. 2017; 9: 76–85.
  35. 35. Acharya N, Das K, Kundu PK. Cattaneo-Christov intensity of magnetized upper-convected Maxwell nanofluid flow over an inclined stretching sheet: A generalised Fourier and Fick’s perspective. Int J Mech Sci. 2017; https://doi.org/10.1016/j.ijmecsci.2017.05.043.
  36. 36. Sajid M, Ali N, Javed T, Abbas Z. Stretching a curved surface in a viscous fluid. Chin. Phys. Lett. 2010; 27: 0247–03.
  37. 37. Rasool G, Shafiq A, Khalique CM, Zhang T. Magnetohydrodynamic Darcy Forchheimer nanofluid flow over nonlinear stretching sheet. Physica Scripta. 2019; Accepted.
  38. 38. Alarifi IM, Abokhalil AG, Osman M, Lund LA, Ayed MB, Belmabrouk H, et al. MHD flow and heat transfer over vertical stretching sheet with heat sink or source effect. Symmetry. 2019; 11(3): 297.
  39. 39. Khan NS, Islam S, Gul T, Khan I, Khan W, Ali L. Thin film flow of a second grade fluid in a porous medium past a stretching sheet with heat transfer. Alexandria Engineering Journal. 2017; 57(2).
  40. 40. Khattak NSK, Gul T, Islam S, Khan W. Thermophoresis and thermal radiation with heat and mass transfer in a magnetohydrodynamic thin-film second-grade fluid of variable properties past a stretching sheet. European Physical Journal Plus. 2017; 132(1).
  41. 41. Palwasha Z, Islam S, Khan NS, Ayaz H. Non-Newtonian nanoliquids thin-film flow through a porous medium with magnetotactic microorganisms. Applied Nanoscience. 2018; 8(6): 1523–1544.
  42. 42. Khan NS, Gul T, Khan MA, Bonyah E, Islam S. Mixed convection in gravity-driven thin film non-Newtonian nanofluids flow with gyrotactic microorganisms. Results in Physics. 2017; 7: 4033–4049.
  43. 43. Zuhra S, Khan NS, Khan MA, Islam S, Khan W, Bonyah E. Flow and heat transfer in water based liquid film fluids dispensed with graphene nanoparticles. Results in Physics. 2018; 8: 1143–1157.
  44. 44. Palwasha Z, Khan NS, Shah Z, Islam S, Bonyah E. Study of two-dimensional boundary layer thin film fluid flow with variable thermo-physical properties in three dimensions space. AIP Advances. 2018; 8: 105318.
  45. 45. Khan NS, Zuhra S, Shah Z, Bonyah E, Khan W, Islam S, et al. Hall current and thermophoresis effects on magnetohydrodynamic mixed convective heat and mass transfer thin film flow. J. Phys. Commun. 2019; 3: 035009.
  46. 46. Khan NS, Zuhra S, Shah Q. Entropy generation in two phase model for simulating flow and heat transfer of carbon nanotubes between rotating stretchable disks with cubic auto-catalysis chemical reaction. Q. Appl Nanosci. 2019; https://doi.org/10.1007/s13204-019-01017-1.
  47. 47. Sharma PR, Singh G. Effects of variable thermal conductivity and heat source/sink on MHD flow near a stagnation point on a linearly stretching sheet. J. Appl. Fluid Mech. 2009; 2: 13–21.
  48. 48. Pal D, Chatterjee S. Soret and Dufour effects on MHD convective heat and mass transfer of a power-law fluid over an inclined plate with variable thermal conductivity in a porous medium. Appl. Math. Comp. 2013; 219: 7556–7574.
  49. 49. Lund LA, Omar Z, Khan I. Analysis of dual solution for MHD flow of Williamson fluid with slippage. Heliyon. 2019; 5(3): e01345. pmid:30949601
  50. 50. Khan NS, Gul T, Islam S, Khan A, Shah Z. Brownian motion and Thermophoresis effects on MHD mixed convective thin film second-grade nanofluid flow with hall effect and heat transfer past a stretching sheet. Journal of Nanofluids. 2017; 6(5): 812–829.
  51. 51. Khan NS, Gul T, Islam S, Khan I, Alqahtani AM, Alshomrani AS. Magnetohydrodynamic nanoliquid thin film sprayed on a stretching cylinder with heat transfer. Appl. Sci. 2017; 7: 271.
  52. 52. Zuhra S, Khan NS, Islam S. Magnetohydrodynamic second-grade nanofluid flow containing nanoparticles and gyrotactic microorganisms. Computational and Applied Mathematics. 2018; 37(5): 6332–6358.
  53. 53. Khan NS, Gul T, Kumam P, Shah Z, Islam S, Khan W, et al. Influence of inclined magnetic field on Carreau nanoliquid thin film flow and heat transfer with Graphene nanoparticles. Energies. 2019; 12(8): 1459.
  54. 54. Dogonchi AS, Selimefendigil F, Ganji DD. Magneto-hydrodynamic natural convection of CuO-water nanofluid in complex shaped enclosure considering various nanoparticle shapes. International Journal of Numerical Methods for Heat Fluid Flow. 2019; 29(5): 1663–1679.
  55. 55. Dogonchi AS, Sheremet MA, Pop I, Ganji DD. MHD natural convection of CuH2O nanofluid in a horizontal semi-cylinder with a local triangular heater. International Journal of Numerical Methods for Heat Fluid Flow. 2019; 28(12): 2979–2996.
  56. 56. Rasool G, Zhang T. Characteristics of chemical reaction and convective boundary conditions in Powell-Eyring nanofluid flow along a radiative Riga plate. Heliyon. 2019; 5: e01479. pmid:30997431