Изучение гидравлики псевдоожиженного слоя

МИНИСТЕРСТВО ОБРАЗОВАНИЯ И НАУКИ РОССИЙСКОЙ ФЕДЕРАЦИИ ГОСУДАРСТВЕННОЕ ОБРАЗОВАТЕЛЬНОЕ УЧРЕЖДЕНИЕ

ВЫСШЕГО ПРОФЕССИОНАЛЬНОГО ОБРАЗОВАНИЯ

«САМАРСКИЙ ГОСУДАРСТВЕННЫЙ УНИВЕРСИТЕТ»

Химический факультет

Отчет по лабораторной работе

«Изучение гидравлики псевдоожиженного слоя»

Выполнили студенты

курса, 441 группы

Золотарёв Павел

Проверила Журавлева И.И.

Самара 2012

Введение

Гидравлика - наука о законах движения (гидродинамика <#"807743.files/image001.gif">                               

Рис. 1. Неподвижный слой            Рис. 2. Взвешенный слой

Цель работы:

1. Экспериментально изучить зависимость гидравлического сопротивления слоя (∆pcл) от фиктивной скорости газа (ωо).

2. Определить критические скорости газа: скорость псевдоожижения (ωпс) и скорость свободного витания (ωсв).

. Определить эквивалентный диаметр частиц монодисперсного слоя (dэ).

1.     
Основное оборудование и реактивы

а) лабораторная установка (см. рис.1 стр. 2)

б) барометр анероид БАММ-1

Цена деления шкалы: p = 0,1·103 Па; ∆p = ± 0,05·103 Па

в) технические весы марки SCL-150

Цена деления шкалы: m = 5*10-3 г; ∆ m = ± 2,5·10-3г = ±2,5·10-6 кг

г) термометр (0 − 100ºС): Цена деления шкалы: t = 1ºС;∆t = ± 0,5ºС

д) дифманометры (водяные):

Цена деления шкалы: p = 1 мм. вод. ст; ∆p = ± 1 мм. вод. ст.

е) расходомеры (H2SO4 конц.):

Цена деления шкалы: V = 0,1 л/мин = 1,667∙10-6 м3/c; ∆V = ± 1,667∙10-6м3/c

ж) цилиндры (2 шт.): V = 25 мл;

Цена деления шкалы: V = 0,5 см3;∆V = ± 0,25см3 = ±0,25·10-6м3

з) воздуходувка

и) колонка с активированным углем

к) ацетон

В стеклянной колонке (6), внутренним диаметром 37.0±1.0 мм, на металлической сетке находится монодисперсный слой твёрдых частиц активированного угля цилиндрической формы. Воздух, подаваемый воздуходувкой (1), проходит через два расходомера (5, 5А), соединенных параллельно и подаётся в нижнюю часть колонки. Вентиль (3) служит для выпуска воздуха в атмосферу и для регулировки расхода воздуха, подаваемого в колонку. Падение давления воздуха в колонке измеряют дифманометром (4А). Изменение давления в колонке, фиксируемое по дифманометру (4А), представляет собой суммарное сопротивление сетки и собственно слоя. Поскольку сопротивление сетки мало по сравнению с сопротивлением слоя, то падение давления, которое показывает дифманометр (4А), с достаточной степенью точности может быть приравнено к гидравлическому сопротивлению слоя. Температуру, относительное давление воздуха и гидравлическое сопротивление слоя частиц угля, определяют по показаниям термометра (2), дифманометра (4) и дифманометра (4А), соответственно.

Рис.1. Лабораторная установка гидравлики псевдоожиженного слоя

2.  Ход работы

.        Изучили экспериментальную установку.

2.      Полностью открыли вентиль 3, связывающий воздуходувку с атмосферой, после этого закрыли вентили 7 и включили воздуходувку 1. Установили начальный расход воздуха в системе 0 л/мин.

3.      При данном расходе измерили гидравлическое сопротивление и высоту слоя (h0), давление и температуру воздуха.Увеличивая расход воздуха на 2 л/мин (до наиболее возможного), сняли все выше перечисленные показатели.

.        Определили насыпную плотность слоя (ρнас). Для определения плотности частиц (ρчаст) взвесили определенное количество частиц угля; поместив их в мерный цилиндр, измерили объем насыпного слоя. Строго определенное количество ацетона поместили в мерный цилиндр с частицами угля и после полного заполнения пустот ацетоном, измерили суммарный объем (V∑) и объем вытесненной жидкости.

.        Результаты, полученные в ходе работы, справочные данные и рассчитанные величины были сведены в таблицы 1 и 2. На основании полученных данных были построены графические зависимости вида h = f(lg(ωо)) иlg(∆pсл)= f(lg(ωо)).

