Numerical Treatment of Initial-Boundary Value Problems with Mixed Boundary Conditions

In this paper, we extend the reliable modification of the Adomian Decomposition Method coupled to the Lesnic’s approach to solve boundary value problems and initial boundary value problems with mixed boundary conditions for linear and nonlinear partial differential equations. The method is applied to different forms of heat and wave equations as illustrative examples to exhibit the effectiveness of the method. The method provides the solution in a rapidly convergent series with components that can be computed iteratively. The numerical results for the illustrative examples obtained show remarkable agreement with the exact solutions. We also provide some graphical representations for clear-cut comparisons between the solutions using Maple software.


Introduction
Mixed boundary value problems are characterized by a combination of Dirichlet and Neumann conditions along at least one boundary condition.They occur in a wide range of engineering and applied mathematics applications [1] [2] [3] [4].
These applications include the classic electrical potential and electric field conditions on a disk [5], stress and strain conditions around a punch pressing on an elastic surface [6] as well as some applications in porous media problems such as the infiltration and seepage among others [7] [8].Historically, only very few of these problems could be solved using analytic methods.In view of this, many researchers obtained the solutions of initial and boundary value problems by using either initial or boundary condition(s).In recent years, there have been significant developments in the use of various semi-analytical methods for partial differential equations such as the homotopy perturbation method [9] and Adomian Decomposition Method (ADM) [10].Duan and Rach [11] developed a new resolution method of Boundary Value Problems (BVPs) for nonlinear ordinary differential equations using the ADM.It is also well-known that the ADM provides approximate analytic solutions without using the Green function concept, which greatly facilitates analytic approximations and numerical computations.Several different resolution techniques for solving BVPs for nonlinear ordinary differential equations by using the ADM were considered by Adomian and Rach [12]- [18], Adomian [19], and Wazwaz [20]- [26].Also, for a two-point BVP for second-order nonlinear differential equations, Adomian and Rach [17] [18] proposed the double decomposition method in order to avoid solving such nonlinear algebraic equations, and Jang [27] and Ebaid [28] introduced different modified inverse linear operators.Adomian [29] suggested a modified method for the hyperbolic, parabolic and elliptic partial differential equations with initial and boundary conditions by using two equations for u, one inverting the t L operator and the other inverting the x L operator, then, adding them and di- viding by two.Further, with regards to the mixed value problems, Lesnic and Elliot [30] proposed the inverse operator defined by are known functions.However, in this paper, we will present a modified recursion scheme based on the reliable modification of the ADM with new structure of the inverse operator applied to the (BVPs) with mixed boundary conditions using Lesnic's approach and Ebaid's method.The proposed operator allows the appearance of all the conditions in the solution thereby making the solution more realistic.The paper is arranged in the following manner: in Section 2, we analyze the ADM; Section 3 presents the modified ADM suggested by Wazwaz; in Section 4, Lesnic's approach is used to approximate solutions of some problems; the implementation of this new method to some test problems is presented in Section 5; finally, a brief conclusion is given in Section 6.

Analysis of the Adomian Decomposition Method with Mixed Conditions
Nonlinear partial differential equations models in mathematics and physics play an important role in theoretical sciences.The understanding of these nonlinear partial differential equations is also crucial to many applied areas such as meteorology, oceanography, and aerospace industry.Nonlinear partial differential equations are the most fundamental models in studying nonlinear phenomena.
Consider the nonlinear partial differential equation given in an operator form ( ) , Nu x t is a nonlinear operator which is assumed to be analytic function, and ( ) , g x t is the input function that is assumed to be con- tinuous function.The solutions for ( ) , u x t obtained from the operator equa- tions x L u in x-direction and t L u in t-direction are called partial solutions.
We further give the following illustrations:

Boundary Value Problems
Consider the general form of the single second-order nonlinear inhomogeneous temporal-spatial partial differential equation: subject to the mixed boundary conditions , , , where, , , , , , . We consider the x partial solution as Applying the two-fold indefinite integration inverse operator where ( ) ( ) where n A 's are the Adomian polynomials determined from the definitional formula Substituting Equation (2.6) into Equation (2.5), yields the following recursion scheme ( ) The n-term approximation of the solution is ( ) ϕ is now determined.Since 0 u and 1 ϕ are now completely known, we form the next term

