Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Кудряшов Методы нелинейной математич 2008

.pdf
Скачиваний:
511
Добавлен:
16.08.2013
Размер:
1.5 Mб
Скачать

3.13. Применение алгоритма Конта—Форди—Пикеринга 171

u =

u0

+

u1

 

 

(3.13.10)

 

 

x2

 

x3

 

 

в уравнении (3.13.7), получим, что оно удовлетворяется.

 

Подставляя

 

 

v1

 

 

v(x) =

v0

+

 

(3.13.11)

x5

x4

 

 

 

 

 

в уравнение (3.13.8) и принимая во внимание (3.13.10), находим

 

 

 

 

v

=

3

 

u

u

 

 

v

=

 

1

u2.

(3.13.12)

 

 

 

 

 

 

 

 

 

 

 

1

2

 

0

 

1

 

0

2

 

0

 

Полагая

 

 

 

 

y(x) =

y0

+

y1

 

 

 

 

(3.13.13)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x6

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x7

 

 

 

 

в уравнении (3.13.9), имеем:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y1

=

5

 

u02u1

y0

=

 

1

u03.

(3.13.14)

 

 

 

 

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

4

 

 

 

 

Результаты вычислений позволяют написать разложение ре-

шения исходного уравнения в виде:

 

 

 

 

 

 

 

 

 

2

 

+ ε

u0

 

 

u1

+ ε2

 

 

u02

 

3u0u1

+

w(x) =

 

 

 

 

+

 

 

 

+

 

 

x

x3

x2

2x5

2x4

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

+ . . .

 

 

(3.13.15)

 

 

 

 

u03

 

5u02u1

 

 

 

 

 

 

 

+

 

 

 

 

 

 

 

 

4x7

4x6

 

 

 

Из (3.13.15) следует, что решение имеет две произвольные постоянные u0 и u1, соответствующие отрицательным индексам Фукса. Приближённое решение (3.13.15) может быть просуммировано [93], и в результате находится общее решение уравнения (3.12.1).

172 Глава 3. Аналитические свойства нелинейных уравнений

3.14.Преобразование Миуры и пара Лакса для уравнения Кортевега—де Вриза

После открытия солитона [155] многие исследователи предприняли попытку найти преобразование, связывающее решения уравнения Кортевега—де Вриза с решениями других точно решаемых уравнений. Это удалось сделать Р. Миуре, который установил, что преобразование

u = vx − v2

(3.14.1)

позволяет находить решения уравнения Кортевега—де Вриза

K [u] = ut + 6uux + uxxx = 0

(3.14.2)

при известном решении модифицированного уравнения Кортевега—де Вриза

M [v] = vt 6v2vx + vxxx = 0.

(3.14.3)

Миура нашел, что для уравнений (3.14.2) и (3.14.3) справедливо соотношение

K [u] =

2v M [v] .

(3.14.4)

∂ x

Однако само по себе преобразование (3.14.1) не упрощает решение уравнения Кортевега—де Вриза, поскольку связывает решения двух нелинейных уравнений. Тем не менее, именно это преобразование стало ключом к открытию метода обратной задачи рассеяния.

Заметим, что модифицированное уравнение Кортевега—

де Вриза (3.14.3) можно представить в виде

vt 6v2vx + vxxx = vt + ∂ ∂ + 2v vx − v2 = ∂ x ∂ x

∂ ∂

= vt + ∂ x ∂ x + 2v u = 0 ,

3.14. Преобразование Миуры и пара Лакса

173

и поэтому уравнения (3.14.1) и (3.14.3) можно записать как систему уравнений:

vx = u + v2,

(3.14.5)

∂ ∂

vt = ∂ x ∂ x + 2v u .

При этом условие совместности для уравнений (3.14.5)

(vx)t = (vt)x

(3.14.6)

дает уравнение Кортевега—де Вриза (3.14.2). Уравнение (3.14.2) допускает группу преобразований Галилея, поэтому в преобразование Миуры (3.14.1) можно добавить постоянное число λ. В этом случае система уравнений (3.14.5) примет вид

 

vx = u + λ + v2,

(3.14.7)

vt =

 

+ 2v (u − 2λ) .

(3.14.8)

∂ x

∂ x

Если ввести замену

 

 

 

v =

Ψx

(3.14.9)

 

 

 

,

 

 

Ψ

то систему уравнений (3.14.7) можно привести к системе линейных уравнений [155]

Ψxx + (u + λ) Ψ = 0 ,

(3.14.10)

Ψt = (c + ux) Ψ 2 (u − 2λ) Ψx,

где c (t) — находится при интегрировании (3.14.8) по x. Система уравнений (3.14.10) может быть представлена в виде

[155]

&

(3.14.11)

LΨ = λΨ;

 

174 Глава 3. Аналитические свойства нелинейных уравнений

 

Ψt = &

 

