Research Article  Open Access
HighAccuracy Multiscale Simulation of ThreeDimensional Squeezing Carbon NanotubeBased Flow inside a Rotating Stretching Channel
Abstract
Enhancing the heat transfer rate using nanofluids is of great interest to engineers and scientists. This research aims to study the heat and mass transfer analysis of threedimensional squeezing carbon nanotube (CNT) based nanofluid flow inside a rotating stretching channel. The upper wall of the channel is assumed to have a reciprocating movement, and the lower wall is assumed to be stationary and permeable. Also, radiative effects are taken into account using the Taylor series approximation. The momentum and energy equations are transformed into a coupled system of nonlinear ordinary differential equations utilizing similarity solutions. A new multiscale and accurate method was developed to solve the achieved nonlinear systems of equations. Water is chosen as the base fluid; singlewall carbon nanotubes (SWCNTs) and multiwall carbon nanotubes (MWCNTs) are added to it, and then two types of nanofluids were created. The effect of different variables such as the concentration of nanotubes, nanotube’s type, suction parameter, rotation parameter, squeezing number, Eckert number, and radiation parameter on the velocity and temperature profiles is investigated. Our results reveal that the temperature profile is an increasing function of the squeezing number, suction, rotation, and radiation parameters when the upper wall moves towards the lower one.
1. Introduction
Fluid squeezing is a process in which two boundaries approach each other and pressurize the trapped flow between them. This process has numerous applications in science and industry such as hydraulic systems, manufacturing of fibers and papers, lubrication of bearings, cooling systems, and biomedical engineering [1]. Wang was the first one who developed squeezing flow between parallel plates in 1976 [2]. After him, Hamza and Macdonald [3] used a finite difference solution to study the twodimensional squeezed flow in 1981. Later on, a study of magnetohydrodynamic and moving boundaries effect on a squeezed fluid film was done [4, 5]. Islam et al. analyzed the squeezing fluid flow in a porous medium channel [6], and Sherwood studied squeezing flow in the gap between nonparallel circular plates [7]. By the advent of nanotechnology and its pervasive influence in all aspects of engineering, nanofluid and its remarkable features in enhanced heat transfer were introduced by Choi and Eastman [8]. Accordingly, scholars studied the effect of adding nanoparticles to different base fluids both numerically and experimentally [9–16]. One of the most important cases that has several engineering applications is the nanofluid flow between the two walls. In our previous studies, we analyzed the heat transfer enhancement between two vertical flat plates (stationary and moving) by a high precision multiscale approach based on the Galerkin method [17–19]. Adding copper and silver nanoparticles to the base fluid helped to increase the heat transfer rate. Dib et al. [20] investigated the squeezing nanofluid flow by the Adomian decomposition method (ADM). They found out that the type of nanofluid plays an important role in heat transfer enhancement. Gupta and Saha Ray [21] developed a numerical method based on Chebyshev wavelet expansion and concluded that the Nusselt number has a direct relationship with the concentration of the nanoparticles when the fluid is squeezing. Effect of changing the base fluid was studied in [22] for the copper nanoparticles by variation of the parameters method (VPM). Xue [23] proposed a model for thermal conductivity of carbon nanotube (CNT) composites which provide higher thermal conductivity than conventional nanoparticles such as aluminium oxide, titanium oxide, and silver. The study of squeezing nanofluid flow in a rotating channel and a lower stretching porous wall was done in [24, 25] for different nanoparticles in water. Khan et al. [26] simulated the same mentioned case by use of the Runge–Kutta–Fehlberg (RKF) method for CNTs. Ghadikolaei et al. [27, 28] investigated it with the nanofluid flow of and CNTs in ethylene glycol as base fluid.
One of the best methods to solve differential equations is the Galerkin method. The Galerkin method is based on minimizing the residual with respect to the functional basis. There are some improvement methods over the standard Galerkin method; the wavelet Galerkin method (WGM) is one of the most appealing and powerful ones, thanks to the special features of the wavelet basis. One of the notable characteristics of the wavelet functions is their multiresolution analysis (MRA), and this capability provides a multiscale solution when these functions are used as a basis. On the contrary, most of the science and engineering problems are multiscale in nature and using multiscale methods in the simulation will help to a better understanding of the physics behind them.
The aspiration of the current analysis is to study the threedimensional squeezing nanofluid flow (CNTs in water) in a rotating channel with a lower permeable stretching wall by considering radiation heat transfer effect in different scales. To do so, a robust Galerkinbased multiscale method implemented. Our proposed method in this study can provide an initial coarse solution to save time and computational resources, and accordingly, the initial coarse answer can be improved by solving the problem in the higher scales. To the best of our knowledge, there is no data for the mentioned problem. To fill out the gap in the literature, we used a novel, precise, and multiscale method for the numerical simulation of the mentioned problem.
2. Problem Description
An unsteady threedimensional CNTbased nanofluid flow inside a rotating channel with a lower wall in stretching and permeable mode was considered. The lower plate of the channel is located at and is stretching in only the x direction with a velocity of . The permeation velocity at the lower wall is . The upper wall is assumed to have a reciprocating movement and squeezes the nanofluid with the velocity of . Figure 1 shows the schematic of the problem. We assumed that the fluid and nanotubes are incompressible and are trapped between the two infinite horizontal plates, no chemical reaction takes place between them, there is a thermal equilibrium between the particles and base fluid, the CNTs are suspended in the base fluid, and the noslip condition occurs between them. Singlewall carbon nanotubes (SWCNTs) and the multiwall carbon nanotubes (MWCNTs) are considered in this study.
The relationship between a flow’s kinetic energy and the boundary layer enthalpy difference can be described by the use of the Eckert number (), which is a dimensionless number. In our problem, would be showing the heat dissipation effect. Friction in the shear layer at high Eckert numbers will generate heat and will be addressed by the viscous dissipation effect. Pure water was chosen as base fluid, and then SWCNTs and MWCNTs were added to it as nanoparticles. Thermophysical properties of the base fluid and nanoparticles are tabulated in Table 1.

