Advances in Mathematical Physics
Volume 2018 (2018), Article ID 3586514, 17 pages
https://doi.org/10.1155/2018/3586514

Research Article

A Newton-Type Approach to Approximate Travelling Wave Solutions of a Schrödinger-Benjamin-Ono System

Juan Carlos Muñoz Grajales

Departamento de Matemáticas, Universidad del Valle, Calle 13 Nro. 100-00, Cali, Colombia

Correspondence should be addressed to Juan Carlos Muñoz Grajales; jcarlmz@yahoo.com

Received 1 October 2017; Revised 1 March 2018; Accepted 13 March 2018; Published 19 April 2018

Academic Editor: Remi Léandre

Copyright © 2018 Juan Carlos Muñoz Grajales. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

We introduce a Newton’s iterative method to approximate periodic and nonperiodic travelling wave solutions of the Schrödinger-Benjamin-Ono system derived by M. Funakoshi and M. Oikawa. We analyze numerically the influence of the model’s parameters on these solutions and illustrate the collision of two unequal-amplitude solitary waves propagating with different speeds computed by using the proposed numerical scheme.

1. Introduction

This paper is concerned with the nonlinear one-dimensional system (hereafter called the Schrödinger-Benjamin-Ono system (SBO))

for , with periodic spatial boundary conditions. This nonlinear dispersive system was derived by Funakoshi and Oikawa [1] to describe the motion of two fluids with different densities under capillary-gravity waves in a deep water regime. The functions and denote the short and long wave terms, respectively. The parameters , are real numbers, and denotes the Hilbert transform

where stands for the integration in the principal value sense. Other physical scenarios, where system (1) with arises, are the sonic-Langmuir wave interaction in plasma physics (Karpman [2]), the capillary-gravity interaction waves (Djordjevic and Redekopp [3], Grimshaw [4]), and water-wave interaction in a nonlinear medium (see, e.g., [5, 6]).

Due to a balance between dispersive and nonlinear effects, system (1) admits the so-called travelling-wave solutions in the form

where and are periodic real-valued functions or smooth real-valued functions such that, for each , and , as . In this last case, these solutions are called solitary waves. In the literature there are some previous analytical results about periodic and nonperiodic travelling wave solutions to the SBO system. For instance, Angulo and Montenegro [7] established the existence of even solitary wave solutions by employing the concentration compactness method (Lions [8, 9]). Existence and stability of a new family of solitary waves was established in [10] for system (1) and a coupled Schrödinger-KdV model. On the other hand, when , the nonperiodic initial value problem corresponding to the SBO system has been considered by Bekiranov et al. [11], who proved well-posedness in the Sobolev space , with . In the case of , Pecher [12] showed local well-posedness for , and Angulo et al. [13] proved global well-posedness for and . However, there are only a few results known for the spatial periodic problem. For instance, assuming that , Angulo et al. [13] showed that system (1) is locally well posed in the Sobolev space for and they further established existence of a smooth branch of periodic travelling wave solutions of the SBO system for and close to 0, and , where is the fundamental period.

