Расчёт и проектирование установки для получения жидкого кислорода

Санкт-Петербургский государственный Университет низкотемпературных и пищевых технологий.

Кафедра криогенной техники.

Курсовой проект.

по дисциплине «Установки ожижения и разделения газовых смесей».

Расчёт и проектирование установки для получения жидкого кислорода.

Работу выполнил.

студент 452 группы.

Денисов Сергей.

Работу принял.

Пахомов О. В.

Санкт — Петербург 2003 год. Оглавление. Задание на расчёт…3 1. Выбор типа установки и его обоснование…3 2. Краткое описание установки…3 3. Общие энергетические и материальные балансы…4 4. Расчёт узловых точек установки…4 5. Расчёт основного теплообменника…7 6. Расчёт блока очистки…17 7. Определение общих энергетических затрат установки…20 8. Расчёт процесса ректификации…20 9. Расчёт конденсатора — испарителя…20 10. Подбор оборудования…21 11.

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

…22.

Задание на расчёт. Рассчитать и спроектировать установку для получения газообразного кислорода с чистотой 99,5%, производительностью 320 м3/ч, расположенную в городе Владивостоке.

1. Выбор типа установки и его обоснование. В качестве прототипа выбираем установку К — 0,4, т. к. установка предназначена для получения жидкого и газообразного кислорода чистотой 99,5%, а также жидкого азота. Также установка имеет относительно несложную схему.

2. Краткое описание работы установки. Воздух из окружающей среды, имеющий параметры Т = 300 К и Р = 0,1 МПа, поступает в компрессорную станцию в точке 1. В компрессоре он сжимается до давления 4,5 МПа и охлаждается в водяной ванне до температуры 310 К. Повышение температуры обусловлено потерями от несовершенства системы охлаждения. После сжатия в компрессоре воздух направляется в теплообменник — ожижитель, где охлаждается до температуры 275 К, в результате чего большая часть содержащейся в ней влаги конденсируется и поступает в отделитель жидкости, откуда выводится в окружающую среду. После теплообменника — ожижителя сжатый воздух поступает в блок комплексной очистки и осушки, где происходит его окончательная очистка от содержащихся в нём влаги и СО2. В результате прохождения через блок очистки воздух нагревается до температуры 280 К. После этого поток сжатого воздуха направляется в основной теплообменник, где охлаждается до температуры начала дросселирования, затем дросселируется до давления Р = 0,65 МПа. В основном теплообменнике поток разделяется. Часть его выводится из аппарата и поступает в детандер, где расширяется до давления Р = 0,65 МПа и поступает в нижнюю часть нижней колонны. Поток из дросселя поступает в середину нижней колонны. Начинается процесс ректификации. Кубовая жидкость (поток R, содержание N2 равно 68%) из низа нижней колонны поступает в переохладитель, где переохлаждается на 5 К, затем дросселируется до давления 0,13 МПа и поступает в середину верхней колонны. Азотная флегма (поток D, концентрация N2 равна 97%) забирается из верхней части нижней колонны, пропускается через переохладитель, где также охлаждается на 5К, затем дросселируется до давления 0,13 МПа и поступает в верхнюю часть верхней колонны. В верхней колонне происходит окончательная ректификация, внизу верхней колонны собирается жидкий кислород, откуда он направляется в переохладитель, где переохлаждается на 8 — 10 К. Далее поток кислорода направляется в жидкостной насос, где его давление поднимается до 10 МПа, и обратным потоком направляется в основной теплообменник. Затем он направляется в теплообменник — ожижитель, откуда выходит к потребителю с температурой 295 К. Азот из верхней части колонны последовательно проходит обратным потоком переохладитель азотной флегмы и кубовой жидкости, оснновной теплообменник и теплообменник — ожижитель. На выходе из теплообменника — ожижителя азот будет иметь температуру 295 К.

3. Общие энергетические и материальные балансы. V = K + A 0,79V = 0,005K + 0,97A МV? i1B — 2B + Vдетhад? адМ = МVq3 + Мк K? i2K — 3K + V? i3В — 4 В М М — молярная масса воздуха. Мк — молярная масса кислорода.

Принимаем V = 1 моль К + А = 1 К = 1 — А 0,79 = 0,005(1 — А) + 0,97А, А = 0,813 К = 1 — 0,813 = 0,187 Определяем теоретическую производительнсть компрессора. (1/0,187) = х/320 => х = 320/0,187 = 1711 м3/ч = 2207,5 кг/ч.

4. Расчёт узловых точек установки Принимаем: Давление воздуха на входе в компрессор… [pic] Давление воздуха на выходе из компрессора… Рвыхк = 4,5 МПА Температура воздуха на входе в компрессор… pic] Температура воздуха на выходе из компрессора… pic] Температура воздуха на выходе из теплообменника — ожижителя… pic] Температура воздуха на выходе из блока очистки…[pic] Давление в верхней колонне… [pic] Давление в нижней колонне…[pic] Концентрация азота в кубовой жидкости … pic] Концентрация азота в азотной флегме… [pic] Температурный перепад азотной флегмы и кубовой жидкости при прохождении через переохладитель… pic] Температура кубовой жидкости… pic] Температура азотной флегмы…[pic] Температура отходящего азота… pic] Температура жидкого кислорода… pic] Разность температур на тёплом конце теплообменника — ожижителя… pic] Температура азота на выходе из установки… pic] Температурный перепад кислорода??? Т1К — 2К = 10 К На начальной стадии расчёта принимаем: [pic].

Составляем балансы теплообменных аппаратов: а) Баланс теплообменника — ожижителя.

КСр к? Т4К — 5К + АСрА? Т3А — 4А = VCpv? T2В — 3 В б) Балансы переохладителя: [pic] [pic] [pic]находим из номограммы [pic]для смеси азот — кислород.

