数理解析研究所講究録 1271 巻 2002 年 125-134

125
Modffied KP 方程式のソリトン解の二次元的相互作用(Hidekazu TSUJI)

九大・応力研辻英一

九大・応力研及川正行 (Masayuki OIKAWA)

Research Institute for Applied Mechanics, Kyushu Univ.

1

はじめに水の波・プラズマなどの基礎方程式がら導かれるソリトン方程式につぃてこれまで多くの研究がなされてきたが, それらは空間 1 次元を仮定して導かれた物が多く , 表現できるソリトンの相互作用は 1 次元 (1 方向) の場合に限られる. 実際の空間 2 次元系で, 峰が直線的に伸ひてぃる
1 次元的ソリトンが角度を持って相互作用する, すなわち進行方向が互いに異なる 1 次元的なソリトンが同時に存在する場合, 支配方程式は当然 2 次元性を持っものでなければいけないが, 相互作用がきちんと解として書ける事は少ない.

進行方向の角度がある程度異なる場合には, 摂動論を使って相互作用が調べられる場合がある (例えば,一層が深い二層流体についての研究 $[1, 2]$ など). しかしながら, 孤立波がほとんど同じ方向に伝ゎる場合には摂動論が使えなくなる. この場合, 主に伝わる方向に対して直交する方向への変化が低いオーダーになる (弱 2 次元的)と仮定をした時基礎方程式系がら近似的に 2 次元的な方程式が導出される場合がある.

そのような例として我々はこれまでに, 二層流体中の波動を表す支配方程式である “Extended Kadomtsev-Petviashvili(EKP)” Equation [3]$\frac{\partial}{\partial x}[\frac{\partial u}{\partial t}+6Pu\frac{\partial u}{\partial x}-6Qu^{2}\frac{\partial u}{\partial x}+\frac{\partial^{3}u}{\partial x^{3}}]+\frac{1}{2}\frac{\partial^{2}u}{\partial y^{2}}=0$$\mathrm{f}$” $2$-dimensional Benjammin-Ono$(2\mathrm{d}\mathrm{B}\mathrm{O})$$(P>0, Q>0)$ ,”Equation [4]$\frac{\partial}{\partial x}(\frac{\partial u}{\partial t}+u\frac{\partial u}{\partial x}+\mathcal{H}[\frac{\partial^{2}u}{\partial x^{2}}])+\frac{\partial^{2}u}{\partial y^{2}}=0$,$\mathcal{H}[f(x)]\equiv\frac{1}{\pi}\mathrm{P}\int_{-\infty}^{\infty},\frac{f(x’)}{x-x}dx’$,などを数値的に調べてきた.

今回我々は負イオンを含むプラズマのある種の振る舞いを表す (弱 2 次元的) モデル方程式である “ModifimlKP (MKP)” Equation,(1)$\frac{\partial}{\partial X}[\frac{\partial u}{\partial T}+6u^{2}\frac{\partial u}{\partial X}+\frac{\partial^{3}u}{\partial X^{3}}]+\frac{\partial^{2}u}{\partial \mathrm{Y}^{2}}=0$を取り上げる 1.

この方程式は,1 次元性 ($\mathrm{Y}$方向の変化を無視) を仮定すると, ソリトン方程式として知られてぃる

Mo 市 fied Korteweg-de Vries(MKdV) 方程式になり, また, 非線形項が 2 次

分な Kadomtsev-Petviashvili(KP) 方程式 [5] となる.

$\mathrm{K}\mathrm{P}$$(u\partial u/\partial X)$

ならばやはり可積

方程式ではソリトン解が解析的に求まり, 相互作用する角度によって相互作用の性質が変わるというソリトン共鳴現象が起こることがゎかってぃる.

$\mathrm{K}\mathrm{P}$

方程式と異なるタイプの非線形項を持っ MKP 方程式において上のような角度依存性があるかどぅかに注目しながら,1 次元 MKdV ソリトンの 2 次元的相互作用を調べる. また, 今までに調べた他の 2 次元方程式の結果との比較も試みる. なお, しばしば Konopelchenko にょって提出された可積分方程式を MKP 方程式と呼ぶ場合があるが, ここでは式 (1) を指す事にする.

13 次の非線形項の符号は, EKP 方程式と}g 違 V , 正である.
$\backslash$

126
2定式化

