Muhammad Ramzan^1,2, Nomana Abid¹, Dianchen Lu³, Iskander Tlili^,4,51Department of Computer Science, Bahria University, 44000, Islamabad, Pakistan
²Department of Mechanical Engineering, Sejong University, Seoul 143-747, Republic of Korea
³Department of Mathematics, Faculty of Science, Jiangsu University, Zhenjiang 212013, China
⁴Department for Management of Science and Technology Development, Ton Duc Thang University, Ho Chi Minh City, Vietnam
⁵Faculty of Applied Sciences, Ton Duc Thang University, Ho Chi Minh City, Vietnam

Received:2020-01-11Revised:2020-04-11Accepted:2020-04-11Online:2020-07-20


Abstract
This study aims to investigate the time-dependent squeezing of nanofluid flow, comprising carbon nanotubes of dual nature, e.g. single-walled carbon nanotubes, and multi-walled carbon nanotubes, between two parallel disks. Numerical simulations of the proposed novel model are conducted, accompanied by Cattaneo-Christov heat flux in a Darcy-Forchheimer permeable media. Additional impacts of homogeneous–heterogeneous reactions are also noted, including melting heat. A relevant transformation procedure is implemented for the transition of partial differential equations to the ordinary variety. A computer software-based MATLAB function, bvp4c, is implemented to handle the envisioned mathematical model. Sketches portraying impacts on radial velocity, temperature, and concentration of the included parameters are given, and deliberated upon. Skin friction coefficient and local Nusselt number are evaluated via graphical illustrations. It is observed that the local inertia coefficient has an opposite impact on radial velocity and temperature field. It is further perceived that melting and radiation parameters demonstrate a retarding effect on temperature profile.
Keywords： melting heat transfer;Darcy-Forchheimer porous media;Cattaneo-Christov heat flux;carbon nanotubes;squeezing flow;homogeneous–heterogeneous reactions


PDF (1870KB)MetadataMetricsRelated articlesExportEndNote|Ris|BibtexFavorite
Cite this article
Muhammad Ramzan, Nomana Abid, Dianchen Lu, Iskander Tlili. Impact of melting heat transfer in the time-dependent squeezing nanofluid flow containing carbon nanotubes in a Darcy-Forchheimer porous media with Cattaneo-Christov heat flux. Communications in Theoretical Physics, 2020, 72(8): 085801- doi:10.1088/1572-9494/ab8a2c

Acronyms $\quad$Description
MWCNTs $\quad$Multi-walled carbon nanotubes
MHD $\quad$Magnetohydrodynamics
SWCNTs $\quad$Single-walled carbon nanotubes
Symbols $\quad$Description
$b\left(t\right)$ $\quad$Strength of magnetic field
$c$ $\quad$Dimensional constant
${C}_{p}$ $\quad$Capacity of specific heat
${c}_{s}$ $\quad$Heat capacity of solid surface
${\left({C}_{p}\right)}_{{\rm{nf}}}$ $\quad$Specific heat capacity of nanofluid
${D}_{m},{D}_{n}$ $\quad$Diffusion coefficients
${\left({C}_{p}\right)}_{{\rm{f}}}$ $\quad$Specific heat capacity of fluid
${F}^{* }$ $\quad$Forchheimer parameter
$Ha$ $\quad$Hartmann number
${k}_{c}$ $\quad$Thermal conductivity of carbon nanotubes
${K}^{* }$ $\quad$Permeability of spongy media
${k}_{1}$ $\quad$Homogeneous reaction parameter
${k}_{2}$ $\quad$Heterogeneous reaction parameter
$L$ $\quad$Latent heat
$m,n,{K}_{i},{K}_{j}$ $\quad$Concentrations of chemical species
$F$ $\quad$Local inertia coefficient
$k$ $\quad$Thermal conductivity
$m* ,n* $ $\quad$Chemical species
$Me$ $\quad$Melting heat coefficient
$P$ $\quad$Pressure
${q}_{{\rm{rd}}}$ $\quad$Radiative heat flux
$Pr$ $\quad$Prandtl number
${R}_{d}$ $\quad$Thermal radiation coefficient
$R{e}_{L}$ $\quad$Local squeezed Reynolds number
$Re$ $\quad$Reynolds parameter
$Sc$ $\quad$Schmidt parameter
$Sq$ $\quad$Squeeze parameter
$T$ $\quad$Temperature
${T}_{0}$ $\quad$Temperature of solid
${T}_{H}$ $\quad$Temperature at upper disk
${T}_{M}$ $\quad$Surface temperature
$\left(U* ,W* \right)$ $\quad$Components of velocities
$K$ $\quad$Kelvin
${\rm{W}}$ $\quad$Watt
Greek symbols $\quad$Description
$\chi $ $\quad$Similarity transformation variable
$\theta \left(\chi \right)$ $\quad$Dimensionless temperature
$f^{\prime} \left(\chi \right)$ $\quad$Radial velocity (dimensionless velocity)
${k}_{c}$ $\quad$Carbon nanotubes thermal conductivity
${\mu }_{{\rm{f}}}$ $\quad$fluid dynamic viscosity
${k}_{{\rm{nf}}}$ $\quad$Thermal conductivity of nanofluid
$\rho $ $\quad$Density
${\rho }_{{\rm{nf}}}$ $\quad$Nanofluid density
${\sigma }_{{\rm{nf}}}$ $\quad$Nanofluid electrical conductivity
${\nu }_{{\rm{f}}}$ $\quad$Kinematic viscosity of the fluid
$\sigma * $ $\quad$Stefan Boltzmann constant
$\gamma $ $\quad$Thermal relaxation coefficient
${\theta }_{w}$ $\quad$Temperature ratio parameter
${\rho }_{c}$ $\quad$Density of carbon nanotubes
$f\left(\chi \right)$ $\quad$Axial velocity (dimensionless velocity)
$H\left(\chi \right),G\left(\chi \right)$ $\quad$Dimensionless concentrations
${k}_{{\rm{f}}}$ $\quad$Base fluid (water) thermal conductivity
$\mu $ $\quad$dynamic viscosity
${\mu }_{{\rm{nf}}}$ $\quad$Nanofluid dynamic viscosity
${\rho }_{{\rm{f}}}$ $\quad$Density of fluid
${\nu }_{{\rm{nf}}}$ $\quad$Nanofluid kinematic viscosity
$\phi $ $\quad$Volume fraction of nanoparticles
$\nu $ $\quad$Kinematic viscosity
$k* $ $\quad$Mean absorption coefficient
$\lambda $ $\quad$Porosity parameter
Subscripts $\quad$Description
0 $\quad$For solid
1 $\quad$For homogeneous
2 $\quad$For heterogeneous
$c$ $\quad$Carbon nanotubes
${\rm{f}}$ $\quad$Fluid
$H$ $\quad$For upper disk
$L$ $\quad$For local
$M$ $\quad$For lower disk
${\rm{nf}}$ $\quad$Nanofluid
$p$ $\quad$Constant pressure
${\rm{rd}}$ $\quad$Radiative
$s$ $\quad$Solid surface

Introduction