.        Рассчитали эквивалентный диаметр частиц слоя (dэ) и скорость свободного витания (ωсв) через критические значения критериев Лященко (Lyкр) и Архимеда (Arкр).

.        Провели статистическую обработку полученных результатов.

.        Выводы.

.        Список литературы.

. Определение основных характеристик воздуха

. Абсолютное давление воздуха находим по формуле II,9 (стр. 25):

pабс. = pатм + pотн, (1)

где pатм = 1,036 · 105 Па;

pотн = pизб;

мм. вод. ст. = 9,81 Па (поз. 3 стр. 13).

. Плотность воздуха при давлении р и температуре Трассчитываем по формуле 1.5 (поз. 2, стр. 13):

, (2)

где - плотность воздуха при н.у.,

= 1,293 кг/м3(табл.V, стр. 513);

То = 273 К; Т = (273 + t), К;

pо = 1,013 · 105 Па;p = pабс.

.Вязкость воздуха при давлении р и температуре Т определяем по формуле 1.13 (поз.7,стр.15):

, (3)

где с - постоянная Сатерленда,с = 124;

μ - динамический коэффициент вязкости; μо = 17,3 ·10-6 Па · с [табл. V, стр. 513].

.Фиктивную скорость воздухарассчитываем по уравнению объемного расхода (формула 1.17, поз. 10, стр.16):

V = ωо · f, откудаωо = , (4)

где f - площадь поперечного сечения потока, м2; f = 0,785 · d2вн;

dвн - диаметр колонки, м; dвн = 37 · 10-3 м.

Таблица 1

Характеристики воздуха

4. Определение характеристик слоя твердых частиц активированного угля

а) Насыпную плотность частиц рассчитали по формуле:

rнас=, (5)

где mчаст - масса частиц активированного угля; mчаст = 5 г,

Vсл - объем насыпного слоя, занимаемый частицами; Vсл = 9,5 мл,

тогда rнас= =526,3 кг/м3

б) Плотность частицопределяли по формуле:

rчаст =  , (6)

где Vчаст = V∑ − Vац;

VΣ- суммарный объем, занимаемый частицами и ацетоном;

VΣ = 18 мл; Vац = 15 мл,

тогда rчаст = = 1666,7 кг/м3

в)Порозностьнеподвижного слоя частицεоопределяли по формуле3.44а (поз.29,cтр.104):

εо = (7)

εо =

г) Порозность взвешенного слоя частиц ε выражали из формулы 3.52 (поз. 32, стр. 106):

ε =  , (8)

где εо - порозность неподвижного слоя,

hо − высота неподвижного слоя, м; hо = 0,042 м,

h − высота взвешенного слоя.

Таблица 2

Характеристики слоя твёрдых частиц

На основании полученных данных (табл. 1,2) построены графические зависимости вида h = f(lg(ωо)) и ∆pсл = f(lg(ωо)).Анализируя кривые, которые представлены на рис. 2 (а,б), определили скорость воздуха, необходимую для начала образования взвешенного слоя частиц активированного угля (ωпс).

ωпс = 0,2482 м/с,

тогда плотность и динамический коэффициент вязкости воздуха при ωпс и соответствующей ей температуре (t = 32°С) будут составлять:

ρвозд = 1,210 кг/м3;

µвозд. = 18,905·10-6 Па·с

а

б

Рис. 2. Зависимости высоты зернистого слояа и его гидравлического сопротивлениябот скорости потока воздуха

.       
Расчет эквивалентного диаметра частиц слоя (dэ) и скорости свободного витания (ωсв)

Для нахожденияdэчастиц слоя активированного угля провели следующую цепочку расчетов:

Lyпс→ Ar→ dэ

)Рассчитываемкритическое значение критерия Лященко (Lyпс) по упрощенной формуле 3.5а (поз. 3, стр. 94), так как средой является воздух:

,

где ωпс - скорость псевдоожижения;

ρвозд - плотность воздуха при псевдоожижении;

μвозд - динамический коэффициент вязкости при скорости псевдоожижения

) Определяем критическое значение критерия Архимеда(Ar) по графику зависимости вида Ly=f(Ar) при порозности неподвижного слоя εо = 0,684 (рис. 3.8, поз.35 стр. 108) и Lyпс=0,07

Ar =500

3) Рассчитываем эквивалентный диаметр частиц слоя по упрощенной формуле для Ar (ф. 3.3 поз. 3 стр. 93), так как средой является воздух:

, откуда

м

