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

книги из ГПНТБ / Теория ветвления и нелинейные задачи на собственные значения

..pdf
Скачиваний:
8
Добавлен:
24.10.2023
Размер:
9.31 Mб
Скачать

230

Г. Ф. ВАЙНБЕРГЕР

Мы предполагаем, что цилиндры бесконечные. Граничные условия прилипания требуют, чтобы возмущения к и п удовлетворяли равен­ ствам

 

 

и = ди/дг = и — 0 при г = г),

1.

(2.5)

(Условие ди/дг = 0 при г = т), 1

вытекает из

(2.3)

и того факта,

что w =

0 при г =

Г], 1.)

отыскании

функций и (г, z, t),

Наша

краевая

задача состоит в

V (г, z, t), удовлетворяющих уравнениям (2.2Ь), (2.4) и (2.5). Будем разыскивать такие решения в пространстве Ь2-

Применив к этим уравнениям преобразование Фурье по z, получим

( L - Я 2— ^-)

(L — А.2) й - Ѵ Т ( ^ — К ) V= 0,

(2.6)

ü + ( L - X a- 4

- ) o =

0,

 

(2.7)

где использованы обозначения

 

 

 

 

 

 

 

 

 

с о

 

 

и (г,

t, X) =

j eiXzu (г, z, t) dz,

 

 

 

 

 

 

—oo

 

 

 

 

 

с о

 

 

 

 

о (г,

t,X )=

j

eiXzü(r,

z, t)dz,

 

 

 

 

— OO

 

 

 

т _ д *

 

[

д

1

 

( 2. 8)

L — Qr% +

 

r

dr

А

'

 

 

T = — 4AB,

 

К =

— A/B.

(2.9)

Сайндж [1] показал, что рассматриваемое течение всегда устойчиво

при

(J, >

г]2. Поэтому *) будем рассматривать

случай,

когда

р <С ц2.

Из

(1.4)

следует, что А < 0, а поэтому Т >

0. Кроме того,

/С < 1,

так

что

1/г2 — К > 0. Мы интересуемся, существуют

ли при этих

обстоятельствах начальные данные, вызывающие неустойчивость. Возьмем теперь преобразование Лапласа от (2.6) и (2.7) по t. Полу­

чаем

 

 

(L - X2- s ) ( L - X 2) t / - r r ( - ^ - - K ) V

= /,

(2.10)

U + (L—Я,2 —s) V = J.

 

(2.11)

Здесь s — переменная преобразования Лапласа,

U и

V — преоб­

разования от и и о, а / и J — заданные начальные значения.

Мы ищем условия, при которых решения уравнений (2.6) и (2.7) растут, т. е. условия, при которых собственные значения ^несамосопряженной системы (2.10) и (2.11) имеют положительные действитель­

*) Нетрудно показать, что и при р = ту* имеет место асимптотическая устой­ чивость.— Прим. ред.

XIV. И З М Е Н Е Н И Е УСТОЙЧИВОСТИ В Т Е Ч Е Н И И КУЭТТА

231

ные части. Для достаточно малых Т, как известно, эти собственные значения имеют отрицательные действительные части. Поэтому неустой­ чивость возникает, когда число Т ровно настолько велико, чтобы s было чисто мнимым. Однако экспериментально установлено, что воз­ никновению неустойчивости соответствует s = 0. Это побуждает нас сформулировать следующий

П р и н ц и п и з м е н е н и я у с т о й ч и в о с т и д л я т е ч е н и я К у этт а. Если действительный положительный луч содержится в резольвентном множестве оператора, соответствующего однородной

Собственные значения

Рис. 2. Плоскость s.

1

системе (2.10), (2.11), то и правая полуплоскость также содержится в резольвентном множестве (см. рис. 2).

Остальная часть этой лекции посвящена развитию методов мате­ матического обоснования этого принципа.

3. Неотрицательные операторы

Пусть А — матрица, все элементы которой неотрицательны. Назо­ вем такую матрицу неотрицательной. Резольвентный оператор мат­ рицы А допускает формальное разложение

{А — s l y 1= — s - ^ / Als)~1 = —s-1(I+A/s + Äl/s‘i + • • •)• (3.1)

Последнее выражение в (3.1) есть ряд по степеням 1/s, имеющий неотрицательные коэффициенты. Привлечем теперь следующую тео­ рему Принсгейма: если степенной ряд с неотрицательными коэффи­ циентами имеет радиус сходимости R, то R есть особая точка функции, представленной этим рядом. Отсюда следует, что имеется действитель­ ное положительное собственное значение st матрицы А, такое, что множество

