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

21869_4f7d312f91bac43d4530a115cb0453bf

.pdf
Скачиваний:
4
Добавлен:
17.03.2016
Размер:
2.98 Mб
Скачать

С другой стороны jz E. Чтобы поставить знак равенства необходимо ввести коэффициент

jz = σq E. (1)

Этот коэффициент пропорциональности σq называется удельной электрической проводимостью вещества. Заметим аналогию

N = n<v>S/6 jz = Nq = q nq <v>др (N= nq<v>др). (2)

В выражении (2) отсутствует численный коэффициент 1/6, так как движение всего заряда упорядочено и направлено вдоль одной координаты. Скорость <v>др в данном случае называют дрейфовой. Это скорость дрейфа зарядов в слабом электрическом поле.

Запишем уравнение движения частицы согласно 2-го закона Ньютона

mdv/dt = qE, после интегрирования - v = qEt/m + v0.

Произведем усреднение. Тогда v0 = 0 , как скорость в начальный момент времени, t = τ - среднему времени свободного пробега, а v заменяется на <v>др

<v>др = qEτ/m.

Приравняем правые части (1) и (2) и подставим туда значение дрейфовой скорости

σq E = q nq <v>др = q2nqEτ / m.

Отсюда получим значение удельной проводимости, выраженное через микроскопические параметры

σq = q2nqτ/m.

Так как

τ = λ/<v>др = 1/2πdэфф2n<v>др , то

σq = q2nq / m2πdэфф2n<v>др,

где n - концентрация всех частиц, обусловленных столкновениями, а nq - средняя концентрация заряженных частиц. Они могут совпадать.

180

Глава 5 Гидродинамика

§1 Понятие о гидродинамике

1.1 Модель сплошной среды

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

рассматривается при этом как сплошная среда. Если в гидродинамике говорят о смещении некоторой частицы, то речь идет не о смещении отдельной молекулы, а о смещении целого элемента континуального объема. При этом, сколь маленький объем ни взять, в нем частичек предполагается актуально много (столько, сколько нам нужно), чем и хороша такая идеализация частичек, которые и сами в свою очередь состоят из воображаемых частичек всегда в достаточном количестве. Элемент такого объема можно рассматривать как точку, имеющую массу. Характеристиками жидкости здесь являются

скорость v = v(x,y,z,t), давление P = Р(x,y,z,t), плотность ρ = ρ(x,y,z,t).

Подчеркнем, что скорость (давление, плотность) рассматриваются в каждой данной точке пространства с координатами (x,y,z) в момент времени t и относится не к частицам жидкости, а к точкам пространства. Температуру можно считать как неизменной, так и любой удобной для данного рассмотрения.

Плотность ρ считаем неизменной вдоль всего объема жидкости и в течение всего времени движения. Иначе говоря, жидкость у нас несжимаемая. Сжи-

181

маемостью называется способность вещества (тела) изменять свой объем под действием всестороннего давления.

клапан

Сжимаемость = (V/P) - изотермическая, адиабатическая. Пользуются понятием сжимаемости в виде β = - V/VP, Па-1. Для примера в таблице приведены коэффициенты сжимаемости некоторых жидкостей

Вещество

t°C

P, атм

β, 10-6 атм-1

H2SO4

0

1-16

302,5

C2H5OH

20

1-50

112

Hg

20

1-10

3,91

1.2 Уравнение непрерывности

Пусть имеем объем V0, тогда масса жидкости в нем

m = ∫ ρdV. (1) V0

Пусть f - площадь поверхности, ограничивающей объем V0, а df - векторный элемент этой поверхности (с направлением вне (+) или во внутрь (-) ).

182

v

d f

Здесь v - длина в единицу времени, тогда

dV0 = v df, dm= ρdV0 = ρv df, m= ∫ ρ v df. (2) f

Здесь m- масса вытекающей из объема V0 жидкости или втекающей в него, а f - замкнутая поверхность, ограничивающая объем V0. Выражение (2) характеризует массу жидкости, заключенную в данном объеме.

Получим из (1) массу, меняющуюся со временем в данном объеме жидкости

dm/dt = d[∫ ρdV]/dt.

Тогда

d[ ∫ ρdV]/dt + ∫ ρ v df = 0

 

V0

 

f

 

На основании теоремы

Остроградского-Гаусса преобразуем интеграл по замк-

нутой поверхности в интеграл по объему

 

 

 

3

( g ds = div g dV), где div g = ∑ ∂g i /dxi, g i - компонента данного

s

V

 

i =1

вектора,

xi - компонента радиус вектора.)

 

 

∫ ρ v df = div (ρv) dv

 

f

V0

 

 

d[∫ ρdV] /dt + div (ρv) dV = [∂ρ/t + div (ρv)] dV = 0

V0

 

V0

V0

Так как V0 - произвольный объем, то и

∂ρ/t + div (ρv) = 0.

Получено так называемое уравнение непрерывности. В данном случае оно характеризует закон сохранения вещества. Часто используют ρv = j - плотность потока, [j] = кг/м2с. В нашем случае плотность потока массы жидкости (а может быть: газа, частиц и т.п.).

183

1.3 Об уравнении Эйлера

Имеем определение давления

P = dF/dS

продифференцируем обе части по оставшейся координате | d/dx

dP/dx = dF/dV = Fед. об..

Здесь Fед.об. - сила, приходящаяся на единичный объем, учтем ее векторный характер

Fед.об. = Fxi + Fyj + Fzk = (P/x) i + (P/y) j + (P/z) k = grad P.

В состоянии равновесия

Fед.об. = grad P.

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

ρdv/dt = Fед.об. - grad P,

где v - скорость единицы объема, а если учесть и силу веса жидкости, имеем:

ρdv/dt = Fед.об. - grad P + ρg.

Полученное выражение носит название уравнения Эйлера.

1.4 Теорема неразрывности струй

Картина тока жидкости представляется полем вектора скорости. Каждая линия тока является касательной к вектору скорости в данной точке. Густота линий тока пропорциональна величине скорости. Часть жидкости, ограниченная линиями тока называется трубкой жидкости. Вектора скорости не пересекают стенок трубки тока как касательные к ним.

184

S1 v1

S2 v2

Определение теоремы неразрывности струй:

S v = cst

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

Следовательно

S1v1 = S2v2 = ... = Sivi = ...

§ 2 Уравнение Бернулли

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

Рассмотрим трубку со стационарным течением. Пусть трубка ограничена стенками и сечениями S1 и S2 . За время t сечения переместятся на длины

l1 и l2.

185

l1

S1

l2

h1

S2 h2

Здесь h1 и h2 - высоты центров масс элементов объемов трубки тока над заданным уровнем. Согласно теореме о неразрывности струй

S1v1 = S2v2, S1l1/t1 = S2l2/t2,

но так как t1 = t2 S1l1 = S2l2 V1 = V2.

То есть элементарные объемы жидкости, образующиеся около сечений S1,S2 за один и тот же промежуток времени равны друг другу.

1. Рассчитаем кинетическую энергию элементарных объемов

K1 = m1v12/2 = V1 ρv12/2, K2 = V2ρv22/2

2. Рассчитаем потенциальную энергию относительно заданного уровня

W1 = m1gh1 = ρ∆V1gh1, W2 = ρ∆V2gh2.

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

E = E2 - E1 = V2ρv22/2 + ρ∆V2gh2 - V1ρv12 - ρ∆V1gh2.

Данное приращение энергии должно равняться совершаемой над объемом работе при его перемещении. Работа давления на боковые стенки равна 0 (здесь перемещения нет). Остается работа за счет разницы давлений на торцах.

A = P1S1l1 - P2S2l2 = P1V - P2V, E = A ρv22/2 + ρgh2 - ρv12/2 - ρgh1 = P1 - P2.

Произведено сокращение на элементарный объем, который в данном случае остается одинаковым.

186

ρv22/2 + ρgh2 + P2 = ρv12/2 + ρgh1 + P1.

Так как сечения выбирались произвольно, то

ρv2/2 + ρgh + P = cst. (1)

Это выражение будет точным при S 0. Определение

В стационарно текущей жидкости в отсутствие внутреннего трения вдоль любой линии тока справедливо уравнение (1), которое называется уравнением Бернулли.

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

При h1 = h2 (случай горизонтальных линий тока) уравнение имеет вид

ρv12/2 + P1 = ρv22/2 + P2.

Тогда, если v1 > v2, то P1<P2.

Пример: Водоструйный насос

жидкость

здесь

P<Pатмосф

здесь

образуется подсос

в атмосферу

Водоструйным насосом достигается уменьшение давления до 100 мм. рт.ст, что в 7,6 раза меньше атмосферного давления.

§ 3 Ламинарное и турбулентное течения

Рассмотрим течение жидкости по трубе. Ламинарное (слоистое):

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

Турбулентное:

187

жидкость энергично перемешивается во всех направлениях (Р. Осборн, английский физик 1842-1912).

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

Re = ρ v l /η, [η] = Па с = кг/м с.

ρ - плотность, v - средняя по сечению трубы скорость потока, η - коэффициент вязкости жидкости, l - характерный размер для трубы.

Начиная с некоторого критического значения числа Re, течение из ламинарного (для данной конфигурации сечения трубы) переходит в турбулентное. Для круглого сечения это Re = 1000. Характер течения различных жидкостей или (газов) будет одинаков для одинаковых чисел Рейнольдса.

Замечание. Иногда пользуются понятием кинематической вязкости

ν = η/ρ Re = vl/η, [ν] = м 2/с.

§ 4 Формула Пуазейля

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

R

r

 

 

 

l

 

(1)

(2)

Выделим в трубе с текущей жидкостью воображаемый цилиндр. Так как движение жидкости равномерное, то сумма внешних сил, приложенных к этому цилиндру должна быть равна нулю.

Рассчитаем силы, действующие на цилиндр.

1. Основания цилиндра. Составим разность сил давления на основания

(P1 - P2) πr2

188

2.Боковая поверхность. На боковую поверхность действуют силы внутреннего трения (вязкости).

F = K/t = η |dv/dr| Sбок = - η 2π rl dv/dr.

Знак минус означает убывание скорости с ростом r . Сложим силы и приравняем их к 0 (силы трения препятствуют перемещению, которому способствуют силы давления).

(P1 - P2) πr2 + η 2πrl dv/dr = 0

dv/dr = - (P1 - P2) r/2ηl.

Интегрируем полученное дифференциальное уравнение.

dv = - (P1 - P2)/2 ηl rdr,

v = C - (P1 - P2) r2/4ηl (*)

Найдем константу интегрирования C по граничному условию

v = (при r=R) = 0, то есть скорость движения у стенки жидкости равна 0, тогда

C = (P1 - P2)R2/4η l

(1).

Подставим (1) в (*)

и преобразуем.

v = (1 - r2/R2)(P1 - P2)R2/4η l.

То есть v r2 , следовательно, v(r) - парабола. В центре трубы при r = 0 скорость максимальна

vмакс = v(0) = (P1 - P2) R2/4ηl v( r ) = vмакс(1 - r2/R2). (2)

Вычислим поток жидкости (Q) - объем жидкости, протекающей за единицу времени через поперечное сечение трубы S.

189

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]