In this paper, we introduce a numerical scheme for approximating the periodic and nonperiodic travelling wave solutions of system (1), with whose existence was established in previous works ([7, 10, 13]). This scheme combines a Fourier spectral discretization together with a Newton-type iteration which is initialized by means of well-known exact solutions for the case . By using the spectral Crank–Nicolson numerical scheme introduced in [14] for solving system (1), we validate the approximate travelling waves computed and also illustrate the phenomenon of collision of two right-running unequal-amplitude solitary waves propagating with different speeds. The behavior of the colliding solitary waves is found to be similar to solitary wave solutions of other nonintegrable dispersive-type equations, such as the scalar Benjamin equation [15] and the regularized Benjamin-Ono equation [16]. Such numerical study on the travelling wave solutions of the SBO system has not been considered in previous works to the best knowledge of the author. Furthermore, we obtain numerically new periodic travelling wave solutions which are not included in the analytical theory presented in the work by Angulo et al. [13]. These are the main contributions of the present study. Other contribution is to explain the application of a numerical scheme which belongs to the spectral collocation methods to approximate the solutions of integrodifferential equations, such as system (1). These numerical methods have been intensively used in the solution of dispersive-type equations in the last years. We point out that exact solutions for system (1) are not known when . Therefore, a numerical strategy is very important in order to investigate the properties of the solution space, such as to establish the parameter regime for existence of periodic and nonperiodic travelling waves, orbital stability under small initial disturbances, and interactions among these solutions, for example. We point out that the constant in the SBO system is positive for physical meaning [1]. However, as we mentioned above, in previous works [7, 10–12], the SBO system has been studied from a mathematical point of view, establishing well-posedness and that travelling wave solutions for the SBO system are possible when the parameter takes both signs. This is an interesting mathematical fact, and thus we will conduct some numerical experiments which illustrate the family of travelling wave solutions for both and .

This paper is organized as follows. In Section 2, we present a brief review of well-known exact solutions of the SBO system for . In Section 3, we explain the numerical methodology used to find approximations to travelling wave solutions of the SBO system. In Section 4, we compute some periodic and nonperiodic travelling wave solutions of the system, analyze the influence of model’s parameters on the geometry of periodic and nonperiodic travelling waves, and explore the collision of two right-going solitary waves of the SBO system. Finally, Section 5 contains the conclusions of our work.

2. Computing Exact Travelling Wave Solutions

For the sake of completeness, in this section, we present a brief review of a class of exact travelling wave solutions of system (1) (periodic and non-periodic) which are well known in the literature when (see [7, 10, 13] for more details and analytical results on existence of travelling wave solutions for ). These analytical special cases are important because they can be used as initial points for computing branches of solutions when , through Newton’s iteration and numerical continuation, for example.

In the nonperiodic case and , we can construct a family of exact travelling wave solutions in the form (3) of system (1). In first place, note that, for , the real-valued functions must satisfy the system

where the tildes denote differentiation with respect to the variable .

We recall that and their derivatives decay to zero at infinity. Thus we can integrate (5) to get

After substituting the value of into (4), we obtain

By multiplying the previous equation by , we get

Therefore,

Thus, by integrating (9) and using again the decaying properties of the functions , we obtain

Looking for solutions of (11) in the form , where are real constants to be determined, we have that

or, equivalently,

Using the identity , we arrive at

From this equation we can conclude that

On the other hand, evaluating (15) at ,

Substituting the value of the constant in the previous equation, we arrive at

Therefore,

On the other hand, in the periodic case and with , we also have that must also satisfy system (4)-(5). By integrating (5), we get

where is a constant. Here we assume that . Thus, after substituting this expression for into (4), we arrive at

Therefore, we obtain

Then integrating the last equation, we arrive at

where is a constant. Suppose that , are positive real numbers and is negative. Furthermore, let , (,  ) be the roots of the polynomial of the right side of the previous equations; then we can write

and thus

Therefore,

Without less of generality, we may assume that , and consider the change of variable . Substituting (24), we get

Let . Thus

Now let us define . Then

and thus

Substituting this expression into (28), we obtain

We can simplify the previous equation to get

Assuming that , then and integrating the previous equation

where . Using from the definition of the Jacobian elliptic function , we get and hence . In terms of the original variables , we obtain the dnoidal wave solutions corresponding to system (11)-(12):

where

Since the fundamental period of the dnoidal function is , where

it follows that has fundamental period given by

We point out that, due to the factor in the component of the travelling wave solution (3), we have the condition

for some .

3. Methodology for Approximation of Travelling Wave Solutions of the Full SBO System

For , we recall that if is a travelling wave solution in the form

of system (1), then are real-valued functions which must satisfy the following equations:

In first place, we are interested in finding approximations to even solutions with period , of system (1); let us introduce truncated cosine expansions for and :

where

and let us introduce analogous expressions for the coefficients . By substituting expressions (41) into (40), evaluating them at the collocation points

and using the property of the Hilbert transform

we obtain the system of nonlinear equations,

for , which can be written in the form

where the coefficients are the unknowns. We point out that the property (44) of the Hilbert transform makes the numerical solution of system (40) with the spectral collocation proposed in the present paper easier. This is the motivation for the choice of this numerical strategy.

Nonlinear system (46) can be solved by Newton’s iteration. Computation of the cosine series in (41) and the integrals in (42) is performed using the FFT (Fast Fourier Transform) algorithm. The Jacobian of the vector field is approximated by the second-order accurate formula

where and . We stop Newton’s iteration when the relative error between two successive approximations and the value of the vector field are smaller than . We use the travelling wave solutions given in Section 2, when , for initializing the Newton’s iterative scheme explained above.

4. Description and Discussion about the Numerical Results

In the following, we present some numerical experiments performed by using the numerical scheme described in the previous section.

Experiment Set 1 (Periodic Waves). In the first numerical experiment we set , , and , and the fundamental period is . The iterative Newton’s procedure is initialized with

The result of this computer simulation is presented in Figure 1. To check that we have computed really a periodic travelling wave of the SBO system (1), we run the spectral Crank-Nicolson numerical solver to approximate the solutions to this system, introduced by the author in [14], with time step , and FFT points; the spatial computational domain is the interval with , and the initial values are

with being the profiles displayed in Figure 1. The result of this computer simulation at is displayed in Figure 2, superimposed with the expected position of the travelling wave given by (39). We observe a good accordance with a maximum error of between the profiles of the modulus of components and , corroborating that, in fact, we have an approximation of a travelling wave solution to system (1). The same verification was performed successfully for other periodic travelling waves computed for different values of the parameter , displayed in Figures 3, 4, and 5 for , , and , respectively.

Figure 1: Approximations of the functions in the travelling wave solution (39) with 9 Newton’s iterations. The model’s parameters are , , and wave speed and the fundamental period is .

Figure 2: Evolution of the travelling wave solution (39) (with as in Figure 1) of the SBO system at . In solid line, computed with the numerical solver introduced in [14]. In pointed line, for the approximate travelling wave solution (39).

Figure 3: Approximations of the functions in the travelling wave solution (39) with 10 Newton’s iterations. The model’s parameters are , , and wave speed and the fundamental period is .

Figure 4: Approximations of the functions in the travelling wave solution (39) with 11 Newton’s iterations. The model’s parameters are , , and wave speed and the fundamental period is .

Figure 5: Approximations of the functions in the travelling wave solution (39) with 23 Newton’s iterations. The model’s parameters are , , and wave speed and the fundamental period is .

We remark that Angulo et al. in [13] established, by using the concentration compactness method (Lions [8, 9]), the existence of periodic travelling wave solutions of system (1), provided that , where is the fundamental period. Thus, the numerical experiment in Figure 5 with , , that is, , shows that travelling wave solutions of the SBO system may also exist outside the parameter regime considered in [13]. In Figure 6, we display another numerical simulation with , , where we also obtain an approximate travelling wave solution for , by initializing the Newton’s iteration with

In Figure 7, we show new periodic travelling waves of system (1) for computed by using as initial step for Newton’s iteration the periodic solutions obtained in (34), corresponding to the case . In this experiment we use the model’s parameters , , , and . Given the parameters , the values of the parameters are approximated numerically satisfying the relationships obtained in Section 2:

In this case, the results are , , and . The resulting fundamental period of the travelling wave is .

Figure 6: Approximations of the functions in the travelling wave solution (39) with 18 Newton’s iterations. The model’s parameters are , , and wave speed and the fundamental period is .

Figure 7: Approximations of the functions in the travelling wave solution (39) with 8 Newton’s iterations. The model’s parameters are , , and wave speed and the fundamental period is .

