Дифференциальные уравнения движения жидкости

Дифференциальные уравнения движения жидкости - student2.ru Дифференциальные уравнения движения жидкости - student2.ru Дифференциальные уравнения движения жидкости - student2.ru Это уравнение векторное и в проекциях на оси выбранной системы координат дает три уравнения. Вывод дифференциального уравнения движения в общем случае требует громоздких математических выкладок. Поэтому для упрощения вывода рассмотрим одномерное течение жидкости. Выделим в потоке вязкой жидкости, как показано на рис. 2.3, элементарный объем с размерами ребер dx, dy и dz.

 
  Дифференциальные уравнения движения жидкости - student2.ru

Рис. 2.3

Скорость в потоке изменяется только в направлении оси Оу. Силы, действующие на элемент жидкости, можно разделить на массовые (или объемные) и поверхностные. К массовым силам относятся сила тяжести, центробежная сила и электромагнитные силы. Мы в дальнейшем будем учитывать только силу тяжести. К поверхностным силам относятся силы давления и силы трения. Найдем проекции этих сил на ось Ох.

Сила тяжести равна произведению массы элемента на проекцию ускорения свободного падения на ось Ох

Дифференциальные уравнения движения жидкости - student2.ru

Равнодействующая сил давления в проекции на ось Ох равна

Дифференциальные уравнения движения жидкости - student2.ru

Равнодействующая сил трения с учетом уравнения (1.11) в проекции на ось Ох составляет

Дифференциальные уравнения движения жидкости - student2.ru

Суммируя Дифференциальные уравнения движения жидкости - student2.ru получим проекцию равнодействующей всех сил на ось Ох, приложенных к объему

Дифференциальные уравнения движения жидкости - student2.ru (а)

Согласно второму закону механики эта равнодействующая равна произведению массы элемента на проекцию ее ускорения на ось Ох и учитывает силы инерции

Дифференциальные уравнения движения жидкости - student2.ru (б)

Приравнивая правые части (а) и (б) и производя сокращения, получим уравнение движения вдоль оси Ох

Дифференциальные уравнения движения жидкости - student2.ru

В общем случае трехмерного движения несжимаемой жидкости уравнение движения в проекциях на оси Ox, Oy и Oz соответственно имеет вид:

Дифференциальные уравнения движения жидкости - student2.ru (2.5)

Дифференциальные уравнения движения жидкости - student2.ru (2.6)

Дифференциальные уравнения движения жидкости - student2.ru (2.7)

Или в векторном виде

Дифференциальные уравнения движения жидкости - student2.ru (2.8)

Уравнения (2.5) – (2.8) называют уравнениями Навье-Стокса. Член, стоящий в левой части уравнений, представляет собой полную производную от скорости по времени

Дифференциальные уравнения движения жидкости - student2.ru (2.9)

где i соответственно x, y, z.

Первые слагаемые в правой части (2.9) характеризуют локальное изменение скорости во времени в какой-либо точке жидкости. Остальные три слагаемых в правой части характеризуют изменение скорости при переходе от точки к точке. Такая полная производная называется субстанциональной производной. Уравнения движения получены при постоянных теплофизических свойствах жидкости. В то же время свободное движение жидкости (естественная конвекция) определяется разностью плотностей холодных и нагретых частиц жидкости. В общем случае при r¹const необходимо учитывать и энергию деформации жидкости. Поэтому ограничимся приближенным учетом переменности плотности в слагаемом, связанным с силой тяжести в уравнениях движения. Пусть плотность линейно зависит от температуры

Дифференциальные уравнения движения жидкости - student2.ru

где r и r0 – плотности, соответствующие температурам t и t0;

J=t-t0; t0 – некоторая фиксированная температура (точка отсчета).

Подставляя это значение плотности в первое слагаемое правой части (2.8), получим приближенно

Дифференциальные уравнения движения жидкости - student2.ru

Первое слагаемое правой части Дифференциальные уравнения движения жидкости - student2.ru можно трактовать как сумму силы тяжести Дифференциальные уравнения движения жидкости - student2.ru , взятой при определенной плотности, и подъемной (архимедовой) силы Дифференциальные уравнения движения жидкости - student2.ru . Член Дифференциальные уравнения движения жидкости - student2.ru можно представить как градиент гидростатического давления р0 в покоящейся жидкости с плотностью r0. Тогда вместо Дифференциальные уравнения движения жидкости - student2.ru можно написать grad p1, где p1=p-p0. При такой замене приближенное векторное уравнение движения будет описывать и естественное движение жидкости (естественную конвекцию)

Дифференциальные уравнения движения жидкости - student2.ru (2.10)

Таким образом, для задач теплообмена система дифференциальных уравнений сохранения массы (уравнение неразрывности или сплошности), энергии и движения в проекциях на координатные оси оказывается замкнутой. Эта система уравнений в принципе позволяет определить в движущейся жидкости поле температуры T=T(x, y, z, t), поле давлений p=p(x, y, z, t) и поля проекций скоростей wx=wx(x, y, z, t), wy=wy(x, y, z, t), wz=wz(x, y, z, t).

Для задач массообмена, не осложненных теплообменом (в изотермических условиях), уравнение энергии в этой системе заменяется уравнением сохранения массы i-го компонента смеси. Вывод этого уравнения, которое называют уравнением массообмена, аналогичен выводу дифференциального уравнения сохранения энергии при qV=0 и имеет вид

Дифференциальные уравнения движения жидкости - student2.ru (2.11)

где mi=Ci/r – относительная массовая концентрация i-го компонента.

В случаях, когда массообмен осложнен теплообменом (в неизотермических условиях), кроме уравнения (2.11), необходимо уравнение сохранения энергии. Однако вывод этого уравнения с учетом (1.13) усложняется. Для двух компонентной (бинарной) смеси оно имеет вид

Дифференциальные уравнения движения жидкости - student2.ru

Дифференциальные уравнения движения жидкости - student2.ru Дифференциальные уравнения движения жидкости - student2.ru (2.12)

Из этого уравнения видно, что, если удельные изобарные теплоемкости компонентов смеси равны ср1р2, то результирующий перенос энтальпии отсутствует, и это уравнение переходит в ранее полученное уравнение (2.4).

Наши рекомендации