A carbon nanotube is a large, stretched, thin, and tube-shaped molecule of pure carbon of around 1 to 3 nanometers (1 to 3 billionth of a meter) in breadth (diameter), and 100 to 1000 s of nanometers in length. Iijima [1] introduced the theory of carbon nanotubes at the end of the 20th century, when he discovered some potential applications of CNTs for solar cells, radar-absorbing coating, gas storage, composites, semiconductor devices, ultra-capacitors, etc [2]. CNTs are classified as SWCNTs and MWCNTs. A SWCNT has a regular straw shape with only one layer. A MWCNT is a set of nested SWCNTs of increasing diameters. Din and Khan [3] studied the squeezing flow of Casson fluid with non-linear thermal radiation between parallel disks. Haq et al [4] examined MHD nanofluid squeezed flow based on water with CNTs between 2-parallel disks, and concluded that temperature and velocity profiles increase with high nanoparticle volume fraction. Melting heat in the radiative flow of CNTs with homogeneous–heterogeneous reactions was scrutinized by Hayat et al [5]. They found that the Nusselt number increases for large values of nanoparticle volume fraction.

Recently, squeezed flow between two parallel disks has garnered a great deal of attention, with its vast number of potential applications in technological and industrial systems. Many devices such as stirring pistons in engines, hydraulic brakes, and chocolate filler are based on the flow principle between squeezing regions. Stefan [6] proposed the idea of squeezing flow in 1874. Since then, many researchers have explored the problems associated with squeezing flow. The theoretical investigation regarding squeezing flow between parallel disks is presented by Leider and Bird [7]. Qayyum et al [8] discussed time-dependent squeezing Jeffery fluid flow between two parallel disks. Hayat et al [9] discussed squeezed nanofluid flow-based CNTs, and the impact of thermal radiations on Darcy-Forchheimer spongy media. They posited that an augmentation in nanoparticle volume fraction causes a reduction in both the velocity and the temperature of the fluid. Hashmi et al [10] investigated the analytical simulations for squeezing nanofluid flow amidst parallel disks. It should be noted that these studies are discuss the use of two parallel disks, but do not refer to Cattaneo-Christov heat flux.

The Fourier law of heat conduction has been a criterion benchmark in many practical industries for estimating the behavior of heat transmission. Nevertheless, because of the parabolic-heat equation due to an initial disorder, this system suffers a great deal. Cattaneo [11] tackled this drawback of the Fourier model via the addition of thermal time relaxation. This modification has created a hyperbolic heat equation for the temperature field. Also, within finite speeds, heat transmission is permitted to circulate through thermal waves. Tibullo and Zampoli [12] have worked on innumerable practicable applications i.e. to nanofluid flow, applying the Cattaneo-Christov heat conduction model. Christov [13] posited a modification of the Maxwell-Cattaneo model, which is known as the Cattaneo-Christov thermal flux model. Radiative nanofluid flow with the Cattaneo-Christov heat flux model between parallel disks is studied by Dogonchi et al [14]. Lu et al [15] discussed the mathematical model of unsteady fluid flow containing SWCNTs and MWCNTs under conditions of Cattaneo-Christov heat flux, and homogeneous–heterogeneous reactions between two parallel disks, and determined that temperature rises with an increasing thermal relaxation parameter. Zubair et al [16] discussed the 3D Darcy-Forchheimer squeezing nanofluid flow with Cattaneo-Christov heat flux, using four distinct types of nanoparticles, via the analysis of entropy generation.

Henry Darcy [17] determined the fluid flow over a permeable surface, based on the outcomes of water flow experiments over cribs of sand, and hydro-geology. He defined his idea of fluid flow over a spongy media in 1856. Due to its limitations of small velocity with weaker permeability, Philipps Forchheimer [18] modified the momentum equation by velocity square ${v}^{2}$ within Darcian velocity. This became known as the Forchheimer term, as designated by Muskat [19]. Nasir et al [20] scrutinized the radiative 3D Darcy-Forchheimer flow of carbon nanotubes past a stretchable rotating disk. They posited that fluid velocity falls with an upsurge in inertia and porosity parameters. Jha and Kaurangini [21] presented analytic solutions for Darcy-Forchheimer-based spongy media relations. Kaladhar [22] discussed mixed convection flow with dual stratification effects within a Darcy-Forchheimer medium, finding that velocity is reduced for higher estimates of Darcy-Forchheimer number. A numerical solution for a second law analysis of ferrofluid within a spongy semi annulus is investigated by Sheikholeslami et al [23]. Alamri et al [24] proposed the model of a radiative plane Poiseuille flow of nanofluid, using slip conditions past a porous medium. Some recent investigations discussing porous media are referenced in [25–28], among many others.

Based on the above discussion, it can be noted that there is as yet no study in which the influences of Cattaneo-Christov heat flux with carbon nanotubes in a Darcy-Forchheimer porous media are examined between two parallel disks. Comparatively little research has been conducted in the area of carbon nanotubes to date. For this reason, this paper aims to examine melting heat transfer effects in carbon nanotubes- (SWCNTs and MWCNTs) based nanofluid unsteady flow, in a non-linear Darcy-Forchheimer permeable media, between two parallel disks, with Cattaneo-Christov heat flux and homogeneous–heterogeneous reactions. The impact of prominent parameters on surface drag force, and Nusselt number, are portrayed via graphic illustrations. The numerical solution of the present work is obtained by adopting the Finite difference method, this being the default in the bvp4c built-in function of the MATLAB scheme. The layout of this paper consists of: Section one – introduction. Section two is devoted to mathematical modeling, with all required equations, having employed the boundary layer theory to partial differential equations. Section three is a detailed elaboration of the numerical methods applied to the problem. Section four covers results and discussion. The paper concludes with Final remarks and summing up.

Mathematical modeling

Consider an incompressible, time-dependent 2D MHD nanofluid flow, containing CNTs within a Darcy-Forchheimer spongy media under conditions of non-linear thermal radiation and melting heat transfer, between two infinite parallel disks of length(1)$ \begin{eqnarray}Z=x(t)=h\sqrt{1-ct},\end{eqnarray}$with applied magnetic strength $b(t)=\tfrac{{b}_{0}}{\sqrt{1-ct}}$ which is normal to the disks (figure 1). Here, SWCNTs and MWCNTs, along with water (base fluid), are considered. Moreover, the upper disk $Z=x\left(t\right)$ moves up and down with a velocity $\tfrac{{\rm{d}}}{{\rm{dt}}}(Z)$ from the fixed and porous lower disk $Z=0.$ The induced magnetic field is neglected here, because we are using a small Reynolds number. We also consider the cylindrical coordinate system $(R,\alpha ,Z).$ The velocity component $V$ vanishes identically due to rotational flow symmetry $\left(\tfrac{\partial }{\partial \alpha }=0\right).$

Figure 1.

New window|Download| PPT slide
Figure 1.Fluid geometry.


We assume a model, designed by Chaudhary and Merkin [29] for homogeneous–heterogeneous chemical reactions defined as:(2)$ \begin{eqnarray}m* +2n* \to 3n* ,\,{\rm{rate}}\,{K}_{i}m{n}^{2},\end{eqnarray}$(3)$ \begin{eqnarray}m* \to n* ,\,{\rm{rate}}\,{K}_{j}m.\end{eqnarray}$