An Alternative Combination of the Initial and Boundary Conditions
Adomian [29] suggested a modified method for the partial differential equations with initial and boundary conditions by using two equations for u, one inverting the t L operator and the other inverting the x L operator, then, adding them and dividing by two.
, , 0 , 0 , , , Secondly, we consider the x partial solution as in Equation (2.5).Next, we averaged the partial solutions, i.e. add two partial solutions in Equation (2.5) and Equation (2.11) and divide by two, we obtain To illustrate this method for coupled linear and nonlinear partial differential equations, we take two examples in the following section.
Consider the linear homogeneous heat equation 0, 0 1, 0, Rewriting the heat equation in the operator form as 0 Applying the inverse operator , , .
This introduces the recursive relations ( ) ( ) The first approximant is . Applying the x conditions to 1 ϕ , it is clear that ( ) . Thus, if the one-term approximant 1 ϕ were sufficient, the solution would be ( ) The next term is Applying the condition at Continuing in a similar way 2 3 , , , n u u u  are obtained for some n, then we get the approximate solution which converged to the exact solution ( ) Table 1 shows the comparison between the absolute error of the exact and approximate solutions for various values of t. Figure 1 and Figure 2 give the plots for the exact and approximate solutions by using ADM for 0.5, 1.0 t t = = and 0 1 x ≤ ≤ , respectively.American Journal of Computational Mathematics , and the boundary conditions ( ) ( ) To solve initial-boundary value problem, firstly, we consider the t partial solution as Applying the inverse operator This introduces the recursive relations where the nonlinear term 2 u can be expressed by an infinite series of the Adomian polynomials n A given by: ( ) Secondly, we consider the x partial solution Applying the inverse operator This gives the recursive relations ( ) Thus, if the one-term approximant 1 ϕ were sufficient, the solution would be is given by: ( ) ( ) Continuing in a similar way 2 3 , , , n u u u  are obtained for some n, then we get the approximate solution which is converge to the exact solution Table 2 shows the comparison between the absolute error of exact and approximate solutions for various values of t. Figure 3 and Figure 4 give the plots for the exact solution and the approximate solution by application of (ADM) for 0.5, 1.0 t t = =

Modified Adomian Decomposition Method with Mixed Conditions
Wazwaz in [31], suggested a modification to the recurrence relations of the Standard Adomian Decomposition Method (ADM).He assumed that the first   ϕ ϕ ϕ at 1.0 t = .American Journal of Computational Mathematics term u 0 can be set as function f and divided it into two parts, namely f 1 and f 2 .
Under this assumption, we set f f f = + .Also, based on this, we propose a slight variation only on the components u 0 and u 1 , the variation is that only the part f 1 is assigned to the zeroth component u 0 , whereas the remaining part f 2 is combined with the other terms to define u 1 .This reduction in the number of terms of u 0 will result in reduction of the computational work and will accelerate the convergence.Further, this slight variation in the definition of the components u 0 and u 1 may provide the solution by using two iterations only.Furthermore, the calculations below will show that sometimes there is no need to evaluate the so-called Adomian polynomials required for the nonlinear differential equations.An important observation that can be made here is that the success of this method depends mainly on the proper choice of the parts f 1 and f 2 .We have been unable to establish any criterion to judge what forms of f 1 and f 2 can be assumed to yield the acceleration demanded.It appears that trials are the only criteria that can be applied so far.

Boundary Value Problems
Based on the recurrence relation in Equation (2.8) We can set ( ) , then we divide it into two parts, so that we formulate the modified recursive algorithm as follows: Comparing the recursive scheme in Equation (2.8) of the ADM with the scheme in Equation (3.1) of the modified technique leads to the conclusion that in Equation (2.8) the zeroth component was defined by the function f, whereas in Equation (3.1), the zeroth component u 0 was defined by only a part f 1 of f, the remaining part f 2 of f is added to the definition of the component u 1 in (3.1).

An alternative Combination of the Initial and Boundary Conditions
Based on the recurrence relation in Equation (2.13), we can set Then we divide it into two parts, so that we formulate the modified recursive algorithm as follows: Comparing the recursive scheme in Equation.(2.13) of the ADM with the recursive scheme in Equation (3.2) of the modified technique leads to the conclusion that in (2.13) the zeroth component was defined by the function f, whereas in (3.2), the zeroth component u 0 is defined only by a part f 1 of f, the remaining part f 2 of f is added to the definition of the component u 1 in (2.13).
In order to demonstrate the efficiency and applicability of the proposed method, we study example of nonlinear partial differential equations here.

