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

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

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

3.20. Метод Вайса—Табора—Карневейля

191

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

3.20.Метод Вайса—Табора—Карневейля для анализа нелинейных уравнений

Применение гипотезы Абловица—Рамани—Сигура для анализа нелинейных уравнений в частных производных является трудоемкой процедурой. Это связано с тем, что, во-первых, требуется найти все преобразования, с помощью которых решения исходного нелинейного уравнения в частных производных можно найти через решения обыкновенного дифференциального уравнения. Для этой цели обычно используется групповой анализ дифференциальных уравнений, открытый Софусом Ли. Во-вторых, требуется исследовать на свойство Пенлеве все обыкновенные дифференциальные уравнения, полученные из исходного уравнения в частных производных, которых может оказаться достаточно много.

Этот недостаток был устранен в обобщении теста Пенлеве, предложенного в 1983 году в работе Дж. Вайса, М. Табора и Г. Карневейля [150, 151, 152].

Идея метода Вайса—Табора—Карневейля состоит в следующем. Пусть задано нелинейное уравнение в частных производных:

E (ux, ut,

. . . , x, t) = 0.

(3.20.1)

Его решение ищется в виде разложения

 

 

 

u (x, t) =

uj zj−p,

(3.20.2)

j=0

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

где z (x, t) — новая функция, p — число, которое определяется при подстановке

u = u0z−p

(3.20.3)

введущие члены уравнения (3.20.1). Из сравнения выражений с наименьшими степенями находятся семейства решений с конкретными p и u0. При этом коэффициенты uj в (3.20.2) зависят от производных функции z (x, t).

Метод, предложенный Вайсом, Табором и Карневейлем, во многом напоминает алгоритм Ковалевской, но отличается тем, что вместо переменной ξ = z − z0 используется функция z (x, t) и вместо постоянного коэффициента a0 предполагается, что коэффициент uj зависит от производных функции z (x, t). Новая функция z (x, t) является сингулярным многообразием, и поэтому метод Вайса—Табора—Карневейля часто называется методом сингулярных многообразий. Как и в случае алгоритма Ковалевской, на первом этапе в методе Вайса—Табора—Карневейля предполагается, что уравнение (3.20.1) проходит тест Пенлеве, если p

в(3.20.2) — целое.

На втором этапе ищутся индексы Фукса точно так же, как для обыкновенных дифференциальных уравнений. С этой целью в ведущие члены уравнения (3.20.1) подставляются выражение

u = u0z−p + ur zr−p.

(3.20.4)

Третий этап исследования состоит в определении произвольных функций для коэффициентов разложения (3.20.2) с номерами, равными индексам Фукса.

Как и в случае обыкновенного дифференциального уравнения, считается, что (3.20.1) проходит тест Пенлеве, если индексы Фукса — целые и число произвольных функций совпадает с порядком исходного уравнения. Проверка исходного уравнения на свойство Пенлеве, как и в случае обыкновенного дифференциального уравнения, позволяет сделать вывод лишь о необходимом условии интегрируемости уравнения.

3.21. Пенлеве-анализ уравнения Бюргерса

193

Однако Вайс, Табор и Корневейль пошли дальше и для целого ряда нелинейных уравнений получили достаточное условие интегрируемости. С этой целью они предположили, что uj 0 при j > p в (3.20.2) и в результате пришли к преобразованию:

u =

u0

+

u1

+ . . . + up,

(3.20.5)

zp

zp−1

 

 

 

 

позволяющему найти ряд свойств, характерных для точно решаемых нелинейных уравнений в частных производных.

3.21.Пенлеве-анализ уравнения Бюргерса методом Вайса—Табора—Карневейля

Проанализируем уравнение Бюргерса:

ut + uux νuxx = 0

(3.21.1)

методом Вайса—Табора—Карневейля [150]. Ведущими членами этого уравнения являются второе и третье слагаемые. Поэтому, подставляя

 

 

u =

u0

, ux

=

u0,x

pu0zx

,

 

 

 

zp

 

zp

