San José State University |
---|

applet-magic.comThayer WatkinsSilicon Valley & Tornado Alley USA |
---|

The Creation of Secondary Solitary Waves
From the Collision of Solitary Wave Solutions to the Regularized LongWave Equation |

One major fascination of solitons is their character as waves and particles. Solitons collide, interact, and emerge unchanged except for a phase shift. The analogy to subatomic particles is compelling. But more interesting than solitons is the behavior of the solitary solutions to the Regularized Long-Wave (RLW) equation in which solitary solutions collide and may create secondary solitary waves or sinusoidal solutions. This corresponds even better to sub atomic phenomena in that collisions of particles create other particles and/or radiation. The purpose of this paper is to study the process for the creation of the secondary waves. This provides the opportunity to investigate the creation process for possible insights into the corresponding processes of particle physics.

The RLW equation has the form

whereas the KdV equation may be expressed as

The RLW equation was formulated by Peregrine (1966) as an alternative to KdV equation for studying soliton phenomenon. It was proposed because it would not have the same limitations for the size of the time step in numerical solution that the KdV has.

Some aspects of the solutions to the KdV and RLW equations may be derived from analysis. A
traveling wave solution is of the form u(x-vt-x_{0}). Letting z=x-vt-x_{0}
the RLWE becomes:

Integrating this once with respect to z gives

where c_{0} is an arbritrary constant.

Multiplying the above equation by u_{z} gives and equation that can be integrated
with respect to z to give

