Файл: Зак М.А. Неклассические проблемы механики сплошных сред.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 с. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|