• Views 2,213
• Citations 5
• ePub 24
• PDF 851
`Abstract and Applied AnalysisVolume 2013 (2013), Article ID 380484, 9 pageshttp://dx.doi.org/10.1155/2013/380484`
Research Article

## Analytical Solutions of Boundary Values Problem of 2D and 3D Poisson and Biharmonic Equations by Homotopy Decomposition Method

1Institute for Groundwater Studies, Faculty of Natural and Agricultural Sciences, University of the Free State, Bloemfontein 9300, South Africa
2Department of Mathematics and Institute for Mathematical Research, University of Putra Malaysia, 43400 Serdang, Selangor, Malaysia

Received 13 June 2013; Accepted 18 August 2013

#### Abstract

The homotopy decomposition method, a relatively new analytical method, is used to solve the 2D and 3D Poisson equations and biharmonic equations. The method is chosen because it does not require the linearization or assumptions of weak nonlinearity, the solutions are generated in the form of general solution, and it is more realistic compared to the method of simplifying the physical problems. The method does not require any corrected function or any Lagrange multiplier and it avoids repeated terms in the series solutions compared to the existing decomposition method including the variational iteration method, the Adomian decomposition method, and Homotopy perturbation method. The approximated solutions obtained converge to the exact solution as tends to infinity.

#### 1. Introduction

The numerical solution of Poisson equations and biharmonic equations is an important problem in numerical analysis. A vast arrangement of investigating effort has been published on the development of numerical solution of Poisson equations and biharmonic equations. The finite difference schemes of second and fourth order for the solution of Poisson’s equation in polar coordinates have been derived by Mittal and Gahlaut [1]. A numerical method to interpolate the source terms of Poisson’s equation by using B-spline approximation has been devised by Perrey-Debain and ter Morsche [2]. Sutmann and Steffen [3] proposed compact approximation schemes for the Laplace operator of fourth and sixth order; the schemes are based on a Padé approximation of the Taylor expansion for the discretized Laplace operator. Ge [4] used fourth-order compact difference discretization scheme with unequal mesh sizes in different coordinate directions to solve a 3D Poisson equation on a cubic domain. Gumerov and Duraiswami [5] developed a complete translation theory for the biharmonic equation in three dimensions. Khattar et al. [6] derived a fourth-order finite difference approximation based on arithmetic average discretization for the solution of three-dimensional nonlinear biharmonic partial differential equations on a 19-point compact stencil using coupled approach. Altas et al. [7] used multigrid and preconditioned Krylov iterative methods to solve three-dimensional nonlinear biharmonic partial differential equations. Jeon [8] derived scalar boundary integral equation formulas for both interior and exterior biharmonic equations with the Dirichlet boundary data. A spectral collocation method for numerically solving two-dimensional biharmonic boundary-value problems has been reported in [9]. An indirect radial-basis-function collocation method for numerically solving biharmonic boundary-value problems has been reported in [10]. A high-order boundary integral equation method for the solution of biharmonic equations has been presented in [11]. A Galerkin boundary node method for solving biharmonic problems was developed in [12]. An integral collocation approach based on Chebyshev polynomials for numerically solving biharmonic equations for the case of irregularly shaped domains has been developed by Mai-Duy et al. [13]. A numerical method, based on neural-network-based functions, for solving partial differential equations has been in [14]. Mai-Duy and Tanner [15] presented a collocation method based on a Cartesian grid and a 1D integrated radial basis function scheme for numerically solving partial differential equations in rectangular domains and Haar wavelet presented in [16]. The aim of this paper is to solve these problems via the homotopy decomposition method.

#### 2. Method

In this study we follow the method of [1720]. In order to illustrate the basic idea of this method we consider a general nonlinear nonhomogeneous partial differential equation with initial conditions of the following form subject to the initial conditions where is a known function, is the general nonlinear differential operator, and represents a linear differential operator. The method’s first step here is to apply the inverse operator of on both sides (1) to obtain The multi-integrals in (3) can be transformed to So that (3) can be reformulated as Using the homotopy scheme the solution of the previous integral equation is given in a series form as and the nonlinear term can be decomposed as where is an embedding parameter. is He’s polynomials [21] that can be generated by The homotopy decomposition method is obtained by the graceful coupling of decomposition method with He’s polynomials and is given by with Comparing the terms of the same power of gives the solutions of various orders. The initial guess of the approximation is . Some further related results can be seen in [2225].

Lemma 1 (see [17]). The complexity of the homotopy decomposition method is of order .

Proof. The number of computations including product, addition, subtraction, and division are as follows.
In step 2: 0 because it is obtained directly from the initial conditions : 3: 3.
Now in step 4 the total number of computations is equal to .

#### 3. Solutions of the Main Problems