(3.14.12)

 

 

 

AΨ,

 

 

L

 

 

 

 

 

 

где & — оператор Шредингера

2

 

 

 

L =

 

u

(3.14.13)

∂ x2

&

 

 

&

и A — вспомогательный эволюционный оператор, имеющий вид

A = 2 (2λ − u)

+ (c + ux) .

(3.14.14)

 

 

∂ x

Используя условие&

совместности

 

 

 

xx)t = (Ψt)xx ,

(3.14.15)

находим, что уравнение Кортевега—де Вриза эквивалентно операторному уравнению [1, 39, 60]:

 

&

 

&

&

(3.14.16)

 

Lt +

L, A = 0,

где

L, A — коммутатор.

 

 

 

 

Если задано нелинейное уравнение в частных производных:

 

 

 

 

 

 

 

& &

 

 

 

 

 

ut = E (u,

ux,

. . . , x, t) ,

(3.14.17)

то можно дать определение.

Определение 3.14. Система уравнений (3.14.11),(3.14.12) называется парой Лакса для уравнения (3.14.17), если операторное уравнение (3.14.16) эквивалентно уравнению (3.14.17).

Уравнение Кортевега—де Вриза относится к классу нелинейных уравнений в частных производных, которые не могут быть сведены к одному линейному уравнению. Однако оно эквивалентно системе линейных уравнений относительно новой функции

Ψ (x, t, λ).

Используя систему уравнений (3.14.11) и (3.14.12), можно построить решение задачи Коши уравнения Кортевега—де Вриза

3.15. Законы сохранения

175

методом, получившим название метода обратной задачи рассеяния.

В этом методе при заданном начальном условии для уравнения (3.14.2)

u (x, t = 0) = f (x) ,

(3.14.18)

таком, что

(1 + |x|) |f (x) | dx < ∞,

(3.14.19)

−∞

находятся данные рассеяния путем решения стационарного уравнения Шредингера (3.14.11). Затем, из уравнения (3.14.12) определяется временная зависимость данных рассеяния. После чего решается обратная задача: определяется потенциал, который как раз и является решением уравнения Кортевега—де Вриза.

3.15.Законы сохранения для уравнения Кортевега—де Вриза

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

ρ

+

u) = 0,

(3.15.1)

∂ t

∂x

 

 

 

где ρ — плотность газа; u — его скорость; x — координата; t — время.

176 Глава 3. Аналитические свойства нелинейных уравнений

Интегрируя это уравнение по x получаем

 

 

 

d

 

ρ dx = ρ u|−∞.

(3.15.2)

dt

−∞

Предполагая, что потоки массы газа при x → ± ∞ равны нулю, из (3.15.2) находим, что интеграл от плотности:

ρ dx = const

(3.15.3)

−∞

не зависит от времени и является сохраняющейся величиной. Аналогичные соображения могут быть высказаны относительно любого уравнения математической физики, которое представляется в виде:

 

∂ T

+

∂ X

= 0.

(3.15.4)

 

∂ t

 

 

 

∂ x

 

Определение 3.15. Говорят, что уравнение в частных про-

изводных

 

 

 

 

E (u, ux, ut, . . . , x, t) = 0

(3.15.5)

имеет n законов сохранения, если оно может быть представлено в виде

∂ Ti

+

∂ Xi

= 0, i = 1, . . . , n,

(3.15.6)

∂ t

∂ x

 

 

 

где Ti — плотность; Xi — поток; причем Ti и Xi зависят от u и ее производных по x и t.

Открытие солитонов поставило вопрос о законах сохранения

для уравнения Кортевега—де Вриза:

 

 

ut + 6uux + uxxx = 0.

(3.15.7)

Очевидно, что это уравнение может быть представлено в виде

 

∂ u

 

 

 

 

+

 

uxx + 3u2 = 0,

(3.15.8)

 

∂ t

∂ x

3.15. Законы сохранения

177

и поэтому для этого уравнения можно взять

 

T1 = u,

X1 = uxx + 3u2,

(3.15.9)

что дает

 

 

u dx = const.

(3.15.10)

−∞

Были найдены также два других закона сохранения:

T2 = u2,

 

X2 = 2uuxx − ux2 + 4u3,

(3.15.11)

 

1

 

 

 

 

 

T3 = u3

 

ux2 ,

 

 

 

 

2

 

 

(3.15.12)

X3 = 3u2uxx +

9

u4 6uux2 + uxuxxx +

1

uxx2 .

 

 

2

2

 

Однако новые законы сохранения для уравнения Кортевега— де Вриза оказалось найти сложнее. К. Гарднер предложил для этой цели использовать преобразование Миуры в виде [1, 60]

 

u = w − ε wx ε2w2,

(3.15.13)

позволяющее получить соотношение между:

 

wt + 6

w − ε2w2 wx + wxxx = 0