в) Баланс переохладителя кислорода. КCpK? T1К — 2К = RCpR? T2R — 3R.

Принимаем ?T1К — 2К = 10 К.

?T2R — 3R = 0,128*1,686*10/6,621*1,448 = 2,4 Т3R = Т2R + ?T2R — 3R = 74 + 2,4 = 76,4 К i3R = 998,2.

г) Баланс основного теплообменнка. Для определения параметров в точках 3А и 4К разобьём основной теплообменник на 2 трёхпоточных теплообменника: [pic] Истинное значение Vдет вычислим из баланса установки: Vдет = [VMq3 + KMk? i2K — 3K + VM? i4B — 3B — VM? i1B — 2B]/Mhад?ад = [1*29*8 + 0,187*32*(352,8 — 349,9) + 1*29*(522,32 — 516,8) — 1*29*(563,82 — 553,75)]/29*(394,5 — 367,5)*0,7 = 0,2 Vдет = 0,2V = 0,2*1711 = 342 м3/ч Составляем балансы этих теплообменников: I VCpV? T4B — 6B = KCpK? T3K' - 4K + ACpA? T2A' - 3A II (V — Vд) CpV?T6B-5B = KCpK? T3K — 3K' + ACpA? T2A' - 2A Добавим к ним баланс теплообменника — ожижителя. Получим систему из 3 уравнений. III КСр к? Т4К — 5К + АСрА? Т3А — 4А = VCpv? T2В — 3 В Вычтем уравнение II из уравнения I: VCpV? T4B — 6B — (V — Vд) CpV?T6B-5B = KCpK? T3K' - 4K — KCpK? T3K — 3K' + ACpA? T2A' - 3A — ACpA? T2A' - 2A Получаем систему из двух уравнений: I VCpV (T4B — 2T6B + T5B) + VдCpV (T6B — T5B) = KCpK (T4K — T3K) + ACpA? T3A — 2A II КСр к? Т4К — 5К + АСрА? Т3А — 4А = VCpv? T2В — 3 В I 1*1,012(280 — 2*173 + 138) + 0,387*1,093(173 — 138) = 0,128*1,831(T4K — 88) +0,872*1,048(T3А-85) II 1*1,012*(310 — 275) = 0,128*1,093(295 — T4K) + 0,872*1,041(295 — T3А) T4K = 248,4 К T3А = 197,7 К Для удобства расчёта полученные данные по давлениям, температурам и энтальпиям в узловых точках сведём в таблицу:

|№ |1 В |2 В |3 В |4 В |5 В |5 |6 В |7 В |1R |2R |3R | |i, |553,7|563,8|516,8|522,3|319,2|319,2|419,1|367,5|1350 |1131,|1243| |кДж/ | | | | | | | | | |2 | | |кг | | | | | | | | | | | | |Р, |0,1 |4,5 |4,5 |4,5 |4,5 |0,65 |4,5 |4,5 |0,65 |0,65 |0,65| |МПа | | | | | | | | | | | | |Т, К |300 |310 |275 |280 |138 |80 |188 |125 |79 |74 |76,4| |№ |1D |2D |1К |2К |3К |4К |5К |1А |2А |3А |4А | |i, |1015 |2465 |354,3|349,9|352,8|467,9|519,5|328,3|333,5|454,6|553,| |кДж/ | | | | | | | | | | | | |кг | | | | | | | | | | | | |Р, |0,65 |0,65 |0,13 |0,12 |10 |10 |10 |0,13 |0,13 |0,13 |0,13| |МПа | | | | | | | | | | | | |Т, К |79 |74 |93 |84 |88 |248,4|295 |80 |85 |197,7|295 |.

ПРИМЕЧАНИЕ. 1. Значения энтальпий для точек 1R, 2R, 3R, 1D, 2D взяты из номограммы Т — i — P — x — y для смеси азот — кислород. 2. Прочие значения энтальпий взяты из [2].

5. Расчёт основного теплообменника. Ввиду сложности конструкции теплообменного аппарата разобьём его на 4 двухпоточных теплообменника. [pic].