(3.2)

содержится в резольвентном множестве матрицы А, т. е. спектр А обя­ зан лежать в круге | s | ^ si плоскости s (см. рис. 3). Этот результат представляет собой часть теоремы Перрона — Фробениуса о неот­ рицательных матрицах (см. [2]).

232

Г. Ф. В А Й Н Б Е Р Г Е Р

Аналогично, если матрица В такова, что —Л-1 неотрицательна, то

(В — si)-1 = в -1 II - (— s) ( - В - 1)] = В -1 II + (— s) (—В -1) +

+ (- s)2 ( - Л - 1)2 + . . . ] . (3.3)

Так как степенной ряд в (3.3) имеет неотрицательные коэффициен­ ты, существует действительное собственное значение st < 0, такое, что круг

{s: I s I < —Sj}

(3.4)

содержится в резольвентном множестве матрицы В (см. рис. 4). Пусть теперь а — некоторое заданное действительное число, а С

 

Рис. 3. Спектр лежит

Рис.

4. Спектр лежит

 

 

в

заштрихованной

в заштрихованной

 

 

 

 

 

области.

 

области.

 

 

 

матрица, такая,

что —(С — s0/) _1 неотрицательна для

всех действи­

тельных s0 > а.

Тогда

 

 

 

 

 

(С -

si)-1 =

— s0/)e-1 {/ -

(s0-

s) [ -(С -

So/)-1]}'1. (3.5)

При помощи тех же рассуждений, что и ранее,

мы убеждаемся

в существовании

действительного

собственного значения

Si <С s0,

такого, что

круг

 

 

 

 

 

 

 

 

 

 

{s: I s0 — s I <

s0 — s j

 

 

(3.6)

содержится в резольвентном множестве матрицы В (см.

рис.

5).

Так как этот вывод остается в силе для всех s0 >

а, мы можем

сдвигать s0 вправо, оставляя собственное значение Si фиксированным. При So — оо отсюда следует, что полуплоскость

{s:Res>Sj}

(3.7)

содержится в резольвентном множестве матрицы С.

Естественным обобщением матричного оператора с неотрицатель­ ными элементами является интегральный оператор с неотрицатель­ ным ядром. Ентч [3] и Карлин [4] применяли теорему Принсгейма описанным выше способом для доказательства того факта, что теорема Перрона — Фробениуса остается в силе для таких интегральных опера­

XIV.

И З М Е Н Е Н И Е УСТОЙЧИВОСТИ

В Т Е Ч Е Н И И

КУЭТТА

23 3

торов. (По поводу

абстрактной

трактовки таких операторов см.

[4J

и [5] х).)

(3.7)

представляет

собой

утверждение

как раз такого

Результат

рода, какое нужно нам, чтобы установить принцип изменения устой­

чивости. Таким образом, если мы сумеем представить

(2.10) и (2.11)

в форме

 

si] и = /,

(3.8)

где С имеет описанные выше свойства, наша задача будет решена. Если бы С был дифференциальным оператором, подчиненным принци­ пу максимума, можно было бы доказать, что он обладает свойствами, нужными для вывода (3.7). К сожалению, данная задача не попадает в этот класс. Поэтому нам требуются более тонкие результаты.

4. Специальные результаты

 

Вернемся к уравнениям (2.10) и (2.11). Положим

 

r(s) =

[L — a*/ — si]-1,

(4.1)

G (s) =

[(L — %2I) (L X21 si)]-1.

(4.2)

Однородную задачу, соответствующую системе (2.10) и (2.11),

можно записать в виде

U К (s) U = 0,

 

(4.3)

где

 

 

 

 

 

MTG (s) [(Mr2 -

 

Г (s)];

(4.4)

 

 

К (S) =

К )

К (s) зависит от $

нелинейно,

но аналитически.

 

Рассмотрим

резольвенту

оператора

К (s):

 

[ I - К (s)]-1 =

{

/ - [ /

К (So)]"1 IK (s) -

К (so)]}-1 II — К (So)]"1-

 

 

 

 

 

 

(4.5)

Если теперь для всех действительных s0, больших некоторого а, (і) оператор [I К (s0)]-1 неотрицателен,

