(2+1)-мерные нелинейные интегрируемые уравнения: асимптотика и возмущения

О.М.Киселев, ok@ufanet.ru
Институт математики, Уфа

В качестве основных представителей (2+1)-мерных уравнений будут рассматриваться уравнения Деви-Стюартсона (ДС):

iAt+x2A+a2y2A+2k|A|2A+Ap = 0,
x2p-a2y2p = -4kx2(|A|2),   k = ±1;
(1)
Кадомцева-Петвиашвили (КП):
x(t u+6ux u+3x u) = -3a22y u;
(2)
а также Ишимори:
t ®
S
 
+ ®
S
 
×(x2 ®
S
 
+a2y2 ®
S
 
)+x wy ®
S
 
+y wx ®
S
 
= 0
x2 w-a2y2 w+2a2 ®
S
 
(x ®
S
 
×y ®
S
 
) = 0,
(3)
®
S
 
= (S1,S2,S3),   ®
S
 
2
 
= 1.
Во всех этих уравнениях a2 = ±1.

Основное свойство этих уравнений - интегрируемость. Это позволяет исследовать асимптотики решений на больших временах и строить асимптотики решений возмущенных уравнений.

Цели докдада:

1.Рассказать о структурной неустойчивости алгебраического солитона уравнения ДС-2 по отношению к возмущению начального условия.

2. Построить асимптотику при t®Ґ бессолитонного решения уравнения ДС-2.

3. Сформулировать результат об асимптотике решения уравнений Ишимори-1 и КП-2.

4. Построить формальное асимптотическое решение возмущенного уравнения ДС-1.

1   Алгоритм интегрирования ДС-2

Рассмотрим уравнения ДС-2 в виде:

it q+2(2z+2[`z])q+(g+
g
 
)q = 0,
[`z] g = z |q|2,    z = x+iy.

Уравнение ДС-2 в методе обратной задачи связано с двумерной системой Дирака ( Fokas, Ablowitz (1984)):

ж
з
и
[`z]
0
0
z
ц
ч
ш
f = 1
2
ж
з
з
з
з
и
0
q(z,t)
-
q(z,t)
 
0
ц
ч
ч
ч
ч
ш
f,
(4)
ж
з
и
exp(kz)
0
0
exp(
kz
 
)
ц
ч
ш
f(k,z)||z|®Ґ = ж
з
и
1
0
0
1
ц
ч
ш
.
Непрерывная часть данных рассеяния для этой задачи:
b(k) = -i
4p
у
х


C 
dzЩd _
z
 
 q(z)f22exp(-kz),
Пусть существует k = k0 такая, что краевая задача неразрешима. Тогда добавляют еще и дискретную часть:
{k0,n,m}.
Если функция q удовлетворяет уравнениям ДС-2, тогда данные рассеяния очень просто зависят от времени:
t b = 2i( _
k
 
2
 
+k2)b,   t n = 2i( _
k
 
2
0 
-k02)n,   
m(t) = m0+4ik0t.
(5)
Наряду с эллиптической системой Дирака матрица y удовлетворяет еще и краевой задаче по k:
ж
з
и
[`k]
0
0
k
ц
ч
ш
fT = ж
з
з
з
з
и
0
- _
b
 
(k,t)
b(k,t)
0
ц
ч
ч
ч
ч
ш
fT,
ж
з
и
exp(-kz)
0
0
exp(
-kz
 
)
ц
ч
ш
fT(k,z)||z|®Ґ = ж
з
и
1
0
0
1
ц
ч
ш
.
Зная зависимость данных рассеяния от времени и решив систему по k, можно получить q(z,t):
q(z,t) = -i
p
у
х


C 
dpЩd
p
 
 b(p;0)exp(2it(p2+ _
p
 
2
 
) -
pz
 
) f11
-
A(1)2(z;t)exp(-k0z)
 
.
Здесь A(1)2(z;t) - элемент матрицы - решения системы Дирака с нулевым граничным условием при k = k0.

Если система Дирака разрешима при всех k - решение бессолитонное; если существует k0 - солитонное.

2   Неустойчивость солитона

Если в формуле для q принять b є 0, тогда:

q(z,t) =
2
n
 
exp(-it(k02+ _
k
 
2
0 
)+k0z-
k0z
 
)

|z+4ik0t+m|2+|n|2
.

Это решение (солитон) было найдено Аркадьевым, Погребковым и Поливановым в 1989 году.

Рассмотрим уравнение ДС-2 с начальным условием в виде возмущенного солитона:

qe(z,0) = q(z,0)+eq1(z),    1 >> e > 0,    q1(z) О C0Ґ.
Вопрос: cохранится ли солитонная структура данных рассеяния при возмущении? Для ответа необходимо исследовать задачу о возмущении нулевого собственного значения.

Краевую задачу для системы Дирака перепишем в виде:

(I-G[q,k])f = E(kz),
(6)
где E(kz) = diag(exp(kz),exp([`kz])), интегральный оператор G[q,k]:
G[q,k]f = 1
4ip
у
х
у
х


C 
dxЩd
x
 
×                        
× ж
з
з
з
з
з
з
з
з
з
з
и
0
q(x)exp(k(z-x))
z-x
k _
q
 
(x)exp(
k(z-x)
 
)


z-x
 
0
ц
ч
ч
ч
ч
ч
ч
ч
ч
ч
ч
ш
f(x,k).
(7)
Пусть Xў -пространство непрерывных по z, [`z] ограниченных матриц с нормой:
||f|| =
sup
j,z 
(|Ejj(-kz)fj1(z,k)|+ |Ejj(-kz)fj2(z,k)|).
Здесь Elm(-kz) и flm(z,k) - элементы матриц E(-kz) и f соответственно.

Лемма 1 Оператор G[q,k], с непрерывным потенциалом, убывающим как |z|-2 при |z|®Ґ, является компактным в Xў.

Из альтернативы Фредгольма следует, что интегральное уравнение неразрешимо при некотором k = k0, если выполнены условия:

1) существует нетривиальное однородное решение краевой задачи

(I-G*[q,k0])B = 0,       E(zk0)B||z|®Ґ = 0,
(8)
где G*[q,k] - формально сопряженный к G[q,k] оператор относительно полуторалинейной формы:
(f(i),y(j))f = у
х
у
х


C 
dzЩd _
z
 
  ж
и
_
f
 
(z) f1i(z,k)
y1j
 
(z,k) +
+f(z) f2i(z,k)
y2j
 
(z,k) ц
ш
;

2) существует пара чисел (l,m) такая, что

(E(l)(kz),B(m))q 0,
где l = 1,2, m = 1,2,...,N, N - число линейно независимых решений B(m) интегрального уравнения (8).

Оба этих условия выполнены для солитонного потенциала q0(z). Однородное интегральное уравнение (8) имеет два линейно независимых решения B(1)(z) и B(2)(z), таких что

(E(1)(k0z),B(1))q0 = 0,   (E(1)(k0z),B(2))q0 = 4ip,
(E(2)(k0z),B(1))q0 = 4ip,   (E(2)(k0z),B(2))q0 = 0.
Эти решения имеют вид:
B(1)(z) = 1
|z+m|2+|v|2
ж
з
з
з
з
и
-
(z+m)
 
exp(-zk0)
-
n
 
exp(-
zk0
 
)
ц
ч
ч
ч
ч
ш
;
B(2)(z) = 1
|z+m|2+|v|2
ж
з
з
з
з
з
и
- 1
k
nexp(-zk0)
-(z+m)exp(-
zk0
 
)
ц
ч
ч
ч
ч
ч
ш
.

Если рассмотреть задачу о возмущении нулевого собственного значения

(I-G[qe,k])B = lB,
(9)
тогда окажется, что формальные асимптотики возмущенных собственных значений и собственных векторов имеют вид:
~
l
 
(e,k) = e 1
l
 
(k)+o(e),
(10)
~
B
 
(z,k,e) = 0
B
 
(z,k)+e 1
B
 
(z,k)+o(e),
(11)
где k = (k-k0)e-1. Коэффициенты в разложении по e cобственных векторов удовлетворяют уравнениям:
(I-G*[q0,k0]) 0
B
 
= 0,
(12)
(I-G*)[q0,k0] 1
B
 
= 1
F
 
,
(13)
Главный член разложения строится в виде:
0
B
 
(z,k) = a(k)B(1)(z)+b(k)B(2)(z)
(14)
Параметры a и b выбираются так, чтобы уравнение для B было разрешимо. В результате, оказывается:
1
l
 


1 
(k) 1
l
 


2 
(k) =
|4ip
k
 
+Q2|2+|Q1|2

|n|2
.
(15)

Здесь Q1 и Q2 линейные функционалы, зависящие только от возмущения q1(z).

Две ветви собственных значений в зависимости от комплексного параметра k локально, вблизи k0 имеют вид двуполостного параболоида:

Гадыльшин, Киселев 1999:

Утверждение 1 [Неустойчивость солитона]    Возмущение потенциала в системе Дирака приводит к расщеплению двукратного полупростого нулевого собственного значения на два простых, ненулевых. Следовательно, солитонная структура данных рассеяния неустойчива по отношению к возмущению потенциала q(z).

Решение системы Дирака с возмущенным солитонным потенциалом.

Формальную асимптотику fe будем строить в виде

fe(z) = f0(z)+ef1(z)+...,
(16)
где f0 - первый столбец матрицы - решения невозмущенной задачи рассеяния (4) при q(z,0) = q0(z).

Используя формулу Коши-Грина, запишем граничную задачу (4), в форме интегрального уравнения

(I-G[qe,k])fe = E1.
(17)

Приравнивая коэффициенты при одинаковых степенях e, получаем уравнения для f0 и f1:

(I-G[q0,k])f0 = E1,
(18)
(I-G[q0,k])f1 = G[q1,k]f0,
(19)
Эти уравнения разрешимы всюду кроме точки k0.

При |k-k0| < 2e1/2 асимптотическое решение задачи (4) будем искать в виде

fe(k,z) = E(kz)(e-1F-1(k,z)+F0(k,z)+...),
(20)
где
F-1(k,z) ~ - 1
k
exp{-k0z}A(1),      k®Ґ,
(21)
F0(k,z) = ж
з
и
1
0
ц
ч
ш
+O(z-1),       z®Ґ,    k®Ґ.
(22)

Обозначим

G[h,k] = E(-kz)G[ ^
h
 
,k],       ^
h
 
(z,k) = h(z) exp{
kz
 
}.
(23)
Подставляя (20) в (17) и учитывая краевые условия (21) и (22), получаем следующие уравнения:
(I- G[q0,k0])F-1 = 0,
(24)
(I- G[q0,k0])F0 = I1+ G[q1,k0]F-1+
kk G[q0,k0]F-1+
k
 
[`k] G[q0,k0]F-1,
(25)
где I1 - первый столбец единичной матрицы. Принимая во внимание (23), легко показать, что функция
F-1(k,z) = E(-k0z)A(z)a(k),
(26)
где A - матрица со столбцами A(j), j = 1,2, является решением однородного уравнения (24) для любого вектора a(k).

Условия разрешимости для уравнения

(I- G[q0,k0])W = F
имеют вид
(E(k0z) F,B(i))q0 = 0,       i = 1,2.
(27)
Вектор a определяется из условия разрешимости (27) для (25). Заметим, что из определений G[h,k] и G[h,k] следует, что
kk G[h,k0]+
k
 
[`k]G[h,k0] = G[ ~
h
 
,k0],   
(28)
где    ~
h
 
(z,k) = h(z)(
zk
 
-zk).
Для компонент a получим:
a(1)(k)
= -
(
Q2+k
 
)

|Q1|2+|Q2+k|2
,
a(2)(k)
= - Q1
|Q1|2+|Q2+k|2
.
(29)
В результате, функция F-1, определенная в (26) и (29), удовлетворяет условиям согласования (21) и не имеет сингулярностей по k.

Таким образом, для невырожденного возмущения (Q1 0) асимптотическое решение имеет вид (16) при |k-k0| > e1/2 и имеет вид (20), (26), (29) при |k-k0| < 2e1/2.

В результате можно получить данные рассеяния ( Гадыльшин, Киселев, 1996):

be(k) = -i
4p
у
х


C 
dzЩd _
z
 
q(z)f22exp(-kz),
be(k) ~ e-1B-1 ж
з
и
k-k0
e
ц
ч
ш
+B0 ж
з
и
k-k0
e
ц
ч
ш
 
при |k-k0| < 2e1/2,
             B-1(k) = -

Q1
 

| Q1|2+|Q2+k|2
,
Q1,2 = const 0.
                         be(k) ~ eb1(k)   при |k-k0| > e1/2.

Следующий шаг - решение [`D]-задачи для возмущенных данных рассеяния:

ж
з
и
[`k]
0
0
k
ц
ч
ш
fT = ж
з
з
з
з
и
0
k _
b
 

e 
(k,t)
be(k,t)
0
ц
ч
ч
ч
ч
ш
fT,
ж
з
и
exp(-kz)
0
0
exp(
-kz
 
)
ц
ч
ш
fT(k,z)||z|®Ґ = ж
з
и
1
0
0
1
ц
ч
ш
.
Формальная асимптотика этой задачи при e®0, равномерно при t = O(e-1) может быть построена похожим образом. В результате решение можно представить в виде интеграла типа Фурье:
qe(z,t) = -i
p
у
х


C 
dpЩd
p
 
 be(p)exp(2it(p2+ _
p
 
2
 
) -
pz
 
) f11.

Утверждение 2 [О возмущении солитона]    ( Гадыльшин, Киселев, 1999) Возмущение приводит к бессолитонной структуре данных рассеяния. Однако, асимптотическое решение сохраняет солитонную форму главного члена с модулированным параметром: me = m0+e2tpQ2 при t = O(e-1).

3  Асимптотика бессолитонного решения ДС-2 при больших значениях времени

Сформулируем утверждение об асимптотике решения уравнений ДС-2:

Теорема 1 Пусть данные рассеяния задачи состоят из функции b(k) такой, что она сама и все ее производные по k и [`k] до второго порядка - гладкие интегрируемые в R2, функция k b(k) интегрируема в R2, и b(k) удовлетворяет требованию:

1
p

max
z О C 
у
х
у
х


R2 
dldm  |b(l+im)|
|l+im-z|
  < 1,
(30)
тогда асимптотика при t®Ґ решения ДС-2 имеет вид:
q(z,t) = t-1  1
2
 b ж
и
iz
4t
ц
ш
 exp ж
з
и
it
8
(x2+
x
 
2
 
) ц
ч
ш
 + O(t-5/4),
(31)
g = -t-2 V.P. у
х
у
х


C 
dvЩd _
v
 
 |b ж
и
iv
4t
ц
ш
|2

2ip(v-z)2
 + O(t-9/4),
(32)
где x = z/t, равномерно по z О C.

Замечание 1 Если функция b(k) не удовлетворяет условию (30), тогда удается построить формальное асимптотическое при t®Ґ решение ДС-2:

~
q
 
(z,t)
= t-1  1
2
 b ж
и
iz
4t
ц
ш
 exp ж
з
и
it
8
(x2+
x
 
2
 
) ц
ч
ш
,
~
g
 
(z,t)
= -t-2 V.P. у
х
у
х


C 
dvЩd _
v
 
 |b ж
и
iv
4t
ц
ш
|2

(v-z)2
.
Функции [q\tilde] и [g\tilde] удовлетворяют первому уравнению в ДС-2 с точностью O(t-3/2), второму - с точностью O(t-5/2).

Следствие 1 На больших временах, t = O(e-g), g > 1 возмущенный солитон диспергирует

qe ~ (te)-1B-1 ж
з
и
iz
4t
ц
ч
ш
exp ж
з
и
i(z2+ _
z
 
2
 
)

8t
ц
ч
ш
.

В этом разделе построено асимптотическое решение обратной задачи при больших значениях t методом согласования асимптотических разложений.

В решении удобно выделить осциллирующие экспоненты. Для этого примем обозначения:

f11(k,z,t)
= A(k,z,t)exp(kz),
f12(k,z,t)
= B(k,z,t)exp(
kz
 
).
(33)

Перепишем обратную задачу для функций A и B. В результате получим систему уравнений

[`k]A = - _
b
 
exp(-itS)B,
k B = bexp(itS)A,
    ж
з
и
A
B
ц
ч
ш
к
к
к


|k|®Ґ 
= ж
з
и
1
0
ц
ч
ш
(34)
с фазовой функцией специального вида: S = 2(k2+[`k]2)-i(kz-[`kz])/t.

Структура формального асимптотического решения задачи (34) определяется фазой S осциллирующей экспоненты в правой части (34). Асимптотическое по t при t®Ґ решение (34) имеет составной характер. Вне некоторой окрестности стационарной точки фазы S формальное асимптотическое решение строится по обратным степеням t (внешнее разложение).

В малой при t®Ґ окрестности стационарной точки построено так называемое внутреннее разложение. Здесь в качестве калибровочной последовательности используется tn/2,    n = 1,2,....

Существует некоторая область в окрестности стационарной точки, где пригодны оба асимптотических разложения, как внутреннее, так и внешнее. Это общий факт в методе согласования асимптотических разложений. Здесь он используется для однозначного определения коэффициентов разложения.

В этом разделе доказано утверждение:

Теорема 2 Пусть b(k) - гладкая по k,  [`k]  и такая, что b(k) и все ее производные по k  и [`k] до второго порядка гладкие интегрируемые в R2 функции. Тогда формальное асимптотическое по mod(o(t-1)) решение задачи (34) при t®Ґ имеет составной характер. При |k-[(ix)/4]| і Ct-1/2+g, для "g  > 0, где x = z/t:

ж
з
з
з
з
и
~
A
 
~
B
 
ц
ч
ч
ч
ч
ш
= ж
з
и
1
0
ц
ч
ш
+ t-1 ж
з
з
з
з
з
и
1
A
 
1
B
 
ц
ч
ч
ч
ч
ч
ш
,
(35)
где B = B0([`k],x)+B1(k,x)exp(itS), поправки A, B0 и B1 определены в (), (), ().

При | k-[(ix)/4]| Ј Ct-d для "d > 0:

ж
з
з
з
з
и
~
A
 
~
B
 
ц
ч
ч
ч
ч
ш
= ж
з
и
1
0
ц
ч
ш
+t-1/2 ж
з
з
з
з
з
и
0
1
b
 
(l,x)exp(itS)
ц
ч
ч
ч
ч
ч
ш
+
+ t-1 ж
з
з
з
з
з
и
2
a
 
(l,x)
2
b
 
(l,x)exp(itS)
ц
ч
ч
ч
ч
ч
ш
,
где l = Цt (k-[(ix)/4]), поправки ab,    n = 1,2 определены в (), () и ().

Доказательство теоремы 3 состоит из двух частей - построения, так называемого внешнего формального асимптотического по t при t®Ґ решения (33) в области вне малой при t®Ґ окрестности точки k = ix/4 и построения внутреннего разложения (34) в окрестности точки k = ix/4.

Решение задачи (34) будем искать в виде:

A(k,x,t) = 1 +  t-1  1
A
 
(k,x)  +             
+t-2 2
A
 
(k,x)exp(-itS) + ...,
(36)
B(k,x,t) = t-1  1
B
 
 + t-2 2
B
 
(k,x)exp(itS) + ... .
(37)

Подставим (36), (37) в (34), приравняем коэффициенты при одинаковых степенях t-1. В результате получим:

ik 1
B
 


1 
= b(k),    1
B
 


1 
||k|®Ґ = 0;
(38)
-i 2
A
 
exp(-itS)[`k]S  + [`k] 1
A
 
=
= - _
b
 
  1
B
 


1 
 -  _
b
 
1
B
 


0 
exp(-itS),
(39)
1
A
 
||k|®Ґ = 0;
(40)
ik 2
B
 
+k 1
B
 


1 
= b(k) 1
A
 
,    2
B
 
||k|®Ґ = 0.
(41)
Уравнение (38) разрешимо, когда k S 0. Из представления S легко видеть, что k S = 0 при k = [(ix)/4]. Тогда, с точностью до убывающего при |k|®Ґ решения однородного уравнения для B из (34), решение уравнения (38) при k [(ix)/ 4] имеет вид:
1
B
 


1 
= b(k)
-ik S
.
(42)

Не определенное здесь решение однородного уравнения для B - аналитическая функция от [`k] - получено ниже, после анализа асимптотического разложения в окрестности точки k = [(ix)/4].

Уравнение для A определяется неосциллирующими членами уравнения (38). С помощью формулы Коши-Грина получим:

1
A
 
(k,x) = -1
8ip
у
х
у
х


C 
dpЩd _
p
 

k-p
|b(p)|2
p- ix
4
.
(43)

Поправка A определяется осциллирующими членами (38):

2
A
 
=
_
b
 
(k) 1
B
 


0 

i[`k]S
.

В этой формуле присутствует B0, поэтому A может быть построена после анализа асимптотики в окрестности k = [(ix)/4].

Установим область, в которой пригодно асимптотическое разложение (36), (37). В этой области для слагаемых в формуле (37) выполняется условие:

t-2 2
B
 
(k,x) = o(t-1 1
B
 


1 
(k,x)).

Из формулы (42) легко видеть, что функция B1 в точке k = ix/4 имеет полюс первого порядка по k. Решение уравнения (41) - функция B - в точке k = ix/4 имеет полюс третьего порядка по k. Поэтому асимптотическое разложение (35) пригодно при |k-[(ix)/4]| і C t-1/2+g, для "g > 0.

Построим внутреннее разложение - справедливое в малой окрестности точки k = [(ix)/4].

Перейдем к растянутым координатам:

l = Цt(k- ix
4
).
В этих координатах фаза экспоненты в системе (34) имеет вид:
itS = 2i(l2+ _
l
 
2
 
)+itS0.
Здесь
S0 = S|k = [(ix)/4] = 1
8
(x2+
x
 
2
 
).
Оператор дифференцирования по k в системе уравнений (34) перепишем через дифференцирование по l. В результате система уравнений (34) примет вид:
Цt[`l] A
= - _
b
 
exp(-2i(l2+ _
l
 
2
 
)-itS0)B,
Цtl B
= b exp(2i(l2+ _
l
 
2
 
)+itS0)A.
(44)

Формальное асимптотическое решение системы (44) будем искать в виде:

A = 1+t-1/2 1
a
 
(l,x) + t-1 2
a
 
(l,x)+...,      
B = ж
и
t-1/2 1
b
 
(l,x) + t-1 2
b
 
(l,x)+... ц
ш
×
exp(2i(l2+ _
l
 
2
 
)+itS0).
(45)

Подставим (45) в (44). Приравняем коэффициенты при одинаковых степенях t, в результате получим уравнения для ab:


_
l
 
  1
a
 
= 0,      

l
  1
b
 
 + 4il  1
b
 
= b ж
и
ix
4
ц
ш
,
(46)
и уравнения для ab:

_
l
 
  2
a
 
= -
b ж
и
ix
4
ц
ш
 
  1
b
 
,            

l
  2
b
 
 + 4il  2
b
 
= b ж
и
ix
4
ц
ш
  1
a
 
+ l b1+ _
l
 
 b2.
(47)
Здесь b1 = k b|k = [(ix)/4],    b2 = [`k]b|k = [(ix)/4].

Область |l| Ј const t1/2-d"d > 0, где справедливо внутреннее разложение, при d < 1/2 пересекается с областью |k-[(ix)/4]| і t-1/2+gconst  при g < 1/2, в которой пригодно внешнее разложение. Следуя методу согласования асимптотических разложений, асимптотики при |l|®Ґ для (46), (47) можно получить из согласования внутреннего и внешнего асимптотических разложений в области, где пригодны оба разложения.

Асимптотика B1 при k®[(ix)/4]:

1
B
 


1 
= -1
i(4k-ix)
b ж
з
и
ix
4
ц
ч
ш
 -  1
4i
 b1  - 
(k- ix
4
)

4i(k- ix
4
)
 
 b2 + 
+ O(|k- ix
4
|).
(48)

Вычислим асимптотику A  при k®[(ix)/4].

1
A
 
= -1
k- ix
4
  1
8ip
у
х
у
х


C 
dpЩd _
p
 
|b(p)|2 й
к
к
к
к
к
л
1
k-p
- 1
p- ix
4
щ
ъ
ъ
ъ
ъ
ъ
ы
.
Подынтегральная функция в интеграле для A представлена в виде разности. Рассмотрим этот интеграл как разность интегралов от соответствующих функций. В первом из получившихся интегралов сделаем замену переменной r = p-(k-[(ix)/4]). Представим числитель подынтегральной функции в первом интеграле в виде отрезка ряда Тейлора в окрестности точки r. В результате получим:
1
A
 
= -1
2ip
у
х
у
х


C 
drЩd _
r
 

ix
4
-r
1
4 i
(r|b(r)|2 + 

k- ix
4
 

k- ix
4
 [`r]|b(r)|2) =
=   1
4i
 

k- ix
4
 

k- ix
4
 |b ж
и
ix
4
ц
ш
|2  + C ж
и
ix
4
ц
ш
 + O(|k- ix
4
|),
(49)
где
C ж
и
ix
4
ц
ш
= -1
2ip
у
х
у
х


C 
drЩd _
r
 

ix
4
-r
  1
4i
 r|b(r)|2.

Из требования согласования асимптотических разложений (35) и (36) и формул (48), (49) получаются краевые условия:

1
a
 
= 0,    при    |l|®Ґ,
1
b
 
= O(l-1),    при   |l|®Ґ;
(50)
2
a
 
=
_
l
 

4il
 |b ж
и
ix
4
ц
ш
|2 + C ж
и
ix
4
ц
ш
+O(|l|-1),     при   |l|®Ґ
2
b
 
= 1
4i
 b1 - 
_
l
 

4il
 b2 + O(|l|-1) ,    при   |l|®Ґ.
(51)

Решение первого из уравнений (46) совместно с краевым условием из (50) имеет вид:

1
a
 
= 0.
(52)

Ограниченное решение уравнения (46) с краевым условием из (50) получается с помощью формулы Коши-Грина:

1
b
 
= exp(-2i(l2+ _
l
 
2
 
)) b ж
и
ix
4
ц
ш
×
-1
2ip
у
х
у
х


C 
dnЩd _
n
 


n-l
 
exp(2i(n2+ _
n
 
2
 
)).
(53)

Асимптотика b при |l|®Ґ имеет вид:

1
b
 
||l|®Ґ =
b( ix
4
)exp(-2i(l2+ _
l
 
2
 
))

2 _
l
 
 - 
b( ix
4
)

4il
+O(|l|2).
(54)
Из согласования асимптотик (48), (54) и (35) получим:
1
B
 


0 
( _
k
 
, ix
4
) =
b( ix
4
)

2
(k- ix
4
)
 
.
(55)

Уравнения (47) с краевыми условиями (51) решаются с помощью формулы Коши-Грина:

2
a
 
= C ж
и
ix
4
ц
ш
 + 
b ж
и
ix
4
ц
ш
 
  -1
2ip
у
х
у
х


C 
dnЩd _
n
 

n-l
1
b
 
(n, _
n
 
),
2
b
 
= exp(-2i(l2+ _
l
 
2
 
))  -1
2ip
у
х
у
х


C 
dnЩd _
n
 


n-l
 
×
×exp(2i(n2+ _
n
 
2
 
)) (nb1+ _
n
 
b2).
(56)

Теорема 3 об асимптотике решения задачи (34) доказана.

3.0.1   Оценка остатка асимптотики

Можно доказать, что формулы (35), (36) определяют асимптотику решения задачи (34). Для этого рассматривается система интегральных уравнений, соответствующая краевой задаче (34). При ограничении (30) на функцию b(k) интегральный оператор в этой системе оказывается сжимающим. С помощью формул (35), (36) строится составное асимптотическое разложение (35). Малость остатка асимптотики следует из свойства интегрального оператора и малости невязки, которая получается при подстановке составного разложения в систему интегральных уравнений.

Справедливо утверждение:

Теорема 3 Пусть b(k) - гладкая функция такая, что все ее производные по k и [`k] до второго порядка - гладкие интегрируемые в R2, и выполнено условие:

1
p

max
z О C 
у
х
у
х


R2 
dldm  |b(l+im)|
|l+im-z|
  < 1,
тогда решение [`D]-задачи существует. Его асимптотика при t®Ґ определяется формулами (35), (36) с точностью O(t-5/4). Производные по z и [`z] остатка асимптотики - непрерывные равномерно ограниченные при k,z О C - равны O(t-5/4).

3.0.2   Асимптотика решения ДС-2

Подставим асимптотику решения [`D]-задачи в формулу формулу для решения ДС-2.

q(z,t) =
-i
p
у
х
у
х


C 
dpЩd
p
 
 b(p)exp(itS)×            
× ж
з
и
1+t-1 ~
A
 

1 
(p,x,t)+t-5/4 f1(p,x,t) ц
ч
ш
,
(57)
где
~
A
 

1 
(p,x,t) = c  1
A
 
 + (1-c 2
a
 
.

Представим интеграл (57) в виде суммы интегралов.

q(z,t)
= -i
p
у
х
у
х


C 
dpЩd
p
 
 b(p)exp(itS)-
- t-1 i
p
у
х
у
х


C 
dpЩd
p
 
 b(p)exp(itS) ~
A
 

1 
(p,x,t)-
-t-5/4 i
p
у
х
у
х


C 
dpЩd
p
 
 b(p)exp(itS)f1(p,x,t).
Последнее слагаемое равно O(t-5/4). Интеграл во втором слагаемом разобьем на два. При |p-[(ix)/4]| і 2t-1/4 этот интеграл оценивается с помощью интегрирования по частям и имеет порядок t-1. При |p-[(ix)/4]| < 2t-1/4 легко видеть, что этот интеграл равен O(t-1/2). В результате получим:
q(z,t) = -i
p
у
х
у
х


C 
dpЩd
p
 
 b(p)exp(itS) + O(t-5/4).
Асимптотика по t интеграла в этой формуле вычисляется методом стационарной фазы. В результате получается формула (31).

3.1   Решение уравнений Ишимори-1

Для уравнений Ишимори-1 (в (3) a2 = 1) метод обратной задачи рассеяния изложен, например, в работе Конопельченко и Маткаримова. Здесь воспользуемся формулой, связывающей решения уравнений Ишимори и решения вспомогательной линейной задачи (4) для уравнений ДС-2.

Обозначим [S\tilde] = [S\vec][(s)\vec], где [(s)\vec] = (s1,s2,s3), а s1,2,3 - матрицы Паули. Тогда

~
S
 
= -f-1s3f,
(58)
где f - решение (4) с потенциалом, q.

Асимптотика решения уравнений Ишимори-1 при t®Ґ с краевым условием [S\tilde]||z|®Ґ = s3 может быть получена из формулы (58). После подстановки асимптотики решения краевой задачи (4), построенной в предыдущем разделе, получим:

~
S
 
= ж
з
з
з
з
з
и
-1
2
t1/2
1
b
 
2
t1/2
1
b
 
1
ц
ч
ч
ч
ч
ч
ш
+o(t-1/2),
(59)
где b(l)- решение краевой задачи:
l 1
b
 
+4il 1
b
 
= b(ix/4),    1
b
 
||l|®Ґ = 0.
Функция b(k) удовлетворяет условиям теоремы 3

4   Асимптотика уравнения КП-2

Уравнение КП-2 связано с [`D]-задачей:

[`k]f = y b(- _
k
 
)exp(itS),
ky = -f b(k)exp(-itS);
ж
з
и
f
y
ц
ч
ш
||k|®Ґ = ж
з
и
1
1
ц
ч
ш
.
(60)
Здесь S = 4(k3+[`k]3)+(k+[`k])x-i(k2-[`k]2)h,   x = x/t, h = y/t, функция b(k) - неаналитическая функция комплексной переменной k О C.

Решение уравнения КП-2 ( Ablowitz, Bar Yaakov, Fokas):

u(x,y,t) = x у
х
у
х


C 
dkЩd _
k
 
 b(k)y(k,x,y,t)exp(itS).

Исследование решения [`D]-задачи для уравнения КП-2 при t®Ґ более сложно, чем проведенное выше для уравнения ДС-2. Это связано, во-первых, - с тем, что [`D]-задача для класса гладких и достаточно быстро убывающих по (x,y) начальных условий связана с разрывной на мнимой оси функцией b(k). Во-вторых, - структура быстро осциллирующей фазы S в коэффициентах эллиптической системы для f,y более сложная. Наряду с простыми, она содержит и вырожденную стационарную точку.

Теорема 4 Пусть (1+|k|)b(k) О L1ЗC, ab О L1ЗC при В(k) 0, |a| Ј 2 и выполнено условие


sup
z О C 
у
х
у
х


R2 
dkdl к
к
к
b(k+il)
k+il-z
к
к
к
  <  p,
тогда решение задачи Коши для уравнения КП-2 существует при "t > 0. Его асимптотика при t®Ґ носит составной характер:

при 12x+h2 < 0 и t1/3|12x+h2| >> 1:

u(x,y,t) =                                     
= -4t-1 p
12i   ж
Ц

-h2-12x
 
f0 ж
з
и
1
2
  ж
Ц

-h2-12x
 
+ ih
12
ц
ч
ш
×
×exp ж
з
и
-11it   ж
Ц

- y2
t2
-12 x
t
 
ц
ч
ш
+c.c.+o(1).
при 12x+h2 > 0 и t1/3|12x+h2| >> 1:
u = o(t-1);
при |12x+12h2| << 1:              u(x,y,t) =
= 8it-1   ж
Ц

p
 
f(ih/12) ж
з
и
у
х
Ґ

0 
dp1   __
Цp1
 
cos ж
з
и
8p13-z p1 ц
ч
ш
+
+ у
х
Ґ

0 
dp1   __
Цp1
 
sin ж
з
и
8p13-z p1 ц
ч
ш
ц
ч
ш
 + o(t-1).
Здесь x = x/t, h = y/t, z = 8([(y2)/(12t4/3)]+[x/(t1/3)]), f(k) = [1/(2p)]ттR2dx dy  u0(x,y)f(x,y,k,0)exp(-i(k+[`k])x-(k2-[`k]2)y).

5   Возмущение солитона уравнений ДС-1

Солитонное решение уравнения было найдено в 1988 году Boiti, Leon, Martina и Pempinelli. Это решение осциллирует по времени и экспоненциально убывает в пространственных направлениях.

q(x,h,t;r) = rlmexp(it(l2+m2))
2cosh(mx)cosh(lh)
×
1
(1-s|r|2 ml
16
(1+tanh(lh)) (1+tanh(mx)))
,
(61)
g1|x®-Ґ є l2
2cosh2(lh)
,   g2|h®-Ґ є m2
2cosh2(mx)
,
где lm положительные постоянные, определяющиеся из граничных условий при h®-Ґ и x®-Ґ; r - произвольный комлексный параметр.

Построим асимптотическое решение возмущенного уравнения ДС-1:

it Q+ 1
2
(x2+h2)Q+(G1+G2)Q = eiF,
x G1 = - s
2
h |Q|2,    hG2 = - s
2
x |Q|2.
(62)
Здесь e - малый параметр, s = ±1 соответствует так называемым фокусирующему и дефокусирующему уравнению ДС-1.

Возмущение может появиться, например, если, в основном следуя Davey, Stewartson и Djorjevic, Redekopp, рассмотреть поверхностные волны над неровным дном. В этом случае возмущение будет иметь вид: F є AQ. Здесь A - вещественная постоянная, ее знак соответствует уменьшению или увеличению глубины в направлении x.

Будем искать асимптотическое решение в виде:

Q(x,h,t,e) = q(x,h,t;t)+eU(x,h,t,t),
G1(x,h,t,e) = g1(x,h,t,t)+eV1(x,h,t,t),
G2(x,h,t,e) = g2(x,h,t,t)+eV2(x,h,t,t),
(63)

Для построения решения понадобится решать линеаризованное на солитоне уравнение ДС-1:

it U+(x2+h2)U+(g1+g2)U+(V1+V2)q = iF
x V1 = - s
2
h(q _
U
 
+ _
q
 
U),   hV2 = -s
2
x(q _
U
 
+ _
q
 
U).
(64)
Для построения базисного набора решений линеаризованного уравнения понадобятся некоторые сведения из метода обратной задачи для ДС-1. С уравнением ДС-1 связана двумерная гиперболическая система Дирака (см.Нижник, Fokas и Ablowitz, Fokas и Santini):
ж
з
и
x
0
0
h
ц
ч
ш
y = - 1
2
ж
з
и
0
q
s _
q
 
0
ц
ч
ш
y.
(65)

Обозначим через y+ и y- матричные решения задачи Гурса (следуяFokas и Santini):

y+11|x®-Ґ = exp(ikh),
y+12|x®-Ґ = 0,
y+21|h®Ґ = 0,
y+22|h®-Ґ = exp(-ikx);
y-11|x®-Ґ = exp(ikh),
y-12|x®Ґ = 0,
y-21|h®-Ґ = 0,
y-22|h®-Ґ = exp(-ikx).
(66)

Пусть y+(j), j = 1,2, - столбцы матрицы y+. Введем две билинейные формы - аналоги прямого и обратного преобразований Фурье:

(c,m)f = у
х
Ґ

-Ґ 
у
х
Ґ

-Ґ 
dxdh(c1m1 s _
f
 
 + c2m2 f);
(67)
здесь ci и mi - элементы столбцов c и m.
бcmсs = у
х


R2 
dk dl (c1(l)m1(k)s _
s
 
(k,l)+c2(l)m2(k) s(k,l)),
(68)
где cj - элемент строки c.

Обозначим f(j), j = 1,2, строку сопряженную y(j) = [y-j1,y+j2] относительно второй билинейной формы.

Теорема 5 [О разложении] Пусть Q такая, чтоaQ О L1ЗC1 для |a| Ј 3, если функция f такая, что af(x,h) О L1ЗC1 для |a| Ј 4, тогда она может быть представлена в виде:

f = -1
p
бy(1)(x,h,l),f(1)(x,h,k)с[^f],
где
^
f
 
= 1
4p
(y+(1)(x,h,k),f(1)(x,h,l))f.

Теорема 6 [Эволюция коэффициентов] Пусть U - решение первого из линеаризованных уравнений ДС-1 - функция гладкая и интегрируемая по x и h, aU О L1ЗC1 и aF О L1ЗC1, для |a| Ј 4 и t О [0,T0]. Тогда

t ^
U
 
= i(k2+l2) ^
U
 
+ у
х
Ґ

-Ґ 
dkў ^
U
 
(k-kў,l,t)c(kў)+
у
х
Ґ

-Ґ 
dlў ^
U
 
(k,l-lў,t)k(lў)+ ^
F
 
.
(69)

Обозначим g(t) = 1-s[(ml)/4]|r(t)|2 и g0 = g(0). Конечную формулу для асимптотического решения возмущенного уравнения ДС-1 дает

Теорема 7 [О модуляции параметра] Если

g(t) = g0exp(2At),    Arg(r(t)) є const,
где g0 > 1 при s = -1 и 0 < g0 < 1 при s = 1, тогда асимптотическое решение (63) по mod(O(e2)) равномерно пригодно при t = O(e-1).

Следствие 2 [О больших временах] Если переменная t больше, чем e-1, а именно, t << e-1log(log(e-1)), тогда формулы (63) дают асимптотическое решение для (62) по mod(o(1)).

Предположение 1 [О сингулярности] Приведенные выше формулы пригодны для решений (61) без особенностей. Если s = 1, тогда [(ml)/4]|r(t)|2 < 1. Если в возмущении A > 0, тогда |r| возрастает и может появиться особенность в главном члене решения при t®Ґ. Вообще говоря, появление особенности в решении неинтегрируемого уравнения ДС - известное явление (см., например, Papanicolaou, C.Sulem, P.L.Sulem and Wang (1994)).

Список литературы

Гадыльшин Р.Р., Киселев О.М.. О бессолитонной структуре возмущенного солитонного решения уравнения Деви-Стюартсона. ТМФ, 1996, т.106, n2, с.167-173.

Киселев О.М.. Асимтпотика бессолитонного решения уравнения Деви-Стюартсона II. Дифференциальные уравнения, 1997, т.33, n6, c.812-819.

Gadyl'shin R.R., Kiselev O.M.. Asymptotics of perturbed soliton solution for the Davey-Stewartson II equation. http://xxx.lanl.gov/solv-int/9801014

Kiselev O.M.. Perturbation theory for the Dirac equation in two-dimensional space. J.Math.Phys., 1998, v.39, p.2333-2345.

Киселев О.М.. Асимптотика решения задачи Коши для уравнения Деви-Стюартсона I.Теор.Матем.Физ., 1998, т.114, n1, с.104-114.

Киселев О.М.. Базисные функции, связанные с двумерной системой Дирака. Функц. анализ и его приложения,1998, т.32, n1, с.56-59.

Гадыльшин Р.Р., Киселев О.М.. Структурная неустойчивость солитонауравнения Деви-Стюартсона II. ТМФ, 1999, т.118, n3, с.354-361. с.812-819.

Киселев О.М.. Асимптотика решения двумерной системы Дирака с быстро осциллирующими коэффициентами. Матем. сборник, 1999, т190, n2, с.71-92.

Kiselev O.M.. Dromion Perturbation for the Davey-Stewartson-1 Equations. Journal of Nonlinear Mathematical Physics, 2000, v.7, n4, p.411-422.

Kiselev O.M.. Asymptotic behaviour of a solution for the Kadomtsev-Petviashvili-2 equation. http://xxx.lanl.gov/math-ph/0003014


File translated from TEX by TTH, version 1.57.