Double Sumudu Transform Iterative Method for One- Dimensional Nonlinear Coupled Sine-Gordon Equation

In this paper, the combined double Sumudu transform with iterative method is successfully implemented to obtain the approximate analytical solution of the one-dimensional coupled nonlinear sine-Gordon equation (NLSGE) subject to the appropriate initial and boundary conditions which cannot be solved by applying double Sumudu transform only. The solution of the nonlinear part of this equation was solved by a successive iterative method, the proposed technique has the advantage of producing an exact solution, and it is easily applied to the given problems analytically. Two test problems from mathematical physics were taken to show the liability, accuracy, convergence, and efficiency of the proposed method. Furthermore, the results indicate that the introduced method is promising for solving other types of systems of NLPDEs.


Introduction
The system of coupled sine-Gordon equations in the form was introduced by Khusnutdinova and Pelinovsky [1]. The coupled sine-Gordon equations generalize the Frenkel-Kontorova dislocation model [2,3]. System (1) with α = 1 was also proposed to describe the open states in the DNA double helix model [4]. This type of equation has been receiving an enormous amount of attention due to the presence of soliton solutions, and it has applications in the propagation of fluxons in Josephson junctions between two superconductors [5], the motion of a rigid pendulum attached to a stretched wire [6], etc.
Recently, many researchers used various methods to get the solution of system (1). Li et al. [7] presented a lattice Boltzmann model for the two-component system of coupled sine-Gordon equations by using the coupled mesoscopic Boltzmann equations and the Chapman-Enskog multiscale expansion, while the macroscopical governing evolution system can be recovered correctly by selecting suitable discrete equilibrium distribution functions and the amending functions. In paper [8], the authors used a rational exponential ansatz to derive the exact solutions of a coupled sine-Gordon equation. The simplest equation method has been used for finding the exact solutions of coupled sine-Gordon equations by Zhao [9]. In papers [10][11][12], Zhao et al. obtained some new solutions of coupled sine-Gordon equations including Jacobi elliptic function solutions, hyperbolic function solutions, and trigonometric function solutions by the Jacobi elliptic function expansion method, the hyperbolic auxiliary function method, and the symbolic computation method. Sadighi et al. [13] employed the homotopy perturbation method (HPM) for solving both sine-Gordon and coupled sine-Gordon equations. In paper [14], the authors obtained numerical solution of inhomogeneous systems of sine-Gordon equation (1) by finite difference method with fixed-point iteration.
In 2018, the authors of [15] studied the coupled sine-Gordon equations in nonlinear optics, which describe the propagation of an optical pulse in fiber waveguide, and they derived the new exact solutions of the problem through the use of the well-organized modified Kudryashov method. For more related research about solving the sine-Gordon equations, interested readers may refer to [16][17][18].
Since many application problems in science and engineering are formulated by nonlinear partial differential equations, new techniques are required to solve those problems efficiently. The Sumudu transform is an integral transform similar to the Laplace transform, introduced in the early 1990s by Watugala [19] to solve differential equations and control engineering problems. Recently, this method is used to solve differential equations like nonlinear boundary value problems [20] and the Cauchy problem for the wave equation in one dimensional space [21]. In [22], the Sumudu transform is applied to arbitrary powers Dumont bimodular Jacobi elliptic functions for arbitrary powers. Furthermore, Watugala [23] has extended the Sumudu transform to two variables with emphasis on solutions to partial differential equations. Ahmeda et al. [24] studied the convergence properties of double Sumudu transformation and used it to find the exact solution of the Volterra integropartial differential equation. In the paper [25], the authors' applied the double Sumudu transform linked with the Adomian decomposition method or with variational iteration method to find the analytical solution of nonlinear fractional partial differential equations. Kiwne and Sonawane [26] present the proofs of some theorems of double Sumudu transform with application in solving partial differential equations like one-dimensional heat equation, Poisson equation, and second-order hyperbolic equations of two variables. The advantage of the double Sumudu transform is it gives a rapid convergence of the exact solution without any restrictive assumption of the solution compared to other known methods (see [27]). Unfortunately, this transform fails to solve nonlinear partial differential equations like other integral transforms; to solve this problem, this transform is often combined with other methods like the variational method [28], reduced differential transform method [29], double integral transform (Laplace-Sumudu transform) method [30], Adomian decomposition method [25,31], and homotopy perturbation method [20,32].
The main aim of this work is to apply the double Sumudu transform method coupled with the new iterative method (NIT) proposed by Daftardar-Gejji and Jafari in [33] to find an exact/approximate solution of the nonlinear coupled sine-Gordon equation. The new iterative method (NIM) has been extensively used by many researchers for the treatment of linear and nonlinear ordinary and partial differential equations of integer and fractional order (see [33][34][35][36][37]). The method converges to the exact solution if it exists through successive approximations. The nonlinear term in the equation is expanded in terms of Daftardar-Gejji and Jafari polynomials (see [38,39]).
In the present study, we consider the nonlinear initial boundary value problem (IBVP) for the one-dimensional coupled sine-Gordon equation [14] in the following form: subject to the initial conditions and boundary conditions (Cauchy-type BCs) within the problem domain of Ω = fx : a ≤ x ≤ bg and t > 0.
The remaining parts of this paper are structured as follows. In Section 2, we begin with some basic definitions, properties, and theorems of the double Sumudu transform method. Section 3 illustrates the details of the new iterative method and its convergence. Section 4 presents how to obtain an approximate analytical solution of the coupled SGE using the DSTIM. In Section 5, we apply the proposed method to two illustrative examples in order to show its liability, convergence, and efficiency. Finally, concluding remarks are given in Section 6.