Considering the mentioned conditions and assumptions, the continuity, momentum, and energy equations are as follows:where u, v, and denote the velocity in the x, y and z directions, respectively. T, p, , , , and represent the temperature, pressure, effective density, effective dynamic viscosity, effective heat capacity, and radiative heat flux fluid for the nanofluid, respectively. The Rosseland approximation for the optically thick media can be used:where and are showing the Stefan–Boltzman constant and Rosseland mean spectral absorption coefficient, respectively. We assume that there is a small temperature difference in the flow and by use of Taylor series approximation about the temperature , we have
By replacing equations (6) and (7) in equation (5), we have
The relative boundary conditions arewhere means the lower wall of the channel is stretching. The nanofluid effective density and heat capacities are defined by
For effective dynamic viscosity of nanofluid, we use the approximation given by Brinkman [29] for spherical particles:
Effective thermal conductivity can be characterized by the use of different approximations. In this project, we used the model presented by Xue [23] for the carbon nanotubes:where ϕ is the volume fraction of nanoparticles and is showing the thermal conductivity of the carbon nanotubes.
For simplifying the governing equations to a nondimensional set of ordinary differential equations (ODEs), the following transformations are given:where η is the local similarity variable, and represent the velocity along the y and x axes, and and correspond to the rotational velocity and temperature profile, respectively.
By using and substituting equation (14) in equations (1)–(4) and considering equations (8)–(13), we have
The coefficient originated from approximation given by Brinkman in equation (7).
Similar transforms change the boundary conditions of the problem to a more simplified form:
The nondimensional variables in the above equations are defined as follows:where is the rotational parameter; A is the suction parameter; S is the squeezing number, and positive one means the upper wall is moving toward the lower stationary one (and vice versa); is the Prandtl number; is the Eckert number; is the modified Eckert number; and R is the radiation parameter. When there is no movement in the upper wall (), the unsteady flow will turn to a steady flow. We are using the definition of [28] for defining the Nusselt number and skin friction coefficient:
Nondimensional format of the above equations can be written as follows:where the term is called the local Reynolds number.
3. Alpert’s Multiwavelets
Alpert et al. [30, 31] introduced an orthogonal basis for and using multiresolution analysis (MRA). These bases are indexed by a parameter and built based on Legendre polynomials of degree less than r. Over the last decade, Alpert’s multiwavelets have been used in numerous areas of mathematics, engineering, computer science, statistics, physics, etc [30, 32–35].
Let , where and for . Also, suppose that and are the dilation and the translation operators, respectively, such that and .
According to the characteristics of , one can introduce a set of scaling functions for subspaces using the dilation and the translation operators, i.e.,
Based on the definition of this wavelets, it is obvious that is the collection of piecewise polynomials which are of degree less than r on the interval , . As the subspaces are nested, there exist complementary orthogonal subspaces such thatwhere represents orthogonal sumations. This gives rise to an orthogonal decomposition of , i.e.,
So there exists a family of other basis that generate the complementary subspaces , namely,
The functions are called multiwavelets. In other words, from (21), the multiwavelets , , , form a family of bases for the orthogonal complement space . On the contrary, it can be shown that where δ is the Kronecker delta function. Therefore, this family of basis is orthonormal with respect to the inner product.
Due to the fact that and , vector functions and satisfy a matrix refinement equation:where and and for are matrices with constant elements that are introduced in [30].
3.1. Multiscale Transformation
From (21), for a fixed integer , it is easy to find multiscale decomposition . Using this decomposition, we are able to approximate any function by the multiscaling functions of the coarse space and multiwavelets of the higher levels , . Now we introduce multiscale operator that maps onto viawhere and for are the orthonormal projection operators that map onto and , respectively. The multiscale operator defines the multiscale transformation, i.e.,where the coarse scale coefficients and multiwavelets coefficients are determined by the orthonormal property of Alpert’s multiwavelets:
Assume that where . In fact, for a fixed number r and J, is a vector function that includes the scaling functions in the coarse space and the highlevel spaces , . By this introduction, we are able to writewhere is an dimensional vector with entries and for , , and .
It should be mentioned that the coefficients and for , , and must be obtained to find the multiscaling transformation (27). One can calculate the singlescale coefficients, , according to the interpolating property of interpolating scaling functions [30]. But for evaluating the multiwavelets coefficients , this property does not exist. One remedy for this problem is numerically calculating the integrals. To avoid this action, we introduce a new matrix that is known as the wavelet transform matrix and is obtained using twoscale refinement relations (24) and (25). In fact, the wavelet transform matrix helps us to obtain the multiwavelets by using the scaling functions, i.e.,where is an () matrix which is obtained by the following scheme. Generally, the refinement equation of multiscaling functions for neighboring scales J and is given by where , , and is the identity matrix of size . Let the vector function satisfies for where and . Now we are able to determine the wavelet transform matrix as
One can compute the elements of matrices , , , and easily using the interpolation property and orthonormality of Alpert’s multiwavelets [30]. There is also a reconstruction formula [30]:
Relations (24), (25), and (32) yield algorithms for the transition between different scales of the multiresolution analysis.
4. Wavelet Galerkin Method (WGM)
Wavelet Galerkin methods provides an approximated solution to the differential equations by the formalism to generate discrete algorithms. Let us write the largest derivative of the unknown solution of equation (15) in terms of the basis of . Let ; by successive integration from both sides of this equation and using the operational matrix of integration which is introduced in [18, 33], we can writewhere super index N has been used for the largest derivatives and for is a known constant determined by boundary conditions. To use the wavelet Galerkin method, we expand using Alpert’s multiwavelets:
According to (33) and (34), we have
In the same way, we can find the closed form for f, , and θ and their derivatives, namely,
Using wavelet transform matrix and interpolating property of this type of wavelets, we can approximate the nonlinear terms of (15) as
We can also expand the functions , , and in the same way:where for are the vectors. Let us to look back into equation (15). Using equations (36)–(38), we have
Applying the Wavelet Galerkin method (multiplying both side in and integrating over ), we get a system of nonlinear algebraic equations that we solve it using the Newton method to find the unknown coefficients and numerical solutions.
5. Results and Discussion
Implementation of the discussed method on the obtained governing equations provides an analysis of each input parameters effect. To show the multiscale solution for the problem and as an accuracy test, Table 2 and Figure 2 are provided. As it can be seen, results are shown in different scale levels (J) and by using different orders of polynomials (r) in comparison to the reference method, fourthorder Runge–Kutta (RK method). Increasing the order of polynomials and scale levels should help to converge the results to the precise answer. Table 2 represents the accuracy of the proposed method. One can see that, by using the sixth order polynomials and just in the secondscale level, the WGM method can provide results with very high accuracy, and Figure 2 shows the effect of increasing the J and r values on the error for and in the logarithmic scale.