х) Укажем здесь еще на основополагающую в данном круге вопросов статьюКрейна и Рутмана [6]; см. также [7], [8].— Прим. ред.

234

Г. Ф. В А Й Н Б Е РГ Е Р

 

 

(іі) К (s) в окрестности s0 разлагается в ряд по степеням

(s0 — s)

с неотрицательными коэффициентами, т. е. для всех k

 

 

( — зг ) * * & ) > °»

<4-6)

то правая часть равенства (4.5) разлагается в ряд по степеням (s0 — s) с неотрицательными коэффициентами, и, применяя методы § 3, мы при­ ходим к выводу о существований действительного собственного значе­ ния Sj, такого, что спектр оператора К (s) содержится в множестве {s: Re s ^ Sj}. Этот результат эквивалентен принципу изменения устойчивости.

Справедливость условия (і) для нашего оператора устанавливается

немедленно. Чтобы проверить условие (іі),

заметим сначала, что

Г (s) есть интегральный оператор с ядром Г (г,

р; s), представляющим

■собой] преобразование Лапласа функции Грина Г* (г, р, t) следующей

.краевой задачи:

и (г), /) = и (1, f) —и {г, 0) = 0.

Легко видеть, что для этой задачи выполняется принцип максимума, так что / ^ 0 влечет за собой и 0. Следовательно, Г* (г, р, t) ^ 0. Так как

о о

Г (г, р; s) = j

e~str* (г, р, t) dt,

о

 

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

 

 

о о

( ~ ~ і г ) кт (г' р; s) =

J the-stT*(r, р; t)dt> 0

для всех k и всех s > —

о

 

Теперь мы должны рассмотреть оператор G (s). Его можно пред­ ставить в виде интегрального оператора, ядро которого G (k , р; s)

есть функция Грина задачи (L — Я2) (L X2 — s) и = /

с соответст­

вующими краевыми условиями. Далее, имеем

 

ОО

 

G (г, р; s)= j e~stG* (г, р; t) dt,

(4.7)

о

 

где G* — функция Грина нестационарной задачи

 

(.L - X 2) ( L - X 2— Ir )u = f ,

 

и =

= 0

при г — г] и г — 1,

(4.8)

 

и = 0

при ^ = 0.

 

XIV . И З М Е Н Е Н И Е УСТОЙЧИВОСТИ В Т Е Ч Е Н И И 1КУЭТТА

235

Как и выше, неравенства (—d/ds)hG (s) > 0 будут доказаны,

если

будет установлено, что функция Грина G* задачи (4.8) положительна. Если это уже проделано, то при помощи правила дифференцирования произведения легко убедиться в том, что оператор К (s) = К2TG (s) X

X ^ (-^------К) Г (s)] обладает свойством (іі).

Функцию Грина G (г, р; s) можно выписать явно с помощью бессе­ левых функций. Это выражение, однако, слишком сложно, чтобы можно было проверить условие (іі) непосредственно. Тем не менее, используя асимптотические разложения при s-> + °о , мы можем доказать существование такого > 0, что неравенство G* > 0 выпол­ няется для t <H\ .

Функцию G* (г, р, t) можно разложить в ряд по собственным

функциям задачи

(L -

X2) (L + у2) w =

О,

 

 

w =

xso' = 0 при г =

г] и г = 1.

Это разложение имеет вид

 

 

 

° °

2

 

G* = 2 wn (г) Wn (р) е ~ (Уп+%Ці.

 

 

1

 

Функции wn опять-таки выражаются через бесселевы функции;

в частности, можно доказать, что wt >

0. Отсюда следует существова­

ние такого tг, что G* >

0 для t > t2-

 

Вычисление показывает, что для каждого фиксированного X суще­

ствует т]0) достаточно

близкое к 1 и такое, что t2< ^ при rj > т)0-

Таким образом, G* > 0 для всех t, если зазор достаточно мал х).

(Однако при т]

1

критическое волновое число X растет до беско­

нечности, и нужны дополнительные вычисления для доказательства того факта, что условие G* > 0 выполняется равномерно по X для достаточно больших ту)

Можно доказать, что G* при т] -э- 0 стремится к функции Грина краевой задачи

( L - X 2)

и(1, 0 = -§г(1. t) = u(r, 0 ) = 0

для всего цилиндра 0 ^ 1. Можно доказать, что для этой задачи G* > 0 . Так как критическое'волновое число остается здесь конечным, можно утверждать, что принцип изменения устойчивости выполняется, когда радиус внутреннего цилиндра достаточно мал а).

і) При малых l - і і и 1 — р. принцип изменения устойчивости^ следует