4) Рассчитываемскорость свободного витания (ωсв), при которой разрушается взвешенный слой частиц и начинается массовый унос частиц потоком воздуха, по упрощенной формуле для Lyсв (ф. 3.4а, поз. 3, стр. 93), так как средой является воздух:

,

Плотность воздуха при свободном витании(ρвозд) оцениваем при tсв, которую определяем графическим методом по зависимости видаε=f(t)путём экстраполяции на ε=1.

Анализируя полученную зависимость, находим, чтосвободное витание должно происходить приtсв=800С.

Применив формулу 1.13(поз.7, стр.15), находим вязкость воздуха при tсв.

= 21,170·10-6Па·с

μ=21,170·10-6Па·с

Рис.3. Зависимость ε взвешенного слоя от температуры потока воздуха

Плотность воздуха при tсвнаходим по зависимости видаρ=f(t), путём экстраполяции на t=800С. (рис.4)ρвозд= 1,015 кг/м3при температуре свободного витания.

Lyсв− критическое значение критерия Лященконаходим порис. 3.8 (стр. 108)

при e =1 и Ar=500;

Lyсв = 6

Рис.4. Зависимость плотности от температуры потока воздуха.

м/с

. Статистическая обработка результатов

Расчет систематических погрешностей

.Насыпная плотность слоя частиц

rнас = lnmчаст - lnVсл

drнас

drнас

Drнас = rнас×drнас = 526,3 × 0,0268 = 14,113 кг/м3

rнас = 526,3 ±14,113 кг/м3.

. Плотность частиц

lnrчаст = lnmчаст - ln(Vå - Vац)

drчаст

drчаст =

Drчаст = rчаст×drчаст = 1666,7 ×0,1505 = 250,83 кг/м3

rчаст = 1666,7 ±250,83 кг/м3.

Все последующие расчеты проводили при скорости псевдоожижения (ωпс).

.Абсолютное давление воздуха

pабс = ратм + ротн,

ln (рабс) = ln(ратм + ротн)

dрабс =

dрабс =

Dрабс = рабс×dрабс =1,05974 · 5,643·10-4= 59,801 Па

рабс = 1,05974·105±59,801 Па.

4. Плотность воздуха

drвозд =

drвозд =

Drвозд = rвозд×drвозд = 1,21 × 0,002111 = 2,55·10-3 кг/м3

rвозд = 1,21 ±2,55·10-3кг/м3.

. Динамический коэффициент вязкости воздуха

,

dm =

dm =

Dm = m×dm = 18,905 × 10-6× 0,0036 = 0,068 × 10-6 Па×с

m = (18,905 ± 0,068) × 10-6 Па.

. Фиктивная скорость воздуха

ωо =

lnωо = lnV − ln(0,785 ×dвн2)

dωо =

dωо =

Dωо = dωо×ωо = 0,2482×= 1,553·10-3

ωо = 0,2482±1,553·10-3м/с

. Порозность неподвижного слоя

lneо = ln(rчаст − rнас) - ln(rчаст)

Deо = deо×eо =0,358× 0,684=0,244

eо = 0,684±0,244

. Порозность взвешенного слоя твердых частиц

e ==

lnε = ln(h - hо + hоεо) - ln(h)

De = de×e = 0,0188× 0,740 = 0,0139

e = 0,740± 0,0139

9. Критерий Lyпс

DLyпс = Lyпс×dLyпс = 0,07× 0,177 = 0,0124

Lyпс = 0,07±0,0124

. Эквивалентный диаметр частиц угля

,

Ddэ = dэ×ddэ = 2,08×10-4× 0,0865 = 1,8·10-5м,

dэ = (20,8 ± 1,8)×10-5м

. Скорость свободного витания

Dωсв = ωсв×dωсв=1,263×0,0855= 0,108м/с

ωсв =1,263±0,108м/с

газ скорость гидравлика псевдоожижение

9. Выводы

1. Изучили зависимость гидравлического сопротивления (∆pсл) слоя от фиктивной скорости воздуха (ωо).

2. Определили критическую скорость воздуха, необходимую для образования взвешенного слоя частиц - скорость псевдоожижения (ωпс), рассчитали вторую критическую скорость воздуха (ωсв), при которой Архимедова сила уравновешивается силой тяжести частиц слоя, через Lyсв и Arкр.

3. Определили эквивалентный диаметр частиц зернистого слоя (dэ).

10. Список литературы

1.       Касаткин А.Г. Основные процессы и аппараты химической технологии. М., «Альянс», 2004 г. - 752 с.

2.      Павлов К.Ф., Романков П.Г., Носков А.А. Примеры и задачи по курсу процессов и аппаратов химической технологии. М.: ООО «РусМедиаКонсалт», 2004 г. - 576 с.