**American Journal of Computational Mathematics**

Vol.08 No.02(2018), Article ID:85604,22 pages

10.4236/ajcm.2018.82012

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

Nawal Abdullah Alzaid, Huda Omar Bakodah^{ }

Department of Mathematics, King Abdulaziz University, Jeddah, KSA

Copyright © 2018 by authors and Scientific Research Publishing Inc.

This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).

http://creativecommons.org/licenses/by/4.0/

Received: April 11, 2018; Accepted: June 25, 2018; Published: June 28, 2018

ABSTRACT

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.

**Keywords:**

Decomposition Method, Modified Adomian Decomposition Method, Linear and Nonlinear Partial Differential Equations, Mixed Boundary Conditions, Initial-Boundary Value Problem

1. 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 ${L}_{t}$ operator and the other inverting the ${L}_{x}$ operator, then, adding them and dividing by two. Further, with regards to the mixed value problems, Lesnic and Elliot

[30] proposed the inverse operator defined by ${L}_{x}^{-1}={\displaystyle \underset{{x}_{0}}{\overset{x}{\int}}\text{d}{x}^{\prime}}{\displaystyle \underset{1}{\overset{x}{\int}}\text{d}{x}^{\u2033}}$ to solve the

linear homogeneous heat equation ${u}_{t}={u}_{xx},\text{\hspace{0.17em}}{x}_{0}<x<1,\text{\hspace{0.17em}}t>0$ subject to the mixed boundary conditions $u\left(x,t\right)={h}_{1}\left(t\right)$ , ${u}_{x}\left(1,t\right)={h}_{2}\left(t\right)$ , where ${\beta}_{i}=\left(t\right),i=1,2$ 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.

2. 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

${L}_{x}u\left(x,t\right)+{L}_{t}u\left(x,t\right)+Ru\left(x,t\right)+Nu\left(x,t\right)=g\left(x,t\right)$ (2.1)

where ${L}_{x}$ and ${L}_{t}$ are the linear operators to be inverted, which are usually the highest order differential operators in x, and t respectively; ${L}_{x}$ , R is the linear remainder operator; $Nu\left(x,t\right)$ is a nonlinear operator which is assumed to be analytic function, and $g\left(x,t\right)$ is the input function that is assumed to be continuous function. The solutions for $u\left(x,t\right)$ obtained from the operator equations ${L}_{x}u$ in x-direction and ${L}_{t}u$ in t-direction are called partial solutions. We further give the following illustrations:

2.1. Boundary Value Problems

Consider the general form of the single second-order nonlinear inhomogeneous temporal-spatial partial differential equation:

${L}_{xx}u\left(x,t\right)+{L}_{tt}u\left(x,t\right)+Nu\left(x,t\right)=g\left(x,t\right),\text{\hspace{0.17em}}a\le x\le b,\text{\hspace{0.17em}}t>0,$ (2.2)

subject to the mixed boundary conditions

$u\left(a,t\right)={h}_{1}\left(t\right),\text{\hspace{0.17em}}{u}_{x}\left(b,t\right)={h}_{2}\left(t\right)$ (2.3)

where, ${L}_{xx}u\left(x,t\right)=\frac{{\partial}^{2}}{\partial {x}^{2}}u\left(x,t\right),{L}_{tt}u\left(x,t\right)=\frac{{\partial}^{2}}{\partial {t}^{2}}u\left(x,t\right)$ . We consider the x partial solution as

${L}_{xx}u\left(x,t\right)=g\left(x,t\right)-{L}_{tt}u\left(x,t\right)-Nu\left(x,t\right)$ (2.4)

Applying the two-fold indefinite integration inverse operator ${L}_{xx}^{-1}={\displaystyle \int \text{d}x}{\displaystyle \int \text{d}{x}^{\prime}}$ to both sides of Equation (2.4), gives

$u\left(x,t\right)={\Phi}_{x}+{L}_{xx}^{-1}g\left(x,t\right)-{L}_{xx}^{-1}{L}_{tt}u\left(x,t\right)-{L}_{xx}^{-1}Nu\left(x,t\right),$ (2.5)

where ${\Phi}_{x}={c}_{1}\left(t\right)+{c}_{2}\left(t\right)x$ and the constants of integrations ${c}_{1}\left(t\right)$ and ${c}_{2}\left(t\right)$ , are determined from the boundary conditions. We now decompose the following ${\Phi}_{x}$ , the linear and nonlinear terms u and Nu based on ADM as follows:

${\Phi}_{x}={\displaystyle \sum _{n=0}^{\infty}{\Phi}_{x,n}=}{\displaystyle \sum _{n=0}^{\infty}\left({c}_{1,n}\left(t\right)+{c}_{2,n}\left(t\right)x\right)},$

$u\left(x,t\right)={\displaystyle \sum _{n=0}^{\infty}{u}_{n}\left(x,t\right),}$

$Nu\left(x,t\right)={\displaystyle \sum _{n=0}^{\infty}{A}_{n}\left({u}_{0},{u}_{1},{u}_{2},\cdots ,{u}_{n}\right)}$ (2.6)

where ${A}_{n}$ ’s are the Adomian polynomials determined from the definitional formula

${A}_{n}={\frac{1}{n!}\frac{{\text{d}}^{n}}{\text{d}{\lambda}^{n}}N\left({\displaystyle \sum _{k=0}^{\infty}{u}_{k}{\lambda}^{k}}\right)|}_{\lambda =0},\text{\hspace{0.17em}}n\ge 0.$ (2.7)

Substituting Equation (2.6) into Equation (2.5), yields the following recursion scheme

${u}_{0}={\Phi}_{x,0}+{L}_{xx}^{-1}g\left(x,t\right),$

${u}_{n+1}={\Phi}_{x,n+1}-{L}_{xx}^{-1}{L}_{tt}{u}_{n}-{L}_{xx}^{-1}{A}_{n},\text{\hspace{0.17em}}\text{\hspace{0.17em}}n\ge 0.$ (2.8)

The n-term approximation of the solution is ${\phi}_{n}={\displaystyle \sum _{i=0}^{n-1}{u}_{i}\left(x,t\right)}$ for $n>0$ .

Thus, ${\phi}_{1}={u}_{0},{\phi}_{2}={\phi}_{1}+{u}_{1},{\phi}_{3}={\phi}_{2}+{u}_{2}$ , etc., and all ${\phi}_{n}$ ‘s must satisfy the boundary conditions.

The first approximate ${\phi}_{1}={u}_{0}={c}_{1,0}\left(t\right)+{c}_{2,0}\left(t\right)x+{L}_{xx}^{-1}g\left(x,t\right)$ , where the values ${c}_{1}\left(t\right)$ and ${c}_{2}\left(t\right)$ can be evaluated by using the boundary conditions in Equation (2.2)

${{\phi}_{1}|}_{x=a}={h}_{1}\left(t\right),\text{\hspace{0.17em}}\text{\hspace{0.17em}}{{{\phi}^{\prime}}_{1}|}_{x=b}={h}_{2}\left(t\right),$

which results in

${{c}_{1,0}\left(t\right)+{c}_{2,0}\left(t\right)x+{L}_{xx}^{-1}g\left(x,t\right)|}_{x=a}={h}_{1}\left(t\right),$

${{c}_{2,0}\left(t\right)+{\left({L}_{xx}^{-1}g\left(x,t\right)\right)}^{\prime}|}_{x=b}={h}_{2}\left(t\right).$

Thus ${\phi}_{1}$ is now determined. Since ${u}_{0}$ and ${\phi}_{1}$ are now completely known, we form the next term ${u}_{1}={\Phi}_{x,1}-{L}_{xx}^{-1}{L}_{tt}{u}_{0}-{L}_{xx}^{-1}{A}_{0}$ , then ${\phi}_{2}={\phi}_{1}+{u}_{1}$ , and we continue in the same manner to obtain ${u}_{2},{u}_{3},\cdots ,{u}_{n}$ for some $n\ge 0$ . Substituting all these values in $u\left(x,t\right)={\displaystyle \sum _{n=0}^{\infty}{u}_{n}\left(x,t\right)}$ , we get the solution of Equation (2.1).

2.2. 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 ${L}_{t}$ operator and the other inverting the ${L}_{x}$ operator, then, adding them and dividing by two. To convey the basic idea for treatment of initial and boundary conditions by ADM for solving initial boundary value problems, we consider Equations (2.2)-(2.3) with the initial conditions

$u\left(x,0\right)={p}_{1}\left(x\right),\text{\hspace{0.17em}}\text{\hspace{0.17em}}{u}_{t}\left(x,0\right)={p}_{2}\left(x\right).$ (2.9)

where ${p}_{i}\left(x\right),i=1,2$ are known functions.

Firstly, we consider the t partial solution as

${L}_{tt}u\left(x,t\right)=g\left(x,t\right)-{L}_{xx}u\left(x,t\right)-Nu\left(x,t\right)$ (2.10)

Applying the inverse operator ${L}_{tt}^{-1}$ defined by ${L}_{tt}^{-1}={\displaystyle {\int}_{0}^{t}{\displaystyle {\int}_{0}^{t}\text{d}t\text{\hspace{0.05em}}\text{d}t}}$ to both sides of Equation (2.10) and using the initial conditions, gives

$u\left(x,t\right)=u\left(x,0\right)+t{u}_{t}\left(x,0\right)+{L}_{tt}^{-1}g\left(x,t\right)-{L}_{tt}^{-1}{L}_{xx}u\left(x,t\right)-{L}_{tt}^{-1}Nu\left(x,t\right)$ (2.11)

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

$\begin{array}{l}u\left(x,t\right)=\frac{1}{2}(u\left(x,0\right)+t{u}_{t}\left(x,0\right)+{\Phi}_{x}+{L}_{tt}^{-1}g\left(x,t\right)+{L}_{xx}^{-1}g\left(x,t\right)\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}-{L}_{tt}^{-1}{L}_{xx}u\left(x,t\right)-{L}_{xx}^{-1}{L}_{tt}u\left(x,t\right)-{L}_{tt}^{-1}Nu\left(x,t\right)-{L}_{xx}^{-1}Nu\left(x,t\right))\end{array}$ (2.12)

Substituting Equation (2.6) into Equation (2.12), yields the following recursion scheme

${u}_{0}=\frac{1}{2}\left(u\left(x,0\right)+t{u}_{t}\left(x,0\right)+{\Phi}_{x,0}+{L}_{tt}^{-1}g\left(x,t\right)+{L}_{xx}^{-1}g\left(x,t\right)\right)$

${u}_{n+1}=\frac{1}{2}\left({\Phi}_{x,n}-{L}_{tt}^{-1}{L}_{xx}{u}_{n}\left(x,t\right)-{L}_{xx}^{-1}{L}_{tt}{u}_{n}\left(x,t\right)-{L}_{tt}^{-1}{A}_{n}-{L}_{xx}^{-1}{A}_{n}\right),\text{\hspace{0.17em}}n\ge 0.$ (2.13)

To illustrate this method for coupled linear and nonlinear partial differential equations, we take two examples in the following section.

2.3. Numerical Experiments

Example 1.

Consider the linear homogeneous heat equation

${u}_{t}-{u}_{xx}=0,\text{\hspace{0.17em}}\text{\hspace{0.17em}}0\le x\le 1,\text{\hspace{0.17em}}\text{\hspace{0.17em}}t>0,$

