# Interaction of Solitons for Sine-Gordon-Type Equations

Authors: Georgii A. Omel’yanov; Israel Segundo-Caballero
Journal: Journal of Mathematics
 (2013)

---

## Abstract

The subject of our consideration is a family of semilinear wave equations with a small parameter and nonlinearities which provide the existence of kink-type solutions (solitons). Using asymptotic analysis and numerical simulation, we demonstrate that solitons of the same type (kinks or antikinks) interact in the same manner as for the sine-Gordon equation. However, solitons of the different type preserve the shape after the interaction only in the case of two or three waves, and, moreover, under some additional conditions.


---

## Body

## 1. Introduction

We consider the semilinear wave equation(1)ε2(utt-uxx)+F′(u)=0,x∈ℝ1,t>0,
with some smooth nonlinearities F(u) and the parameter ε→0. It is well known [1] that the unique completely integrable representative of the family (1) is the sine-Gordon equation, that is, (1) with F′(u)=sin(u) (see, e.g., [2, 3]).At the same time, there are many nonlinearitiesF(u) such that (1) admits exact traveling wave solutions of the soliton type:
(2)u(x,t,ε)=ω(±βx-Vtε),β=(1-V2)-1/2,ω(η)∈𝒞∞(ℝ),ω(η)⟶0forη⟶-∞,ω(η)⟶1forη⟶+∞.
Traditionally, solution (2) with the sign “plus” is called “kink,” whereas (2) with the sign “minus” is called “antikink.” It is easy to check that the conditions (A)
F(z)∈𝒞∞(ℝ),F(z)>0 for z∈(0,1),(B)
F(i)(z0)=0,i=0,1,…,k,F(k+1)(z0)>0,wherez0=0 and z0=1, and k=1 or k=3, are sufficient for the existence of kink/antikink solutions such that
(3)|ω′(η)|≤cη-2asη⟶±∞.Moreover, under the periodicity condition(C)
F(z+1)=F(z)any combination of kink-antikink waves, that is,(4)uΣ=∑i=1N‍ω(±βix-Vit-xi0ε),xi+10-xi0>1,0<t≪1,will approximately sufficiently well the exact solution of the corresponding Cauchy problem. This brings up the question about the character of the interaction between the entities (2).The first results of this topic have been obtained in [4–6], where the two-wave asymptotics have been constructed and then simulated numerically. In the present paper, we continue the cited investigations considering the interaction of three or more waves. The point is that there is a hypothesis (Danilov and Subochev [7], Boris Dubrovin’s, private communication) that there are sufficiently many equations with sine-Gordon scenario of two solitary waves interaction, but three waves can interact in the same manner for the completely integrable equations only. In fact, the situation is more complicated. Our main result consists in the conclusion that, in the leading term with respect to ε, two solitary waves (2) interact like sine-Gordon solitons but the stability of three or more waves depends on their parameters. Namely, structures with one kink and two antikinks or with one antikink and two kinks remain stable for special choice of the velocities, whereas interaction of four or more waves does not change the structure only for the entities of the same type (all of them are kinks or antikinks). Our main tool here is the numerical simulation.The contents of the paper is the following. In Section2 we present the asymptotics for the interaction of two solitary waves; the description of the finite difference scheme is contained in Section 3; in Section 4 we consider the numerical results.


## 2. Asymptotic Solution

