Mathematical Problems in Engineering

Volume 2013 (2013), Article ID 816262, 9 pages

http://dx.doi.org/10.1155/2013/816262

Research Article

Approximation of First Grade MHD Squeezing Fluid Flow with Slip Boundary Condition Using DTM and OHAM

Inayat Ullah, Hamid Khan, and M. T. Rahim

Department of Mathematics, National University of Computer and Emerging Sciences, FAST Peshawar Campus, Peshawar 25000, Pakistan

Received 27 December 2012; Accepted 24 April 2013

Academic Editor: Safa Bozkurt Coskun

Copyright © 2013 Inayat Ullah 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.

A steady two-dimensional axisymmetric flow of an incompressible viscous fluid under the influence of a uniform transverse magnetic field with slip boundary condition is studied. An ordinary nonlinear differential equation is formed by transforming the Navier-Stokes equations using the transformation . Differential transform and optimal homotopy analysis methods have been used to obtain the solutions by varying pertinent flow parameters. By using residuals in each case, the validity of solutions is established. Excellent results are obtained by using the proposed schemes. The influence of different parameters on the flow is shown through graphs.

1. Introduction

Squeezing flows induce when normal stresses or vertical velocities are applied externally by means of a mobile boundary [1]. Stefan [2] carried out the fundamental research in this field. Engineers studied the analysis of Newtonian fluid squeezed between two infinite planar plates [3, 4]. There are many applications of squeezing flow in food industry, especially in chemical engineering [4]. Polymer processing, compression, and injection modeling are the practical examples of squeezing flow. Thin Newtonian liquid films squeezing between two plates were studied by Grimm [5]. Squeezing flow under the influence of magnetic field is widely applied to bearing with liquid-metal lubrication [6–9]. Islam et al. [10] studied squeezing fluid flow between the two infinite parallel plates in a porous medium channel.

In case of many polymeric liquids when the weight of molecule is high, then they show slip at the boundary. The no-slip boundary condition is not applicable in this case. In many cases such as thin film problems, rarefied fluid problems, fluids containing concentrated suspensions, and flow on multiple interfaces, the no-slip boundary condition fails to work. Navier [11], for the first time, proposed the general boundary condition which demonstrates the fluid slip at the surface. The difference of fluid velocity and velocity of the boundary is proportional to the shear stress at that boundary. The proportionality constant is named the slip parameter having length as its dimension.

The slip condition is of great importance especially when fluids with elastic character are under consideration [12]. Newtonian fluid was considered by Ebaid [13] to study the effects of magnetic field and wall slip conditions on the peristaltic transport in an asymmetric channel. It has great importance in medical sciences, particularly in polishing artificial heart valves and internal cavities in many manufactured parts achieved by embedding such fluids with abrasives [14]. The influence of slip on the peristaltic motion of third-order fluid in asymmetric channel is studied by Hayat et al. [15]. The effects of slip condition on the rotating flow of a third grade fluid in a nonporous medium are investigated by Hayat and Abelman [16]. The extension of the work in [16] to a porous medium and obtaining the numerical solutions for the steady magnetohydrodynamics flow of a third grade fluid in a rotating frame is presented by Abelman et al. [17].

Various perturbation techniques can be applied to solve nonlinear differential equations [18]. The limitations of these methods are based on the assumption of small parameter and there is no proper way to select this parameter. The idea of homotopy was combined with perturbation in the last decade. Liao and He did the fundamental work on it. Homotopy analysis method (HAM) was proposed by Liao [19–21]. Homotopy perturbation method was introduced by He [22–24]. Marinca et al. [25, 26] introduced OHAM for approximate solution of nonlinear problems of thin film flow of a fourth-grade fluid down a vertical cylinder and for the study of the behavior of nonlinear mechanical vibration of electrical machines. It is observed that HPM and HAM are the special cases of OHAM [27, 28].

Differential transform method (DTM) was initially introduced by Zhou in 1986 [29]. Çatal [30–33] studied solution of free vibration equations of beam on elastic soil, buckling analysis of partially embedded pile in elastic soil, analysis of free vibration of beam on elastic soil, and response of forced Euler-Bernoulli beams by using DTM. Ayaz [34] applied this method to differential algebraic equations (DAEs) of index-1. Liu and Song [35] investigated that DTM is effective in case of index-2 DAEs but not in index-3. Ayaz [36] studied the applications of two-dimensional DTM in case of partial differential equations.

The objective of this research paper is to use OHAM and DTM for studying MHD squeezing flow with slip boundary condition between two infinite plates approaching each other slowly. In Section 2, the basic equations are derived and the model is reduced to a nonlinear boundary value problem. In Sections 3 and 4, the basic ideas of DTM and its use for our problem are presented while the idea of OHAM and its application to our problem are in Sections 5 and 6, respectively.