These reactions are presumed to be isothermal. The governing system with boundary layer equations are represented as:(4)$ \begin{eqnarray}\displaystyle \frac{\partial U* }{\partial R}+\displaystyle \frac{U* }{R}+\displaystyle \frac{\partial W* }{\partial Z}=0,\end{eqnarray}$(5)$ \begin{eqnarray}\begin{array}{l}\displaystyle \frac{\partial U* }{\partial t}+U\displaystyle \frac{\partial U* }{\partial R}+W\displaystyle \frac{\partial U* }{\partial Z}=-\displaystyle \frac{1}{{\rho }_{{\rm{nf}}}}\displaystyle \frac{\partial P}{\partial R}+\displaystyle \frac{{\mu }_{{\rm{nf}}}}{{\rho }_{{\rm{nf}}}}\left(\displaystyle \frac{{\partial }^{2}U* }{\partial {R}^{2}}\right.\\ \,+\,\left.\displaystyle \frac{1}{R}\displaystyle \frac{\partial U* }{\partial R}-\displaystyle \frac{U* }{{R}^{2}}+\displaystyle \frac{{\partial }^{2}U* }{\partial {Z}^{2}}\right)-\displaystyle \frac{{\sigma }_{{\rm{f}}}}{{\rho }_{{\rm{nf}}}}{b}^{2}\left(t\right)U* \\ \,-\,\displaystyle \frac{{\nu }_{{\rm{nf}}}U* }{K}-F* {U}^{* 2},\end{array}\end{eqnarray}$(6)$ \begin{eqnarray}\begin{array}{l}\displaystyle \frac{\partial W* }{\partial t}+U\displaystyle \frac{\partial W* }{\partial R}+W\displaystyle \frac{\partial W* }{\partial Z}=-\displaystyle \frac{1}{{\rho }_{{\rm{nf}}}}\displaystyle \frac{\partial P}{\partial Z}+\displaystyle \frac{{\mu }_{{\rm{nf}}}}{{\rho }_{{\rm{nf}}}}\left(\displaystyle \frac{{\partial }^{2}W* }{\partial {R}^{2}}\right.\\ \,+\,\left.\displaystyle \frac{1}{R}\displaystyle \frac{\partial W* }{\partial R}+\displaystyle \frac{{\partial }^{2}W* }{\partial {Z}^{2}}\right)-\displaystyle \frac{{\nu }_{{\rm{nf}}}W* }{K}-F* {W}^{* 2},\end{array}\end{eqnarray}$(7)$ \begin{eqnarray}\begin{array}{l}\displaystyle \frac{\partial T}{\partial t}+U* \displaystyle \frac{\partial T}{\partial R}+W* \displaystyle \frac{\partial T}{\partial Z}\\ \,+\,\varepsilon \left(\begin{array}{l}\displaystyle \frac{{\partial }^{2}T}{\partial {t}^{2}}+\displaystyle \frac{\partial U* }{\partial t}\displaystyle \frac{\partial T}{\partial R}+2U* \displaystyle \frac{{\partial }^{2}T}{\partial t\partial R}+2W* \displaystyle \frac{{\partial }^{2}T}{\partial t\partial Z}\\ +\displaystyle \frac{\partial W* }{\partial t}\displaystyle \frac{\partial T}{\partial Z}+U* \displaystyle \frac{\partial U* }{\partial R}\displaystyle \frac{\partial T}{\partial R}+W* \displaystyle \frac{\partial W* }{\partial Z}\displaystyle \frac{\partial T}{\partial Z}+U* \displaystyle \frac{\partial W* }{\partial R}\displaystyle \frac{\partial T}{\partial Z}\\ +W* \displaystyle \frac{\partial U* }{\partial Z}\displaystyle \frac{\partial T}{\partial R}+2U* W* \displaystyle \frac{{\partial }^{2}T}{\partial R\partial Z}+{U}^{* 2}\displaystyle \frac{{\partial }^{2}T}{\partial {R}^{2}}+{W}^{* 2}\displaystyle \frac{{\partial }^{2}T}{\partial {Z}^{2}}\end{array}\right)\\ \,=\,-\displaystyle \frac{{k}_{{\rm{nf}}}}{{\left(\rho {C}_{p}\right)}_{{\rm{nf}}}}\left(\displaystyle \frac{{\partial }^{2}T}{\partial {R}^{2}}+\displaystyle \frac{1}{R}\displaystyle \frac{\partial T}{\partial R}+\displaystyle \frac{{\partial }^{2}T}{\partial {Z}^{2}}\right)-\displaystyle \frac{1}{{\left(\rho {C}_{p}\right)}_{{\rm{nf}}}}\displaystyle \frac{\partial {q}_{{\rm{rd}}}}{\partial Z},\end{array}\end{eqnarray}$(8)$ \begin{eqnarray}\begin{array}{l}\displaystyle \frac{\partial m}{\partial t}+U* \displaystyle \frac{\partial m}{\partial R}+W* \displaystyle \frac{\partial m}{\partial Z}={D}_{m}\left(\displaystyle \frac{{\partial }^{2}m}{\partial {R}^{2}}+\displaystyle \frac{1}{R}\displaystyle \frac{\partial m}{\partial R}+\displaystyle \frac{{\partial }^{2}m}{\partial {Z}^{2}}\right)\\ \,-\,{K}_{i}m{n}^{2},\end{array}\end{eqnarray}$(9)$ \begin{eqnarray}\begin{array}{l}\displaystyle \frac{\partial n}{\partial t}+U* \displaystyle \frac{\partial n}{\partial R}+W* \displaystyle \frac{\partial n}{\partial Z}={D}_{n}\left(\displaystyle \frac{{\partial }^{2}n}{\partial {R}^{2}}+\displaystyle \frac{1}{R}\displaystyle \frac{\partial n}{\partial R}+\displaystyle \frac{{\partial }^{2}n}{\partial {Z}^{2}}\right)\\ \,+\,{K}_{i}m{n}^{2}.\end{array}\end{eqnarray}$With boundary conditions(10)$ \begin{eqnarray}\begin{array}{l}U* =0,\,W* =\displaystyle \frac{d}{dt}(Z),\,T={T}_{H},\,m\to {m}_{0},\,n\to 0,\\ {\rm{at}}\,Z=x\left(t\right),\end{array}\end{eqnarray}$(11)$ \begin{eqnarray}\begin{array}{l}U* =0,\,{k}_{nf}{\left(\displaystyle \frac{\partial T}{\partial Z}\right)}_{Z=0}={\rho }_{nf}\left({c}_{s}\left({T}_{M}-{T}_{0}\right)+L\right)W* \left(R,0\right),\\ T={T}_{M},\,{D}_{m}\displaystyle \frac{\partial m}{\partial Z}=m{K}_{j},\,{D}_{n}\displaystyle \frac{\partial n}{\partial Z}=-m{K}_{j},\,{\rm{at}}\,Z=0.\end{array}\end{eqnarray}$