In Figure 8, we show the computation of other periodic travelling wave of the SBO system with , using as initial step in Newton’s scheme the periodic solution (34) with the parameters , , and . Again, the parameters are approximated numerically and the results are , , and . Thus the resulting fundamental period of the travelling wave is . We point out that, in this numerical simulation, the travelling wave is faster () than the one computed in the previous numerical experiment ().

Figure 8: Approximations of the functions in the travelling wave solution (39) with 6 Newton’s iterations. The model’s parameters are , , and wave speed and the fundamental period is .

It is important to remark that the first purpose of the previous numerical experiments is to show that the proposed Newton-pseudospectral scheme can be used to approximate with good accuracy and low computational cost, travelling-wave solutions of the SBO system. This fact was corroborated for every travelling wave computed (see, e.g., Figure 2), where we see that the corresponding wave components translate with the expected speed without changing their shape, and neither dissipation nor dispersion was observed during the time evolution. In second place, we obtained numerically new periodic travelling wave solutions of the SBO system in a parameter’s regime not covered by the theoretical results in previous works, suggesting that the regime of existence of such type of solutions of the SBO system can be extended to the case .

Experiment Set 2 (Solitary Waves). In second place, we wish to compute approximations of solitary wave solutions of the SBO system. The following numerical simulations have the goal of showing that the proposed numerical solver is not restricted to the periodic case, and it can also be further employed to obtain approximations to solitary wave solutions of the SBO system on the real line, provided that the spatial computational domain is chosen large enough so that the solution does not reach the computational boundaries.

In this case, the functions and and their derivatives must decay to zero at infinity. Thus, we run Newton’s iteration to approximate the cero of the vector field given in (46) by using as initial step the solitary wave (19) corresponding to . Here we set , , , , and FFT points. We point out that Newton’s scheme explained in the previous section is designed for searching spatial periodic solutions of the SBO system. However, due to the decay properties at infinity of a solitary wave solution, we can also use Newton’s procedure adopting a spatial computational domain as the interval with large enough. In this case we set . The result of this computer simulation is presented in Figure 9 and the corresponding verification that it corresponds really to a solitary wave solution of the SBO system is presented in Figure 10. In this validation we use the numerical parameters and FFT points for the spatial discretization in the spectral numerical scheme introduced in [14]. The maximum error observed at between the two wave profiles is roughly . Other solitary wave solutions of the SBO system were computed for different values of the parameter and the results are displayed in Figures 11, 12, and 13, for , , and , respectively. The parameters are the same as before. In the experiments in Figure 11, we use , , and for the experiments in Figures 12 and 13, we set , , and , , respectively. We recall that the wave component of a travelling solution of the SBO system is given here by . Thus, observe that, in all solitary waves computed in this paper, the component is negative. This is because solution (19) (corresponding to ) was adopted as the starting point for Newton’s iteration in the case of solitary waves.

Figure 9: Approximations of the functions in the solitary wave solution (39) with 6 Newton’s iterations. The model’s parameters are , , , and wave speed .

Figure 10: Evolution of a solitary wave solution (39) (with as in Figure 9) of the SBO system at . In solid line, computed with the numerical solver introduced in [14]. In pointed line, for the approximate solitary wave solution (39).

Figure 11: Approximations of in the solitary wave solution (39) with 6 Newton’s iterations. The model’s parameters are , , , and wave speed .

Figure 12: Approximations of the functions in the solitary wave solution (39) with 5 Newton’s iterations. The model’s parameters are , , , and wave speed .

Figure 13: Approximations of the functions in the solitary wave solution (39) with 5 Newton’s iterations. The model’s parameters are , , , and wave speed .

Experiment Set 3 (Influence of Model’s Parameters). The previous numerical experiments illustrated the shape of some travelling wave solutions of the SBO system, not computed in previous works, to the knowledge of the author. In the next set of experiments, we wish to initiate the analysis of the dependence of the these travelling wave solutions (periodic and non periodic case) on the parameters appearing in the SBO system.