Истинное значение Vдет вычислим из баланса установки: Vдет = [VMq3 + KMk? i2K — 3K + VM? i4B — 3B — VM? i1B — 2B]/Mhад?ад = [1*29*8 + 0,128*32*(352,8 — 349,9) + 1*29*(522,32 — 516,8) — 1*29*(563,82 — 553,75)]/29*(394,5 — 367,5)*0,7 = 0,2 Vдет = 0,2V = 0,2* = 342,2 м3/ч Составляем балансы каждого из четырёх теплообменников: I VA (i4B — i1) + Vq3 = A (i3A — i3) II VK (i4B — i2) + Vq3 = K (i4K — i4) III (VA — Vда)(i1 — i5B) + Vq3 = A (i3 — i2A) IV (VК — Vдк)(i2 — i5B) + Vq3 = К (i4 — i2К) Здесь VA + VК = V, Vда + Vдк = Vд Параметры в точках i1 и i2 будут теми же, что в точке 6 В Температуру в точке 5 В задаём: Т5 В = 138 К Р5 В = 4,5 МПа i5 В = 319,22 кДж/кг = 9257,38 кДж/кмоль Принимаем VA = А = 0,813, VК = К = 0,187, Vдк = Vда = 0,1, q3 = 1 кДж/кг для всех аппаратов. Тогда из уравнения I VA (i4B — i6В) + Vq3 = A (i3A — i3) 0,813(522,32 — 419,1) + 1 = 0,813(454,6 — i3) i3 = (394,6 — 112,5)/0,813 = 324,7 кДж/кг Т3 = 140 К Проверяем полученное значение i3 с помощью уравнения III: (0,872 — 0,1)(394,5 — 319,22) + 1 = 0,872(i3 — 333,5) 59,1 = 0,872i3 — 290,8 i3 = (290,8 + 59,1)/0,872 = 401,3 кДж/кг Уменьшим VА до 0,54: 0,54(522,32 — 419,1) + 1 = 0,872(454,6 — i3) i3 = (394,6 — 70,023)/0,872 = 372,2 кДж/кг Проверяем полученное значение i3 с помощью уравнения III: (0,54 — 0,1)(394,5 — 319,22) + 1 = 0,872(i3 — 333,5) i3 = (290,8 + 34,123)/0,872 = 372,6 кДж/кг Т3 = 123 К Тогда из уравнения II: VK (i4B — i6В) + Vq3 = K (i4K — i4) 0,56(522,32 — 419,1) + 1 = 0,128(467,9 — i4) 72,6 = 59,9 — 0,128 i4 i4 = (72,6 — 59,9)/0,128 = 332 кДж/кг Т4 = 140 К Рассчитываем среднеинтегральную разность температур для каждого из четырёх теплообменников. а) Материальный баланс теплообменника I: VA (i4B — i1) + Vq3 = A (i3A — i3) Из баланса расчитываем истинное значение теплопритоков из окружающей среды: 0,54*1,15(280 — 173) + 1*q3 = 0,872*1,99(197,7 — 123) q3 = 121,9 — 66,4 = 55,5 кДж/кг Рассчитываем коэффициенты В и D: VA (i4B — i6В) + Vq3 = A (i3A — i3) VA? iB + Vq3 = A? iA ?iB = A? iA/ VA — V q3/VA | ?iA/ ?iA ?iB = A? iA/ VA — Vq3* ?iA/ ?iA В = A/VA = 0,872/0,54 = 1,645 D = V q3/VA ?iA = 1*55,5/0,54*(197,7 — 123) = 0,376 ?iB = В? iA — D? iA = С? iA = (1,635 — 0,376) ?iA = 1,259 ?iA.

Составляем таблицу:

|№ |ТВ, К |iв, |?iВ |ТА, К |iА, |?iА | | | |кДж/кг | | |кДж/кг | | |0 — 0 |280 |522,32 |0 |197,7 |454,6 |0 | |1 — 1 |272 |512,0 |10,324 |190,23 |- |8,2 | |2 — 2 |261 |501,7 |20,648 |182,76 |- |16,4 | |3 — 3 |254 |491,3 |30,971 |175,29 |- |24,6 | |4 — 4 |245 |481,0 |41,295 |167,82 |- |32,8 | |5 — 5 |235 |470,7 |51,619 |160,35 |- |41 | |6 — 6 |225 |460,4 |61,943 |152,88 |- |49,2 | |7 — 7 |218 |450,1 |72,267 |145,41 |- |57,4 | |8 — 8 |210 |439,73 |82,59 |137,94 |- |65,6 | |9 — 9 |199 |429,4 |92,914 |130,47 |- |73,8 | |10 — 10|188 |419,12 |103,2 |123 |372,6 |82 |.

Строим температурные кривые:

?Тсринт = n/?(1/?Тср) |№ |?Тср |1/?Тср | |1 |82 |0,012 | |2 |82 |0,012 | |3 |78 |0,0128 | |4 |79 |0,0127 | |5 |77 |0,013 | |6 |72 |0,0139 | |7 |73 |0,0137 | |8 |72 |0,0139 | |9 |69 |0,0145 | |10 |65 |0,0154 |.

?(1/?Тср) = 0,1339 ?Тср = 10/0,1339 = 54,7 К б) Материальный баланс теплообменника II: VK (i4B — i6В) + Vq3 = K (i4K — i4) Из баланса расчитываем истинное значение теплопритоков из окружающей среды: 0,56*1,15(280 — 173) + 1*q3 = 0,187*1,684(248,4 — 140) q3 = 23,4 — 68,9 = -45,5 кДж/кг Рассчитываем коэффициенты В и D: VК (i4B — i6В) + Vq3 = K (i4K — i4) VК? iB + Vq3 = К? iК ?iB = К? iК/ VК — V q3/VК | ?iК/ ?iК ?iB = К? iК/ VК — Vq3* ?iК/ ?iК В = К/VК = 0,128/0,56 = 0,029 D = V q3/VК ?iК = -1*45,5/0,56*(248,4 — 140) = -0,75 ?iB = В? iК — D? iК = С? iК = (0,029 + 0,75) ?iК = 0,779 ?iК Составляем таблицу: |№ |ТВ, К |iв, кДж/кг|?iВ |ТК, К |iК, |?iК | | | | | | |кДж/кг | | |0 — 0 |280 |522,32 |0 |248,4 |332 |0 | |1 — 1 |272 |511,7 |10,589 |237,56 |- |13,593| |2 — 2 |261 |501,1 |21,178 |226,72 |- |27,186| |3 — 3 |254 |490,6 |31,767 |215,88 |- |40,779| |4 — 4 |245 |480 |42,356 |205,04 |- |54,372| |5 — 5 |235 |469,3 |52,973 |194,2 |- |67,975| |6 — 6 |225 |458.8 |63,534 |183,36 |- |81,558| |7 — 7 |218 |448,2 |74,123 |172,52 |- |95,151| |8 — 8 |210 |437,6 |84,735 |161,68 |- |108,77| |9 — 9 |199 |427 |95,301 |150,84 |- |122,33| |10 — 10 |188 |419,12 |105,9 |140 |467,93 |135,93|.

?Тсринт = n/?(1/?Тср) |№ |?Тср |1/?Тср | |1 |32 |0,3 125 | |2 |34 |0,0294 | |3 |34 |0,0294 | |4 |40 |0,025 | |5 |41 |0,0244 | |6 |42 |0,0238 | |7 |45 |0,0222 | |8 |48 |0,0208 | |9 |48 |0,0208 | |10 |48 |0,0208 |.