из теории возмущений, так как в предельном случае 1 — И = 1 — И' = 0 полу­ чается (самосопряженная) задача о конвекции в слое жидкости 19]. Оценке допустимых значений 1 — т) и 1 — (А посвящена работа Дэвиса [11]. Прим.

^ 2) Устойчивость в этом предельном случае исследована в работе Овчинни­ ковой [10].— Прим. ред.

236

Г. Ф. В А Й Н Б Е Р Г Е Р

5.

Заключение

Для интегрального оператора К (s) с ядром К (s; х, у)

условия (і)

и (іі)

(§ 4) можно заменить следующими более слабыми

условиями:

(i)

оператор [/ — К (s)]“1 неотрицателен для всех достаточно боль­

ших действительных

s;

 

 

(ii) К (s) и —

(s) неотрицательны для действительных s ^

0;

(iii) I К (s + іт; X, у )\ < К (s; х, у) для всех действительных

т

и действительных s ^

0.

 

 

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

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

ЛИТЕРАТУРА

[1]Synge J. L., On the stability of a viscous liquid between rotating coaxial cylinders, Proc. Roy. Soc. (London) A, 167 (1938), 250—256.

[2]Гантмахер Ф. P., Теория матриц, «Наука», M., 1967.

[3]Jentsch R., Ober Integralgleichungen mit positivem Kern, Crelle’s J., 141 (1912), 235_244.

[4]Karlin S., Positive operators, J. Math. Mech., 8 (1959), 907—937.

[5]Красносельский M. А., Положительные решения операторных уравнений,

Физматгиз, М., 1962.

[6*]Крейн М. Г. и Рутман М. А., Линейные операторы, оставляющие инва­ риантным конус в пространстве Банаха, УМН, 3, вып. 1 (1948), 3—95.

[7*]Шефер X., Топологические векторные пространства, «Мир», М., 1971. [8*]Красносельский М. А ., Вайникко Г. М., Забрейко П. П., Рутицкий Я- Б.,

Стеценко В. Я-, Приближенные методы решения операторных уравнений,

«Наука», М., 1969.

[9*]Линь Цзя-цзяо, Теория гидродинамической устойчивости, ИЛ, М., 1958. [Ю*]Овчинникова С. Н., Устойчивость течения Куэтта в случае широкого зазора

между вращающимися цилиндрами, ПММ, 34, (1970), 302—307. [11*]Дэвис С. Г., О принципе изменения устойчивости, сб. Механика, 4 (122),

(1970), 56—73.

XV

ВОЗМУЩЕНИЕ РЕШЕНИЙ НЕКОТОРЫХ НЕЛИНЕЙНЫХ КРАЕВЫХ ЗАДАЧ

М. Мильман

1.Введение

Вэтой лекции рассматриваются три нелинейные задачи на собственные

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

2. Установившийся температурный режим в теле, содержащем нелинейный источник тепла

Пусть Т (х) — установившееся распределение температуры в трех­ мерной области D, содержащей источник тепла величины S (Т). Если область D ограничена непроводящей поверхностью В, то Т является решением следующей задачи:

АT = XS(T)

в D,

(2.1)

4 ^ = 0

на ß .

(2.2)

on

 

 

Здесь S (Т) пропорционально локальной скорости теплообразования, а X — параметр, характеризующий интенсивность источника. Нужно решить задачу (2.1), (2.2) относительно Т (х), когда 5 (Т) — нели­ нейная функция.

Предположим, что функция S (Т) имеет один и тот же знак при всех Т. Тогда, поскольку граница В теплонепроницаема, в нестацио­ нарном случае следует ожидать, что Т будет монотонной функцией времени. Поэтому, чтобы существовал установившийся режим, мы предполагаем, что S (Т) должна менять знак. Следовательно, мы считаем, что существует по крайней мере одно значение темпера­

туры Т о, такое, что

(2.3)

S (То) = 0.

Тогда при произвольном X задача (2.1), (2.2) имеет решение

 

Т(х) = Т0.

(2-4)

Будем искать теперь решения, отличные от (2.4). (Другой

подход

к этой задаче дан в лекции VIII.)

 

238

М. М И ЛЬМ А Н