In Figure 14, we compare the computed periodic travelling waves for the SBO system for . The rest of model’s parameters are fixed at , , and . From these computer simulations, we can see that the amplitudes of both components of the computed travelling wave increase when the parameter decreases. In contrast, in Figure 15, we explore the effect of the frequency on the shape of the computed travelling wave. The parameters , , and are left fixed. Notice that, in this case, the amplitude of the wave components results in being proportional to the value of the frequency . In the case of solitary waves of the SBO system, we conducted similar numerical experiments to analyze the effect of the model’s parameters on the amplitude of the computed profiles. The results are shown in Figure 16 (with ,  , , and  ) and Figure 17 (with ,  , , and ). In contrast to the case of periodic waves, in this case the behavior of the solitary wave components differs from each other. For instance, notice that the amplitude of the component decreases when increases. Meanwhile, the component decreases its amplitude when decreases. Furthermore, the amplitude of both components is directly proportional to the frequency .

Figure 14: Comparison of periodic travelling wave solutions of the SBO system for different values of the parameter .

Figure 15: Comparison of periodic travelling wave solutions of the SBO system for different values of the frequency .

Figure 16: Comparison of solitary wave solutions of the SBO system for different values of the parameter .

Figure 17: Comparison of solitary wave solutions of the SBO system for different values of the frequency .

To analyze the effect of the parameter on the shape of travelling waves of system SBO, we conduct additional numerical experiments setting , , , and . The resulting profiles computed with the Newton procedure are displayed in Figure 18. Observe that the amplitude of the computed travelling waves is inversely proportional to the parameter . On the other hand, fixing , , and , we study in Figure 19 the influence of the parameter . In these numerical experiments, we take and see that the corresponding travelling wave solution is more sensitive to changes of the parameter (which alters the first derivative term ), in contrast to previous simulations varying the other parameters in the SBO system.

Figure 18: Comparison of periodic travelling wave solutions of the SBO system for different values of the parameter .

Figure 19: Comparison of periodic travelling wave solutions of the SBO system for different values of the parameter .

Experiment Set 4 (Solitary Wave Collisions). In this set of experiments we wish to explore the interaction between two solitary waves of system (1) with different speeds of propagation and amplitudes. As mentioned above, this phenomenon has not been considered in previous works on the SBO system, to the best of the author’s knowledge. To do this, we compute the evolution of the initial profile

where are the solutions of (40) with wave speed and are the solutions corresponding to wave speed . The model’s parameters for computing the profiles ,   are taken as , , and , and the spatial domain is the interval with . We point out that the profiles , have been moved to the right a distance of 20 units, so that the initial separation between the solitary waves is 20 units. Furthermore, the profile with a higher velocity is located behind the one with the lower velocity . The numerical parameters are and FFT points in this computer simulation. The results of the interaction of the two solitary waves are displayed in Figures 20 and 21, and in the contour plots shown in Figures 22 and 23, we can see the trajectory of each solitary wave in the plane before and after the collision. When the larger solitary wave begins to shrink and the smaller solitary wave begins to grow. This process continues until the wave profiles have swapped identities, after which they separate. We always observe two local maxima in both components and during the entire interaction. When , the two solitary waves start to separate from each other. From these numerical experiments, we see that, before collision, the solitary wave with higher velocity moves towards the one with slower velocity , with the expected speeds, and then one passes through another. We point out that, after collision, some small dispersive trailing waves are generated around the pulses. Then these tails separate from the solitary waves, lengthen, and become weaker progressively as time evolves. Further observe that the tails are the strongest for the wave component , which is affected by the Hilbert operator in the second equation in the SBO system. We obtained similar results in other numerical simulations performed with different pairs of solitary wave solutions of the SBO system. This behavior is very different from that observed in other well-known integrable dispersive models, such as the Korteweg-de Vries equation [17], the Benjamin-Bona-Mahony equation [18], and the Benjamin-Ono equation [19], where the collision of two superimposed solitary-waves is elastic (i.e., clean, without forming dispersive tails). For the Benjamin-Ono equation, see also the numerical study of solitary wave collisions in the work by Thomeé and Vasudeva Murthy [20] and Pelloni and Dougalis [21]. We remark that this phenomenon has already been observed by Kalisch and Bona [22], who studied the interaction of solitary-wave solutions of the Benjamin-Ono equation, the regularized Benjamin-Ono equation, and the Benjamin equation, which are dispersive models where the Hilbert transform is also present. In accordance with that work, this behavior suggests that the SBO system is not integrable; that is, it does not have exact solitary-wave solutions, similar to the case of the Benjamin equation and the regularized Benjamin-Ono equation.