For essentially nonintegrable interaction problems it is impossible to construct either explicit solutions (classical or weak) or asymptotics in the classical sense. However, it is possible to construct an asymptotic solution in the weak sense (see, e.g., [4, 5, 8–11] and references therein). The main advantage of this approach is the possibility of reducing the problem of describing nonlinear waves interaction to a qualitative analysis of some ordinary differential equations (instead of partial differential equations). This method takes into account the fact that kinks (as well as solitons [9, 10]) which are smooth for ε>0 become nonsmooth in the limit as ε→0. So it is possible to treat such solutions as a mapping 𝒞∞((0,T);𝒞∞(ℝx1)) for ε=const>0 and only as 𝒞((0,T);𝒟′(ℝx1)) uniformly in ε≥0. Accordingly, the remainder should be small in the weak sense. This rather trivial observation allowed to reach a progress for some old problems about nonlinear wave interaction for nonintegrable equations. As for the equations of the form (1), it should be noted that there is an obstacle to apply the standard 𝒟′ construction. Indeed, in the 𝒟′ sense, the differential terms of (1) are subordinated to the nonlinear term. Moreover, the left-hand side of (1) is of the value O(ε2) in the weak sense for any u of the form (4) and t≪1. Obviously, this prevents the construction of the correct asymptotics for the Cauchy problem. To overcome this obstacle, in [4] has been constructed a new definition of asymptotic solutions, which involves in the leading term the derivatives of u with arguments x/ε and t/ε.Definition 1.
A sequenceu(t,x,ε), belonging to 𝒞∞((0,T);𝒞∞(ℝx1)) for ε=const>0 and belonging to 𝒞((0,T);𝒟′(ℝx1)) uniformly in ε, is called a weak asymptotic mod O𝒟′(ε2) solution of (1) if the relation
(5)2ddt∫-∞∞‍ε2utuxψdx+∫-∞∞‍{(εut)2+(εux)2-2F(u)}ψxdx=O(ε2)
holds uniformly in t for any test function ψ=ψ(x)∈𝒟(ℝ1).Here the right-hand side is a𝒞∞-function for ε=const>0 and a piecewise continuous function uniformly in ε≥0. The estimate is understood in the 𝒞(0,T) sense:
(6)g(t,ε)=O(εk)⟷maxt∈[0,T]|g(t,ε)|≤cεk.
The left-hand side of (5) is the result of multiplication of (1) by ψ(x)ux and integration by parts in the case of smooth u. Therefore, the relation (5) is satisfied automatically for any exact solution. On the other hand, the relation (5) is just the orthogonality condition that appears for single-phase asymptotics [12, 13]. This condition both guarantees the first correction existence and allows to find an equation for the distorted kink’s front motion.Definition 2.
A functionv(t,x,ε) is said to be of the value O𝒟′(εk) if the relation
(7)∫-∞∞‍v(t,x,ε)ψ(x)dx=O(εk)
holds uniformly in t for any test function ψ∈𝒟(ℝx1).Let us consider the interaction of two kinks,(8)u|t=0=∑i=12‍ω(βix-xi0ε),ε∂u∂t|t=0=-∑i=12‍βiViω′(βix-xi0ε),
where βi=1/1-Vi2, |Vi|∈(0,1), and the initial front positions xi0 are such that x20-x10>1. Obviously, it is assumed that the trajectories x=Vit+xi0 have a joint point x=x* at a time instant t=t*.The asymptotic ansatz for the problem (1), (8) has the following form:
(9)u=∑i=12‍{ω(βix-Φi(t,τ,ε)ε)+Ai(τ)U(βix-Φi(t,τ,ε)μ2ε)}.
Here Φi=ϕi0(t)+εϕi1(τ),ϕi0=Vit+xi0 are the trajectories of the noninteracting kinks, τ=ψ0(t)/ε denotes the “fast time,” ψ0(t)=ϕ20(t)-ϕ10(t). The phase corrections ϕi1 are smooth functions such that
(10)ϕi1⟶0asτ⟶-∞,ϕi1⟶ϕi1∞=constiasτ⟶+∞,
with a rate not less than 1/|τ|. Furthermore, Ai(τ)∈𝒞∞ are exponentially vanishing as |τ|→∞ functions, μ is a sufficiently small parameter, ε<μ≪1, and
(11)U(η)=dmU0(η)dηm,
where m≥1 is an arbitrary number and U0(η)∈𝒞∞ is a sufficiently fast vanishing function as |η|→∞.The main result, which is known for the problem (1), (8), is the following.Theorem 3 (see [4]).
Let the assumptions (A)–(C) hold. Set the additional assumptions(D)
F(1/2+z)=F(1/2-z),(E)
let the functionF(z) be such that the inequality
(12)∫-∞∞‍F(ω(η)+ω(θη))dη≤∫-∞∞‍{F(ω(η))+F(ω(θη))}2dη
holds uniformly inθ∈(0,∞). Then the interaction of kinks in the problem (1), (8) preserves the sine-Gordon scenario with accuracy O𝒟′(ε2) in the sense of Definition 1. The weak asymptotic solution of (1), (8) has the form (9) with a special choice of the amplitudes Ai and of the parameter μ.Remark 4.
The symmetry (D) has been assumed to simplify the asymptotic analysis and it is not very important.Remark 5.
The sense of assumption (E) is the following. The phase correctionsϕi1 are solutions of a 2×2-dynamical system with a singularity whose support divides the phase plane into two parts with the possible exception of the point (0,0). Assumptions (10) are satisfied (consequently, the sine-Gordon scenario takes place) if and only if there exists a specific trajectory which goes from one half-plane to the other one through the point (0,0). When Ai in (9) are equal to zero, the existence of the trajectory implies the appearance of an additional very complicated assumption. This condition can be made more coarse and transformed into the simplest form (12). Such version can be treated as an admissible one since it is satisfied for the sine-Gordon equation for any velocities Vi,i=1,2. The same is true for the nonlinearity
(13)F(u)=sin4(πu).
Taking into account a freedom in the choice of the amplitudes Ai,i=1,2, assumption (12) can be made weaker. However, the dynamical system with Ai≠0,i=1,2, is very complicated and its complete analysis remains undone.Obviously, all stated above remains true for the antikink-antikink interaction.Let us focus our attention in the kink-antikink interaction, that is, in (1) with initial data
(14)u|t=0=∑i=12‍ω(Siβix-xi0ε),ε∂u∂t|t=0=-∑i=12‍SiβiViω′(Siβix-xi0ε),
where S1=1,S2=-1, and the notation βi, Vi, xi0 is the same as in (8).The asymptotic ansatz for the solution of the problem (1), (14) differs a little bit from (9), namely,
(15)u=∑i=12‍{ω(Siβix-Φi(t,τ,ε)ε)+Ai(τ)U(Siμβix-Φi(t,τ,ε)ε)},
with the same notation and assumption (10).Technically, the construction of (15) is similar to the kink-kink case. However, the resulting dynamical system for the phase corrections becomes much more complicated. Moreover, it is impossible to simplify the additional assumption, which appears here also, without loss of the adequacy. For this reason we do not present the explicit form of the additional condition but state only the existence of the weak asymptotics (15) under some restrictions for V1 and V2. We refer the readers to [5] for the explicit statement.Finally we note that there is a correspondence between weak asymptotic solutions and energy relations for (1).Theorem 6.
Let the assumptions of Theorem3 hold. Then two kinks (9) preserve mod O𝒟′(ε2) their forms after the interaction if and only if they satisfy the conservation law
(16)ddt∫-∞∞‍utuxdx=0
and the energy relation
(17)2ddt∫-∞∞‍xε2utuxdx+∫-∞∞{(εut)2+(εux)2-2F(u)}dx=0.‍Similar conclusion is true for the kink-antikink pair (15) [5].