Mathematically, thermophysical properties are shown as:(12)$ \begin{eqnarray}\begin{array}{l}{\mu }_{{\rm{nf}}}=\displaystyle \frac{{\mu }_{{\rm{f}}}}{{\left(1-\phi \right)}^{2.5}},\,{\nu }_{nf}=\displaystyle \frac{{\mu }_{nf}}{{\rho }_{nf}},\,{\alpha }_{{\rm{nf}}}=\displaystyle \frac{{k}_{{\rm{nf}}}}{{\left(\rho {C}_{p}\right)}_{{\rm{nf}}}},\\ {\nu }_{{\rm{f}}}=\displaystyle \frac{{\mu }_{{\rm{f}}}}{{\rho }_{{\rm{f}}}},\,\displaystyle \frac{{k}_{{\rm{nf}}}}{{k}_{{\rm{f}}}}=\displaystyle \frac{\left(1-\phi \right)+2\phi \displaystyle \frac{{k}_{C}}{{k}_{C}-{k}_{{\rm{f}}}}\,\mathrm{ln}\,\displaystyle \frac{{k}_{C}+{k}_{{\rm{f}}}}{2{k}_{f}}}{\left(1-\phi \right)+2\phi \displaystyle \frac{{k}_{C}}{{k}_{C}-{k}_{{\rm{f}}}}\,\mathrm{ln}\,\displaystyle \frac{{k}_{C}+{k}_{{\rm{f}}}}{2{k}_{{\rm{f}}}}},\\ {\rho }_{{\rm{nf}}}=\left(1-\phi \right){\rho }_{{\rm{f}}}+\phi {\rho }_{C},\,{\left(\rho {C}_{p}\right)}_{{\rm{nf}}}=\left(1-\phi \right){\left(\rho {C}_{p}\right)}_{{\rm{f}}}\\ \,+\,\phi {\left(\rho {C}_{p}\right)}_{C}.\end{array}\end{eqnarray}$

The thermophysical features of water and CNTs are appended in table 1.


Table 1.
Table 1.Thermophysical characteristics of base fluid, SWCNTs and MWCNTs [30].

Thermophysical traits	${C}_{p}\left(\displaystyle \frac{{\rm{J}}}{{\rm{kgK}}}\right)$	$\rho \left(\displaystyle \frac{{\rm{kg}}}{{{\rm{m}}}^{{\rm{3}}}}\right)$	$k\left(\displaystyle \frac{{\rm{W}}}{{\rm{mK}}}\right)$
Water (base fluid)	4179.00	997.100	0.613 00
Nanoparticles (SWCNTs)	425	2600	6600
Nanoparticles (MWCNTs)	796	1600	3000

New window|CSV

From equation (7), by utilizing the Roseland thermal radiation approximation [31], we obtain the value of ${q}_{{\rm{rd}}}$ as:(13)$ \begin{eqnarray}{q}_{{\rm{rd}}}=-\displaystyle \frac{4\sigma * }{3k* }\displaystyle \frac{\partial {T}^{4}}{\partial Z}=-\displaystyle \frac{16\sigma * }{3k* }{T}^{3}\displaystyle \frac{\partial T}{\partial Z},\end{eqnarray}$with(14)$ \begin{eqnarray}T={T}_{M}\left[1+\left({\theta }_{w}-1\right)\theta \right].\end{eqnarray}$

Similarity transformation

Similarity transformations are defined as:(15)$ \begin{eqnarray}\begin{array}{l}U* \,=\,\displaystyle \frac{cR}{2\left(1-ct\right)}f^{\prime} \left(\chi \right),\,W* \,=\,-\displaystyle \frac{ch}{{\left(1-ct\right)}^{\displaystyle \frac{1}{2}}}f\left(\chi \right),\\ m={m}_{0}H\left(\chi \right),\,n={m}_{0}G\left(\chi \right),\,b\left(t\right)=\displaystyle \frac{{b}_{0}}{{\left(1-ct\right)}^{\displaystyle \frac{1}{2}}},\\ \chi =\displaystyle \frac{Z}{h{\left(1-ct\right)}^{\displaystyle \frac{1}{2}}},\,\theta =\displaystyle \frac{T-{T}_{M}}{{T}_{H}-{T}_{M}}.\end{array}\end{eqnarray}$

By means of the above transformation, equation (4) is satisfied, and equations (5)–(9) are transformed into:(16)$ \begin{eqnarray}\begin{array}{l}\displaystyle \frac{1}{{\left(1-\phi \right)}^{2.5}}f^{\prime\prime\prime\prime} \left(\chi \right)-\left[\left(1-\phi \right)+\phi \displaystyle \frac{{\rho }_{c}}{{\rho }_{f}}\right]\\ \,\times \,\left[\begin{array}{l}3Sqf^{\prime\prime} \left(\chi \right)+Sq\chi f^{\prime\prime\prime} \left(\chi \right)\\ -2Sqf\left(\chi \right)f^{\prime\prime\prime} \left(\chi \right)+\mathrm{Re}Ff^{\prime} \left(\chi \right)f^{\prime\prime} \left(\chi \right)\end{array}\right]\\ \,-\,\left[\displaystyle \frac{\lambda \mathrm{Re}}{{\left(1-\phi \right)}^{2.5}}+Ha\right]f^{\prime\prime} \left(\chi \right)=0,\end{array}\end{eqnarray}$(17)$ \begin{eqnarray}\begin{array}{l}\displaystyle \frac{{k}_{{\rm{nf}}}}{{k}_{{\rm{f}}}}\left\{\left[1+{R}_{d}{\left(1+\left({\theta }_{w}-1\right)\theta \right)}^{3}\right]\theta ^{\prime} \right\}^{\prime} \\ \,+\,SqPr\left[\left(1-\phi \right)+\phi \displaystyle \frac{{\left(\rho {C}_{p}\right)}_{c}}{{\left(\rho {C}_{p}\right)}_{{\rm{f}}}}\right]\\ \left\{\left(2f\left(\chi \right)-\chi \right)\theta ^{\prime} \left(\chi \right)-\gamma \left[\begin{array}{l}\left(\begin{array}{l}{\chi }^{2}-4\chi f\left(\chi \right)\\ +4{f}^{2}\left(\chi \right)\end{array}\right)\theta ^{\prime\prime} \left(\chi \right)\\ +\left(\begin{array}{l}4f\left(\chi \right)f^{\prime} \left(\chi \right)+3\chi \\ -2\chi f^{\prime} \left(\chi \right)-6f\left(\chi \right)\end{array}\right)\theta ^{\prime} \left(\chi \right)\end{array}\right]\right\}\\ \,=\,0,\end{array}\end{eqnarray}$(18)$ \begin{eqnarray}\displaystyle \frac{1}{Sc}H^{\prime\prime} \left(\chi \right)+f\left(\chi \right)H^{\prime} \left(\chi \right)-\displaystyle \frac{\chi }{2}H^{\prime} \left(\chi \right)-{k}_{1}H\left(\chi \right){G}^{2}\left(\chi \right)=0,\end{eqnarray}$(19)$ \begin{eqnarray}\displaystyle \frac{\delta }{Sc}G^{\prime\prime} \left(\chi \right)+f\left(\chi \right)G^{\prime} \left(\chi \right)-\displaystyle \frac{\chi }{2}G^{\prime} \left(\chi \right)+{k}_{1}H\left(\chi \right){G}^{2}\left(\chi \right)=0,\end{eqnarray}$with(20)$ \begin{eqnarray}\delta =\displaystyle \frac{{D}_{n}}{{D}_{m}}.\end{eqnarray}$

Based on equations (2) and (3), chemical species $m* $ and $n* $ cannot be analogous, but both can be identical in magnitude, provided $\delta =1.$ Thus, from equation (20), presuming that ${D}_{n}$ and ${D}_{m}$ are identical (i.e., $\delta =1$), we obtain(21)$ \begin{eqnarray}H\left(\chi \right)+G\left(\chi \right)=1,\end{eqnarray}$using the above property, equations (18) and (19) become(22)$ \begin{eqnarray}\begin{array}{l}\displaystyle \frac{1}{Sc}H^{\prime\prime} \left(\chi \right)+f\left(\chi \right)H^{\prime} \left(\chi \right)-\displaystyle \frac{\chi }{2}H^{\prime} \left(\chi \right)\\ \,-\,{k}_{1}H\left(\chi \right){\left[1-H\left(\chi \right)\right]}^{2}=0,\end{array}\end{eqnarray}$

