Mathematical Problems in Engineering
Volume 2019 (2019), Article ID 9890626, 18 pages
https://doi.org/10.1155/2019/9890626

Research Article

High-Accuracy Multiscale Simulation of Three-Dimensional Squeezing Carbon Nanotube-Based Flow inside a Rotating Stretching Channel

S. Hadi Seyedi,¹ Behzad Nemati Saray,² and Ali Ramazani³

¹Department of Mechanical Engineering, Wayne State University, 5050 Anthony Wayne Drive, Detroit, MI 48202, USA
²Department of Mathematics, Institute for Advanced Studies in Basic Sciences (IASBS), Zanjan 45137-66731, Iran
³Department of Mechanical Engineering, Massachusetts Institute of Technology, 77 Massachusetts Ave, Cambridge, MA 02139, USA

Correspondence should be addressed to Ali Ramazani; ramazani@mit.edu

Received 30 April 2019; Accepted 24 July 2019; Published 26 August 2019

Academic Editor: Xiao-Qiao He

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.

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 three-dimensional 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; single-wall 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 two-dimensional 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 three-dimensional 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 Galerkin-based 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 three-dimensional CNT-based 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 no-slip condition occurs between them. Single-wall carbon nanotubes (SWCNTs) and the multiwall carbon nanotubes (MWCNTs) are considered in this study.

Figure 1: Geometry of the problem: nanofluid in a rotating stretching channel.

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.

Table 1: Thermophysical properties of base fluid and carbon nanotubes [26].

Physical properties
Base fluid
SWCNT
MWCNT

4179
425
796

997.1
2600
1600

0.613
6600
3000

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 are

where 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 that

where 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 via

where 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 high-level spaces , . By this introduction, we are able to write

where 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 single-scale 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 two-scale 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 write

where 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, fourth-order 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 second-scale 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.

Table 2: Comparison of numerical values for and between the different methods taking and .

[36]

A
S

0.5
0.5
0.0
−7.419696
−0.237555
−7.419696
−0.237555

1.0
−4.445146
0.075782
−4.445146
0.075782

2.0
−1.297045
0.300313
−1.297045
0.300313

3.0
2.018006
0.469820
2.018006
0.469820

0.5
0.3
2.0
0.025724
0.372682
0.025724
0.372682

0.5
−1.297045
0.300313
−1.297045
0.300313

0.8
−3.369658
0.172293
−3.369658
0.172293

1.0
−4.814071
0.072293
−4.814071
0.072293

1.2
−6.311701
−0.040879
−6.311701
−0.040879

1.0
0.5
2
−1.291472
0.600474
−1.291472
0.600474

3.0
−1.231657
1.796681
−1.231657
1.796681

5.0
−1.110254
2.979537
−1.110254
2.979539

7.0
−0.924676
4.142754
−0.924676
4.142759

10.0
−0.522086
5.844396
−0.522085
5.844413

(a)

(b)

Figure 2: Effect of parameters J and r on the error.

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 multi-wall 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)

Figure 3: Effects of the nanoparticle volume fraction on the velocity and temperature profiles when  = 5,  = 0.1,  = 6.2, A = S = 0.5, and R = 0, taking and .

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 x-directional 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)

Figure 4: Effects of the rotation parameter, on the velocity and temperature profiles when ϕ =  = 0.1,  = 6.2, A = S = 0.5, and R = 0, taking and .

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 y-directional 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 x-directional velocity profile, there is a reverse relationship between the A and values. Moreover, a decrease in the x-directional 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)

Figure 5: Effects of the suction parameter on the velocity and temperature profiles when ϕ =  = 0.1,  = 6.2,  = 5, S = 0.5, and R = 0, taking and .

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 y-directional velocity for values by changing the squeezing number; however, for the x-directional 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)

Figure 6: Effects of the squeezing number (S) on the velocity and temperature profiles when ϕ =  = 0.1,  = 6.2,  = 5, A = 0.5, and R = 0, taking and .

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)

Figure 7: Effects of the R value on the temperature profile and Nusselt number when ϕ =  = 0.1,  = 6.2,  = 5, A = 0.5, and S = 0.5, taking and .

(a)

(b)

(c)

(d)

Figure 8: Effects of the squeezing number and concentration on the Nusselt number and skin friction coefficient when  = 0.1,  = 6.2,  = 5, and A = 0.5, taking and .

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)

Figure 9: Effects of the Eckert number on the Nusselt number in the presence of S and R values when ϕ =  = 0.1,  = 6.2,  = 5, and A = 0.5, taking and .

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)

Figure 10: Effects of the modified Eckert number on the Nusselt number in the presence of S and R values when ϕ =  = 0.1,  = 6.2,  = 5, and A = 0.5, taking and .

6. Conclusion

Three-dimensional CNT-based 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 Galerkin-based 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 Non-Newtonian Fluid Mechanics, vol. 132, no. 1–3, pp. 1–27, 2005.