with specified conditions $u\left(0,t\right)=0,{u}_{x}\left(1,t\right)=\mathrm{cos}\left(1\right){\text{e}}^{-t}$ .

Rewriting the heat equation in the operator form as ${L}_{t}u-{L}_{xx}u=0$ .

Applying the inverse operator ${L}_{xx}^{-1}$ defined by ${L}_{xx}^{-1}={\displaystyle \int \text{d}x}{\displaystyle \int \text{d}{x}^{\prime}}$ , gives

$u\left(x,t\right)={c}_{1}\left(t\right)+{c}_{2}\left(t\right)x+{L}_{xx}^{-1}{L}_{t}u\left(x,t\right).$

This introduces the recursive relations

${u}_{0}={c}_{1}\left(t\right)+{c}_{2}\left(t\right)x,$

${u}_{n+1}={c}_{1,n+1}\left(t\right)+{c}_{2,n+1}\left(t\right)x+{L}_{xx}^{-1}{L}_{t}{u}_{n},\text{\hspace{0.17em}}n\ge 0.$

The first approximant is ${\phi}_{1}={u}_{0}={c}_{1,0}\left(t\right)+{c}_{2,0}\left(t\right)x$ . Applying the x conditions to ${\phi}_{1}$ , it is clear that ${c}_{1,0}\left(t\right)=0$ and ${c}_{2,0}\left(t\right)=\mathrm{cos}\left(1\right){\text{e}}^{-t}$ . Thus, if the one-term approximant ${\phi}_{1}$ were sufficient, the solution would be ${\phi}_{1}={u}_{0}=\mathrm{cos}\left(1\right){\text{e}}^{-t}x$ , The next term is

$\begin{array}{c}{u}_{1}={c}_{1,1}\left(t\right)+{c}_{2,1}\left(t\right)x+{L}_{xx}^{-1}{L}_{t}{u}_{0}\\ ={c}_{1,1}\left(t\right)+{c}_{2,1}\left(t\right)x+{L}_{xx}^{-1}{L}_{t}\left(\mathrm{cos}\left(1\right){\text{e}}^{-t}x\right)\end{array}$

Then ${\phi}_{2}={u}_{0}+{u}_{1}$ is given by

${\phi}_{2}=\mathrm{cos}\left(1\right){\text{e}}^{-t}x+{c}_{1,1}\left(t\right)+{c}_{2,1}\left(t\right)x-\frac{1}{6}\mathrm{cos}\left(1\right){\text{e}}^{-t}{x}^{3}$

Applying the condition at $x=0$ , we have ${c}_{1,1}\left(t\right)=0$ . From the condition on x at 1, we get

${c}_{2,1}\left(t\right)=\frac{1}{2}\mathrm{cos}\left(1\right){\text{e}}^{-t}$ , thus ${u}_{1}=\frac{1}{2}\mathrm{cos}\left(1\right){\text{e}}^{-t}x-\frac{1}{6}\mathrm{cos}\left(1\right){\text{e}}^{-t}{x}^{3}$

Continuing in a similar way ${u}_{2},{u}_{3},\cdots ,{u}_{n}$ are obtained for some n, then we get the approximate solution ${u}_{ap}={\displaystyle \sum _{n=0}^{\infty}{u}_{n}\left(x,t\right)}$ which converged to the exact solution ${u}_{ex}=\mathrm{cos}\left(x\right){\text{e}}^{-t}$ .

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 $t=0.5,t=1.0$ and $0\le x\le 1$ , respectively.

Table 1. Absolute errors using ADM at $t=0.5,t=1.0$ and $0\le x\le 1$ .

Figure 1. The exact solution and the approximate solution using ADM for ${\phi}_{3},{\phi}_{5},{\phi}_{10}$ at $t=0.5$ and $0\le x\le 1$ .

Figure 2. The exact solution and the approximate solution using ADM for ${\phi}_{3},{\phi}_{5},{\phi}_{10}$ at $t=1.0$ and $0\le x\le 1$ .

Example 2

Consider the nonlinear inhomogeneous wave equation

${u}_{tt}-{u}_{xx}-u+{u}^{2}=xt+{x}^{2}{t}^{2},\text{\hspace{0.17em}}\text{\hspace{0.17em}}t>0,$

with specified initial conditions $u\left(x,0\right)=1,{u}_{t}\left(x,0\right)=x$ , and the boundary conditions $u\left(0,t\right)=1,{u}_{x}\left(0,t\right)=t$ . Rewriting the wave equation in the operator form as

${L}_{tt}u-{L}_{xx}u-u+{u}^{2}=xt+{x}^{2}{t}^{2}.$

To solve initial-boundary value problem, firstly, we consider the t partial solution as ${L}_{tt}u={L}_{xx}u+u-{u}^{2}+xt+{x}^{2}{t}^{2}$ .

Applying the inverse operator ${L}_{tt}^{-1}$ defined by ${L}_{tt}^{-1}={\displaystyle {\int}_{0}^{t}{\displaystyle {\int}_{0}^{t}\text{d}t\text{\hspace{0.05em}}\text{d}t}}$ , gives

$u\left(x,t\right)=u\left(x,0\right)+t{u}_{t}\left(x,0\right)+{L}_{tt}^{-1}\left(xt+{x}^{2}{t}^{2}\right)+{L}_{tt}^{-1}{L}_{xx}u\left(x,t\right)+{L}_{tt}^{-1}u-{L}_{tt}^{-1}{u}^{2}$ .

This introduces the recursive relations

${u}_{0}=u\left(x,0\right)+t{u}_{t}\left(x,0\right)+{L}_{tt}^{-1}\left(xt+{x}^{2}{t}^{2}\right)$

${u}_{n+1}={L}_{tt}^{-1}{L}_{xx}{u}_{n}+{L}_{tt}^{-1}{u}_{n}-{L}_{tt}^{-1}{A}_{n},\text{\hspace{0.17em}}\text{\hspace{0.17em}}n\ge 0.$

where the nonlinear term ${u}^{2}$ can be expressed by an infinite series of the Adomian polynomials ${A}_{n}$ given by:

${A}_{0}={\left({u}_{0}\right)}^{2},$

${A}_{1}=2{u}_{0}{u}_{1},$

${A}_{2}={\left({u}_{1}\right)}^{2}+2{u}_{0}{u}_{2},$

${A}_{3}=2{u}_{1}{u}_{2}+2{u}_{0}{u}_{3},$

${A}_{4}={\left({u}_{2}\right)}^{2}+2{u}_{1}{u}_{3}+2{u}_{0}{u}_{4},$

$\vdots $

So that first two terms are

${u}_{0}=1+xt+\frac{1}{6}x{t}^{3}+\frac{1}{12}{x}^{2}{t}^{4},$

$\begin{array}{l}{u}_{1}=\frac{1}{180}{t}^{6}-\frac{1}{12960}{x}^{4}{t}^{10}-\frac{1}{2529}{x}^{3}{t}^{9}-\frac{1}{2016}{x}^{2}{t}^{8}-\frac{1}{252}{x}^{3}{t}^{7}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}-\frac{1}{72}{x}^{2}{t}^{6}-\frac{1}{120}x{t}^{5}-\frac{1}{12}{x}^{2}{t}^{4}-\frac{1}{6}x{t}^{3},\end{array}$

Secondly, we consider the x partial solution ${L}_{xx}u={L}_{tt}u-u+{u}^{2}-xt-{x}^{2}{t}^{2}$ .

Applying the inverse operator ${L}_{xx}^{-1}$ defined by ${L}_{xx}^{-1}={\displaystyle \int \text{d}x{\displaystyle \int \text{d}{x}^{\prime}}}$ , gives

$u\left(x,t\right)={\Phi}_{x}+{L}_{xx}^{-1}\left(-xt-{x}^{2}{t}^{2}\right)+{L}_{xx}^{-1}{L}_{tt}u\left(x,t\right)-{L}_{xx}^{-1}u\left(x,t\right)+{L}_{xx}^{-1}{u}^{2}\left(x,t\right),$

This gives the recursive relations

${u}_{0}={\Phi}_{x,0}+{L}_{xx}^{-1}\left(-xt-{x}^{2}{t}^{2}\right),$

${u}_{n+1}={\Phi}_{x,n}+{L}_{xx}^{-1}{L}_{tt}{u}_{n}-{L}_{xx}^{-1}{u}_{n}+{L}_{xx}^{-1}{A}_{n},\text{\hspace{0.17em}}n\ge 0.$

So that

$\begin{array}{c}{u}_{0}={c}_{1,0}\left(t\right)+{c}_{2,0}\left(t\right)x+{L}_{xx}^{-1}\left(-xt-{x}^{2}{t}^{2}\right)\\ ={c}_{1,0}\left(t\right)+{c}_{2,0}\left(t\right)x-\frac{1}{12}{x}^{4}{t}^{2}-\frac{1}{6}{x}^{3}t,\end{array}$

The first approximant is ${\phi}_{1}={u}_{0}$ . Applying the x conditions to ${\phi}_{1}$ it is clear that ${c}_{1,0}\left(t\right)=1$ and ${c}_{2,0}\left(t\right)=t$ . Thus, if the one-term approximant ${\phi}_{1}$ were

sufficient, the solution would be ${\phi}_{1}={u}_{0}=1+xt-\frac{1}{12}{x}^{4}{t}^{2}-\frac{1}{6}{x}^{3}t$ . The next term is

${u}_{1}={c}_{1,1}\left(t\right)+{c}_{2,1}\left(t\right)x+{L}_{xx}^{-1}{L}_{tt}{u}_{0}-{L}_{xx}^{-1}{u}_{0}+{L}_{xx}^{-1}{A}_{0}.$

Then ${\phi}_{2}={u}_{0}+{u}_{1}$ is given by:

$\begin{array}{l}{\phi}_{2}=1+xt+{c}_{1,1}\left(t\right)+{c}_{2,1}\left(t\right)x-\frac{1}{180}{x}^{6}+\frac{1}{12960}{t}^{4}{x}^{10}+\frac{1}{2592}{t}^{3}{x}^{9}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}+\frac{1}{2016}{t}^{2}{x}^{8}-\frac{1}{252}{t}^{3}{x}^{7}-\frac{1}{72}{t}^{2}{x}^{6}-\frac{1}{120}t{x}^{5},\end{array}$

Applying the condition at $x=0$ , we have ${c}_{1,1}\left(t\right)=0$ . From the condition on x at 0, we get ${c}_{2,1}\left(t\right)=0$ , thus

$\begin{array}{c}{u}_{1}=-\frac{1}{180}{x}^{6}+\frac{1}{12960}{t}^{4}{x}^{10}+\frac{1}{2592}{t}^{3}{x}^{9}+\frac{1}{2016}{t}^{2}{x}^{8}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}-\frac{1}{252}{t}^{3}{x}^{7}-\frac{1}{72}{t}^{2}{x}^{6}-\frac{1}{120}t{x}^{5}+\frac{1}{12}{t}^{2}{x}^{4}+\frac{1}{6}t{x}^{3},\end{array}$

Next, we average the partial solutions, i.e. add two partial solutions and divide by two, so we obtain

${u}_{0}=\frac{1}{2}\left(2+2xt+\frac{1}{12}{x}^{2}{t}^{4}+\frac{1}{6}x{t}^{3}-\frac{1}{12}{x}^{4}{t}^{2}-\frac{1}{6}{x}^{3}t\right)$