Figure 20: View in the plane of the collision of two solitary waves of the SBO system with velocities and . Observe that the colliding waves have different initial amplitude.

Figure 21: View in the plane of the collision of two solitary waves of the SBO system with velocities and . Observe that the colliding waves have different initial amplitude.

Figure 22: Wave-crest trajectory of collision of the two solitary waves of different speed and amplitude shown in Figure 20.

Figure 23: Wave-crest trajectory of collision of the two solitary waves of different speed and amplitude shown in Figure 21.

We point out that the parameter space for system (1) is large and thus, in this paper, we have presented only a limited number of numerical experiments to illustrate the geometry of travelling wave solutions of the SBO system, in a parameter regime where no exact solution is known, and to initiate the study of interactions between these permanent wave states. Further numerical and analytical studies are necessary to determine the complete parameter space for existence of travelling waves to study important analytical properties, such as orbital stability, and to develop an exhaustive discussion of solitary wave interactions for the SBO system.

5. Conclusions

In this paper, we introduced a numerical scheme for approximating the travelling wave solutions (periodic and nonperiodic) of system (1), whose existence was established in previous works [7, 10, 13]. This scheme combines a Fourier spectral discretization together with a Newton-type iteration, using as initial step well-known exact solutions when . Using the spectral Crank-Nicholson numerical scheme introduced by the author in [14] for approximating time evolution of solutions of the SBO system, we also illustrated the collision of two right-running unequal-amplitude solitary waves propagating with different speeds. The behavior of the colliding solitary waves was found to be similar to solitary waves of other well-known nonintegrable dispersive-type equations, such as the scalar regularized Benjamin-Ono equation and the Benjamin equation, for instance. As mentioned in the work by Kalisch and Bona [22], this phenomenon suggests that the SBO system is not integrable. Furthermore, we obtained numerically new periodic travelling wave solutions which are not included in the analytical theory presented in the work by Angulo et al. [13].

Conflicts of Interest

The author declares that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This research was supported by Universidad del Valle, Calle 13, Nro. 100-00, Cali, Colombia, under Research Project C.I. 71020.

References

M. Funakoshi and M. Oikawa, “The resonant interaction between a long internal gravity wave and a surface gravity wave packet,” Journal of the Physical Society of Japan, vol. 52, no. 6, pp. 1982–1995, 1983.

V. I. Karpman, “On the dynamics of sonic-langmuir solitons,” Physica Scripta, vol. 11, no. 5, pp. 263–265, 1975.

V. D. Djordjevic and L. G. Redekopp, “On two-dimensional packets of capillary-gravity waves,” Journal of Fluid Mechanics, vol. 79, no. 4, pp. 703–714, 1977.

R. H. J. Grimshaw, “The Modulation of an Internal Gravity-Wave Packet, and the Resonance with the Mean Motion,” Studies in Applied Mathematics, vol. 56, no. 3, pp. 241–266, 1977.

