Лабораторная работа № 15
ОПРЕДЕЛЕНИЕ КОЭФФИЦИЕНТА ДИНАМИЧЕСКОЙ ВЯЗКОСТИ ВОЗДУХА ПРИ РАЗЛИЧНЫХ ТЕМПЕРАТУРАХ
Цель работы изучение явления внутреннего трения (вязкости) в газах, экспериментальное определение коэффициента динамической вяз-
кости воздуха методом истечения через капилляр (метод Пуазейля), изу- |
|||
чение зависимости коэффициента динамической вязкости от температу- |
|||
О |
|
|
|
ры. |
|
|
|
|
1. Метод измерения и расчетные соотношения |
||
Ф |
|
Вязкость (или внутреннее трение) |
|
и |
|
— свойство газов и жидкостей, про- |
|
|
являющееся в возникновении сил со- |
||
|
Я |
противления (трения) между их от- |
|
|
дельными слоями. При ламинарном |
||
|
С |
(слоистом) течении жидкости или |
|
|
|
|
|
|
|
|
газа в трубе скорости u различных |
|
|
|
слоев образуют определенный про- |
|
|
Н |
|
|
|
|
филь по сечению трубы (рис. 1). |
|
|
|
аибольшей скоростью обладает осе- |
Рис. 1. Профиль скорости жидкости |
вой слой жидкости. Скорость слоев, |
||
(газа) при ламинарном течении в |
соприкасающихся со стенкой канала |
||
|
трубе |
Иравна нулю (пристеночный слой |
|
|
|
|
«прилипает» к стенке). |
За счет разной скорости отдельных слоевУжидкости или газа между ни- |
|||
|
|
|
М |
ми возникают силы трения. При этом быстро движущиеся слои тормозятся, а медленно движущиеся — ускоряются. Постоянство профиля скоро-
сти поддерживается постоянным перепадом давления (р1 – р2), который |
||
обеспечивается внешними условиями. |
Э |
|
|
|
|
Экспериментально установлено, что сила трения F, действующая меж- |
||
|
|
И |
ду отдельными слоями жидкости, направлена по касательной к скорости |
||
слоев и пропорциональна модулю градиента скорости du |
, характеризую- |
|
|
dr |
|
щему быстроту изменения скорости в направлении r перпендикулярном скорости слоев (1):
F η |
|
du |
|
S , |
(1) |
|
|
||||
|
|
dr |
|
|
|
где S – площадь поверхности между взаимодействующими слоями;— размерный коэффициент пропорциональности, называемый коэффициентом динамической вязкости (внутреннего трения) жидкости или газа. Единицей динамической вязкости является Па∙с.
Выражение (1) известно как закон вязкого трения Ньютона.
С точки зрения молекулярно-кинетической теории газов вязкость определяется переносом импульса молекул из слоя в слой в направлении, перпендикулярном скорости упорядоченного движения газа. В результате теплового (хаотического) движения молекулы переходят из слоя газа с большей скоростью u в соседний слой с меньшей скоростью и ускоряют этот слой. Наоборот, молекулы из «медленного» слоя переходят в более быстрый соседний слой и замедляют его.
Молекулярно-кинетическая теория идеального газа позволяет получить
выражение для коэффициента динамической вязкости |
|
||||||||||
О |
η 1 |
ρ |
|
|
|
, |
|
(2) |
|||
v |
l |
|
|||||||||
где — плотность газа; v |
3 |
|
|
|
|
|
|
|
|
||
— средняя скорость теплового движения моле- |
|||||||||||
Ф |
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
кул; l — средняя длина свободного пробега молекул газа. |
|
||||||||||
Средняя скорость теплового движения молекул идеального газа опре- |
|||||||||||
деляетсявыражением |
|
|
|
|
|
|
|
|
|
||
|
|
Я v |
|
8RT |
, |
(3) |
|||||
|
|
|
|
||||||||
|
|
|
|
|
|
πμ |
|
||||
где R — универсальная газовая постоянная, Т — абсолютная температура, |
|||||||||||
|
|
С |
|
|
|
|
|
|
|
|
|
— молярная масса газа. |
|
|
|
|
|
|
|
|
|
Поскольку средняя длина свободного пробега молекул идеального
газа не зависит от температуры, то, в соответствии с (2) и (3) можно счи- |
|
Н |
|
тать, что коэффициент динамической вязкости такого газа с ростом тем- |
|
пературы увеличивается по закону |
|
И |
(4) |
η const T |
|
У |
|
Анализируя течение газа через канал малого диаметра (капилляр),
можно найти соотношение между коэффициентом динамической вязко-
сти исследуемого газа, объемным расходом, перепадом давления на кон-
цах капилляра и его геометрическими характеристиками (диаметром и |
|
длиной). Это соотношение является основойМметода капилляра (метод |
|
|
Э |
Пуазейля), используемого для измерения коэффициента динамической |
|
вязкости. |
И |
В дальнейшем изложении для кратности будем говорить о течении в канале газа, хотя метод Пуазейля применим для исследования вязкости как жидкой, так и газообразной фаз.
Рис.2. К выводу формулы Пуазейля
Рассмотрим ламинарное течение газа через цилиндрический канал малого диаметра (капилляр). Пусть внутренний радиус капилляра r0. Наличие внутреннего трения приводит к возникновению градиента скорости упорядоченного движения. Как уже отмечено выше, для обеспечения стационарного течения газа необходимо поддерживать постоянным во вре-
мени перепад давления (р1 – р2) между входом в капилляр и выходом из |
||||||||||||||||||
него. Для определения зависимости u u (r ) мысленно выделим в газе |
||||||||||||||||||
цилиндрический объем радиуса r и длины L (рис. 2). На торцах этого ци- |
||||||||||||||||||
О |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
линдра поддерживается разность давлений p1 – p2 |
, что соответствует |
|
||||||||||||||||
результирующей силе |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
Ф |
F p p |
2 |
r 2 . |
|
|
|
(5) |
|||||||||||
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
На боковую поверхность цилиндра со стороны окружающего слоя газа |
||||||||||||||||||
и |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
действует с |
ла внутреннего трения, которая согласно закону вязкого тре- |
|||||||||||||||||
ния Ньютона (1) равна |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
Я |
|
du |
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
F(r) |
η |
2πrL , |
|
|
(6) |
|||||||
где 2πrL |
|
|
|
|
С |
|
dr |
|
|
|
|
|
|
|
|
|||
|
— площадь боковой поверхности выделенного цилиндра. Так |
|||||||||||||||||
как газ движется равномерно (стационарное течение), то векторная сумма |
||||||||||||||||||
сил, действующих на выделенный цилиндр, и сумма их проекций на ось z |
||||||||||||||||||
равна нулю (рис.2): |
|
|
|
|
|
du |
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
p |
|
И |
. |
(7) |
|||||||||
|
|
|
|
|
Нp πr η |
2πrL 0 |
||||||||||||
|
|
|
1 |
2 |
|
|
|
|
dr |
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Из (7) с учетом убывания скорости u при удалении от оси симметрии z |
||||||||||||||||||
du |
|
du |
|
|
du |
|
|
|
|
|
У |
|
||||||
|
|
|
|
|
|
|
|
|||||||||||
( dr 0 и |
|
dr |
|
|
dr ) имеем |
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
p1 p2 |
|
|
|
М |
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
du |
|
rdr. |
|
(8) |
||||||||
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
2ηL |
|
|
|
|
|
Э |
|
||
Интегрирование выражения (8) дает |
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
И |
||||||||||||
|
|
|
|
|
|
u |
p1 p2 |
r |
2 |
С. |
||||||||
|
|
|
|
|
|
4ηL |
|
|
|
Постоянную C определяем из граничного условия: при r = r0 u = 0, где r0 — внутренний радиус капилляра. Отсюда
Сp14ηLp2 r02
иокончательно для скорости получаем:
u |
p1 p2 |
r02 |
r2 . |
(9) |
|
4ηL |
|||||
|
|
|
|
Из (9) следует, что при ламинарном течении жидкости или газа профиль скорости, изображенный на рис. 1 и рис. 2 носит параболический характер.
С целью определения объема V газа, протекающего через поперечное
сечение капилляра за время , выделим в газе цилиндрический слой с внутренним радиусом r и внешним r dr (рис. 2), считая во всех точках Оэтого слоя скорость газа u одной и той же. За время из выделенного
нами слоя вытечет объем газа dV uτ2πrdr . Используя (9) и интегрируя
это выражение по всей площади сечения капилляра (при этом r изменяет-
ся от 0 до r0), получим выражение для объема газа вытекающего из ка- |
|||||||||||||
|
и |
|
|
|
|
|
|
|
|
|
|
||
пилляраФза время |
|
πr 4 |
|
|
|
|
|
|
|
||||
|
|
Я |
V |
p p |
|
τ. |
|
||||||
|
|
|
0 |
|
|
||||||||
|
|
|
|
|
|
||||||||
|
|
|
8ηL |
|
1 |
|
2 |
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
||
Очевидно, что за единицу времени через капилляр проходит объем |
|||||||||||||
газа |
|
|
|
|
|
πr 4 |
|
|
|
|
|
|
|
|
|
|
|
G |
|
|
p p , |
(10) |
|||||
|
|
|
|
|
0 |
|
|||||||
|
|
|
СV |
|
|
|
1 |
|
|
2 |
|
||
|
V |
|
|
8ηL |
|
|
|
||||||
где G |
τ |
называется объемным расходом жидкости или газа, который |
|||||||||||
V |
|
|
И |
|
|
||||||||
может быть достаточно точноНизмерен в эксперименте. Соотношение (10) |
|||||||||||||
известно как формула Пуазейля. |
|
|
|
|
|
|
|
|
|
||||
В соответствии с (10), зная геометрию капилляра, измерив в опыте |
|||||||||||||
|
|
|
|
|
|
|
У |
|
перепад давлений на входе и выходе из него, а также измерив объемный расход, можно экспериментально определить коэффициент динамической вязкости жидкости или газа
|
πr4 |
p p |
2 |
|
|
М |
|
|
η |
0 |
1 |
|
. |
(11) |
|||
|
8G L |
|
|
|||||
|
|
V |
|
|
|
|
Э |
|
|
|
|
|
|
|
|
|
|
2. Описание схемы установки |
И |
|||||||
|
|
|
|
|
|
|
||
Схема экспериментальной установки представлена на рис. 3. |
|
Микрокомпрессор 1 создает необходимое давление и прокачивает воздух сначала через измеритель объемного расхода газа (ротаметр) 2, а затем по магистральной трубке – в медный термостат 3, внутри которого находится измерительный участок – капилляр 4. Из капилляра воздух сбрасывается в атмосферу. Термостат предназначен для нагрева и поддержания необходимой температуры воздуха при его входе в капилляр, и тем самым позволяет проводить измерения вязкости воздуха при различных температурах. Температура воздуха измеряется термопарой 5, спай которой находится при входе воздуха в капилляр и с помощью термоэлектродов 7 подключена к показывающему прибору – измерителю температуры
8. Медный термостат 3 помещен в керамическую трубку, на внешней поверхности которой смонтирован электрический нагреватель 9, питаемый от регулятора мощности 10 с ручками управления 15 и 16, расположенными на лицевой панели (рис. 4). Перепад давления на капилляре измеряется U-образным водяным манометром 6 по разности высот столбов жидкости в левом и правом коленах. Расчет перепада давлений проводится по формуле
новки в целом. |
Н |
|
|
Расположение и назначение тумблеров и ручек управления данной |
|||
установкой приводится в инструкции, размещенной на установке. |
|
||
|
И |
|
|
|
p1 p2 |
У |
|
|
М |
|
|
|
Э |
|
|
|
И |
||
|
ρвg h1 h2 , |
(12) |
|
где в — плотность |
воды при |
комнатной температуре; g — ускорение |
|
свободного падения; |
h1 h2 — |
разность уровней воды в трубках мано- |
|
О |
|
|
|
метра. |
|
|
|
Экспер ментальная установка смонтирована в типовом модуле (рис. 4) |
|||
наФлицевой |
которого выведены тумблеры и ручки включения, |
||
управления, а также необходимые измерительные приборы. На лицевую |
панели панель выведеныЯ: трубка ротаметра 2, предназначенного для измерения
объемного расхода воздуха; трубки U-образного манометра 6 с миллиметровой шкалой; лицевая панель измерителя температуры 8, который включается тумблером 13; тумблеры 12, 14 и 11, служащие соответсвенно
для включения иСвыключения микрокомпрессора, вентиллятора и уста-
Рис. 3. Схема экспериментальной |
|
установки |
Рис. 4. Лицевая панель модуля |
|
|
|
установки |
3.Порядок выполнения работы
1.Заполните табл.1 спецификации измерительных приборов. Запишите данные установки, которые приведены на лицевой панели. Подключите установку к электрической сети.
|
|
|
|
|
|
|
|
|
|
|
|
Таблица 1 |
|||
|
|
|
|
Спецификация измерительных приборов |
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
Название прибора |
Пределы |
Цена |
Инструментальная |
|||||||||||
О |
и его тип |
измерения |
деления |
погрешность |
|
||||||||||
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Ф |
|
|
|
|
|
|
|
|
|
|
|
||||
|
и |
Данные установки: |
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|||||
Длина кап лляра |
L = |
мм, |
L = |
мм. |
|
|
|
|
|
|
|||||
Диаметр капилляраЯd = мм, |
d = |
мм. |
|
|
|
|
|
|
|||||||
|
|
|
|
|
С |
|
|
|
|
|
|
|
|||
2. Включив компрессор, проведите измерения расхода воздуха через |
|||||||||||||||
капилляр в делениях шкалы ротаметра и высоты столбов воды h1 и h2 в |
|||||||||||||||
|
|
|
|
|
|
|
Н |
|
|
|
|
|
|
||
коленах манометра при комнатной температуре. Точное значение темпе- |
|||||||||||||||
ратуры воздуха определите по показаниям измерителя 8. Результаты из- |
|||||||||||||||
мерений запишите в табл. 2. |
И |
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
У |
||||||
|
|
|
|
|
|
|
Результаты измерений |
Таблица 2 |
|||||||
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
М |
|||||
|
|
|
№ |
|
t, oС |
h1, |
h2, |
(p1 – p2), |
|
GV |
|
||||
|
|
|
3 |
|
|
|
мм |
мм |
Па |
|
|
дел. |
м3/с |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Э |
|||
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
И |
||
|
|
|
4 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
5 |
|
|
|
|
|
|
|
|
|
|
|
|
3.Включите блок нагревателя 9 и установите регулятор температуры
вположение 2 (примерно 100o С). Через 10-15 минут убедитесь в том, что режим течения и температурный режим стабилизировались. Измерьте температуру на входе в капилляр, объемный расход воздуха (в делениях) и перепад давления на капилляре (h1 и h2). Результаты запишите в табл. 2.
4.Повторите действия, описанные в п. 2 для последующих темпера-
турных режимов (не менее трех).
Замечание. Максимальный температурный режим определяется конструктивными особенностями конкретной установки.
5. После окончания измерений выключите блок нагревателя и включите вентилятор для охлаждения установки.
6. При достижении комнатной температуры по показаниям измерителя 8 (рис. 4) отключите вентилятор, компрессор и установку в целом.
4. Обработка результатов измерений
О1. Рассчитайте перепад давлений воздуха на капилляре для всех температурных режимов течения, используя формулу (12). Результаты расчетаФзанесите в табл.2.
2. Используя паспортные данные ротаметра, установленного на конкретнойиустановке, произведите перерасчет показаний ротаметра в делениях в объемный расход, выраженный в м3/с.
3. По формулеЯ(11) рассчитайте значения коэффициента динамической вязкости воздуха для всех температурных режимов. Результаты расчета занесите в табл.3.
ную зависимость (Т) с соотношением (4). Сделайте выводы.
4. ПостройтеСграфик зависимости (Т). Сравните эксперименталь-
5. Для одного из температурных режимов рассчитайте погрешность |
|||||||||||||||||||||||
|
|
|
|
|
|
Н |
|
|
|
|
|
|
|
|
|
|
|
||||||
косвенного измерения . Запишите окончательный результат в стандарт- |
|||||||||||||||||||||||
ной форме. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
η = |
|
16 |
|
r0 |
2 |
2 |
h1 2 |
|
ρ 2 |
|
GV 2 |
|
L 2 |
|||||||||
η |
|||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
r |
|
|
|
|
|
У |
|
|
|
L |
|||||||
|
|
|
|
|
|
|
Иh ρ |
|
G |
|
|||||||||||||
|
|
|
|
|
0 |
|
|
1 |
|
|
|
|
V |
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Таблица 3 |
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
М |
|
|
||||||
Экспериментальные значения коэффициента динамической вязкости при |
|||||||||||||||||||||||
|
|
|
|
|
|
различных температурах |
|
|
Э |
||||||||||||||
|
t, oС |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
η·107, |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
И |
||||
|
Па∙с |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
5.Контрольные вопросы
1.Какое течение жидкости или газа называется ламинарным?
2.Поясните качественно, от каких параметров зависит объёмный расход жидкости или газа при течении в трубе круглого сечения.
3.Объясните возникновение внутреннего трения при течении газа в трубе на основе молекулярно-кинетической теории.
4.Как проявляется внутреннее трение (вязкость) в жидкости и газе?
5. Запишите уравнение, выражающее закон вязкого трения Ньютона. Дайте определение величинам, входящим в это уравнение.
6. Опишите содержание метода Пуазейля для определения коэффициента динамической вязкости. Нарисуйте принципиальную схему установки.
7. Как и почему изменяется перепад давления на каппиляре при изменении температуры газа?
8. Каков порядок проведения эксперимента. Какие физические величины измеряются непосредственно в опыте?
О9. Чем объяснить увеличение вязкости газов с ростом температуры?
10. Назовите основные источники погрешности при эксперимен-
тальном определении коэффициента динамической вязкости газа методом
Пуазейля. Выведете формулу для расчета абсолютной погрешности кос- |
||
венногоФ |
я η. |
|
Я |
||
измерен |
6. Рекомендуемая литература |
|
1. Савельев И.В. Курс общей физики: Учебное пособие. В 3-х тт. Т.1. |
||
|
С |
|
Механика. Молекулярная физика. 5-е изд., стер. – СПб.: Изд-во «Лань». |
||
2006 г. С. 252 – 258. |
|
|
2. Д.А. Иванов, |
И.В. Иванова, А.Н. Седов, А.В. Славов. Механика. |
|
|
|
Н |
Молекулярная физика и термодинамика: Конспект лекций/ Под ред. А.В. |
||
Кириченко. – М.: Издательство МЭИ. 2003. С.151 – 154, 160– 166. |
||
|
|
И |
|
|
У |
|
|
М |
|
|
Э |
|
|
И |