zp+1

 

 

u

0, xx

 

 

2pu

z

 

 

pu z

 

 

u z2

uxx =

 

 

 

 

0,x x

0

xx

+ p (p + 1)

0 x

 

zp

 

 

zp+1

 

zp+1

zp+2

в уравнение (3.21.1) и приравнивая выражения нулю при наименьшей степени z (x, t), получаем соотношение:

 

pu2zx

 

νp (p + 1) u0z2

 

 

0

 

x

= 0.

(3.21.2)

 

z2p+1

zp+2

 

Из равенства (3.21.2) следует, что если u0

= 0 и zx = 0, то,

во-первых, 2p + 1 = p + 2, откуда p = 1 и, во-вторых,

 

 

u0 = zx.

 

 

(3.21.3)

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

Таким образом, для уравнения Бюргерса имеем одно семейство решений:

p = 1, u0 = zx.

(3.21.4)

Поскольку p = 1, то выполняется первое необходимое условие для того, чтобы уравнение Бюргерса проходило тест Пенлеве.

На втором этапе в ведущие члены уравнения (3.21.1) подставляется соотношение:

u =

zx

+ uj zj−1.

(3.21.5)

z

Поскольку

 

zxx

 

z2

+ uj (j − 1) zj−2zx + . . . ,

ux =

 

 

+

 

x

z

 

z2

 

 

 

z

 

 

 

z z

 

3 z

 

uxx =

 

xxx

+

x xx

x

+

 

z

 

 

z2

z3

+uj (j − 1) (j − 2) zx2zj−3 + uj (j − 1) zj−2zxx + . . . ,

то после подстановки u, ux и uxx в ведущие члены уравнения (3.21.1) и приравнивания выражений при первой степени uj , имеем соотношение:

u

νz2 zj−1

[2

2 (j

1)

(j

1) (j

2)] = 0,

(3.21.6)

j

x

 

 

 

 

 

 

откуда приходим к уравнению для индексов Фукса:

j2 − j − 2 = 0,

(3.21.7)

имеющему корни j1 = 1, j2 = 2.

Таким образом второе необходимое условие теста Пенлеве для уравнения Бюргерса также выполняется.

3.21. Пенлеве-анализ уравнения Бюргерса

195

Индекс j1 = 1 соответствует произвольно выбранной функции z (x, t), а коэффициент u2 в (3.20.2) должен быть произвольным.

Подстановка

u =

zx

+ u1 + u2z

(3.21.8)

z

в уравнение (3.21.1) и приравнивание выражений при одинаковых степенях z (x, t) дает следующую цепочку уравнений:

 

zt + u1zx νzxx = 0,

(3.21.9)

 

 

 

 

(zt + u1zx νzxx) = 0,

(3.21.10)

 

∂x

u1,t + u1u1,x νu1xx − u2zxx 2u2,xzx = 0 .

(3.21.11)

Очевидно, что если (3.21.9) удовлетворяется, то (3.21.10) переходит в тождество. Коэффициент u2 из (3.21.10) не определяется, и поэтому его можно взять произвольным. Таким образом, уравнение Бюргерса проходит тест Пенлеве. Если коэффициент u2 (x, t) взять равным нулю, то уравнение (3.21.11) примет вид

u1,t + u1u1,x νu1,xx = 0.

(3.21.12)

Пусть u1 является решением уравнения (3.21.12). В этом случае из решения линейного уравнения (3.21.9) можно найти z (x, t). Тогда формула (3.21.8) при u2 = 0 примет вид

u =

 

∂x ln z + u1.

(3.21.13)

Эта формула снова дает решение уравнения Бюргерса.

В частности, поскольку u1 = 0 является тривиальным решением уравнения (3.21.1), то, принимая во внимание (3.21.13), имеем хорошо известное преобразование Коула—Хопфа для решений уравнения Бюргерса [91, 101]:

u =

 

∂x ln z.

(3.21.14)

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

Это преобразование позволяет по решению уравнения теплопроводности

zt ν zxx = 0