Мы ищем решение Т (х , е) и X (е) задачи (2.1), (2.2), зависящее дифференцируемым образом от параметра е и сводящееся к (2.4) при 8 = 0. Представим Т (х, е) и X (е) в окрестности е = 0 рядами Тейлора по е:

Т(Х, 8) =

Г 0 + 8 Т ( х ) +

4 - 8 2Г ( х ) + . . . ,

(2.5)

Я (е) =

Яо -|~вЯ -[—2*

~f" *• • *

(2.6)

Уравнения для коэффициентов этих разложений можно получить повторным дифференцированием (2.1) и (2.2) по е, полагая затем е = 0.

Дифференцируя (2.1) и (2.2) по е один раз и полагая 8 = 0, полу­ чаем

[Д -Я

05 '(Г 0)]Т = 0,

(2.7)

- £ - = 0 на В.

(2.8)

дп

 

 

Предположим, что 5'(Т0) ф 0; тогда уравнение (2.7) с граничным усло­ вием (2.8) представляет собой линейную задачу на собственные значе­

ния относительно Х0 к Т (х), которую можно получить линеаризацией уравнения (2.1). Эта задача имеет систему ортонормированных собст­ венных функций (х), соответствующих собственным значениям Х0п, п = 0, 1, 2, . . . . Если — простое собственное значение,, то при

Ьо =

Ьо», п = 0, 1, 2, . . . .

(2-9)

задача (2.7), (2.8) имеет нетривиальные решения

 

 

Т = Ап<рп,

(2.10)

где Ап — неопределенная

постоянная.

 

Нужно отметить, что собственное значение

 

 

^оо = 9

(2.11)

всегда простое; ему соответствует собственная функция х)

Фо = У"1/3,

(2-12)

где V — объем области D. Случай кратного собственного значения рассматривается ниже.

Дифференцируя (2.1), (2.2) дважды по е и полагая е = 0, получаем

[Д - Х0S' (То)] Т = X0S" (То) Г2 + 2XS' (Г0) Т,

(2.13)

- £ - = 0 на В.

(2.14)

г) Нормированная.— Прим, перев.

XV. В О ЗМ У Щ Е Н И Е Р Е Ш Е Н И Й Н Е Л И Н Е Й Н Ы Х К Р А Е В Ы Х ЗА Д А Ч

2 3 9

Однородной задачей, соответствующей (2.13), (2.14), является все та же задача (2.7), (2.8). Поэтому краевая задача (2.13), (2.14) имеет решение в том и только том случае, если правая часть (2.13) ортого­ нальна к — соответствующей собственной функции самосопряжен­ ной однородной задачи (2.7), (2.8). Из условия ортогональности и (2.9)„ (2.10) получаем

І п = ---- b f l g ™

j cpndx.

(2.15)

 

D

 

Можно подставить теперь это значение і п в уравнение (2.13) и решить

его относительно Тп. Если, переопределяя соответствующим образом *)

параметр е, исключить из Тп произвольное слагаемое, кратное фпу

то Тп будет пропорционально А\.

 

 

 

Предположим теперь, что собственное значение Х0п имеет крат­

ность k и ему соответствуют собственные функции фп1,

. . ., фп?г,

так

что

к

 

 

 

 

 

 

 

 

Т’п = 3 ^ п у Ф пЬ

(2.16>

 

 

3=1

 

 

где

Anj — постоянные.

Соответствующие

(2.16) k условий

разреши­

мости уравнения (2.ІЗ)

имеют вид

 

 

 

 

к

 

 

 

0 = K0nS

(То) 2 ChnijAniAn] -f- 2XS' (То) Anm,

 

 

 

і, 3

 

 

где

 

m*=lf .

k.

(2.17)

 

 

 

 

 

 

= ^ фппгфпгфпу dx»

(2.18)

 

 

D

 

 

Уравнения (2.17) определяют Anj и І с точностью до произвольного общего множителя Мп, и решение, возможно, не единственно. Запи­ шем это решение в виде

= A47lp.7j) Anj = МПСП] , j = 1, . .. , k\

(2.19)

тогда

 

к

 

T n ~ M n 2 ^п^фпу'

(2.20)

3=1

 

Определение постоянных CnJ в этом случае аналогично задаче нахож­ дения «правильной» линейной комбинации невозмущенных собст­ венных состояний в теории возмущений собственного значения крат­ ности k в квантовой механике.

См. пример на стр. 241.— Прим, перев.

Соседние файлы в папке книги из ГПНТБ