# Solutions of Two-Dimensional Nonlinear Sine-Gordon Equation via Triple Laplace Transform Coupled with Iterative Method

Authors: Alemayehu Tamirie Deresse; Yesuf Obsie Mussa; Ademe Kebede Gizaw
Journal: Journal of Applied Mathematics
 (2021)

---

## Abstract

This article presents triple Laplace transform coupled with iterative method to obtain the exact solution of two-dimensional nonlinear sine-Gordon equation (NLSGE) subject to the appropriate initial and boundary conditions. The noise term in this equation is vanished by successive iterative method. The proposed technique has the advantage of producing exact solution, and it is easily applied to the given problems analytically. Four test problems from mathematical physics are taken to show the accuracy, convergence, and the efficiency of the proposed method. Furthermore, the results indicate that the introduced method is promising for solving other type systems of NLPDEs.


---

## Body

## 1. Introduction

The sine-Gordon (SG) equation is a nonlinear hyperbolic PDE, which was originally considered in the nineteenth century in the course of study of surfaces of constant negative curvature often used to describe and simulate the physical phenomena in a variety of fields of engineering and science, such as nonlinear waves, propagation of fluxons, and dislocation of metals [1–4]. Because sine-Gordon equation leads to different types of soliton solutions, it has been receiving an enormous amount of attention. Soliton solution travels without experiencing any deformation through the medium, even when it collides with another soliton. The solitons, identified in many wave and particle systems, are of importance in the theory of nonlinear differential equations. As one of the crucial equations in nonlinear science, the sine-Gordon equation has been constantly investigated and solved numerically and analytically in recent years [5–8]. For the one-dimensional sine-Gordon equation, Maitama and Hamza [5] introduced analytical method called the Natural Decomposition Method (NDM) for solving nonlinear Sine-Gordon equation. Fayadh and Faraj [9] also applied combined Laplace transform method and VIM to get the approximate solution of the one-dimensional sine-Gordon equation.For the two-dimensional sine-Gordon equation, Su [6] obtained numerical solution of two-dimensional nonlinear sine-Gordon equation using localized method of approximate particular solutions. In [10], the author developed and analyzed an energy-conserving local discontinuous Galerkin method for the two-dimensional SGE on Cartesian grids. Duan et al. [11] proposed a numerical model based on lattice Boltzmann method to obtain the numerical solutions of two-dimensional generalized sine-Gordon equation, and the method was extended to solve the other two-dimensional wave equations, such as nonlinear hyperbolic telegraph equation as indicated in [12].In 2020, [13] developed a local Kriging meshless solution to the nonlinear (2+1)-dimensional sine-Gordon equation. The meshless shape function is constructed by Kriging interpolation method to have Kronecker delta function property for the two-dimensional field function, which leads to convenient implementation of imposing essential boundary conditions.In paper [14], the authors constructed kink wave solutions and traveling wave solutions of the (2+1)-dimensional sine-Gordon equation from the well-known AKNS system. For more related research about solving the sine-Gordon equations, readers may refer to [15–20].The Laplace transform method (LTM) is one of the integral transform methods that have been intensively used to solve linear and nonlinear equations [21]. The Laplace transform method is used frequently in engineering and physics; the output of a linear time invariant system can be calculated by convolving its unit impulse response with the input signal. Performing this calculation in Laplace space turns the convolution into a multiplication; the latter is easier to solve because of its algebraic form. The Laplace transform can also be used to solve differential equations and is used extensively in electrical engineering [22–25]. The Laplace transform reduces a linear differential equation to an algebraic equation, which can then be solved by the formal rules of algebra. The original differential equation can then be solved by applying the inverse Laplace transform.Recently, the concept of single Laplace transform is extended to double Laplace transform to solve some kind of differential equations and fractional differential equations such as linear/nonlinear space-time fractional telegraph equations, functional, integral, and partial differential equations [26–28]. Dhunde and Waghmare [29] applied the double Laplace transform method for solving a one-dimensional boundary value problems. Through this method, the boundary value problem is solved without converting it into ordinary differential equation; therefore, there is no need to find complete solution of ordinary differential equation. This is the biggest advantage of the proposed method.Furthermore, different scholars extended the double Laplace transform method to triple Laplace transform (TLT) to solve two-dimensional nonlinear partial differential equations arising in various natural phenomena. In [30], the authors used triple Laplace transform method to the solution of fractional-order partial differential equations by using Caputo fractional derivative. Through [31–34], the triple Laplace transform method was applied to obtain the solution of fractional-order telegraph equation in two dimensions, linear Volterra integro-differential equations in three dimensions, third-order Mboctara equations, and the proof of some of its properties like linearity property, change of scale property, first shifting property, and convolution theorem property, and differential property and triple integral property are given.Similar to the Laplace transform method, an iteration method (IM) is a fascinating task in an applied scientific branches to find the solution of nonlinear differential equation. The iterative procedure of the proposed method leads to a series, which can be summed up to find an analytical formula, or it can form a suitable approximation [35]. The error of the approximation can be controlled by properly truncating the series [36]. More surprisingly, an IM has showed effective and more rapid convergent series solution (see [37]).The purpose of this paper, is to apply triple Laplace transform (TLT) and iterative method (IM) developed in [38] to find the exact solution of two-dimensional nonlinear sine-Gordon equation (NLSGE) subject to appropriate initial and boundary conditions. Dhunde and Waghmare in [37, 39] applied double Laplace transform iteration method (TLTIM) to solve nonlinear Klein-Gordon and telegraph equations. By this method, the noise terms disappear in the iteration process, and single iteration gives the exact solution.In the present study, we are interested in the following two-dimensional sine-Gordon equation [40, 41].
(1)uttx,y,t+βutx,y,t=αuxxx,y,t+uyyx,y,t−ϕx,ysinux,y,t+hx,y,t,subject to the initial conditions
(2)ux,y,0=ϕ1x,y,utx,y,0=ϕ2x,y,x,y∈Ω,and boundary conditions (Cauchy type BCs)
(3)u0,y,t=g1y,t,ux0,y,t=g2y,t,ux,0,t=g3x,t,uyx,0,t=g4x,t,within the problem domain of Ω=x,y:a≤x≤b,c≤x≤d for t>0. In Equation (1), the parameter β≥0 is the damping factor for the dissipative term. When β=0, this equation will be reduced to an undamped sine-Gordon equation and when β>0 to the damped one. The function ϕx,y represents Josephson current density, and the functions ϕ1x,y and ϕ2x,y are specified wave modes and velocity.The remaining parts of this paper is structured as follows. In Section2, we begin with some basic definitions, properties, and theorems of triple Laplace transform method. Section 3 illustrates the details of the new iterative method and its convergence. Section 4 presents the description of the model; how the approximate analytical solutions of the given SGE equations is obtained using triple Laplace transform method coupled with iterative method. In Section 5, we apply the proposed method to four illustrative examples in order to show its liability, convergence, and efficiency. Finally, concluding remarks are given in Section 6.


