- •ВВЕДЕНИЕ
- •1. ИСХОДНЫЕ ПОНЯТИЯ
- •1.1. О методе размерностей
- •1.2. Некоторые сведения из кинематики жидкости
- •1.3. Дифференцирование по времени интеграла, взятого по объему, который движется вместе с жидкостью
- •1.4. О записи уравнений движения в различных системах координат
- •2. СИСТЕМА УРАВНЕНИЙ ДВИЖЕНИЯ ЖИДКОСТИ КАК СПЛОШНОЙ СРЕДЫ
- •2.1. Уравнение закона сохранения массы
- •2.2. Уравнение количества движения
- •2.3. Уравнение энергии
- •2.4. Уравнения термодинамического состояния
- •2. 5. Вязкие напряжения и теплопроводность
- •2.6. Некоторые дополнительные сведения из термодинамики
- •3. ОСНОВЫ ДИНАМИКИ ИДЕАЛЬНОЙ ЖИДКОСТИ
- •3.1. Интегралы уравнений движения идеальной жидкости
- •3.2. Скорость звука
- •3.3. Различные формы записи уравнения Бернулли
- •3.4. Одномерное установившееся течение газа в канале переменного сечения
- •3.5. Некоторые примеры применения уравнения Бернулли
- •3.6. Установившееся истечение газа из сосуда через отверстие
- •3.7. Неустановившееся истечение газа из сосуда (опорожнение сосуда)
- •4. ВОЛНЫ В СВЕРХЗВУКОВОМ ТЕЧЕНИИ
- •4.1. О качественном различии дозвуковых и сверхзвуковых течений. Течение разрежения. Ударные волны
- •4.2. Течение Прандтля - Майера
- •4.3. Основы теории ударных волн
- •4.5. Конус в сверхзвуковом потоке
- •4.6. Теория Ньютона
- •4.7. Метод характеристик
- •Библиографический список
p =Cργ |
или |
p |
=C . |
(2.71) |
|
ργ |
|||||
|
|
|
|
||
Второе из соотношений (2.71) означает, что величина |
p / ργ |
в каждой движущейся элементарной частице не изменяется. Напомним, что этот результат получается при определенных допущениях, а именно: газ невязкий и нетеплопроводный и нет объемного поглощения или излучения тепла. Последние два допущения фактически означают, что движение каждой частицы происходит без теплообмена с окружаюшими ее другими частицами и с внешними телами. Течение, в котором справедливы эти допущения, называется адиабатическим. В этой связи соотношения (2.71) можно рассматривать как интеграл системы уравнений идеального совершенного газа при адиабатическом течении. Этот интеграл называется адиабатой Лапласа–Пуассона. Из второго из полученных соотношений (2.71) и определения энтропии (2.65) следует, что в каждой частице при ее движении не изменяется энтропия s . С другой стороны, из (2.71) следует, что в рассматриваемом случае плотность в движущейся частице зависит только от давления. Такое свойство называется баротропностью. В данном случае мы имеем баротропность для каждой элементарной частицы газа. Если константа C в (2.71) одна и та же для всех элементарных частиц или, другими словами, во всем поле течения, то вся среда (газ) называется баротропной. Очевидно, в этом случае мы имеем изоэнтропное течение.
3. ОСНОВЫ ДИНАМИКИ ИДЕАЛЬНОЙ ЖИДКОСТИ
Идеальная жидкость – это жидкость, в которой вязкие напряжения полагаются равными нулю. В такой жидкости в каждой точке действуют только нормальные напряжения, которые сводятся к давлению. Поскольку молекулярные механизмы явлений вязкости и теплопроводности одинаковы, то, не учитывая вязкость, следует пренебречь также и теплопроводностью. Для большинства задач гидромеханики, аэродинамики и газовой динамики модель идеальной жидкости – это и модель нетеплопроводной жидкости (исключение составляют задачи, связанные с течением жидких металлов, некоторые задачи динамики плазмы
163
и т.п.). Таким образом, под идеальной жидкостью мы будем понимать и нетеплопроводную жидкость. Рассматривая идеальную среду, будем также считать, что в ней отсутствует теплообмен за счет объемного поглощения или излучения энергии. Вместе с предыдущими допущениями это означает, что мы рассматриваем течения, в которых каждая жидкая частица не обменивается теплом с другими частицами и внешними телами, т.е. адиабатические течения.
3.1. Интегралы уравнений движения идеальной жидкости
Интеграл Бернулли. В данном случае ставится задача упростить исходную систему уравнений движения жидкости, что позволит не только избежать необходимости решать уравнения в частных производных или обыкновенные дифференциальные уравнения, но и получить в конечном итоге алгебраические соотношения, связывающие изучаемые величины. Пойдем по пути отыскания такого соотношения, производя при этом необходимые допущения. Допущения эти, конечно, должны иметь четкое физическое истолкование и не упускать основные свойства изучаемой среды. В результате получим соотношения, являющиеся самыми простыми в гидромеханике, вместе с тем отражающие основные закономерности движения жидкости и охватывающие широкий класс течений.
Рассмотрим сначала установившееся движение жидкости, т.е.
∂∂t = 0 . Запишем уравнение количества движения в форме Гро-
меки – Лэмба:
|
Φ + |
v2 |
|
− |
p |
. |
(3.1) |
Ω ×v = − |
|
|
|
||||
|
|
2 |
|
|
ρ |
|
|
|
|
|
|
|
|
||
Здесь сделано предположение о том, что массовые силы |
|||||||
имеют потенциал, т.е. |
|
|
|
|
|
|
|
F = − Φ . |
|
|
|
|
|
(3.2) |
Умножим скалярно векторное уравнение (3.1) на некоторое перемещение dr = dxi + dyj + dzk . В результате получим
164
|
|
2 |
|
|
dp |
= −(Ω×v )dr = − |
dx |
dy |
dz |
|
|
|
|
v |
|
|
|
|
|
|
|||
d |
Φ + |
|
|
+ |
ρ |
Ωx |
Ωy |
Ωz |
. (3.3) |
||
2 |
|||||||||||
|
|
|
|
|
vx |
vy |
vz |
|
|||
|
|
|
|
|
|
|
|
Для получения простейшего результата желательно в левой части последнего уравнения иметь полный дифференциал. Поэтому допустим, что существует некоторая функция P , удовлетворяющая условию
dP = |
dp |
. |
(3.4) |
|
|||
|
ρ |
|
Для того чтобы соотношение (3.4) имело место, необходимо выполнение условия
ρ = ρ(p). |
(3.5) |
Напомним, что среда, в которой плотность есть функция только давления, называется баротропной. Для баротропных жидкостей уравнение неразрывности и векторное уравнение количества движения, записанное в проекциях на три координатные оси, образуют замкнутую систему, так как эти четыре уравнения содержат как раз четыре искомые функции, ибо, используя (3.5), можно исключить ρ, оставив в качестве неизвестных только vx ,
vy , vz и p .
Простейшим примером баротропной среды является несжимаемая жидкость, в которой плотность вообще не зависит от давления:
ρ = const . |
(3.6) |
В этом случае легко найти выражение для функции давления P . Действительно,
P = ∫dp |
= |
p |
+ const . |
(3.7) |
|
ρ |
|||||
ρ |
|
|
|
Для совершенного газа свойство баротропности обсуждалось в конце предыдущего раздела. Было показано, что гладкое изоэнтропное течение (в котором нет скачков уплотнения и других разрывов) является одновременно и баротропным с одной и той же константой в интеграле Лапласа – Пуассона (2.71):
ϑ = |
p |
= const . |
(3.8) |
|
ργ |
||||
|
|
|
||
|
|
165 |
|
В этом случае выражение для функции P принимает вид
P = ∫ |
dp |
= ∫ |
d (ϑργ ) |
= |
γϑργ−1 |
+const = |
γ |
p |
+const . |
(3.9) |
|||
|
|
|
|
|
|
|
|||||||
ρ |
ρ |
|
γ −1 |
γ −1 |
ρ |
||||||||
|
|
|
|
|
|
В дальнейшем мы рассмотрим оба этих важнейших случая: течение несжимаемой жидкости, когда имеет место (3.6), и особенно подробно случай изоэнтропного течения совершенного газа, для которого выполняется соотношение (3.8).
Вернемся к уравнению (3.3). С учетом (3.4) его можно переписать в виде
|
|
v 2 |
|
dx |
dy |
dz |
|
|
||||
|
|
|
|
|
|
|
|
|
|
|||
d |
Φ + P + |
|
|
= − |
Ω |
x |
Ω |
y |
Ω |
z |
. |
(3.10) |
|
||||||||||||
|
|
2 |
|
|
|
|
|
|
|
|||
|
|
|
|
vx |
vy |
vz |
|
|
||||
|
|
|
|
|
|
|
Уравнение (3.10) можно проинтегрировать, если положить правую часть равной нулю. В этом случае получим знаменитый
интеграл Бернулли: |
v 2 |
|
|
|
Φ + P + |
= const . |
(3.11) |
||
2 |
||||
|
|
|
Напомним, что соотношение (3.11) было получено при следующих допущениях: жидкость идеальная, течение установившееся, массовые силы имеют потенциал, жидкость баротропная, а также при условии, что правая часть уравнения (3.10) равна нулю. Что это за условия? В соответствии с теоремой линейной алгебры определитель в правой части (3.10) равен нулю, когда:
а) |
dx = dy = dz или, что то же, dr ||v ; в этом случае (3.11) |
|||||
|
vx |
vy |
vz |
|
|
|
выполняется вдоль любой линии тока; |
||||||
б) |
dx |
= dy |
= dz |
или, что то же dr || Ω ; в этом случае |
||
|
Ω |
x |
Ω |
Ω |
z |
|
|
|
y |
|
|
(3.11) выполняется вдоль любой вихревой линии;
в) Ω||v , т.е, если вектор вихря в каждой точке параллелен вектору скорости; в этом случае (3.11) выполняется всюду в потоке;
г) Ωx = 0, Ωy = 0, Ωz = 0 или, что то же, Ω = 0 , т.е. течение
безвихревое, а значит, потенциальное; в этом случае (3.11) выполняется во всем потоке;
166
д) vx = 0, vy = 0, |
vz = 0 , т.е. |
|
когда |
жидкость |
покоится; в |
||||||||
этом случае (3.11) выполняется всюду в потоке. |
|
|
|||||||||||
Интеграл Лагранжа–Коши. |
Рассмотрим теперь неустано- |
||||||||||||
вившееся движение жидкости, т.е. когда |
∂ |
|
≠ 0 |
. Запишем урав- |
|||||||||
∂t |
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|||
нение количества движения в форме Громеки – Лэмба (2.17): |
|||||||||||||
|
∂v |
|
|
|
v 2 |
|
p |
|
|
|
|||
|
|
+ Ω×v = − Φ + |
|
− |
|
|
|
. |
|
(3.12) |
|||
|
|
|
|
|
|
|
|||||||
|
∂t |
|
|
|
2 |
|
ρ |
|
|
|
|
||
|
|
|
|
|
|
|
|
|
Мы считаем, как и при |
выводе интеграла Бернулли, что |
||
F = − Φ |
, P = |
p |
, ρ = ρ(p). Кроме того, будем по- |
|
|
ρ |
|
лагать течение безвихревым ( Ω = 0 ). В последнем случае суще-
ствует функция (потенциал скорости) ϕ такая, что v = ϕ |
. То- |
|||||||||||||
гда |
|
∂v |
= |
∂ ϕ |
= |
∂ϕ |
, и с учетом приведенных выше соотно- |
|||||||
|
∂t |
∂t |
∂t |
|||||||||||
|
|
|
|
|
|
|
|
|
|
|||||
шений |
|
уравнение |
(3.12) |
можно |
переписать |
так: |
||||||||
( |
∂ϕ |
+Φ +P + |
v2 |
) = 0 . Отсюда |
|
|
|
|||||||
∂t |
2 |
|
|
|
||||||||||
|
|
|
|
|
|
|
|
+ Φ + P + v 2 |
|
|
||||
|
|
|
|
|
|
|
|
|
|
∂ϕ |
= f (t). |
(3.13) |
||
|
|
|
|
|
|
|
|
|
|
∂t |
||||
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
Поскольку потенциал скорости ϕ задается с точностью до
произвольной функции времени, то уравнение (3.13) сводится к следующему:
∂ϕ |
+Φ +P + |
v2 |
= const . |
(3.14) |
|
∂t |
2 |
||||
|
|
|
Последнее соотношение называется интегралом Лагранжа – Коши.
В качестве примера его применения укажем задачу о неустановшемся течении несжимаемой жидкости ( ρ = const ) в поле сил
тяжести. В этом случае Φ = gz и P = ρp и уравнение (3.14)
примет вид
167