and the boundary Equations (10) and (11) become(23)$ \begin{eqnarray}\begin{array}{l}f^{\prime} \left(0\right)=0,\,\theta \left(0\right)=0,\,\displaystyle \frac{{k}_{{\rm{nf}}}}{{k}_{{\rm{f}}}}{M}_{e}\theta ^{\prime} \left(0\right)+PrRe\left[\Space{0ex}{3.20ex}{0ex}\left(1-\phi \right)\right.\\ \,+\left.\phi \displaystyle \frac{{\rho }_{c}}{{\rho }_{{\rm{f}}}}\right]f\left(0\right)=0,\,H^{\prime} \left(0\right)={k}_{2}H\left(0\right),\,{\rm{at}}\,\chi =0,\end{array}\end{eqnarray}$(24)$ \begin{eqnarray}f(1)=\displaystyle \frac{1}{2},\,f^{\prime} (1)=0,\,\theta (1)=1,\,H\left(1\right)\to 1,\,{\rm{at}}\,\chi =1.\end{eqnarray}$

Based on equation (23), ${M}_{e}$ is the melting heat coefficient as:(25)$ \begin{eqnarray}{M}_{e}=\left[\displaystyle \frac{{c}_{f}\left({T}_{H}-{T}_{M}\right)}{L+{c}_{s}\left({T}_{M}-{T}_{0}\right)}\right],\end{eqnarray}$which is the amalgamation of two numbers, ${c}_{f}\left({T}_{H}-{T}_{M}\right)/L$, and ${c}_{s}\left({T}_{M}-{T}_{0}\right)/L$, known as Stefan numbers, for solid and liquid states. In the above equations, non-dimensional coefficients are defined as:(26)$ \begin{eqnarray}\begin{array}{l}Sq=\displaystyle \frac{c{h}^{2}}{2{\nu }_{f}},\,Ha={\left(\displaystyle \frac{{\sigma }_{f}{{b}_{0}}^{2}{h}^{2}}{{\mu }_{f}}\right)}^{\displaystyle \frac{1}{2}},\,Pr=\displaystyle \frac{{\mu }_{f}{\left(\rho {C}_{p}\right)}_{c}}{{\rho }_{f}{k}_{f}},\\ Re=\displaystyle \frac{c{h}^{2}}{{\nu }_{{\rm{f}}}},\,F=F* R=\displaystyle \frac{{C}_{b}R}{{\left({K}^{* }\right)}^{\displaystyle \frac{1}{2}}},\,\lambda =\displaystyle \frac{{\nu }_{{\rm{f}}}}{cK* },\\ {R}_{d}=\displaystyle \frac{16\sigma * {{T}_{M}}^{3}}{3k* {k}_{{\rm{f}}}},\,{k}_{1}=\displaystyle \frac{{K}_{i}{{A}_{0}}^{2}\left(1-ct\right)}{c},\\ \gamma =\displaystyle \frac{\varepsilon c}{2\left(1-ct\right)},Sc=\displaystyle \frac{c{h}^{2}}{{D}_{m}},{k}_{2}=\displaystyle \frac{{K}_{j}h{\left(1-ct\right)}^{\displaystyle \frac{1}{2}}}{{D}_{m}}.\end{array}\end{eqnarray}$

Surface drag force and rate of heat flux are classified by:(27)$ \begin{eqnarray}\begin{array}{l}{C}_{fr}=\displaystyle \frac{{\mu }_{{\rm{nf}}}{\left(\displaystyle \frac{\partial U* }{\partial Z}+\displaystyle \frac{\partial W* }{\partial R}\right)}_{Z=x\left(t\right)}}{{\rho }_{{\rm{f}}}{\left(-\displaystyle \frac{ch}{2\sqrt{1-ct}}\right)}^{2}},\\ Nu=\displaystyle \frac{{k}_{{\rm{nf}}}h}{{k}_{{\rm{f}}}\left({T}_{H}-{T}_{M}\right)}\left[{\left({q}_{r}\right)}_{w}-{\left(\displaystyle \frac{\partial T}{\partial Z}\right)}_{Z=x\left(t\right)}\right],\end{array}\end{eqnarray}$using equation (15), we get(28)$ \begin{eqnarray}\begin{array}{l}\displaystyle \frac{{h}^{2}}{{R}^{2}}R{e}_{L}{C}_{fr}=\left[\displaystyle \frac{1}{{\left(1-\phi \right)}^{2.5}\left[\left(1-\phi \right)+\phi \displaystyle \frac{{\rho }_{c}}{{\rho }_{{\rm{f}}}}\right]}\right]f^{\prime\prime} \left(1\right),\\ \sqrt{1-ct}Nu=-\left[\displaystyle \frac{{k}_{{\rm{nf}}}}{{k}_{{\rm{f}}}}\right.\\ \,+\,\left.{R}_{d}{\left(1+\left({\theta }_{W}-1\right)\theta \left(1\right)\right)}^{3}\right]\theta ^{\prime} \left(1\right),\end{array}\end{eqnarray}$where(29)$ \begin{eqnarray}R{e}_{L}=\displaystyle \frac{{\rho }_{f}Rch\sqrt{1-ct}}{2{\mu }_{{\rm{f}}}}.\end{eqnarray}$

Solution methodology

