On the Study of Oscillating Viscous Flows by Using the Adomian-Padé Approximation
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.
The Adomian decomposition method has been widely studied and applied to solve mathematical problems . 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 . 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 . 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 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.
The present method can play an important role in practical cases. For example, if the boundary conditions at the plate can be measured by current meters or other facilities, velocity profile at any elevation can be mathematically calculated without using the equipment to measure the flow speed at the desired place.
3. Pressure-Driven Pulsating Flow
Another oscillating flow, the pressure-driven pulsating flow, is studied in this section. This flow describes that a fluid bounded by two parallel plates located at is driven by a time-dependent pressure gradient. The governing equation iswhere is the pressure andis given. By introducing the following dimensionless variables:(17) can be rewritten aswhere and are assigned hereafter for the purpose of simplification. The boundary conditions at both plates areThe exact solution for (20) and (21) is (see  for details)where denotes the real part of the complex. The velocity and its gradient at the plane will be adopted as the beginning conditions. They arewhere the former is obtained from (22) and the latter indicates that the flow is symmetric with respect to . It is also noted that boundary conditions at both plates will be seen as a condition for verifying the accuracy of the derived solution. Applying (8) to (20) with the help of (23) leads toNow we substitute (11) into (24); the result isUsing defined in previous section and balancing it with , the coefficients of are solved:Figure 3 presents a comparison of different approximation solutions for the upper fluid domain (). Ratios of (solid lines) and (dash lines) to the exact solution are plotted for times (thin lines) and (bold lines). It is seen that the polynomial solution behaves slightly better than the Padé solution; namely, the former approximation is slightly closer to its exact solution when approaches unity. The reason is that, in the range of , the Padé solution may not give a more excellent simulation than the truncated polynomial solution . This implies that higher-order Padé approximant may be not required to improve the accuracy of the polynomial solution for problems considered in a finite domain.
This paper presents the examination of Stokes’ second problem and the pressure-driven pulsating flow by applying the Adomian-Padé technique. For Stokes’ second problem, higher-order Adomian-Padé solution behaves very well in comparison with the exact solution while the spatial parameter grows. For pressure-driven pulsating flow, the Adomian approximation provides satisfactory results and the application of the Padé approximant may be unnecessary for the case of the finite domain. The above results demonstrate that the method used in this paper can be applied to solve more complicated problems with the help of the Fourier analysis and powerful computers.
Conflict of Interests
The author declares that there is no conflict of interests regarding the publication of this paper.
Financial support from the Ministry of Science and Technology of Taiwan through Grants NSC 101-2221-E-270-001-MY2 and MOST 104-2911-I-006-301 is acknowledged.
A. H. Bokhari, G. Mohammad, M. T. Mustafa, and F. D. Zaman, “Adomian decomposition method for a nonlinear heat equation with temperature dependent thermal properties,” Mathematical Problems in Engineering, vol. 2009, Article ID 926086, 12 pages, 2009.View at: Publisher Site | Google Scholar | MathSciNet
A.-M. Wazwaz, R. Rach, and J.-S. Duan, “Adomian decomposition method for solving the Volterra integral form of the Lane-Emden equations with initial values and boundary conditions,” Applied Mathematics and Computation, vol. 219, no. 10, pp. 5004–5019, 2013.View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet
G. A. Baker and P. Graves-Morris, Padé Approximants, Cambridge University Press, 2nd edition, 1996.
I. G. Currie, Fundamental Mechanics of Fluids, McGraw-Hill, New York, NY, USA, 2nd edition, 1993.