2. Problem Formulation

Let us consider, in the presence of a magnetic field, a squeezing flow of an incompressible Newtonian fluid with constant density and viscosity , squeezed between two large planar parallel plates separated by a small distance approaching each other with a low constant velocity (Figure 1). Assume that the flow is quasisteady [1, 3], and the Navier-Stokes equations governing such flow when inertial terms are retained are

where denotes the material time derivative, is the Cauchy stress tensor given by with . is the total magnetic field given by . and represent the imposed and induced magnetic fields, respectively. The modified Ohm’s law and Maxwell’s equations (see [37] and the references therein), in the absence of displacement currents, are

Here is the electric current density, represents the electrical conductivity, the electric field, and the magnetic permeability. If , , and are constant, is negligible as compared to , is perpendicular to so that the Reynold number is small with no electric field in the fluid flow region and then the magneto hydrodynamic force involved can be written as

Assuming that the plates are nonconducting and the magnetic field is applied along the -axis. The gap distance between the plates changes slowly with time for small values of the velocity so that it can be taken as constant. An axisymmetric flow in cylindrical coordinates with -axis perpendicular to plates and at the plates. For axial symmetry, is represented by . In view of negligible body forces with no tangential velocity, Navier-Stokes equations [1, 6, 10] in cylindrical coordinates are

where

Introducing the stream function , we have,

Eliminating from (4), we have

where

Using the transformation , (7) can be written as

subject to the slip boundary conditions

The nondimensional parameters are , , and . Omitting the *, (9) and (10) become

with and , are the Reynolds and Hartmann numbers, respectively.

Figure 1: A steady squeezing axisymmetric fluid flow between two parallel plates.

3. Differential Transform Method

One-dimensional differential transform of a function is defined as follows [29–35]:

is also called the -function. The inverse transform of is is defined as

Combining (12) and (13), we can write

Usually we express in a finite series, that is,

which implies that

is negligibly small.

Some fundamental theorems on one-dimensional differential transform are as follows.

Theorem 1. If , then .

Theorem 2. If , then .

Theorem 3. If , then .

Theorem 4. If , then

4. Analysis of Differential Transform Method

For fourth-order boundary value problem [10],

with boundary conditions

The differential transform is given by

with transformed boundary conditions

Using (20) and (21), for is evaluated to obtain the solution up to as follows:

5. Application of DTM to Our Problem

Here, we consider an incompressible Newtonian fluid, squeezed between two large planar, parallel plates which are separated by a small distance and moving towards each other with velocity . For small values of , the gap distance between the plates changes slowly with time , so that it can be taken as constant, and the flow is steady as

with slip boundary conditions

The differential transform of (23) is given by

with transformed boundary conditions

where and are unknowns to be determined later [38]. Using (25) and (26), the values of , are

As , have been nondimensionalized, so let us fix them one to find the values of and from the following system:

which leads us to the following values:

Using (22), the approximate solution of the BVP is as follows:

This is a series solution and we note that the terms approach to zero if we go forward which is the beauty of DTM and other series solutions.

6. Basic Idea of OHAM

Let us consider a boundary value problem [25–27]

along with boundary conditions

where denotes the independent variable, , are unknown and known functions and , , are linear, nonlinear, and boundary operators, respectively.

According to OHAM a homotopy, satisfying

for and , is constructed with nonzero auxiliary function for with and for an unknown function . For the equalities and hold. So it is concluded that the solution approaches from to as varies from to . For ,

We choose the auxiliary function such that

where are the convergence controlling constants to be determined. Expanding in Taylor’s series about to get the approximation

substituting (36) into (33), and equating the coefficients of like powers of , we get the following.

Zeroth-order problem is

First-order problem is

Second-order problem is

Generally,

with boundary conditions

where is the coefficient of in the expansion of about as

The convergence of the series (36) depends upon . For convergence at , the th-order approximation is

Substituting (43) into (31), the expression for the residual is

will be the exact solution if , but in case of nonlinearity, it does not happen generally. To find the optimal values of the constants , there are different methods to apply. We apply the method of least square as follows:

Minimizing this function, we have

where and are chosen in the domain of the problem for locating the suitable . The approximate solution is well determined for these values of .

7. Application of OHAM to Our Problem

Using , , and in (11), we have the following.

Zeroth-order problem is

First-order problem is

Second-order problem is

Third-order problem is

Fourth-order problem is

Fifth-order problem is

By considering the fifth-order solution, we have

The residual of the problem is

In order to find , , we apply the method of least square as follows:

Solving (56) for , , we get

Using these values of , , the approximate solution is given by:

As , the components of vanish when going onward. In Table 1, the solutions obtained by OHAM and DTM along with their residual are presented. Table 2 is constructed to discuss residuals for different values of .