D. J. Benney, “Significant interactions between small and large scale surface waves,” Studies in Applied Mathematics, vol. 55, no. 2, pp. 93–106, 1976.

D. J. Benney, “General theory for interactions between short and long waves,” Studies in Applied Mathematics, vol. 56, no. 1, pp. 81–94, 1977.

J. Angulo and J. F. Montenegro, “Existence and evenness of solitary-wave solutions for an equation of short and long dispersive waves,” Nonlinearity, vol. 13, no. 5, pp. 1595–1611, 2000.

P.-L. Lions, “The concentration-compactness principle in the calculus of variations. The locally compact case. Part 1,” Annales De L Institut Henri Poincaré-Analyse Non Linéare, vol. 1, no. 2, pp. 109–145, 1984.

P.-L. Lions, “The concentration-compactness principle in the calculus of variations. The locally compact case. II,” Annales de l'Institut Henri Poincaré (C) Analyse Non Linéaire, vol. 1, no. 4, pp. 223–283, 1984.

J. Angulo Pava, “Stability of solitary wave solutions for equations of short and long dispersive waves,” Electronic Journal of Differential Equations, No. 72, 18 pages, 2006.

D. Bekiranov, T. Ogawa, and G. Ponce, “Interaction equations for short and long dispersive waves,” Journal of Functional Analysis, vol. 158, no. 2, pp. 357–388, 1998.

H. Pecher, “Rough solutions of a Schrödinger-Benjamin-Ono system,” Differential and Integral Equations. An International Journal for Theory & Applications, vol. 19, no. 5, pp. 517–535, 2006.

J. Angulo, C. Matheus, and D. Pilod, “Global well-posedness and non-linear stability of periodic traveling waves for a Schrödinger-Benjamin-Ono system,” Communications on Pure and Applied Analysis, vol. 8, no. 3, pp. 815–844, 2009.

J. C. Muñoz Grajales, “Error analysis of an implicit spectral scheme applied to the Schrödinger-Benjamin-Ono system,” International Journal of Differential Equations, Article ID 6930758, Art. ID 6930758, 12 pages, 2016.

T. B. Benjamin, “Solitary and periodic waves of a new kind,” Philosophical Transactions of the Royal Society of London. Series A. Mathematical, Physical Sciences and Engineering, vol. 354, no. 1713, pp. 1775–1806, 1996.

J. Angulo, M. Scialom, and C. Banquet, “The regularized Benjamin-Ono and {BBM} equations: well-posedness and nonlinear stability,” Journal of Differential Equations, vol. 250, no. 11, pp. 4011–4036, 2011.

T. R. Marchant and N. F. Smyth, “Soliton interaction for the extended Korteweg-de Vries equation,” IMA Journal of Applied Mathematics, vol. 56, no. 2, pp. 157–176, 1996.

T. B. Benjamin, J. L. Bona, and J. J. Mahony, “Model equations for long waves in nonlinear dispersive systems,” Philosophical Transactions of the Royal Society A: Mathematical, Physical & Engineering Sciences, vol. 272, no. 1220, pp. 47–78, 1972.

W.-R. Sun, B. Tian, H. Zhong, and R.-X. Liu, “Breather and double-pole solutions for the Benjamin-Ono equation in a stratified fluid,” Waves in random and complex media : propagation, scattering and imaging., vol. 26, no. 2, pp. 168–175, 2016.

V. Thomeé and A. S. Vasudeva Murthy, “A numerical method for the Benjamin-Ono equation,” BIT Numerical Mathematics, vol. 38, no. 3, pp. 597–611, 1998.

B. Pelloni and V. A. Dougalis, “Numerical solution of some nonlocal, nonlinear dispersive wave equations,” Journal of Nonlinear Science, vol. 10, no. 1, pp. 1–22, 2000.

H. Kalisch and J. L. Bona, “Models for internal waves in deep water,” Discrete and Continuous Dynamical Systems - Series A, vol. 6, no. 1, pp. 1–20, 2000.