Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Лекции гидродинамика Тф-14 / ВВЕДЕНИЕ, ГЛАВЫ 1- 2.ppt
Скачиваний:
181
Добавлен:
26.05.2014
Размер:
3.79 Mб
Скачать

Уравнения идеальной жидкости содержат пять скалярных функций:

три компоненты скорости;

плотность;

давление.

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

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

vn 0, если граница неподвижна;

vn vГ , если скорость границы равна vГ

Если имеются две несмешивающиеся жидкости (1 и 2), то на поверхности их раздела равны их нормальные компоненты скорости и выполняется равенство давлений:

v n1 v n2

p1 p2

2.3. Гидростатика

Если жидкость покоится, то из уравнения Эйлера следует, что

 

p F

 

 

В однородном поле тяжести F

 

 

g Откуда:

F

p g

 

 

 

Последнее уравнение описывает механическое равновесие жидкости в поле сил тяжести. Если внешние силы вообще отсутствуют, то

p 0

 

p const

 

Задача 1: показать, что при const давление в жидкости есть

p gz const

Задача 2: определить форму поверхности жидкости в поле тяжести в цилиндрическом сосуде, вращающемся вокруг своей оси с постоянной угловой скоростью

2.4. Стационарные течения, уравнение Бернулли

Важным случаем гидродинамических течений являются стационарные

(установившиеся) течения. В таких течениях в

каждой точке

пространства

скорость течения остается постоянной во времени.

В этом случае

v/ t 0

Тогда:

 

p

 

(v )v -

 

Воспользуемся теперь изоэнтропийностью течений. В этом случае можно записать для тепловой функции (энтальпии):

dw Tds Vdp Vdp (ds 0)

Отсюда:

dw Vdp (1/ )dp

(1/ ) p w

 

Поэтому:

(v )v - w

 

Используя соотношение векторного анализа:

(1/2) (v2 ) v, rotv (v )v

Получаем:

v, rotv (w v2 /2)

 

Откуда:

(1/2) (v2 )

v, rotv w

 

Введем понятие о линиях тока:

Даниил

Бернулли

 

 

(1667-1748)

v

Линии тока – линии, касательные к которым указывают направление вектора скорости в точке касания в данный момент времени – они определятся системой дифференциальных уравнений:

dx dy dz

vx

vy

vz

Если уравнение умножить на единичный вектор касательный к линии тока в каждой точке, то, поскольку градиент на некоторое направление равен производной, взятой по этому направлению, то (проекция v, rotv равна нулю!)

( / l ) v2 / 2 w 0

Отсюда следует, что

v2 w const 2

Последнее уравнение носит название уравнения Бернулли (1738г.)

Если течение жидкости происходит в поле сил тяжести (сила тяжести действует вдоль направления z, то уравнение Бернулли будет иметь вид (задача: показать !):

Даниил

Бернулли (1667-1748)

v2 w gz const 2

2.5. Поток энергии Вычислим изменение со временем энергию единичного объема жидкости

2

 

 

 

- сумма кинетической и внутренней энергии жидкости

 

 

 

 

 

 

 

 

Изменение со временем энергии есть:

 

 

 

 

 

 

r

 

 

 

 

 

 

 

 

 

 

v

2

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

v

2

 

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

v

 

 

 

 

 

 

 

 

 

 

 

 

 

v

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

v

 

 

 

 

 

 

 

 

t 2

 

t

2

 

 

t

 

 

2

t

t

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Используя уравнения неразрывности и движения, получаем:

 

 

 

 

 

 

 

 

 

 

 

 

v2

 

v2

 

 

r r

 

 

 

 

 

r r r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

div( v) - v

- v(v )v

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

2

2

 

 

 

 

 

 

 

 

 

Изменение внутренней энергии со временем можно найти из

 

термодинамического соотношения:

 

 

 

 

 

 

 

 

p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d Tds pdV Tds

 

 

d

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Но энтальпия (тепловая функция) есть

w = e+ p/ r = pV

 

Откуда:

 

 

d ( ) d d wd Tds

Отсюда можно записать (используя соотношение адиабатичности ds/dt=0):

( )

 

 

 

s

r

r

t

w

t

T

t

wdiv( v) Tv s

Окончательно можно записать:

 

 

 

 

 

 

 

2

 

 

v

2

r

r

 

v

2

 

v

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(w

 

 

)div( v) (v )

w

 

 

t

2

 

 

 

2

 

 

 

2

 

 

 

 

 

 

 

Проинтегрируем полученное выражение по объему:

 

 

r

 

 

div

v

 

 

 

 

 

 

 

2

 

 

 

 

 

v

 

 

 

 

w

 

 

2

 

 

 

 

 

 

 

 

 

v2t 2

V

q

 

 

 

 

 

 

2

 

 

 

 

 

 

 

2

r

 

 

 

r

 

 

 

 

r

v

 

 

 

 

 

v

 

 

 

Ñ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

w

 

v

 

 

 

 

w dS

dV

 

div v

2

dV

 

 

2

 

 

V

 

 

 

 

S

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

r v2

 

 

q v

 

 

w

- вектор плотности потока энергии

 

2

 

 

 

 

Абсолютная величина вектора плотности потока энергии есть количество энергии, протекающей в единицу времени через единичную поверхность, расположенную перпендикулярно к направлению скорости

Информация к размышлению:

Вектор потока энергии показывает, что каждая частица жидкости переносит при своем движении энергию w+v2/2. Факт, что переносится не величина внутренней энергии, а тепловая функция связано со следующим обстоятельством. Если подставить w p /

то полный поток энергии через замкнутую поверхность есть:

 

 

r

Ñ

r

v2

 

r

 

Ñ

r

r

 

 

 

 

 

 

 

 

Q

v

 

 

 

dS

 

pvdS

2

 

S

 

 

 

 

S

 

 

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

2.6. Поток импульса

Вычислим изменение со временем импульса p v единичного объема

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

 

 

 

 

 

 

v

 

 

vi

v

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

i

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

t

 

 

 

 

 

 

( vk )

 

 

 

 

 

Запишем уравнение неразрывности в форме

 

 

 

 

 

 

 

 

 

t

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

vk x

 

 

x

 

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Уравнение Эйлера в форме

 

 

vi

 

 

 

 

vi

 

 

 

1 p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k

 

 

 

 

 

 

 

 

i

 

 

 

 

 

 

 

 

Тогда получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

v

 

v

vi

p

v

( vk )

p

 

 

 

 

( v

v

 

)

 

 

 

 

 

x

 

x

 

 

 

t

i

 

 

k x

 

 

x

 

i

x

k

 

 

 

 

 

 

 

i

 

k

 

 

 

 

 

 

 

 

k

 

 

 

 

i

 

 

 

 

 

 

 

 

 

 

 

i

 

 

 

k

 

 

 

Можно записать

 

 

 

 

 

 

p

 

 

 

 

p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

ik x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

 

 

 

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Тогда получим

 

 

vi

ik

 

 

t

xk

 

 

 

Тензор Пik определяется так:

 

 

ik

p ik vi vk

Тензор Пik является симметричным, т.е.

ik p ik vi vk ki p ki vk vi

Чтобы выяснить физический смысл тензора Пik, проинтегрируем по некоторому объему:

 

 

vi dV

 

ik

dV

ik dSk

t

 

 

 

x

Ñ

 

V

 

V

k

 

Слева стоит изменение в единицу времени i –компоненты импульса в объеме, а справа – количество импульса, протекающего в единицу времени через поверхность