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

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

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

2.6. Интегрирование дифференциальных уравнений

111

(2.6.3) будет также интеграл XF = C2. Однако поскольку уравнение (2.6.2) является уравнением 1-го порядка, то эти решения функционально зависимы, и поэтому

ξ

∂F

+ η

∂F

= f (F ) .

(2.6.5)

∂x

∂y

 

 

 

 

Учитывая (2.6.3) и (2.6.5) имеем

∂F

=

Qf

 

,

∂F

=

P f

 

.

 

 

 

 

 

 

 

 

∂x

ξQ

ηP

∂y

ξQ

ηP

 

 

 

 

 

 

 

 

 

 

Откуда приходим к выражению

Qdx − P dy

=

dF

,

(2.6.6)

ξQ − ηP

f (F )

 

 

 

являющемуся полным дифференциалом.

Соотношение (2.6.6) приводит к результату Софуса Ли: оператор (2.6.4) допускается уравнением (2.6.2) тогда и только тогда, когда функция

μ =

1

(2.6.7)

ξQ − ηP

является его интегрирующим множителем. Поэтому, зная какойлибо оператор, допускаемый уравнением 1-го порядка по формуле (2.6.7), можно найти интегрирующий множитель.

Используя этот подход, проинтегрируем уравнение Рикатти, имеющее вид [32]:

dy

+ y2 =

2

.

(2.6.8)

dx

 

 

x2

 

Подставляя x = ax, y = by, находим, что уравнение (2.6.8) будет инвариантным, если b = 1/a, и поэтому оно допускает группу преобразований растяжения

x = xeα, y = yeα.

(2.6.9)

112

Глава 2. Элементы группового анализа

 

Инфинитезимальный оператор, соответствующий преобразо-

ваниям (2.6.9) имеет вид

 

 

 

 

 

X = x

− y

(2.6.10)

 

 

 

.

 

∂x

∂y

Уравнение (2.6.8) можно представить в виде

2

 

dx = 0.

 

dy + y2

 

 

(2.6.11)

x2

Интегрирующий множитель для (2.6.11) находим по формуле

(2.6.7)

 

 

x

 

 

 

 

μ =

 

 

 

.

(2.6.12)

 

 

 

x2y2 − xy − 2

Он приводит к решению уравнения (2.6.11) в виде:

 

 

xy − 2

=

c

.

(2.6.13)

 

 

 

 

xy + 1

x3

 

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

хода основано на следующей теореме [30, 32].

Теорема 2.5. Всякая однопараметрическая локальная группа преобразований x i = f i (x, a) в RN с помощью некоторой невырожденной замены переменных zi = zi (x) может быть

приведена к группе переноса.

 

 

 

 

 

Доказательство.

Пусть

 

группе

преобразований

x = f i (x, a) (i = 1, . . . ,

N ) соответствует оператор

 

X = ξi (x)

 

.

(2.6.14)

 

∂xi

 

 

 

 

 

Этот оператор может быть представлен в виде

X˜ = X zi

 

 

,

(2.6.15)

 

 

 

∂zi

2.6. Интегрирование дифференциальных уравнений

113

если использовать соотношение

 

 

 

 

 

 

ξi

=

ξi

∂zj

 

.

(2.6.16)

∂xi

∂xi

∂zj

 

 

 

 

 

Выбирая в качестве новых переменных любой набор N − 1

функционально независимых инвариантов

 

Z1 = I1 (x) , Z2 = I2 (x) , . . . , ZN −1 = IN −1 (x) ,

 

а в качестве ZN переменную, которая находится из уравнения

 

X (z) = 1,

 

 

(2.6.17)

получим набор переменных, определяющих группу преобразований вдоль оси ZN .

Рассмотрим интегрирование уравнения (2.6.8) с помощью соответствующей замены переменных. Уравнение (2.6.8) допускает оператор (2.6.10), который в соответствии с вышедоказанной теоремой приводит к группе преобразований переноса после замены

t = t (x, y) ,

u = u (x, y) ,

(2.6.18)

которая находится из системы уравнений

 

x

∂u

− y

∂u

= 0,

(2.6.19)

 

 

 

 