(a)
(b)
Effect of the concentration of nanoparticles ϕ on the heat transfer and velocity profiles is one of the most important surveys in any nanofluid study. For our case, pure water was chosen as the base fluid, and then two carbon nanotubes with different thermophysical properties, which are depicted in Table 1, were added to it. We increased the amount of ϕ, and Figure 3 shows this effect on the velocity and temperature profiles. It can be seen in Figure 3(a) that there is a direct relationship between the nanoparticle volume fraction and the velocity in the y direction () for both single and multiwall carbon nanotubes; also, this effect is more remarkable in the first half of the channel. In Figure 3(b), one can see that when η is between zero and 0.35, by increasing the nanoparticle concentration, we will have higher velocity in the x direction; however, for the rest of the η values, there is an opposite trend for velocity (). Figure 3(c) depicts the behavior of rotational velocity by increasing the nanoparticle volume fraction. One can see a reverse flow in the center of the channel. It is possible to reduce the mentioned reverse flow by adding more nanotubes to the water. One important point about the discussed figure is that, in higher concentrations, the MWCNT’s rotational velocity will change more than SWCNT’s velocity. The physical reason can be related to the fact of lower densities for SWCNTs. Effect of volume fraction on temperature profiles can be found in Figure 3(d), and the opposite trend of temperature profiles with ϕ is obvious. By adding more nanoparticles to the base fluid, the higher thermal conductivity of the nanoparticles causes an increase in the heat transfer. Accordingly, temperature profiles will drop (Figure 3(d)). By zooming in on the temperature profile, we can see higher values of θ for the MWCNTs which is due to lower thermal conductivity of these types of particles. Another important conclusion from the temperature profiles is that when we compare the differences in heat transfer amount between the pure water and 10 percent nanofluid, we see a huge difference. However, adding more nanoparticles (20 percent case) would not provide that much difference in comparison to the lower percentage nanofluid flow. In some practical applications, this issue can be useful when the nanoparticles are expensive or when adding more of them is causing some side problems.
(a)
(b)
(c)
(d)
Effect of rotational parameter on the y axis velocity is highlighted in Figure 4(a). When we raise the rotational parameters, there is a dual behavior for the . Near the lower wall, increasing the rotational parameter causes a decrease in the y component velocity. However, there is a reverse trend for the vicinity of the upper wall. Figure 4(b) depicts the change in velocity in the x direction by increasing the value. There is a direct relationship between the increasing value and the value of xdirectional velocity in the core (far from the two walls). On the contrary, one can see the reverse trend for and . Moreover, as it can be expected, the absolute values of rotational velocity increase by raising the rotational parameter (Figure 4(c)). By evaluating all the three mentioned regions in Figures 4(b) and 4(c), it can be concluded that effect of the rotational parameter on the SWCNTs are higher than MWCNTs. In Figure 4(d), we see the direct effect of the rotational parameter on the temperature profile, while the upper plate is moving toward the lower one.
(a)
(b)
(c)
(d)
To study the effect of the suction parameter A on the velocity and temperature profiles, the value of the suction parameter increased from 0.5 to 2 by an increment of 0.5. Based on the obtained results in Figure 5, as it can be predicted easily, one can see that an increase in A increases the velocity in the y direction. It is crystal clear that the increase in ydirectional velocity near the stationary wall would be maximum due to the suction at the lower wall. By getting closer to the moving wall, the effect of the suction parameter would be recessive. For the xdirectional velocity profile, there is a reverse relationship between the A and values. Moreover, a decrease in the xdirectional velocity can lead to flow reversal in the center of the domain. Figure 5(c) illustrates the effect of the suction parameter on the rotational velocity. As it can be seen, by increasing the A value from 0.5 to 2, the absolute value of rotational velocity increased and amount of increase in the SWCNTs is higher than MWCNTs due to their lower densities. The last picture of the Figure 5 resulted by changing the A value and studying its effect on the temperature profile by keeping other parameters fixed. As it can be predicted easily, near the lower plate the temperature profile would be higher because of the suction source that we have in the bottom of the geometry. By moving toward the upper plate, the effect of the suction parameter will be weakened and the temperature profile will drop. The temperature profile for the solution which contains SWCNT will be higher compared to the MWCNTs nanofluid due to the higher thermal conductivity of the SWCNTs.
(a)
(b)
(c)
(d)
For the analysis of the variation of the squeezing number and its effect on the velocity and temperature profile, a series of Figure 6 are plotted. Based on the results in Figures 6(a)–6(d), there is a direct relationship between the squeeze number S and the velocity profile in the y direction when we compress the fluid between the plates (Figure 6(a)). The mentioned trend is opposite in the decompression period (Figure 6(b)). We are witnessing the same behavior as mentioned for ydirectional velocity for values by changing the squeezing number; however, for the xdirectional velocities, maximum changes occur close to the center of the channel. Effect of the squeezing number on the rotational velocity is shown in Figures 6(e) and 6(f) for the compression and decompression cases. As we saw in previous parts, the effect of the squeezing number on the SWCNTs is higher than MWCNTs. It is worth to mention that when the upper wall is moving apart (), the effect of reverse flow in the center of the channel can be seen (Figures 6(d)–6(f)). Finally, one of the important analysis that we can have by the variation of the squeezing number is the change in the temperature profile. It is interesting that, by compression or decompression of the nanofluid, the temperature profile between them will increase. In the compression case, the solution with MWCNT is experiencing a relatively higher change; however, in the decompression mode, more changes in the temperature profile in the fluid with SWCNTs (Figures 6(g) and 6(h)) are obvious.
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)
As it is obvious and discussed previously, by increasing the heat transfer rate, the temperature profile will drop. Different methods exist to raise the rate of heat transfer including adding nanoparticles and changing the squeezing number. One of the other methods of increasing the heat transfer rate for our geometry is increasing the radiation parameter. Figure 7(a) shows this matter when the upper wall is moving toward the stationary lower one. As it has shown, increasing the radiation parameter has helped to drop the temperature profile which is an indicator of more heat transfer. By taking a closer look into the obtained results, we see that, for the SWCNT case, we have more temperature drop which has caused higher thermal conductivity of this type of nanotubes. Also, Figure 7(b) illustrates the change of the Nusselt number by the R value for both lower and upper plates, and the same mentioned trend is visible for the Nusselt number. Figure 8 compares the effect of different kinds of nanotubes, concentration, and squeezing number on the Nusselt number and skin friction coefficient. One can see in Figures 8(a) and 8(b) that increasing the concentration helps to have more heat transfer rate and bigger Nusselt number. Also, the variation of squeezing number changes the rate of heat transfer and Nusselt number. At the lower wall, increasing the value of S causes an increase in the Nusselt number; however, this relationship is reversed for the upper wall. Moreover, the SWCNTs are mostly affected more with these variations of S and ϕ values due to their higher thermal conductivities. Analysis of the deviation in the skin friction coefficient can be found in Figures 8(c) and 8(d). For the upper wall in the geometry, increasing the nanoparticle volume fraction yields higher skin friction coefficient for both types of particles; however, squeezing number plays an opposite role and causes a downfall in the values. We see a completely different trend for the values of the skin friction coefficient in the lower plate by changing the ϕ and S parameters. By scrutinizing the Figures 8(c) and 8(d), it can be concluded that values corresponding to MWCNT’s are higher on the upper wall and SWCNTs are taking higher values for the lower plate due to the difference in their densities.
(a)
(b)
(a)
(b)
(c)
(d)
The Eckert number shows the relationship between kinetic energy of the flow and enthalpy difference across the thermal boundary layer. For characterizing the heat dissipation phenomena, the effects of changing the Eckert number in the presence of the radiation parameter are shown in Figures 9(a) and 9(b) for both upper and lower plates. In the upper plate, increasing the Eckert number has a direct effect on the rate of heat transfer and Nusselt number for both the compression and decompression cases. On the contrary, increasing the values causes a downfall in the Nusselt number at the lower plate of the channel for both and cases. There is an interesting point about the obtained results, and the effect of changing the nanoparticle type is more remarkable in the decompression period and upper plate.
(a)
(b)
At the end of this section, effects of the change in the modified Eckert number and radiation parameter R are depicted in Figures 10(a) and 10(b). The results of the simulations in the mentioned scale show that there is a direct relationship between the modified Eckert number and heat transfer rate at the upper plate. For the lower plate, one can see an opposite relationship between the and .
(a)
(b)
6. Conclusion
Threedimensional CNTbased nanofluid flow in the presence of radiation heat transfer was analyzed by a new multiscale method. The 3D squeezing flow modeled by a set of nonlinear ordinary differential equations was obtained by implementing similarity solutions on the governing equations. The acquired equations were solved by a Galerkinbased spectral element method with Wavelet basis in different scales. Values of and were compared with a verified method from the literature for different scenarios to validate the algorithm and code. After verification of the method, the effect of different parameters on the velocity and temperature profiles is analyzed. The results of the simulations show that one can reach higher values in the temperature profile using MWCNTs instead of SWCNTs when the walls are moving toward each other. Moreover, the velocity components in the compression mode increase with the nanoparticle volume fraction in both solutions.
Data Availability
The raw/processed data are available from the corresponding author on a reasonable request.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
References
 J. Engmann, C. Servais, and A. S. Burbidge, “Squeeze flow theory and applications to rheometry: a review,” Journal of NonNewtonian Fluid Mechanics, vol. 132, no. 1–3, pp. 1–27, 2005. View at: Publisher Site  Google Scholar
 C.Y. Wang, “The squeezing of a fluid between two plates,” Journal of Applied Mechanics, vol. 43, no. 4, pp. 579–583, 1976. View at: Publisher Site  Google Scholar
 E. A. Hamza and D. A. Macdonald, “A fluid film squeezed between two parallel plane surfaces,” Journal of Fluid Mechanics, vol. 109, pp. 147–160, 1981. View at: Publisher Site  Google Scholar
 E. A. Hamza, “The magnetohydrodynamic effects on a fluid film squeezed between two rotating surfaces,” Journal of Physics D: Applied Physics, vol. 24, no. 4, pp. 547–554, 1991. View at: Publisher Site  Google Scholar
 P. Singh, V. Radhakrishnan, and K. A. Narayan, “Squeezing flow between parallel plates,” IngenieurArchiv, vol. 60, no. 4, pp. 274–281, 1990. View at: Publisher Site  Google Scholar
 S. Islam, H. Khan, I. A. Shah, and G. Zaman, “An axisymmetric squeezing fluid flow between the two infinite parallel plates in a porous medium channel,” Mathematical Problems in Engineering, vol. 2011, Article ID 349803, p. 10, 2011. View at: Publisher Site  Google Scholar
 J. D. Sherwood, “Squeeze flow of a powerlaw fluid between nonparallel plates,” Journal of NonNewtonian Fluid Mechanics, vol. 166, no. 56, pp. 289–296, 2011. View at: Publisher Site  Google Scholar
 S. U. S. Choi and J. A. Eastman, “Enhancing thermal conductivity of fluids with nanoparticles,” in Proceedings of the ASME International Mechanical Engineering Congress and Exposition, pp. 99–105, ASME, San Francisco, USA, November 1995. View at: Google Scholar
 S. K. Das, S. U. S. Choi, W. Yu, and T. Pradet, Nanofluids: Science and Technology, Wiley, Hoboken, NJ, USA, 2007.
 S. Kakaç and A. Pramuanjaroenkij, “Review of convective heat transfer enhancement with nanofluids,” International Journal of Heat and Mass Transfer, vol. 52, no. 1314, pp. 3187–3196, 2009. View at: Publisher Site  Google Scholar
 X.Q. Wang and A. S. Mujumdar, “A review on nanofluids—part II: experiments and applications,” Brazilian Journal of Chemical Engineering, vol. 25, no. 4, pp. 631–648, 2008. View at: Publisher Site  Google Scholar
 W. Daungthongsuk and S. Wongwises, “A critical review of convective heat transfer of nanofluids,” Renewable and Sustainable Energy Reviews, vol. 11, no. 5, pp. 797–817, 2007. View at: Publisher Site  Google Scholar
 A. Mokmeli and M. SaffarAvval, “Prediction of nanofluid convective heat transfer using the dispersion model,” International Journal of Thermal Sciences, vol. 49, no. 3, pp. 471–478, 2010. View at: Publisher Site  Google Scholar
 H. R. Rayatzadeh, M. SaffarAvval, M. Mansourkiaei, and A. Abbassi, “Effects of continuous sonication on laminar convective heat transfer inside a tube using waterTiO_{2} nanofluid,” Experimental Thermal and Fluid Science, vol. 48, pp. 8–14, 2013. View at: Publisher Site  Google Scholar
 K. Khanafer, K. Vafai, and M. Lightstone, “Buoyancydriven heat transfer enhancement in a twodimensional enclosure utilizing nanofluids,” International Journal of Heat and Mass Transfer, vol. 46, no. 19, pp. 3639–3653, 2003. View at: Publisher Site  Google Scholar
 S. S. Nourazar, M. Hatami, D. D. Ganji, and M. Khazayinejad, “Thermalflow boundary layer analysis of nanofluid over a porous stretching cylinder under the magnetic field effect,” Powder Technology, vol. 317, pp. 310–319, 2017. View at: Publisher Site  Google Scholar
 S. H. Seyedi, B. N. Saray, and M. R. H. Nobari, “Using interpolation scaling functions based on Galerkin method for solving nonNewtonian fluid flow between two vertical flat plates,” Applied Mathematics and Computation, vol. 269, pp. 488–496, 2015. View at: Publisher Site  Google Scholar
 S. H. Seyedi, B. Nemati Saray, and A. Ramazani, “On the multiscale simulation of squeezing nanofluid flow by a highprecision scheme,” Powder Technology, vol. 340, pp. 264–273, 2018. View at: Publisher Site  Google Scholar
 S. Hadi Seyedi, “Multiresolution simulation of flow in porous media using wavelet analysis,” Amirkabir University of Technology—Tehran Polytechnic, Tehran, Iran, 2013, Master thesis. View at: Google Scholar
 A. Dib, A. Haiahem, and B. Bousaid, “Approximate analytical solution of squeezing unsteady nanofluid flow,” Powder Technology, vol. 269, pp. 193–199, 2015. View at: Publisher Site  Google Scholar
 A. K. Gupta and S. Saha Ray, “Numerical treatment for investigation of squeezing unsteady nanofluid flow between two parallel plates,” Powder Technology, vol. 279, pp. 282–289, 2015. View at: Publisher Site  Google Scholar
 U. Khan, N. Ahmed, M. Asadullah, and S. Tauseef Mohyuddin, “Effects of viscous dissipation and slip velocity on twodimensional and axisymmetric squeezing flow of Cuwater and Cukerosene nanofluids,” Propulsion and Power Research, vol. 4, no. 1, pp. 40–49, 2015. View at: Publisher Site  Google Scholar
 Q. Z. Xue, “Model for thermal conductivity of carbon nanotubebased composites,” Physica B: Condensed Matter, vol. 368, no. 1–4, pp. 302–307, 2005. View at: Publisher Site  Google Scholar
 N. Freidoonimehr, B. Rostami, M. M. Rashidi, and E. Momoniat, “Analytical modelling of threedimensional squeezing nanofluid flow in a rotating channel on a lower stretching porous wall,” Mathematical Problems in Engineering, vol. 2014, Article ID 692728, 14 pages, 2014. View at: Publisher Site  Google Scholar
 H. Shahmohamadi and M. M. Rashidi, “VIM solution of squeezing MHD nanofluid flow in a rotating channel with lower stretching porous surface,” Advanced Powder Technology, vol. 27, no. 1, pp. 171–178, 2016. View at: Publisher Site  Google Scholar
 U. Khan, N. Ahmed, and S. T. MohyudDin, “Numerical investigation for three dimensional squeezing flow of nanofluid in a rotating channel with lower stretching wall suspended by carbon nanotubes,” Applied Thermal Engineering, vol. 113, pp. 1107–1117, 2017. View at: Publisher Site  Google Scholar
 S. S. Ghadikolaei, K. Hosseinzadeh, and D. D. Ganji, “Investigation on three dimensional squeezing flow of mixture base fluid (ethylene glycolwater) suspended by hybrid nanoparticle (Fe_{3}O_{4}Ag) dependent on shape factor,” Journal of Molecular Liquids, vol. 262, pp. 376–388, 2018. View at: Publisher Site  Google Scholar
 S. S. Ghadikolaei, K. Hosseinzadeh, M. Hatami, D. D. Ganji, and M. Armin, “Investigation for squeezing flow of ethylene glycol (C_{2}H_{6}O_{2}) carbon nanotubes (CNTs) in rotating stretching channel with nonlinear thermal radiation,” Journal of Molecular Liquids, vol. 263, pp. 10–21, 2018. View at: Publisher Site  Google Scholar
 H. C. Brinkman, “The viscosity of concentrated suspensions and solutions,” The Journal of Chemical Physics, vol. 20, no. 4, p. 571, 1952. View at: Publisher Site  Google Scholar
 B. Alpert, G. Beylkin, D. Gines, and L. Vozovoi, “Adaptive solution of partial differential equations in multiwavelet bases,” Journal of Computational Physics, vol. 182, no. 1, pp. 149–190, 2002. View at: Publisher Site  Google Scholar
 B. Alpert, G. Beylkin, R. Coifman, and V. Rokhlin, “Waveletlike bases for the fast solution of secondkind integral equations,” SIAM Journal on Scientific Computing, vol. 14, no. 1, pp. 159–184, 1993. View at: Publisher Site  Google Scholar
 B. N. Saray and J. Manafian, “Sparse representation of delay differential equation of pantograph type using multiwavelets galerkin method,” Engineering Computations, vol. 35, no. 2, pp. 887–903, 2018. View at: Google Scholar
 B. N. saray, “An efficient algorithm for solving Volterra integrodifferential equations based on Alpert’s multiwavelets Galerkin method,” Journal of Computational and Applied Mathematics, vol. 348, pp. 453–465, 2019. View at: Publisher Site  Google Scholar
 B. N. Saray, M. Lakestani, and C. Cattani, “Evaluation of mixed Crank–Nicolson scheme and Tau method for the solution of Klein–Gordon equation,” Applied Mathematics and Computation, vol. 331, pp. 169–181, 2018. View at: Publisher Site  Google Scholar
 J. S. Geronimo and F. Marcellán, “On Alpert multiwavelets,” Proceedings of the American Mathematical Society, vol. 143, no. 6, pp. 2479–2494, 2015. View at: Publisher Site  Google Scholar
 A. S. Butt and A. Ali, “Analysis of entropy generation effects in unsteady squeezing flow in a rotating channel with lower stretching permeable wall,” Journal of the Taiwan Institute of Chemical Engineers, vol. 48, pp. 8–17, 2015. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2019 S. Hadi Seyedi et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.