$\begin{array}{c}{u}_{1}=\frac{1}{2}(\frac{1}{180}{t}^{6}-\frac{1}{12960}{x}^{4}{t}^{10}-\frac{1}{2592}{x}^{3}{t}^{9}-\frac{1}{2016}{x}^{2}{t}^{8}-\frac{1}{252}{x}^{3}{t}^{7}-\frac{1}{72}{x}^{2}{t}^{6}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}-\frac{1}{120}x{t}^{5}-\frac{1}{12}{x}^{2}{t}^{4}-\frac{1}{6}x{t}^{3}-\frac{1}{180}{x}^{6}+\frac{1}{12960}{t}^{4}{x}^{10}+\frac{1}{2592}{t}^{3}{x}^{9}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}+\frac{1}{2016}{t}^{2}{x}^{8}-\frac{1}{252}{t}^{3}{x}^{7}-\frac{1}{72}{t}^{2}{x}^{6}-\frac{1}{120}t{x}^{5}+\frac{1}{12}{t}^{2}{x}^{4}+\frac{1}{6}t{x}^{3}),\end{array}$

Continuing in a similar way ${u}_{2},{u}_{3},\cdots ,{u}_{n}$ are obtained for some n, then we get the approximate solution ${u}_{ap}={\displaystyle \sum _{n=0}^{\infty}{u}_{n}\left(x,t\right)}$ which is converge to the exact solution ${u}_{ex}=1+xt$ .

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 $t=0.5,t=1.0$

3. 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

Table 2. Absolute errors using ADM at $t=0.5,t=1.0$ and $0\le x\le 1$ .

Figure 3. The exact solution and the approximate solution using ADM for ${\phi}_{3},{\phi}_{5},{\phi}_{10}$ at $t=0.5$ .

Figure 4. The exact solution and the approximate solution using ADM for ${\phi}_{3},{\phi}_{5},{\phi}_{10}$ at $t=1.0$ .

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}_{1}+{f}_{2}$
. 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.

3.1. Boundary Value Problems

Based on the recurrence relation in Equation (2.8)

${u}_{0}={\Phi}_{x,0}+{L}_{xx}^{-1}g\left(x,t\right),$

${u}_{n+1}={\Phi}_{x,n+1}-{L}_{xx}^{-1}{L}_{tt}{u}_{n}-{L}_{xx}^{-1}{A}_{n},\text{\hspace{0.17em}}\text{\hspace{0.17em}}n\ge 0.$

We can set ${u}_{0}={\Phi}_{x,0}+{L}_{xx}^{-1}g\left(x,t\right)=f$ , then we divide it into two parts, so that we formulate the modified recursive algorithm as follows:

${u}_{0}={f}_{1},$

${u}_{1}={f}_{2}+{\Phi}_{x,1}-{L}_{xx}^{-1}{L}_{tt}{u}_{0}-{L}_{xx}^{-1}{A}_{0},$

${u}_{n+1}={\Phi}_{x,n+1}-{L}_{xx}^{-1}{L}_{tt}{u}_{n}-{L}_{xx}^{-1}{A}_{n},\text{\hspace{0.17em}}\text{\hspace{0.17em}}n\ge 1.$ (3.1)

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).

3.2. An alternative Combination of the Initial and Boundary Conditions

Based on the recurrence relation in Equation (2.13), we can set

${u}_{0}=\frac{1}{2}\left(u\left(x,0\right)+t{u}_{t}\left(x,0\right)+{\Phi}_{x,0}+{L}_{tt}^{-1}g\left(x,t\right)+{L}_{xx}^{-1}g\left(x,t\right)\right)=f$

Then we divide it into two parts, so that we formulate the modified recursive algorithm as follows:

${u}_{0}={f}_{1},$

${u}_{1}={f}_{2}+\frac{1}{2}\left({\Phi}_{x,1}-{L}_{tt}^{-1}{L}_{xx}{u}_{0}-{L}_{xx}^{-1}{L}_{tt}{u}_{0}-{L}_{tt}^{-1}{A}_{0}-{L}_{xx}^{-1}{A}_{0}\right),$

${u}_{n+1}=\frac{1}{2}\left({\Phi}_{x,n+1}-{L}_{tt}^{-1}{L}_{xx}{u}_{n}-{L}_{xx}^{-1}{L}_{tt}{u}_{n}-{L}_{tt}^{-1}{A}_{n}-{L}_{xx}^{-1}{A}_{n}\right),\text{\hspace{0.17em}}\text{\hspace{0.17em}}n\ge 1.$ (3.2)

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.

3.3. Numerical Experiments

Example 3

Consider the nonlinear inhomogeneous heat equation

${u}_{t}={u}_{xx}+{u}^{2}-{x}^{2}{t}^{2}+x,\text{\hspace{0.17em}}\text{\hspace{0.17em}}t>0,$

with specified conditions $u\left(x,0\right)=0,u\left(0,t\right)=0,{u}_{x}\left(0,t\right)=t$ .

Rewriting the heat equation in the operator form as

${L}_{t}u={L}_{xx}u+{u}^{2}-{x}^{2}{t}^{2}+x.$

Applying the inverse operator ${L}_{xx}^{-1}$ defined by ${L}_{xx}^{-1}={\displaystyle \int \text{d}x{\displaystyle \int \text{d}{x}^{\prime}}}$ , we get

$u\left(x,t\right)={\Phi}_{x}+{L}_{xx}^{-1}\left({x}^{2}{t}^{2}-x\right)+{L}_{xx}^{-1}{L}_{t}u\left(x,t\right)-{L}_{xx}^{-1}{u}^{2}\left(x,t\right),$

This gives the recursive relations

${u}_{0}={\Phi}_{x,0}+{L}_{xx}^{-1}\left({x}^{2}{t}^{2}-x\right),$

${u}_{n+1}={\Phi}_{x,n+1}+{L}_{xx}^{-1}{L}_{t}{u}_{n}-{L}_{xx}^{-1}{A}_{n},\text{\hspace{0.17em}}\text{\hspace{0.17em}}n\ge 0.$

So that

$\begin{array}{c}{u}_{0}={c}_{1,0}\left(t\right)+{c}_{2,0}\left(t\right)x+{L}_{xx}^{-1}\left({x}^{2}{t}^{2}-x\right)\\ ={c}_{1,0}\left(t\right)+{c}_{2,0}\left(t\right)x+\frac{1}{12}{x}^{4}{t}^{2}-\frac{1}{6}{x}^{3},\end{array}$

The first approximant is ${\phi}_{1}={u}_{0}$ Applying the x conditions to ${\phi}_{1}$ it is clear that ${c}_{1,0}\left(t\right)=0$ and ${c}_{2,0}\left(t\right)=t$ . Thus, if the one-term approximant ${\phi}_{1}$ were

sufficient, the “solution” will be ${\phi}_{1}={u}_{0}=xt+\frac{1}{12}{x}^{4}{t}^{2}-\frac{1}{6}{x}^{3}$ . We set ${u}_{0}=f$ ,

then we divide it into two parts, so that

${u}_{0}={f}_{1}=xt,$

$\begin{array}{c}{u}_{1}={f}_{2}+{c}_{1,1}\left(t\right)+{c}_{2,1}\left(t\right)x+{L}_{xx}^{-1}{L}_{t}{u}_{0}-{L}_{xx}^{-1}{A}_{0}\\ =\left(\frac{1}{12}{x}^{4}{t}^{2}-\frac{1}{6}{x}^{3}\right)+{c}_{1,1}\left(t\right)+{c}_{2,1}\left(t\right)x+{L}_{xx}^{-1}{L}_{t}{u}_{0}-{L}_{xx}^{-1}{A}_{0}.\end{array}$

Then ${\phi}_{2}={u}_{0}+{u}_{1}$ is given by:

${\phi}_{2}=xt+\frac{1}{12}{x}^{4}{t}^{2}-\frac{1}{6}{x}^{3}+{c}_{1,1}\left(t\right)+{c}_{2,1}\left(t\right)x+{L}_{xx}^{-1}{L}_{t}{u}_{0}-{L}_{xx}^{-1}{A}_{0}$

Applying the condition at $x=0$ , we have ${c}_{1,1}(t)=0$ . From the condition on x at 0, we get ${c}_{2,1}\left(t\right)=0$ , thus ${u}_{1}=0$ , and ${u}_{n+1}=0,n\ge 1$ . Thus, the approximate solution is $u={u}_{0}=xt$ , which is the exact solution.

4. 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, Neumann or mixed boundary conditions. In [30] the authors defined the operator ${L}_{t}^{-1}$ and as definite integrals given by

${L}_{t}^{-1}={\displaystyle {\int}_{0}^{t}\text{d}{t}^{\prime}},\text{\hspace{0.17em}}\text{\hspace{0.17em}}{L}_{x}^{-1}={\displaystyle {\int}_{{x}_{0}}^{x}\text{d}{x}^{\prime}}{\displaystyle {\int}_{1}^{{x}^{\prime}}\text{d}{x}^{\u2033}}$ (4.1)

to solve the linear homogeneous heat equation ${u}_{t}={u}_{xx},\text{\hspace{0.17em}}{x}_{0}<x<1,\text{\hspace{0.17em}}t>0$ , subject to the mixed boundary conditions $u\left({x}_{0},t\right)={h}_{1}\left(t\right),{u}_{x}\left(1,t\right)={h}_{2}\left(t\right)$ .

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.

4.1. Boundary Value Problems

Consider Equations. (2.2)-(2.3), and applying the inverse operator ${L}_{xx}^{-1}$ defined by

${L}_{xx}^{-1}={\displaystyle {\int}_{a}^{x}\text{d}{x}^{\prime}}{\displaystyle {\int}_{b}^{{x}^{\prime}}\text{d}{x}^{\u2033}}$ to both sides of Equation (2.4), gives

$u\left(x,t\right)=u\left(a,t\right)+\left(x-a\right)\frac{\partial u\left(b,t\right)}{\partial x}+{L}_{xx}^{-1}g\left(x,t\right)-{L}_{xx}^{-1}{L}_{tt}u\left(x,t\right)-{L}_{xx}^{-1}Nu\left(x,t\right)$ (4.2)

i.e. the boundary conditions can be used directly to solve the boundary value problem in x-direction.

Substituting $u\left(x,t\right)={\displaystyle \sum _{n=0}^{\infty}{u}_{n}\left(x,t\right)}$ , and $Nu\left(x,t\right)={\displaystyle \sum _{n=0}^{\infty}{A}_{n}\left({u}_{0},{u}_{1},{u}_{2},\cdots ,{u}_{n}\right)}$ into Equation (4.2), gives

$\begin{array}{c}{\displaystyle \sum _{n=0}^{\infty}{u}_{n}\left(x,t\right)}=u\left(a,t\right)+\left(x-a\right)\frac{\partial u\left(b,t\right)}{\partial x}+{L}_{xx}^{-1}g\left(x,t\right)\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}-{L}_{xx}^{-1}{L}_{tt}{\displaystyle \sum _{n=0}^{\infty}{u}_{n}\left(x,t\right)}-{L}_{xx}^{-1}{\displaystyle \sum _{n=0}^{\infty}{A}_{n}\left(x,t\right)}\end{array}$

This yields the recursive relations