C.-Y. Wang, “The squeezing of a fluid between two plates,” Journal of Applied Mechanics, vol. 43, no. 4, pp. 579–583, 1976.

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.

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.

P. Singh, V. Radhakrishnan, and K. A. Narayan, “Squeezing flow between parallel plates,” Ingenieur-Archiv, vol. 60, no. 4, pp. 274–281, 1990.

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.

J. D. Sherwood, “Squeeze flow of a power-law fluid between non-parallel plates,” Journal of Non-Newtonian Fluid Mechanics, vol. 166, no. 5-6, pp. 289–296, 2011.

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.

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. 13-14, pp. 3187–3196, 2009.

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.

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.

A. Mokmeli and M. Saffar-Avval, “Prediction of nanofluid convective heat transfer using the dispersion model,” International Journal of Thermal Sciences, vol. 49, no. 3, pp. 471–478, 2010.

H. R. Rayatzadeh, M. Saffar-Avval, M. Mansourkiaei, and A. Abbassi, “Effects of continuous sonication on laminar convective heat transfer inside a tube using water-TiO₂ nanofluid,” Experimental Thermal and Fluid Science, vol. 48, pp. 8–14, 2013.

K. Khanafer, K. Vafai, and M. Lightstone, “Buoyancy-driven heat transfer enhancement in a two-dimensional enclosure utilizing nanofluids,” International Journal of Heat and Mass Transfer, vol. 46, no. 19, pp. 3639–3653, 2003.

S. S. Nourazar, M. Hatami, D. D. Ganji, and M. Khazayinejad, “Thermal-flow boundary layer analysis of nanofluid over a porous stretching cylinder under the magnetic field effect,” Powder Technology, vol. 317, pp. 310–319, 2017.

S. H. Seyedi, B. N. Saray, and M. R. H. Nobari, “Using interpolation scaling functions based on Galerkin method for solving non-Newtonian fluid flow between two vertical flat plates,” Applied Mathematics and Computation, vol. 269, pp. 488–496, 2015.

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.

S. Hadi Seyedi, “Multi-resolution simulation of flow in porous media using wavelet analysis,” Amirkabir University of Technology—Tehran Polytechnic, Tehran, Iran, 2013, Master thesis.

A. Dib, A. Haiahem, and B. Bou-said, “Approximate analytical solution of squeezing unsteady nanofluid flow,” Powder Technology, vol. 269, pp. 193–199, 2015.

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.

U. Khan, N. Ahmed, M. Asadullah, and S. Tauseef Mohyud-din, “Effects of viscous dissipation and slip velocity on two-dimensional and axisymmetric squeezing flow of Cu-water and Cu-kerosene nanofluids,” Propulsion and Power Research, vol. 4, no. 1, pp. 40–49, 2015.

Q. Z. Xue, “Model for thermal conductivity of carbon nanotube-based composites,” Physica B: Condensed Matter, vol. 368, no. 1–4, pp. 302–307, 2005.

N. Freidoonimehr, B. Rostami, M. M. Rashidi, and E. Momoniat, “Analytical modelling of three-dimensional 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.

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.

U. Khan, N. Ahmed, and S. T. Mohyud-Din, “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.

S. S. Ghadikolaei, K. Hosseinzadeh, and D. D. Ganji, “Investigation on three dimensional squeezing flow of mixture base fluid (ethylene glycol-water) suspended by hybrid nanoparticle (Fe₃O₄-Ag) dependent on shape factor,” Journal of Molecular Liquids, vol. 262, pp. 376–388, 2018.

S. S. Ghadikolaei, K. Hosseinzadeh, M. Hatami, D. D. Ganji, and M. Armin, “Investigation for squeezing flow of ethylene glycol (C₂H₆O₂) carbon nanotubes (CNTs) in rotating stretching channel with nonlinear thermal radiation,” Journal of Molecular Liquids, vol. 263, pp. 10–21, 2018.

H. C. Brinkman, “The viscosity of concentrated suspensions and solutions,” The Journal of Chemical Physics, vol. 20, no. 4, p. 571, 1952.

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.

B. Alpert, G. Beylkin, R. Coifman, and V. Rokhlin, “Wavelet-like bases for the fast solution of second-kind integral equations,” SIAM Journal on Scientific Computing, vol. 14, no. 1, pp. 159–184, 1993.

B. N. Saray and J. Manafian, “Sparse representation of delay differential equation of pantograph type using multi-wavelets galerkin method,” Engineering Computations, vol. 35, no. 2, pp. 887–903, 2018.

B. N. saray, “An efficient algorithm for solving Volterra integro-differential equations based on Alpert’s multi-wavelets Galerkin method,” Journal of Computational and Applied Mathematics, vol. 348, pp. 453–465, 2019.

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.

J. S. Geronimo and F. Marcellán, “On Alpert multiwavelets,” Proceedings of the American Mathematical Society, vol. 143, no. 6, pp. 2479–2494, 2015.

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.