## 2. Definitions and Properties of Triple Laplace Transform Method

In this section, we give some essential definitions, properties, and theorems of triple Laplace transform of partial differential equation, which should be used in the present study.Definition 1 (see [31]).
Letfx,y,t be a function of three variables x,y, and t defined in the positive xyt-plane. The triple Laplace transform of the function fx,y,t is defined by
(4)Lxytfx,y,t=Fk,p,s=∫0∞e−kx∫0∞e−py∫0∞e−stfx,y,tdtdydx,whenever the integral exist. Here, Lxytfx,y,t denotesLxLyLtfx,y,t, and k,p, and s are complex numbers. From this definition, we deduce that:
(5)Lxytfxgyht=FkGpHs=LxfxLygyLtht.Definition 2 (see [34]).
The inverse triple Laplace transform ofFk,p,s,Lxyt−1Fk,p,s=fx,y,t is given by the complex triple integral formula
(6)Lxyt−1Fk,p,s=fx,y,t=12πi∫q−i∞q+i∞e−kx12πi∫r−i∞r+i∞e−py12πi∫z−i∞z+i∞e−stFk,p,sdsdpdk,where Lxyt−1Fk,p,s denotes Lx−1Ly−1Lt−1Fk,p,s, and Fk,p,s must be an analytic function for all k,p, and s in the region defined by the inequality Rek≥q,Rep≥r, and Res≥z.Property 3 (see [34]).
The triple Laplace transform of second-order partial derivatives are given by(7)Lxyt∂2∂x2fx,y,t=k2Fk,p,s−kF0,p,s−Fx0,p,s,Lxyt∂2∂y2fx,y,t=p2Fk,p,s−pFk,0,s−Fyk,0,s,Lxyt∂2∂t2fx,y,t=s2Fk,p,s−sFk,p,0−Ftk,p,0,Lxyt∂2∂x∂yfx,y,t=kpFk,p,s−F0,p,s−Fk,0,s.Furthermore, the triple Laplace transform of first-order partial derivatives are given by(8)Lxyt∂∂xfx,y,t=sFk,p,s−F0,p,s,Lxyt∂∂yfx,y,t=sFk,p,s−Fk,0,s,Lxyt∂∂tfx,y,t=sFk,p,s−Fk,p,0.


### 2.1. Existence and Uniqueness of the Triple Laplace Transform

Theorem 4 (Existence).
Letfx,y,t be a continuous function on the interval 0,∞ which is of exponential order, that is, for some a,b,cϵℝ. Consider
(9)supx,y,t>0fx,y,texpax+by+ct<0.
Under this condition, the triple transform,Fk,p,s=∫0∞∫0∞∫0∞e−kxe−pye−stfx,y,tdxdydt, exists for all p>a,s>bandk>c and is in actuality infinitely differentiable with respect to p>a,s>bandk>c.Theorem 5 (Uniqueness).
Letfx,y,t and gx,y,t be continuous functions defined for x,y,t≥0 and having Laplace transforms, Fk,p,s and Gk,p,s, respectively. If Fk,p,s=Gk,p,s, then fx,y,t=gx,y,t.For the proof, see [34].


### 2.2. Some Properties of Triple Laplace Transform

