The present article investigates the squeezing flow of two types of nanofluids such as Cu-water and Cu-kerosene between two parallel plates in the presence of magnetic field. The governing non-linear partial differential equations are transformed into ordinary differential equations by applying suitable similarity transformation and then solved numerically using RK-4 method with shooting technique and analytically using differential transformation method (DTM). The influence of arising relevant parameters on flow characteristics has been discussed through graphs and tables. A comparative study has been taken into account between existing results and present work and it is found to be in excellent harmony.
Squeezing flow; Nanofluid; Magnetic field; Differential transformation method
The squeezing flow of Newtonian and non-Newtonian fluids continues to stimulate significant interest by the researchers owing to increasing application in different fields of engineering, technology such as polymer processing, transient loading of mechanical components, compression, the squeezed film in power transmission, moving pistons, and chocolate filler. Such flows are performed between two moving parallel plates. Classical work in this sector was first reported by Stefan [1] considering the squeezing flow using lubrication approach in 1874. After that Reynolds [2] considered the problems for elliptic plates in 1886. Archibald [3] analyzed the squeezing flow through rectangular plates in 1956. Gradually day by day considerable efforts have been devoted by the researchers to avail the squeezing flow easier and reliable. Rashidi et al. [4] discussed the two dimensional axisymmetric squeezing flow between parallel plates. Effects of magnetic field in the squeezing flow between infinite parallel plates were reported by Siddiqui et al. [5]. Hamdan and Baron [6] investigated the squeezing flow of dusty fluid between parallel disks and discussed the squeezing effects on the velocity profiles. Domairry and Aziz [7] delivered the approximate analytic solution for the squeezing flow of viscous fluid between parallel disks with suction or blowing. Hayat et al. [8] extended the work of Domairry and Aziz [7] to analyze the squeezing flow of non-Newtonian fluids by taking second grade fluids.
Recently the study of heat transfer characteristics of squeezing flow of viscous fluid has gained considerable attention owing to its various applications in many branches of science and engineering. The squeezing flow through a porous surface has been addressed by Mahmood et al. [9]. Investigation reveals that magnitude of local Nusselt number increases with Prandtl number. The results are satisfied with Mustafa et al. [10]. Heat transfer characteristics in a squeezing flow between parallel disks has been studied by Duwairi et al. [11]. Khaled and Vafai [12] analyzed the hydromagnetic effects on flow and heat transfer over a horizontal surface placed in an externally squeezed free stream. They found that Nusselt number and wall shear stress are both increasing functions of the magnetic parameter. Analytical investigation of the unsteady squeezing flow of viscous Jeffery fluid between parallel disks was performed by Qayyum et al. [13] and discussed the porosity and squeezing effects on the velocity profiles.
The modern age is rightly called the age of science. Whenever we cast our eyes, we see that ongoing advanced technology requires further improvement of heat transfer from energy saving point of view. Since the conventional heat transfer fluids such as water, kerosene, and ethylene glycol have low thermal conductivity, modern science bestowed blessing in the form of nanofluid in which nano-sized particles are added to the base fluid to increase the heat transfer capabilities of the base fluid. It should be noticed that there have been published several recent papers [14], [15], [16], [17] and [18] on the mathematical and numerical modeling of convective heat transfer mechanism in nanofluids. Mandy [19] investigated mixed convection flow and heat transfer of nanofluids due to an unsteady stretching sheet. These models have some advantage over experimental studies due to many factors that influence nanofluid properties. Recently, nanofluid flow and heat transfer characteristics under different flow configuration were discussed by Hatami et al. [20] and [21]. Analytical investigation of unsteady squeezing nanofluid flow has been discussed by Pourmehran et al. [22]. They showed that highest value of Nusselt number can be obtained by selecting silver as nanoparticle. All the abovementioned models are single phase model. Recently a new dimension has been added by [23] and [24] in terms of complex geometry or two phase model for simulating nanofluid.
Nowadays modern science and engineering fields are blessed with so many numerical and analytical methods to obtain accurate approximate solutions from non-linear equations. The methods are variation of parameter method (VPM) [25] and [26], differential transformation method (DTM) [27], [28], [29], [30] and [31], adomian decomposition method (ADM) [32] and [33], homotopy analysis method (HAM) [34] and [35], homotopy perturbation method (HPM) [36], etc. In this article we have applied relatively a novel analytic technique, say differential transformation method (DTM) which can be said as a modified version of Taylor series method. Zhou [37] was the first to initiate the differential transformation method (DTM). It is the method in which we determine the coefficient of Taylor series of the function by solving the recursive equation from the given differential equation. While the Taylor series method requires very time consuming symbolic computation and derivatives of the data function for higher order equations, DTM plays a powerful authority to obtain the solution in a purely numerical way providing a smooth, functional form within a very few steps. Another essential advantage of this DTM method is that it reduces the size of computational work providing the desired accuracy together with fast convergence rate. More or less every semi analytical technique is helpful to obtain solutions but DTM is free from any linearization, perturbation, discretization or restrictive assumption like HAM, HPM, etc. That’s why it is not affected due to any round off errors. Hence it can be applied directly to linear and non-linear ordinary differential equations as a promising tool.
Motivated by the above investigations the present paper deals with the squeezing flow of Cu-water and Cu-kerosene nanofluid considering the presence of external applied magnetic field. Similarity transformation has been used to obtain ordinary differential equations from the governing equations. The reduced ordinary differential equations have been solved by DTM as well as by RK-4 method with shooting technique.
Consider the viscous incompressible nanofluid flow and heat transfer in a two dimensional co-ordinate system. The co-ordinate system is chosen in such a way that x-axis is measured along the plate and y axis is normal to the plate as shown in Fig. 1. The squeezing flow has been performed through a system having two parallel plates situated at Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle h(t)=H{\left(1-\alpha t\right)}^{1/2}}
distance apart where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \alpha >0} refers to the squeezing movement of both plates with velocity Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle v(t)=\frac{dh}{dt}} until they touch each other at Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle t=\frac{1}{\alpha }}
. Also Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \alpha <0}
refers to movement of the plate away to each other and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \alpha } is called the characteristic parameter having dimension of time inverse. H is the initial position of the plate at time t = 0. A uniform magnetic field of strength Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle B(t)=B_0{\left(1-\alpha t\right)}^{-1/2}} is applied normal to the plate where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle B_0} is the initial intensity of the magnetic field. In the mathematical formulation scheme we proceed with the 0 following assumption that there is no chemical reaction, radiative heat transfer, nanoparticles and base fluid are in thermal equilibrium and no slip occurs between them. All body forces are assumed to be neglected. Here we have considered two types of nanofluids cu-water and cu-kerosene. The thermophysical properties of the nanofluid are given in Table 1.
|
Figure 1. Physical model of the problem. |
Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \rho }
(kg/m3) |
Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle C_p}
(J/kg K) |
Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \kappa }
(W/m K) | |
---|---|---|---|
Pure water | 997.1 | 4179 | 0.613 |
Kerosene | 783 | 2090 | 0.145 |
Copper | 8933 | 385 | 401 |
Under the above stated situation the governing equations are as follows:
|
(1) |
|
(2) |
|
(3) |
|
(4) |
where u, v are the velocity components in x and y directions respectively. Here Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle T\mbox{,}p\mbox{,}{\rho }_{nf}\mbox{,}{\mu }_{nf}\mbox{,}{\left(\rho C_p\right)}_{nf}\mbox{,}k_{nf}}
represents the temperature, pressure, effective density, effective dynamic viscosity, effective heat capacity and effective thermal conductivity of the nanofluids respectively. Now the relations between base fluid and nanoparticles are given by
|
(5) |
|
(6) |
|
(7) |
|
(8) |
The relevant boundary conditions for the present problem are as follows:
|
(9) |
Now in order to convert the partial differential equation into ordinary differential equations, let us introduce the dimensionless functions Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle f\mbox{,}\theta }
and similarity variable Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \eta } as follows:
|
(10) |
where prime denotes the differentiation with respect to Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \eta } . Using the above transformation and simplifying Eqs. (2) and (3) by eliminating pressure terms we have
|
(11) |
|
(12) |
with the association of the boundary conditions
|
(13) |
where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle A_1=(1-\phi )+\phi \frac{{\rho }_s}{{\rho }_f}\mbox{,}A_2=} Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): (1-\phi )+\phi \frac{{\left(\rho C_p\right)}_s}{{\left(\rho C_p\right)}_f}\mbox{,}A_3= Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): \frac{k_s+2k_f-2\phi (k_f-k_s)}{k_s+2k_f+\phi (k_f-k_s)}\mbox{,} Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle S=\frac{\alpha H^2}{2{\nu }_f}}
denotes the squeeze number having the property S > 0 corresponds to the plates moving apart and S < 0 corresponds to the plates moving to each other, Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle Pr=\frac{{\mu }_f{\left(\rho C_p\right)}_f}{{\rho }_fk_f}} is the Prandtl number, Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle Ec=\frac{{\rho }_f}{{\left(\rho C_p\right)}_fT_H}{\left[\frac{\alpha x}{2(1-\alpha t)}\right]}^2} is the modified Eckert number, Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle M=\frac{{\sigma }_{nf}}{{\mu }_{nf}}H^2\alpha B_0^2} is the Hartman number, and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \delta =\frac{H(1-\alpha t)}{x}} is the dimensionless length.
For practical purpose the skin friction coefficient and Nusselt number can be defined respectively as follows:
|
(14) |
In terms of dimensionless variables the reduced skin friction coefficient and the reduced Nusselt number are respectively given by
|
(15) |
where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle {Re}_x=\frac{x\alpha H}{2{\nu }_f}}
is the local Reynolds number.
We can define a function Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle g(\lambda )}
in differential transformation method as follows:
|
(16) |
where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle G(k)}
is called the T-function or simply transformed form of Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle g(\lambda )}
. The inverse transformation of Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle G(k)}
is defined as
|
(17) |
Now putting the value of Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle G(k)}
in (17) as we defined earlier in (16) we can rewrite (17) as
|
(18) |
In the upcoming solution procedure profile we have used some of the basic operations of DTM which we have listed in Table 2.
Original function | T-function |
---|---|
Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle h(\lambda )=\frac{d^mg(\lambda )}{d{\lambda }^m}} | Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle H[k]=\frac{(k+m)!}{k!}G[k+m]} |
Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle h(\lambda )=g(\lambda ).p(\lambda )} | Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle H[k]={\sum }_{n=0}^kG[n]P[k-n]} |
Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle h(\lambda )={\lambda }^n} | Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle H[k]=\delta (k-n)=\begin{array}{l} 1\mbox{,}k=n\\ 0\mbox{,}k\not =n \end{array}} |
Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle h(\lambda )=\alpha g(\lambda )\pm \beta p(\lambda )} | Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle H[k]=\alpha G[k]\pm \beta P[k]} |
Now applying DTM in Eq. (11), we can obtain
|
(19) |
where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \Delta (m)=\begin{array}{l} 1\mbox{,}m=1\\ 0\mbox{,}m\quad \not =\quad 1 \end{array}}
The corresponding boundary conditions are transformed as
|
(20) |
Similarly, we have from (12)
|
(21) |
where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \Delta (m)=\begin{array}{l} 1\mbox{,}m=1\\ 0\mbox{,}m\quad \not =\quad 1 \end{array}}
The relevant boundary conditions become
|
(22) |
where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle F[k]\mbox{,}\Theta [k]}
are the transformed form of Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle f(\eta )\mbox{,}\theta (\eta )} and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle k_1\mbox{,}k_2\mbox{,}k_3} are the constants which can be obtained using the boundary condition. Now the iteration schemes are as follows:
|
(23) |
and
|
(24) |
Substituting the above Eqs. (23) and (24) into the main equation based on DTM we can obtain the following solutions:
|
(25) |
|
(26) |
Considering Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \eta =1\mbox{,}S=0.5\mbox{,}Pr=0.71\mbox{,}Ec=0.5\mbox{,}M=} Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): 0.5\mbox{,}\delta =0.1
one can easily get Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle k_1=1.460028967\mbox{,}k_2=-0.420057934\mbox{,}k_3=}
Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): 1.016183698\mbox{.}
Now we proceed to solve numerically the Eqs. (11) and (12) together with the boundary conditions (13) by switching them to an initial value problem. The Eqs. (11) and (12) can be written in the form of first order differential equations as follows:
|
(27) |
and
|
(28) |
Since the Eqs. (27) and (28) are initial value problem; hence, to integrate them we need a value for Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle z(0)}
i.e., Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle f^{{'}}(0)\mbox{,}}
Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle p(0)}
i.e., Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle f^{{''}}(0)\mbox{,}}
Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle q(0)}
i.e., Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle f^{{'''}}(0)} and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle r} i.e., Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle {\theta }^{{'}}(0)\mbox{.}} We know the value of Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle {\theta }^{{'}}(0)} and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle f^{{''}}(0)} from the boundary condition. But no such values for Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle f^{{'}}(0)} and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle f^{{'''}}(0)} are given in the boundary conditions. So we make some proper guess of Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle f^{{'}}(0)} and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle f^{{'''}}(0)\mbox{.}} After that integration is carried out and compared with the values of Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle f(1)\mbox{,}f^{{'}}(1)\mbox{,}\theta (1)\mbox{.}} We adjust the value of Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle f^{{'}}(0)\mbox{,}}
Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle f^{{'''}}(0)}
to obtain better approximation of the solution. We take the series of values for Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle f^{{'}}(0)\mbox{,}}
Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle f^{{'''}}(0)}
and apply RK-4 method with different values of step-sizes such as Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle h=0.01\mbox{,}0.001\mbox{,}\ldots \mbox{.}} etc. The inner iteration is done with the convergence criterion of Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle {10}^{-6}} in all cases.
To check the efficiency of RK-4, a comparison between analytical solution obtained by DTM and numerical solution obtained by RK-4 is tabulated in Table 3. Here we have calculated the values of Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle f^{{''}}(1)}
of Cu-water for various values of squeeze number S. The results obtained are in good agreement. On the other hand to check the validity of our present work we consider a situation of Cu-water nanofluid in the absence of magnetic field when two plates are moving apart i.e., S > 0. We arrange the situation in mathematical equation by letting the values of parameter as Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle S=0.5\mbox{,}\delta =0.1\mbox{,}M=0.0\mbox{,}\phi =}
Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): 0.0 . We have calculated the values of Nusselt number Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle -{\theta }^{{'}}(1)}
for various values of Pr and Ec and tabulated these values in Table 4. We see that the values are in excellent agreement with Mustafa et al. [10] and Pourmehran et al. [22].
S | Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle f^{{''}}(1)} | |
---|---|---|
DTM | RK-4 | |
0.5 | −3.413019 | −3.413020 |
1.0 | −3.736838 | −3.736839 |
1.5 | −4.030501 | −4.030502 |
2.0 | −4.030037 | −4.030038 |
2.5 | −4.551056 | −4.551057 |
Pr | Ec | Mustafa et al. [10] | Pourmehran et al. [22] | Present work |
---|---|---|---|---|
0.5 | 1.0 | 1.522368 | 1.518859607 | 1.522367498 |
1.0 | – | 3.026324 | 3.019545607 | 3.026323559 |
2.0 | – | 5.98053 | 5.967887511 | 5.980530398 |
5.0 | – | 14.43941 | 14.41394678 | 14.43941324 |
1.0 | 0.5 | 1.513162 | 1.509772834 | 1.513161807 |
– | 1.2 | 3.631588 | 3.623454726 | 3.631588269 |
– | 2.0 | 6.052647 | 6.039091204 | 6.052647108 |
– | 5.0 | 15.13162 | 15.09772808 | 15.13161784 |
In this section we will discuss the effects of various physical parameters on velocity distribution and temperature distribution. The whole discussion has been performed through graphs and tables. In the simulation the default values of the parameters are taken as Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle Pr=0.71\mbox{,}Ec=0.5\mbox{,}M=0.5\mbox{,}S=0.5\mbox{,}\delta =} Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): 0.1\mbox{,}\phi =0.15
unless otherwise specified.
Fig. 2 demonstrates the influence of the squeezing parameter S on nanofluid velocity. It reveals that for Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle 0\leqslant \eta \leqslant 0.42}
(not accurately determined) the velocity profile decreases with increasing S but for Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \eta >0.42} (not accurately determined) it starts increasing as shown in the figure. Physically this can be explained that the decrease in the fluid velocity near the wall region causes to increase velocity gradient there. Since the mass flow rate is kept conservative hence decrease in the fluid velocity near the wall region will be compensated by the increasing fluid velocity near the central region. That’s why at Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \eta =0.42} (not accurately determined) we found a point of separation and backflow occurs. Also the rate of decrease in the fluid velocity for Cu-water is quicker than Cu-kerosene within the region Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle 0\leqslant \eta \leqslant 0.42}
. But the effect is opposite for Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \eta >0.42\mbox{.}}
From Table 5 we see that the reduced skin friction coefficient Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle {Cf}_r} decreases as S increases for both Cu-water and Cu-kerosene.
|
Figure 2. Effects of squeezing parameter S on velocity. |
S | Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle {Cf}_r} | Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle {Nu}_r} | ||
---|---|---|---|---|
Cu-water | Cu-kerosene | Cu-water | Cu-kerosene | |
−1.0 | −1.938344 | −1.819084 | 1.225071 | 1.232899 |
0.5 | −3.096871 | −2.983827 | 1.137759 | 1.138204 |
1.0 | −3.390695 | −3.275727 | 1.119464 | 1.111067 |
Fig. 3 shows that the nanofluid temperature decreases as S increases. Since Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle S=\frac{\alpha H^2}{2{\nu }_f}}
hence increase in the squeezing parameter S represents the decrease in kinematic viscosity as well as the speed at which the plates move apart. Consequently the friction between the boundary surface of the plates and nanoparticles reduces and so the temperature of the nanofluid. This result is found to be identical with Mustafa et al. [10]. Table 5 shows that Nusselt number Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle {Nu}_r} decreases with increasing S for both Cu-water and Cu-kerosene.
|
Figure 3. Effects of squeezing parameter S on temperature. |
Figure 4(a) and Figure 4(b) depict the effect nanoparticle volume fraction φ on nanofluid velocity for S > 0 and S < 0. It has been found that for Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle 0\leqslant \eta \leqslant 0.42}
(not accurately determined) velocity increases as φ increases but the effect is opposite for Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \eta >0.42} (not accurately determined). The nature of the effect is same for both for S > 0 and S < 0. But when two plates are moving apart i.e., S > 0 then the velocity of Cu-kerosene is less as compared to Cu-water within the region Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle 0\leqslant \eta \leqslant 0.42\mbox{.}} However the effect is opposite and prominent when S < 0 i.e., so-called squeezing flow. From Table 6 it is seen that the influence of φ increases the reduced skin friction coefficient Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle {Cf}_r} for both Cu-water and Cu-kerosene.
|
Figure 4(a). Effect of nanoparticle volume fraction on velocity when S < 0. |
|
Figure 4(b). Effect of nanoparticle volume fraction on velocity when S > 0. |
Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \phi } | M | Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle {Cf}_r} | |||
---|---|---|---|---|---|
Cu-water | Cu-kerosene | ||||
S = 0.5 | S = −0.5 | S = 0.5 | S = −0.5 | ||
0.02 | 0.5 | −3.096871 | −2.388320 | −2.983826 | −2.274793 |
0.04 | – | −2.889905 | −2.180419 | −2.711130 | −2.000632 |
0.06 | – | −2.736134 | −2.025792 | −2.515986 | −1.804083 |
0.02 | 1.5 | −3.402609 | −2.773329 | −3.275986 | −2.646386 |
– | 3 | −4.277192 | −3.811202 | −4.112513 | −3.646416 |
– | 4.5 | −5.414585 | −5.074220 | −5.201616 | −4.861222 |
We have represented the temperature profiles against η for various values of nanoparticle volume fraction φ in Figure 5(a) and Figure 5(b) for both S < 0 and S > 0 respectively. Fig. 5(a) shows that the fluid temperature possesses high value because of increasing values of φ for S < 0 i.e., the squeezing flow. The reason behind the phenomena is that when two plates move to each other and nanoparticle volume increases then there are more collision between nanoparticles and particles with the boundary surface of the plate and consequently the resulting friction gives rise to increase the temperature within the fluid near the boundary region. It should be noted that fluid temperature of Cu-kerosene is significantly increased as compared to Cu-water. Higher values of Prandtl number of kerosene are responsible for this phenomenon. But the result is totally opposite for S > 0 i.e., when two plates move apart. On the other hand Nusselt number increases with φ for both nanofluids as shown in Table 7.
|
Figure 5(a). Effect of nanoparticle volume fraction on temperature when S < 0. |
|
Figure 5(b). Effect of nanoparticle volume fraction on temperature when S > 0. |
Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \phi } | M | Ec | Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle {Nu}_r} | |||
---|---|---|---|---|---|---|
Cu-water | Cu-kerosene | |||||
S = 0.5 | S = −0.5 | S = 0.5 | S = −0.5 | |||
0.02 | 0.5 | 0.5 | 1.137759 | 1.160594 | 1.138203 | 1.162025 |
0.04 | – | – | 1.199843 | 1.223263 | 1.200806 | 1.226410 |
0.06 | – | – | 1.266466 | 1.290603 | 1.267997 | 1.295712 |
0.02 | 1.5 | – | 1.157179 | 1.149260 | 1.157937 | 1.149597 |
– | 3.0 | – | 1.251637 | 1.206060 | 1.252779 | 1.205117 |
– | 4.5 | – | 1.426942 | 1.376159 | 1.428062 | 1.375016 |
0.02 | 0.5 | 0.5 | 2.275519 | 2.321188 | 2.276407 | 2.324051 |
– | – | 1.5 | 3.413278 | 3.481782 | 3.414611 | 3.486075 |
– | – | 2.5 | 4.551038 | 4.642375 | 4.552815 | 4.648101 |
The impact of Hartmann number M on velocity distribution has been portrayed in Figure 6(a) and Figure 6(b) for both S > 0 and S < 0. It reveals that as M increases the velocity decreases for Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle 0\leqslant \eta \leqslant 0.42}
(not accurately determined) and it is true because for electrically conducting fluid in the presence of magnetic field there will always be Lorentz force which slows the motion of fluid in the boundary layer region. But for Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \eta >0.42} (not accurately determined) the reverse effect is seen because of the same mass flow reason as discussed in the previous Section 4.1. Moreover when S < 0 two plates are very close to each other, then the situation together with retarding Lorentz force creates adverse pressure gradient. Whenever such forces act over a long time then there might be a point of separation and back flow occurs. For S > 0 the reason is slightly different. When two plates move apart then a vacant space occurs and fluid in that region goes with high velocity so that mass flow conservation will not be violated. That’s why for Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \eta >0.42} we suddenly find an accelerated flow. For Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle 0\leqslant \eta \leqslant 0.42} velocity of Cu-water is less as compared to Cu-kerosene for S < 0 whereas the opposite effect occurs for S > 0. It should be noted from Table 6 that reduced skin friction coefficient Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle {Cf}_r} decreases as M goes high.
|
Figure 6(a). Effects of Hartmann number M on velocity when S < 0. |
|
Figure 6(b). Effects of Hartmann number M on velocity when S > 0. |
Figure 7(a) and Figure 7(b) exhibits that the temperature of the nanofluid decreases as Hartmann number increases for both S < 0 and S > 0. For S < 0 the effect is more prominent and temperature of Cu-kerosene is higher than Cu-water nanofluid. The effect is opposite for S > 0. It is worth mentioning that for S > 0 we see a crossover of temperature profile at η = 0.82 (not accurately determined). But after that region the effect is negligible. From Table 7 we see that Nusselt number increases as M increases.
|
Figure 7(a). Effects of Hartmann number M on temperature when S < 0. |
|
Figure 7(b). Effects of Hartmann number M on temperature when S > 0. |
The influence of Eckert number Ec on nanofluid temperature profile has been portrayed in Figure 8(a) and Figure 8(b) for both S > 0 and S < 0. It reveals that the temperature significantly increases with the increasing value of Ec by means of frictional heating. Also the Nusselt number increases as Ec increases for both S > 0 and S < 0 as shown in Table 7. It should be noted that Cu-kerosene possess higher temperature than Cu-water for S < 0. But the reverse effect is seen for S > 0.
|
Figure 8(a). Effects of Eckert number Ec on temperature when S < 0. |
|
Figure 8(b). Effects of Eckert number Ec on temperature when S > 0. |
In the present paper the influence of magnetic field on squeezing flow of Cu-water and Cu-kerosene nanofluids has been discussed. The model is transformed and rendered into dimensionless form and then solved analytically using differential transformation method (DTM) as well as numerically using RK-4 together with comparison to show the efficiency of DTM. The numerical discussion has been performed through graphs and tables to illustrate the details of the flow characteristics. Based on the whole discussion the main conclusions of our investigation are as follows.
It is observed that fluid velocity decreases with the influence of squeezing parameter and Hartmann number within the region Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle 0\leqslant \eta \leqslant 0.42\mbox{.}}
But fluid velocity increases within the same region due to nanoparticle volume fraction φ whereas opposite effect is seen in both cases for Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \eta >0.42\mbox{.}}
.
The authors wish to express their sincere thanks to the honorable reviewers for the valuable comments and suggestions to improve the quality of the paper.
Published on 11/04/17
Licence: Other
Are you one of the authors of this document?