Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги / Общая физика. Гидродинамика и теплообмен

.pdf
Скачиваний:
2
Добавлен:
12.11.2023
Размер:
2.01 Mб
Скачать

Для прояснения ситуации рассмотрим следующий пример. Пусть необходимо измерить температуру воды в реке (ϕ = Т). Если мы измеряем температуру, стоя на мосту, то под Tt мы понимаем из-

менение Т со временем в конкретной точке A с координатами (xA, yA, zA). Если же мы проводим измерения, находясь на плоту, который плывет по течению со скоростью реки, то субстанциональная производная dTdt дает нам скорость изменения температуры какой-то

жидкой частицы, которая перемещается вместе с рекой по течению. Выражению (3.3) можно придать более компактный вид. Для этого рассмотрим скалярное произведение вектора скорости и опера-

тора «набла»:

r

r

r

r

 

r

 

r

 

r

 

 

 

v

= (vxi

+ vy j

+ vz k )

 

i

+

 

j

+

 

k

= vx

 

+ vy

 

+ vz

 

,

x

y

 

x

y

 

 

 

 

 

 

 

 

 

 

z

 

 

 

z

так как для ортов выполняются правила скалярного произведения:

r r

r r

r r

r r

r

r

r

r

= 0 .

i i = 1;

j

j

= 1; k

k = 1; i

j

= 0; i

k = 0;

j

k

Тогда субстанциональная или полная производная (3.3) примет

компактный вид

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dϕ

=

∂ϕ + (vr )ϕ .

 

 

 

(3.4)

 

 

 

 

 

 

 

 

 

 

 

 

dt

t

 

 

 

 

 

 

 

3.2. Уравнение энергии (теплопереноса)

 

 

Уравнением

энергии

является

дифференциальное

уравнение

в частных производных, описывающее температурное поле в движущейся среде. Вывод этого уравнения аналогичен выводу уравнения теплопроводности в твердом теле (см. подразд. 1.2). Отличие состоит в том, что теперь элементарный объем dV движется вместе с жидкостью со скоростью потока в той точке, где находится этот объем dV. Поэтому в полученном ранее уравнении (1.12) частную производную по времени Tt нужно заменить субстанциональной про-

изводной dTdt , в результате чего получим

41

dT = a T + qV . dt ρc

Если внутренние источники тепла отсутствуют, т.е. qV = 0 , то

 

T + (vr )T = a T ,

 

 

 

 

(3.5)

 

t

 

 

 

 

 

 

 

 

 

или в развернутом виде

 

 

 

 

 

 

 

 

 

 

T + vx

T + vy

T

+ vz

T

 

2

2

+

2

 

 

= a

T2 +

T2

T2

.

(3.6)

t

x

y

 

z

 

x

y

 

z

 

 

В уравнение (3.6) вошли три компоненты скорости vx , vy , vz , по-

этому необходимо получить уравнения движения, которые бы описывали изменение скорости во времени и в пространстве.

3.3. Уравнения движения

Вывод уравнений движения жидкости не является строгим, но зато он нагляден [2]. Выделим в потоке вязкой жидкости элементарный объем dV = dx dy dz . В основе вывода – второй закон Ньютона

F = ma : сила равна массе, умноженной на ускорение. Силы, действующие на dV, разделим на массовые или объемные (сила тяжести) и поверхностные (силы давления и трения).

Рассмотрим составляющую движения вдоль оси 0х. Сначала бу-

дем предполагать, что скорость

жидкости vx

= vx ( y,t)

изменяется

только в направлении оси 0y, потом это ограничение снимем.

Сила тяжести, приложенная к центру масс, вычисляется сле-

дующим образом:

 

 

 

 

df1 = ρgxdV = ρgxdxdydz .

 

(3.7)

Равнодействующая сил давления на верхнюю и нижнюю грани

(рис. 3.1)

 

 

 

 

 

p

 

p

(3.8)

df2 = pdydz p +

dx

dydz = −

dV.

 

x

 

x

 

42

 

 

 

 

Определим равнодействующую сил трения df3 . Согласно рис. 3.1

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

−τdxdz , (3.9)

вто время как около правой грани, наоборот, скорость движения жидкости больше, чем в элементе dV, поэтому

всечении y + dy сила трения направленав сторону движенияи равна

z dy

p y

dx

 

p + (dp /dx)dx

x

vx(y)

Рис. 3.1

 

 

 

 

τ +

∂τ dy dxdz .

(3.10)

 

 

 

 

 

 

y

 

 

По закону вязкого трения Ньютона

 

 

 

 

 

 

τ = µ

dvx

.

 

 

(3.11)

 

 

 

 

 

 

 

 

 

 

 

 

dy

 

Из (3.9)–(3.11) находим равнодействующую сил трения

df

 

=

dτ

dydxdz = µ

d2vx

dV .

(3.12)

3

 

 

 

 

dy

 

 

 

dy2

 

Такое сравнительно простое выражение получилось при предпо-

ложении, что vx = vx ( y,t) ,

однако в общем случае vx

меняется по

всем трем направлениям. В этом случае проекция на ось 0х равнодействующих сил трения

 

 