∂x

∂y

x

∂t

− y

∂t

= 1.

(2.6.20)

 

 

 

∂x

∂y

В качестве одной из переменных можно взять инвариант

 

 

u = xy,

 

(2.6.21)

а в качестве второй переменной получаем

 

 

t = ln x.

 

(2.6.22)

114

Глава 2. Элементы группового анализа

Эта замена приводит группу преобразований растяжений к группе преобразований переноса вдоль оси x, поскольку

t = ln x = ln x + a = t + a ,

 

 

u

=

y

=

y

= u .

 

x 2

 

 

 

 

 

 

x2

 

Подставляя (2.6.21) в уравнение (2.6.8), приходим к уравне-

нию

 

 

 

 

 

 

 

du

+ u2 − u − 2 = 0,

(2.6.23)

 

dt

которое в результате интегрирования дает решение

 

u + 1

 

ln

 

 

= 3 t + c,

 

u − 2

 

имеющее в старых переменных вид

 

y =

2x3 + c

(2.6.24)

 

.

x (x3 − c)

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

2.7.Группы преобразований для линейного уравнения теплопроводности

Найдем группы преобразований, которые допускаются линейным уравнением теплопроводности

ut = a2uxx.

(2.7.1)

2.7. Группы для линейного уравнения теплопроводности 115

Здесь u — температура cреды, которая зависит от x и t, a2 — коэффициент температуропроводности.

Общий вид инфинитезимального оператора можно представить формулой

1

2

 

 

X = ξ

 

+ ξ

 

+ η

 

,

(2.7.2)

∂x

∂t

 

 

 

 

∂u

 

где ξ1, ξ2 и η — координаты касательного векторного поля, которые требуется найти. Оператор, продолженный на вторые производные, можно записать в виде:

X = X + ξ1

+ ξ2

+ ξ11

.

(2.7.3)

∂ux

 

 

2

 

∂ut

∂uxx

 

Действуя оператором (2.7.3) на уравнение (2.6.1), получаем соотношение

ξ2 − a2ξ11 = 0,

(2.7.4)

в котором производные ut и uxx связаны уравнением (2.6.1). Находим выражения для ξ1, ξ2 и ξ11 через ξ1, ξ2и η.

Поскольку

ξ1 = Dx (η) − uxDx ξ1 − utDx ξ2 ,

ξ2 = Dt (η) − uxDt ξ1 − utDt ξ2 ,

ξ11 = Dx 1) − uxxDx ξ1 − uxtDx ξ2 ,

то

ξ1 = ηx + uxηu − uxξ1x (ux)2 ξ1u − a2uxxξ2x − a2uxuxxξ2u, (2.7.5)

ξ2 = ηt + a2uxxηu − uxξ1t − a2uxuxxξ1u

(2.7.6)

−a2uxxξ2t − a4 (uxx)2 ξ2u ,

116

Глава 2. Элементы группового анализа

ξ11 = ξ1x + uxξ1u − uxxξ1x − uxuxxξ1u − a2uxxxξ2x − a2uxuxxxξ2u = = ηxx + uxxηu + uxηux − uxxξ1x − uxξ1xx 2uxuxxξ1u

2 (ux)2 ξux1 − a2uxxxξx2 − a2uxxξxx2 − a2 (uxx)2 ξu2

(2.7.7)

 

−a2uxuxxxξu2 − a2uxuxxξux2 + uxηxu + (ux)2 ηuu

 

(ux)2 ξxu1 (ux)3 ξuu1 − a2uxuxxξxu2 − a2 (ux)2 uxxξuu2

 

−uxxξx1 − uxuxxξu1 − a2uxxxξx2 − a2uxuxxxξu2 .

 

Подставляя (2.7.6) и (2.7.7) в определяющее уравнение (2.7.4)

и приводя подобные члены, получим уравнение:

 

 

 

 

ηt − a2ηxx + ux ξt1 + a2ξxx1 2a2ηux +

 

 

 

)2

 

 

2ξ1

2

η

 

 

3 2ξ1 +

 

 

 

 

