Application of Generalized Integral Method (GIRM) to Numerical Evaluations of SolitontoSoliton and SolitontoBottom Interactions
Gantulga Tsedendorj^{1}, Hiroshi Isshiki^{2}, Rinchinbazar Ravsal^{3}
^{1}Department of Mathematics, National University of Mongolia, Ulaanbaatar, Mongolia
^{2}IMA, Institute of Mathematical Analysis, Osaka, Japan
^{3}Department of Administration, National University of Mongolia, Ulaanbaatar, Mongolia
Email address:
To cite this article:
Gantulga Tsedendorj, Hiroshi Isshiki, Rinchinbazar Ravsal. Application of Generalized Integral Method (GIRM) to Numerical Evaluations of SolitontoSoliton and SolitontoBottom Interactions. Applied and Computational Mathematics. Special Issue: Integral Representation Method and its Generalization. Vol. 4, No. 31, 2015, pp. 7886. doi: 10.11648/j.acm.s.2015040301.16
Abstract: Numerical evaluations of solitonsoliton and solitontobottom interaction have many applications in various fields. On the other hand, Generalized Integral Representation Method (GIRM) is known as a convenient numerical method for solving Initial and Boundary Value Problem of differential equations such as advective diffusion. In this work, we apply onestep GIRM to numerical evaluations of propagation of a single soliton, solitontosoliton interaction and solitontobottom interaction. Firstly, in case of a single soliton, the bottom is considered to be constant in order to understand the behavior of the soliton propagation as it travels in the middle of the sea. Next, in case of solitontobottom, we study behavior of a single soliton propagation when the bottom has different geometries. Finally, we evaluate interaction of two different i.e., big and small solitons. To carry out with the studies, we derive and implement GIRM to numerically solve the Kortewegde Vries (KdV) equation. In order to verify the theory, numerical experiments are conducted and accurate approximate solutions are obtained in each case of the soliton interactions.
1. Introduction
The propagation of solitary waves has many applications in various fields in particular, in tsunami studies. Solitary waves or solitons are propagating through a certain medium with constant velocity and describe a variety of nonlinear wave phenomena in one dimension. This phenomenon was first observed, described and reproduced by J.S. Russell [1].
It is wellknown that the nonlinear wave equation of Korteweg and de Vries (KdV) [2] with the appropriate initial conditions admits solitons as its solutions. A modified version of the KdV equation with the variabledepth is studied in Ref [3]. Some derivations of the variabledepth KdV equation that allow a bottom profile are discussed in Ref [3] and [5]. Another approach, spectral and pseudospectral numerical schemes for the KdV equation are discussed in Ref. [4].
In this study, we develop and implement a high precision numerical scheme based on GIRM, in order to study and evaluate a single soliton propagation, solitontobottom and solitontosoliton interactions. GIRM is known as a convenient alternative method for solving Initial and Boundary Value Problem of differential equations such as advective diffusion. In Ref. [6], GIRM is applied to one and twodimensional diffusion problems and twodimensional Burgers’ equation. It is also shown that the determination of a fundamental solution for GIRM is always possible in advance. Effects of Generalized Fundamental Solution (GFS) on GIRM is discussed in Ref. [7]. Further, GIRM is applied to fluid dynamic motion of gas in Ref. [8] and to tidal wave propagation in Ref [9] in order to obtain the accurate numerical solutions. Full computer codes of One and Twostep GIRMs written in widely used computational languages such as Matlab, C, and FORTRAN are given in Ref [10].
2. Theory of 1Step GIRM
2.1. KdV Equation
Let and be the coordinate and time is the surface elevation. The Kortwegde Vries or KdV equation is given by
(1)
with
and (2)
where and are the depth of water and speed of the linear wave, respectively. KdV equation was obtained under the assumption that the water depth is constant. However, if the water depth changes slowly, KdV equation could reflect this change. Hence, we consider and as functions of and i.e. and . The linear wave speed is given by
(3)
with is the gravitational acceleration. If we rewrite Eq. (1) in a more general form, we have
(4)
In Eq. (4), diffusion term is added, where is the diffusion coefficient. The wave form and the speed of a soliton obtained analytically are given by
(5)
(6)
2.2. OneStep GIRM for KdV Equation
First, we derive 1step Generalized Integral Representation (GIR) for Eq. (4). Since KdV equation expresses a wave traveling to the positive direction of , we consider it in the region , where is large enough. We also assume function and its derivatives tend to zero, as tends to and .
Multiplying both sides of Eq. (4) by function of and , and integrating over , we obtain
(7)
Next, transforming all the spatial derivatives of in Eq. (7) and taking into account the boundary conditions, yield
(8)
Finally, rearranging and exchanging variables and , in Eq. (8), we have
(9)
where is a Generalized Fundamental Solution (GFS) chosen properly. The determination of a fundamental solution for GIRM is always possible in advance [610]. In our case, we take the Gaussian GFS:
(10)
Eq. (9) is a GIR of Eq. (4). This integral representation is applied to numerical solution of Eq. (4). If is known in , then Eq. (9) is an integral equation with unknowns in , where is the kernel function of the integral equation. Namely, we can obtain numerically, if we use, for instance, the following procedure:
Let be known at time
Obtain from Eq. (9)
Then
Add to Repeat process. (11)
In the present paper, an iteration is used to obtain (→ Appendix A).
3. Numerical Method and Results of 1Step GIRM
In numerical experiments, we take for simplicity, the coefficients , , and to change slowly with respect to . In this case, Eq. (9) would be approximated by
(12)
To begin with, we introduce a uniform mesh as follows:
(13)
and
,, ,
etc. (14)
We approximate each term in Eq. (12) as follows:
(15a)
(15b)
(15c)
(15d)
Hence Eq. (12) can be discretized as
(16)
where we denote
(17a)
and
(17b)
Eq. (16) is satisfied at the points and the unknowns are , . Hence, we have N equations for N unknowns.
3.1. Propagation of a Single Soliton
First, we numerically evaluate propagation of a single soliton as if it moves in an infinite open sea. Numerical experiment is conducted straightforward by using Eq. (16) along with Eq. (17a) and Eq. (17b). For comparison, we also evaluate it by Finite Difference Method using second order RungeKutta scheme. Numerical results are shown in Fig. 1. The accuracy of the numerical results by 1step GIRM using implicit time evolution is high. Values of the parameters used in the numerical experiments are:
(18)
and the initial condition is
(19)
We provide Matlab code for propagation of a single soliton in Appendix A.
3.2. SolitontoSoliton Interaction
We consider an interaction of big and small solitons. From Eq. (6), the big one moves faster than the smaller one. Hence, the big one catches up the small one. We introduce a frame moving with mean speed of the individual solitons:
,
and , then Eq. (1) becomes
(20)
Therefore, a GIR of Eq. (20) is
(21)
with
and (22)
Fig. 2 shows the interaction of two solitons calculated by 1Step GIRM. Both solitons move to the right and the speed of the bigger one is faster than the speed of the small one (see Eq. (6)). The calculation was conducted in a frame moving with mean speed of the individual solitons. The accuracy of the numerical results obtained by 1step GIRM is very high. Initial condition is given by
(23)
and values of the parameters are
(24)
3.3. SolitontoBottom Interaction
We consider an interaction of a soliton to the bottom with three different geometries: bottom geometry 1.
(25a)
bottom geometry 2
(25b)
bottom geometry 3
(25c)
and initial condition is given by
(26)
Values of parameters used in the experiments are
, , , , , , , , (27a)
for bottom geometries 1 and 2 and
, , , , , , , , (27b)
for bottom geometry 3, respectively.
Interactions of a single soliton to the bottom with different geometries are shown in Fig. 3. The tendency of the numerical results by 1step GIRM seems reasonable.
4. Conclusion
We developed an accurate numerical scheme based on GIRM for solitary wave phenomena including a single soliton propagation, solitontobottom interaction and solitontosoliton interaction. In case of a single soliton, the bottom is considered constant in order to understand the behavior of the soliton propagation as it travels in the middle of the sea. In case of solitontobottom interaction, we study behavior of a single soliton propagation when the bottom has different geometries. Finally, we evaluate interaction of two different i.e., big and small solitons. In order to verify the theory, numerical experiments are conducted and accurate approximate solutions are obtained in admissible time in each case of the soliton interactions.
Appendix A. Matlab Code for Propagation of a Single Soliton

References