Problem 1. Consider the following equation The exact solution of the previous equation is given as In the view of the homotopy decomposition method, (11) can be first transformed to Following the decomposition techniques, we obtain the following equation Comparing the terms of the same power of leads to The following solutions are obtained: In the same manner one can obtain the rest of the components. But for eight terms were computed and the asymptotic solution is given by Therefore in general for any we have This is the exact solution of the problem. Figures 1 and 2 show the comparison of the exact solution and the approximated one for . The approximate solution and the exact solution are compared in Figures 1 and 2, respectively.
The numerical errors for are evaluated in Table 1.

Table 1: Evaluation of numerical errors for .
Figure 1: Exact solution.
Figure 2: Approximated solution for the 4 first terms.

Problem 2. Consider 3D Poisson equation: Following the discussion presented earlier we obtain the following set of integral equations: The following solutions are obtained: In the same manner one can obtain the rest of the components. But for six terms were computed and the asymptotic solution is given by Therefore, for any , the partial sum is given as Thus And this is the exact solution to the problem. One can evaluate error committed by choosing the first terms in the series solutions, in the same manner as in Table 1. The accuracy of the results is estimated by error function

Problem 3. Let us consider the following biharmonic equation for which the exact solution is The aim of this part is to compare the numerical results obtained via HDM and the method used in [26].
Applying the steps involved in the HDM, we arrive at the following: In the same manner, one can obtain the remaining term by using the following recursive formula: In this paper we consider only the first six terms of the series solution as follows: To access the accuracy of the method used in paper, we compare in Table 2 the numerical results of the above equation, the solution obtained in [26] with the exact solution.

Table 2: Comparison of the HDM and [1] results with the exaction solution for .

Problem 4. We consider the 2D biharmonic equation subject to the initial conditions: In the view of the homotopy decomposition method, the following integral equations are obtained: It is worth noting that if the zeroth component is defined, then the remaining components can be completely determined such that each term is determined by using the previous terms, and the series solutions are thus entirely determined. Finally, the solution is approximated for : Therefore for any we have the following: Thus The exact solution of (31) is given by Figures 3 and 4 are the graphical representation of the previous solution. We have plotted the solution for (31) in Figure 3 and showed absolute value of the solution in Figure 4.

Figure 3: Analytical solution.
Figure 4: Absolute value of the solution.

Theorem 2. Let be a nonzero natural number and let ; then two dimensional biharmonic equation of form with along the boundaries has an exact solution as follows

Proof. Use the step of the homotopy decomposition method.

Problem 5. We consider the 3D biharmonic equation: In the view of the homotopy decomposition method, the following integral equations are obtained: Solving the previous integral equations, the series solutions for the first terms are given as Therefore taking the limit at tending to infinity we obtained

#### 4. Conclusion

In this paper the recent homotopy decomposition [1821] is used to solve the 2D and 3D Poisson equations and biharmonic equations. The method is chosen because it does not require the linearization or assumptions of weak nonlinearity, the solutions are generated in the form of general solution, and it is more realistic compared to the method of simplifying the physical problems. The method does not require any corrected function any Lagrange multiplier and it avoids repeated terms in the series solutions compared to the existing decomposition method including the variational iteration method and the Adomian decomposition method. The approximated solutions obtained converge to the exact solution as tends to infinity. The numerical values are presented in Table 1 shows that the method is very efficient and accurate.

#### Acknowledgment

The authors would like to thank referee(s) for very useful comments regarding the details and their remarks that improved the presentation and the contents of the paper.

#### References