?(1/?Тср) = 0,245 ?Тср = 10/0,245 = 40,3 К в) Материальный баланс теплообменника III: (VA — Vда)(i6В — i5B) + Vq3 = A (i3 — i2A) Из баланса расчитываем истинное значение теплопритоков из окружающей среды: (0,54 — 0,1)*2,204(188 — 138) + 1*q3 = 0,813*1,684(123 — 85) q3 = 55,8 — 33,9 = 21,9 кДж/кг Рассчитываем коэффициенты В и D: (VA — Vда)(i6В — i5B) + Vq3 = A (i3 — i2A) (VА — Vда) ?iB + Vq3 = А? iА ?iB = А? iА/ (VА — Vда) — V q3/VА | ?iА/ ?iА ?iB = А? iА/ (VА — Vда) — Vq3* ?iА/ ?iА В =А/(VА — Vда) = 0,813/0,44 = 1,982 D = V q3/(VА — Vда) ?iА = 1*21,9/0,44*(372,6 — 333,5) = 0,057 ?iB = В? iА — D? iА = С? iА = (1,982 — 0,057) ?iА = 1,925 ?iА Составляем таблицу: |№ |ТВ, К |iв, кДж/кг|?iВ |ТА, К |iА, |?iА | | | | | | |кДж/кг | | |0 — 0 |188 |394,5 |0 |123 |372,6 |0 | |1 — 1 |175 |387 |7,527 |119,2 |- |3,91 | |2 — 2 |168 |379,4 |15,1 |115,4 |- |7,82 | |3 — 3 |162 |371,92 |22,58 |111,6 |- |11,73 | |4 — 4 |158 |364,4 |30,1 |107,8 |- |15,64 | |5 — 5 |155 |356,9 |37,6 |104 |- |19,55 | |6 — 6 |152 |349,3 |45,2 |100,2 |- |23,46 | |7 — 7 |149 |341,8 |52,7 |96,4 |- |27,37 | |8 — 8 |145 |334,3 |60,2 |92,6 |- |31,28 | |9 — 9 |141 |326,8 |67,741 |88,8 |- |35,19 | |10 — 10 |138 |319,22 |75,28 |85 |333,5 |39,1 |.

?Тсринт = n/?(1/?Тср) |№ |?Тср |1/?Тср | |1 |56 |0,0179 | |2 |53 |0,0189 | |3 |50 |0,02 | |4 |50 |0,02 | |5 |51 |0,0196 | |6 |52 |0,0192 | |7 |53 |0,0189 | |8 |52 |0,0192 | |9 |52 |0,0192 | |10 |53 |0,0189 |.

?(1/?Тср) = 0,192 ?Тср = 10/0,245 = 52 К г) Материальный баланс теплообменника IV: (VК — Vдк)(i6В — i5B) + Vq3 = К (i4 — i2К) Из баланса расчитываем истинное значение теплопритоков из окружающей среды: (0,56 — 0,1)*2,204(188 — 138) + 1*q3 = 0,128*1,742(123 — 88) q3 = 7,804 — 50,7 = - 42,9 кДж/кг Рассчитываем коэффициенты В и D: (VК — Vдк)(i6В — i5B) + Vq3 = К (i4 — i2К) (Vк — Vдк) ?iB + Vq3 = К? iк ?iB = К? iк/ (VК — Vдк) — V q3/VК | ?iК/ ?iК ?iB = К? iК/ (VК — Vдк) — Vq3* ?iК/ ?iК В =К/(VК — Vдк) = 0,128/0,46 = 0,278 D = V q3/(VК — Vдк) ?iк = -1*42,9/0,46*(372,6 — 332) = - 1,297 ?iB = В? iК — D? iК = С? iк = (0,278 + 1,297) ?iК = 1,488 ?iК Составляем таблицу: |№ |ТВ, К |iв, кДж/кг|?iВ |ТК, К |iК, |?iК | | | | | | |кДж/кг | | |0 — 0 |188 |394,5 |0 |140 |332 |0 | |1 — 1 |174 |387,17 |7,33 |134,8 |- |5,06 | |2 — 2 |167 |379,8 |14,7 |129,6 |- |10,12 | |3 — 3 |162 |371,6 |22,9 |124,4 |- |15,18 | |4 — 4 |158 |365,2 |29,3 |119,2 |- |20,24 | |5 — 5 |155 |357,9 |36,6 |114 |- |25,3 | |6 — 6 |152 |350,5 |44 |108,8 |- |30,36 | |7 — 7 |149 |343,2 |51,3 |103,6 |- |35,42 | |8 — 8 |146 |335,9 |58,6 |98,4 |- |40,48 | |9 — 9 |143 |328,6 |65,9 |93,2 |- |45,54 | |10 — 10 |138 |319,22 |75,28 |88 |372,6 |50,6 |.

?Тсринт = n/?(1/?Тср).

|№ |?Тср |1/?Тср | |1 |40 |0,025 | |2 |37 |0,027 | |3 |38 |0,026 | |4 |39 |0,0256 | |5 |41 |0,0244 | |6 |43 |0,0233 | |7 |45 |0,0222 | |8 |47 |0,0213 | |9 |50 |0,02 | |10 |50 |0,02 |.

