Файл: Зак М.А. Неклассические проблемы механики сплошных сред.pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 26.07.2024
Просмотров: 99
Скачиваний: 0
Fu= —pdv/dt. Появление нового вектора v требует привлечения дополнительных соотношений, а именно кинематических соотно шений совместности перемещений и поворотов (I—2—41), (I—2—43). Таким образом, к замкнутой системе уравнений ста тики § 2, п. 3 добавляются двенадцать скалярных уравнений, со держащих двенадцать новых переменных Vi, TV. Тем самым замкнутость уравнений динамики выполняется. Граничные усло вия полученной системы с частными производными первого по
рядка по координатам и времени сохраняются в |
форме (2—5). |
|
2. |
Модель ламинарной вязкой окидкости. |
Уравнения движ |
ния ламинарной вязкой жидкости, строго говоря, нельзя полу чить из уравнений равновесия (2—4), так как в процессе движе
ния здесь сопутствующая система координат |
детерминирована |
|
лишь в бесконечно малых интервалах времени. |
Поэтому |
обра |
тимся к принципу наименьшего принуждения |
(§ 1, п. 1). |
Запи |
шем уравнения связи (I—3—18) в виде S[(dv*/dxj) • (dv*/dXj)] = 0;
умножим этЮ выражение на неопределенные непрерывные ска
лярные функции Тц'. |
Очевидно, |
что |
образуют см. (2—2а) |
координаты некоторого тензора |
Т. |
Поэтому можно записать |
|
|
|
|
(11-3-1) |
Проинтегрируем (4) |
по У и сложим с бФа= 0: |
Проделав здесь те же преобразования, что и при выводе (2—4), получим
P4 f = d i v r + F , |
(11-3-2) |
причем граничные условия по-прежнему |
даются соотношением |
(2—5). Уравнения (2) замкнуты в силу |
уравнений состояния, |
связывающих компоненты тензора напряжений Т и тензора ско ростей деформаций С, которые выражаются через производные по координатам от проекций v в соответствии с (I—3—10), а также с учетом соотношений, связывающих плотность и давле ние: ф=ф(р,р). Эти уравнения позволяют определить поле ско ростей и поле тензора напряжений ламинарной вязкой жидкости. Если в уравнение состояния входит температура, то следует до бавить еще уравнение баланса энергии и уравнение теплопровод ности.
Заметим, что в уравнениях (2) скорость v детерминирована в пространстве D, но не в пространстве D0, т. е. функция v (г0, f) может оказаться всюду разрывной. Поэтому ускорение а в этом случае следует вычислять по схеме (I—3—14).
50
Положим, |
что |
рассматриваемая |
жидкость |
несжимаема. |
|
В этом случае все предыдущие выводы остаются в силе с той |
|||||
лишь разницей, что ввиду односторонности связи |
(I—3—20) сле |
||||
дует считать р^О . |
В частности, если div v^O , то величина divv |
||||
перестает быть недетерминированной и р==0; |
|
||||
3. |
Модель идеальной жидкости. |
Уравнения движения идеаль |
|||
ной жидкости, |
вообще говоря, нельзя получить из уравнеций рав |
новесия (2—4) путем добавления сил инерции, так как в идеаль ной жидкости сопутствующая система не определена; поэтому не справедливы вариационные формулировки типа (1—8). Обра тимся к основной функциональной формулировке принципа наи меньшего принуждения, записанного в форме (1—5). Будем ис ходить из представления v* = vc*+vr*, v+= v+c-j-v+r, где vc*, vr* —
возможные значения скоростей |
vc и vr (g), которые |
вводятся |
||
в соответствии с гл. |
I, § 3, п. |
3. |
То же относится и к |
скорости |
свободного движения v+. |
|
|
|
|
Обратимся к функционалу |
(1—7). Для случая V,= V можно |
|||
записать |
1 |
|
|
|
|
|
|
|
|
Фа = |
ffp|a*(?)— F [2д Ы V -> min. |
(II—3^—3) |
||
|
ко |
|
|
|
Минимизация функционала (3) |
должна проводитНся |
с учетом |
||
условия |
|
|
|
|
1 |
|
|
1 |
|
div J v (?) |
= ---- -l |
^ |
+ J v (?) d? • grad Pj |
(11—3—4) |
о |
|
|
0 |
|
и уравнения состояния. Если последнее содержит температуру, то необходимо добавить уравнения теплопроводности'и теплово го баланса.
Результатом минимизации должно явиться нахождение ис
тинного поля скоростей v (г, |
t, g). Напомним, что в соответствии |
с введенной в гл. I, § 3, п. 3 |
моделью идеальной жидкости каж |
дой точке пространства сопоставляется не одно значение скоро сти, а функция v (g).
Из условий (3) и (4) можно перейти к дифференциальным уравнениям для идеальной жидкости. Действителыщ, пользуясь
непрерывностью V? (г), из |
(3) полуйаем |
1 |
|
JJ(pa — F)-8Ev * d ? d I/= 0 . |
|
ко |
; |
Учитывая, что вариационная формулировка связи для идеаль ной жидкости может быть записана в виде равенства
1 |
1 |
8? div J v (I) d? = |
| div 8£ v (?) d? — 0, |
4* 51
умножим это равенство на некоторый скаляр [—р (г, f)]. Прини мая во внимание, что
I р div 8v* d V — \ div(/?8v*) d V — j (gradp)-b\*dV,
v v v
fdiV (jOov*) r f l/= fpn-8v*rfo,
V a
где n — нормаль к поверхности a, ограничивающей объем V, и полагая, что на этой поверхности S$v* = 0, получаем
1 |
grad/?)-8e\* d U V = 0 , |
|
j J (ра — F + |
|
|
/ о |
5 |
|
откуда следуют уравнения движения идеальной жидкости*: |
||
p(r, t) а (г, t, |
£)= F — grad р (г, t). |
(II—3—5) |
Здесь р — гидростатическое давление.
Эти уравнения замыкаются добавлением уравнения баланса массы (4) и уравнения состояния. Если положить 5v/5£=0, то из
(5) получим уравнения Эйлера для идеальной жидкости:
Р (г, t) а (г, t) — F — g rad /;(г, t). |
(II—3—6) |
Использование уравнений (5) вместо уравнений (6) должно |
|
привести к новым результатам при определении |
сопротивления |
тел, движущихся в идеальной жидкости, по следующим сообра жениям. Если вычислить кинетическую энергию идеальной жид кости, то можно получить соотношение
2ЭК= J J ръЧ Ы У = Ч J p |v c + v r |a rf£rfl/=
V0 ио
= J 'pv2ed V + U pI v \d 'd V ,
|
V |
V 0 . |
|
причем |
1 |
|
1 |
\ \ p \ c- \ rdbdV— \ |
p4cA \ rd'idV = 0. |
||
|
но |
H |
O |
Из этого соотношения следует, что кинетическая энергия со |
|||
держит дополнительное слагаемое, |
отражающее наличие неде |
терминированных составляющих поля скоростей, на возникнове ние которых должна затрачиваться работа.
С формально математической точки зрения различие между
уравнениями (5) и |
(6) состоит в том, что уравнения (6) |
миними |
|
зируют функционал |
(3) в классе дифференцируемых |
функций |
|
v (г), а уравнения |
|
(5) •—в классе дифференцируемых |
функций |
divv. Последний класс функций является более широким, и, сле
довательно, решения уравнений (5) доставляют |
более «глубо |
|
кий» минимум функционалу (3), нежели решения уравнений |
(6). |
|
Остановимся, наконец, на математических предпосылках воз |
||
никновения зависимости v от £. Из (5) следует, |
что такая зави |
|
* Так как v определяется тремя компонентами о,-, то |
в дальнейшем, |
не |
меняя обозначений, под g будем понимать совокупность трех независимых чи сел g = {gi, | 2, |з>, O s £ ii< l, полагая, что =
52
симость может возникнуть или в силу специфических граничных, или в силу специфических начальных условий. Первый случай будет проиллюстрирован в данном пункте. Второй случай, свя занный с потерей устойчивости решения в классе дифференциру емых функций v (г), будет рассмотрен в гл. III.
Итак, классические уравнения Эйлера для идеальной жидко сти недостаточны для полного описания движения, так как не дают возможности определения недетерминированных составля ющих поля скоростей. Попытаемся дать физическую интерпрета цию появления недетерминиров.анных составляющих скоростей в идеальной жидкости.
Согласно классической механике сплошных сред идеальная жидкость характеризуется отсутствием Dev Г; но тогда в силу принципа освобождаемости от связей неизбежно должны возник нуть дополнительные степени свободы, проявляющиеся в возник новении недетерминированных, т. е. разрывных на трехмерных континуумах, составляющих скоростей vr. Поэтому уравнения идеальной жидкости, записанные в общепринятой форме (6), не являются полными и, строго говоря, могут использоваться лишь при условии, что вязкость, а следовательно, и компоненты Dev Т бесконечно малы, но не в точности равны нулю и получаются предельным переходом из (2). Отмеченное обстоятельство хоро шо иллюстрируется, например, тем, что уравнения идеальной жидкости с разрывными начальными условиями", записанные в форме (6), не имеют единственного решения; решение стано вится единственным лишь при введении бесконечно малой вязко сти.
Все приведенные выше рассуждения остаются справедливы
ми и для несжимаемой жидкости при условии, что |
0, |
так |
как |
|
связь (I—3—22) является односторонней. |
Если |
div v > 0 , |
то |
|
div v перестает быть детерминированной |
и р = 0; |
среда |
превра |
щается в совокупность свободных точек, в которой нельзя про вести не только индивидуализацию кривых и поверхностей, но также и индивидуализацию объемов (разбрызгивание, кавита ция) . Для описания движения в этом случае следует обратиться к обыкновенным дифференциальным, уравнениям движения сво бодных точек (1—2).
В заключение данного пункта приведем пример, иллюстри рующий возможность изучения неклассических течений жидко сти посредством принципа наименьшего принуждения.
П р и м е р. Рассмотрим плоское стационарное течение идеаль ной несжимаемой жидкости без массовых сил с граничными ус ловиями v = eXr при х = ± у , где е — постоянный вектор, направ ленный вдоль оси z, г — радиус-вектор точек жидкости, лежащий в плоскости х, у.
Если искать решение в классе дифференцируемых функций и, следовательно, исходить из уравнений Эйлера, то комплексную
53