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