?(1/?Тср) = 0,235 ?Тср = 10/0,245 = 42,6 К д) Расчёт основного теплообменника. Для расчёта теплообменника разбиваем его на 2 трёхпоточных. Для удобства расчёта исходные данные сводим в таблицу. |Поток |Рср, |Тср, К |Ср, |Уд. Объём |?, |?, | | |ат. | |кДж/кгК |v, м3/кг |кг*с/м2 |Вт/мК, | | | | | | |*107 |*103 | |Прямой |45 |226,5 |1,187 |0,005 |18,8 |23,6 | |(воздух) | | | | | | | |Обратный |100 |190 |2,4 |0,106 |108 |15 | |(О2 под | | | | | | | |дав) | | | | | | | |Обратный |1,3 |155 |1,047 |0,286 |9,75 |35,04 | |(N2 низ | | | | | | | |дав) | | | | | | |.

Прямой поток. 1) Скорость потока принимаем? = 1 м/с 2) Секундный расход Vсек = V*v/3600 = 1711*0,005/3600 = 2,43*10−3 м3/с 3) Выбираем тубку ф 12×1,5 мм 4) Число трубок n = Vсек/0,785dвн? = 0,243/0,785*0,0092*1 = 39 шт Эквивалентный диаметр dэкв = 9 — 5 = 4 мм 5) Критерий Рейнольдса Re =? dвн?/g? = 1*0,004*85,4/9,81*18,8*10−7 = 32 413 6) Критерий Прандтля Pr = 0,802 (см. [2]) 7) Критерий Нуссельта: Nu = 0,023 Re0,8 Pr0,33 = 0,015*324 130,8*0,8020,33 = 63,5 8) Коэффициент теплоотдачи: ?В = Nu?/dвн = 63,5*23,6*10−3/0,007 = 214,1 Вт/м2К Обратный поток (кислород под давлением): 1) Скорость потока принимаем? = 1 м/с 2) Секундный расход Vсек = V*v/3600 = 320*0,0011/3600 = 9,8*10−5 м3/с 3) Выбираем тубку ф 5×0,5 мм гладкую. 4) Критерий Рейнольдса Re =? dвн?/g? = 1*0,007*330,1/9,81*106*10−7 = 21 810 5) Критерий Прандтля Pr = 1,521 (см. [2]) 6) Критерий Нуссельта: Nu = 0,023 Re0,8 Pr0,4 = 0,015*218 100,8*1,5210,33 = 80,3 7) Коэффициент теплоотдачи: ?В = Nu?/dвн = 80,3*15*10−3/0,007 = 172 Вт/м2К Обратный поток (азот низкого давления) 1) Скорость потока принимаем? = 15 м/с 2) Секундный расход Vсек = V*v/3600 = 1391*0,286/3600 = 0,11 м3/с 3) Живое сечение для прохода обратного потока: Fж = Vсек/? = 0,11/15 = 0,0074 м² 4) Диаметр сердечника принимаем Dc = 0,1 м 4) Критерий Рейнольдса Re =? dвн?/g? = 15*0,004*2,188/9,81*9,75*10−7 = 34 313 5) Критерий Нуссельта: Nu = 0,0418 Re0,85 = 0,0418*343 130,85=299,4 7) Коэффициент теплоотдачи: ?В = Nu?/dвн = 299,4*35,04*10−3/0,01 = 1049 Вт/м2К Параметры всего аппарата: 1) Тепловая нагрузка азотной секции QA = A? iA/3600 = 1391*(454,6 — 381,33)/3600 = 28,3 кВт 2) Среднеинтегральная разность температур? Тср = 54,7 К 3) Коэффициент теплопередачи КА = 1/[(1/?в)*(Dн/Dвн) + (1/?А)] = 1/[(1/214,1)*(0,012/0,009) + (1/1049)] = 131,1 Вт/м2 К 4) Площадь теплопередающей поверхности FA = QA/KA ?Тср = 28 300/131,1*54,7 = 3,95 м² 5) Средняя длина трубки с 20% запасом lА = 1,2FA /3,14DHn = 1,2*3,95/3,14*0,012*32 = 3,93 м 6) Тепловая нагрузка кислородной секции QК = К? iA/3600 = 0,183*(467,93 — 332)/3600 = 15,1 кВт 7) Коэффициент теплопередачи КК = 1/[(1/?в) + (1/?К) *(Dн/Dвн)] = 1/[(1/214,1) + (1/172) *(0,01/0,007)]=77 Вт/м2 К 8) Площадь теплопередающей поверхности FК = QК/KК ?Тср = 15 100/77*25 = 7,8 м² 9) Средняя длина трубки с 20% запасом lК = 1,2FК /3,14DHn = 1,2*7,8/3,14*0,01*55 = 5,42 м Принимаем l = 5,42 м. 10) Теоретическая высота навивки. Н = lt2/?Dср = 17*0,0122/3,14*0,286 = 0,43 м.

Второй теплообменник.

|Поток |Рср, |Тср, К |Ср, |Уд. Объём |?, кг*с/м2|?, | | |ат. | |кДж/кгК |v, м3/кг | |Вт/мК, | | | | | | |*107 |*103 | |Прямой |45 |155,5 |2,328 |0,007 |142,62 |23,73 | |(воздух) | | | | | | | |Обратный |100 |132,5 |1,831 |0,104 |943,3 |106,8 | |(О2 под | | | | | | | |дав) | | | | | | | |Обратный |1,3 |112,5 |1,061 |0,32 |75,25 |10,9 | |(N2 низ | | | | | | | |дав) | | | | | | |.

Прямой поток. 1) Скорость потока принимаем? = 1 м/с 2) Секундный расход Vсек = V*v/3600 = 1875*0,007/3600 = 2,6*10−3 м3/с.