For non-linear systems of ODEs (16), (17), and (22), with boundary conditions (23) and (24), we employ the finite difference default method of the bvp4c built-in function of the MATLAB scheme, which is fourth order accurate, and a grid size of 0.01 is taken with the tolerance ${10}^{-7}.$ Using the following numerical code, we obtain first order ODEs as:(30)$ \begin{eqnarray}\begin{array}{l}{y}_{\left(1\right)}=f\left(\chi \right),\,{y}_{\left(2\right)}=f^{\prime} \left(\chi \right),\,{y}_{\left(3\right)}=f^{\prime\prime} \left(\chi \right),\\ {y}_{\left(4\right)}=f^{\prime\prime} \left(\chi \right),\\ {f}^{////}\left(\chi \right)=\left[{y^{\prime} }_{\left(4\right)}\right]=y{y}_{\left(1\right)}={\left(1-\phi \right)}^{2.5}\left\{\Space{0ex}{3.40ex}{0ex}\left(1-\phi \right)\right.\\ \,+\left.\phi \displaystyle \frac{{\rho }_{c}}{{\rho }_{f}}\right\}\left\{\begin{array}{l}3Sq{y}_{\left(3\right)}+Sq\chi {y}_{\left(4\right)}\\ -2Sq{y}_{\left(1\right)}{y}_{\left(4\right)}+ReF{y}_{\left(2\right)}{y}_{\left(3\right)}\end{array}\right\}\\ \,+\left\{\lambda Re+Ha{\left(1-\phi \right)}^{2.5}\right\}{y}_{\left(3\right)},\end{array}\end{eqnarray}$(31)$ \begin{eqnarray}\begin{array}{l}{y}_{\left(5\right)}=\theta \left(\eta \right),\,{y}_{\left(6\right)}=\theta ^{\prime} \left(\chi \right),\\ \theta ^{\prime\prime} \left(\chi \right)=\left[{y^{\prime} }_{\left(6\right)}\right]=y{y}_{\left(2\right)}\\ \,=\,\left(\displaystyle \frac{\begin{array}{l}-SqPr\left\{\left(1-\phi \right)+\phi \displaystyle \frac{{\left(\rho {C}_{p}\right)}_{c}}{{\left(\rho {C}_{p}\right)}_{f}}\right\}\left[\begin{array}{l}\left(2{y}_{(1)}-\chi \right){y}_{\left(6\right)}-\\ \gamma \left\{\begin{array}{l}4{y}_{\left(1\right)}{y}_{\left(2\right)}+3\chi \\ -2\chi {y}_{\left(2\right)}-6{y}_{\left(1\right)}\end{array}\right\}{y}_{\left(6\right)}\end{array}\right]\\ -\displaystyle \frac{{k}_{nf}}{{k}_{f}}3{R}_{d}{\left\{1+\left({\theta }_{w}-1\right){y}_{\left(5\right)}\right\}}^{2}{\left({y}_{\left(6\right)}\right)}^{2}\end{array}}{\begin{array}{l}SqPr\gamma \left\{\left(1-\phi \right)+\phi \displaystyle \frac{{\left(\rho {C}_{p}\right)}_{c}}{{\left(\rho {C}_{p}\right)}_{f}}\right\}\left(\begin{array}{l}{\chi }^{2}-4\chi {y}_{\left(1\right)}\\ +4{({y}_{(1)})}^{2}\end{array}\right)\\ +\displaystyle \frac{{k}_{nf}}{{k}_{f}}\left[1+{R}_{d}{\left\{1+\left({\theta }_{w}-1\right){y}_{\left(5\right)}\right\}}^{3}\right]\end{array}}\right);\\ {y}_{\left(7\right)}=H\left(\chi \right),\,{y}_{\left(8\right)}=H^{\prime} \left(\chi \right),\end{array}\end{eqnarray}$(32)$ \begin{eqnarray}\begin{array}{l}H^{\prime\prime} \left(\chi \right)=\left[{y^{\prime} }_{\left(8\right)}\right]=y{y}_{\left(3\right)}=Sc\left[\displaystyle \frac{\chi }{2}{y}_{(8)}-{y}_{\left(1\right)}{y}_{\left(8\right)}\right.\\ \,+\left.{k}_{1}{y}_{\left(7\right)}{\left[1-{y}_{\left(7\right)}\right]}^{2}\Space{0ex}{3.0ex}{0ex}\right].\end{array}\end{eqnarray}$With the boundary conditions(33)$ \begin{eqnarray}\begin{array}{l}{y}_{\left(2\right)}(0)=0,{y}_{\left(5\right)}(0)=0,\displaystyle \frac{{k}_{nf}}{{k}_{f}}{M}_{e}{y}_{\left(6\right)}(0)+\left\{\Space{0ex}{3.25ex}{0ex}\left(1-\phi \right)\right.\\ \,+\,\left.\phi \displaystyle \frac{{\rho }_{c}}{{\rho }_{f}}\right\}PrRe{y}_{(1)}(0)=0,{y}_{\left(8\right)}(0)={k}_{2}{y}_{\left(7\right)}(0),\end{array}\end{eqnarray}$(34)$ \begin{eqnarray}{y}_{1}(\infty )-\displaystyle \frac{1}{2};{y}_{2}(\infty );{y}_{5}(\infty )-1;{y}_{7}(\infty )-1.\end{eqnarray}$

Results and discussion

This segment examines the outcomes of dimensionless velocity $f^{\prime} (\chi ),$ temperature $\theta (\chi )$, and concentration fields $H(\chi )$ for numerous arising parameters, so as to reflect the behavior of fluid flow, and heat and mass transport. The subsequent discussion is presented for both SWCNTs and MWCNTs. Parameters used in this analysis are $Pr=6.2,Re=0.4,\lambda =0.2,Ha=0.5,{M}_{e}\,=\,0.01,\gamma \,=\,0.5,Sq=Sc=1.0,{k}_{1}={k}_{2}\,=0.7,\phi =F=0.1,{R}_{d}=0.9,{\theta }_{w}=1.1.$ Figure 2(a) illustrates how positive and negative values of squeezed number ${S}_{q}$ affect velocity $f^{\prime} (\chi ).$ The graph illustrates how, for both single and multi-walled CNTs, velocity profile $f^{\prime} (\chi )$ increases with the contraction of disks, i.e., negative values, whereas for positive values where upper and lower plates are driven further apart, the opposite behavior can be observed. In the case of contraction $\left({S}_{q}=-1,-2,-\mathrm{3...}\right),$ the fluid is exposed to a squeezed force, which causes it to move with increased velocity. Hence, velocity is augmented. Nevertheless, for ${S}_{q}=1,2,\mathrm{3...}$ when both disks move away from each other, a gap is produced between the disks. The fluid moves in a reverse direction to fill this gap; thus, velocity reduces. In figure 2(b) effect of the ${S}_{q}$ squeezing parameter is exhibited versus $\theta (\chi ).$ In the case of the contraction of disks $\left({S}_{q}=-1,-2,-\mathrm{3...}\right),$ the temperature profile $\theta (\chi )$ establishes diminishing behavior, whereas an opposite trend is seen when the disks are driven apart from each other i.e. ${S}_{q}=1,2,\mathrm{3...}.$ This is because movement of fluid increases when disks are driving away from each other, thus the temperature increases. An inverse impact of the melting heat transfer parameter ${M}_{e}$ can be observed in figures 3(a) and (b), for velocity $f^{\prime} (\chi )$ and temperature $\theta (\chi )$ profiles, respectively. As the molecular motion enhances due to melting heat transit, owing to the fact that ${M}_{e}$ leads to increased molecular motion from hot fluid toward cold surface, which in turn causes an increase in velocity. In contrast, with the temperature field, convective flow causes heat transfer to the melting surface more promptly, which results in decreasing temperature $\theta (\chi ).$ Figures 4(a) and (b) depict the effect of the local inertia coefficient $F$ on velocity $f^{\prime} (\chi )$ and temperature $\theta (\chi )$ of fluid flow. It can be observed that local inertia coefficient $F$ has an inverse influence on both fields. Here, velocity falls whereas temperature increases. As porous media cause resistance in a fluid flow, the result is a reduction in dimensionless velocity $f^{\prime} (\chi ).$ Figure 5 illustrates the behavior of the temperature ratio parameter ${\theta }_{w}$ on fluid temperature. A retarding effect of ${\theta }_{w}$ on dimensionless temperature can be observed. The effects of porosity parameter $\lambda $ on radial velocity are shown in figure 6, where velocity falls for higher values of $\lambda .$ The porosity of spongy media causes high resistivity to fluid flow; hence, velocity declines. Figures 7(a) and (b) show the impact of the nanoparticle volume fraction parameter $\phi $ on the radial velocity $f^{\prime} (\chi )$, and temperature $\theta (\chi )$ of the fluid. An increase in the quantity of nanoparticles in the base fluid (water) leads to a thickening of the fluid. Due to this, velocity reduces (figure 7(a)). On the other hand, in figure 7(b), for ordinary fluid (i.e., $\phi =0,$ in the absence of volume proportion) the temperature field is highest when the disks are driven further from each other (i.e., $Sq=1,2,\mathrm{3...}$), while augmentation of volume fraction $\phi $ causes a reduction in the temperature of the fluid. This is because the thermal conductivity of nanofluid increases by using a small concentration of nanoparticles, and when we increase the volume of nanoparticles, the thermal conductivity of the nanofluid decreases; hence, the temperature decreases. Figure 8 demonstrates the effects of the radiation parameter ${R}_{d}$ on the temperature $\theta (\chi ).$ The temperature drops for larger values of ${R}_{d}.$ The transfer of energy to the fluid declines owing to higher estimates of ${R}_{d},$ thus decreasing the fluid temperature. In figure 9, the impact of the thermal relaxation time coefficient $\gamma $ is shown. Higher values of $\gamma $ cause an increment in dimensionless temperature $\theta (\chi )$ for both types of CNTs. Figure 10 indicates that concentration $H(\chi )$ decreases for higher values of ${k}_{1}.$ The same outcome can be detected in figure 11 for ${k}_{2}.$ It is therefore deduced that concentration eventually reduces as the reactants are used throughout homogeneous–heterogeneous reactions. The effects of Schmidt number $Sc$ are portrayed in figure 12. The concentration profile is reduced for increasing values of $Sc.$ As $Sc$ is the ratio of momentum to mass diffusivity, greater $Sc$ estimates indicate lower mass diffusivity, which causes a reduction in fluid concentration $H(\chi ).$ In figure 13, a retardation effect of porosity parameter $\lambda $ on local inertia coefficient $F$ can be observed for surface drag force. Influences of melting parameter ${M}_{e}$ and squeezing parameter ${S}_{q}$ on the rate of heat transfer are depicted in figure 14. Here, it can be observed that augmentation in the melting coefficient ${M}_{e}$ causes an increment in the rate of heat transfer. Since molecular motion increases with high melting heat transmission, hence the rate of heat is increased when we augment the melting parameter. Figure 15 illustrates the impact of radiation coefficient ${R}_{d}$ and temperature ratio coefficient ${\theta }_{w}$ on the Nusselt number. The rate of heat transfer falls for larger estimates of ${R}_{d}.$ This is because the energy from radiation phenomena is being used in the melting .process; thus, a significant decay in the Nusselt number can be seen.

Figure 2.

New window|Download| PPT slide
Figure 2.(a). Variations of squeezing parameter ${S}_{q}$ on radial velocity $f^{\prime} (\chi ).$ (b). Variations of squeezing parameter ${S}_{q}$ on temperature distribution $\theta (\chi ).$

Figure 3.

New window|Download| PPT slide
Figure 3.(a). Variations of melting parameter ${M}_{e}$ on radial velocity $f^{\prime} (\chi ).$ (b). Variations of melting parameter ${M}_{e}$ on temperature distribution $\theta (\chi ).$

Figure 4.

New window|Download| PPT slide
Figure 4.(a). Variations of local inertia coefficient $F$ on radial velocity $f^{\prime} (\chi ).$ (b). Variations of local inertia coefficient $F$ on temperature distribution $\theta (\chi ).$

Figure 5.

New window|Download| PPT slide
Figure 5.Variations of temperature ratio parameter ${\theta }_{w}$ on temperature distribution $\theta (\chi ).$

Figure 6.

New window|Download| PPT slide
Figure 6.Variations of porosity parameter $\lambda $ on radial velocity $f^{\prime} (\chi ).$

Figure 7.

New window|Download| PPT slide
Figure 7.(a). Variations of nanoparticle volume fraction $\phi $ on radial velocity $f^{\prime} (\chi ).$ (b). Variations of nanoparticle volume fraction $\phi $ on temperature distribution $\theta (\chi ).$

Figure 8.

New window|Download| PPT slide
Figure 8.Variations of radiation coefficient ${R}_{d}$ on temperature distribution $\theta (\chi ).$

Figure 9.

New window|Download| PPT slide
Figure 9.Variations of thermal relaxation parameter $\gamma $ on temperature distribution $\theta (\chi ).$

Figure 10.

New window|Download| PPT slide
Figure 10.Variations of homogeneous reaction Coefficient ${k}_{1}$ on concentration distribution $H(\chi ).$

Figure 11.

New window|Download| PPT slide
Figure 11.Variations of heterogeneous reaction coefficient ${k}_{2}$ on concentration distribution $H(\chi ).$

Figure 12.

New window|Download| PPT slide
Figure 12.Variations of Schmidt number ${S}_{c}$ on concentration distribution $H(\chi ).$

Figure 13.

New window|Download| PPT slide
Figure 13.Variations of porosity parameter $\lambda $ and local inertia parameter $F$ on surface drag force.

Figure 14.

New window|Download| PPT slide
Figure 14.Variations of melting parameter ${M}_{e}$ and squeezing parameter ${S}_{q}$ on Nusselt number.

Figure 15.

New window|Download| PPT slide
Figure 15.Variations of radiation parameter ${R}_{d}$ and temperature ratio parameter ${\theta }_{w}$ on Nusselt number.


Table 2 depicts the comparative estimates of surface drag force for various values of squeezing parameter and Hartmann number, with Lu et al [15] as limiting case. A remarkable resemblance is achieved in this regard.


Table 2.
Table 2.Comparative results of surface drag force for distinct values of squeezing parameter and Hartmann number with Lu et al [15].

${S}_{q}$	$Ha$	Lu et al [15]	Present results
0.5	0.0	−3.146 1941	−3.134 6178
	1.0	−3.194 0816	−3.183 3610
	2.0	−3.241 3602	−3.236 7145
−1.0	1.0	−2.759 6174	−2.747 1465
0.0		−3.049 6468	−3.031 1268
1.0		−3.338 1297	−3.327 8263

New window|CSV

Concluding remarks

In the presented model, 2D time-dependent magnetohydrodynamic squeezing nanofluid flow between two parallel disks with suspended carbon nanotubes is discussed. The analysis is performed under conditions of non-linear thermal radiation, amalgamated with melting heat, and homogeneous–heterogeneous reactions. Cattaneo-Christov heat flux is engaged in place of the conventional Fourier law of heat conduction. The proposed model is transformed into a non-linear form and processed by means of the default Finite difference method of the bvp4c built-in function of the MATLAB scheme. The key findings of the presented model are:•The local inertia coefficient has an opposite impact on the radial velocity and temperature field.
•Porous media creates resistance in a fluid flow, resulting in a reduction in fluid velocity and an augmentation in the temperature of the fluid.
•The melting parameter has a retarding effect on temperature, whereas radial velocity increases.
•An increase in the radiation coefficient leads to a reduction in the temperature of the fluid.
•Radiative energy is being used in the melting process; consequently, the temperature of the fluid decreases.
•The porosity coefficient has a retarding influence on radial velocity and surface drag force.

Acknowledgments

This work is supported by Bahria University, Islamabad, Pakistan.

Reference View Option By original order
By published year
By cited within times
By Impact factor

[1]Iijima S 1991 Helical microtubules of graphitic carbon
Nature 354 56
DOI:10.1038/354056a0 [Cited within: 1]

[2]Endo M 2004 Applications of carbon nanotubes in the twenty-first century
Phil. Trans. R. Soc. A 362 1823
DOI:10.1098/rsta.2004.1437 [Cited within: 1]

[3]Mohyud-Din S T Khan S I 2016 Nonlinear radiation effects on squeezing flow of a Casson fluid between parallel disks
Aerosp. Sci. Technol. 48 186 192
DOI:10.1016/j.ast.2015.10.019 [Cited within: 1]

[4]Haq R U Hammouch Z Khan W A 2016 Water-based squeezing flow in the presence of carbon nanotubes between two parallel disks
Therm. Sci. 20 1973 1981
DOI:10.2298/TSCI141102148H [Cited within: 1]

[5]Hayat T Muhammad K Muhammad T Alsaedi A 2018 Melting heat in radiative flow of carbon nanotubes with homogeneous-heterogeneous reactions
Commun. Theor. Phys. 69 441
DOI:10.1088/0253-6102/69/4/441 [Cited within: 1]

[6]Stefan M J 1874 ‘Versuch uber die scheinbare adhesion', Akademie Der Wissenschaften in Wien
Math. Natur. 69 713 721
[Cited within: 1]

[7]Leider P J Bird R B 1974 Squeezing flow between parallel disks, I. Theoretical analysis
Ind. Eng. Chem. Fund. 13 336 341
DOI:10.1021/i160052a007 [Cited within: 1]

[8]Qayyum A Awais M Alsaedi A Hayat T 2012 Unsteady squeezing flow of Jeffrey fluid between two parallel disks
Chin. Phys. Lett. 29 034701
DOI:10.1088/0256-307X/29/3/034701 [Cited within: 1]

[9]Hayat T Haider F Muhammad T Alsaedi A 2018 Darcy–Forchheimer squeezed flow of carbon nanotubes with thermal radiation
J. Phys. Chem. Solids 120 79 86
DOI:10.1016/j.jpcs.2018.04.016 [Cited within: 1]

[10]Hashmi M M Hayat T Alsaedi A 2012 On the analytic solutions for squeezing flow of nanofluid between parallel disks
Nonlinear Anal. Model. Control 17 418 430
DOI:10.15388/NA.17.4.14048 [Cited within: 1]

[11]Cattaneo C 1948 Sulla conduzione del calore
Atti Sem Mat Fis Univ Modena 3 83 101
[Cited within: 1]

[12]Tibullo V Zampoli V 2011 A uniqueness result for the Cattaneo– Christov heat conduction model applied to incompressible fluids
Mech. Res. Commun. 38 77 79
DOI:10.1016/j.mechrescom.2010.10.008 [Cited within: 1]

[13]Christov C I 2009 On frame indifferent formulation of the Maxwell–Cattaneo model of finite-speed heat conduction
Mech. Res. Commun. 36 481 486
DOI:10.1016/j.mechrescom.2008.11.003 [Cited within: 1]

[14]Dogonchi A S Chamkha A J Seyyedi S M Ganji D D 2018 Radiative nanofluid flow and heat transfer between parallel disks with penetrable and stretchable walls considering Cattaneo–Christov heat flux model
Heat Transf - Asian Res. 47 735 753
DOI:10.1002/htj.21339 [Cited within: 1]

[15]Lu D Li Z Ramzan M Shafee A Chung J D 2019 Unsteady squeezing carbon nanotubes-based nano-liquid flow with Cattaneo–Christov heat flux and homogeneous–heterogeneous reactions
Appl. Nanosci. 9 169 178
DOI:10.1007/s13204-018-0899-1 [Cited within: 4]

[16]Zubair M Shah Z Islam S Khan W Dawar A 2019 Study of three-dimensional Darcy–Forchheimer squeezing nanofluid flow with Cattaneo–Christov heat flux based on four different types of nanoparticles through entropy generation analysis
Adv. Mech. Eng. 11 1 17
DOI:10.1177/1687814019851308 [Cited within: 1]

[17]Darcy H 1856 Les Fontaines Publiques de la ville de Dijon: Exposition et ApplicationVictor Dalmont
[Cited within: 1]

[18]Forchheimer P 1901 Wasserbewegung durch boden
Z. Ver. Deutsch, Ing. 45 1782 1788
[Cited within: 1]

[19]Muskat M 1938 The flow of homogeneous fluids through porous media
Soil Sci. 46 169
DOI:10.1097/00010694-193808000-00008 [Cited within: 1]

[20]Nasir S vShah Z Islam S Khan W Bonyah E Ayaz M Khan A 2019 Three dimensional Darcy-Forchheimer radiated flow of single and multiwall carbon nanotubes over a rotating stretchable disk with convective heat generation and absorption
AIP Adv. 9 035031
DOI:10.1063/1.5087468 [Cited within: 1]

[21]Jha B K Kaurangini M L 2011 Approximate analytical solutions for the nonlinear Brinkman-Forchheimer-extended Darcy flow model
Appl. Math. 2 1432 1436
DOI:10.4236/am.2011.212202 [Cited within: 1]

[22]Kaladhar K 2017 Double stratification effects on mixed convection flow of couple stress fluid in a non-Darcy porous medium with heat and mass fluxes
Comput. Appl. Math. 36 611 626
DOI:10.1007/s40314-015-0248-x [Cited within: 1]

[23]Sheikholeslami M Ellahi R Shafee A Li Z 2019 Numerical investigation for second law analysis of ferrofluid inside a porous semi annulus
Int. J. Numer. Method. H 29 1079 1102
DOI:10.1108/HFF-10-2018-0606 [Cited within: 1]

[24]Alamri S Z Ellahi R Shehzad N Zeeshan A 2019 Convective radiative plane Poiseuille flow of nanofluid through porous medium with slip: an application of Stefan blowing
J. Mol. Liq. 273 292 304
DOI:10.1016/j.molliq.2018.10.038 [Cited within: 1]

[25]Sheikholeslami M 2019 New computational approach for exergy and entropy analysis of nanofluid under the impact of Lorentz force through a porous media
Comput. Method. Appl. M 344 319 333
DOI:10.1016/j.cma.2018.09.044 [Cited within: 1]

[26]Li Z Shafee A Ramzan M Rokni H B Al-Mdallal Q M 2019 Simulation of natural convection of Fe₃O₄-water ferrofluid in a circular porous cavity in the presence of a magnetic field
Eur. Phys. J Plus. 134 77
DOI:10.1140/epjp/i2019-12433-5

[27]Bhatti M M Zeeshan A Ellahi R Bég O A Kadir A 2019 Effects of coagulation on the two-phase peristaltic pumping of magnetized Prandtl biofluid through an endoscopic annular geometry containing a porous medium
Chin. J. Phys. 58 222 234
DOI:10.1016/j.cjph.2019.02.004

[28]Li Z Ramzan M Shafee A Saleem S Al-Mdallal Q M Chamkha A J 2019 Numerical approach for nanofluid transportation due to electric force in a porous enclosure
Microsyst. Technol. 25 2501 2514
DOI:10.1007/s00542-018-4153-2 [Cited within: 1]

[29]Chaudhary M A Merkin J H 1995 A simple isothermal model for homogeneous heterogeneous reactions in boundary-layer flow. I Equal diffusivities
Fluid Dyn. Res. 16 311
DOI:10.1016/0169-5983(95)00015-6 [Cited within: 1]

[30]Khan U Ahmed N Mohyud-Din S T 2017 Numerical investigation for three-dimensional squeezing flow of nanofluid in a rotating channel with lower stretching wall suspended by carbon nanotubes
Appl. Therm. Eng. 113 1107 1117
DOI:10.1016/j.applthermaleng.2016.11.104 [Cited within: 1]

[31]Rosseland S 1931 Astrophysics and Nuclear-Theoretical Foundations41 44
[Cited within: 1]