2vx

 

2vx

 

2vx

 

 

df3

= µ

 

 

+

 

 

+

 

 

dV = µ ∆vxdV .

(3.13)

x

2

y

2

z

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

43

Запишем теперь в проекции на ось 0х второй закон Ньютона для элемента dV:

ρdvx dV = df1 + df2 + df3 , dt

где слева стоит субстанциональная производная. Подставим выражения (3.7), (3.8) и (3.13) в последнее уравнение и сократим левую и правую части на dV:

vx

r

 

= ρgx

p

+ µ ∆vx .

(3.14)

ρ

t

+ (v )vx

x

 

 

 

 

 

 

 

В общем случае трехмерного движения аналогичные уравнения можно получить для vy и vz :

v

ρyt

ρvzt

 

 

= ρgy

p + µ ∆vy ,

 

+ (vr )vy

(3.15)

 

 

 

 

y

 

 

r

 

= ρgz

p

+ µ ∆vz .

(3.16)

+ (v

)vz

z

 

 

 

 

 

 

Выражения (3.14)–(3.16) представляют собой уравнения движения несжимаемой жидкости – уравнения Навье−Стокса. Можно объединить их, записав в векторной форме:

vr

+ (vr )vr = −

1

p + ν∆vr

+ gr.

(3.17)

t

ρ

 

 

 

 

Уравнение (3.17) есть основное уравнение гидродинамики вязкой несжимаемой жидкости в поле силы тяжести.

В уравнениях движения появилась еще одна неизвестная величина – давление р. Теперь число неизвестных (T, vx , vy , vz , p) пре-

вышает число уравнений ((3.6), (3.14), (3.15), (3.16)). Для замыкания системы необходимо добавить еще одно уравнение.

44

3.4. Уравнение неразрывности (сплошности)

Уравнение неразрывности выводится на основе закона сохранения массы. Выделим в потоке жидкости неподвижный элементарный объем в форме параллелепипеда dV = dxdydz (рис. 3.2) и посчитаем

массу жидкости, протекающую через него за время dt. В направлении оси 0х через

грань y0z втекает масса жидкости dM x = ρvxdtdydz .

Через противоположную грань вытекает масса

dM

 

v

 

 

(ρvx )

dx

dydzdt .

= ρ

x

+

 

x

 

x

 

 

 

 

 

 

 

Излишек массы жидкости, выте-

Рис. 3.2

кающей из dV в направлении оси 0х,

 

 

 

dM

dM

x

=

(ρvx ) dV dt .

 

x

 

x

 

 

 

 

 

 

 

Аналогичные выражения получаются для двух других направлений:

 

dM

dM

 

=

(ρvy )

dV dt ;

 

 

y

 

 

 

 

y

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

dM

dM

z

=

(ρvz ) dV dt .

 

 

 

z

 

 

z

 

 

 

 

 

 

 

 

 

 

 

 

Полный избыток массы жидкости, вытекающей из dV,

 

(ρvx )

 

 

(ρvy )

 

 

(ρvz )

r

dM ′ − dM =

 

+

 

 

 

 

 

+

 

 

dV dt = div(ρv)dV dt ,

x

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

z

 

равен уменьшению плотности жидкости в этом объеме за время dt:

∂ρ dV dt .

t

45

Приравняв оба выражения и поделив их на dVdt, получим уравнение неразрывности для сжимаемых жидкостей (газов):

∂ρ + div(ρvr) = 0 .

(3.18)

t

 

Длянесжимаемых жидкостей ρ = const, поэтому (3.18) упрощается:

 

 

 

divvr = 0

(3.19)

или в развернутом виде

vx

+

vy

+

vz

= 0 .

x

y

z

 

 

 

 

3.5.Краевые условия

Вслучае, когда движение жидкости является неустановившимся, необходимо задавать как начальные, так и граничные условия. На твердой поверхности выполняется так называемое условие прилипания жидкости

vr S = 0

или

vn S = 0; vτ S = 0 ,

где vn и vτ – нормальная и касательная составляющие скорости со-

ответственно.

Первое условие означает непроницаемость твердой стенки, второе – отсутствие проскальзывания жидкости вдоль стенки.

Измерения скоростей показали, что толщина неподвижного слоя жидкости весьма мала: она составляет несколько молекулярных слоев, т.е. не превышает 5 107 см.

В зависимости от рассматриваемой задачи возможны другие граничные условия: на границе двух несмешивающихся жидкостей, на свободной поверхности жидкости и т.д.

46

Контрольные вопросы

1.Дать определение конвективного теплообмена.

2.Дать определение субстанциональной производной.

3.Сравнить уравнение теплопереноса дляжидкостей и твердых тел.

4.Вывести уравнения Навье−Стокса.

5.Получить уравнение сплошности.

6.Привести краевые условия в задачах движения жидкости.

47

4. НЕКОТОРЫЕ ТОЧНЫЕ РЕШЕНИЯ УРАВНЕНИЙ ДВИЖЕНИЯ ВЯЗКОЙ НЕСЖИМАЕМОЙ ЖИДКОСТИ