${u}_{0}=u\left(a,t\right)+\left(x-a\right)\frac{\partial u\left(b,t\right)}{\partial x}+{L}_{xx}^{-1}g\left(x,t\right),$

${u}_{n+1}=-{L}_{xx}^{-1}{L}_{tt}{u}_{n}-{L}_{xx}^{-1}{A}_{n},\text{\hspace{0.17em}}n\ge 0.$ (4.3)

4.2. 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 partial solution as in Equation (4.2) to get

$\begin{array}{c}u\left(x,t\right)=\frac{1}{2}(u\left(x,0\right)+{u}_{t}\left(x,0\right)t+u\left(a,t\right)+\left(x-a\right)\frac{\partial u\left(b,t\right)}{\partial x}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}+{L}_{tt}^{-1}g\left(x,t\right)+{L}_{xx}^{-1}g\left(x,t\right)-{L}_{tt}^{-1}{L}_{xx}u\left(x,t\right)-{L}_{xx}^{-1}{L}_{tt}u\left(x,t\right)\\ \text{\hspace{0.05em}}\text{\hspace{0.05em}}\begin{array}{c}\text{\hspace{0.05em}}\\ \text{\hspace{0.05em}}\end{array}-{L}_{tt}^{-1}Nu\left(x,t\right)-{L}_{xx}^{-1}Nu\left(x,t\right))\end{array}$ (4.4)

Substituting $u\left(x,t\right)={\displaystyle \sum _{n=0}^{\infty}{u}_{n}\left(x,t\right)}$ , and $Nu\left(x,t\right)={\displaystyle \sum _{n=0}^{\infty}{A}_{n}\left({u}_{0},{u}_{1},{u}_{2},\cdots ,{u}_{n}\right)}$ into Equation (4.3), gives

$\begin{array}{c}{\displaystyle \sum _{n=0}^{\infty}{u}_{n}\left(x,t\right)}=\frac{1}{2}(u\left(x,0\right)+{u}_{t}\left(x,0\right)t+u\left(a,t\right)+\left(x-a\right)\frac{\partial u\left(b,t\right)}{\partial x}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}+{L}_{tt}^{-1}g\left(x,t\right)+{L}_{xx}^{-1}g\left(x,t\right)-{L}_{tt}^{-1}{L}_{xx}{\displaystyle \sum _{n=0}^{\infty}{u}_{n}\left(x,t\right)}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}-{L}_{xx}^{-1}{L}_{tt}{\displaystyle \sum _{n=0}^{\infty}{u}_{n}\left(x,t\right)}-{L}_{tt}^{-1}{\displaystyle \sum _{n=0}^{\infty}{A}_{n}\left(x,t\right)}-{L}_{xx}^{-1}{\displaystyle \sum _{n=0}^{\infty}{A}_{n}\left(x,t\right)})\end{array}$

so that the recurrence relations are

${u}_{0}=\frac{1}{2}\left(u\left(x,0\right)+{u}_{t}\left(x,0\right)t+u\left(a,t\right)+\left(x-a\right)\frac{\partial u\left(b,t\right)}{\partial x}+{L}_{tt}^{-1}g\left(x,t\right)+{L}_{xx}^{-1}g\left(x,t\right)\right)$

${u}_{n+1}=\frac{1}{2}\left(-{L}_{tt}^{-1}{L}_{xx}{u}_{n}-{L}_{xx}^{-1}{L}_{tt}{u}_{n}-{L}_{tt}^{-1}{A}_{n}-{L}_{xx}^{-1}{A}_{n}\right),\text{\hspace{0.17em}}\text{\hspace{0.17em}}n\ge 0.$ (4.5)

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.

4.3. Numerical Experiments

Example 4

Consider the linear homogeneous heat equation

${u}_{t}-{u}_{xx}=0,\text{\hspace{0.17em}}\text{\hspace{0.17em}}0\le x\le 1,\text{\hspace{0.17em}}\text{\hspace{0.17em}}t>0,$

with specified conditions $u\left(x,0\right)={\text{e}}^{x},\text{\hspace{0.17em}}u\left(0,t\right)={\text{e}}^{t},\text{\hspace{0.17em}}{u}_{x}\left(1,t\right)={\text{e}}^{1+t}$ .

Rewriting the heat equation in the operator form as ${L}_{t}u={L}_{xx}u$ . Applying the

inverse operator ${L}_{xx}^{-1}$ defined by ${L}_{xx}^{-1}={\displaystyle {\int}_{0}^{x}\text{d}{x}^{\prime}}{\displaystyle {\int}_{1}^{{x}^{\prime}}\text{d}{x}^{\u2033}}$ , we get the recursive relations

${u}_{0}=u\left(0,t\right)+x\frac{\partial u\left(1,t\right)}{\partial x},$

${u}_{n+1}={L}_{xx}^{-1}{L}_{t}{u}_{n},\text{\hspace{0.17em}}\text{\hspace{0.17em}}n\ge 0.$

So that

${u}_{0}={\text{e}}^{t}+x{\text{e}}^{1+t},$

${u}_{1}=\frac{1}{6}{\text{e}}^{1+t}{x}^{3}+\frac{1}{2}{\text{e}}^{t}{x}^{2}-{\text{e}}^{t}x-\frac{1}{2}{\text{e}}^{1+t}x,$

${u}_{2}=\frac{1}{120}{\text{e}}^{1+t}{x}^{5}+\frac{1}{24}{\text{e}}^{t}{x}^{4}+\frac{1}{3}\left(\frac{-1}{2}{\text{e}}^{t}-\frac{1}{4}{\text{e}}^{1+t}\right){x}^{3}+\frac{5}{24}x{\text{e}}^{1+t}+\frac{1}{3}{\text{e}}^{t}x,$

Continuing in a similar way ${u}_{2},{u}_{3},\cdots ,{u}_{n}$ are obtained for some n, then we get the approximate solution ${u}_{ap}={\displaystyle \sum _{n=0}^{\infty}{u}_{n}\left(x,t\right)}$ which converged to the exact solution ${u}_{ex}={\text{e}}^{x+t}$ .

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 $t=0.5,t=1.0$ .

Example 5

Consider the linear homogeneous heat equation

${u}_{t}={u}_{xx}+u,\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}t>0,$

with specified conditions $u\left(0,t\right)={\text{e}}^{\left(1-{\text{\pi}}^{2}\right)t},\text{\hspace{0.17em}}{u}_{x}\left(0,t\right)=0$ .

Rewriting the heat equation in the operator form a ${L}_{t}u={L}_{xx}u+u$ .

Applying the inverse operator ${L}_{xx}^{-1}$ defined by ${L}_{xx}^{-1}={\displaystyle {\int}_{0}^{x}\text{d}{x}^{\prime}}{\displaystyle {\int}_{1}^{{x}^{\prime}}\text{d}{x}^{\u2033}}$ , we get the recursive relations

${u}_{0}=u\left(0,t\right)+x\frac{\partial u\left(0,t\right)}{\partial x},$

${u}_{n+1}={L}_{xx}^{-1}{L}_{t}{u}_{n}-{L}_{xx}^{-1}{u}_{n},\text{\hspace{0.17em}}n\ge 0.$

So that

${u}_{0}={\text{e}}^{\left(1-{\text{\pi}}^{2}\right)t},$

${u}_{1}=\frac{1}{2}\left(\left(1-{\text{\pi}}^{2}\right){\text{e}}^{\left(1-{\text{\pi}}^{2}\right)t}-{\text{e}}^{\left(1-{\text{\pi}}^{2}\right)t}\right){x}^{2},$

${u}_{2}=\frac{1}{4}\left(\frac{1}{6}{\left(1-{\text{\pi}}^{2}\right)}^{2}{\text{e}}^{\left(1-{\text{\pi}}^{2}\right)t}-\frac{1}{3}\left(1-{\text{\pi}}^{2}\right){\text{e}}^{\left(1-{\text{\pi}}^{2}\right)t}+\frac{1}{6}{\text{e}}^{\left(1-{\text{\pi}}^{2}\right)t}\right){x}^{4},$

Continuing in a similar way ${u}_{3},{u}_{4},\cdots ,{u}_{n}$ are obtained for some n, then we get the approximate solution ${u}_{ap}={\displaystyle \sum _{n=0}^{\infty}{u}_{n}\left(x,t\right)}$ which converged to the exact solution ${u}_{ex}={\text{e}}^{\left(1-{\text{\pi}}^{2}\right)t}\mathrm{cos}\left(\text{\pi}x\right)$ .

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 $t=0.5,t=1.0$ .

Example 6

Consider the nonlinear inhomogeneous wave equation

${u}_{tt}-{u}_{xx}+{u}^{2}=2{x}^{2}-2{t}^{2}+{x}^{4}{t}^{4},\text{\hspace{0.17em}}\text{\hspace{0.17em}}t>0,$

with specified initial condition $u\left(x,0\right)=0,{u}_{t}\left(x,0\right)=0$ and the boundary conditions $u\left(0,t\right)=0,{u}_{x}\left(0,t\right)=0$ .

Rewriting the wave equation in the operator form as

${L}_{tt}u-{L}_{xx}u+{u}^{2}=2{x}^{2}-2{t}^{2}+{x}^{4}{t}^{4}.$

Table 3. Absolute errors of Lesnic’s approach at $t=0.5,t=1.0$ .

Table 4. Absolute errors of Lesnic’s approach at $t=0.5,t=1.0$ .

Figure 5. The exact solution and the approximate solution of Lesnic’s approach for ${\phi}_{3},{\phi}_{5},{\phi}_{10}$ at $t=0.5$ .

Figure 6. The exact solution and the approximate solution of Lesnic’s approach for ${\phi}_{3},{\phi}_{5},{\phi}_{10}$ at $t=1.0$ .

Figure 7. The exact solution and the approximate solution of Lesnic’s approach for ${\phi}_{3},{\phi}_{5},{\phi}_{10}$ at $t=0.5$ .

Figure 8. The exact solution and the approximate solution of Lesnic’s approach for ${\phi}_{3},{\phi}_{5},{\phi}_{10}$ at $t=1.0$ .

Firstly, we consider the t partial solution as

${L}_{tt}u={L}_{xx}u-{u}^{2}+2{x}^{2}-2{t}^{2}+{x}^{4}{t}^{4}.$

Applying the inverse operator ${L}_{tt}^{-1}$ defined by ${L}_{tt}^{-1}={\displaystyle {\int}_{0}^{t}{\displaystyle {\int}_{0}^{t}\text{d}t\text{d}t}}$ , gives the recursive relations

${u}_{0}=u\left(x,0\right)+t{u}_{t}\left(x,0\right)+{L}_{tt}^{-1}\left(2{x}^{2}-2{t}^{2}+{x}^{4}{t}^{4}\right),$

${u}_{n+1}={L}_{tt}^{-1}{L}_{xx}{u}_{n}-{L}_{tt}^{-1}{A}_{n},\text{\hspace{0.17em}}\text{\hspace{0.17em}}n\ge 0.$

So that

${u}_{0}=\frac{1}{30}{x}^{4}{t}^{6}-\frac{1}{6}{t}^{4}+{x}^{2}{t}^{2},$