MKP 方程式は, プラズマの基礎方程式から近似的に導出されるが, 簡単にそれにつ て示す. 空間 1次元の場合の MKdV 方程式は Watanabe [6] によって導出された. 2 次元方程式はその拡張として導出さ$\mathrm{A}\mathrm{a}$れる.

基礎方程式は, 正と負のイオンについての運動方程式と連続の式, 電子のボルツマン方程式と静電ポテンシャルについての式である.
$\frac{\partial n_{\epsilon}}{\partial t}+\nabla\cdot(n_{\epsilon}u_{\epsilon})=0$,$\partial^{2}\phi$$n_{\mathrm{e}}=\exp\phi$

ここで,$n$は密度,$\mathrm{u}=(u, v)$

イオン, 電子を表し,は速度,,$\phi$である.

$s=\alpha \mathrm{o}\mathrm{r}\beta$$(C_{\alpha}=-1,$$\frac{\partial \mathrm{u}_{\epsilon}}{\partial t}+(u_{\epsilon}\cdot\nabla)\mathrm{u}_{\epsilon}=C_{\epsilon}\nabla\phi$.

$\partial^{2}\phi$$\frac{-\tau}{\partial x^{2}}+’=n_{\mathrm{e}}+n\rho-n_{\alpha}\overline{\partial y^{2}}$は静電ボテンシャルである. 添字$Q$は質量比 m\beta /m。であり,量は, 摂動を受けていない時の電子密度, イオンプラズマ振動数やデ’ くる.

ここで変動を表す微小量$\epsilon$.

$\alpha,\beta,e$はそれぞれ正イオン, 負$\nabla=(\partial/\partial x,\partial/\partial y)$$\text{イ}$,である. 全ての長などによって無次元化されて$\mathrm{A}\mathrm{a}$によって展開する.

$n\text{。}=1+\epsilon n_{\mathrm{e}1}+\epsilon^{2}n_{\mathrm{e}2}+\epsilon^{3}n_{\mathrm{e}3}+\cdots$,$u_{\epsilon}=\epsilon u_{\alpha 1}+\epsilon^{2}u_{\alpha 2}+\epsilon^{3}u_{\alpha 3}+\cdots$$n_{\epsilon}=n_{\epsilon 0}+\epsilon n_{\epsilon 1}+\epsilon^{2}n_{\epsilon 2}+\epsilon^{3}n_{\epsilon 3}+\cdots$,$v_{\epsilon}=\epsilon^{2}v_{\epsilon 1}+\epsilon^{3}v_{\epsilon 2}+\cdots$$\phi=\epsilon\phi_{1}+\epsilon^{2}\phi_{2}+\epsilon^{3}\phi_{3}+\cdots$

また,$C \rho=\frac{1}{Q})$

,

,

.

座標は以 T のような変換を行う.
$\xi=\epsilon(x-Vt)$

弱 2 次元性の仮定が程式は$\mathrm{K}\mathrm{P}$,$y$$\eta=\epsilon^{2}y$,$\tau=\epsilon^{3}Vt$. だし $V$ は$(^{}$$V^{2}+ \frac{3n\rho 0}{Q^{2}V^{2}}-\frac{3n_{\alpha 0}}{V^{2}}=0$

方向のスケーリングに入っている. なお,$V$を満たす).

についての条件が無$\mathrm{A}\mathrm{a}$場合はモデノレ方方程式となる. 計算すると,,$\frac{\partial}{\partial\xi}[\frac{\partial\phi_{1}}{\partial\tau}+\frac{1}{4}(\frac{15n\rho 0}{Q^{3}V^{6}}-\frac{15n_{\alpha 0}}{V^{6}}-1)\phi_{1}^{2}\frac{\partial\phi_{1}}{\partial\xi}+\frac{1}{2}\frac{\partial^{3}\phi_{1}}{\partial\xi^{3}}]+\frac{\partial^{2}\phi_{1}}{M}=0$と MKP 方程式が導かれる.

3

適当な変数変換をする事によって以下式 (1) の形を取り扱う.

