Journal of Applied Mathematics

Volume 2015, Article ID 864190, 5 pages

http://dx.doi.org/10.1155/2015/864190

## On the Study of Oscillating Viscous Flows by Using the Adomian-Padé Approximation

^{1}Division of Mathematics, General Education Center, Chienkuo Technology University, Changhua City 500, Taiwan^{2}International Wave Dynamics Research Center, National Cheng Kung University, Tainan 701, Taiwan

Received 8 November 2014; Accepted 2 April 2015

Academic Editor: Charalampos Tsitouras

Copyright © 2015 Chi-Min Liu. 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

The Adomian-Padé technique is applied to examine two oscillating viscous flows, the Stokes’ second problem and the pressure-driven pulsating flow. Main purposes for studying oscillating flows are not only to verify the accuracy of the approximation solution, but also to provide a basis for analyzing more problems by the present method with the help of Fourier analysis. Results show that the Adomian-Padé approximation presents a very excellent behavior in comparison with the exact solution of Stokes’ second problem. For the pulsating flow, only the Adomian decomposition method is required to perform the calculation as the fluid domain is finite where the Padé approximant may not provide a better solution. Based on present results, more problems can be mathematically solved by using the Adomian-Padé technique, the Fourier analysis, and powerful computers.

#### 1. Introduction

The Adomian decomposition method has been widely studied and applied to solve mathematical problems [1]. It provides an efficient way to study a rather wide class of nonlinear as well as stochastic equations without linearization, perturbation, closure approximation, or discretization methods [2]. The main idea of the method is to decompose the key variable, say, , into an infinite series of () without requirements of weak nonlinearity and small perturbation. The first decomposed term is determined by the given initial or boundary conditions. Then other decomposed terms can be calculated by applying the integral operators to the governing equation with the help of lower-order terms. In recent years, there has been a great amount of efforts on applying the Adomian decomposition method have been applied to many fields [3–8] and some improvements of algorithms were also presented [9–12].

As all decomposed terms except the first one are calculated by integration, the solution is usually expressed in a form of a polynomial of the variable of integration (the temporal or spatial variable). Therefore, the applicable range of the Adomian decomposed solution will diverge quickly while this integration variable grows. To overcome the weakness, the Padé approximation is adopted to improve the accuracy of the solution. The main idea of Padé approximant is to transfer the original polynomial into a rational function of the order approximant [13]. In comparison of the form of Taylor series, Padé approximant usually gives a larger valid range of the derived solution. The Padé approximation has been widely adopted to improve the accuracy of the Adomian decomposed solution [14–16].

In this paper two cases of oscillating viscous flows will be investigated by using the Adomian-Padé approximation. The main reason why we study oscillating flows is that the present study can provide a basis for analyzing more problems while applying the Fourier analysis. The organization of this paper is as follows. In Section 2, Stokes’ second problem describing the flow induced by an oscillating plate below is examined. Boundary conditions required for calculating the velocity evolution of the whole domain are the velocity and its gradient at the plate. The derived solution is compared with the exact solution for verifying its validity. In Section 3 the pressure-driven pulsating flow between two infinite parallel plates is calculated by giving the oscillating pressure gradient along the flow direction. The velocity at the center plane is given as the beginning condition. Results in Sections 2 and 3 show very excellent behaviors in comparison with their exact solutions. Conclusion remarks are made in Section 4.

#### 2. Stokes’ Second Problem

A viscous flow generated by an oscillating plate below is the well-known Stokes’ second problem. The plate is located at the plane and oscillates with the velocity in its own plane. The kinematic viscosity of the fluid in the domain is . The corresponding dimensionless variables are assigned bywhere is the velocity, the spatial parameter, and the time. Hence the problem can be described with dimensionless governing equationand boundary conditionsThe exact solution for (1) to (3) is (see [17, 18] for details)It is noted that if the initial condition, , is considered in this problem, the corresponding solution will possess both steady-state and transient parts shown as [17]where denotes the complementary error function with the complex argument.

Now the Adomian-Padé method is applied to solve the problem. First we define the operatorHence (2) can be represented asApplying the inverse operatorto (7) yieldsThe first two terms in the right-hand side of (9) contain two boundary conditions at the plate (). They play an important role while applying the beginning conditions of the Adomian approximation. They are given by (4) and its derivativeNote that another boundary condition which is not used in the Adomian analysis will be regarded as a condition for examining the accuracy of the derived solution. Now the velocity is assumed to beSubstituting (11) into (9) gives the following relations:By calculating (13) term by term, we have The coefficients of Padé solution can be determined by balancing and The results areBy following similar processes shown above, higher-order Padé approximants can be readily obtained with the help of commercial software (e.g., Mathematica) and the details are neglected herein. Figure 1 shows velocity profiles calculated by the Adomian-Padé , , , and solutions and by the exact solution (4) for the case . Solid and dash lines indicate the Padé solutions and the exact solution, respectively. It is seen that the higher-order Adomian-Padé solution will result in a more accurate behavior in comparison with the exact solution. Due to rapid and great progress of calculating ability of computers, the applicable domain in will be efficiently expanded by much higher-order Padé solutions. Solutions for are plotted in Figure 2. Results are similar to those demonstrated for Figure 1.