## 3. Finite Differences Scheme

The actual numerical simulation for (1) is realized for a finite x-interval, x∈[0,L]. For this reason we simulate the Cauchy problem by the following mixed problem:
(18)ε2(utt-uxx)+F′(u)=0,x∈(0,L),t∈(0,T),u|x=0=νℓ,u|x=L=νr,u|t=0=u0(xε),ε∂u∂t|t=0=u1(xε),
where u0 is a combination of kinks and antikinks of the form (4) and u1 denotes its time derivative calculated at t=0, νℓ=u0|x=0,νr=u0|x=L. To simulate by (18) the interaction phenomena, we assume that L, T, and the initial front positions xi0,i=1,2, are such that the intersection point of the solitary wave fronts belongs to QT=(0,T)×(0,L). Furthermore, let L, T, and xi0 be such that uniformly in t≤T(19)|uΣ|x∈[0,δ]-νℓ|≤cε2,|uΣ|x∈[L-δ,L]-νr|≤cε2
for some sufficiently small δ>0. Since it is impossible to create any finite difference scheme for the problem (18), which remains stable uniformly in ε→0 and t∈(0,T),T=const, we will treat ε as a small but fixed constant. However, we will fix any relation between ε and finite differences scheme parameters.To create a finite differences scheme for (18) we should choose appropriate approximations for the differential terms and for the nonlinear term. Let us do it separately.


### 3.1. Preliminary Nonlinear “Scheme”