(3.21.15)

находить решения уравнения (3.21.1). Полагая u1 = z в (3.21.12), находим, что если z (x, t) удовлетворяет уравнению Бюргерса, то формула (3.21.13) запишется в виде

u =

 

∂x ln z + z.

(3.21.16)

При известном решении уравнения Бюргерса эта формула позволяет находить другие решения этого уравнения.

Таким образом, анализ уравнения Бюргерса на свойство Пенлеве методом Вайса—Табора—Карневейля позволил, во-первых, установить, что это уравнение проходит тест Пенлеве, и, вовторых, найти формулу Коула—Хопфа, с помощью которой исходное уравнение может быть сведено к линейному уравнению теплопроводности.

3.22. Решение задачи Коши для уравнения Бюргерса

Покажем, что уравнение Бюргерса

 

ut + uux = μ uxx.

(3.22.1)

с помощью преобразования [91, 101]

 

u =

ln z

(3.22.2)

 

 

∂ x

можно свести к линейному уравнению.

Действительно, учитывая (3.22.2), имеем соотношение:

1

 

ut + uux μ uxx = x z (zt μ zxx) ,

(3.22.3)

3.22. Решение задачи Коши для уравнения Бюргерса

197

из которого следует, что решение уравнения теплопроводности

zt μ zxx = 0

(3.22.4)

по формуле (3.22.2) преобразуется в решение уравнения (3.22.1). Преобразование Коула—Хопфа (3.22.2) позволяет найти ре-

шение задачи Коши для уравнения Бюргерса. Пусть в начальный момент имеем

 

 

 

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

 

(3.22.5)

тогда из (3.22.2) получаем

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

0

 

z (x, t = 0) = Ψ (x) = exp

1

ϕ (ξ) dξ .

(3.22.6)

 

 

 

 

 

 

 

 

 

Решение задачи Коши для линейного уравнения теплопровод-

ности находится с помощью функции Грина [4, 64]:

 

z (x, t) =

1

Ψ (θ) exp

,

(x − θ)2

dθ.

(3.22.7)

 

 

 

 

2πμ t

 

t

-

−∞

Используя (3.22.7), получаем решение задачи Коши для уравнения Бюргерса в виде

 

.

θ)

 

 

 

G,x,t)

 

 

 

(x

 

 

 

 

u (x, t) =

 

−∞

t

 

exp #

 

$ dθ

,

(3.22.8)

.

 

 

 

 

 

 

 

 

 

 

 

 

G,x,t)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

−∞ exp #

 

 

$ dθ

 

 

 

 

где

 

 

 

 

 

 

 

 

 

 

 

 

 

G , x, t) =

(x − θ)2

+

θ

ϕ (ξ) dξ.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(3.22.9)

 

 

 

 

 

2t

 

 

 

 

 

0

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

Рассмотрим решение (3.22.8) при фиксированных x и t и при μ 0. Основной вклад в интегралы в (3.22.8) дают окрестности

стационарных точек функции G , x, t), таких, что

 

 

∂ G

= ϕ (θ)

x − θ

= 0.

(3.22.10)

 

θ

 

 

t

 

Пусть θ = θ (x, t) — решение этого уравнения, тогда интеграл вблизи θ = θ вычисляется приближенно методом перевала [4] в соответствии с формулой:

g (θ) exp G (

2,μ

 

dθ ≈ g (θ ) G (θ )

1

G.

 

 

exp

 

 

 

 

 

 

θ x, t)

 

 

 

 

 

 

4πμ

2

 

(θ )

−∞

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В этом случае из (3.22.8) имеем

 

 

 

 

 

 

 

 

 

 

 

u (x, t)

x − θ ,

 

(3.22.11)

 

 

 

 

 

 

 

 

 

 

 

t

 

 

где θ — корень уравнения

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ϕ (θ )

 

x − θ = 0.

 

(3.22.12)

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

Таким образом, решение (3.22.8) при μ 0 можно предста-

вить в виде

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u (x, t)

ϕ (θ ) ,

