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.

1. 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 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 [1012], 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 [1618].

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 [3337]). 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 and .

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.

2. A Brief Introduction of Double Sumudu Transform (DST)

Definition 1 (see [24, 40]). Sumudu transform of a function is defined for all real numbers as the function , given by

Definition 2 (see [30]). The inverse Sumudu transform of a function is denoted by the symbol and is defined by the Bromwich contour integral

Definition 3 (see [24, 40, 41]). Let be a function that can be expressed as a convergent infinite series and let , then the double Sumudu transform of the function in the positive quadrant of the -plane is given by where and are the transform variables for x and t, respectively, whenever the improper integral is convergent.

Definition 4 (see [42, 43]). The inverse double Sumudu transform is defined by the following form:

Remark 5 (see [26, 40]). If 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 is given by Double Sumudu transform for the first partial derivative with respect to is given by Double Sumudu transform for the second partial derivative with respect to is given by Double Sumudu transform for the second partial derivative with respect to is given by

2.1. Convergence Theorem of Double Sumudu Transform

Lemma 6 (see [24]). Let be a function of two variables continuous in the positive quadrant of the -plane. If the integral converges at , , then the integral converges for ,

2.2. Existence and Uniqueness of the Double Sumudu Transform

If is an exponential order, then and as and if there exists a positive constant such that for all then

And we write as Or, equivalently,

And we write this as as .

The function is called an exponential order as , and clearly, it does not grow faster than

Theorem 7 (existence). If a function is continuous in every finite intervals and and of exponential order , then the double Sumudu transform of exists for all provided that and .

For the proof, see [42, 43].

Theorem 8 (uniqueness). Let and be continuous functions defined for and having double Sumudu transforms and , respectively. If , then

2.3. Some Properties of Double Sumudu Transform

Property 9 (linearity of TLT, [42, 43]). If and are two functions of and such that and , then where and are constants.

Property 10 (change of scale property, [44]). If and are two functions of and such that and , then (i)(ii)where and are nonzero constants.

Property 11 (first shifting property). If then

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]). If then where is the Heaviside unit step function defined by

3. The New Iterative Method

Daftardar-Gejji and Jafari [33] have considered the following functional equation: where is a nonlinear operator in a Banach space such that and is a given element of the Banach space . Here, is assumed to be the solution of equation (23) whose series expansion is given by

The nonlinear operation 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 term approximate solution of equation (24) is given by

3.1. Convergence of the New Iterative Method

Theorem 13. If is a continuously differentiable functional in a neighborhood of and for all , then the series is absolutely convergent [33].

4. Double Sumudu Transform Coupled with Iterative Method

To solve the system of sine-Gordon equation (2), first, we decompose the source function into and . The part with the terms in equation (2) always leads to the simple algebraic expression while applying the inverse double Sumudu transform. The portion is combined with the nonlinear term of equation (2) to avoid noise terms in the iteration process. Similarly, .

Applying the double Sumudu transform to both sides of the system of equation (2) and by using Table 1 and Definition 4, we get where , , ,, , , and are Sumudu transforms of functions of two variables.

Multiplying both sides of each equation in system (31) by yields

Next, applying a single Sumudu transform to the systems of equations (3) and (4), we get

Here, and () 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

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

Now, by (25), the nonlinear term 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.

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 following IBVP for one-dimensional coupled sine-Gordon equation on : 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 (46) by , we obtain

Next, applying the single Sumudu transform to the systems of equations (44) and (45), we get

By substituting the system of equations (49) and (48) into the system of equation (47) and simplifying, we obtain

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 , 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 , we have