As usually, we define a meshQT,τ,h={(xi,tj)=(ih,jτ),i=0,…,I,j=0,…,J} over QT and denote
(20)yij=u(xi,tj),yitj=yij+1-yijτ,yit-j=yij-yij-1τ,yixj=yi+1j-yijh,yix-j=yij-yi-1jh,yitt-j=(yitj)t-,yixx-j=(yixj)x-.
Let us consider the following system of nonlinear equations:
(21)ε2(yitt-j-yixx-j+1)+F′(yij+1)=0,i=1,…,I-1,j=2,3,…,y0j=νℓ,yIj=νr,j=0,1,…,yi0=u0(xiε),εyit0=u~1(xiε,τ),i=0,…,I,
where u~1(xi/ε,τ) is such that last equality in (18) is approximated with accuracy O(τ2). Obviously, the local approximation accuracy of (21) is O(τ2+h2).To simplify the notation, we will write(22)y:=yij,y^:=yij+1,yˇ:=yij-1.So the short form of (21) is the following:
(23)ε2(ytt--y^xx-)+F′(y^)=0.Our first result consists in obtaining the boundedness condition for the problem (21) solution.Lemma 7 (see [6]).
Letε be a sufficiently small constant and let
(24)τε2≤const.
Suppose that system (21) is solvable for any j=2,…,J. Then uniformly in j(25)∥εyt∥2+∥εy^x∥2+2∥F(y^)∥2+τε2{|∥ε2ytt-∥|2(j)+|∥ε2yxx-∥|2(j)}≤{∥εyt0∥2+∥εyx1∥2+2∥F(y1)∥2}ectjτ/ε2(1+O(τ)),
where ∥·∥ and |∥·∥|(j) are the ℒ2 norms, namely,
(26)∥f∥2=h∑i=1I-1‍|fi|2,|∥f∥|2(j)=τ∑k=1j‍∥fk∥2.
Here and in what follows c denotes a const > 0 which does not depend on h, τ, and ε.As a consequence of this lemma and the identity(27)yij=yi0+τ∑k=0j-1‍yitk,
we obtain the inequality
(28)∥yj∥2≤2∥y0∥2+2tjτ∑k=0j-1‍∥ytk∥2≤2∥y0∥2+2tj2ε2c0,
where c0>0 denotes the right-hand side in (25). Obviously, this estimate is very rough. However, it can be improved a little for the specific initial data (8) and (14).Lemma 8 (see [6]).
Let the assumptions of Lemma7 be satisfied. Then for the initial data u0(xi/ε), u~1(xi/ε,τ), which approximate the Cauchy data (8) or (14), the following estimate holds uniformly in j:
(29)ε{∥ytj∥+∥yxj∥+∥yj∥}≤c.


### 3.2. Linearization

Now we should verify the solvability of (21) for any fixed j≥1, that is, of the equation
(30)yj+1-τ2yxx-j+1+τ2ε2F′(yj+1)=Gj,Gj=yj+τyt-j,
as well as select a way to linearize the nonlinearity. To this aim let us construct the sequence of functions φ(s):={φ0(s),…,φI(s)}, s≥0, such that φ(0)=yj and φ(s) for s≥1 satisfies the equation
(31)φ(s)-τ2φxx-(s)+τ2ε2{F′(φ(s-1))+F′′(φ(s-1))×(φ(s)-φ(s-1))F′(φ(s-1))+F′′(φ(s-1))}=Gj,φ0(s)=νℓ,φI(s)=νr.The solvability of the algebraic system (31) is obvious for sufficiently small τ and τ/ε2≤const. To simplify the notation we write φ:=φ(s),φ-:=φ(s-1),φ--:=φ(s-2). Let also define
(32)w:=φ-φ-,w-:=φ--φ--.Lemma 9 (see [6]).
Let assumption (24) be satisfied and let τ be sufficiently small. Then
(33)∥φ∥2+τ2ε2∥εφx∥2≤(1+cτ){∥yj∥2+cτ(∥yj∥2+∥εytj∥2)}+c{∥w∥2+τ2∥wx∥2}2,(34)g(1)≤cτ3/2,g(s)≤cτg2(s-1)fors>1,
where
(35)g(s):=∥w(s)∥2+τ2∥wx(s)∥2.Thus we immediately conclude that the terms of thew-sequence vanish very rapidly,
(36)∥w(1)∥2≤cτ3/2,∥w(2)∥2≤cτ4,∥w(3)∥2≤cτ9,….
By virtue of (33), the terms of the φ-sequence are bounded uniformly in s(37)∥φ(s)∥2≤∥yj∥2(1+O(τ)).
Furthermore, for any n>0(38)∥φ(s+n)-φ(s)∥≤∑i=1n‍∥ws+i∥≤∥ws+1∥∑i=1∞‍∥ws+i∥∥ws+1∥≤c∥ws+1∥.
This implies the main statement of this subsection.Theorem 10 (see [6]).
Let assumption (24) be satisfied and ε = const. Then for sufficiently small τ the sequence φ converges in the ℒh2 sense to the solution of (30). Moreover,
(39)∥yj+1-φ(2)∥≤cτ9/2.


### 3.3. Algorithm for the Numerical Simulation

