Файл: Зак М.А. Неклассические проблемы механики сплошных сред.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 26.07.2024

Просмотров: 87

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

правление распространения волн, X— скорость распространения волн. Добавляя к этим соотношениям вытекающее из теории мелкой воды условие

dp

+ ' д(pV)

dt

L dx

J

получаем формулу для скорости X:

 

^ = v c* ± ] / g ( r i +

h) + j (vr*f d \ . (Ill—6—48)

Если игж=0, то'приходим к классическому варианту [15]. Од­

нако получающееся в этом случае движение неустойчиво (ввиду

зависимости X от т] происходит опрокидывание волн, сопровож­ дающееся турбулизацией их гребней, т. е. возникновением ско­ ростей vr). В то же время движение с турбулизовавшимися греб­ нями становится устойчивым, так как падение X из-за уменьше­ ния высоты волн компенсируется возрастанием X вследствие по­

явления пульсационных скоростей vr в соответствии с (48).

В заключение рассмотрим модель

турбулентной

жидкости,

в которой учитываются электромагнитные эффекты.

Для этого

достаточно к силеFc уравнения (42)

добавить силу Лоренца

Fn = yrotB X B , где В ^индукция магнитного поля. Останавли­

ваясь на случае бесконечной электропроводимости, запишем

[9]

-45- =

rot (vc X В), divB = 0.

(Ill—6—49)

Рассматривая (42),

(43),

(49) и уравнение баланса массы,

нетрудно придти к формулам

(45) и (47). Однако для скорости

распространения волн разрыва rot vc получается

формула,

от­

личная от (46). Действительно, проектируя (42)

(с учетом

Fn)