Property 6 (Linearity Property of TLT).
Iffx,y,t and gx,y,t are two functions of x,y, and t such that Lxytfx,y,t=Fk,p,s and Lxytgx,y,t=Gk,p,s, then
(10)Lxytαfx,y,t+βgx,y,t=αLxytfx,y,t+βLxytgx,y,t=αFk,p,s+βGk,p,s,where α and β are constants.For the proof, see [31–33].Property 7 (Change of Scale Property).
LetLxytfx,y,t=Fk,p,s,then
(11)Lxytfax,by,ct=1abcFka,pb,sc,where a,b, and c are nonzero constants.For the proof, see [31, 33, 34].Property 8(First Shifting Property).
IfLxytfx,y,t=Fk,p,s, then
(12)Lxyteax+by+ctfx,y,t=Fk−a,p−b,s−c.For the proof, see [31, 33, 34].Property 9 (Second Shifting Property).
IfLxytfx,y,t=Fk,p,s, then
(13)Lxytfx−a,y−b,t−cHx−a,y−b,t−c=e−ak−bp−csFk,p,s,where Hx,y,t is the Heaviside unit step function defined by
(14)Hx,y,t=Hx−a,y−b,t−c:x>a,y>b,t>c,Hx−a,y−b,t−c:x<a,y<b,t<c.For the proof, see [31, 33].Property 10.
IfLxytfx,y,t=Fk,p,s, then Lxyteax+by+ct=1/k−ap−bs−c, for all x,y, and t.For the proof, see [31].


## 2.1. Existence and Uniqueness of the Triple Laplace Transform

Theorem 4 (Existence).
Letfx,y,t be a continuous function on the interval 0,∞ which is of exponential order, that is, for some a,b,cϵℝ. Consider
(9)supx,y,t>0fx,y,texpax+by+ct<0.
Under this condition, the triple transform,Fk,p,s=∫0∞∫0∞∫0∞e−kxe−pye−stfx,y,tdxdydt, exists for all p>a,s>bandk>c and is in actuality infinitely differentiable with respect to p>a,s>bandk>c.Theorem 5 (Uniqueness).
Letfx,y,t and gx,y,t be continuous functions defined for x,y,t≥0 and having Laplace transforms, Fk,p,s and Gk,p,s, respectively. If Fk,p,s=Gk,p,s, then fx,y,t=gx,y,t.For the proof, see [34].


## 2.2. Some Properties of Triple Laplace Transform

Property 6 (Linearity Property of TLT).
Iffx,y,t and gx,y,t are two functions of x,y, and t such that Lxytfx,y,t=Fk,p,s and Lxytgx,y,t=Gk,p,s, then
(10)Lxytαfx,y,t+βgx,y,t=αLxytfx,y,t+βLxytgx,y,t=αFk,p,s+βGk,p,s,where α and β are constants.For the proof, see [31–33].Property 7 (Change of Scale Property).
LetLxytfx,y,t=Fk,p,s,then
(11)Lxytfax,by,ct=1abcFka,pb,sc,where a,b, and c are nonzero constants.For the proof, see [31, 33, 34].Property 8(First Shifting Property).
IfLxytfx,y,t=Fk,p,s, then
(12)Lxyteax+by+ctfx,y,t=Fk−a,p−b,s−c.For the proof, see [31, 33, 34].Property 9 (Second Shifting Property).
IfLxytfx,y,t=Fk,p,s, then
(13)Lxytfx−a,y−b,t−cHx−a,y−b,t−c=e−ak−bp−csFk,p,s,where Hx,y,t is the Heaviside unit step function defined by
(14)Hx,y,t=Hx−a,y−b,t−c:x>a,y>b,t>c,Hx−a,y−b,t−c:x<a,y<b,t<c.For the proof, see [31, 33].Property 10.
IfLxytfx,y,t=Fk,p,s, then Lxyteax+by+ct=1/k−ap−bs−c, for all x,y, and t.For the proof, see [31].


## 3. The New Iterative Method

Consider the following functional equation [38]:
(15)ux¯=Nux¯+fx¯,where N is a nonlinear operator in a Banach space such that N:B⟶B, and f is a known function.x¯=x1,x2,x3,⋯,xn, and u is assumed to be the solution of Equation (15) having the series form:
(16)ux¯=∑i=0∞uix¯.The nonlinear operationN can then be decomposed as
(17)N∑i=0∞ui=Nu0+∑i=1∞N∑r=0iur−N∑r=0i−1ur.Using Equations (16) and (17), Equation (15) is equivalent to
(18)∑i=0∞ui=f+Nu0+∑i=1∞N∑r=0iur−N∑r=0i−1ur.From Equation (18), we define the following recurrence relation:
(19)u0=f,u1=Nu0,u2=Nu0+u1−Nu0,un+1=Nu0+⋯+un−Nu0+⋯+un−1,n=1,2,⋯.Thus,(20)u1+⋯+un+1=Nu0+⋯+un,n=1,2,⋯,and hence,
(21)∑i=0∞ui=f+N∑i=0∞ui.Therefore, thenth term approximate solution of Equation (16) is given by
(22)u=u0+∑i=1n−1ui=u0+u1+u2+⋯+un−1,n>1.Theorem 11.
IfN is a continuously differentiable functional in a neighborhood of u0 and Nnu0≤M≤e−1 for all n, then the series ∑n=o∞un+1 is absolutely convergent [35].


## 4. Description of the Method