Finally the equation can be solved for u_{z} and integrated.
It
is readily found that the RLWE
has solitary wave solutions of the form a·sech^{2}(b(x-vt-x_{0}).
That is to say, if at t=0,
u=a·sech^{2}(bx) then the solution for all future times will be of that form.
This is a form found frequently in
soliton solutions. When v is positive the waveform moves to the right and when it is negative the wave form moves to
the left.

Collisions of solitary solutions are created by taking initial conditions which are
the sum of two or more
solutions. The collision occurs at the value of x and t such that the argument
(x-vt-x_{0}) is zero for both waves.
It is shown elsewhere that the parameters a, b, and v must satisfy the following
conditions:

and

a=12b

These conditions are equivalent to

and

v = l+a/3.

It is easily seen that for b to be real a must be greater than zero or less than -3. This forbidden interval of (-3,0) for a is of importance. It means that v will be greater than 1 or less than 0. Note that the larger the magnitude of the amplitude the closer the value of the parameter b is to (1/2) and consequently the larger the speed of the wave. The parameter b indicates how sharp of a peak the wave has. The width of the wave, as measured by a e-ratio or standard deviation, is inversely proportional to b. Thus the larger the magnitude of the amplitude the narrower the wave. The minimum width is 2.

Initially it was thought that because of the similarity of the RLW equation to the KdV equation that it would also have n-soliton solutions. Eilbeck and McGuire found in 1975 that after the collision of two solitary waves the same waveforms reappeared with amplitudes within 3/10 of 1 percent of the originals. It was concluded that the RLW equation had soliton solutions, but Abdulloev (1976) found that the discrepancy could not be accounted for by numerical inaccuracies.

Santarelli in 1978 found that secondary solitary waves are created by the collision of large amplitude solitary waves. The results of Santarelli's study are shown in Figures 1 through 9.

In Figure 1 the positive amplitude wave is to the left of the negative amplitude wave.

In Figure 2 the waves are colliding.

In Figure 3 the waves are superimposed but higher magnitude of the negative amplitude wave results in it predominating. There is a slight change of scale for Figure 3 to allow the display of more detail. It must be emphasize that these results are not simply from the addition of the two initial waves. The two initial waves are used as a initial conditions for the numerical solution of the RWL equation.

In Figure 4 the waves have effectively passed through each other, but with diminished amplitude.

In Figure 5 the secondary waves can be seen developing. There is a notch on the negative amplitude wave that will develop into the secondary negative amplitude wave. The secondary positive amplitude is also developing but it is less obvious at this stage.

In Figures 6 through 9 the secondary waves have developed and continue to move; the positive amplitude to the right and the negative amplitude waves to the left.

Lewis and Tjon (1979) using even
larger amplitude solitary waves found that multiple secondary waves are created, always
in pairs with positive and
negative amplitudes. Furthermore, Lewis and Tjon found resonance effects in which the
relationships between the
parameters of the solitary waves became important.
They also
found cases in which the the positive and negative solitary waves annihilate each other
and create a sinusoidal
residual solutions called appropriately *radiation*.

According to Lewis and Tjon the known conserved quantities of the RLW equation are:

where all of the integrations are over the interval [−∞,+∞].
For ^{2}(bz)_{0}

so the first invariant is equivalent to
(a_{i}/b_{i}) for the waves before
the collision being equal to
(a_{j}/b_{j} after the collision. The second invariant for a solitary
wave ^{2}(bz)

which in turn is equal to

Since
∫sech^{6}ydy/∫sech^{4}ydy = 4/5 the second invariant is
proportional to a^{2}[(1/b)+(16/5)b] Because the integrand of
the second invariant is nonlinear it cannot be asserted that its
value for u_{1}+u_{2} is the sum of its values for u_{1} and
u_{2}.

However if the waves for u_{1} and u_{2} are widely separated so
that u_{1}u_{2} is virtually zero then additivity is closely
approximated. Thus a^{2}[(1/b_{i})+(16/5)b_{i}] long before the
collision is equal to its value long after the collision.

The third invariant for a solitary wave is proportional to
a^{2}[(4a/15b) + (4/5)b] and so the sum of these quantities before the
collision is equal to the sum for the
waves long after the collision.

Table 1 gives the parameters of the incident waves in the collision studies and the secondary waves. Unfortunately the studies do not give the parameters of the primary waves after the collision so one cannot use the results to illustrate the invariances. The changes in the primary waves were small so they were neglected.

Table 1: Paramters of the Solitary Wave Solutions of the RLW Equation Involved in Studies of Collisions |
|||
---|---|---|---|

Amplitude | Shape | Velocity | |

a | b | v | |

Primary | |||

5.3333 | 0.4000 | 2.7777 | |

-9.8181 | 0.6000 | -2.2727 | |

Secondary | |||

1.074 | 0.2567 | 1.3579 | |

-3.825 | 1.0770 | -0.2748 | |

Primary | |||

35.2653 | 0.4800 | 12.7551 | |

-39.7647 | 0.6000 | -12.2549 | |

Secondary | |||

17.6326 | 0.4622 | 6.8777 | |

-23.8588 | 0.5347 | -6.9630 | |

8.8163 | 0.4319 | 3.9388 | |

-14.3153 | 0.5624 | -3.7718 |

It is possible to make some interesting computations from the results of the studies. In both studies the magnitude of a/b for the left-moving secondary wave is 46 percent of the combined magnitude of a/b for both secondary waves. The value of the second invariant for the negative amplitude wave is 64.01 whereas its value for the positive amplitude secondary wave is 5.44. Thus there is no equipartition principle involved in the creation of the secondary solitary waves. The values of the third invariant for the negative and positive amplitude secondary waves are -1.25 and 1.53, respectively.

A collision of two primary waves in which two secondary waves are produced involves four unknown parameters, the amplitudes of the primary waves after collision and the amplitudes of the secondary waves. There are three invariances but they are nonlinear so it is not possible to say whether the unknown parameters are determined uniquely by the invariances; i.e., the conservation laws. Actually the after-collision waves cannot be a unique solution of the invariances because the there is always the solution of the after-collision waves being equal to the before-collision waves.

Let u_{1} and u_{2} be two solutions to the RLW equation and
let u_{3} be a solution with the
initial conditions equal to u_{1}+u_{2}
at t=0.

The RLW equation may then be put into the form

Now let w=u_{1}−u_{2}.
Since u_{1} and u_{1},
satisfy the RLW equation it
follows that w satisfies the equation

(w-w

This is the equation which describes the formation of secondary solitary waves. The initial
condition is w(x)=0 for all x. A collision in where (u_{1}u_{2}) is nonzero and
such that (u_{1}u_{2})_{x} is nonzero. This quantity, the gradient
of (u_{1}u_{2}), is what drives the creation of the secondary waves.
From the numerical studies it is
reasonable to conclude that
w will consist of two separated parts; one associated with the positive amplitude
solitary wave and one with the
negative amplitude wave.

A rearrangement puts the equation into the form

(w

A form for numerical solution is

(Δw

where Δw is w(t+Δt)-w(t).

This is an equation of the form

which has the solution

Therefore the numerical solution for the secondary wave amplitude w is given by

w(x,t+Δt) = w(x,t) + Δt∫

where f(x) is given by

f(x) = (w+½w

or, in a more computationally efficient form

f(x) = [w((1+½w)+u

Since the forcing term f(x) is a derivative it is best to make use of integration by parts
to eliminate the need to approximate the derivative. Thus let f(x)=F_{x}(x) then

The limit of sinh(x-z)F(z) as z→−∞ needs to be investigated but for now it is presumed to be zero. Since sinh(0)=0, the estimating equation for computation is

w(x,t+Δt) = w(x,t) − Δt∫

where F(x) is given by

F(x) = (w+½w

or in the computationally more efficient form

F(x) = w(1 + ½w + u

Suppose that u_{1} and u_{1} are nearly Dirac delta functions of opposite
signs that coincide (except for sign) at x=0 at t=0. Initially (t=0) the F function will be zero everywhere
except at x=0 where u_{1} and u_{2} coincide. At the first time step
Δt, then
w(x,Δt) will be -Δt*cosh(x)∫F(z)dz for x>0,

where ∫F(z)dz denotes the integral across the interval about zero. For x<0 the convolution integral will be zero because z≤x and thus the integration never covers the interval for which F(z) is nonzero.

The w function will be zero for x<0 and some positive
value, say w_{0} for x>0. For t≥0 the product u_{1}u_{2} will be
zero everywhere.

Before dealing with the interaction of two solitary waves it is worthwhile to show how such waves procede without interaction, as is shown below. In this diagram the images of past positions are not erased.

So the initial time step Δt produces values for the w function of the form

= -Δt[cosh(x)A] for x≥0

where A=∫F(z)dz.

The value of Δt should be on the order of the time it takes for u_{1}
and _{2} to slip past each other. This would be sech² waves on the order
of (1/b_{1}+1/b_{2})/(|v_{1}|+|v_{2}|). For the next time
step u_{1}u_{2} is zero and the w function is determined by the term w*(u_{1}+u_{2}),
so except for the intervals covered by u_{1} and u_{2} there will not be
any augmentation of the level of w.

Once the level of w over some interval is such that the w+½w² term is
significant then the level of w in that interval grows on its own without any help from
the u_{1} or u_{2} functions.

(To be continued.)