Explosion of Soliton in A Magnetic Field: 128.114.34.22 On: Sun, 30 Nov 2014 00:16:13
Explosion of Soliton in A Magnetic Field: 128.114.34.22 On: Sun, 30 Nov 2014 00:16:13
Explosion of Soliton in A Magnetic Field: 128.114.34.22 On: Sun, 30 Nov 2014 00:16:13
Dispersive properties of a magnetized plasma with a fieldaligned drift and inhomogeneous transverse flow
Phys. Plasmas 3, 3091 (1996); 10.1063/1.871656
This article is copyrighted as indicated in the article. Reuse of AIP content is subject to the terms at: http://scitation.aip.org/termsconditions. Downloaded to IP:
128.114.34.22 On: Sun, 30 Nov 2014 00:16:13
Explosion of soliton in a magnetic field
Katsuhiro Nishinari
Department of Aeronautics. Faculty of Engineering. The University of Tokyo. Hongo 7-3-1. Bunkyo-ku.
Tokyo 113. Japan
Kanji Abe
College of Arts and Sciences. The University of Tokyo. Komaba 3-8-1. Meguroku. Tokyo 153. Japan
Junkichi Satsuma
Department of Mathematical Sciences. The University of Tokyo. Komaba 3-8-1, Meguroku, Tokyo 153.
Japan
(Received 12 July 1994; accepted 9 September 1994)
A dynamics of a solitary pulse of the electrostatic ion cyclotron wave that propagates perpendicular
to an applied magnetic field is considered. It is shown that the solitary wave will be singular in some
range of parameters in the system, such as the plasma density and the magnitude of an applied
magnetic field. This fact shows that there is a possibility of controlling the place where explosion
of the solitary wave occurs. © 1994 American Institute of Physics.
3728 Phys. Plasmas 1 (12), December 1994 1070-664X/94/1 (12)/3728/3/$6.00 © 1994 American Institute of Physics
This article is copyrighted as indicated in the article. Reuse of AIP content is subject to the terms at: http://scitation.aip.org/termsconditions. Downloaded to IP:
128.114.34.22 On: Sun, 30 Nov 2014 00:16:13
and similarly for <P and v. We calculate, up to the third order (a3+ s / a 4)
of E, and finally obtains
f3= albi/(2alP)2+a2b~/(2a2q)2' (16)
iA r + aIAU+ a2A7]7]+ a31AI2A + a4AQ=O, (6) The equation (15) is generally solved by the elliptic function.
2 For the homoclinic orbit of this system, the solution is ex-
Q~e+ 01 Q 7]7]+ 02(IAI )u=o, (7)
pressed by
where A = n\l) and Q = nb 2
) + 03V~~). The coefficients aj,
i = I, 2, 4 and 0;, i = I, 2 are given by 8ad exp( rat)
R= f3+8ad 2 exp(2Ta?)' (17)
I a2wl2 I a2wl
0'1 =2 ak ,0'2=2 ak2 , where d is an arbitrary constant. If a>O and f3>O, (17) gives
x kx=O Y kx=O a regular solitary wave. The solution (17) will be singular if
f3<O. Existence of this type of solitary wave solution indi-
cates that there is a possibility to realize an explosion. That
is, if a is kept positive and f3 becomes negative, then the
See Ref. 5 for coefficients 0'4 and b':l, which are complicated solitary wave will explode. This is easily imaged by the fact
functions of the wave number of the carrier wave, the plasma that the maximum amplitude of R, J2 al f3, goes to infinity
density, and the magnitude of magnetic field. The system (6) as f3'-,,0.
and (7) is the well-known DS equations. It must be noted here that a and f3 given in (16) are
We now consider the solitary wave solution of the DS functions of the plasma density, and the magnitude of the
equations. Let us assume that the solution is of the form magnetic field, both of which can be controlled from outside
the system. Using this fact, a scenario of the explosion is
A =R(x,y,t)exp i(px+qy+ rt), (9) described as follows: First, consider a situation where a soli-
tary wave is excited and propagates in the region where a>O
where R(x,y,t) is a real function and p, q, and r are real
and f3>O. We expect to have a stable traveling wave in this
constants. We substitute (9) into (6) to obtain
region. Assume that the system is inhomogeneous and the
R t + 2alpRx+2a2qRy=O, (10) magnitude of the plasma density or the magnetic field varies
locally. Then, due to a change in the value of a, we have a
from the imaginary part. A solution of (10) is
change that f3 becomes negative while a is kept positive. The
b1x b 2y solitary wave explodes at that place. Of course, we can arti-
R =R(~), ~= -2-+ -2--(b l +b 2)t. (11) ficially give such an inhomogeneity to the system, and con-
alP a2q
trol the place where explosion occurs. We think that this idea
In (11), b l and b 2 are constants that determine the direction will be useful for an energy transport in a plasma system.
of the wave propagation. Substituting (11) into (7), we find Perturbation analysis may not be applicable to the singu-
that lar case. Moreover, the existence of inhomogeneity may
cause other phenomena, such as deformation and fission of
(12) soliton.9 We consider, however, that the tendency for explo-
sion may appear in a real system. The reason of this is as in
where the prime means the derivative with respect to ~. In- the following. Let us consider "collapse," which was first
tegration of (12) gives discussed by Zakharov. 10 The explosion noted above may
have some resemblance to the collapse. The collapse is ex-
Q=sIR2+CO~+CI , (13) plained by the singularity of a localized solution in the mul-
tidimensional NLS equation, which is derived by the same
where
perturbation analysis as ours. Since the tendency for collapse
is observed in real systems such as plasma, optics, etc.,l1 we
may say that the explosion will also be expected to occur.
The physical mechanism of explosion, however, is dif-
and co' C I are integration constants. The constant Co can be ferent from that of collapse. In the case of collapse, the am-
taken to be zero from a physical point of view. plitude goes to infinity within a finite time, keeping its total
The real part of (6) gives energy. It is apparent, however, that in our explosion the total
energy goes to infinity. This is easily seen by considering the
- rR+ 0'1 (R xx - p2R) + a2(R yy -q 2R)+ a3R3+ a4RQ=0.
(14) area of cross section of the soliton given by f R d ~
= 'TrJ21 f3, which goes to infinity as f3'-,,0. The soliton is
We substitute (II) and (13) into (14) to obtain sustained by the balance between the nonlinear effect of the
R"-aR+f3R 3 =O, (15) ponderomotive force and the dispersion effect. From our re-
sults, we can artificially alter the balance and make the soli-
where ton taller by varying the parameter a. As its cross section of
it represents the wave energy, we can then conclude that the
explosion is caused by injecting energy into the soliton from
outside through a suitable change of the parameter a.
-0.15
-0.2
FIG. I. Dependence of the sign of ex on a and k, . The shaded parts and {3 becomes negative with a>O, which are shown by three
unshaded parts indicate the region where ex is negative and positive, respec- arrows in Fig. 2. Let us analyze one of these in detail. Figure
tively. On the thick curve, ex goes to infinity.
3 shows the a dependence of {3 for k y = 0.8. From Fig. 1 and
Fig. 3, we find a regular solitary wave exists in the region
Finally, we give an estimation of whether the above situ- 0.359<a. If a becomes smaller than 0.359, {3 becomes
ation really occurs or not in this system. This is done by negative with a>O.
studying the signs of a and {3 using (8). Since there are many In conclusion, we have proposed in this Letter a new
parameters in a and {3, we only give some examples. In the type of explosion, which is induced by changing physical
following calculation, all parameters except a and ky are parameters in the system. We think that this type of explo-
fixed as p=2, q= t, m =n = I, r= -2, and Cl =0. The signs sion occurs also in other fields, such as fluid dynamics, op-
of a and {3 are illustrated in Fig. 1 and Fig. 2, respectively. tics, and so on. Investigation of the possibility will be a
From these figures, we see that there exist three cases where subject of future work.
ACKNOWLEDGMENTS
We are grateful to Dr. P. M. Santini, Dr. M. Tanaka, Dr.
T. Yajima, and Mr. K. Matsui for helpful suggestions and
ky comments. Thanks are due to Ms. Yumiko Nojima for her
1.2
assistance in writing figures.
This work is partially supported by a Grant-in-Aid for
the Encouragement of Young Scientists from the Ministry of
Education, Science. and Culture (No. 2426).
1A. Davey and K. Stewartson, Proe. R. Soc. London Ser. A 338. 101
(1974).
2A. S. Fokas and M. J. Ablowitz. J. Math. Phys. 25. 2494 (I 984}.
3J. Satsuma and M. J. Ablowitz, J. Math. Phys. 20, 1496 (1979).
4M. BoW. J. J.-P. Leon. L. Martina. and F. Pempinelli, Phys. Lett. A 132,
432 (1988).
5K. Nishinari, K. Abe, and J. Satsuma, Phys. Plasmas 1, 2559 (1994).
o 0.2 0.4 0.6 0.8 1.2 1.4 6K. Shimizu and Y. H. Ichikawa, J. Phys. Soc. Jpn. 33, 789 (1972).
a 7K. Nishinari. K. Abe, and J. Satsuma, J. Phys. Soc. Jpn. 62, 2021 (1993).
8B. Malomed, K. Matsui. and J. Satsuma, "Solitons and collapse in the
generalized deri vati ve nonlinear Schriidinger equation," submitted to J.
FIG. 2. Dependence of the sign of f3 on a and k y . The shaded parts and Phys. Soc. Jpn.
unshaded parts indicate the region where f3 is negative and positive, respec- 9y. S. Kivshar and B. A. Malomed, Rev. Mod. Phys. 61. 763 (1989).
tively. On thick curves, f3 goes to infinity. Three arrows show ways of 1Ov. E. Zakharov. Zh. Eksp. Teor. Fiz. 62. 1745 (1972).
changing the sign of {3. which cause the explosion. 11M. V. Goldman. Rev. Mod. Phys. 56, 709 (1984).