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

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

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

220

 

 

П. Г. Р АБ ИН О ВИ Ч

 

 

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Рп -+ в Я 4.

15—01285

226 П. Г. РА Б И Н О В И Ч

Так как, согласно (5.7), мы имеем оценки для высших производных от 61І„, 6ѲП, 8Рп, то, используя интерполяционное неравенство [22]

 

| ф Ь - < « | ф ] і - ^

| ф ] І /й,

 

(5.10)

заключаем, что

би„ -> би,

6ѲП

в Нт+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),

101111.

[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 ) +

 

 

 

 

 

 

 

 

 

г

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

п

В \ д*ѵ

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