$\begin{array}{c}{u}_{1}=\frac{11}{840}{x}^{2}{t}^{8}+\frac{1}{6}{t}^{4}-\frac{1}{163800}{x}^{8}{t}^{14}+\frac{1}{11880}{x}^{4}{t}^{12}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}-\frac{1}{10}\left(\frac{1}{135}{x}^{6}+\frac{1}{324}\right){t}^{10}-\frac{1}{30}{x}^{4}{t}^{6},\end{array}$

Secondly, we consider the x partial solution

${L}_{xx}u={L}_{tt}u+{u}^{2}-2{x}^{2}+2{t}^{2}-{x}^{4}{t}^{4}.$

Applying the inverse operator ${L}_{xx}^{-1}$ defined by ${L}_{xx}^{-1}={\displaystyle {\int}_{0}^{x}\text{d}{x}^{\prime}}{\displaystyle {\int}_{0}^{{x}^{\prime}}\text{d}{x}^{\u2033}}$ , gives the recursive relations

${u}_{0}=u\left(0,t\right)+x\frac{\partial u\left(0,t\right)}{\partial x}+{L}_{xx}^{-1}\left(-2{x}^{2}+2{t}^{2}-{x}^{4}{t}^{4}\right),$

${u}_{n+1}={L}_{xx}^{-1}{L}_{t}{u}_{n}-{L}_{xx}^{-1}{A}_{n},\text{\hspace{0.17em}}\text{\hspace{0.17em}}n\ge 0.$

So that

${u}_{0}=-\frac{1}{30}{t}^{4}{x}^{6}-\frac{1}{6}{x}^{4}+{x}^{2}{t}^{2},$

$\begin{array}{c}{u}_{1}=-\frac{11}{840}{t}^{2}{x}^{8}+\frac{1}{6}{x}^{4}+\frac{1}{163800}{t}^{8}{x}^{14}+\frac{1}{11880}{t}^{4}{x}^{12}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}+\frac{1}{10}\left(\frac{-1}{135}{t}^{6}+\frac{1}{324}\right){x}^{10}+\frac{1}{30}{t}^{4}{x}^{6},\end{array}$

Next, we average the partial solutions, i.e. add two partial solutions and divide by two, so we obtain

${u}_{0}=\frac{1}{2}\left(\frac{1}{30}{x}^{4}{t}^{6}-\frac{1}{6}{t}^{4}+2{x}^{2}{t}^{2}-\frac{1}{30}{t}^{4}{x}^{6}-\frac{1}{6}{x}^{4}\right),$

$\begin{array}{c}{u}_{1}=\frac{1}{2}(\frac{11}{840}{x}^{2}{t}^{8}+\frac{1}{6}{t}^{4}-\frac{1}{163800}{x}^{8}{t}^{14}+\frac{1}{11880}{x}^{4}{t}^{12}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}-\frac{1}{10}\left(\frac{1}{135}{x}^{6}+\frac{1}{324}\right){t}^{10}-\frac{1}{30}{x}^{4}{t}^{6}-\frac{11}{840}{t}^{2}{x}^{8}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}+\frac{1}{6}{x}^{4}+\frac{1}{163800}{t}^{8}{x}^{14}+\frac{1}{11880}{t}^{4}{x}^{12}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}+\frac{1}{10}\left(\frac{-1}{135}{t}^{6}+\frac{1}{324}\right){x}^{10}+\frac{1}{30}{t}^{4}{x}^{6}),\end{array}$

Continuing in a similar way ${u}_{2},{u}_{3},\cdots ,{u}_{n}$ are obtained for some n, then we get the approximate solution ${u}_{ap}={\displaystyle \sum _{n=0}^{\infty}{u}_{n}\left(x,t\right)}$ which converged to the exact solution ${u}_{ex}={x}^{2}{t}^{2}$ .

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 5. Absolute errors of Lesnic’s approach at $t=0.5,t=1.0$ .

Figure 9. The exact solution and the approximate solution of Lesnic’s approach for ${\phi}_{3},{\phi}_{5},{\phi}_{10}$ at $t=0.5$ .

Figure 10. The exact solution and the approximate solution of Lesnic’s approach for ${\phi}_{3},{\phi}_{5},{\phi}_{10}$ at $t=1.0$ .

results for the exact solution and the approximate solution by application of Lesnic’s approach for $t=0.5\text{\hspace{0.17em}},\text{\hspace{0.17em}}t=\mathrm{1.0.}$

Example 7

Consider the nonlinear inhomogeneous heat equation

${u}_{t}={u}_{xx}+{u}^{2}-u-{t}^{2}{\text{e}}^{2x}+{\text{e}}^{x},\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}t>0,$

with specified conditions $u\left(x,0\right)=0,u\left(0,t\right)=t,{u}_{x}\left(0,t\right)=t$ . Rewriting the heat equation in the operator form as

${L}_{t}u={L}_{xx}u+{u}^{2}-u-{t}^{2}{\text{e}}^{2x}+{\text{e}}^{x}.$

Applying the inverse operator ${L}_{xx}^{-1}$ defined by ${L}_{xx}^{-1}={\displaystyle {\int}_{0}^{x}\text{d}{x}^{\prime}}{\displaystyle {\int}_{0}^{{x}^{\prime}}\text{d}{x}^{\u2033}}$ , we get the recursive relations

${u}_{0}=u\left(0,t\right)+x\frac{\partial u\left(0,t\right)}{\partial x}+{L}_{xx}^{-1}\left({t}^{2}{\text{e}}^{2x}-{\text{e}}^{x}\right),$

${u}_{n+1}={L}_{xx}^{-1}{L}_{t}{u}_{n}+{L}_{xx}^{-1}{u}_{n}-{L}_{xx}^{-1}{A}_{n},\text{\hspace{0.17em}}\text{\hspace{0.17em}}n\ge 0.$

So that

${u}_{0}=t+xt+1-\frac{1}{4}{t}^{2}-\frac{1}{2}{t}^{2}x+x+\frac{1}{4}{t}^{2}{\text{e}}^{2x}-{\text{e}}^{x},$

We set ${u}_{0}=f$ then we divide it into two parts, so that

${u}_{0}={f}_{1}=t+xt,$

$\begin{array}{c}{u}_{1}={f}_{2}+{L}_{xx}^{-1}{L}_{t}{u}_{0}+{L}_{xx}^{-1}{u}_{0}-{L}_{xx}^{-1}{A}_{0}\\ =\left(1-\frac{1}{4}{t}^{2}-\frac{1}{2}{t}^{2}x+x+\frac{1}{4}{t}^{2}{\text{e}}^{2x}-{\text{e}}^{x}\right)+\frac{1}{6}{x}^{3}+\frac{1}{2}{x}^{2}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}-\frac{1}{12}{t}^{2}{x}^{4}-\frac{1}{3}{t}^{2}{x}^{3}-\frac{1}{2}{t}^{2}{x}^{2}+\frac{1}{6}t{x}^{3}+\frac{1}{2}t{x}^{2},\end{array}$

Continuing in a similar way ${u}_{2},{u}_{3},\cdots ,{u}_{n}$ are obtained for some n, then we get the approximate solution ${u}_{ap}={\displaystyle \sum _{n=0}^{\infty}{u}_{n}\left(x,t\right)}$ which converged to the exact solution ${u}_{ex}=t{\text{e}}^{x}$ .

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 $t=0.5$ .

5. Conclusion

The main goal of this paper is the development of an efficient and accurate method 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.

Table 6. Absolute errors of Lesnic’s approach and modified of Lesnic’s Approach at $t=0.5$ .

Figure 11. The exact solution and the approximate solution of Lesnic’s approach for ${\phi}_{3},{\phi}_{5},{\phi}_{10}$ at $t=0.5$ .

Figure 12. The exact solution and the approximate solution of Modified of Lesnic’s approach for ${\phi}_{3},{\phi}_{5},{\phi}_{10}$ at $t=0.5$ .

The accuracy and computational efficiency of the proposed method are verified by numerical examples. Thus, the method is highly recommended for varieties of mixed boundary value problems.

Cite this paper

Alzaid, N.A. and Bakodah, H.O. (2018) Numerical Treatment of Initial-Boundary Value Problems with Mixed Boundary Conditions. American Journal of Computational Mathematics, 8, 153-174. https://doi.org/10.4236/ajcm.2018.82012

References