Table 1: DTM and OHAM solutions along with the residuals.

DTM solution
DTM residual
OHAM solution
OHAM residual

0. 0
0.0
0.0
0.0
0.0

0.1
0.075739

0.075739

0.2
0.152935

0.152935

0.3
0.233046

0.233046

0.4
0.317540

0.317540

0.5
0.407893

0.407893

0.6
0.505591

0.505591

0.7
0.612134

0.612134

0.8
0.729034

0.729034

0.9
0.857813

0.857813

1. 0
1.0

1.0

Table 2: Residuals for different values of γ using DTM.

0.0
0.0
0.0
0.0
0.0

0.0

8. Numerical Results and Discussion

Here we analyzed, under a transverse magnetic field, the squeezing fluid flow between two parallel infinite plates. DTM and OHAM are applied to obtain the solution of reduced nonlinear boundary value problem. The convergence region can easily be adjusted using OHAM as compared to DTM in this case. Differential transform method which is the straight forward application of Taylor’s series suffers in general and particularly in initial value problems. Beside this fact, DTM leads easily to closed form solution. OHAM does not require discretization like other numerical and approximate methods. This technique requires less computational effort and the convergence of this method to the exact is fast enough. For computational purposes, we use Mathematica software.

In order to get an insight into the problem, it is important to discuss the influence of different flow parameters on the velocity profile of the fluid.

Figure 2 represents the influence of Reynold’s number on the velocity of the fluid. It is observed from Figure 2 that by increasing the velocity of the fluid decreases. The influence of Hartmann number on the velocity of the fluid is shown in Figure 3. It is concluded that the velocity increases by increasing in case of . However, from Figure 4, the velocity decreases due to an increase in for . This observation is in agreement with the literature (see [39, Figure 5]). Figure 5 shows the influence of the slip parameter on the fluid velocity. By increasing , it is observed that the velocity of the fluid increases.

Figure 2: Plot shows the Graph of for different values of .

Figure 3: Analysis of for different values of when .

Figure 4: Analysis of by varying when .

Figure 5: Analysis of by Varying .

References

T. C. Papanastasiou, G. C. Georgiou, and A. N. Alexandrou, Viscous Fluid Flow, CRC Press, 1994.

M. J. Stefan, “Versuch Uber die scheinbare adhesion,” Akademie Der Wissenschaften in Wien, vol. 69, article 713, 1874.

Q. K. Ghori, M. Ahmed, and A. M. Siddiqui, “Application of homotopy perturbation method to squeezing flow of a newtonian fluid,” International Journal of Nonlinear Sciences and Numerical Simulation, vol. 8, no. 2, pp. 179–184, 2007.

X. J. Ran, Q. Y. Zhu, and Y. Li, “An explicit series solution of the squeezing flow between two infinite parallel plates,” Communications in Nonlinear Science and Numerical Simulation, vol. 14, no. 1, pp. 119–132, 2007.

R. J. Grimm, “Squeezing flows of Newtonian liquid films an analysis including fluid inertia,” Applied Scientific Research, vol. 32, no. 2, pp. 149–166, 1976.

W. F. Hughes and R. A. Elco, “Magnetohydrodynamic lubrication flow between parallel rotating disks,” Journal of Fluid Mechanics, vol. 13, pp. 21–32, 1962.

S. Kamiyama, “Inertia Effects in MHD hydrostatic thrust bearing,” Transactions ASME, vol. 91, pp. 589–596, 1969.

E. A. Hamza, “Magnetohydrodynamic squeeze film,” Journal of Tribology, vol. 110, no. 2, pp. 375–377, 1988.

S. Bhattacharyya and A. Pal, “Unsteady MHD squeezing flow between two parallel rotating discs,” Mechanics Research Communications, vol. 24, no. 6, pp. 615–623, 1997.

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, 10 pages, 2011.

C.-L.-M.-H. Navier, “Sur les lois de l’équilibre et du mouvement des corps solides élastiques,” Bulletin des Sciences par la Société Philomatique de Paris, pp. 177–181, 1823.

C. le Roux, “Existence and uniqueness of the flow of second-grade fluids with slip boundary conditions,” Archive for Rational Mechanics and Analysis, vol. 148, no. 4, pp. 309–356, 1999.

A. Ebaid, “Effects of magnetic field and wall slip conditions on the peristaltic transport of a Newtonian fluid in an asymmetric channel,” Physics Letters A, vol. 372, no. 24, pp. 4493–4499, 2008.

L. J. Rhooades, R. Resnic, T. O' Bradovich, and S. Stegman, “Abrasive flow machining of cylinder heads and its positive effects on performance and cost characteristics,” Tech. Rep., Dearborn, Mich, USA, 1996.

T. Hayat, M. U. Qureshi, and N. Ali, “The inuence of slip on the peristaltic motion of third order fluid in an asymmetric channel,” Physics Letters A, vol. 372, pp. 2653–2664, 2008.