A Brief Introduction of Double Sumudu Transform (DST)
Definition 1 (see [24,40]). Sumudu transform of a function f ðtÞ is defined for all real numbers t > 0 as the function Sð pÞ, given by Definition 2 (see [30]). The inverse Sumudu transform of a function FðpÞ is denoted by the symbol S −1 ½FðpÞ = f ðtÞ and is defined by the Bromwich contour integral Definition 3 (see [24,40,41]). Let f ðx, tÞ be a function that can be expressed as a convergent infinite series and let ðx, t Þ ∈ ℝ + n , then the double Sumudu transform of the function f ðx, tÞ in the positive quadrant of the xt-plane is given by where x, t ≥ 0, and k, p are the transform variables for x and t, respectively, whenever the improper integral is convergent.

Advances in Mathematical Physics
Definition 4 (see [42,43]). The inverse double Sumudu transform S x −1 S t −1 ½Fðk, pÞ = f ðx, tÞ is defined by the following form: Remark 5 (see [26,40]). If f ðx, tÞ is a continuous function having a second partial derivative, then double Sumudu transforms of partial derivative of the first and second are defined as follows: Double Sumudu transform for the first partial derivative with respect to t is given by Double Sumudu transform for the first partial derivative with respect to x is given by Double Sumudu transform for the second partial derivative with respect to t is given by Double Sumudu transform for the second partial derivative with respect to x is given by

Convergence Theorem of Double Sumudu Transform
Lemma 6 (see [24]). Let f ðx, tÞ be a function of two variables continuous in the positive quadrant of the xt-plane. If the integral converges at k = k 0 , p = p 0 , then the integral converges for k < k 0 , p < p 0 :

Existence and Uniqueness of the Double Sumudu
Transform. If f ðx, tÞ is an exponential order, then c and d as x ⟶ ∞,t⟶∞, and if there exists a positive constant K such that for all x > X, t > T, then And we write f ðx, tÞ = Oe cx+dt as ⟶∞, t ⟶ ∞: Or, equivalently, And we write this as f ðx, tÞ = Oðe ax+by+ct Þ as x ⟶ ∞ and t⟶∞.
The function f ðx, tÞ is called an exponential order as x ⟶ ∞, t ⟶ ∞, and clearly, it does not grow faster than Ke cx+dt x ⟶ ∞, t ⟶ ∞: If a function f ðx, tÞ is continuous in every finite intervals ð0, XÞ and ð0, TÞ and of exponential order e cx+dt , then the double Sumudu transform of f ðx, tÞ exists for all 1/k and 1/p provided that Re ð1/kÞ > c and Re ð 1/pÞ > d.

Some Properties of Double Sumudu Transform
Property 9 (linearity of TLT, [42,43]). If f ðx, tÞ and gðx, tÞ are two functions of x and t such that S x S t f f ðx, tÞg = Fðk, pÞ and S x S t fgðx, tÞg = Gðk, pÞ, then where α and β are constants.

Advances in Mathematical Physics
Property 10 (change of scale property, [44]). If f ðx, tÞ and gðx, tÞ are two functions of x and t such that S x S t f f ðx, tÞg = Fðk, pÞ and S x S t fgðx, tÞg = Gðk, pÞ, then where c and d are nonzero constants.
Proof. from (7), we have By using the integration by substitution method for the integral inside the bracket and the integral outside the bracket, we get Similarly, ☐ Property 12 (second shifting property, [42,43]).
where Hðx, tÞ is the Heaviside unit step function defined by