x = θ + ϕ (θ ) t.

 

(3.22.13)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Оно совпадает с решением задачи Коши для простейшего ква-

зилинейного уравнения

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∂ u

+ u

∂ u

= 0,

 

 

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

 

(3.22.14)

 

 

 

 

 

 

 

 

 

∂ t

∂ x

 

 

 

 

 

 

 

 

 

 

 

 

которое получено в п. 1.10 методом характеристик.

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

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

199

3.23.Анализ уравнения Кортевега—де Вриза методом Вайса—Табора—Карневейля

Проверим уравнение Кортевега—де Вриза

 

ut + 6uux + uxxx = 0

(3.23.1)

на свойство Пенлеве, используя метод Вайса—Табора— Карневейля [150].

Ведущими членами уравнения являются второе и третье слагаемые. Подставляя (3.20.3) в эти слагаемые и приравняв нулю выражения при наименьшей степени z (x, t), находим соотноше-

ние

6 p u2zx

 

p (p + 1) (p + 2) u0z3

 

 

 

 

 

 

0

=

x

,

(3.23.2)

z2p+1

zp+3

из которого следует, что p = 2, u0 = 2zx2 . Таким образом, уравнение Кортевега—де Вриза имеет одно семейство решений, у которого

p = 2, u0 = 2zx2.

(3.23.3)

Из (3.23.3) следует, что первое необходимое условие теста Пе-

нлеве для уравнения (3.21.16) выполняется.

 

Подставляя выражение

 

u = 2 zx2 z2 + uj zj−2

(3.23.4)

снова в ведущие члены уравнения и приравнивая нулю слагаемые при первой степени uj , имеем соотношение:

u

z3 zj−5

[(j

2) (j

3) (j

4)

12 (j

2) + 24] = 0,

(3.23.5)

j

x

 

 

 

 

 

 

из которого следует уравнение для индексов Фукса:

j3 9 j2 + 14 j + 24 = 0.

(3.23.6)

Это уравнение имеет корни j1 = 1, j2 = 4, j3 = 6, откуда следует, что уравнение Кортевега—де Вриза на втором этапе также проходит тест Пенлеве.

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

На третьем этапе исследования проверяется произвольность коэффициентов u4 и u6. С этой целью в уравнение (3.23.1) подставляется выражение

u = 2 zx2 z2 + u1 z1 + u2 + u3z + u4z2 + u5z3 + u6z4.

(3.23.7)

Приравнивая

нулю

выражения

при

степенях

z4, z3, . . . , z1 и z0, соответственно, получаем

 

 

 

 

 

u1 = 2zxx,

 

(3.23.8)

 

 

ztzx + 4zxzxxx 3zxx2 + 6u2zx2

= 0,

(3.23.9)

 

zxt + 6u2zxx + zxxxx 6u3zx2 = 0,

 

(3.23.10)

 

 

 

 

 

 

 

zxt

+ zxxxx + 6u2zxx 6u3zx2 = 0.

(3.23.11)

 

∂x

Из (3.23.11) следует, что при выполнении (3.23.10) уравнение (3.23.11) также удовлетворяется, и, следовательно, коэффициент u4 может быть взят произвольным. Продолжая вычисления дальше, устанавливаем, что коэффициент u6 можно также взять произвольным.

Полагая uj = 0 при j 2 , из (3.23.7) имеем

 

 

2z2

2zxx

 

2

 

u =

x

+

 

+ u2 = 2

 

ln z + u2,

(3.23.12)

z2

z

∂x2

где u2 удовлетворяет опять же уравнению Кортевега—де Вриза

u2t + 6u2u2, x + u2, xxx = 0.

(3.23.13)

Таким образом, применение теста Пенлеве для уравнения Кортевега—де Вриза показало, что необходимое условие интегрируемости для него выполняется. В работах Дж. Вайса [150, 151, 152, 153] проиллюстрировано, что уравнение sin-Гордона, нелинейное уравнение Шредингера и ряд других нелинейных уравнений в частных производных также проходят тест Пенлеве.