Механика жидкости и газа. В 4 ч. Ч. 2
.pdf
6. УРАВНЕНИЯ ДВИЖЕНИЯ ВЯЗКОЙ ЖИДКОСТИ (УРАВНЕНИЯ НАВЬЕ–СТОКСА)
Уравнения движения вязкой жидкости можно получить из уравнений движения в напряжениях, выполнив некоторые преобразования. Рассмотрим лишь одну проекцию этих уравнений (на ось Х):
dux |
|
1 |
|
|
|
|
|
|
X |
|
pxx |
|
yx |
zx . |
|||
|
|
|
||||||
dt |
|
|
x |
y |
z |
|
||
|
|
|||||||
Как было показано при рассмотрении модели вязкой жидкости, нормальные напряжения
pxx p 2 ux 2 div u .x 3
Для упрощения задачи будем считать жидкость несжимаемой ( div u 0 ), тогда
p |
|
|
|
|
|
|
|
u |
x |
|
|
p |
|
|
|
|
|
2u |
|
|||||||||
|
|
|
xx |
|
|
|
|
|
p 2 |
|
|
|
|
|
2 |
|
|
|
x . |
|||||||||
|
|
|
|
|
|
|
|
|
|
x |
|
|
|
|
||||||||||||||
x |
x |
|
|
|
|
x |
|
|
|
|
|
|
|
2 |
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
x |
|
|
Касательное напряжение |
|
|
|
|
u |
y |
u |
|
|
; |
|
|
||||||||||||||||
yx |
|
|
|
|
x |
|
|
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
y |
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
x |
|
|
|
|
|
|
|||||
|
yx |
|
|
|
|
u |
y |
u |
x |
|
|
2u |
y |
|
|
|
2u |
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
x , |
||||||||||
|
|
y |
|
y |
y x |
|
|
|
||||||||||||||||||||
y |
|
|
|
x |
|
|
|
|
|
|
|
2 |
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
y |
|
|
аналогично |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
zx |
|
|
|
u |
x |
u |
z |
|
|
2u |
x |
|
2u |
z . |
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
y |
|
|
|
|
2 |
|
|
|
|
||||||||||||||||
z |
|
|
|
|
z |
x |
|
|
|
|
|
z x |
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
z |
|
|
|
|
|
|
|
|
|
(6.1)
(6.2)
(6.3)
41
Суммируя (6.1)–(6.3) и группируя члены, получаем
|
p |
|
2 |
|
|
|
2 |
|
|
|
2 |
|
|
|
|
|
|
|
2 |
ux |
|
|
2 |
|
|
|
2 |
|
|
|||||||
|
|
ux |
|
|
|
|
ux |
|
|
|
ux |
|
|
|
|
|
|
|
uy |
|
uz |
|
. |
|||||||||||||
x |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
y x |
z x |
|||||||||||||
|
|
|
2 |
|
|
|
|
2 |
|
|
|
|
2 |
|
|
|
|
|
|
|
|
2 |
|
|
|
|
||||||||||
|
|
|
x |
|
|
|
y |
|
|
|
z |
|
|
|
|
|
|
|
x |
|
|
|
|
|
|
|
|
|
|
|
||||||
Третий член можно записать в виде |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||
|
|
|
|
|
u |
|
|
u |
y |
|
|
u |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
x |
|
|
|
|
|
z |
|
|
|
|
|
div |
u |
, |
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
x |
y |
|
|
|
z |
|
|
|
|
|
x |
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
x |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
но жидкость несжимаема и div u 0 . Таким образом, получаем
du |
x |
|
1 p |
|
|
2u |
x |
|
2u |
x |
|
2u |
x |
|
|
|
X |
|
|
|
|
|
|
|
|
|
. |
||||
|
|
x |
|
|
|
|
2 |
||||||||
dt |
|
|
2 |
2 |
|
||||||||||
|
|
|
|
|
|
|
x |
|
|
y |
|
|
z |
|
|
Выражение в скобках есть ничто иное, как оператор Лап-
ласа 2 |
|
, а |
|
. Окончательно получаем |
|
|
|
|
|
|
|
|
|
||||||||||||||||||
u |
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
u |
x |
|
|
|
u |
x |
|
|
|
u |
x |
|
|
u |
x |
|
1 p |
|
|
2u |
x |
|
2u |
x |
|
2u |
x |
|
||
|
|
u |
|
|
|
|
u |
|
|
u |
|
|
X |
|
|
|
|
|
|
|
|
|
. |
||||||||
|
|
|
|
|
|
x |
|
|
|
|
z |
x |
|
|
|
|
|
|
|
||||||||||||
|
t |
x |
|
|
|
y |
y |
z |
|
|
2 |
2 |
2 |
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
x |
|
|
y |
|
|
z |
|
|
Аналогично можно расписать и две другие проекции. Полученная система уравнений движения вязкой жидкости носит название системы уравнений Навье–Стокса.
В векторной форме можно записать
a |
|
|
1 |
grad p 2 |
|
. |
|
|
F |
(6.4) |
|||||||
u |
||||||||
|
||||||||
|
|
|
|
|
||||
42
Как видно уравнение (6.6) отличается от уравнения движения идеальной жидкости дополнительным членом ( 2u ),
учитывающим действие сил вязкого трения.
Целью гидродинамического расчета является нахождение полей скоростей и давлений, т. е. в результате расчета должны быть найдены четыре величины: ux , u y , uz и p. Принципи-
ально это оказывается возможным, так как три уравнения Навье–Стокса (в проекциях) плюс уравнение неразрывности образуют замкнутую систему. Плотность и вязкость, входящие в них, считаются известными, а проекции массовых сил (X, Y, Z) задаются условиями конкретной задачи.
С чисто математических позиций уравнения Навье–Стокса относятся к классу нелинейных дифференциальных уравнений в частных производных второго порядка. Одно из наиболее неприятных их свойств – нелинейность, обусловленная наличием конвективных членов ускорения. Следует отметить, что до настоящего времени, вследствие практически непреодолимых математических трудностей, не получено ни одного общего решения уравнений Навье–Стокса в их полном виде, т. е. при сохранении всех конвективных членов и всех членов, учитывающих вязкость. Известны лишь отдельные частные решения.
Одним из основных граничных условий при интегрировании является условие «прилипания», т. е. равенство нулю скорости жидкости на стенке.
43
СОДЕРЖАНИЕ
1. Движение деформируемой жидкой частицы. . . . . . . . . . 3 2. Кинематика вихревого движения. . . . . . . . . . . . . . . . . . . . 13 3. Потенциальное движение жидкости. . . . . . . . . . . . . . . . . 21 4. Преобразования Громеки–Лэмба. . . . . . . . . . . . . . . . . . . . 31 5. Модель вязкой жидкости. . . . . . . . . . . . . . . . . . . . . . . . . . 37 6. Уравнения движения вязкой жидкости
(уравнения Навье–Стокса) . . . . . . . . . . . . . . . . . . . . . . . . . 41
Учебное издание
КАЧАНОВ Игорь Владимирович КУЛЕБЯКИН Виталий Васильевич НЕДБАЛЬСКИЙ Викентий Константинович
МЕХАНИКА ЖИДКОСТИ И ГАЗА
Курс лекций
В4 частях
Ча с т ь 2
Редактор Т.Н. Микулик Компьютерная верстка Н.А. Школьниковой
Подписано в печать 02.02.2011. Формат 60 841/16. Бумага офсетная. Отпечатано на ризографе. Гарнитура Таймс.
Усл. печ. л. 2,56. Уч.-изд. л. 2,00. Тираж 100. Заказ 1159. Издатель и полиграфическое исполнение:
Белорусский национальный технический университет. ЛИ № 02330/0494349 от 16.03.2009.
Проспект Независимости, 65. 220013, Минск.
44
