![](/user_photo/_userpic.png)
книги из ГПНТБ / Теория ветвления и нелинейные задачи на собственные значения
..pdf230 |
Г. Ф. ВАЙНБЕРГЕР |
Мы предполагаем, что цилиндры бесконечные. Граничные условия прилипания требуют, чтобы возмущения к и п удовлетворяли равен ствам
|
|
и = ди/дг = и — 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) = —2А |
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] применяли теорему Принсгейма описанным выше способом для доказательства того факта, что теорема Перрона — Фробениуса остается в силе для таких интегральных опера
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-> + °о , мы можем доказать существование такого t±> 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.) |
|