![](/user_photo/_userpic.png)
книги из ГПНТБ / Теория ветвления и нелинейные задачи на собственные значения
..pdf220 |
|
|
П. Г. Р АБ ИН О ВИ Ч |
|
|
||
b2, |
...,(mb)2 имеют |
взаимно |
непересекающиеся |
области |
значений |
||
в |
(0, |
ш„]. Отсюда |
следует, |
что |
кратное Rc |
может получаться |
|
не более чем для п точек множества Sm. Итак, доказана |
|
||||||
|
Т е о р е м а 2. Всем точкам (а, Ь) на дуге (а2 + |
Ъ2) — соп, а, Ь>0, |
|||||
кроме |
некоторого |
исключительного |
множества, |
которое |
счетно |
и состоит только из изолированных точек, соответствует простое наименьшее собственное значение задачи (3.2).
Обозначим собственное решение, соответствующее |
Rc, через ис, |
|
Ѳс, рс. Введем нормировку |
|
|
(|Ѵис|2 + Яс|ѴѲс12) d x = l . |
(3.20) |
|
Далее рассмотрим неоднородные линейные уравнения |
||
Au — Ѵр + |
ЯсѲе = f, |
|
АѲ + w = |
cp, |
(3.21) |
V-u = 0, |
|
|
где f = (/, g, h) и cp имеют вид (2.4). Добавим краевое |
условие (2.2) |
и будем искать решение, также имеющее вид (2.4). Чтобы оно суще ствовало, вектор (f, ф) должен быть ортогонален к решениям сопря женных однородных уравнений. Последние можно записать в виде
АП — ѴР + Ѳе = |
0, |
АѲ + RCW = 0, |
(3.22) |
v - u = о,
где П, Ѳ, Р тоже имеют вид (2.4) и удовлетворяют краевым усло виям (2.2). Проверка показывает, что можно положить П = ис, Ѳ = RCQC, Р = рс. Таким образом, необходимое условие, которому должны удовлетворять f, ф, принимает вид
j (uc • f + і?сѲсф) äx = 0. |
(3.23) |
т
Как показывает следующая теорема, это условие также и доста точно. Чтобы сформулировать ее, введем некоторые обозначения. Пусть Нт — замыкание множества бесконечно дифференцируемых по X, у, z функций, которые (2я/а)-периодичны по х и (2л/6)-периодич- ны по у в метрике, порождаемой скалярным произведением
(Р, |
Я)т= 2 J DapDaq d x+ J |
|
pqdx. |
|
|
іа|=та<jß |
|
d\°\ |
|
Здесь |сг| = <ті+ а2 + а3, причем au аг, а3> 0 |
, |
|||
a Da= |
||||
Обозначим норму в |
Нт через || ||т . Нетрудно |
|
показать, что || ||Л |
X I I I . Н Е Е Д И Н С Т В Е Н Н О С Т Ь П Р Я М О У Г О Л Ь Н Ы Х Р Е Ш Е Н И Й З А Д А Ч И Б Е Н А Р А 221
эквивалентна (для функций вида (2.4)) норме | |m, определяемой равенством
ОО |
1 |
|
1 |
№ = 2 |
[f l |
n(m) |
(2) j2 dz + (/2a2 + /e2ft2)m ^ p (г)2 dzJ . (3.24) |
Hjk |
|||
ft= 0 |
- 1 |
|
- 1 |
Здесь pffi означает pjh и pjk есть коэффициент Фурье функции р.
Детали мы опускаем.
Т е о р е м а 3. Если функции f, ф 6 Нт имеют вид (2.4) и удов летворяют условию ортогональности (3.23), то неоднородная систе ма (3.21) имеет единственное решение и, Ѳ, р вида (2.4), удовлетворяю щее краевым условиям (2.2). Пара и, Ѳортогональна (в смысле ( , )0) паре Uc Ѳс. Кроме того, и, Ѳ £ Нт+г, р € Нт+І и
и Іт + 2 + [ Ѳ |m+2 + [ р Im+I ^ ОС ([ f \т + I ф [т), |
(3.25) |
где а — постоянная. (Как и раньше, норма | и \т вектора и по опреде лению равна сумме норм его компонент.)
Эту теорему можно вывести из более общей теоремы Дуглиса и Ниренберга [22] или же непосредственно путем подстановки (3.14) в (3.21), сведения задачи к бесконечной системе обыкновенных дифференциаль ных уравнений, последующего доказательства существования и оце нок для этой системы и возвращения к исходной системе. Подроб ности см. в [24].
4.Итерационная схема
Вэтом параграфе мы укажем способ построения решения систе мы (3.21). Сходимость данной схемы будет установлена в § 5.
Мы знаем, что уравнения (2.1) обладают только тривиальным реше нием при R ^ Rc. Это побуждает нас разыскивать решение около
тривиального при R > |
Rc и малой |
разности R — Rc. Пусть х = |
= (X, у, z). Попытаемся найти решение в виде |
||
и (х, |
е) = еис (х) + |
е28и (х, е), |
Ѳ(х, |
е) = еѲс (х) + |
е2бѲ (х, е), |
Р (х, |
е) = грс (х) + |
(4.1) |
е28р (х, е), |
||
R (е) |
= Rc + e28tf |
(е). |
Здесь би, 6Ѳ, Ьр должны иметь вид (2.4), и пара би, 6Ѳ должна быть ортогональна к ис, Ѳс относительно скалярного произведения ( , )0. Параметр е можно считать амплитудой решения. Для сходимости рассматриваемой схемы нам необходимо предположение о малости | е |.
222 |
П. |
Г. Р А Б И Н О В И Ч |
|
|
|
Подставляя (4.1) в (4.2), находим, что |
|
|
|||
Аби — ѵбр + |
R cбѲе = [(ис + |
еби) -ѵ] (ис + |
еби) — |
|
|
|
|
|
— ебR (Ѳс + ебѲ) е, |
(4.2) |
|
|
АбѲ + w = Р [(ис + |
еби)-Ѵ] |
(Ѳс + ебѲ), |
(4.3) |
|
|
|
у б и |
= 0. |
|
(4.4) |
При этом би, 6Ѳ удовлетворяют краевым условиям (2.2).
Введем следующую итерационную схему решения этих уравнений:
би0 = 0, 6Ѳ0 = 0, |
|
|
|
Абип+1 — Ѵбрп+і + |
ДсбѲп+1е = [(uc + |
ебип)-Ѵ] (uc + |
ебил) — |
|
|
- ебRn (Ѳс + ебѲп), (4.5) |
|
ДбѲп+1 + |
бшл+1 = Р [(ис -Ь еби„) -V] (Ѳс + |
ебѲл), |
|
|
У б и п+1 = |
0. |
|
Так как уравнения (4.5) имеют вид (3.21), их нельзя решить, если правая часть не удовлетворяет условию ортогональности (3.23). В то же время мы располагаем свободным параметром бRn. Поэтому опре делим бRn, потребовав, чтобы правая часть системы (4.5) удовлетво ряла условию (3.23). Мы сумеем определить бRn таким способом, если окажется, что при образовании интеграла (3.23) коэффициент при
бRn не равен нулю Но этот коэффициент равен —е \ (Ѳс + ебѲп) wc dx.
Далее, имеем
1 2 я / а 2 п /Ь
-1 0 0
так как D, С > 0 при z 6 (—1, 1)- Оставшееся слагаемое имеет поря док е. В ходе итерации выполняется точечная оценка для 6Ѳ„. Поэтому коэффициент при бRn не равен нулю, если величина | е | достаточно мала. Таким образом, можно предполагать, что правая часть систе мы (4.5) удовлетворяет условию (3.23). Легко видеть, что если би„, 6ѲП, брп имеют вид (2.4), то это верно и для правой части системы (4.5). Поэтому, согласно теореме 3, можно провести итерацию, и бип+1, 6Ѳ„+1, 8рп+1 также будут иметь вид (2.4). Поскольку ис, Ѳс бесконечно дифференцируемы, это же верно и для бип, 6ѲП, брп согласно теоре ме 3 и неравенству Соболева.
В § 5 мы покажем, что данная итерационная схема сходится для достаточно малых е. В этом предположении доказана
Т е о р е м а 4. Существует нетривиальное решение задачи (2.1), (2.2) вида (4.1) для 0 <; е < е0. Это решение бесконечно дифференци руемо.
X I I I . Н Е Е Д И Н С Т В Е Н Н О С Т Ь П Р Я М О У Г О Л Ь Н Ы Х Р Е Ш Е Н И Й З А Д А Ч И Б Е Н А Р А 2 2 3
С л е д с т в и е . 8R (г) > 0.
Д о к а з а т е л ь с т в о . По теореме единственности Джозефа система (2.1), (2.2) имеет только тривиальное решение при R ^ . Rc x)
Прежде чем перейти к § 5, проведем более детальный анализ нашей итерационной схемы. Напомним введенное выше условие нормировки
J (VucI |*2+ .ReIѴѲ |2) dx = l.
Чё
Так как — ис, —Ѳс удовлетворяют этому условию нормировки вме сте с ис, Ѳс, можно найти второе нетривиальное решение задачи (2.1), (2.2) , используя в качестве исходных данных —ис, —Ѳс вместе ис, Ѳс. Таким образом, мы приходим к следующей теореме:
Т е о р е м а 5. Задача (2.1), (2.2) имеет по крайней мере два нетривиальных решения вида (2.4) для 0 < е С е02).
Будем различать эти два решения, обозначая первое из них через
и+(е) = |
еис + |
е2би+(е), Ѳ+(е) = еѲс + |
е2бѲ+(е), д+(е) = |
грс + |
e28p+(s), |
|||||
R+(e) = |
Rc + |
e28R+(e) |
и |
заменяя |
+ на •— для |
второго. |
Мы |
не |
||
станем подробно описывать зависимость решений от х, у, z. |
|
|
||||||||
Пусть 5 означает линейное преобразование от и (х, у, z), ѵ (х, у, z), |
||||||||||
w (х, у, |
г), Ѳ (х, у, z), |
р (х, у, г) |
к |
и (х, у, |
—г), |
ѵ (х, у, —г), |
||||
— W (х, у, —г), —Ѳ (х, у, —z), р (х, у, —г). |
|
|
|
|
||||||
Л е м м а |
9. Если и, |
Ѳ, |
р — решение |
задачи |
(2.1), |
(2.2), |
то |
5и, |
||
5Ѳ, Sp — тоже решение (при том же R). |
|
|
|
|
|
Д о к а з а т е л ь с т в о получается непосредственным подсчетом. Согласно лемме 9, каждое из уже найденных решений задачи (2.1), (2.2) приводит после преобразования 5 к новому решению. Однако на самом деле эти решения не являются новыми, так как они функ
ционально зависят от уже известных. Л е м м а 10.
Su-(e) = u+(s), Se-(e) - Ѳ+(е), Sp-(e) = р+(е),
R-(s) = R+(e).
Доказательство опускается.
З а м е ч а н и е . Так как R (е) — четная функция от е, даже без использования теоремы единственности для системы (2.1) приходим к выводу, что решения системы (2.1) могут ответвляться от нулевого решения только по одну сторону от Rc.
См. примечание на стр. 212. — Прим. ред.
2)В [37] доказано, что их ровно два. Недавно Г. К. Тер-Григорьянц показал (работа печатается в ПММ), что если отказаться от условий четности, то возни кает ровно три новых решения с точностью до сдвигов в горизонтальной пло скости. Это прямоугольная конвекция и два решения типа плоской конвекции.—
Прим. ред.
224 |
П. Г. Р А Б И Н О В И Ч |
5.Доказательство сходимости
Вэтом параграфе мы завершим доказательство существования, показав, что итерационная схема (4.5) сходится для малых | е |. При этом в сочетании с оценками (3.25) линейной теории существенно используются рассуждения, связанные со сжимающими отображе ниями. Мы будем работать с решением и+, но знак + опустим. Дока зательство для решения и- то же самое. По соображениям техниче ского порядка наряду с квадратичными оценками нам будут необхо димы также и точечные оценки. Поэтому положим
|
II Ф (х) 11= II ф II |
= |
sup |
I Ф(х) I, |
|
|
|||||||
|
|
|
|
|
|
|
X |
|
|
|
|
|
|
|
II ф 111 = |
II |
ф II + |
II ф* |
II + |
II |
ф » |
II |
+ |
II фг II. |
|
||
|
II и |
И= |
max |
(И и |
И, |
Ио ||, |
|| w ||). |
|
|
||||
Мы также воспользуемся следующей леммой: |
|
|
|
||||||||||
Л е м м а |
11. | Фф L |
< |
а (|| ф || |
| ф |m |
+ |
II Ф II |
I Ф U ), |
где т > 3 |
|||||
и а — постоянная, зависящая от т. |
|
|
|
|
|
|
|
|
|||||
Д о к а з а т е л ь с т в о |
по |
существу |
такое же, как при выводе |
||||||||||
аналогичного неравенства в [22]. Оно |
использует |
соответствую |
|||||||||||
щие интерполяционные неравенства, и мы его опустим. |
|
||||||||||||
Д о к а з а т е л ь с т в о |
т е о р е м ы |
|
4. |
Пусть |
уже |
доказана |
|||||||
оценка |
|
|
|
|
|
|
|
|
|
|
|
|
|
I би? |m+2 + |
I 6Ѳ; |m+2 + |
I бpj |
|m+1 ^ |
/Cn. |
j |
= |
1, |
• • |
n, |
(5.1) |
где Kn — постоянная и m ^ 3. Покажем, что если взять | s | доста точно малым, то можно выбрать постоянную К п = К, не зависящую от п, так, чтобы последовательности {бип, 6ѲП} и {брп) были огра ничены в Нт+г, Нт+І соответственно. Далее, имеем
SRj — ( 8 j wc(Ѳс + |
ебѲ7-) d x ) 1 j {uc-[(uc + e6u7-)-V] (uc + |
e6u,-) + |
||||
|
|
-f- R$c [Uc-j- e6uj) • V] (Ѳс -f- вбѲу)}d x . |
||||
Положим ß = |
j wcQcdx. Мы можем оценить |
Jдасб02 dx |
сверху вели- |
|||
|
<ё |
I wcI Kn- |
|
|
|
|
ЧИНОЙ I Wc 11 6 0 j К |
|
|
и |
используя |
||
Полагая |
2 max ( || uc ||ь /?с ||Ѳе||, |uc|m+1, Ѳс |m+1) |
|||||
(5.1), мы можем для бR} получить оценку |
|
|
|
|||
|
|
\èRj\ < ± ( 4 у К п + 2гК1). |
|
(5.2) |
||
Из (3.25) |
имеем |
|
|
|
|
|
I б ііп + і |m+2 + I бѲп+ l |m+l + |б р п+1 |m+1^ |
0i {[(uc -f- б б і1п ) -V] (Uc + |
eÖUn) + |
||||
|
sdRn (Ѳс -f* eö9n) -|-Pj |
(uc-)- e 6 u n )• Ѵ(ѲС-f- ббѲп) |m}. (5.3) |
X III . Н Е Е Д И Н С Т В Е Н Н О С Т Ь |
П Р Я М О У Г О Л Ь Н Ы Х |
Р Е Ш Е Н И Й ЗА Д А Ч И |
Б Е Н А Р А 225 |
||||||
Неравенство |
Соболева |
[22] |
утверждает, |
что |
|
|
|
||
|
|
II Ф \\j < |
а |
I ф Іт+/, |
т > 2. |
|
|
(5.4) |
|
Полагая о = а | (uc-V) uc |
\т + Р | (ис-Ѵ) Ѳс \т и используя (5.1), |
||||||||
(5.2), (5.4), |
а также лемму 11 для оценки правой части в (5.3), находим |
||||||||
I би„+1 |m+2 + |
I 6Ѳ„+1 |m+2 + |
I 6pn+1 |m+1 ^ |
|
|
|
||||
|
|
< |
у + |
I 8 I УіКп + |
^ г К І |
+ |
I 6 |3y3M , (5.5) |
||
где Yi, Y2, Тз — полиномы от а, |
ß-1, у. |
|
Уг» |
Уз) < |
1/4 и поло |
||||
Поэтому, если теперь взять |
[ е | шах (2у, Yi» |
||||||||
жить Кі = |
2у и Кп = 2у, то мы придем к неравенству |
|
|||||||
|
I |
lm+2 + |
I ^Ѳп+1 lm+2 “Ь I ^Pn+l Іт+г"'-^ 2у. |
(5-6) |
|||||
Таким |
образом, можно взять Кп+ 1 = 2у, и последовательности |
{8un, 6Ѳ„}, {брп} ограничены в Нт+2, Нт+І соответственно. Последо вательность действительных чисел {б/?п} также ограничена. В силу стандартных теорем о компактности можно выделить подпоследова тельность нашей последовательности, сходящуюся к решению систе мы (2.1). На самом деле вся последовательность сходится к этому реше нию системы (2.1). Чтобы в этом убедиться, положим
6U„ = би„+1 — 6un, |
6Ѳ = 6Ѳ„+1 — 6Ѳ„, бРп = |
8рп+1 — брп. |
|
Из (5.6) следует оценка |
|
|
|
I 6Un lm+2 + |
I 0Ѳ„ |m+2 + | 6Pn |m+l < |
4y. |
(5.7) |
Используя уравнения, |
которым удовлетворяют |
6Un, 8@п, |
бРп, |
и примененную выше технику, нетрудно показать, что выполняется неравенство
I би„ І2 + I бѳ„ |2 + I 5Рп І1 < |
|
|
|
< (I 8 I у4 + 82у5 + |
I 8 I3 у6) (I 6Un_! I, + I 6ѲП_! Ю, |
(5.8) |
|
где у4, у5, ув — полиномы от a, |
ß_1, у, Р, |
Rc. Поэтому, если выбрать |
|
J е I столь малым, что |
|
|
|
М у 4 + е 2у6 + | е |3у6< - ^ - , |
|
||
то будет выполняться неравенство |
|
|
|
I 6ПП|2 + I 6Ѳ„ І2 + I âPn ІІ^м-у ( I |
|і 4“ I бѲп-1 |і)? |
(5-9) |
из которого следует, что 6U„ -+- би, 6Ѳ„ -+ 6Ѳ в Я 2 и 8Рп -+ 8р в Я 4.
15—01285
226 П. Г. РА Б И Н О В И Ч
Так как, согласно (5.7), мы имеем оценки для высших производных от 61І„, 6ѲП, 8Рп, то, используя интерполяционное неравенство [22]
|
| ф Ь - < « | ф ] і - ^ |
| ф ] І /й, |
|
(5.10) |
||
заключаем, что |
би„ -> би, |
6ѲП |
8Ѳ |
в Нт+1, а |
6рп -> бр |
в Нт . |
На самом деле, |
по теореме |
Банаха — Сакса, би, |
6Ѳ 6 Нт + 2 |
и бр £ |
||
6 Нт+І. |
показано, что (би, |
6Ѳ) 6 # m+2> бр £ Нт+1, обычная |
||||
Как только |
процедура последовательного улучшения гладкости приводит к выво ду, что би, 6Ѳ, бр бесконечно дифференцируемы. Точнее, из известных свойств дифференцируемости би, 6Ѳ и лемм 5, 11 следует, что правая часть в (4.2) принадлежит Нт+1, а тогда (би, 6Ѳ) 6 Нт+з> бр € # m+a по теореме 3 и т. д.
Таким образом, доказательство теоремы 4 завершено.
Эта лекция основана на работе [23], где можно найти дальнейшие подробности. По поводу более поздних результатов см. также [24].
ЛИТЕРАТУРА
[1]Krishnamurti R., Finite amplitude thermal convection with changing mean temperature: The stability of hexagonal flows and the possibility of finite amplitude instability, dissertation, UCLA, 1967.
[2]Kosmieder E. L., On convection on a uniformly heated plane, Beiträge zur Physik der Atmosphäre, 39 (1966), 1-—11.
[3]Pearson J. R. A., On convection cells induced by surface tension, J. Fluid Mech., 4 (1958), 489—500.
[4]Chandrasekhar S., Hydrodynamic and hydromagnetic stability, Oxford
University Press, 1961.
[5] Malkus W. V. R. and Veronis G., Finite amplitude cellular convection,
J.Fluid Mech., 4 (1958), 225—260.
[6]Busse F., Das Stabilitätsverhalten der Zellularkonvection bei endlicher Amplitude, dissertation, University of Munich 1962. (Translation from Ger man by S. H. Davis, Rand Corp., Santa Monica, Calif. 1966.)
[7]Segel L. A., The structure of nonlinear cellular solutions of the Boussinesq equations, J. Fluid Mech., 21 (1965), 345—358.
[8]Chorin A. J., Numerical study of thermal convection in a fluid layer heated from below, New York University, Courant Institute of Mathematical Scien ces technical reports, 1966.
[9]Joseph D. D., On the stability of the Boussinesq equations, Arch. Rational Mech. Anal., 20 (1965), 59—71.
[10]Veite W., Stabilitätsverhalten und Verzweigung Stationärer Lösungen der Navier-Stokesschen Gleichungen, Arch. Rational Mech. Anal., 16 (1964),
97—125.
[11]Veite W., Stabilität und Verzweigung stationärer Lösungen der Navier-Sto kesschen Gleichungen beim Taylorproblem, Arch. Ratiohai Mech. Anal., 22
(1966), |
1—14. |
[12]Юдович В. И., Пример рождения вторичного стационарного или периоди ческого течения при потере устойчивости ламинарного течения вязкой жидкости, ПММ, 29, вып. 3 (1965), 453—467.
[13]Ладыженская О. А., Математическая теория вязкой несжимаемой жидко сти, Физматгиз, М., 1961.
[14]Курант Р. и Гильберт Д., Методы математической физики, т. 1, Гостех-
издат, М., 1951.
X III . Н Е Е Д И Н С Т В Е Н Н О С Т Ь П Р Я М О У Г О Л Ь Н Ы Х Р Е Ш Е Н И Й ЗА Д А Ч И Б Е Н А Р А 227
[15] Kirchgässner К-, Die Instabilität der Strömung zwischen zwei rotierenden Zylindern gegenüber Taylor-Wirbeln für beliebige Spaltbreiten, Z.A.M .P., 12 (1961), 14—30.
[16]Karlin S., Total positivity and applications, Stanford University Press, Stanford, Calif., to appear.
[17]Крейн М. Г. и Рутман М. А., Линейные операторы, оставляющие инва
риантным конус в пространстве Банаха, УМН, 3, вып. 1 (1948), 3 —95.
[18]Reid W. Н. and Harris D. L., Some further results on the Bénard problem, Phys. Fluids, 1 (1958), 102—110.
[19]Данфорд H. и Шварц Дж. T., Линейные операторы, Общая теория, т. I,
ИЛ, М., 1962.
[20]Rabinowitz Р. Н., Periodic solutions of nonlinear hyperbolic partial diffe rential equations, Comm. Pure Appl. Math., 20 (1967), 145—205.
[21]Nirenberg L., On elliptic partial differential equations, Ann. Scuola Norm. Super. Pisa, Ser. 3, 13 (1959), 1—48.
[22]Douglis A. and Nirenberg L., Interior estimates for elliptic systems of partial differential equations, Comm. Pure Appl. Math., 8 1955), 503—538.
[23]Rabinowitz P. H., Nonuniqueness of rectangular solutions of the Bénard
problem, Stanford University, technical report. |
rectangular solutions |
|
[24] Rabinowitz P. |
H., Existence and nonuniqueness of |
|
of the Bénard problem, Arch. Rational Mech. Anal., 29 (1968), 32—57. |
||
[25*]Сорокин В. С., |
О стационарных движениях жидкости, |
подогреваемой снизу, |
ПММ, 18, вып. 2 (1954), 197—204.
[26*]Горьков Л. П., Стационарная конвекция в плоском слое жидкости вблизи критического режима теплопередачи, ЖЭТФ, 33, № 2 (8) (1957), 402—411.
[27*]Сорокин В. С., Вариационный метод в теории конвекции, ПММ, 17, вып. 1
(1953).
[28*]Уховский М. Р. и Юдович В. И., Об-уравнениях стационарной конвекции,
ПММ, 27 вып. 2 (1963), 197—204.
[29*]Юдович В. И., О возникновении конвекции, ПММ, 30, вып. 6 (1966). [30*]Юдович В. И., Математические вопросы гидродинамической теории устой
чивости, Тезисы докладов Всесоюзной межвузовской конференции по при менению методов функционального анализа к решению нелинейных задач,
Баку, 1965.
[31*]Юдович В. И., О бифуркации вращательных течений жидкости, ДАН
СССР, 169, № 2 (1966), 306—309.
[32*]Юдович В. И., Вторичные течения и неустойчивость жидкости между вра щающимися цилиндрами, ПММ, 30, вып. 4 (1966), 688—698.
[33*]Иванилов Ю. П. и Яковлев Г. И., О бифуркации течений жидкости между вращающимися цилиндрами, ПММ, 30, вып. 4 (1966), 768—773.
[34*]Ворович И. И., Некоторые вопросы устойчивости оболочек в большом,
ДАН СССР, 122, № 1 (1958), 37—41.
[35*]Ворович И. И., Некоторые оценки числа решений для уравнений Кармана в связи с проблемой устойчивости пластин и оболочек. Проблемы гидроди намики и механики сплошной среды (к 60-летию академика Л. И. Седова), «Наука», М., 1969, стр. 111—118.
[36*]Юдович В. И., Устойчивость конвекционных потоков, ПММ, 31, вып. 2
(1967), 272—281.
[37*]Юдович В. И., Свободная конвекция и ветвление, ПММ, 31, № 1 (1967),
101— 111.
[38*]Ландау Л. Д ., Лифшиц Е. М., Механика сплошных сред, ГИТТЛ, М., 1953.
15*
XIV
ИЗМЕНЕНИЕ УСТОЙЧИВОСТИ В ТЕЧЕНИИ КУЭТТА
Г. Ф. Вайнбергер
1. Введение
Рассмотрим течение вязкой несжимаемой жидкости между двумя вра щающимися цилиндрами. Внутренний цилиндр имеет радиус т] < 1 и угловую скорость й, а внешний цилиндр — радиус 1 и угловую скорость рй, причем р > О (см. рис. 1).
I
В цилиндрических координатах с единицами измерения, выбран ными так, чтобы вязкость ѵ = 1, уравнения Навье — Стокса имеют вид
диг |
|
„ 2 |
|
|
|
|
|
|
|
|
2 âuQ |
ur |
|
|
|
иѳ |
|
|
|
|
|
|
|
|
) , |
(1.1a) |
|||
dt |
(U-V)«r---- і~ = |
âr |
(\ гp' ) + ( д“г ' |
|
ѲѲ |
|
||||||||
|
r |
|
|
|
|
|
||||||||
дщ |
(u - v ) |
итщ |
|
- ! JL |
|
|
|
|
|
|
|
|||
dt |
Uq |
|
|
( f ) + |
|
|
|
|
|
|||||
|
|
|
|
г |
0Ѳ |
2 |
dur |
щ |
|
|
||||
|
|
|
|
|
|
|
+ |
( А^ѳ -b |
|
(1.1b) |
||||
|
|
|
|
|
|
|
r2 |
de |
г2 |
|
||||
|
|
du. |
.(u.V) u z= |
|
— ^ - ( j - ) |
+ Auz, |
|
|
(1.1с) |
|||||
|
|
dt |
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
где |
|
|
|
â |
|
щ |
|
|
|
|
|
|
|
|
|
|
U‘V = u |
|
d |
, |
|
â |
|
|
|
||||
|
|
|
—“Г |
r |
de |
+ UZ |
dz |
|
|
|
||||
|
|
|
|
dr |
|
|
|
|
|
|
||||
|
|
Д = |
ö2 , |
1 |
d |
|
, |
1 |
d2 |
|
a2 |
|
|
|
|
|
a/-2 |
|
âr |
|
|
г2 |
аѳ2 |
|
dz* |
|
|
|
|
Уравнение неразрывности |
имеет вид |
|
|
|
|
|
|
диГ |
. ит . |
1 а«ѳ , |
dr |
г |
аѳ |
duz |
0. |
( 1 . 2) |
|
dz |
|||
|
|
XIV. ИЗМЕНЕНИЕ УСТОЙЧИВОСТИ В ТЕЧЕНИИ КУЭТТА |
2 2 9 |
Здесь ur, uq, tiz — физические компоненты скорости, а р — давле ние. Эти уравнения допускают точное решение
ur = 0, ыѳ = Лг + -^-, uz= 0, |
= j (ul/r)dr. |
(1.3) |
Граничные условия для вязкого течения состоят в том, что ско рость жидкости на твердых стенках должна равняться скорости стенок. В нашей задаче эти граничные условия на цилиндрах г — ц и г = 1 выполняются, если положить
|
b = qV |
' |
(1- 4) |
и |
Если цилиндры имеют конечную длину, |
скажем | г | ^ |
L, то (1.3) |
(1.4) доставляют решение краевой задачи, |
изображенной на рис. 1, |
||
только в том случае, когда щ при г = + L |
задано равенствами (1.3) |
||
и |
(1.4). |
|
|
Течение, заданное равенствами (1.3) и (1.4), называется течением Куэтта. Мы хотим изучить условия, при которых это течение устой
чиво относительно бесконечно малых возмущений. |
|
2. Задача о возмущениях |
|
Положим |
|
цг = е«, иѳ = Аг + В/г+гѵ, uz=aw, |
^ |
(р/р) = j {{Ar + В/г)21г] dr + ел. |
|
Будем разыскивать решения и, v, w, л, обладающие осевой симмет рией, т. е. не зависящие от Ѳ.
Подставим (2.1) в (1.1) и (1.2). Приравнивая нулю коэффициенты при е, получим уравнения, которым подчиняются осесимметричные
возмущения: |
|
|
|
|
|
І ж - ь + ± ] |и —2(Л + В/г*)о = |
W ’ |
(2.2а) |
|||
[ ± |
- ь + ± \ |
ѵ+ 2Аи —0, |
|
(2.2Ь) |
|
дл |
|
|
|||
[ |
т |
- д ] “”= - |
|
(2.2с) |
|
~~~дГ ’ |
|
||||
|
|
|
|
|
|
|
|
|
— 0. |
|
(2.3) |
Действуя на |
|
(2.2а) оператором (д2/дг2), на |
(2.2с) — оператором |
—(д2/дгдг), складывая |
результаты и используя (2.3), получаем урав |
|||
нение четвертого порядка |
|
|
|
|
1 \ г д |
/ . |
1 \ 1 |
п |
В \ д*ѵ |