Вследствие нелинейности системы уравнений движения жидкости точные решения удается найти лишь для небольшого числа сравнительно простых течений. Это возможно, если движение жидкости обладает симметрией, в результате которой дифференциальные уравнения в частных производных сводятся к обыкновенным дифференциальным уравнениям.

4.1. Установившееся движение жидкости между параллельными плоскостями – течение Куэтта

Рассмотрим изотермическое (T = const), установившееся ( t = 0 ) течение несжимаемой (ρ = const) жидкости между двумя параллельными горизонтальными плоскостями y = ±h (рис. 4.1). Течение одномер-

ное, т.е.

 

vy

= 0; vz

= 0 ,

а отлична от

нуля компонента скорости

vx = vx ( y) .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Из уравнения Навье−Стокса

 

 

 

 

 

 

 

 

 

 

 

(3.15) получаем

p y = 0 ,

откуда

 

 

 

 

 

 

 

 

 

 

 

следует,

что

p = const( y) . Анало-

 

 

 

 

 

 

 

 

 

 

 

гично

 

 

из

(3.16)

получаем

 

 

 

 

 

Рис. 4.1

 

 

 

p z = 0

 

и

p = const(z) . Значит,

 

 

 

 

 

 

 

 

давление зависит только от коор-

динаты х:

 

 

 

 

 

 

 

p = p(x) . Уравнение неразрывности (3.19) удовлетворяет-

ся автоматически.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Запишем уравнение (3.14) в развернутом виде:

 

 

 

 

 

 

vx

 

 

vx

 

 

vx

 

vx

 

1 p

 

µ

 

2vx

2vx

2vx

 

 

 

 

+ vx

 

+ vy

 

+ vz

 

= −

 

 

+

 

 

 

2 +

 

2 +

 

 

.

 

t

x

y

z

ρ ∂x

ρ

x

y

z

2

 

 

 

 

 

 

 

 

 

 

 

 

 

48

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Все слагаемые левой части уравнения равны нулю, справа отличны от нуля первое и третье слагаемое. Таким образом, остается

d2vx

=

1

 

dp

,

(4.1)

dy2

 

 

 

µ dx

 

где используются прямые производные, так как остались только

функции одной переменной: vx = vx ( y)

и p = p(x) .

 

Пусть задан постоянный напор жидкости:

 

 

dp

=

p = const ,

 

p

< 0 ,

(4.2)

 

dx

x

 

x

 

тогда уравнение (4.1) упрощается, так как остается только одна неизвестная – скорость потока vx :

d2vx

=

1

p .

(4.3)

 

 

dy2 µ ∆x

 

Запишем граничные условия. Предположим, что верхняя плоскость движется со скоростью v1 , а нижняя – со скоростью v2 , тогда условия прилипания надвижущихсятвердых границах будут иметьвид

vx

y=h = v1; vx

y=− h = v2 .

(4.4)

Обыкновенное дифференциальное уравнение (4.3) нетрудно проинтегрировать. Его общее решение имеет вид

vx

=

1

p

y2

+ C1 y + C2 ,

(4.5)

 

 

 

 

µ ∆x 2

 

 

где С1 и С2 – константы интегрирования, значения которых находим из граничных условий (4.4):

v1

=

1

p

h2

+ C1h + C2 ; v2

=

1

p

h2

C1h + C2 .

 

 

 

 

 

 

µ ∆x 2

 

 

µ ∆x 2

 

49

Решив два последних уравнения, получим

C1

=

v1 v2

;

C2

=

v1 + v2

1

p

h2

.

 

 

 

 

 

 

2h

 

2

 

µ ∆x 2

Подставим значения констант интегрирования в общее решение (4.5) и получим решение поставленной задачи:

vx

= −

(h2 y2 )

p +

v1 v2

y +

v1 + v2

.

(4.6)

 

2h

 

 

 

2µ ∆x

2

 

 

Рассмотрим несколько частных случаев.

1. Если движение безнапорное, т.е. px = 0 , то распределение скорости по сечению потока будет линейным:

vx =

v1 v2

 

y +

v1 + v2

.

(4.7)

2h

 

2

 

 

 

 

 

 

2. Если в безнапорном потоке нижняя плоскость неподвижна

( v2 = 0 ), то из (4.7) далее получим

 

 

 

 

 

 

 

 

 

 

 

v1

 

y

 

vx =

 

 

1+

 

 

 

.

(4.8)

2

 

 

 

 

 

 

 

 

 

h

 

3. Если в напорном потоке обе стенки неподвижны ( v1 = 0; v2

= 0 ),

тоиз решения (4.6) следует, что

 

 

 

 

 

 

 

 

 

 

vx = −

(h2 y2 )

p .

(4.9)

 

 

 

 

 

 

2µ

 

 

 

x

 

Последнее уравнение описывает так называемое течение Пуазейля между неподвижными параллельными плоскостями. Максимальная скорость такого потока жидкости будет наблюдаться при y = 0:

(vx )max = −

h2

p .

(4.10)

2µ ∆x

На рис. 4.2 показан параболи-

ческий профиль скорости течения Рис. 4.2 Пуазейля.

50