Steps to be followed to apply triple Laplace Transform properties as in Table1 coupled with new iterative method are as follows:Step 1.
Applying triple Laplace transform on both sides of Equation (1) and using Property 3, we get
(23)s2u¯k,p,s−su¯k,p,0−u¯tk,p,0+βsu¯k,p,s−u¯k,p,0+α−p2u¯k,p,s+pu¯k,0,s+u¯yk,0,s−k2u¯k,p,s+ku¯0,p,s+u¯x0,p,s+Lxytϕx,ysinux,y,t=Lxythx,y,t,where hx,y,t is the source functioning Equation (1).Step 2.
Now, employing double Laplace transform to Equations (2) and (3), we have
(24)φ1¯k,p=u¯k,p,0,φ2¯k,p=u¯tk,p,0,andg1¯p,s=u¯0,p,s,g2¯p,s=u¯x0,p,s,(25)g3¯k,s=u¯k,0,s,g4¯k,s=u¯yk,0,s,respectively.Table 1
Triple Laplace transform for some functions of three variables [33, 34].
Functionsfx,y,tTriple Laplace transformFk,p,sabcabckpsxyt1k2p2s2xmyntrm!n!r!km+1pn+1sr+1e−ax−by−ct1a+kb+pc+scosxcosycostkps1+k21+p21+s2sinxsinysint11+k21+p21+s2sinx+y+t−1+ps+kp+s1+k21+p21+s2cosx+y+t−k+p+s−kps1+k21+p21+s2coshax+by+ct121k−ap−bs−c+1k+ap+bs+csinhax+by+ct121k−ap−bs−c−1k+ap+bs+cBy substituting Equations (24) and (25) into Equation (23) and simplifying, we obtain
(26)u¯k,p,s=1s2−αp2−αk2+βss+βu¯k,p,0+u¯tk,p,0−αpu¯k,0,s+u¯yk,0,s+ku¯0,p,s+u¯x0,p,s+Lxythx,y,t−ϕx,ysinux,y,t.Step 3.
By implementing the triple inverse Laplace transformation of Equation (26), we obtain
(27)u¯k,p,s=Lxyt−11s2−αp2−αk2+βss+βu¯k,p,0+u¯tk,p,0−αpu¯k,0,s+u¯yk,0,s+ku¯0,p,s+u¯x0,p,s+Lxythx,y,t−ϕx,ysinux,y,t.Step 4.
Assume that(28)ux,y,t=∑i=0∞uix,y,t,is the solution of Equation (1).Then, substituting Equation (28) into Equation (27), we obtain
(29)∑i=0∞uix,y,t=Lxyt−11s2−αp2−αk2+βss+βu¯k,p,0+u¯tk,p,0−αpu¯k,0,s+u¯yk,0,s+ku¯0,p,s+u¯x0,p,s+Lxythx,y,t−ϕx,ysin∑i=0∞uix,y,t.Step 5.
By implementing the new iterative method, the nonlinear termsinux,y,t in Equation (29) is decomposed as
(30)sin∑i=0∞uix,y,t=sinu0x,y,t+∑i=1∞sin∑r=0iurx,y,t−sin∑r=0i−1urx,y,t.Using Equation (30), Equation (29) is equivalent to
(31)∑i=0∞uix,y,t=Lxyt−11s2−αp2−αk2+βss+βu¯k,p,0+u¯tk,p,0−αpu¯k,0,s+u¯yk,0,s+ku¯0,p,s+u¯x0,p,s+Lxythx,y,t−ϕx,ysinu0x,y,t+∑i=1∞sin∑r=0iurx,y,t−sin∑r=0i−1urx,y,t.Step 6.
Using triple Laplace transform coupled with new iterative method, we introduce the recursive relations and get(32)u0x,y,t=Lxyt−11s2−αp2−αk2+βss+βu¯k,p,0+u¯tk,p,0−αpu¯k,0,s+u¯yk,0,s+ku¯0,p,s+u¯x0,p,s,u1x,y,t=Lxyt−11s2−αp2−αk2+βsLxythx,y,t−ϕx,ysinu0x,y,t,un+1x,y,t=−Lxyt−11s2−αp2−αk2+βsLxytϕx,y∑i=1∞sin∑r=0iurx,y,t−sin∑r=0i−1urx,y,t,n≥1.Step 7.
The solution of Equations (1)–(3) in series form is given by
(33)ux,y,t=u0x,y,t+u1x,y,t+u2x,y,t+⋯=u0x,y,t.


## 5. Illustrative Examples

