Application of the Adomian Decomposition Method to Oscillating Viscous Flows
Chi-Min Liu
Division of Mathematics, General Education Center, Chienkuo Technology University, Changhua City, Taiwan
Email address:
To cite this article:
Chi-Min Liu. Application of the Adomian Decomposition Method to Oscillating Viscous Flows. Applied and Computational Mathematics. Vol. 5, No. 3, 2016, pp. 121-132. doi: 10.11648/j.acm.20160503.15
Received: June 3, 2016; Accepted: June 13, 2016; Published: June 29, 2016
Abstract: In this paper three oscillating viscous flows are studied by applying the Adomian decomposition method (ADM). Major improvement is on the choice of the assignment of the first term of the decomposition series. Different from past studies in which the initial velocity profile of the whole domain is assigned as the first term of the decomposition series, the assignment in present study is simply the boundary velocity for Stokes’ second problem and the pressure gradient for pulsatile flows. This improvement demonstrates and implies that ADM is not only good in approaching the known exact solution, but also possesses the practicability in treating realistic problems. The derived approximate solutions accurate up to any order can be obtained after two key parameters are determined. Present results show an excellent agreement with those calculated by the exact solutions. Based on the present results, more periodic problems can be analyzed by ADM with the help of Fourier analysis.
Keywords: Adomian Decomposition Method, Stokes’ Second Problem, Pulsatile Flow, Starting Assignment
1. Introduction
Adomian decomposition method (ADM) [1-2] is a powerful method to analyze mathematical and practical problems. Without the requirements of discretization, linearization and perturbation, ADM can be applied to solve linear or nonlinear, ordinary or partial differential equations and integral equations for either initial-value or boundary-value problems. A great deal of efforts has been devoted to the applications of ADM [3-11], the improvement associated with Padé approximation [12-14], and the examination of the convergent behavior [15-16].
Recently, Liu [17] applied ADM to investigate Stokes’ second problem and planar pulsatile flow in a semi-infinite domain. He assigned the first term of the decomposition series from the known exact solutions [18-19] of the velocity and the shear stress at the plate for Stokes’ second problem, and the pressure gradient for pulsatile flow. Using the differential operator with respect to the spatial coordinate, the derived solution is displayed in terms of powers in spatial variable. Similarly, Farkhadnia et al. [20] also applied ADM and other two methods to study Stokes’ second problem. The starting assignment of the first term of the decomposition series (Eq.(64)) is the velocity profile at which is directly adopted from the exact solution [21-22]. Different from Liu’s solution [17], the solution is shown in terms of powers in time due to the different definition of the differential operator they chose. Though Liu and Farkhadnia et al. successfully validated the efficiency and accuracy of ADM by assigning the first term of the decomposition series acquired from the exact solution, however, the ability of ADM for simulating realistic flows by simply assigning the given boundary velocity (for Stokes’ problem) or the pressure gradient (for pulsatile flow) as the first term of the series is not well examined.
Due to aforementioned demonstration, Stokes’ second problem in a finite-depth domain, pulsatile flows between plates and in a circular tube will be studied by applying ADM. Results for three cases are respectively shown in Sections 2 to 4. Each Section includes the derivation processes of the ADM solution and the examination of convergent behavior of two key parameters. The derived solution will be compared with the exact solution to verify the efficiency and accuracy of ADM. Conclusions are made in Section 5.
2. Stokes’ Second Problem in A Finite-Depth Domain
Stokes’ problem [23], also known as Rayleigh problem [24], describes that a viscous fluid which occupies a semi-infinite domain is driven by a moving plate below. The plate moving in a constant speed or in an oscillation type is named as the first or second problem. Analytical solutions to the flow velocity [18-19, 25] and the wall stress [19] were well investigated and discussed. Extended from the semi-infinite domain of fluid, Stokes’ problem for a fluid in a finite-depth domain was studied [26]. The flow considered is bounded by a free surface above and a moving plate below. Solutions of the first and second problems are derived and shown in the form of infinite series. For the second problem, the exact solution can be separated into the steady-state and the transient parts. In addition to analytical investigation, numerical methods based on discretization schemes, which include finite-difference method [27], finite-element method [28] and mesh-free numerical scheme [29], are performed to analyze the semi-infinite Stokes’ problems. In general numerical results show a good agreement with analytical solutions.
In this section, ADM is applied to solve Stokes’ second problem in a finite-depth domain by simply assigning the plate velocity to the first term of the decomposition series instead of those adopted in [17] and [20]. Consider a fluid of kinematic viscosity which is bounded by a plate located at and a free surface at . The fluid is driven by the oscillating plate with the speed . The momentum equation and boundary conditions for the velocity are
(1)
(2)
(3)
Dimensional parameters
, , , (4)
are defined and applied to Eqs.(1) to (3). The results are
(5)
(6)
(7)
The value of will be assigned to be unity henceforward without loss of generality. To apply ADM, a differential operator
(8)
is defined and then Eq.(5) can be shown as
(9)
Applying the inverse operation of Eq.(8)
(10)
to Eq.(9) results in
(11)
The solution of is assumed to be
(12)
By using the boundary condition Eq.(6), the first term is assigned as
(13)
where represents the velocity gradient at which is an unknown function of and will be determined later. The remaining terms at the right-hand side of Eq.(12) can be recursively calculated by
(14)
Now to are shown below for clear comprehension
(15)
where the superscript demotes the i-th differentiation with respect to , and the factorial. Therefore the complete solution can be written as
(16)
Inserting Eq.(16) into Eq.(7) gives
(17)
As the plate oscillates steadily, is assumed to be
(18)
where and are constant parameters to be determined. Substituting Eq.(18) into Eq.(17) results in
(19)
Consider the approximate solution collecting terms up to , we define
(20)
Now parameters and for can be solved by
(21)
For example, if is considered, parameters and can be determined by
(22)
and the corresponding solutions are
(23)
Based on above derivation, the approximate solution Eq.(20) accurate up to any order can be solved after parameters and are determined by Eq.(21).
Based on the derived solution, the convergent behavior of parameters and is firstly examined. The solution pairs of for to by solving Eq.(21) are listed in Table 1. In each column, the decimal is the approximate value and the fraction in the bracket is the exact value solved from Eq.(21). It is clearly found that solutions of and will converge when becomes large.
Next, the approximate solutions , , and the exact solution [26]
, ,(24)
are plotted in Figure 1 for comparison. Velocity distributions along at , , and are respectively displayed while the period of oscillation is . , and are respectively drawn in green, blue and red solid lines and is displayed in black dash line. It shows that provides the best behavior than while is the poorest one. The corresponding quantitative measurement of errors is displayed in Figure 2. Error for each approximate solution is defined by
, .(25)
In this figure , and are displayed in green, blue and red curves, respectively. It is found that is around 0.1 while is smaller than 0.01 for most time domain. is much smaller than and which indicates again the rapid convergence of the present ADM solution.
3. Pulsatile Flow Between Horizontal Plates
A flow with periodic variation is known as pulsatile flow. It is commonly observed in arterial blood circulation, hydraulic machines and motor engines. Exact solutions of pulsatile flows between plates and in a circular tube are well known and are referred to [30] and [31]. In this section, a pressure-driven pulsatile flow bounded by two fixed horizontal plates is studied by ADM. An oscillating pressure gradient is the only driven force for the flow system. Liu [17] studied the same problem by assigning the velocity and its derivative at the middle plane between two plates. For the sake of simulating a real flow, the pressure gradient will be the only assignment in present study.
Consider a fluid of density and kinematic viscosity bounded by two plates at . The momentum equation is
(26)
where the pressure gradient is given by
(27)
As the flow velocity at plates is zero, the corresponding boundary conditions are
(28)
Applying the following dimensionless variables
, , ,(29)
to Eqs.(26) and (28), the results are
(30)
where
(31)
In the following analysis will be assigned without loss of generality. The boundary conditions at both plates become
(32)
Now applying the operator defined in Eq.(8) to Eq.(30) leads to
(33)
Applying the inverse operator of Eq.(10) to above equation results in
(34)
As the flow is symmetrical to , it implies that has to be zero. The velocity at is assumed to be which will be determined later. Now the starting term of the decomposition series is given by
(35)
The remaining term can be recursively calculated by Eq.(14). It gives
(36)
The complete solution is shown as
(37)
Now the boundary condition at the upper plate, , is adopted to calculate . The result is
(38)
As the pressure gradient is oscillatory and so as the flow, is assumed to be of the form
(39)
where and are unknowns to be determined. Substituting Eq.(39) into Eq.(38) gives
(40)
For the approximate solution , parameters and can solved by
(41)
Solution pairs of and for the first four approximate solutions are shown in Table 2. It is also found that and rapidly converge when the order of grows. Figure 3 shows a comparison of , , and the exact solution [30]
, .(42)
As the parameters and converge very quickly, very slight difference exists between , and the exact solution. This is also observed in Figure 4 which displays the quantitative measurement of errors defined by Eq.(25).
(a)
(b)
(c)
(d)
4. Pulsatile Flow in a Circular Tube
The problem studied in this section is almost the same as that in previous section, except the geometry is changed from a planar type to a circular type. Consider a viscous flow in a circular tube of radius driven by a pressure gradient, the momentum equation reads
(43)
where the pressure gradient is given by Eq.(27). As the flow velocity around the tube is zero due to viscosity, the corresponding boundary condition is
(44)
Besides, according to the symmetrical geometry of the tube, it implies another boundary condition
(45)
Applying the following dimensionless variables
(46)
to Eqs.(43) and (45), it results in
(47)
(48)
where
(49)
In the following analysis will be assigned without loss of generality. Now we define the differential operator
(50)
and its inversion
(51)
and then apply the operator to Eq.(50), it gvies
(52)
Applying the inverse operator to above equation results in
(53)
which can be further written as
(54)
where indicates the velocity at . Now the first term of the decomposition series is assigned by
(55)
The remaining terms can be recursively calculated and are shown below
(56)
The complete solution is displayed
(57)
Using the boundary condition gives
(58)
where
(59)
Similar to the idea adopted in previous section, is assumed to be
(60)
where and are unknowns to be determined. Substituting Eq.(60) into Eq.(58) results in
(61)
Now and for can be determined by
(62)
Similar to the analysis in previous section, the convergence behavior of and shown in Table 3 is excellent. Figure 5 displays the velocity profiles along at four different times. Both and are almost the same as the exact solution while behaves slightly poor.
5. Concluding Remarks
In this paper three types of oscillating viscous flows are examined by applying Adomian decomposition method. Different from previous studies for the same problem, the assignment of the first term of the decomposition series is simply the plate velocity for Stokes’ second problem and the pressure gradient for pulsatile flows. In each case approximate solution accurate up to any order can be easily acquired after two constant parameters are determined. It is found that two parameters and the corresponding approximate solution converge very rapidly for all cases. Present results show an excellent agreement with the exact solution which demonstrates the capacity of ADM for simulating realistic flows. Based on the present method, more complicated initial-value and boundary-value problems with periodic variations can be simulated or solved by ADM associated with the Fourier analysis in the future.
Acknowledgements
The financial support from Ministry of Science and Technology of Taiwan with the grants MOST-104-2221-E-270-003-MY2 and NSC-102-2911-I-006-302 is acknowledged.
References