и (43) на нормаль п к поверхности разрыва с касательной г и учитывая, что в силу уравнения баланса массы и второго равен­ ства в (48) [ducVdrJ^O, [д5т/дт]=0, а в силу кинематических со­ отношений на фронте разрыва [<3/д/г]=0, получим

1

Отсюда

' = 'Ос ±

J (О2&

Сравнивая эту формулу с (46), молено констатировать, что магнитная индукция увеличивает упругость волн разрыва rotvc

115


и далее сохраняет их упругость при и,-=0. Полученный результат хорошо согласуется с известным из магнитной гидродинамики фактом стабилизирующего влияния магнитного поля, повышаю­ щего устойчивость ламинарных течений лшдкости.

З а м е ч а н и е . Уравнения (42), (43) и баланса массы при физически оправданных граничных условиях могут все лее иметь бесчисленное миолеество решений, как и уравнение Эйлера (см. п. 10), если рассматривать, например, задачу об обтекании твер­ дого тела. Такая ситуация возникает из-за того, что на границе обтекаемого тела молено однозначно задать лишь нормальную составляющую скорости, в то время как касательная составляю­ щая v^° остается неопределённой. Так как произвольное задание v 1° не противоречит ни граничным условиям, ни исходной систе­ ме уравнений, появление мнолеества решений на мнолеестве {vT0} неизбелено. Для выделения из этого мнолеества ' единственного решения следует обратиться вновь к принципу наименьшего принуледения, выбирая такое vx+°, которое минимизирует принуледение Ф на {vt0}, т. е.

®(v0T+) = inf ®(v°t).

Напомним, что в п. 10 именно таким путем выделялось един­ ственное решение для плоского безвихревого обтекания, так как в этом случае неопределенность Vt° проявлялась в неопределен­ ности циркуляции Г.

Отметим также, что при решении уравнений (42), (43) для задач, в которых причиной появления vr является потеря устой­ чивости однозначного поля скоростей (например, возникновение, конвекции, размывание поверхности тангенциального разрыва), начальные и граничные значения vr имеют вероятностную приро­ ду. Для решения таких задач в детерминированной постановке мол<но рассуждать следующим образом: появление vn в соответ­ ствии с приведенными выше результатами, стабилизирует двилсение лшдкости, делая его устойчивым в классе многозначных фун­ кций v(g); пусть {vr*} — мнолсество начальных и граничных ус­ ловий, при которых во всех точках исследуемой области эта устойчивость имеет место; тогда истинный вектор vr+* молено найти из минимизации принуждения

Ф(vr+*) = inf Ф (vr*).

15.Упрощенная .турбулентная модель жидкости. Эта модель может быть получена на базе уравнений (42), (43) в результате

приближенной аппроксимации: vr(g) = v r° sign (g— 0,5), где vr° не зависит от g, т. е. замены многозначного поля скоростей дву­ значным. Тогда

~jf~ + (vc ’ V) vc + (v°r• v ) vr° = Fc -

- y grad p,

+ (vc• V) vr° + (v,° • V) vc = -

x vr°.

116


В этом приближении движение жидкости характеризуется двумя векторными полями: vc й vr. Подчеркнем, что векторы скоростей vc и vr можно складывать только «энергетически» (i\ е. vc2 + vr2 =

= 2 Эц/р).

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

кающие из сохранения массы, количества движения

и энергии

при переходе через эту поверхность [16]:

 

[рХ'] =0,рХ' [ч)п\ = \р], [г>2]/ 2 = [да].

(111— 6 — 5 0 )

Здесь n" = v-n — проекция скорости v на внешнюю нормаль п к поверхности разрыва, w — тепловая функция, А'=А— vn, где А— скорость перемещения поверхности разрыва.

Интегрируя (50) по £ от 0 до 1, получим

[ рУ ] = 0 , РУ [ v *\ = [ р ] + р v rn ° [г»,л° ] ,

РХ1

У с 2+

У г™

I______

] _ _

(111—6—51)

 

2

 

' к

— 1

 

 

г д е

Vcn =

vc-ri,

vrn°

 

 

Следствием равенств (50), (51) является условие

 

 

 

У К „ ° ] = ® гп ° К 1 -

 

Подставляя его в (50), найдем

 

 

 

Р У К " ] = [р] + Р К , : 0) 2 М / У -

 

или, учитывая,

что [А]=0 и вводя обозначение [х\ = х+Х-,

пере­

пишем

 

 

 

 

 

 

р+[ху — (У

 

} ( х у - х у ) = - ( р + -/7_).

 

Но р+Х/с4- = р_Х'с-,

поэтому последнее выражение можно

запи­

сать так:

 

 

 

+Хр—+хр

lpx+] l'cJ’

 

 

 

[ У ]

 

 

 

 

 

[ ] [АСН

[p ip - п /1

 

т. е.

 

 

 

 

 

 

 

y ,2- ( ^

V

+ р^.

 

 

 

 

 

 

р+

 

 

Для ненулевых разрывов, т. е. при [Ас']#0, имеем

V = ± / ^ . M + ( ^ V -

(III—6—52)

Эта формула отличается от классической добавочным слагае­ мым, исчезающим при t>r=0.

117


Нетрудно проверить, что при [/?], [р]->-0 формула (52) перехо­ дит в (45), если последнюю отнести к рассматриваемой упрощен­ ной модели.

Наконец, привлечение третьего равенства в (50) позволяет найти соотношение между давлениями и плотностями на фронте разрыва, которое отличается от классического и дается аналогом адиабаты Гюгонио.

Действительно, введем плотность потока газа через поверх­ ность разрыва / = р+о+=р_о_, тогда из (52) молено получить

Р = -£±£^ L +

Р+2 K

, V

= ^

+ Р-2 М

Л

(III—6—53)

а третье равенство в (50)

записать так:

 

 

J

^ ’ + [К „ °)2] = - 2 И -

 

(III- 6 - 5 4 )

Подставляя (53) в (54), найдем

 

 

 

?4Р~р! - 7 - ',~)

(_±_ _

 

) + 2 ( V n V _ 2 м

=

_ 2 [« ],

(-£ - +

Ш = 2 I [®1 +

К ^.ТП -

('III—6 - 55}

При иг=0 (55) переходит в классическую ударную адиабату. Выведем аналог ударной поляры для косой ударной волны.

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

tg ? =

^ - - ^ cos*-,

(III—6—56)

ь т

vc+ sin у.

4

'

где х — угол между скоростью газа vc+ за фронтом волны с каса­ тельной к этому фронту.

Из (53) с учетом свойств идеального газа следует

(vc+ cos X sin у — vc+ sin х cos <р)2 + (гуп0^ 2 _

 

(у с- sin о)2 - f (vrrp ) J

 

6— 1

,

2М]2

(III- 6 - 5 7 )

k + 1

(k + 1) sin2 tf

 

где Mi — число Маха перед фронтом волны, k — показатель ади­ абаты, v°rn — проекция vr° на нормаль к фронту волны. Исклю­ чая из (57) посредством (56) угол ф, придем к зависимости для угла х, или к зависимости между uc+cosx и oc+sinx, переходя­ щей в классическую ударную поляру при vr=0.


 

 

 

 

 

 

 

 

 

Л И Т Е Р А Т У Р А

 

 

 

 

 

 

 

1. С о б о л е в

С.

 

Л.

Уравнения

математической

физики. М.,

Гостехиздат,

1.954. 444

с.

О

потере

устойчивости формы идеально гибкой нити. —

2. 3 а к

М.

А.

«Прикладная

математика и

механика»,

1968, № 6, с. 1092— 1096.

 

 

3. 3 а к

М.

А.

Разрывные

движения

в механике

сплошной

среды. — «Ме­

ханика твердого тела»,

1971, № 3, с. 199—200.

 

множеств. М.,

Гос­

4. А л е к с а н д р о в

П. С.

Введение

в общую

теорию

техиздат, 1948. 411 с.

упругости. М., «Наука»,

1970. 939

с.

 

 

 

5. Л у р ь е

А.

И.

Теория

 

 

 

6. Б у л г а к о в

 

Б.

В.

Колебания. М., Гостехиздат, 1954. 891 с.

 

 

M..t

7. Р а ш е в с к и й

 

П.

К.

Риманова

геометрия

и

тензорный

анализ.

«Наука», 1967. 664 с.

 

Курс

дифференциальной

геометрии. М., Гостех­

8. Р а ш е в с к и й

 

П.

К.

 

издат, 1956. 420 с.

 

 

 

 

 

 

 

 

 

 

 

492

с.

9. С е д о в

Л.

И. Механика сплошной среды. М., «Наука», 1970.

10. С о к о л о в с к и й

В.

В.

Статика сыпучей среды. М., Гостехиздат,

1954.

275 с.

 

 

 

И.

Г.

Лекции об

уравнениях с

частными производными.

11. П е т р о в с к и й

М., Физматгиз, 1961. 400 с.

 

 

А. А. Уравнения

математической

фи­

12. Т и х о н о в

А.

Н.,

С а м а р с к и й

зики. М., Гостехиздат, 1953. 679 с.

 

 

Судпромгиз, 1958.

370

с.

13. Н о в о ж и л о в

В.

В.

Теория

упругости. Л.,

14. 3 а к

М.

А.

Единственность

и

устойчивость

решения задачи

о

малых

возмущениях

гибкой нити со

свободным концом. — «Прикладная

меха­

ника»,

1970, № 6, с. 1048— 1052.

М.

Механика

сплошных

сред.

М.г

15. Л а н д а у

 

Л.

 

Д.,

Л и ф ш и ц

Е.

Гостехиздат, 2954. 795 с.

 

 

 

 

 

 

 

 

 

 

1957.

16. Л о й ц я н с к и й

Л.

Г. Механика жидкости и газа. М., Гостехиздат,

784 с.