Numerical Experiments
Example 3 Consider the nonlinear inhomogeneous heat equation , 0, Rewriting the heat equation in the operator form as , , , , This gives the recursive relations ( ) Thus, if the one-term approximant 1 ϕ were sufficient, the "solution" will be then we divide it into two parts, so that is given by: ( ) ( ) Applying the condition at 0 x = , we have 1,1 ( ) 0 c t = .From the condition on American Journal of Computational Mathematics x at 0, we get ( )

Improvement of the Inverse Operator with Mixed Boundary Conditions
In 1999, Lesnic and Elliott [30] employed ADM for solving some inverse boundary value problems in heat conduction, also they applied the modification method to deal with noisy input data and obtain a stable approximate solution.In [32] Lesnic investigated the application of the decomposition method involving computational algebra for solving more complicated problems with Dirichlet, or mixed boundary conditions.In [30] the authors defined the operator 1 t L − and as definite integrals given by to solve the linear homogeneous heat equation 0 , 1, 0 to the mixed boundary conditions ( ) ( ) ( ) ( ) x u x t h t u t h t = = .
We will make an extension to the inverse operator (4.1) given in [30] to all cases of problems, so that we consider in this section, two types of problems: boundary value problems and initial-boundary value problems.

Boundary Value Problems
Consider Equations.(2.2)-( 2.3), and applying the inverse operator to both sides of Equation (2.4), gives i.e. the boundary conditions can be used directly to solve the boundary value problem in x-direction.

An Alternative Combination of the Initial and Boundary Conditions
Here, we take the average of the t partial solution as in Equation (2.11), and the x so that the recurrence relations are To give a clear overview of these methods, we have chosen several differential equations.The examples will be approached by the Lesnic's Approach and the modified technique for comparison reasons.We also compare the approximate solution with the exact solution.

Numerical Experiments
Example 4 Consider the linear homogeneous heat equation 0, 0 1, 0, ,0 e , 0, e , 1, e Rewriting the heat equation in the operator form as t xx L u L u = .Applying the inverse operator  Table 3 shows the comparison between the absolute error of exact and approximate solutions for various values of t. Figure 5 and Figure 6 show the results for the exact solution and the approximate solution by application of Lesnic's approach for 0.5, 1.0 Consider the linear homogeneous heat equation , 0, Rewriting the heat equation in the operator form a Applying the inverse operator Continuing in a similar way 3 4 , , , n u u u  are obtained for some n, then we get the approximate solution Table 4 shows the comparison between the absolute error of exact and approximate solutions for various values of t. Figure 7 and Figure 8 show the results for the exact solution and the approximate solution by application of Lesni's approach for 0.5, 1.0 t t = = .
Example 6 Consider the nonlinear inhomogeneous wave equation and the boundary conditions ( ) ( ) Rewriting the wave equation in the operator form as Continuing in a similar way Table 5 shows the comparison between the absolute error of exact and approximate solutions for various values of t. Figure 9 and Figure 10 show the    Table 6 shows the comparison between the absolute error of exact and approximate solutions Lesnic's approach and modified of Lesnic's approach for various values of t. Figure 11 and Figure 12 show the results for the exact solution and the approximate solution by application of Lesnic's approach and modified of Lesnic's Approach for 0.5 t = .

Conclusion
The main goal of this paper is the development of an efficient and accurate me-

Figure 4 .
Figure 4.The exact solution and the approximate solution using ADM for 3 5 10, , the exact solution.
thod to solve linear and nonlinear initial-boundary value problems with mixed boundary conditions.The method developed provides the solution in the form of rapidly convergent series relative to the Adomian decomposition method.It also introduces a slight change in the formulation of Adomian recursive relation.The method presented is applied to several examples of boundary-value problems and initial boundary-value problems with mixed type boundary conditions.
., and all n ,

Table 1 .
Absolute errors using ADM at

Table 2 .
Absolute errors using ADM at

Table 3 .
Absolute errors of Lesnic's approach at

Table 4 .
Absolute errors of Lesnic's approach at

Table 5 .
Absolute errors of Lesnic's approach at American Journal of Computational Mathematics results for the exact solution and the approximate solution by application of , , ϕ ϕ ϕ at 1.0 t =

Table 6 .
Absolute errors of Lesnic's approach and modified of Lesnic's Approach at