解析についてMKP 方程式の 2 次元的な解は解析的に求められていない. 例えば, 双線形形式による解法を考えることができるが, これについては Chen らが EKP 方程式について行った解析 [7] がほとんど同じ今のところ,形で適用できる. それによると 1 ソリトン解 (本質的に MKdV ソリトンである) は簡単に求められるが,2 ソリトン解を求めようとするとその 2 つのソリトンは同じ方向を向くという条件を満たす必要があるた
め, 結局 MKdV 方程式の 2 ソリトン解が出ることになり, 求めたい 2 次元的な解は得られな .
$\mathrm{A}\mathrm{a}$このため,ここでは数値的に調べる事を考える. 具体的には, 方向の異なる 2 つの MKdV ソリトンを置いた初期条件を与え, どのように波が時間発展していくかを調べる.

127
数値計算の方法は以前に$2\mathrm{d}\mathrm{B}\mathrm{O}$

方程式を調べた [4] 方法とほぼ同じである. 空間微分は擬スペクトル法を用い, 時間発展は, 非線形の陰的部分を反復で解$\text{く}$Crank-Nicolson 法を用いる. 初期条件は以下のように置く,$u=A\mathrm{s}\mathrm{e}\mathrm{c}\mathrm{h}$$u=A\mathrm{s}\mathrm{e}\mathrm{c}\mathrm{h}$ここで$L_{\mathrm{Y}}$は$\mathrm{Y}$$[A((X-X0+\Omega \mathrm{Y})]$for $0<\mathrm{Y}<L_{Y}/2$ ,for$[A((X-X\mathit{0}-\Omega \mathrm{Y})]$$L_{\mathrm{Y}}/2<\mathrm{Y}<L_{\mathrm{Y}}$.

方向の領域の大きさである. パラメーター $A$ が孤立波の振幅を 2,分をそれぞれ表す.

$\Omega$$\Omega$が進行方向のが大きいほど, 進行方向の違いが大きい事になる. 時間が進むにっれて,$\mathrm{Y}$成

2 っの孤立波は $X$ について正の方向に進むので, 相互作用する場所も 2 っの孤立波の対称軸の所で正の方向へと進んでいく.

境界条件は $X,$$\mathrm{Y}$方向ともに周期境界条件とするが,$\mathrm{Y}$方向については, 孤立波が外から入ってきてぃるとみなせるように領域の端の値をある時間単位で設定する, 計算で用いたパラメーターは $A=1,0.75,0.5$ ,$\Omega$については 0.25 から 4 の間の値である.

である. 領域の点数は $X$ 方向に 4096 点,条件の対称性により,4$\mathrm{Y}$

時間空間のメッシュ幅は $dX=0.1$ ,$\mathrm{Y}$$d\mathrm{Y}=0.1$,$dT=0.005$

方向に 1024 点を取る (実際は領域の周期性と方程式・初期方向は半分の点数で計算できる).

計算結果
まず, 系の時間発展について 2 つの典型的な例を示す. この例は振幅 $A$ を一定にして$\Omega$を変化させた場合である.

41

進行方向の成す角度が大きい (進行方向の成す角度が大きい ($\Omega$$\Omega$が大きい) 場合の結果が大きい) 場合の結果を Fig.l に示す.

Fig .l(a) が初期状態であり, 時間発展するにつれて領域の中心で 2 っのソリトンが出会って相互作用を行う. この場合, 相互作用が進むにつれて, 峰の伸ひる方向が元の波とほぼ等しい 2 つの孤立波が形成され伝播していく. また領域中央付近の相互作用する場所には, 相互作用前のソリトンの 2 倍をわすかに超える盛上がりが徐々に形成され, ある程度の時間で定常な状態に達する.
ある程度相互作用をした後の様子が Fig.

$1(\mathrm{b})$に示されている. この段階では. 新しく出来た孤立波がこれから相互作用する孤立波の延長上に伸ひているように見える. 実際, 衝突後にできた孤立波をもう少し詳しく見てみると, Fig .l(c) に峰に沿った振幅を示すが, 峰に沿ってほぼ振幅が一定であり, 生成された孤立波が 1 次元的である事が分かる. また, Fig .l(d) は峰の伸ひる方向と直角の方向に断面を取って孤立波の外形を見たものであるが, その波形が MKd好愁螢肇鵑犯鷯錣砲茲

一致する事がわかる. すなわち, 相互作用の十分進んだ定常状態では, 初期に置いた 2 つのソリトンは独立に伝播しているように見える.
このように, 進行方向の成す角度が大きい場合には 2 っのソリトンは相互作用によってあまり影響を受けず, 定常状態では 2 つの無限に伸ひる 1 次元的ソリトンが重ね合わされたパターンが安定に伝播するという結果が得られた.

2 本研究では 2 つの孤立波の振幅を同じに取る.

128$X$(a)

$x$

(b)
Fig 1(cont)

129

054
052
0.5

$\approx$048046210220

(。)

230

240

250

1020$X$

0.6
0.5
0.4

0.3
$\approx$

0.2

0.1
0-0.1
(d) -20-100$S$

Fig 1:Development of the wave interaction for $A=0.5,$ $\Omega=4.0$ . (a) Initial condition $T=0$ .

profile at $T=2.925$ .

$(\mathrm{c})$Amplitude of the solitary wave after the interaction at $T=2.925$ .

the solitary wave after the interaction at $T=2.925$ . Axis
$S=\mathrm{O}$$S$$(\mathrm{d})$$(\mathrm{b})$Wave

Proffle ofis indicated in the Fig 1(b) and the origin

is located on $(X, \mathrm{Y})=(230.0,44.9)$ . Solid line: Numerical result, Dashed line: MKdV Soliton

13042

進行方向の成す角度が小さい ($\Omega$が小さい) 場合の結果

一方, Fig 2 は $\Omega=1$ , つまり進行方向のずれが小さい場合である. この場合,

Fig .2(b) に示すように,相互作用する場所には大きな振幅を持つ 1 次元的な孤立波 (以後 stem と呼ぶ) が形成される.

この stem の振幅は相互作用が進むにつれて大きくなり, やがてある一定値を取るようになる. 一方 $\Omega=4$ の場合と同様,相互作用する前の波に対応する孤立波が後方 ( $X$ が負の方向) に生成される.の場合に比べて, 振幅が小さ$\text{く}$しカル, それらは $\Omega=4$, 波の進行方向が対応する孤立波に比べてずれており,stem の存在による位相ずれも見られる. さらに詳しく見てみると, 振幅は峰に沿って大きく変化しており (Fig .2(c)), また, その形は MKdV ソリトンの形からずれている (Fig .2(d)). これらの事から, 後方の孤立波は 1 次元的 MKdVソリトンとみなす事ができない. 一方, stem についてであるが, Fig .2(e) にあるように, 振幅が一定となった所での stem は定常に伝わる MKdV ソリトンと見なす事ができる.

$\mathrm{K}\mathrm{P}$方程式などのソリトン共鳴では,1 点から伸ひる 3 つのソリトンが波数ベクトルに関してある共鳴条件を満たすという特徴があり, この場合にもそのような関係があるかどうかには興味がある. しかし, ソリトンの伸ひている点 (stem の端) 後方
ではさまざまな波が励起されているため, なにか本質的な関係が存在するかどうかを判断する事は難しく,今後の課題である.

43stem の性質に関する解析結果
進行方向の角度が小さい場合に相互作用によって形成される stem の性質について示す.

次に,Fig.3 に, さまざまな$\Omega$について波の進行に伴う stem の$\mathrm{Y}$方向の伸ひ$L_{\epsilon}$を示す. ここで伸ひは, 領域の対称線から stem の端までの長さとし, stem の端は, 相互作用する孤立波と繋がっていて大きさが $A$になった場所で決めた.

$\Omega=4$

では stem はほとんど形成されず,$\Omega$が大きいほど,stem の伸ひ方が大きい事がわかる.
しかしながら,ざまな$\Omega$stem の振幅については違った傾向がある. Fig.4 はある程度孤立波が進んだ場所でのさまに対する stem の振幅 um。を示す. ここでは方程式が変換 $(u, X, \mathrm{Y}, T)arrow(u/a, X/a, \mathrm{Y}/a^{2},T/a^{3})$に対して不変である事から,ケーリングで見えなくなる.の孤立波があろう.

$\overline{\Omega}$$\mathrm{Y}$$\overline{\Omega}\equiv\Omega/A$$\overline{\Omega}$と $\overline{u}_{\max}\equiv u_{\max}/A$ を比較する事にする.

が小さい時には, u-m。は軸に平行にある状態なので,がある程度大きくなると,$\overline{\Omega}$大きくなるほど弱まり, ある程度大きな$\overline{\Omega}arrow 0$$\overline{\Omega}$

このため $A$ の違いはスに比例して大きくなる. なお,$\overline{\Omega}=0$は1つで u-m。が 1 に近づいているように見えるのは自然での増加に対して u-。へは減少するようになるが, その傾向は$\overline{\Omega}$

では, 一定値を取るようになる.

この値が$u- \mathrm{m}\text{。}=2$$\overline{\Omega}$がである事と先に示した $\Omega=4$ の数値計算の結果から, 非常に角度が大きい場合の 2 つの孤立波の相互作用はほとん

ど変化の無い重ね合わせでないかとの予想が得られる.
また,

u-m。の変化が増加から減少に変わる部分は非常に興味深い.

可積分である$\mathrm{K}\mathrm{P}$

方程式でも同様なstem の角度依存性が見られるが, さらに最大値 (初期孤立波の 4 倍) 及ひその時の角度が厳密に求められている. 振幅の大きい stem を安定に数値計算する事が難しいため,しいが, 例えば$\mathrm{K}\mathrm{P}$

この問題について明らかにするのは難の最大値を超えるかどうかなどの視点から今後調べる予定である.

131$X$(a)$X$(b)
Fig 2(cont.)

132$\approx$

1.4
1.2
1$\approx$

0.8
0.6
0.4
0.20——————————–-0.2(e) 280

290

310

300

320

$X$

Wave
Fig.2 :Development of the wave interaction for $A=0.5,$ $\Omega=1.0$ . (a) Initial condition $T=0$ .
profile at $T=25.005$ . (c) Amplitude of the solitary wave after the interaction at $T=25.005$ . (d) Profile$(\mathrm{b})$of solitary wave after the interaction at $T=25.005$ . Axis is indicated in the Fig.2(b) and the origin
is located on $(X,\mathrm{Y})=(280.0, 22.3)$ . Solid line: Numerical result, Dashed line: MKdV Soliton$S$$S=\mathrm{O}$solution. (e) Proffle of the stem wave at $T=25.05$ . Solid line: Numerical raeult, Dashed line: MKdVsoliton solution.

133

30

2520
$\mathrm{d}^{\mathrm{O}\mathrm{l}}15$

1050010

30

20

50

40

60

70

$X-X_{0}$Fig 3:Development of stem length for various $\Omega.A=0.5$ .

4350.50.75$\mathrm{F}$$+$$\fbox$

310$\mathrm{o}^{\copyright}$$\mathrm{Q}$$+$$\ovalbox{\tt\small REJECT}^{\triangleleft}\mathrm{s}$$-$$2.52$

ffl.

$+$

$+$

$\theta$

1.5$\bigoplus_{-}$

1050246810$\alpha A$

Fig.4 :Dependence of the stem amplitude at $X-X_{0}=25.0$ on$\Omega$.

1345

まとめ
本研究では,

MKP 方程式を用いて MKdV ソリトンの 2 次元的相互作用を調べた. その結果, 進行方向の角度の違いによって相互作用が非常に小さい場合と, 初期のソリトンより高い振幅を持つ第 3 のソリトンが生成される場合とがある事が明らかになった. ソリトンの出会う角度が近づくにつれて stem の長さはより長くなるが, 定常振幅は, ある程度の範囲に収まる. これらの傾向は今まで調べた空間 2 次元の方程式での結果とある程度類似している. この事を統一的に説明できるかどうかは今後の課題である.
今回, 物理系が弱 2 次元という仮定の元に孤立波の相互作用を調べたが, スケーリングについて仮定の無い場合について基礎方程式を摂動論によって調べる事は興味深い. その場合には個々の孤立波が独立に進むという予想が本研究の一つの結果としてあるが,これが正しいのかを解析的に明らかにできないかと考えている.

参考文献
1) 及川正行: 九州大学応用力学研究所所報 60 (1984) 467.
2) Y. Matsuno: Proc. R. Soc. Lond. A454 (1998) 835.
3) H. Tsuji and M. Oikawa: J. Phys. Soc.

$\mathrm{J}\mathrm{p}\mathrm{n}$.

$62$(1993) 3881.

4) H. Tsuji and M. Oikawa . Fluid Dyn. Ilae. 29 (2001) 251.
$\cdot$5) B. B. Kadomtsev and$\mathrm{V}.\mathrm{I}$. Petviashvfli Sov. Phys. Dokl. 15 (1970) 539.

6) S. Watambe: J. Phys. Soc.

$\mathrm{J}\mathrm{p}\mathrm{n}$. 53 (1984) 950.

7) Chen Yongze and P. L.-F. Liu: Wave Motion, 27 (1998) 321.