1. R. C. Mittal and S. Gahlaut, “High-order finite-difference schemes to solve Poisson's equation in polar coordinates,” IMA Journal of Numerical Analysis, vol. 11, no. 2, pp. 261–270, 1991.
2. E. Perrey-Debain and H. G. ter Morsche, “B-spline approximation and fast wavelet transform for an efficient evaluation of particular solutions for Poisson's equation,” Engineering Analysis with Boundary Elements, vol. 26, no. 1, pp. 1–13, 2002.
3. G. Sutmann and B. Steffen, “High-order compact solvers for the three-dimensional Poisson equation,” Journal of Computational and Applied Mathematics, vol. 187, no. 2, pp. 142–170, 2006.
4. Y. Ge, “Multigrid method and fourth-order compact difference discretization scheme with unequal meshsizes for 3D Poisson equation,” Journal of Computational Physics, vol. 229, no. 18, pp. 6381–6391, 2010.
5. N. A. Gumerov and R. Duraiswami, “Fast multipole method for the biharmonic equation in three dimensions,” Journal of Computational Physics, vol. 215, no. 1, pp. 363–383, 2006.
6. D. Khattar, S. Singh, and R. K. Mohanty, “A new coupled approach high accuracy numerical method for the solution of 3D non-linear biharmonic equations,” Applied Mathematics and Computation, vol. 215, no. 8, pp. 3036–3044, 2009.
7. I. Altas, J. Erhel, and M. M. Gupta, “High accuracy solution of three-dimensional biharmonic equations,” Numerical Algorithms, vol. 29, no. 1–3, pp. 1–19, 2002.
8. Y. Jeon, “New indirect scalar boundary integral equation formulas for the biharmonic equation,” Journal of Computational and Applied Mathematics, vol. 135, no. 2, pp. 313–324, 2001.
9. N. Mai-Duy and R. I. Tanner, “A spectral collocation method based on integrated Chebyshev polynomials for two-dimensional biharmonic boundary-value problems,” Journal of Computational and Applied Mathematics, vol. 201, no. 1, pp. 30–47, 2007.
10. N. Mai-Duy and T. Tran-Cong, “Solving biharmonic problems with scattered-point discretization using indirect radial-basis-function networks,” Engineering Analysis with Boundary Elements, vol. 30, no. 2, pp. 77–87, 2006.
11. N. Mai-Duy and R. I. Tanner, “An effective high order interpolation scheme in BIEM for biharmonic boundary value problems,” Engineering Analysis with Boundary Elements, vol. 29, no. 3, pp. 210–223, 2005.
12. X. Li and J. Zhu, “A Galerkin boundary node method for biharmonic problems,” Engineering Analysis with Boundary Elements, vol. 33, no. 6, pp. 858–865, 2009.
13. N. Mai-Duy, H. See, and T. Tran-Cong, “A spectral collocation technique based on integrated Chebyshev polynomials for biharmonic problems in irregular domains,” Applied Mathematical Modelling, vol. 33, no. 1, pp. 284–299, 2009.
14. M. W. M. G. Dissanayake and N. Phan-Thien, “Neural-network-based approximations for solving partial differential equations,” communications in Numerical Methods in Engineering, vol. 10, no. 3, pp. 195–201, 1994.
15. N. Mai-Duy and R. I. Tanner, “A collocation method based on one-dimensional RBF interpolation scheme for solving PDEs,” International Journal of Numerical Methods for Heat & Fluid Flow, vol. 17, no. 2, pp. 165–186, 2007.
16. Z. Shi, Y.-Y. Cao, and Q.-J. Chen, “Solving 2D and 3D Poisson equations and biharmonic equations by the Haar wavelet method,” Applied Mathematical Modelling, vol. 36, no. 11, pp. 5143–5161, 2012.
17. A. Atangana and J. F. Botha, “Analytical solution of the ground water flow equation obtained via homotopy decomposition method,” Journal of Earth Science & Climatic Change, vol. 3, p. 115, 2012.
18. A. Atangana and A. Secer, “The time-fractional coupled-Korteweg-de-Vries equations,” Abstract and Applied Analysis, vol. 2013, Article ID 947986, 8 pages, 2013.
19. A. Atangana and A. Ahmed, “A generalized version of a low velocity impact between a rigid sphere and a transversely isotropic strain-hardening plate supported by a rigid substrate using the concept of non-integer derivatives,” Abstract Applied Analysisaccepted, vol. 2013, Article ID 671321, 9 pages, 2013.
20. A. Atangana and E. Alabaraoye, “Solving system of fractional partial differential equations arisen in the model of HIV infection of CD4+ cells and attractor one-dimensional Keller-Segel equations,” Advance in Difference Equation, vol. 2013, no. 94, 2013.
21. S. T. Mohyud-Din, A. Yildirim, and M. M. Hosseini, “Variational iteration method for initial and boundary value problems using He's polynomials,” International Journal of Differential Equations, vol. 2010, Article ID 426213, 28 pages, 2010.
22. A. A. Elbeleze, A. Kılıçman, and B. M. Taib, “Homotopy perturbation method for fractional Black-Scholes European option pricing equations using Sumudu transform,” Mathematical Problems in Engineering, vol. 2013, Article ID 524852, 7 pages, 2013.
23. A. A. Elbeleze, A. Kılıçman, and B. M. Taib, “Application of homotopy perturbation and variational iteration method for Fredholm integro-differential equation of fractional order,” Abstract and Applied Analysis, vol. 2012, Article ID 763139, 14 pages, 2012.
24. A. M. A. El-Sayed and D. Hammad, “A reliable treatment of homotopy perturbation method for solving the nonlinear Klein-Gordon equation of arbitrary (fractional) orders),” Journal of Applied Mathematics, vol. 2012, Article ID 581481, 13 pages, 2012.
25. M. T. Kajani, M. Ghasemi, and E. Babolian, “Comparison between the homotopy perturbation method and the sine-cosine wavelet method for solving linear integro-differential equations,” Computers & Mathematics with Applications, vol. 54, no. 7-8, pp. 1162–1168, 2007.
26. A. K. Khalifa, “The decomposition method for one dimensional biharmonic equations,” International Journal of Simulation and Process Modelling, vol. 2, no. 1-2, 2006.