In order to show the validity and effectiveness of the method under consideration, some examples are presented here.Example 1.
Consider the two-dimensional NLSGE1 with β=0,α=1,ϕx,y=1, hx,y,t=sinex+y−2t+2ex+y−2t on the domain Ω=0,22,t≥0. That is,
(34)utt=uxx+uyy−sinux,y,t+sinex+y−2t+2ex+y−2t,with initial conditions
(35)ux,y,0=ex+y,utx,y,0=−2ex+y,and boundary conditions
(36)u0,y,t=ey−2t,ux0,y,t=ey−2t,ux,0,t=ex−2t,uyx,y,0=ex−2t.Solution: Applying properties of triple Laplace transform on both sides of Equation (34), we get
(37)s2u¯k,p,s−su¯k,p,0−u¯tk,p,0−p2u¯k,p,s+pu¯k,0,s+u¯yk,0,s−k2u¯k,p,s+ku¯0,p,s+u¯x0,p,s+Lxytsinux,y,t=Lxytsinex+y−2t+2k−1p−1s+2.Applying double Laplace transform to Equations (35) and (36), we obtain
(38)u¯k,p,0=1k−1p−1,u¯tk,p,0=−2k−1p−1,u¯0,p,s=1p−1s+2,u¯x0,p,s=1p−1s+2,(39)u¯k,0,s=1k−1s+2,u¯yk,0,s=1k−1s+2,respectively.By substituting Equations (38) and (39) into Equation (37), we get
(40)s2−p2−k2u¯k,p,s=2k−1p−1s+2+s−2k−1p−1−p+1k−1s+2−k+1p−1s+2+Lxytsinex+y−2t−sinux,y,t.Simplifying (40) gives us
(41)u¯k,p,s=1k−1p−1s+2+1s2−p2−k2Lxytsinex+y−2t−sinux,y,t.Applying inverse triple Laplace transform to Equation (41), we get
(42)ux,y,t=ex+y−2t+Lxyt−11s2−p2−k2Lxytsinex+y−2t−sinux,y,t.Now, applying the new iterative method to Equation (42), we obtain the components of the solution as follows:
(43)u0x,y,t=ex+y−2t,(44)u1x,y,t=Lxyt−11s2−p2−k2Lxytsinex+y−2t−sinex+y−2t=0.(45)un+1ux,y,t=−Lxyt−11s2−p2−k2Lxytsin∑r=0nurx,y,t−sin∑r=0n−1urx,y,t,n≥1.Now, we define the recurrence relation from Equation (45) for n≥1 as follows:
(46)u2x,y,t=Lxyt−11s2−p2−k2Lxytsinu0x,y,t+u1x,y,t−sinu0x,y,t=0,u3x,y,t=Lxyt−11s2−p2−k2Lxytsinu0x,y,t+u1x,y,t+u2x,y,t−sinu0x,y,t+u1x,y,t=0.In the same way, we obtain u4x,y,t=u5x,y,t=0 and so on.Therefore, the solution of Example1 by using Equation (33) is
(47)ux,y,t=ex+y−2t.Let us now test the convergence of the obtained series solution. From Equation (42), we have
(48)u0x,y,t=ex+y−2t,Nux,y,t=Lxyt−11s2−p2−k2Lxytsinex+y−2t−sinux,y,t.Thus, for allx,y,t≥0, we have
(49)Nu0x,y,t=Lytx−11s2−p2−k2Lxytsinex+y−2t−sinu0x,y,t=Lxyt−11s2−p2−k2Lxytsinex+y−2t−sinex+y−2t=0,Therefore,Nu0x,y,t=0=0<1/e.
(50)N′ux,y,t=Lytx−11s2−p2−k2LxLyLtex+y−2t′cosex+y−2t−cosux,y,tu′x,y,t,where N′ux,y,t represents the partial derivatives ∂/∂xux,y,t or ∂/∂yux,y,t or ∂/∂tux,y,t.(51)N∂∂tux,y,t=Lxyt−11s2−p2−k2Lxyt−2ex+y−2tcosex+y−2t−cosux,y,t∂∂tux,y,t.Then,(52)N∂∂tu0x,y,t=Lxyt−11s2−p2−k2Lxyt−2ex+y−2tcosex+y−2t+2ex+y−2tcosex+y−2t=0,N∂∂xu0x,y,t=N∂∂yu0x,y,t=Lxyt−11s2−p2−k2Lxytex+y−2tcosex+y−2t−ex+y−2tcosex+y−2t=0.Therefore,Nu0′x,y,t=0=0<1/e.
(53)N∂2∂t2u0x,y,t=Lxyt−11s2−p2−k2Lxyt4ex+y−2tcosex+y−2t−4e2x+2y−4tsinex+y−2t−cosu0x,y,t∂2∂t2u0x,y,t+sinux,y,t∂∂tu0x,y,t2=Lxyt−11s2−p2−k2Lxyt4ex+y−2tcosex+y−2t−4e2x+2y−4tsinex+y−2t−4ex+y−2tcosex+y−2t+4e2x+2y−4tsinex+y−2t=0,N∂2∂x2u0x,y,t=N∂2∂y2u0x,y,t=N∂2∂x∂yu0x,y,t=Lxyt−11s2−p2−k2Lxytex+y−2tcosex+y−2t−e2x+2y−4tsinex+y−2t−ex+y−2tcosex+y−2t+e2x+2y−4tsinex+y−2t=0,N∂2∂x∂tu0x,y,t=N∂2∂y∂tu0x,y,t=Lxyt−11s2−p2−k2Lxyt−2ex+y−2tcosex+y−2t+2e2x+2y−4tsinex+y−2t+2ex+y−2tcosex+y−2t−2e2x+2y−4tsinex+y−2t=0.Therefore,N′′u0x,y,t=0=0<1/e,whereN′′u0x,y,t represents all the second-order partial derivatives ofu0x,y,t.Similarly, by principle of Mathematical induction, we haveN3u0x,y,t=N4u0x,y,t=⋯=Nku0x,y,t=0<1/e, for all k≥0.As the condition of Theorem11 are satisfied, the series solution ux,y,t=∑i=0∞uix,y,t is absolutely convergent, and hence, the solution obtained by the new iterative method is convergent on the domain of interest.Example 2.
Consider Equation (1) with β=0,α=1,ϕx,y=1, hx,y,t=sinsinx,y,t+sinx,y,t on the domain Ω=0,2π2,t≥0. That is,
(54)uttx,y,t=uxxx,y,t+uyyx,y,t−sinux,y,t+sinsinx+y+t+sinx+y+t,with initial conditions
(55)ux,y,0=sinx+y,utx,y,0=cosx+y,and boundary conditions
(56)u0,y,t=siny+t,ux0,y,t=cosy+t,ux,0,t=sinx+t,uyx,0,t=cosx+t.Solution: Applying properties of triple Laplace transform on both sides of Equation (54), we obtain
(57)s2u¯k,p,s−su¯k,p,0−u¯tk,p,0−p2u¯k,p,s+pu¯k,0,s+u¯yk,0,s−k2u¯k,p,s+ku¯0,p,s+u¯x0,p,s+Lxytsinux,y,t=Lxytsinsinx+y+t+−1+ps+kp+s1+k21+p21+s2.Applying double Laplace transform to Equations (55) and (56), we obtain
(58)u¯k,p,0=k+p1+p21+k2,u¯tk,p,0=kp−11+p21+k2,(59)u¯0,p,s=k+s1+p21+s2,u¯x0,p,s=ps−11+p21+s2,u¯k,0,s=k+s1+k21+s2,u¯yk,0,s=ks−11+k21+s2,respectively.By substituting Equations (58) and (59) in to Equation (57) and simplifying, we get
(60)u¯k,p,s=−1+ps+kp+s1+k21+p21+s2+1s2−p2−k2Lxytsinsinx+y+t−sinux,y,t.Applying inverse triple Laplace transform to Equation (60), we get
(61)ux,y,t=sinx+y+t+Lxyt−11s2−p2−k2Lxytsinsinx+y+t−sinux,y,t.Now, applying the new iterative method to Equation (61), we obtain the components of the solution as follows:Substituting Equation (16) into Equation (61), we obtain the components of the solution as follows:
(62)u0x,y,t=sinx+y+t,u1x,y,t=Lxyt−11s2−p2−k2Lxytsinsinx+y+t−sinsinx+y+t=0,u2x,y,t=Lxyt−11s2−p2−k2Lxytsinu0x,y,t+u1x,y,t−sinu0x,y,t=0,u3x,y,t=Lxyt−11s2−p2−k2Lxytsinu0x,y,t+u1x,y,t+u2x,y,t−sinu0x,y,t+u1x,y,t=0,and so on.Therefore, the solution of Example2 by using equation (33) is
(63)ux,y,t=sinx+y+t,as indicated in Baccouch [10] and Deresse et al. [41].Next, we test the convergence of the obtained series solution. From Equation (61) we have,
(64)u0x,y,t=sinx+y+tandNux,y,t=Lxyt−11s2−p2−k2Lxytsinsinx+y+t−sinux,y,t.Thus, for allx,y,t≥0, we have
(65)Nu0x,y,t=Lxyt−11s2−p2−k2Lxytsinsinx+y+t−sinu0x,y,t=Lxyt−11s2−p2−k2Lxytsinsinx+y+t−sinsinx+y+t=0.Therefore,Nu0x,y,t=0=0<1/e.
(66)N∂∂tux,y,t=Lxyt−11s2−p2−k2Lxytcossinx+y+tcosx+y+t−cosux,y,t∂∂tux,y,t,N∂∂tu0x,y,t=Lxyt−11s2−p2−k2Lxytcossinx+y+tcosx+y+t−cossinx+y+tcosx+y+t=0.Similarly,N∂/∂xu0x,y,t=N∂/∂yu0x,y,t=0.Therefore,Nu0′x,y,t=0=0<1/e.
(67)N∂2∂t2u0x,y,t=Lxyt−11s2−p2−k2Lxyt−sinsinx+y+tcos2x+y+t+cossinx+y+tsinx+y+t−cosu0x,y,t∂2∂t2u0x,y,t−sinux,y,t∂∂tu0x,y,t2=Lxyt−11s2−p2−k2Lxyt−sinsinx+y+tcos2x+y+t−cossinx+y+tsinx+y+tsinsinx+y+tcos2x+y+t+cossinx+y+tsinx+y+t=0.Similarly,(68)N∂2∂x2u0x,y,t=N∂2∂y2u0x,y,t=N∂2∂x∂yu0x,y,t=N∂2∂x∂tu0x,y,t=N∂2∂y∂tu0x,y,t=0.Therefore,N′′u0x,y,t=0=0<1/e,whereN′′u0x,y,t represents all the second-order partial derivatives of u0x,y,t.Similarly, by principle of Mathematical induction, we haveN3u0x,y,t=N4u0x,y,t=⋯=Nku0x,y,t=0<1/e, for all k≥0.As the condition of Theorem11 are satisfied, the series solution ux,y,t=∑i=0∞uix,y,t is absolutely convergent, and hence, the solution obtained by the new iterative method is convergent on the domain of interest.Example 3.
Consider Equation (1) on the domain Ω=−1/2,1/22 with β=0,α=1/2π2,ϕx,y=1, and hx,y,t=sincosπxcosπycost+sinx,y,t, That is,
(69)uttx,y,t=12π2uxxx,y,t+uyyx,y,t−sinux,y,t+sincosπxcosπycost,with initial conditions
(70)ux,y,0=cosπxcosπy,utx,y,0=0,and boundary conditions
(71)u0,y,t=cosπycost,ux0,y,t=0,ux,0,t=cosπxcost,uyx,0,t=0.Solution: Applying properties of triple Laplace transform on both sides of Equation (69), we get
(72)s2u¯k,p,s−su¯k,p,0−u¯tk,p,0+12π2−p2u¯k,p,s+pu¯k,0,s+u¯yk,0,s−k2u¯k,p,s+ku¯0,p,s+u¯x0,p,s+LxLyLtsinux,y,t=LxLyLtsincosπxcosπycost.Applying double Laplace transform to Equations (70) and (71), we obtain
(73)u¯k,p,0=kpπ2+p2π2+k2,u¯tk,p,0=0,(74)u¯k,0,s=ksπ2+k21+s2,u¯yk,0,s=0,u¯0,p,s=psπ2+p21+s2,u¯x0,p,s=0,respectively.By substituting Equations (73) and (74) into Equation (72) and simplifying, we get
(75)u¯k,p,s=kpsπ2+p2π2+k21+s2+2π22π2s2−p2−k2Lxytsincosπxcosπycost−sinux,y,t.Applying inverse triple Laplace transform to Equation (75), we obtain
(76)ux,y,t=cosπxcosπycost+Lxyt−12π22π2s2−p2−k2Lxytsincosπxcosπycost−sinux,y,t.Now, applying the new iterative method to Equation (76), we obtain the components of the solution as follows:
(77)u0x,y,t=cosπxcosπycost,u1x,y,t=Lxyt−12π22π2s2−p2−k2Lxytsincosπxcosπycost−sinu0x,y,t,=Lxyt−12π22π2s2−p2−k2Lxytsincosπxcosπycost−sincosπxcosπycost=0,(78)u2x,y,t=Lxyt−12π22π2s2−p2−k2Lxytsinu0x,y,t+u1x,y,t−sinu0x,y,t=0,u3x,y,t=Lxyt−12π22π2s2−p2−k2Lxytsinu0x,y,t+u1x,y,t+u2x,y,t−sinu0x,y,t+u1x,y,t=0,and so on.Therefore, the solution of Example3 by using Equation (33) is
(79)ux,y,t=cosπxcosπycost.This is exactly the same as the result obtained by Kang et al. [16] and Deresse et al. [41].We now test the convergence of the obtained series solution of Example3.From Equation (61), we have u0x,y,t=cosπxcosπycost and
(80)Nux,y,t=Lxyt−12π22π2s2−p2−k2Lxytsincosπxcosπycost−sinux,y,t.Thus, for allx,y,t≥0, we have
(81)Nu0x,y,t=Lxyt−12π22π2s2−p2−k2Lxytsincosπxcosπycost−sinu0x,y,t=Lxyt−12π22π2s2−p2−k2Lxytsincosπxcosπycost−sincosπxcosπycost=0.Therefore,Nu0x,y,t=0=0<1/e.
(82)N∂∂tux,y,t=Lxyt−12π22π2s2−p2−k2Lxyt−coscosπxcosπycostcosπxcosπysint−cosux,y,t∂∂tux,y,t,N∂∂tu0x,y,t=Lxyt−12π22π2s2−p2−k2Lxyt−coscosπxcosπycostcosπxcosπysint+coscosπxcosπycostcosπxcosπysint=0.Similarly,N∂/∂xu0x,y,t=N∂/∂yu0x,y,t=0.Therefore,Nu0′x,y,t=0=0<1/e.(83)N∂2∂t2u0x,y,t=Lxyt−12π22π2s2−p2−k2Lxyt−sincosπxcosπycostcos2πxcos2πysin2t−coscosπxcosπycostcosπxcosπycost−cosu0x,y,t∂2∂t2u0x,y,t−sinux,y,t∂∂tu0x,y,t2=Lxyt−12π22π2s2−p2−k2Lxyt−sincosπxcosπycostcos2πxcos2πysin2t−coscosπxcosπycostcosπxcosπycost+sincosπxcosπycostcos2πxcos2πysin2t+coscosπxcosπycostcosπxcosπycost=0.Similarly,(84)N∂2∂x2u0x,y,t=N∂2∂y2u0x,y,t=N∂2∂x∂yu0x,y,t=N∂2∂x∂tu0x,y,t=N∂2∂y∂tu0x,y,t=0,Therefore,Nu0′′x,y,t=0=0<1/e,whereN′′u0x,y,t represents all the second order partial derivatives ofu0x,y,t.Similarly, by principle of Mathematical induction, we haveN3u0x,y,t=N4u0x,y,t=⋯=Nku0x,y,t=0<1/e, for all k≥0.As the condition of Theorem11 are satisfied, the solution series ux,y,t=∑i=0∞uix,y,t is absolutely convergent, and hence, the solution obtained by the new iterative method is convergent on the domain of interest.Example 4.
Consider Equation (1) with β=1,α=1,ϕx,y=2, hx,y,t=2sine−t1−cosπx1−cosπy−π2e−tcosπx+cosπy−2cosπxcosπy on the domain Ω=0,22,t≥0. That is,
(85)uttx,y,t+utx,y,t=uxxx,y,t+uyyx,y,t−2sinux,y,t+2sine−t1−cosπx1−cosπy−π2e−tcosπx+cosπy−2cosπxcosπy,with initial conditions
(86)ux,y,0=1−cosπx1−cosπy,utx,y,0=−1−cosπx1−cosπy,and boundary conditions
(87)u0,y,t=0,ux0,y,t=0,ux,0,t=0,uyx,0,t=0.Solution: Applying properties of triple Laplace transform on both sides of Equation (85), we obtain
(88)s2u¯k,p,s−su¯k,p,0−u¯tk,p,0+su¯k,p,s−u¯k,p,0−p2u¯k,p,s+pu¯k,0,s+u¯yk,0,s−k2u¯k,p,s+ku¯0,p,s+u¯x0,p,s+2Lxytsinux,y,t=2Lxytsine−t1−cosπx1−cosπy−π2s+1kk2+π2+pp2+π2−2kpk2+π2p2+π2.Applying double Laplace transform to Equations (86) and (87), we obtain
(89)u¯k,p,0=1s−kk2+π2−pp2+π2+kpk2+π2p2+π2,u¯tk,p,0=−1s−kk2+π2−pp2+π2+kpk2+π2p2+π2,and
(90)u¯k,0,s=0,u¯yk,0,s=0,u¯0,p,s=0,u¯x0,p,s=0,respectively.By substituting Equations (89) and (90) into Equation (88) and simplifying, we get
(91)u¯k,p,s=1s+11s−kk2+π2−pp2+π2+kpk2+π2p2+π2+2s2+s−p2−k2Lxytsine−t1−cosπx1−cosπy−sinux,y,t.Applying inverse triple Laplace transform to Equation (91), we obtain
(92)ux,y,t=e−t1−cosπx1−cosπy+Lxyt−12s2+s−p2−k2Lxytsine−t1−cosπx1−cosπy−sinux,y,t.Now, we apply the iteration process.Substituting Equation (16) into Equation (92), we obtain the components of the solution as follows:
(93)u0x,y,t=e−t1−cosπx1−cosπy,u1x,y,t=Lxyt−12s2+s−p2−k2Lxytsine−t1−cosπx1−cosπy−sinu0x,y,t,=Lxyt−12s2+s−p2−k2Lxytsine−t1−cosπx1−cosπy−sine−t1−cosπx1−cosπy=0,u2x,y,t=Lxyt−12s2+s−p2−k2Lxytsinu0x,y,t+u1x,y,t−sinu0x,y,t=0,u3x,y,t=Lxyt−12s2+s−p2−k2Lxytsinu0x,y,t+u1x,y,t+u2x,y,t−sinu0x,y,t+u1x,y,t=0,and so on.Therefore, the solution of Example4 by using Equation (33) is
(94)ux,y,t=e−t1−cosπx1−cosπy.This result shows an excellent agreement with the one obtained in [12, 40, 41].Now, we test the convergence of the obtained series solution. From Equation (61), we have
(95)u0x,y,t=e−t1−cosπx1−cosπy,Nux,y,t=Lxyt−12s2+s−p2−k2Lxytsine−t1−cosπx1−cosπy−sinux,y,t.Thus, for allx,y,t≥0, we have
(96)Nu0x,y,t=Lxyt−12s2+s−p2−k2Lxytsine−t1−cosπx1−cosπy−sinu0x,y,t=Lxyt−12s2+s−p2−k2Lxytsine−t1−cosπx1−cosπy−sine−t1−cosπx1−cosπy=0.Therefore,Nu0x,y,t=0=0<1/e.
(97)N∂∂tux,y,t=Lxyt−12s2+s−p2−k2Lxyt−cose−t1−cosπx1−cosπye−t1−cosπx1−cosπy−cosux,y,t∂∂tux,y,t,N∂∂tu0x,y,t=Lxyt−12s2+s−p2−k2Lxyt−cose−t1−cosπx1−cosπye−t1−cosπx1−cosπy+cose−t1−cosπx1−cosπye−t1−cosπx1−cosπy=0.Similarly,N∂/∂xu0x,y,t=N∂/∂yu0x,y,t=0.Therefore,Nu0′x,y,t=0=0<1/e.
(98)N∂2∂t2u0x,y,t=Lxyt−12s2+s−p2−k2Lxytsine−t1−cosπx1−cosπye−2t1−cosπx21−cosπy2−cose−t1−cosπx1−cosπye−t1−cosπx1−cosπy−cosu0x,y,t∂2∂t2u0x,y,t−sinux,y,t∂∂tu0x,y,t2=Lxyt−12s2+s−p2−k2Lxyte−t1−cosπx1−cosπye−2t1−cosπx21−cosπy2−cose−t1−cosπx1−cosπye−t1−cosπx1−cosπy−e−t1−cosπx1−cosπye−2t1−cosπx21−cosπy2+cose−t1−cosπx1−cosπye−t1−cosπx1−cosπy=0.Similarly,(99)N∂2∂x2u0x,y,t=N∂2∂y2u0x,y,t=N∂2∂x∂yu0x,y,t=N∂2∂x∂tu0x,y,t=N∂2∂y∂tu0x,y,t=0.Therefore,Nu0′′x,y,t=0=0<1/e,whereN′′u0x,y,t represents all the second order partial derivatives ofu0x,y,t.Similarly, by principle of Mathematical induction, we haveN3u0x,y,t=N4u0x,y,t=⋯=Nku0x,y,t=0<1/e, for all k≥0.As the condition of Theorem11 are satisfied, the series solution ux,y,t=∑i=0∞uix,y,t is absolutely convergent, and hence, the solution obtained by the new iterative method is convergent on the domain of interest.


## 6. Conclusion

In this paper, triple Laplace transform coupled with iterative method is applied to obtain exact solution of two-dimensional nonlinear Sine-Gordon equation subject to initial and boundary conditions. Four illustrative examples are presented to show the validity of the method under consideration. The solutions of Examples1, 2, 3, and 4 obtained by the proposed method are in an excellent agreement with the same problem that has been considered in [6, 10–12, 16, 40, 41], and further, nontrivial problems that are solved using earlier methods become trivial in the sense that the decomposition ux,y,t=∑n=0∞unx,y,t=u0x,y,t+u1x,y,t+u2x,y,t+⋯+umx,y,t+⋯ consists of only one term, i.e.,ux,y,t=u0x,y,t. From this study, we concluded that triple Laplace transform coupled with iterative method finds quite practical analytical results with less computational work.



---

*Source: 9279022-2021-10-04.xml*