Since the accuracyO(τ9/2) is much less than the accuracy of the finite differences scheme (21), we obtain the following algorithm for the numerical simulation of the problem (18) solution:For any fixedj=1,2,…,[T/τ],T=const, we(i)
defineφ(0):=yj,(ii)
calculateφ(s),s=1,2, accordingly with (31),(iii)
defineyj+1:=φ(2), redefine j:=j+1, and go back to (i).By virtue of the estimates (29) and (39), this algorithm allows to calculate bounded in ℒ2(QT,h,τ) numerical solution of the problem (18).Note that this result can be improved. Moreover, it turns out that the algorithm is absolutely stable. To prove this we state firstly the proposition.Lemma 11 (see [6]).
Let assumption (24) be satisfied and ε = const. Then uniformly in s(40)∥εφt(s)∥≤const,∥εφx(s)∥≤const.
Moreover, uniformly in j(41)∥ε2yxtj∥+∥ε2yxx-j+1∥≤cε,∥εytj+1-εφt(s)∥≤cεp(s),∥εyxj+1-εφx(s)∥≤cεp(s),
with some p(s) that tend to infinity as s→∞.An immediate consequence of the Lemmas7–11 is the following.Theorem 12 (see [6]).
Let the assumption (24) be satisfied and ε = const. Then the solution of the above described finite differences scheme converges to the solution of the problem (18) as τ,h→0 in the W21(QT) sense.Finally, in view of the boundedness of the sequenceφj(s), s=1,2,j=2,3,…, it is easy to establish our last statement.Theorem 13 (see [6]).
Under the assumptions of Theorem12 the above described finite differences scheme is stable in the W21(ΩT,τ,h) sense.


## 3.1. Preliminary Nonlinear “Scheme”

As usually, we define a meshQT,τ,h={(xi,tj)=(ih,jτ),i=0,…,I,j=0,…,J} over QT and denote
(20)yij=u(xi,tj),yitj=yij+1-yijτ,yit-j=yij-yij-1τ,yixj=yi+1j-yijh,yix-j=yij-yi-1jh,yitt-j=(yitj)t-,yixx-j=(yixj)x-.
Let us consider the following system of nonlinear equations:
(21)ε2(yitt-j-yixx-j+1)+F′(yij+1)=0,i=1,…,I-1,j=2,3,…,y0j=νℓ,yIj=νr,j=0,1,…,yi0=u0(xiε),εyit0=u~1(xiε,τ),i=0,…,I,
where u~1(xi/ε,τ) is such that last equality in (18) is approximated with accuracy O(τ2). Obviously, the local approximation accuracy of (21) is O(τ2+h2).To simplify the notation, we will write(22)y:=yij,y^:=yij+1,yˇ:=yij-1.So the short form of (21) is the following:
(23)ε2(ytt--y^xx-)+F′(y^)=0.Our first result consists in obtaining the boundedness condition for the problem (21) solution.Lemma 7 (see [6]).
Letε be a sufficiently small constant and let
(24)τε2≤const.
Suppose that system (21) is solvable for any j=2,…,J. Then uniformly in j(25)∥εyt∥2+∥εy^x∥2+2∥F(y^)∥2+τε2{|∥ε2ytt-∥|2(j)+|∥ε2yxx-∥|2(j)}≤{∥εyt0∥2+∥εyx1∥2+2∥F(y1)∥2}ectjτ/ε2(1+O(τ)),
where ∥·∥ and |∥·∥|(j) are the ℒ2 norms, namely,
(26)∥f∥2=h∑i=1I-1‍|fi|2,|∥f∥|2(j)=τ∑k=1j‍∥fk∥2.
Here and in what follows c denotes a const > 0 which does not depend on h, τ, and ε.As a consequence of this lemma and the identity(27)yij=yi0+τ∑k=0j-1‍yitk,
we obtain the inequality
(28)∥yj∥2≤2∥y0∥2+2tjτ∑k=0j-1‍∥ytk∥2≤2∥y0∥2+2tj2ε2c0,
where c0>0 denotes the right-hand side in (25). Obviously, this estimate is very rough. However, it can be improved a little for the specific initial data (8) and (14).Lemma 8 (see [6]).
Let the assumptions of Lemma7 be satisfied. Then for the initial data u0(xi/ε), u~1(xi/ε,τ), which approximate the Cauchy data (8) or (14), the following estimate holds uniformly in j:
(29)ε{∥ytj∥+∥yxj∥+∥yj∥}≤c.


## 3.2. Linearization