The New Iterative Method
Daftardar-Gejji and Jafari [33] have considered the following functional equation: where N is a nonlinear operator in a Banach space such that N : B ⟶ B and f is a given element of the Banach space B.
Here, u is assumed to be the solution of equation (23) whose series expansion is given by The nonlinear operation N can then be decomposed as Using Equations (24) and (25), Equation (23) is equivalent to From Equation (26), we define the following recurrence relation: Thus, and hence, Therefore, the n th term approximate solution of equation (24) is given by

Double Sumudu Transform Coupled with Iterative Method
To solve the system of sine-Gordon equation (2), first, we decompose the source function f ðx, tÞ into f 1 ðx, tÞ and f 2 ðx, tÞ. The part f 1 ðx, tÞ with the terms in equation (2) always leads to the simple algebraic expression while applying the inverse double Sumudu transform. The portion f 2 ðx, tÞ is combined with the nonlinear term of equation (2) to avoid noise terms in the iteration process. Similarly, hðx, tÞ = h 1 ðx, tÞ + h 2 ðx, tÞ.
Applying the double Sumudu transform to both sides of the system of equation (2) and by using Table 1 and Definition 4, we get , v x f 1 , and h 1 are Sumudu transforms of functions of two variables.
Multiplying both sides of each equation in system (31) by k 2 p 2 yields Next, applying a single Sumudu transform to the systems of equations (3) and (4), we get u k, 0 Here, φ i and g i (i = 1, 2, 3, 4) are Sumudu transforms of functions of single variables.
By substituting the systems of equations (34) and (33) into the systems of equation (32) and simplifying, we obtain u k, p ð Þ= Applying the inverse double Sumudu transform to the system of equation (35), we obtain Now, we start an iterative process. We assume that is the solution of the coupled nonlinear sine-Gordon equation (2). Substituting the systems of equation (37) into the systems of equation (36), we obtain

Advances in Mathematical Physics
Now, by (25), the nonlinear term sin ðuðx, tÞ − vðx, tÞÞ in the systems of equation (38) is decomposed as Substituting the system of equation (38) into the systems of equation (37), we get Now, we define the recurrence relation from the systems of equation (40) as follows:    Therefore, the solution of the systems of SGE equations (2)-(4) in series form is given by Finally, Theorem 13 will be applied to demonstrate the series solution obtained from the new iterative method which is convergent.

Illustrative Examples
In order to show the validity and effectiveness of the method under consideration, some examples are presented here. Example 1. Consider the following IBVP for onedimensional coupled sine-Gordon equation on Ω = ½0, π, t > 0: with initial condition and boundary condition Solution. Applying the properties of double Sumudu transform to both sides of the system of equation (43), we get Multiplying both sides of each equation in the system Advances in Mathematical Physics (46) by k 2 p 2 , we obtain Next, applying the single Sumudu transform to the systems of equations (44) and (45), we get u k, 0 By substituting the system of equations (49) and (48) into the system of equation (47) and simplifying, we obtain u k, p ð Þ= pk 1 + k 2 − Applying the inverse double Sumudu transform to (50), we get Now, applying the new iterative method to a system of equations (51), we obtain the components of the solution as follows:  For n = 1, in the same way, we obtain and so on. Therefore, the solution of Example 1 in the sense of (42) is which is the same as the result obtained by Yildirim et al. [14].
Let us now test the convergence of the obtained series solution. From the system of equation (51), we have Thus, for all x, t ≥ 0, we have Therefore,     , Similarly, we have for all k ≥ 0 by the principle of mathematical induction. As the conditions of Theorem 13 are satisfied, the series solution obtained by the new iterative method is convergent on the domain of interest.
Example 2. Consider the following IBVP for onedimensional coupled sine-Gordon equation on Ω = ½0, π, t > 0: with initial condition boundary condition Solution. Applying the properties of the double Sumudu transform to both sides of the system of equations (66) , we get Thus, for all x, t ≥ 0, we have Therefore, Therefore, Similarly, we have for all k ≥ 0 by the principle of mathematical induction. As the conditions of Theorem 13 are satisfied, the series solution obtained by the new iterative method is convergent on the domain of interest.

Conclusion
The double Sumudu transform method has been known to be a powerful tool for solving many fractional partial differential equations, integral equations, and so many other equations. In this paper, we have presented a new method, that is, the amalgamation of double Sumudu transform and a numerical method which is known as an iterative method to obtain the exact solution of one dimensional coupled nonlinear sine-Gordon equation. From the solutions of illustrative Examples 1 and 2, we can see that the decomposition From this study, we concluded that double Sumudu transform linked with an iterative method finds quite practical analytical results with less computational work.

Data Availability
No data were used to support the study.

Conflicts of Interest
The author declares that there is no conflict of interest regarding the preparation of this research paper.