3) Выбираем тубку ф 10×1,5 мм гладкую. 4) Число трубок n = Vсек/0,785dвн? = 0,0026/0,785*0,0072*1 = 45 шт Эквивалентный диаметр dэкв = 9 — 5 = 4 мм 5) Критерий Рейнольдса Re =? dвн?/g? = 1*0,004*169,4/9,81*142,62*10−7 = 83 140 6) Критерий Прандтля Pr =1,392 (см. [2]) 7) Критерий Нуссельта: Nu = 0,023 Re0,8 Pr0,33 = 0,015*831 400,8*1,3920,33 = 145 8) Коэффициент теплоотдачи: ?В = Nu?/dвн = 145*10,9*10−3/0,007 = 225,8 Вт/м2К Обратный поток (кислород под давлением): 1) Скорость потока принимаем? = 1 м/с 2) Секундный расход Vсек = V*v/3600 = 800*0,104/3600 = 1,2*10−4 м3/с 3) Выбираем тубку ф 10×1,5 мм с оребрением из проволоки ф 1,6 мм и шагом оребрения tп = 5,5 мм 4) Критерий Рейнольдса Re =? dвн?/g? = 1*0,007*1067,2/9,81*75,25*10−7 = 101 200 5) Критерий Прандтля Pr = 1,87 (см. [2]) 6) Критерий Нуссельта: Nu = 0,023 Re0,8 Pr0,4 = 0,015*1 012 000,8*1,870,33 = 297,2 7) Коэффициент теплоотдачи: ?В = Nu?/dвн = 297,2*10,9*10−3/0,007 = 462,8 Вт/м2К Обратный поток (азот низкого давления) 1) Скорость потока принимаем? = 15 м/с 2) Секундный расход Vсек = V*v/3600 = 2725*0,32/3600 = 0,242 м3/с 3) Живое сечение для прохода обратного потока: Fж = Vсек/? = 0,242/15 = 0,016 м² 4) Диаметр сердечника принимаем Dc = 0,1 м 4) Критерий Рейнольдса Re =? dвн?/g? = 15*0,01*3,04/9,81*75,25*10−7 = 60 598 5) Критерий Нуссельта: Nu = 0,0418 Re0,85 = 0,0418*605 980,85=485,6 7) Коэффициент теплоотдачи: ?В = Nu?/dвн = 485,6*10,9*10−3/0,01 = 529,3 Вт/м2К Параметры всего аппарата: 1) Тепловая нагрузка азотной секции QA = A? iA/3600 = 2725(391,85 — 333,5)/3600 = 57 кВт 2) Среднеинтегральная разность температур? Тср = 52 К 3) Коэффициент теплопередачи КА = 1/[(1/?в)*(Dн/Dвн) + (1/?А)] = 1/[(1/225,8)*(0,01/0,007) + (1/529,3)] = 121,7 Вт/м2 К 4) Площадь теплопередающей поверхности FA = QA/KA ?Тср = 57 000/121,7*52 = 9 м² 5) Средняя длина трубки с 20% запасом lА = 1,2FA /3,14DHn = 1,2*9/3,14*0,01*45 = 7,717 м 6) Тепловая нагрузка кислородной секции QК = К? iК/3600 = 0,128*(352,8 — 332)/3600 = 4,6 кВт 7) Коэффициент теплопередачи КК = 1/[(1/?в) + (1/?К) *(Dн/Dвн)] = 1/[(1/225,8) + (1/529,3) *(0,01/0,007)] = 140,3 Вт/м2 К 8) Площадь теплопередающей поверхности FК = QК/KК ?Тср = 4600/140*42,6 = 0,77 м² 9) Средняя длина трубки с 20% запасом lК = 1,2FК /3,14DHn = 1,2*0,77/3,14*0,01*45 = 0,654 м Принимаем l = 7,717 м. 10) Теоретическая высота навивки. Н = lt2/?Dср = 7,717*0,0122/3,14*0,286 = 0,33 м.

Окончательный вариант расчёта принимаем на ЭВМ.

6. Расчёт блока очистки. 1) Исходные данные: Количество очищаемого воздуха … V = 2207,5 кг/ч = 1711 м3/ч Давление потока … Р = 4,5 МПа Температура очищаемого воздуха… Т = 275 К Расчётное содержание углекислого газа по объёму… С = 0,03% Адсорбент… NaX Диаметр зёрен … dз = 4 мм Насыпной вес цеолита??? ц = 700 кг/м3 Динамическая ёмкость цеолита по парам СО2… ад = 0,013 м3/кг Принимаем в качестве адсорберов стандартный баллон диаметром Da = 460 мм и высоту слоя засыпки адсорбента L = 1300 мм. 2) Скорость очищаемого воздуха в адсорбере:? = 4Va/n?Da2 n — количество одновременно работающих адсорберов; Vа — расход очищаемого воздуха при условиях адсорбции, т. е. при Р = 4,5 МПа и Тв = 275 К: Va = VTB P/T*PB = 1711*275*1/273*45 = 69,9 кг/ч? = 4*69,9/3*3,14*0,462 = 140,3 кг/ч*м2 Определяем вес цеолита, находящегося в адсорбере: Gц = nVад? ц = L*?*n*?*Da2/4 = 1*3,14*0,462*1,3*700/4 = 453,4 кг Определяем количество СО2, которое способен поглотить цеолит: VCO2 = Gц*aд = 453,4*0,013 = 5,894 м³ Определяем количество СО2, поступающее каждый час в адсорбер: VCO2' = V*Co = 3125*0,0003 = 0,937 м3/ч Время защитного действия адсорбента: ?пр = VCO2/ VCO2' = 5,894/0,937 = 6,29 ч Увеличим число адсорберов до n = 4. Тогда:? = 4*69,9/4*3,14*0,462 = 105,2 кг/ч*м2 Gц = 4*3,14*0,462*1,3*700/4 = 604,6 кг VCO2 = Gc *aд = 604,6*0,013 = 7,86 м³? пр = 7,86/0,937 = 8,388 ч. Выбираем расчётное время защитного действия? пр = 6 ч. с учётом запаса времени.