Now we should verify the solvability of (21) for any fixed j≥1, that is, of the equation
(30)yj+1-τ2yxx-j+1+τ2ε2F′(yj+1)=Gj,Gj=yj+τyt-j,
as well as select a way to linearize the nonlinearity. To this aim let us construct the sequence of functions φ(s):={φ0(s),…,φI(s)}, s≥0, such that φ(0)=yj and φ(s) for s≥1 satisfies the equation
(31)φ(s)-τ2φxx-(s)+τ2ε2{F′(φ(s-1))+F′′(φ(s-1))×(φ(s)-φ(s-1))F′(φ(s-1))+F′′(φ(s-1))}=Gj,φ0(s)=νℓ,φI(s)=νr.The solvability of the algebraic system (31) is obvious for sufficiently small τ and τ/ε2≤const. To simplify the notation we write φ:=φ(s),φ-:=φ(s-1),φ--:=φ(s-2). Let also define
(32)w:=φ-φ-,w-:=φ--φ--.Lemma 9 (see [6]).
Let assumption (24) be satisfied and let τ be sufficiently small. Then
(33)∥φ∥2+τ2ε2∥εφx∥2≤(1+cτ){∥yj∥2+cτ(∥yj∥2+∥εytj∥2)}+c{∥w∥2+τ2∥wx∥2}2,(34)g(1)≤cτ3/2,g(s)≤cτg2(s-1)fors>1,
where
(35)g(s):=∥w(s)∥2+τ2∥wx(s)∥2.Thus we immediately conclude that the terms of thew-sequence vanish very rapidly,
(36)∥w(1)∥2≤cτ3/2,∥w(2)∥2≤cτ4,∥w(3)∥2≤cτ9,….
By virtue of (33), the terms of the φ-sequence are bounded uniformly in s(37)∥φ(s)∥2≤∥yj∥2(1+O(τ)).
Furthermore, for any n>0(38)∥φ(s+n)-φ(s)∥≤∑i=1n‍∥ws+i∥≤∥ws+1∥∑i=1∞‍∥ws+i∥∥ws+1∥≤c∥ws+1∥.
This implies the main statement of this subsection.Theorem 10 (see [6]).
Let assumption (24) be satisfied and ε = const. Then for sufficiently small τ the sequence φ converges in the ℒh2 sense to the solution of (30). Moreover,
(39)∥yj+1-φ(2)∥≤cτ9/2.


## 3.3. Algorithm for the Numerical Simulation

Since the accuracyO(τ9/2) is much less than the accuracy of the finite differences scheme (21), we obtain the following algorithm for the numerical simulation of the problem (18) solution:For any fixedj=1,2,…,[T/τ],T=const, we(i)
defineφ(0):=yj,(ii)
calculateφ(s),s=1,2, accordingly with (31),(iii)
defineyj+1:=φ(2), redefine j:=j+1, and go back to (i).By virtue of the estimates (29) and (39), this algorithm allows to calculate bounded in ℒ2(QT,h,τ) numerical solution of the problem (18).Note that this result can be improved. Moreover, it turns out that the algorithm is absolutely stable. To prove this we state firstly the proposition.Lemma 11 (see [6]).
Let assumption (24) be satisfied and ε = const. Then uniformly in s(40)∥εφt(s)∥≤const,∥εφx(s)∥≤const.
Moreover, uniformly in j(41)∥ε2yxtj∥+∥ε2yxx-j+1∥≤cε,∥εytj+1-εφt(s)∥≤cεp(s),∥εyxj+1-εφx(s)∥≤cεp(s),
with some p(s) that tend to infinity as s→∞.An immediate consequence of the Lemmas7–11 is the following.Theorem 12 (see [6]).
Let the assumption (24) be satisfied and ε = const. Then the solution of the above described finite differences scheme converges to the solution of the problem (18) as τ,h→0 in the W21(QT) sense.Finally, in view of the boundedness of the sequenceφj(s), s=1,2,j=2,3,…, it is easy to establish our last statement.Theorem 13 (see [6]).
Under the assumptions of Theorem12 the above described finite differences scheme is stable in the W21(ΩT,τ,h) sense.


## 4. Results of Numerical Simulation