(3.15.14)

и уравнением Кортевега—де Вриза (3.15.7). Это соотношение имеет вид

ut + 6uux + uxxx =

 

= 1 2w − ε x wt + 6

w − ε2w2 wx + wxxx .

(3.15.15)

178 Глава 3. Аналитические свойства нелинейных уравнений

Из (3.15.15) следует, что u (x, t) не содержит параметр ε. Это возможно, если w = w (x, t, ε). Будем искать w (x, t, ε) в виде разложения по степеням ε:

wn (x, t) εn.

 

w (x, t, ε) =

(3.15.16)

n=0

 

 

Из (3.15.14) получаем

 

 

w (x, t, ε) dx = const,

(3.15.17)

−∞

 

 

и поэтому

 

 

wn (x, t) dx = const,

n = 1, 2, . . .

(3.15.18)

−∞

 

 

Подставляя (3.15.16) в формулу (3.15.13) и приравнивая выражения при одинаковых степенях ε нулю, имеем

w0 = u, w1 = w0, x = ux,

w2 = w1, x + w02 = uxx + u2,

(3.15.19)

w3 = w2,x + 2w0w1 = uxxx + 4uux,

w4 = w3, x + 2w0w2 + w12 = uxxxx + 6uuxx + 5u2x + 2u3.

Все эти выражения дают законы сохранения для уравнения Кортевега—де Вриза. Процедура поиска wn может быть продолжена, что ведет к бесконечному множеству законов сохранения. Это важное свойство типично для точно решаемых уравнений методом обратной задачи рассеяния.

3.16. Преобразования Бэклунда

179

3.16.Преобразования Бэклунда для уравнения Кортевега—де Вриза

Если преобразование Миуры для уравнения Кортевега — де Вриза использовать в виде

u = vx λ − v2,

(3.16.1)

где λ — постоянная, то соотношение Миуры примет вид:

ut + 6uux + uxxx =

+ 2v

 

 

 

 

∂ x

 

vt 6 v2

+ λ

vx + vxxx .

(3.16.2) Соотношение (3.16.2), связывающее модифицированное уравнение Кортевега—де Вриза и уравнение Кортевега—де Вриза, после замены v → −v остается. Однако решение уравнения

Кортевега—де Вриза при этом будет другим

u = λ − vx − v2.

(3.16.3)

Складывая и вычитая (3.16.1), (3.16.3), имеем

u +2 u = λ + v2,

Обозначив

u = ωx,

из (3.16.4) находим

u − u

= vx.

(3.16.4)

2

 

 

u = ωx,

 

(3.16.5)

 

 

1

ω ω 2 = 0.

 

 

ωx + ωx + 2λ +

 

(3.16.6)

 

2

Модифицированное уравнение Кортевега—де Вриза

 

 

vt 6

v2 + λ vx + vxxx = 0,

 

с учетом (3.16.4) можно записать в виде

 

ωt ωt + 3

ωx + ωx ωx ωx + ωxxx ωxxx = 0.

(3.16.7)

180 Глава 3. Аналитические свойства нелинейных уравнений

При этом уравнение, которому удовлетворяют ω и ω , имеет вид

ωt + 3ω2

+ ωxxx = 0.

(3.16.8)

x

 

 

Система уравнений (3.16.6) и (3.16.7) является преобразованиями Бэклунда для уравнения (3.16.8), поскольку при заданном ω , решение ω, найденное из (3.16.6), (3.16.7), является также решением уравнения (3.16.8).

Преобразования Бэклунда (3.16.6), (3.16.7) могут быть использованы при нахождении аналитических решений уравнения Кортевега—де Вриза. Полагая ω = 0 в этой системе, из (3.16.6) имеем уравнение:

ωx 2k2 +

1

ω2 = 0,

λ = −k2,

(3.16.9)

2

с решением в виде

 

 

 

 

ω (x, t) = 2k th [kx + c (t)] .

(3.16.10)

Подставив (3.16.10) в (3.16.7), находим c (t) = 4k2t + ϕ0, где ϕ0 — постоянная. Решение (3.16.10) дает односолитонное решение уравнения Кортевега—де Вриза. Используя итерационный процесс из (3.16.6), (3.16.7), можно найти другие решения уравнения Кортевега—де Вриза.

Однако при нахождении решений удобно воспользоваться алгебраической формулой, которая может быть получена из (3.16.6). Полагая ω = ω1, ω = ω0, 2λ = λ1, из (3.16.6) имеем [39, 60]:

ω1x + ω0x = λ1

1

1 ω0)2 .

(3.16.11)

2

Аналогичным путем получаем соотношения:

 

ω2x + ω0x = λ2

1

2 ω0)2 ,

(3.16.12)

 

2

ω3x + ω1x = λ2

1

3 ω1)2 ,

(3.16.13)

 

2