2) Ориентировочное количество азота для регенерации блока адсорберов: Vрег = 1,2*GH2O /x' ?рег GH2O — количество влаги, поглощённой адсорбентом к моменту регенерации GH2O = GцаН2О = 604,2*0,2 = 120,84 кг? рег — время регенерации, принимаем? рег = 0,5 ?пр = 3 ч. х' - влагосодержание азота при Тср. вых и Р = 105 Па: Тср. вых = (Твых.1 + Твых.2)/2 = (275 + 623)/2 = 449 К х = 240 г/м3 Vрег = 1,2*120,84/0,24*3 = 201,4 м3/ч Проверяем количество регенерирующего газа по тепловому балансу: Vрег *?N2*CpN2*(Твх + Твых. ср)* ?рег = ?Q ?Q = Q1 + Q2 + Q3 + Q4 + Q5 Q1 — количество тепла, затраченное на нагрев металла; Q2 — количество тепла, затраченное на нагрев адсорбента, Q3 — количество тепла, необходимое для десорбции влаги, поглощённой адсорбентом; Q4 — количество тепла, необходимое для нагрева изоляции; Q5 — потери тепла в окружающую среду. Q1 = GмСм (Тср' - Tнач') Gм — вес двух баллонов с коммуникациями; См — теплоёмкость металла, См = 0,503 кДж/кгК Tнач' - температура металла в начале регенерации, Tнач' = 280 К Тср' - средняя температура металла в конце процесса регенерации, Тср' = (Твх' + Твых')/2 = (673 + 623)/2 = 648 К Твх' - температура азота на входе в блок очистки, Твх' = 673 К; Твых' - температура азота на выходе из блока очистки, Твх' = 623 К;