The numerical algorithm has been realized as a program and tested using the sine-Gordon equation in the cases of one, two, three, and four solitary waves. Next we used the nonlinearities (13) and the following:
(42)F(u)=2-cos(2πu)-cos(4πu),(43)F(u)=1+sin(2π(u+54))cos(4π(u+54)).
Obviously, these functions satisfy the assumptions (A)–(D). Note also that functions (13), (42), and (43) for u∈(0,1) have one, three, and five critical points respectively. At the same time, for (13) F(i)(0)=0 if i=0,…,3 and F(4)(0)>0. Therefore, the explicit kink type solution for this case, that is,
(44)ω(η)=1πarccot(-2πη),
tends to zero as η-1 when η→-∞. Conversely, for functions (42), and (43) F(i)(0)=0 if i=0,1 and F(2)(0)>0. Therefore, the kink type solutions for these cases tend to zero with exponential rates when η→-∞. However, their explicit form remains unknown; that is, why we simulate them solving numerically the Cauchy problem
(45)dωdη=2F(ω),η>0,ω|η=0=12
and, by virtue of the condition (C), define ω with negative argument as ω(η)=1-ω(-η). To calculate the solution of (45) we use the Runge-Kutta method of the forth order with the mesh step hη=0.01.


### 4.1. Two-Wave Interaction

In accordance with the asymptotic analysis, two solitons interact preserving the shape; see Figures1 and 2 for the nonlinearity (13) and ε=0.1 (here and in what follows we numerate the waves from the left to the right). Moreover, this remains true independently of the wave parameters V1, V2 for all nonlinearities under consideration. Thus we can conclude that the additional conditions, which appear for the asymptotics, are restrictions of the asymptotic method only.Figure 1
Evolution of the kink-kink pair,V1>0, V2<0.Figure 2
Evolution of the kink-antikink pair,V1>0, V2<0.


### 4.2. Three-Wave Interaction

There are 32 combinations of three solitons (kinks and antikinks) with trajectories which intersect at one instant of time. In view of the symmetryx→-x,t→-t for (1), we can reduce this number to 16, considering only admissible combinations of three kinks and of two kinks and one antikink. Next we reduce the number of combinations to 12 using the symmetry t→-t. All these combinations have been analyzed numerically with the same result for each of the nonlinearities (13), (42), (43): three kinks preserve the shape after the interaction for all admissible velocities, whereas there are only two stable combinations for two kinks and one antikink. The corresponding numerical results are depicted in Figures 3–5 for (43) with ε=0.1. In Figure 6 we present the stable structure of one kink and two antikinks which is dual to Figure 4 in the sense the symmetry x→-x,t→-t. So that, there are 12 stable structures of kinks (K) and antikinks (A):(1)
KKK: (V1>0, V2>0, V3>0), (V1>0, V2>0, V3<0), (V1>0, V2<0, V3<0), (V1<0, V2<0, V3<0);(2)
KKA: (V1>0, V2<0, V3<0);(3)
KAA: (V1>0, V2>0, V3<0);(4)
AAK: (V1>0, V2<0, V3<0);(5)
AKK: (V1>0, V2>0, V3<0);(6)
AAA: (V1>0, V2>0, V3>0), (V1>0, V2>0, V3<0), (V1>0, V2<0, V3<0), (V1<0, V2<0, V3<0).Figure 3
Evolution of the kink-triplet,V1>0, V2<0, V3<0.Figure 4
Evolution of the kink-kink-antikink triplet,V1>0, V2<0, V3<0.Figure 5
Evolution of the antikink-kink-kink triplet,V1>0, V2>0, V3<0.Figure 6
Evolution of the antikink-antikink-kink triplet,V1>0, V2<0, V3<0.


### 4.3. Interaction of Four or More Waves

We obtain again that the interaction of waves of the same type does not destroy the structure; see, for example, Figure7 with 6 kinks, Vi>0,i=1,2,3, and Vj<0, j=4,5,6. As for waves of different types, we did not find any stable combination. In fact, even for four waves the number of possible combinations is too large to simulate each of them, so we checked a part of the combinations only. However, we guess that there are no stable structures for the following reason: each combination of four waves can be considered either as a union of a triplet with an additional wave on its left, or as a union of a triplet with an additional wave on its right. It turns out that one of these triplets should be unstable since the list of stable triplets is too scanty. Thus the total combination is unstable. We illustrate this in Figure 8 for the interaction of two kinks and two antikinks. Indeed, such combination can be treated as the union of the stable triplet depicted in the Figure 4 and the antikink on his right. On the contrary, the same combination is the union of the kink and the kink-antikink-antikink triplet with Vi<0,i=1,2,3, which is unstable.Figure 7
Interaction of 6 kinks for the nonlinearity (43).Figure 8
Evolution of four waves for the nonlinearity (43), V1>0, V2<0, V3<0, V4<0.