T. Hayat and S. Abelman, “A numerical study of the influence of slip boundary condition on rotating flow,” International Journal of Computational Fluid Dynamics, vol. 21, no. 1, pp. 21–27, 2007.

S. Abelman, E. Momoniat, and T. Hayat, “Steady MHD flow of a third grade fluid in a rotating frame and porous space,” Nonlinear Analysis: Real World Applications, vol. 10, no. 6, pp. 3322–3328, 2009.

A. H. Nayfeh, Problems in Perturbation, A Wiley-Interscience Publication, John Wiley & Sons, New York, NY, USA, 1985.

S. J. Liao, The proposed homotopy analysis technique for the solution of nonlinear problems [Ph.D. thesis], Shanghai Jiao Tong University, 1992.

S. Liao, Beyond Perturbation. Introduction to the Homotopy Analysis Method, vol. 2 of CRC Series: Modern Mechanics and Mathematics, Chapman & Hall/CRC, Boca Raton, Fla, USA, 2004.

S. Liao, “On the homotopy analysis method for nonlinear problems,” Applied Mathematics and Computation, vol. 147, no. 2, pp. 499–513, 2004.

J.-H. He, “Homotopy perturbation method for solving boundary value problems,” Physics Letters A, vol. 350, no. 1-2, pp. 87–88, 2006.

J. H. He, “Application of homotopy perturbation method to nonlinear wave equations,” Chaos, Solitons and Fractals, vol. 26, no. 3, pp. 695–700, 2005.

J.-H. He, “Comparison of homotopy perturbation method and homotopy analysis method,” Applied Mathematics and Computation, vol. 156, no. 2, pp. 527–539, 2004.

V. Marinca, N. Herişanu, and I. Nemeş, “Optimal homotopy asymptotic method with application to thin film flow,” Central European Journal of Physics, vol. 6, no. 3, pp. 648–653, 2008.

V. Marinca and N. Herişanu, “Application of optimal homotopy asymptotic Method for solving nonlinear equations arising in heat transfer,” International Communications in Heat and Mass Transfer, vol. 35, no. 6, pp. 710–715, 2008.

J. Ali, S. Islam, S. Islam, and G. Zaman, “The solution of multipoint boundary value problems by the optimal homotopy asymptotic method,” Computers & Mathematics with Applications, vol. 59, no. 6, pp. 2000–2006, 2010.

V. Marinca and N. Herisanu, “Optimal homotopy perturbation method for strongly nonlinear differential equations,” Nonlinear Science Letters A, vol. 1, no. 3, pp. 273–280, 2010.

J. K. Zhou, Deferential Transformation and Its Application for Electrical Circuits, Huazhong University Press, Wuhan, China, 1986.

S. Çatal, “Solution of free vibration equations of beam on elastic soil by using differential transform method,” Applied Mathematical Modelling, vol. 32, no. 9, pp. 1744–1757, 2008.

S. Çatal and H. H. Çatal, “Buckling analysis of partially embedded pile in elastic soil using differential transform method,” Structural Engineering and Mechanics, vol. 24, no. 2, pp. 246–269, 2006.

S. Çatal, “Analysis of free vibration of beam on elastic soil using differential transform method,” Structural Engineering and Mechanics, vol. 24, no. 1, pp. 51–63, 2006.

S. Catal, “Response of forced Euler-Bernoulli beams using differential transform method,” Structural Engineering and Mechanics, vol. 42, no. 1, pp. 95–119, 2012.

F. Ayaz, “Applications of differential transform method to differential-algebraic equations,” Applied Mathematics and Computation, vol. 152, no. 3, pp. 649–657, 2004.

H. Liu and Y. Song, “Differential transform method applied to high index differential-algebraic equations,” Applied Mathematics and Computation, vol. 184, no. 2, pp. 748–753, 2007.

F. Ayaz, “On the two-dimensional differential transform method,” Applied Mathematics and Computation, vol. 143, no. 2-3, pp. 361–374, 2003.

M. R. Mohyuddin and T. Götz, “Resonance behaviour of viscoelastic fluid in Poiseuille flow in the presence of a transversal magnetic field,” International Journal for Numerical Methods in Fluids, vol. 49, no. 8, pp. 837–847, 2005.

S. H. Chang and I. L. Chang, “A new algorithm for calculating one-dimensional differ- ential transform of nonlinear functions,” Applied Mathematics and Computation, vol. 195, no. 2, pp. 799–808, 2008.

S. Islam, M. Ullah, G. Zaman, and M. Idrees, “Approximate solutions to MHD squeezing fluid flow,” Journal of Applied Mathematics & Informatics, vol. 29, no. 5-6, pp. 1081–1096, 2011.