Для определения веса блока очистки определяем массу одного баллона, который имеет следующие геометрические размеры: наружний диаметр… Dн = 510 мм, внутренний диаметр… Dвн = 460 мм, высота общая… Н = 1500 мм, высота цилиндрической части… Нц = 1245 мм. Тогда вес цилиндрической части баллона GM' = (Dн2 — Dвн2) Нц*?м*?/4 = (0,512 — 0,462)*1,245*7,85*103*3,14/4 = 372,1 кг, где? м — удельный вес металла, ?м = 7,85*103 кг/м3. Вес полусферического днища GM'' = [(Dн3/2) — (Dвн3/2)]* ?м*4?/6 = [(0,513/2) — (0,463/2)]*7,85*103*4*3,14/6 = 7,2 кг Вес баллона: GM' + GM'' = 382 + 7,2 = 389,2 кг Вес крышки с коммуникациями принимаем 20% от веса баллона: GM''' = 389,2*0,2 = 77,84 кг Вес четырёх баллонов с коммуникацией: GM = 4(GM' + GM'' + GM''') = 4*(382 + 7,2 + 77,84) = 1868 кг. Тогда: Q1 = 1868*0,503*(648 — 275) = 3,51*105 кДж Количество тепла, затрачиваемое на нагревание адсорбента: Q2 = GцСц (Тср' - Tнач') = 604,6*0,21*(648 — 275) = 47 358 кДж Количество тепла, затрачиваемое на десорбцию влаги: Q3 = GH2OCp (Ткип — Тнач') + GH2O*? = 120,84*1*(373 — 275) + 120,84*2258,2 = 2,8*105 кДж? — теплота десорбции, равная теплоте парообразования воды; Ср — теплоёмкость воды. Количество тепла, затрачиваемое на нагрез изоляции: Q4 = 0,2Vиз ?изСиз (Тиз — Тнач) = 0,2*8,919*100*1,886*(523 — 275) = 8,3*104 кДж Vиз = Vб — 4Vбалл = 1,92*2,1*2,22 — 4*0,20 785*0,512*0,15 = 8,919 м3 — объём изоляции. ?из — объёмный вес шлаковой ваты, ?из = 100 кг/м3 Сиз — средняя теплоёмкость шлаковой ваты, Сиз = 1,886 кДж/кгК Потери тепла в окружающую среду составляют 20% от? Q = Q1 + Q2 + Q4: Q5 = 0,2*(3,51*105 + 47 358 + 8,3*104) = 9.63*104 кДж Определяем количество регенерирующего газа: Vрег = (Q1 + Q2 + Q3 + Q4 + Q5)/ ?N2*CpN2*(Твх + Твых. ср)* ?рег = =(3,51*105 + 47 358 + 2,8*105 + 8,3*104 + 9,63*104)/(1,251*1,048*(673 — 463)*3) = 1038 нм3/ч Проверяем скорость регенерирующего газа, отнесённую к 293 К: ?рег = 4 Vрег*293/600*?*Da2 *n*Tнач = 4*1038*293/600*3,14*0,462*2*275 = 5,546 м/с n — количество одновременно регенерируемых адсорберов, n = 2.

Определяем гидравлическое сопротивление слоя адсорбента при регенерации. ?Р = 2f? L?2/9,8dэх2 где? Р — потери давления, Па; f — коэффициент сопротивления;? — плотность газа, кг/м3; L — длина слоя сорбента, м; dэ — эквивалентный диаметр каналов между зёрнами, м;? — скорость газа по всему сечению адсорбера в рабочих условиях, м/с;? — пористость слоя адсорбента,? = 0,35 м2/м3. Скорость регенерирующего газа при рабочих условиях:? = 4*Vрег*Твых.ср./3600*?*Da2*n*Тнач = 4*1038*463/3600*3,14*0,462*2*275 = 1,5 м/с Эквивалентный диаметр каналов между зёрнами: dэ = 4*?*dз/6*(1 — ?) = 4*0,35*4/6*(1 — 0,35) = 1,44 мм. Для определения коэффициента сопротивления находим численное значение критерия Рейнольдса: Re = ?*dэ*?/?*?*g = 1,5*0,144*0,79*107/0,35*25*9,81 = 198,8 где? — динамическая вязкость,? = 25*10−7 Па*с;? — удельный вес азота при условиях регенерации,? = ?0 *Р*Т0/Р0*Твых.ср = 1,251*1,1*273/1,033*463 = 0,79 кг/м3 По графику в работе [6] по значению критерия Рейнольдса определяем коэффициент сопротивления f = 2,2 Тогда: ?Р = 2*2,2*0,79*1,3*1,52/9,81*0,144*0,352 = 587,5 Па Определяем мощность электроподогревателя: N = 1,3* Vрег*?*Ср*(Твх — Тнач)/860 = 1,3*1038*1,251*0,25(673 — 293)/860 = 70,3 кВт где Ср = 0,25 ккал/кг*К.

7. Определение общих энергетических затрат установки l = [V?в RToc ln (Pk/Pn)]/?из Кж*3600 = 1711*0,287*296,6*ln (4,5/0,1)/0,6*320*3600 = 0,802 кВт где V — полное количество перерабатываемого воздуха, V = 2207,5 кг/ч = 1711 м3/ч ?в — плотность воздуха при нормальных условиях, ?в = 1,29 кг/м3 R — газовая постоянная для воздуха, R = 0,287 кДж/кгК ?из — изотермический КПД, ?из = 0,6 Кж — количество получаемого кислорода, К = 320 м3/ч Тос — температура окружающей среды, принимается равной средне — годовой температуре в городе Владивостоке, Тос = 23,60С = 296,6 К.

8. Расчёт процесса ректификации. Расчёт процесса ректификации производим на ЭВМ (см. распечатки 4 и 5). Вначале проводим расчёт нижней колонны. Исходные данные вводим в виде массива. Седьмая строка массива несёт в себе информацию о входящем в колонну потоке воздуха: принимаем, что в нижнюю часть нижней колонны мы вводим жидкий воздух. 1 — фазовое состояние потока, жидкость; 0,81 — эффективность цикла. Поскольку в установке для ожижения используется цикл Гейландта (х = 0,19), то эффективность установки равна 1 — х = 0,81. 0,7812 — содержание азота в воздухе; 0,0093 — содержание аргона в воздухе; 0,2095 — содержание кислорода в воздухе. Нагрузку конденсатора подбираем таким образом, чтобы нагрузка испарителя стремилась к нулю.

8. Расчёт конденсатора — испарителя. Расчёт конденсатора — испарителя также проводим на ЭВМ с помощью программы, разработанной Е. И. Борзенко.

В результате расчёта получены следующие данные (смотри распечатку 6): Коэффициент телоотдачи в испарителе… ALFA1 = 1130,7 кДж/кгК Коэффициент телоотдачи в конденсаторе… ALFA2 = 2135,2 кДж/кгК Площадь теплопередающей поверхности… F1 = 63,5 м³ Давление в верхней колонне… Р1 = 0,17 МПа.

10. Подбор оборудования. 1. Выбор компрессора. Выбираем 2 компрессора 605ВП16/70. Производительность одного компрессора …16±5% м3/мин Давление всасывания…0,1 МПа Давление нагнетания…7 МПа Потребляемая мощность…192 кВт Установленная мощность электродвигателя…200 кВт 2. Выбор детандера. Выбираем ДТ — 0,¾. Характеристики детандера: Производительность… V = 340 м3/ч Давление на входе… Рвх = 4 МПа Давление на выходе… Рвых = 0.6 МПа Температура на входе… Твх = 188 К Адиабатный КПД??? ад = 0,7 3. Выбор блока очистки. Выбираем стандартный цеолитовый блок осушки и очистки воздуха ЦБ — 2400/64. Характеристика аппарата: Объёмный расход воздуха… V=2400 м3/ч Рабочее давление: максимальное… Рмакс = 6,4 МПа минимальное… Рмин = 3,5 МПа Размеры сосудов…750×4200 мм. Количество сосудов…2 шт. Масса цеолита… М = 2060 кг Список используемой литературы: 1. Акулов Л. А., Холодковский С. В. Методические указания к курсовому проектированию криогенных установок по курсам «Криогенные установки и системы» и «Установки сжижения и разделения газовых смесей» для студентов специальности 1603. — СПб.; СПбТИХП, 1994. — 32 с. 2. Акулов Л. А., Борзенко Е. И., Новотельнов В. Н., Зайцев А. В. Теплофизические свойства криопродуктов. Учебное пособие для ВУЗов. — СПб.: Политехника,.

2001. — 243 с. 3. Архаров А. М. и др. Криогенные системы: Основы теории и расчёта: Учебное пособие для ВУЗов, том 1., — М.: Машиностроение, 1998. — 464 с. 4. Архаров А. М. и др. Криогенные системы: Основы теории и расчёта: Учебное пособие для ВУЗов, том 2., — М.: Машиностроение, 1999. — 720 с. 5. Акулов Л. А., Холодковский С. В. Криогенные установки (атлас технологических схем криогенных установок): Учебное пособие. — СПб.:

СПбГАХПТ, 1995. — 65 с. 6. Кислород. Справочник в двух частях. Под ред. Д. Л. Глизманенко. М., «Металлургия», 1967.

Распечатка 1. Расчёт основного теплообменника.

Распечатка 2. Расчёт теплообменника — ожижителя.

Распечатка 3. Расчёт переохладителя.

Распечатка 4. Расчёт процесса ректификации в нижней колонне.

Распечатка 5. Расчёт процесса ректификации в верхней колонне.

Распечатка 6. Расчёт конденсатора — испарителя.

Распечатка 7. Расчёт переохладителя кислорода. ———————————- [pic].

[pic].

[pic].

[pic].