## 4.1. Two-Wave Interaction

In accordance with the asymptotic analysis, two solitons interact preserving the shape; see Figures1 and 2 for the nonlinearity (13) and ε=0.1 (here and in what follows we numerate the waves from the left to the right). Moreover, this remains true independently of the wave parameters V1, V2 for all nonlinearities under consideration. Thus we can conclude that the additional conditions, which appear for the asymptotics, are restrictions of the asymptotic method only.Figure 1
Evolution of the kink-kink pair,V1>0, V2<0.Figure 2
Evolution of the kink-antikink pair,V1>0, V2<0.


## 4.2. Three-Wave Interaction

There are 32 combinations of three solitons (kinks and antikinks) with trajectories which intersect at one instant of time. In view of the symmetryx→-x,t→-t for (1), we can reduce this number to 16, considering only admissible combinations of three kinks and of two kinks and one antikink. Next we reduce the number of combinations to 12 using the symmetry t→-t. All these combinations have been analyzed numerically with the same result for each of the nonlinearities (13), (42), (43): three kinks preserve the shape after the interaction for all admissible velocities, whereas there are only two stable combinations for two kinks and one antikink. The corresponding numerical results are depicted in Figures 3–5 for (43) with ε=0.1. In Figure 6 we present the stable structure of one kink and two antikinks which is dual to Figure 4 in the sense the symmetry x→-x,t→-t. So that, there are 12 stable structures of kinks (K) and antikinks (A):(1)
KKK: (V1>0, V2>0, V3>0), (V1>0, V2>0, V3<0), (V1>0, V2<0, V3<0), (V1<0, V2<0, V3<0);(2)
KKA: (V1>0, V2<0, V3<0);(3)
KAA: (V1>0, V2>0, V3<0);(4)
AAK: (V1>0, V2<0, V3<0);(5)
AKK: (V1>0, V2>0, V3<0);(6)
AAA: (V1>0, V2>0, V3>0), (V1>0, V2>0, V3<0), (V1>0, V2<0, V3<0), (V1<0, V2<0, V3<0).Figure 3
Evolution of the kink-triplet,V1>0, V2<0, V3<0.Figure 4
Evolution of the kink-kink-antikink triplet,V1>0, V2<0, V3<0.Figure 5
Evolution of the antikink-kink-kink triplet,V1>0, V2>0, V3<0.Figure 6
Evolution of the antikink-antikink-kink triplet,V1>0, V2<0, V3<0.


## 4.3. Interaction of Four or More Waves

We obtain again that the interaction of waves of the same type does not destroy the structure; see, for example, Figure7 with 6 kinks, Vi>0,i=1,2,3, and Vj<0, j=4,5,6. As for waves of different types, we did not find any stable combination. In fact, even for four waves the number of possible combinations is too large to simulate each of them, so we checked a part of the combinations only. However, we guess that there are no stable structures for the following reason: each combination of four waves can be considered either as a union of a triplet with an additional wave on its left, or as a union of a triplet with an additional wave on its right. It turns out that one of these triplets should be unstable since the list of stable triplets is too scanty. Thus the total combination is unstable. We illustrate this in Figure 8 for the interaction of two kinks and two antikinks. Indeed, such combination can be treated as the union of the stable triplet depicted in the Figure 4 and the antikink on his right. On the contrary, the same combination is the union of the kink and the kink-antikink-antikink triplet with Vi<0,i=1,2,3, which is unstable.Figure 7
Interaction of 6 kinks for the nonlinearity (43).Figure 8
Evolution of four waves for the nonlinearity (43), V1>0, V2<0, V3<0, V4<0.


## 5. Conclusion

Summarizing all stated above, we can deduce that there exists a family of nonlinearities such that kink-kink and kink–antikink pairs preserve the sine-Gordon scenario of interaction at least in the leading term in the asymptotic sense. Apparently, this family can be specified by assumptions (A)–(D).As for multiwave interactions, the situation is more complicated. Apparently, a sufficiently large number of solitons of the same type interact preserving the shape, whereas there are only four stable combinations of kinks and antikinks. In fact, this is rather unexpected, since single kink and antikink waves have the same properties. The second strange phenomenon is that there does not appear any perturbation of the radiation type, comparable with the collision of solitons for the KdV-type equations [11]. Let us note also that our algorithm (see Section 3) allows to eliminate “radiation due to discreteness effects” which appears for the trivial linearization; see [14].



---

*Source: 845926-2013-03-20.xml*