+ (ux

 

 

 

 

3a

 

 

 

a

 

uu

+ (ux) a

 

uu

 

 

 

 

 

 

 

 

 

ux

 

 

 

 

 

 

 

 

 

 

 

 

 

a2ξ2

 

 

u

+ a4ξ2

 

+

 

 

 

 

(2.7.8)

 

+uxx

+ 2a2ξ1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

x

xx

 

 

 

 

 

 

 

 

+u

x

u

xx

3a2ξ1

 

a2ξ2

+ 2a4ξ2

 

+

 

 

 

 

 

 

 

 

 

 

 

 

 

u

 

 

 

 

ux

 

 

 

 

 

 

 

4ξ2

+ u

 

u

xxx

2a4

ξ2

2

 

 

a4

ξ2

= 0.

 

+uxxx2a

 

x

 

 

x

 

 

 

u + (ux) uxx

uu

 

 

 

Откуда находим уравнения для определения координат ξ1, ξ2

и η:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ξ2

= ξ2 = 0,

 

 

 

ξ1

= 0,

 

ξ1 = 0,

 

 

η

uu

= 0,

(2.7.9)

u

x

 

 

 

 

 

 

 

u

 

 

 

 

 

xx

 

 

 

 

 

 

 

ηt = a2ηxx,

 

 

ξt2 = 2ξx1 ,

ξt1 = 2a2ηux.

(2.7.10)

Из уравнений (2.7.9) находим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ξ2 = Ψ (t) ,

 

 

 

 

 

 

 

(2.7.11)

 

 

 

 

 

 

 

 

ξ1 = ϕ1 (t) + xϕ2 (t) ,

 

 

 

 

 

(2.7.12)

2.7. Группы для линейного уравнения теплопроводности 117

η = b (x, t) + ub1 (x, t) .

(2.7.13)

Подставив полученные соотношения во второе и третье урав-

нение (2.7.10), имеем

1

 

 

 

 

 

 

 

 

 

 

 

 

 

Ψt,

 

 

 

 

ϕ2 (t) =

 

 

 

 

 

 

 

(2.7.14)

2

 

 

 

 

 

 

 

x2

 

 

 

 

 

 

x2

 

b1 (x, t) = ϕ (t)

 

 

 

 

· ϕ1t

 

 

Ψtt.

2a2

8a2

Принимая во внимание эти выражения, из первого уравнения

(2.7.10) приходим к соотношениям

 

 

 

 

 

 

 

 

 

bt = a2bxx,

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

ϕt =

 

Ψtt,

 

 

 

(2.7.15)

4

 

 

 

ϕ1tt = 0,

 

 

 

Ψttt = 0.

 

Откуда находим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ψ (t) = B2 + 2B3t + B6t2,

 

ϕ1 (t) = B1 + B4t,

 

 

 

(2.7.16)

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ϕ (t) = B5

 

B6t.

 

 

 

 

2

 

 

 

 

Таким образом, координаты касательного векторного поля ξ1,

ξ2и η принимают вид

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ξ1 = B1 + B4t + (B3 + B6t) x,

 

ξ2 = B2 + 2B3t + B6t2,

 

 

 

(2.7.17)

1

 

 

 

 

 

 

B4x

B6x2

 

 

 

 

 

 

 

η = b (x, t) + B5u −

 

B6t +

 

 

+

 

 

u,

2

 

2a2

 

4a2

где B1, . . . , B6 — произвольные постоянные, b (x, t) — произвольная функция, которая удовлетворяет уравнению теплопроводности (2.7.1).

118

Глава 2.

Элементы группового анализа

 

Полагая Bi = 0

(i = 1, . . . , 6)

из (2.7.17) находим беско-

нечное число операторов

 

 

 

 

 

 

 

X

 

= b (x, t)

,

(2.7.18)

 

 

 

 

 

 

 

∂u

 

которые допускаются линейным уравнением теплопроводности (2.7.1). Из оператора (2.7.18) следует семейство преобразований

u (x, t) = a b (x, t) + u (x, t) , x = x, t = t ,

которое выражает свойство линейности для решений уравнения (2.7.1). Полагая b (x, t) = 0 и последовательно одну из постоянных Bi = 1, а остальные равными нулю, найдем еще шесть

операторов:

X1

=

,

 

 

X2 =

 

 

,

 

 

 

 

∂t

 

 

 

∂x

 

 

 

 

 

 

 

 

 

 

 

 

X4 = t

 

xu

 

 

 

 

 

 

 

 

 

 

 

,

 

 

 

∂x

2a2

∂u

 

 

 

 

 

 

 

 

 

X6 = xt

 

+ t2

 

 

 

 

∂x

∂t

 

 

X3

= x

 

+ 2t

,

 

∂x

∂t

 

 

 

 

 

 

 

 

 

 

X5 = u

,

 

(2.7.19)

 

 

 

 

 

 

 

 

 

 

 

∂u

 

 

t

 

x2

 

 

 

+

 

u

 

.

 

 

2

4a2

∂u

 

 

Эти операторы содержат важную информацию о симметрии уравнения теплопроводности.

Используя координаты касательных векторных полей в операторах (2.7.19), можно найти группы преобразований и инварианты, соответствующие этим операторам.

Например, решая уравнения Ли

dx = x , x (a = 0) = x , da

dt = 2t , t (a = 0) = t da

с координатами, соответствующими оператору X3, находим группу преобразований растяжения, которая допускается уравнением (2.6.1):

x = eax, t = e2at, u = u.

(2.7.20)

2.8. Группы для нелинейного уравнения теплопроводности 119

Откуда следуют два инварианта

I1

=

x 2

=

x2

, I2 = u

= u.

(2.7.21)

t

t

 

 

 

 

 

 

Эти выражения приводят к хорошо известному факту, что линейное уравнение теплопроводности (2.7.1) имеет инвариантное решение в виде

u = C1 w(z),

z =

C2 x

,

(2.7.22)

1

 

 

t 2

 

где C1 и C2 — произвольные постоянные.

Аналогично, используя операторы (2.7.19), можно найти и другие группы преобразований, которые соответствуют уравнению (2.7.1).

2.8.Группы преобразований для нелинейного уравнения теплопроводности

Для нелинейного уравнения теплопроводности

∂u

=

un

∂u

(2.8.1)

∂t

∂x

∂x

 

 

 

группа преобразований находится аналогично, как это было сделано для уравнения (2.7.1).

Пусть оператор X представляется в виде

X = ξ1

+ ξ2

+ η

,

(2.8.2)

∂x

∂t

∂u

 

 

 

 

 

где ξ1, ξ2и η — координаты касательного векторного поля, которые следует найти. Предполагая, что

X = X + ξ1

+ ξ2

+ ξ11

,

(2.8.3)

∂ux

 

 

2

 

∂ut

∂uxx

 

120

Глава 2. Элементы группового анализа

и действуя этим оператором на уравнение (2.8.1), получаем соотношение

ξ2 − n (n − 1) (ux)2 un−2η 2nun−1uxξ1

(2.8.4)

−nun−1uxxη − unξ11 = 0.

Вычисляя по формулам продолжения компоненты ξ1, ξ2и ξ11 и подставляя выражения для них в определяющее уравнение (2.8.4), получаем четыре независимых оператора

X1 =

, X2 =

 

,

X3 = 2t

+ x

,

∂t

 

 

 

∂t

∂x

 

 

 

∂x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(2.8.5)

 

 

X4 =

nx ∂

+ u

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

2

∂x

∂u

 

 

 

 

В случае n = 4/3 уравнение (2.8.1) допускает дополнительный к (2.8.5) оператор в виде

2

(2.8.6)

 

 

 

 

 

X5 = x ∂x + 3xu

∂u .

 

Уравнение (2.8.1) при n = 0 переходит в линейное уравнение (2.7.1), имеющее группу преобразований, представленную в предыдущем параграфе.

Заметим, что поскольку операторы группы преобразований линейные, то, кроме “основных” операторов (2.8.5), уравнением (2.8.1) допускается и произвольная сумма их, умноженная на постоянные коэффициенты. Например, уравнение (2.8.1) допускает оператор

X =

1

,

(2.8.7)

c0

 

∂t

∂x

который ведет к инвариантам

ξ = x − c0t, u (x, t) = v,