- 1. Gedye, R., Smith, F., Westaway, K., Ali, H., Baldisera, L., Laberge, L. and Rousell, J. (1986) The Use of Microwave Ovens for Rapid Organic Synthesis. Tetrahedron Letters, 27, 279-282. http://dx.doi.org/10.1016/S0040-4039(00)83996-9
- 2. Giguere, R.J., Bray, T.L. and Duncan, S.M. (1986) Application of Commercial Microwave Oven to Organic Synthesis. Tetrahedron Letters, 27, 4945-4948. http://dx.doi.org/10.1016/S0040-4039(00)85103-5
- 3. Jin, Q.H. (1999) Microwave Chemistry. Science Press, Beijing.
- 4. Kappe, C.O. (2004) Controlled Microwave Heating in Modern Organic Synthesis. Angewandte Chemie International Edition, 43, 6250-6284. http://dx.doi.org/10.1002/anie.200400655
- 5. (a) Pradhan, D.K., Dharamrajan, T.S., Mishra, M.R. and Mishra, A. (2012) An Overview of Microwave Oven in the Field of Synthetic Chemistry. International Journal of Research and Development in Pharmacy & Life Sciences (IJRDPL), 44, 44-50. (b) Taylor, M., Atri, S.S. and Minhas, S. Evalueserve Analysis: Developments in Microwave Chemistry. http://www.rsc.org/images/evalueserve_tcm18-16758.pdf (c) Horeis, A.G., Pichler, S., Stadler, A., Gössler, W. and Kappe, C.O. (2001) Microwave-Assisted Organic Synthesis—Back to the Roots. 5th International Electronic Conference on Synthetic Organic Chemistry (ECSOC-5), Basel, 1-30 September 2001. http://www.mdpi.org/ecsoc-5.htm
- 6. http://shodhganga.inflibnet.ac.in/bitstream/10603/3197/8/08_chapter%201.pdf
- 7. The type of radiation emitted by MW ovens is non-ionizing. This means that unlike X-ray or UV light, it does not have potential to cause cancer. Although it possess possible heat-burn risks. Dr. Spencer himself, despite being literally surrounded by intense microwaves for a long period of his life, lived to the ripe old age of 76, dying apparently of natural causes. (a) Stuerga, D. and Delmotte, M. (2002) Chapter 1. Wave-Material Interactions, Microwave Technology and Equipment. In: Loupy, A., Ed., Microwaves in Organic Synthesis, Wiley-VCH, Weinheim, 1-34. (b) Mingos, M.D.P. (2004) In: Lidstrom, P. and Tierney, J.P., Eds., Microwave-Assisted Organic Synthesis, Blackwell, Oxford. (c) Gabriel, C., Gabriel, S., Grant, E.H., Halstead, B.S. and Mingos, D.M.P. (1998) Dielectric Parameters Relevant to Microwave Dielectric Heating. Chemical Society Reviews, 27, 213-223. (d) Kerr, J.A. (1992) Strengths of Chemical Bonds. In: Lide, D.R., Ed., CRC Handbook of Chemistry and Physics, 76th Edition, CRC Press, Boca Raton, Ann Arbor, London, Tokyo, 51. (e) Atkins, P.W. (1990) Physical Chemistry. Oxford University Press, Oxford, 938. (f) Stuerga, D.A.C. and Gaillard, P. (1996) Microwave Athermal Effects in Chemistry: A Myth’s Autopsy Part I: Historical Background and Fundamentals of Wave-Matter Interaction. The Journal of Microwave Power Electromagnetic Energy, 31, 87-100.
- 8. In specific cases, magnetic field interactions have also been observed, see: (a) Stass, D.V., Woodward, J.R., Timmel, C.R., Hore, P.J. and McLauchlan, K.A. (2000) Radiofrequency Magnetic Field Effects on Chemical Reaction Yields. Chemical Physics Letters, 329, 15-22. http://dx.doi.org/10.1016/S0009-2614(00)00980-5 (b) Timmel, C.R. and Hore, P.J. (1996) Oscillating Magnetic Field Effects on the Yields of Radical Pair Reactions. Chemical Physics Letters, 257, 401-408. http://dx.doi.org/10.1016/0009-2614(96)00466-6 (c) Woodward, J.R., Jackson, R.J., Timmel, C.R., Hore, P.J. and McLauchlan, K.A. (1997) Resonant Radiofrequency Magnetic Field Effects on a Chemical Reaction. Chemical Physics Letters, 272, 376-382. http://dx.doi.org/10.1016/S0009-2614(97)00542-3
- 9. Hayes, B.L. (2002) Microwave Synthesis: Chemistry at the Speed of Light. CEM Publishing, Matthews NC and References Cited Therein.
- 10. Gaba, M. and Dhingra, N. (2011) Microwave Chemistry: General Features and Applications. Indian Journal of Pharmaceutical Education and Research, 45, 175-183.
- 11. http://www.milestonesci.com
- 12. http://www.maos.net
- 13. http://www.cem.com
- 14. Westaway, K.C. and Gedye, R.J. (1995) The Question of Specific Activation of Organic Reactions by Microwaves. Journal of Microwave Power and Electromagnetic Energy, 30, 219-230.
- 15. Kuhnert, N. (2002) Microwave-Assisted Reactions in Organic Synthesis—Are There Any Nonthermal Microwave Effects? Angewandte Chemie International Edition, 41, 1863-1866. http://dx.doi.org/10.1002/1521-3773(20020603)41:11<1863::AID-ANIE1863>3.0.CO;2-L
- 16. Langa, F., De la Cruz, P., De la Hoz, A., Diaz-Ortiz, A. and Diez-Barra, E. (1997) Microwave Irradiation: More than Just a Method for Accelerating Reactions. Contemporary Organic Synthesis, 4, 373-386. http://dx.doi.org/10.1039/co9970400373
- 17. Saillard, R., Poux, M., Berlan, J. and Audhvy-Peaudecert, M. (1995) Microwave Heating of Organic Solvents: Thermal Effects and Field Modeling. Tetrahedron, 51, 4033-4042. http://dx.doi.org/10.1016/0040-4020(95)00144-W
- 18. Stuerga, D. and Gaillard, P. (1996) Microwave Heating as a New Way to Induce Localized Enhancements of Reaction Rate. Non-Isothermal and Heterogeneous Kinetics. Tetrahedron, 52, 5505-5510. http://dx.doi.org/10.1016/0040-4020(96)00241-4
- 19. Caddick, S. (1995) Microwave Assisted Organic Reactions. Tetrahedron, 51, 10403-10432. http://dx.doi.org/10.1016/0040-4020(95)00662-R
- 20. Strauss, C.R. and Trainor, R.W. (1995) Developments in Microwave-Assisted Organic Chemistry. Australian Journal of Chemistry, 48, 1665-1692. http://dx.doi.org/10.1071/CH9951665
- 21. Bose, A.K., Banik, B.K., Lavlinskaia, N., Jayaraman, M. and Manhas, M.S. (1997) MORE Chemistry in a Microwave. Chemtech, 27, 18-124.
- 22. Lidstrom, P., Tierney, J., Wathey, B. and Westman, J. (2001) Microwave Assisted Organic Synthesis—A Review. Tetrahedron, 57, 9225-9283. http://dx.doi.org/10.1016/S0040-4020(01)00906-1
- 23. Nuchter, M., Ondruschka, B. and Lautenschlager, W. (2003) Microwave-Assisted Chemical Reactions. Chemical Engineering & Technology, 26, 1208-1216. http://dx.doi.org/10.1002/ceat.200301836
- 24. Nuchter, M., Ondruschka, B., Bonrath, W. and Gum, A. (2004) Microwave Assisted Synthesis—A Critical Technology Overview. Green Chemistry, 6, 128-141. http://dx.doi.org/10.1039/b310502d
- 25. Abramovitch, R.A. (1991) Applications of Microwave Energy in Organic Chemistry. A Review. Organic Preparations and Procedures International, 23, 683-711. http://dx.doi.org/10.1080/00304949109458244
- 26. Leadbeater, N.E. (2010) In Situ Reaction Monitoring Of Microwave-Mediated Reactions Using IR Spectroscopy. Chemical Communications, 46, 6693-6695. http://dx.doi.org/10.1039/c0cc01921f
- 27. (a)http://www.mrl.ucsb.edu/polymer-characterization-facility/instruments/microwave-reactor-biotage (b)http://www.biotagepathfinder.com/pressurePC.jsp
- 28. For recent reviews on microwave-assisted organic synthesis see: (a) Adam, D. (2003) Microwave Chemistry: Out of the Kitchen. Nature, 421, 571-572. http://dx.doi.org/10.1038/421571a (b) Kappe, C.O. and Dallinger, D. (2006) The Impact of Microwave Synthesis on Drug Discovery. Nature Reviews Drug Discovery, 5, 51-63. http://dx.doi.org/10.1038/nrd1926 (c) Nagariya, A.K., Meena, A.K., Kiran, K., Yadav, A.K., Niranjan, U.S., Pathak, A.K., Singh, B. and Rao, M.M. (2010) Microwave Assisted Organic Reaction as New Tool in Organic Synthesis. Journal of Pharmacy Research, 3, 575-580. (d) Jignasa, K.S., Ketan, T.S., Bhumika, S.P. and Anuradha, K.G. (2010) Microwave Assisted Organic Synthesis: An Alternative Synthetic Strategy. Der Pharma Chemica, 2, 342-353. (e) Bose, A.K., Banik, B.K., Lavlinskaia, N., Jayaraman, M. and Manhas, M.S. (1997) MORE Chemistry in a Microwave. Chemtech, 27, 18-24. (f) Bose, A.K., Manhas, M.S., Ganguly, S.N., Sharma, A.H. and Banik, B.K. (2002) MORE Chemistry for Less Pollution: Applications for Process Development. Synthesis, 11, 1578-1591. http://dx.doi.org/10.1055/s-2002-33344
- 29. Moghaddam, F.M. and Sharifi, A. (1995) Microwave Promoted Acetalization of Aldehydes and Ketones. Synthetic Communications, 25, 2457-2461. http://dx.doi.org/10.1080/00397919508015450
- 30. Mogilaiah, K., Kavita, S. and Babu, H.R. (2003) Microwave Assisted Addition-Elimination Reactions of 3-Phenylimino-2-Indolinones with 2-Hydrazino-3-(p-Chlorophenyl)-1, 8-Naphthyridine: A Simple and Facile Synthesis of 3-(3-p-Chlorophenyl-l,8-Naphthyridin-2-yl Hydrazono)-2-Indolinones. Indian Journal of Chemistry, 42B, 1750-1752.
- 31. Abramovitch, R.A., Shi, Q. and Bogdal, D. (1995) Microwave-Assisted Alkylations of Activated Methylene Groups. Synthetic Communications, 25, 1-8. http://dx.doi.org/10.1080/00397919508010781
- 32. Miljanic, O.S., Volhardt, K.P.C. and Whitener, G.D. (2003) An Alkyne Metathesis-Based Route to ortho-Dehydrobenzannulenes. Synlett, 1, 29-34.
- 33. Motorina, I.A., Parly, F. and Grierson, D.S. (1996) Selective O-Allylation of Amidoalcohols on Solid Support. Synlett, 4, 389-391. http://dx.doi.org/10.1055/s-1996-5416
- 34. Mccarroll, A.J., Sandham, D.A., Titumb, L.R., Dek Lewis, A.K., Cloke, F.G.N., Davies, B.P., Desantand, A.P., Hiller, W. and Caddicks, S. (2003) Studies on High-Temperature Amination Reactions of Aromatic Chlorides Using Discrete Palladium-N-Heterocyclic Carbene (NHC) Complexes and in Situ Palladium/Imidazolium Salt Protocols. Molecular Diversity, 7, 115-123. http://dx.doi.org/10.1023/B:MODI.0000006863.14423.da
- 35. Angrish, C., Kumar, A. and Chouhan, S.M.S. (2005) Microwave Assisted Aromatic Nucleophilic Substitution Reaction of Chloronitrobenzenes with Amines in Ionic Liquids. Indian Journal of Chemistry, 44B, 1515-1518.
- 36. Ermolat’ev, D.S., Giménez, V.N., Babaev, E.V. and Van der Eycken, E. (2006) Efficient Pd(0)-Mediated Microwave-Assisted Arylation of 2-Substituted Imidazo[1,2-a]Pyrimidines. Journal of Combinatorial Chemistry, 8, 659-663.
- 37. Azarifar, D. and Ghasemnejad, H. (2003) Microwave-Assisted Synthesis of Some 3,5-Arylated-2-Pyrazolines. Molecules, 8, 642-648. http://dx.doi.org/10.3390/80800642
- 38. Al-Obeidi, F., Austin, R.E., Okonoya, J.F. and Bond, D.R.S. (2003) Microwave-Assisted Solid-Phase Synthesis (MASS): Parallel and Combinatorial Chemical Library Synthesis. Mini-Reviews in Medicinal Chemistry, 3, 449-460. http://dx.doi.org/10.2174/1389557033488042
- 39. Kim, J.K., Kwon, P.S., Kwon, T.W., Chung, S.K. and Lee, J.W. (1996) Application of Microwave Irradiation Techniques for the Knoevenagel Condensation. Synthetic Communications, 26, 535-542. http://dx.doi.org/10.1080/00397919608003646
- 40. Burten, G., Cao, P., Li, G. and Rivero, R. (2003) Palladium-Catalyzed Intermolecular Coupling of Aryl Chlorides and Sulfonamides under Microwave Irradiation. Organic Letters, 5, 4373-4376. http://dx.doi.org/10.1021/ol035655u
- 41. Arevela, R.K. and Leadbeater, N.E. (2003) Rapid, Easy Cyanation of Aryl Bromides and Chlorides Using Nickel Salts in Conjunction with Microwave Promotion. Journal of Organic Chemistry, 68, 9122-9125. http://dx.doi.org/10.1021/jo0350561
- 42. Crawford, K.R., Bur, S.K., Straub, C.S. and Padwa, A. (2003) Intramolecular Cyclization Reactions of 5-Halo- and 5-Nitro-Substituted Furans. Organic Letters, 5, 3337-3340. http://dx.doi.org/10.1021/ol035233k
- 43. Lerestif, J.M., Perocheav, J., Tonnard, F., Bazareav, J.P. and Hamelin, J. (1995) 1,3-Dipolar Cycloaddition of Imidate Ylides on Imino-Alcohols: Synthesis of New Imidazolones Using Solvent Free Conditions. Tetrahedron, 51, 6757-6774. http://dx.doi.org/10.1016/0040-4020(95)00321-X
- 44. Jayakumar, G., Ajithabai, M.D., Santhosh, B., Veena, C.S. and Nair, M.S. (2003) Microwave Assisted Acetylation and Deacetylation Studies on the Triterpenes Isolated from Dysoxylum malabaricum and Dysoxylum beddomei. Indian Journal of Chemistry, 42B, 429-431.
- 45. Calinescu, I., Calinescu, R., Martin, D.I. and Radoiv, M.T. (2003) Microwave-Enhanced Dechlorination of Chlorobenzene. Research on Chemical Intermediates, 29, 71-81. http://dx.doi.org/10.1163/156856703321328424
- 46. Mavoral, J.A., Cativicla, C., Garcia, J.I., Pires, E., Rovo, A.J. and Figueras, F. (1995) Diels-Alder Reactions of α-Amino Acid Precursors by Heterogeneous Catalysis: Thermal vs. Microwave Activation. Applied Catalysis A, 131, 159-166. http://dx.doi.org/10.1016/0926-860X(95)00131-X
- 47. Santagoda, V., Fiorino, F., Perissuti, E., Severino, B., Terracciano, S., Cirino, G. and Caliendo, G. (2003) A Convenient Strategy Of Dimerization By Microwave Heating And Using 2,5-Diketopiperazine As Scaffold. Tetrahedron Letters, 44, 1145-1148.
- 48. Mazo, P.C. and Ríos, L.A. (2010) Esterification and Transesterification Assisted by Microwaves of Crude Palm Oil: Homogeneous Catalysis. Latin American Applied Research, 40, 343-349.
- 49. Diaz-Ortiz, A., De la Hoz, A., Merrero, M.A., Prieto, P., Sanchez-Migallon, A., Cassio, F.P., Arriela, A., Vivanco, S. and Foces, C. (2003) Enhancing Stereochemical Diversity by Means of Microwave Irradiation in the Absence of Solvent: Synthesis of Highly Substituted Nitroproline Esters via 1,3-Dipolar Reactions. Molecular Diversity, 7, 175-180. http://dx.doi.org/10.1023/B:MODI.0000006799.01924.2e
- 50. Jnagaki, T., Fukuhara, T. and Hara, S. (2003) Effective Fluorination Reaction with Et3N.3HF under Microwave Irradiation. Synthesis, 8, 1157-1159.
- 51. Lechmann, F., Pilotti, A. and Luthman, K. (2003) Efficient Large Scale Microwave Assisted Mannich Reactions Using Substituted Acetophenones. Molecular Diversity, 7, 145-152. http://dx.doi.org/10.1023/B:MODI.0000006809.48284.ed
- 52. Kiasat, A.R., Kazemi, F. and Rafati, M. (2003) Microwave Promoted Rapid Oxidation of Alcohols Using Cobalt Nitrate Hexahydrate Supported on Silica Gel Under Solvent Free Conditions. Synthetic Communications, 33, 601-605. http://dx.doi.org/10.1081/SCC-120015814
- 53. Gospondinova, M., Gredard, A., Jeannin, M., Chitanv, G.C., Carpov, A., Thiery, V. and Besson, T. (2002) Efficient Solvent-Free Microwave Phosphorylation of Microcrystalline Cellulose. Green Chemistry, 4, 220-222. http://dx.doi.org/10.1039/b200295g
- 54. Vu, Z.T., Liu, L.J. and Zhuo, R.X. (2003) Microwave-Improved Polymerization of ζ-Caprolactone Initiated by Carboxylic Acids. Journal of Polymer Science Part A: Polymer Chemistry, 41, 13-21. http://dx.doi.org/10.1002/pola.10546
- 55. Srikrishana, A. and Kumar, P.P. (1995) Naphthalenes Microwave Irradiation Induced Rearrangement on Montmorillonite K-10. Tetrahedron Letters, 36, 6313-6316. http://dx.doi.org/10.1016/0040-4039(95)01219-8
- 56. Danks, T.N. and Wagner, G. (2005) Microwave-Assisted Reductions. In: Tierney, J.P. and Lidström, P., Eds., Microwave Assisted Organic Synthesis, Blackwell Publishing Ltd., Oxford. http://dx.doi.org/10.1002/9781444305548.ch4
- 57. Garbacia, S., Desai, B., Lavastre, O. and Kappe, C.O. (2003) Microwave-Assisted Ring-Closing Metathesis Revisited. On the Question of the Nonthermal Microwave Effect. Journal of Organic Chemistry, 68, 9136-9139. http://dx.doi.org/10.1021/jo035135c
- 58. Laurent, A., Jacquault, P., Di Martino, J.L. and Hamelin, J. (1995) Fast Synthesis of Amino Acid Salts and Lactams without Solvent under Microwave Irradiation. Chem. Commun., 11, 1101.
- 59. Nielsen, T.E. and Schreiber, S.L. (2007) Towards the Optimal Screening Collection: A Synthesis Strategy. Angewandte Chemie International Edition, 47, 48-56. http://dx.doi.org/10.1002/anie.200703073
- 60. (a) Hayes, B.L. (2004) Recent Advances in Microwave-Assisted Synthesis. Aldrichimica Acta, 37, 66-76. (b) Bogdal, D. (2005) Microwave-Assisted Organic Synthesis: One Hundred Reaction Procedures (Tetrahedron Organic Chemistry). Elsevier Science, Amsterdam.
- 61. (a) Seijas, J.A., Vázquez-Tato, M.P. and Carballido-Reboredo, R. (2005) Solvent-Free Synthesis of Functionalized Flavones under Microwave Irradiation. Journal of Organic Chemistry, 70, 2855-2858. (b) Ranu, B.C., Guchhait, S.K., Ghosh, K. and Patra, A. (2000) Construction of Bicyclo[2.2.2]octanone Systems by Microwave-Assisted Solid Phase Michael Addition Followed by Al2O3-Mediated Intramolecular Aldolisation. An Eco-Friendly Approach. Green Chemistry, 2, 5-6. http://dx.doi.org/10.1039/a907689a
- 62. Porcheddu, A., Mura, M.G., De Luca, L., Pizzetti, M. and Taddei, M. (2012) From Alcohols to Indoles: A Tandem Ru Catalyzed Hydrogen-Transfer Fischer Indole Synthesis. Organic Letters, 14, 6112-6115. http://dx.doi.org/10.1021/ol3030956
- 63. Miura, T., Biyajima, T., Fujii, T. and Murakami, M. (2012) Synthesis of α-Amino Ketones from Terminal Alkynes via Rhodium-Catalyzed Denitrogenative Hydration of N-Sulfonyl-1,2,3-Triazoles. Journal of the American Chemical Society, 134, 194-196. http://dx.doi.org/10.1021/ja2104203
- 64. García-Muñoz, Á., Ortega-Arizmendi, A.I., García-Carrillo, M.A., Díaz, E., Gonzalez-Rivas, N. and Cuevas-Yañez, E. (2012) Direct, Metal-Free Synthesis of Benzyl Alcohols and Deuterated Benzyl Alcohols from p-Toluenesulfonylhydrazones Using Water as Solvent. Synthesis, 44, 2237-2242. http://dx.doi.org/10.1055/s-0031-1290372
- 65. Yu, C.W., Chen, G.S., Huang, C.W. and Chern, J.W. (2012) Efficient Microwave-Assisted Pd-Catalyzed Hydroxylation of Aryl Chlorides in the Presence of Carbonate. Organic Letters, 14, 3688-3691. http://dx.doi.org/10.1021/ol301523q
- 66. Collados, J.F., Toledano, E., Guijarro, D. and Yus, M. (2012) Microwave-Assisted Solvent-Free Synthesis of Enantiomerically Pure N-(Tert-Butylsulfinyl)Imines. Journal of Organic Chemistry, 77, 5744-5750. http://dx.doi.org/10.1021/jo300919x
- 67. Nagaraj, M., Boominathan, M., Muthusubramanian, S. and Bhuvanesh, N. (2012) Microwave-Assisted Metal-Free Synthesis of 2,8-Diaryl-6-Aminoimidazo-[1,2-a]Pyridine via Amine-Triggered Benzannulation. Synlett, 23, 1353-1357. http://dx.doi.org/10.1055/s-0031-1290979
- 68. Yang, D., Burugupalli, S., Daniel, D. and Chen, Y. (2012) Microwave-Assisted One-Pot Synthesis of Isoquinolines, Furopyridines, and Thienopyridines by Palladium-Catalyzed Sequential Coupling-Imination-Annulation of 2-Bromoarylaldehydes with Terminal Acetylenes and Ammonium Acetate. Journal of Organic Chemistry, 77, 4466-4472. http://dx.doi.org/10.1021/jo300494a
- 69. Cantillo, D., Gutmann, B. and Kappe, C.O. (2011) Mechanistic Insights on Azide-Nitrile Cycloadditions: On the Dialkyltin Oxide-Trimethylsilyl Azide Route and a New Vilsmeier-Haack-Type Organocatalyst. Journal of the American Chemical Society, 133, 4465-4475. http://dx.doi.org/10.1021/ja109700b
- 70. Baghbanzadeh, M., Pilger, C. and Kappe, C.O. (2011) Palladium-Catalyzed Direct Arylation of Heteroaromatic Compounds: Improved Conditions Utilizing Controlled Microwave Heating. Journal of Organic Chemistry, 76, 8138-8142. http://dx.doi.org/10.1021/jo201516v
- 71. Revelant, G., Dunand, S., Hesse, S. and Kirsch, G. (2011) Microwave-Assisted Synthesis of 5-Substituted 2-Aminothiophenes Starting from Arylacetaldehydes. Synthesis, 2011, 2935-2940.
- 72. Kumar, A. and Rao, V.K. (2011) Microwave-Assisted and Yb(OTf)3-Promoted One-Pot Multicomponent Synthesis of Substituted Quinolines in Ionic Liquid. Synlett, 15, 2157-2162. http://dx.doi.org/10.1055/s-0030-1261200
- 73. Barnett, D.S. and Schaus, S.E. (2011) Asymmetric Propargylation of Ketones Using Allenylboronates Catalyzed by Chiral Biphenols. Organic Letters, 13, 4020-4023. http://dx.doi.org/10.1021/ol201535b
- 74. Rostami, M., Khosropour, A.R., Mirkhani, V., Mohammadpoor-Baltork, I., Moghadam, M. and Tangestaninejad, S. (2011) An Efficient, Simple, and Scaleable Domino Reaction to Diverse N-(1-Oxo-1H-inden-2-yl)benzamides Catalyzed by HPW at Nano-SiO2 under Microwave Irradiation. Synlett, 2011, 1677-1682.
- 75. Talukdar, D., Saikia, L. and Thakur, A.J. (2011) Zirconyl Chloride: An Efficient, Water-Tolerant, and Reusable Catalyst for the Synthesis of N-Methylamides. Synlett, 11, 1597-1601.
- 76. Allam, B.K. and Singh, K.N. (2011) An Efficient Phosphine-Free Heck Reaction in Water Using Pd(L-Proline)2 as the Catalyst Under Microwave Irradiation. Synthesis, 2011, 1125-1131.
- 77. Seijas, J.A., Vazquez-Tato, M.P., Martinez, M.M. and Corredoira, G.N. (1999) Direct Synthesis of Imides from Dicarboxylic Acids Using Microwaves. Journal of Chemical Research, 1999, 420-425.
- 78. Dahmani, Z., Rahmouni, M., Brugidou, R., Bazureau, J.P. and Hamelin, J. (1998) A New Route to α-Hetero β-Enamino Esters Using a Mild and Convenient Solvent-Free Process Assisted by Focused Microwave Irradiation. Tetrahedron Letters, 39, 8453-8456. http://dx.doi.org/10.1016/S0040-4039(98)01894-2
- 79. Bogdal, D., Pielichowski, J. and Jaskot, K. (1997) Remarkable Fast N-Alkylation of Azaheterocycles under Microwave Irradiation in Dry Media. HeteroCycles, 45, 715-722. http://dx.doi.org/10.3987/COM-96-7731
- 80. Azarifar, D. and Ghasemnejad, H. (2003) Microwave-Assisted Synthesis of Some 3,5-Arylated 2-Pyrazolines. Molecules, 8, 642-648. http://dx.doi.org/10.3390/80800642
- 81. Michael, C.J. and Dennis, M.P. (1988) Microwave Heating Apparatus for Laboratory Analyses. CEM Corporation. US Patent No. 4835354.
- 82. Price, W.H. and Earl, E.K. (1999) Pressure Vessel with Composite Sleeve. CEM Limited, L.L.C. US Patent No. 6534140.
- 83. Bernd, O., et al. (2002) Device for Performing Safety Functions in Areas with High Frequency Radiation. Mikrowellen-Labor-SystemeGmbh. US Patent No. 20020096340A1.
- 84. Edward, J.T., Price, W.H. and Earl Edward, K. (1999) Pressure Sensing Reaction Vessel for Microwave Assisted Chemistry. CEM Corporation. US Patent No. 6124582.
- 85. Price, W.H. (1999) Sealing Closure for High-Pressure Vessels in Microwave Assisted Chemistry. CEM Corporation. US Patent No. 6287526.
- 86. Stone-Elander, S. and Elander, N. (2002) Microwave Applications in Radiolabelling with Short-Lived Positron-Emitting Radionuclides. Journal of Labelled Compounds and Radiopharmaceuticals, 45, 715-746. http://dx.doi.org/10.1002/jlcr.593
- 87. Stadler, A., Yousefi, B.H., Dallinger, D., Walla, P., Van der Eycken, E., Kaval, N. and Kappe, C.O. (2003) Scalability of Microwave-Assisted Organic Synthesis. From Single-Mode to Multimode Parallel Batch from Single-Mode to Multimode Parallel Batch Reactors. Organic Process Research Development, 7, 707-716.
- 88. Moseley, J.D. and Kappe, C.O. (2011) A Critical Assessment of the Greenness and Energy Efficiency of Microwave-Assisted Organic Synthesis. Green Chemistry, 13, 794-806. http://dx.doi.org/10.1039/c0gc00823k
- 89. MILESTONES Inc. (2014) Digestion. http://www.milestonesci.com/digestion.php
- 90. MILESTONES Inc. (2014) Extraction. http://www.milestonesci.com/extraction.php
- 91. Developments in Microwave Chemistry (2005) Royal Society of Chemistry, Evalueserve, London, UK.
- 92. Kakita, Y., Kashige, N., Murata, K., Kuroiwa, A., Funatsu, M. and Watanabe, K. (1995) Inactivation of Lactobacillus Bacteriophage PL-1 by Microwave Irradiation. Microbiology and Immunology, 39, 571-576. http://dx.doi.org/10.1111/j.1348-0421.1995.tb02244.x
- 93. Kakita, Y., Funatso, M., Miake, F. and Watanabe, K. (1999) Effects of Microwave Irradiation on Bacteria Attached to the Hospiral White Coats. International Journal of Occupational Medicine and Environmental Health, 12, 123-126.
- 94. Sekhon, B.S. (2010) Microwave-Assisted Pharmaceutical Synthesis: An Overview. International Journal of PharmTech Research, 2, 827-833.
- 95. Nüchter, M., Ondruschka, B., Bonrath, W. and Gum, A. (2004) Microwave Assisted Synthesis—A Critical Technology Overview. Green Chemistry, 6, 128-141.
- 96. Brett, A. and Christopher, R. (2005) Toward Rapid, “Green”, Predictable Microwave-Assisted Synthesis. Accounts of Chemical Research, 38, 653-661. http://dx.doi.org/10.1021/ar040278m
- 97. Carslaw, H.S. and Jaeger, J.C. (1948) Conduction of Heat in Solids. Clarendon Press, Oxford.
- 98. Carrier, G.F., Krook, M. and Pearson, C.E. (1966) Functions of a Complex Variable. McGraw-Hill, New York.
- 99. Courant, R. and Hilbert, D. (1953) Methods of Mathematical Physics. Vol. 1, Interscience Publishers Inc., New York.
- 100. Sneddon, I.N. (1966) Mixed Boundary Value Problems in Potential Theory. Wiley, New York.
- 101. Tranter, C.J. (1951) Integral Transforms in Mathematical Physics. Wiley, New York.
- 102. Fabrikant, V.I. (1991) Mixed Boundary Value Problems of Potential Theory and their Applications in Engineering. Kluwer, Boston.
- 103. Sherwood, J.D. and Stone, H.A. (2001) Leakage through Filtercake into a Fluid Sampling Probe. Physics of Fluids, 13, 1151-1159. https://doi.org/10.1063/1.1360712
- 104. Warrick, A.W., Broadbridge, P. and Lomen, D.O. (1992) Approximations for Diffusion from a Disc Source. Applied Mathematical Modelling, 16, 155-161. https://doi.org/10.1016/0307-904X(92)90067-D
- 105. He, J.H. (2008) An Elementary Introduction to the Homotopy Perturbation Method. Computers & Mathematics with Applications, 57, 410-412. https://doi.org/10.1016/j.camwa.2008.06.003
- 106. Wazwaz, A.M. (2006) The Modified Decomposition Method for Analytic Treatment of Differential Equations. Applied Mathematics and Computation, 173, 165-176. https://doi.org/10.1016/j.amc.2005.02.048
- 107. Duan, J.S. and Rach, R. (2011) A New Modification of the Adomian Decomposition Method for Solving Boundary Value Problems for Higher Order Nonlinear Differential Equations. Applied Mathematics and Computation, 218, 4090-4118. https://doi.org/10.1016/j.amc.2011.09.037
- 108. Adomian, G. and Rach, R. (1983) Inversion of Nonlinear Stochastic Operators. Journal of Mathematical Analysis and Applications, 91, 39-46. https://doi.org/10.1016/0022-247X(83)90090-2
- 109. Adomian, G. (1986) Nonlinear Stochastic Operator Equations. Academic, Orlando.
- 110. Adomian, G. (1989) Nonlinear Stochastic Systems Theory and Applications to Physics. Kluwer Academic, Dordrecht. https://doi.org/10.1007/978-94-009-2569-4
- 111. Adomian, G., Rach, R. and Meyers, R. (1991) An Efficient Methodology for the Physical Sciences. Kybernetes, 20, 24-34. https://doi.org/10.1108/eb005909
- 112. Adomian, G. and Rach, R. (1993) Analytic Solution of Nonlinear Boundary Value Problems in Several Dimensions by Decomposition. Journal of Mathematical Analysis and Applications, 174, 118-137. https://doi.org/10.1006/jmaa.1993.1105
- 113. Adomian, G. and Rach, R. (1993) A New Algorithm for Matching Boundary Conditions in Decomposition Solutions. Applied Mathematics and Computation, 57, 61-68. https://doi.org/10.1016/0096-3003(93)90012-4
- 114. Adomian, G. and Rach, R. (1994) Modified Decomposition Solution of Linear and Nonlinear Boundary-Value Problems. Nonlinear Analysis, 23, 615-619. https://doi.org/10.1016/0362-546X(94)90240-2
- 115. Adomian, G. (1994) Solving Frontier Problems of Physics: The Decomposition Method. Kluwer Academic, Dordrecht. https://doi.org/10.1007/978-94-015-8289-6
- 116. Wazwaz, A.M. (2000) Approximate Solutions to Boundary Value Problems of Higher Order by the Modified Decomposition Method. Computers & Mathematics with Applications, 40, 679-691. https://doi.org/10.1016/S0898-1221(00)00187-5
- 117. Wazwaz, A.M. (2000) The Modified Adomian Decomposition Method for Solving Linear and Nonlinear Boundary Value Problems of 10th-Order and 12th-Order. International Journal of Nonlinear Sciences and Numerical Simulation, 1, 17-24. https://doi.org/10.1515/IJNSNS.2000.1.1.17
- 118. Wazwaz, A.M. (2001) A Reliable Algorithm for Obtaining Positive Solutions for Nonlinear Boundary Value Problems. Computers & Mathematics with Applications, 41, 1237-1244. https://doi.org/10.1016/S0898-1221(01)00094-3
- 119. Wazwaz, A.M. (2001) The Numerical Solution of Fifth-Order Boundary Value Problems by the Decomposition Method. Journal of Computational and Applied Mathematics, 136, 259-270. https://doi.org/10.1016/S0377-0427(00)00618-X
- 120. Wazwaz, A.M. (2001) The Numerical Solution of Sixth-Order Boundary Value Problems by the Modified Decomposition Method. Applied Mathematics and Computation, 118, 311-325. https://doi.org/10.1016/S0096-3003(99)00224-6
- 121. Wazwaz, A.M. (2001) A Reliable Algorithm for Solving Boundary Value Problems for Higher-Order Integro-Differential Equations. Applied Mathematics and Computation, 118, 327-342. https://doi.org/10.1016/S0096-3003(99)00225-8
- 122. Wazwaz, A.M. (2002) The Numerical Solution of Special Fourth-Order Boundary Value Problems by the Modified Decomposition Method. International Journal of Computer Mathematics, 79, 345-356. https://doi.org/10.1080/00207160211928
- 123. Jang, B. (2008) Two-Point Boundary Value Problems by the Extended Adomian Decomposition Method. Journal of Computational and Applied Mathematics, 219, 253-262. https://doi.org/10.1016/j.cam.2007.07.036
- 124. Ebaid, A.E. (2010) Exact Solutions for a Class of Nonlinear Singular Two-Point Boundary Value Problems: The Decomposition Method. Zeitschrift für Naturforschung, 65a, 1-6. https://doi.org/10.1515/zna-2010-0301
- 125. Adomian, G. (1986) A New Approach to the Heat Equation—An Application of the Decomposition Method. Journal of Mathematical Analysis and Applications, 113, 202-209. https://doi.org/10.1016/0022-247X(86)90344-6
- 126. Lesnic, D. and Elliott, L. (1999) The Decomposition Approach to Inverse Heat Conduction. Journal of Mathematical Analysis and Applications, 232, 82-98. https://doi.org/10.1006/jmaa.1998.6243
- 127. Wazwaz, A.M. (1999) A Reliable Modification of Adomian’s Decomposition Method. Applied Mathematics and Computation, 102, 77-86. https://doi.org/10.1016/S0096-3003(98)10024-3
- 128. Lesnic, D. (2001) A Computational Algebraic Investigation of the Decomposition Method for Time-Dependent Problems. Applied Mathematics and Computation, 119, 197-206. https://doi.org/10.1016/S0096-3003(99)00257-X