Моделирование теплообменных аппаратов в стационарном режиме

Размещено на http://www.allbest.ru/

1. Моделирование теплообменных аппаратов в стационарном режиме

теплообменный аппарат химический температурный

Цель работы

1. Cоставить уравнения тепловых балансов потоков в теплообменных аппаратах.

2. Составить алгоритм расчета в соответствии с вариантом заданий.

3. Составить программу вычислений для ЭВМ.

4. Рассчитать температурные профили потоков в теплообменном аппарате.

5. Провести исследования влияния технологических параметров и анализ полученных результатов.

1.1 Моделирование и интенсификация работы теплообменной аппаратуры

Реализация тепловых процессов в промышленности требует установки крупногабаритного теплообменного оборудования с большой площадью поверхности теплопередачи. Например, в агрегатах синтеза аммиака большой единичной мощности (1360 т/с) АМ-70 и АМ-76 из 205 единиц основного оборудования 57 составляют различные типы теплообменных аппаратов с общей поверхностью теплообмена 150000 м2. На изготовление теплообменных аппаратов ежегодно расходуется большое количество остродефицитных труб из нержавеющей стали и титана.

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

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

Наиболее широкое применение в настоящее время находят рекуперативные теплообменники, которые по своим основным конструктивным признакам делятся:

· на теплообменные аппараты, изготавливаемые из труб различной формы и диаметров (кожухотрубчатые, «труба в трубе», змеевиковые погружные, оросительные, витые, воздушного охлаждения с оребренными трубами);

· теплообменные аппараты, изготовленные из листа (пластинчатые: разборные, полуразборные, неразборные, пластинчато-ребристые; спиральные; ламельные и панельные);

· теплообменные аппараты, совмещенные с различными типами химических аппаратов и реакторов.

Теплообменники могут быть изготовлены из различных металлов, графита и фторопластов различных типов.

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

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

Конструктивные особенности теплообменных аппаратов определяют область, в которой они могут быть применены для различных температур и давлений. Наиболее широко применяются трубчатые теплообменные аппараты, работающие в широком диапазоне температур (от -200 до +475 оС) и давлении до 0,6 МПа. Однако эти аппараты имеют низкие коэффициенты теплопередачи (1000-1500 Вт/м2К), высокую металлоемкость (до 37 кг/м2); для их изготовления необходимы остродефицитные из нержавеющей стали бесшовные трубы и значительные трудозатраты. Они имеют низкую степень унификации узлов и деталей - 10-12 %.

Более низкие значения коэффициентов теплопередачи и высокую металлоемкость имеют оросительные теплообменники, теплообменники «труба в трубе», аппараты воздушного охлаждения, которые, однако, могут работать при более высоких давлениях, чем кожухотрубные.

Научно-технический прогресс в химическом машиностроении в последние годы характеризуется созданием большого количества высокопроизводительного оборудования большой единичной мощности, в том числе теплообменных аппаратов для химической, нефтехимической и микробиологической промышленности.

Анализ параметров работы кожухотрубчатых теплообменников в химической и смежных областях промышленности показывает, что около 70 % теплообменников применяется для давлений до 0,1 МПа и температур до 200оС. Увеличение теплосъёма на единицу площади теплообменного оборудования кожухотрубчатого типа обычно сопровождается экстенсивным ростом поверхности теплообмена, размеров, массы и его стоимости.

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



1.2 Основные закономерности теплообмена

Тепловые процессы в химической технологии имеют как самостоятельное значение при сушке, выпаривании, нагревании, охлаждении и т. д., так и сопровождают химические и массообменные процессы [3-5].

Теплообмен обусловлен стремлением системы к тепловому равновесию, т. е. к равенству температур в системе. В случае неоднородности поля температур возникает поток энергии. Связь между градиентом температуры и молекулярным потоком теплоты (qT) определена законом теплопроводности Фурье:

, (1)

где - коэффициент теплопроводности среды, Вт/(мК);

qT, Вт/м2.

При движении в жидкостях и газах происходит конвективный перенос энергии веществом:

, (2)

где u - скорость движения среды, м/с;

- плотность вещества, кг/м3;

I - энтальпия, Дж/кг.

Таким образом, при конвективном теплообмене плотность теплового потока q определяется суммой молекулярной и конвективной составляющих:

. (3)



Этот процесс достаточно сложный, поэтому для удобства расчета переноса теплоты между поверхностью твердого тела и движущейся сплошной средой используют закон теплоотдачи Ньютона - Рихмана:

, (4)

где - коэффициент теплоотдачи, Вт/(м2К);

F - поверхность теплообмена, м2;

Тс - температура стенки;

Тж - температура среды.

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

Теплоотдача определяется не только тепловыми, но и гидродинамическими условиями. Поэтому конвективный теплообмен описывается дифференциальным уравнением Фурье - Кирхгофа

, (5)

где - коэффициент температуропроводности, м2/с;

t - время, с.

Количество тепла, передаваемое от нагретого теплоносителя к холодному теплоносителю, определяется основным уравнением теплопередачи [5, 13]



, (6)

где - коэффициент теплопередачи, Вт/(м2С);

- средняя разность температур между теплоносителями.

При теплопередаче через стенку толщиной с коэффициент тепло-передачи можно рассчитать с помощью уравнения аддитивности термических сопротивлений на пути теплового потока:

, (7)

где 1 и 2 - коэффициенты теплоотдачи от жидкости к стенке и от стенки к другой жидкости соответственно, Вт/(м2К), (Дж/с м2К);

с - теплопроводность материала стенки, Вт/(мК);

rЗ1 и rЗ2 - термические сопротивления слоёв загрязнений с обеих сторон стенки, м2К/Вт.

Это уравнение справедливо для передачи тепла через плоскую или цилиндрическую стенку при условии, что ( и - наружный и внутренний радиусы цилиндра соответственно).

Ориентировочные значения коэффициентов теплопередачи KТ в Вт/(м2К), а также значения теплопроводности загрязнений стенок в Вт/(м2К) приведены ниже.

Вид теплообмена	Коэффициент теплопередачи
	для вынужденного движения	для свободного движения
От газа к газу. От газа к жидкости. От конденсирующегося пара к газу. От жидкости к жидкости: - для воды; - для углеводорода, масел. От конденсирующегося водяного пара к воде. От конденсирующегося водяного пара к органическим жидкостям. От конденсирующегося пара органических жидкостей к воде. От конденсирующегося водяного пара к кипящей жидкости	10-40 10-60 10-60 800-1700 120-270 800-3500 120-340 300-800 -	4-12 6-20 6-12 140-340 30-60 300-1200 60-170 230-460 300-2500
Теплоноситель. Вода: - загрязнённая; - среднего качества; - хорошего качества; - дистиллированная.	1400-1860 1860-2900 2900-5800 11600
Воздух. Нефтепродукты, масла, пары хладагентов. Нефтепродукты сырые. Органические жидкости, рассолы, жидкие хладагенты. Водяной пар, содержащий масла. Пары органических жидкостей	2800 2900 1160 5800 5800 11600

1.3 Модели процессов теплообмена

При построении математических моделей теплообменных аппаратов предварительно проводится структурный анализ по выявлению количества и видов однородных потоков тепловой энергии, имеющих место в аппарате. Для каждого потока записывается математическое описание в виде выражения, характеризующего изменения температуры в потоке теплоносителя во времени, обусловленное движением потока и теплопередачей [3, 5]. Предварительно формулируются допущения.

Если структура потока теплоносителя соответствует модели идеального перемешивания, то для математического описания этого потока можно использовать уравнение (8) с учетом теплопередачи:



, (8)

где V - объем потока идеального перемешивания, м3;

- плотность теплоносителя, кг/м3;

- удельная теплоемкость теплоносителя, Дж/(кгК);

- объемная скорость потока, м3/с;

F - поверхность теплообмена, м2;

- коэффициент теплопередачи, Вт/(м2К);

(Тст - Т) - разность температур стенки и теплоносителя, К;

Твх - температура потока на входе, К;

t - время, с.

Если структура потока соответствует модели идеального вытеснения, то для математического описания потока можно использовать уравнение (8) с учетом теплопередачи:

, (9)

где Sв - площадь поперечного сечения потока, м2;

L - длина зоны идеального вытеснения, м;

l - пространственная координата, изменяющаяся от 0 до L;

Т=Т (l, t) - функция распределения температуры потока теплоносителя по пространственной координате во времени.

Обычно в уравнениях (8) и (9) принимают коэффициент теплоотдачи, плотность и теплоемкость теплоносителя постоянными в исследуемом ограниченном интервале изменения температуры. Предполагается, что объемные скорости потоков остаются постоянными.

Cоставим математические модели некоторых типов теплообменных аппаратов.

Теплообменник типа «перемешивание-перемешивание» (рис. 1)

Примем, что тепло передается от первого потока теплоносителя ко второму. Режим движения потоков - идеальное перемешивание.

Рис. 1 Схематическое изображение теплообменника типа «перемешивание-перемешивание»

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

(10)

где КТ - коэффициент теплопередачи.

Если тепловой емкостью стенки, разделяющей потоки теплоносителя, пренебречь нельзя, то необходимо к уравнениям (10) добавить уравнение изменения температуры:

(11)



где G3 - вес стенки, кг;

С3 - удельная теплоемкость металла стенки, ; Т3 - температура стенки, К; 1, 2 - коэффициенты теплоотдачи, Вт/м2К.

Теплообменник типа «перемешивание - вытеснение» (рис. 2)

Рис. 2 Схематическое изображение теплообменника типа «перемешивание - вытеснение»

Без учета тепловой емкости стенки:

; (12)

с начальным условием Т2(l,0) = T20(l) и граничным условием Т2(0,t)=T20(t). С учётом теплоёмкости стенки:

;

; (13)



с начальным условием Т2(x,0) = T20(x) и граничным условием Т2(0,) = T20().

Теплообменник типа «вытеснение-вытеснение» (рис. 3)

Рис. 3 Схематическое изображение теплообменника типа «вытеснение-вытеснение»

· Без учета тепловой емкости стенки:

; (14)

с начальными условиями Т1(l,0) = T10(l) и Т2(l,0) = T20(l) и граничными условиями: для прямотока Т1(0,t) = T10(t), Т2(0,t) = T20(t);

для противотока Т1(0,t) = T10(t), Т2(L,t) = T20(t).

· С учётом теплоёмкости стенки:

;

; (15)



с начальными условиями и граничными условиями, заданными для системы (14). Если необходимо учесть влияние теплоемкости внешней стенки, то добавляем еще одно уравнение, соответствующее теплообмену одной среды с внешней стенкой, и система общих уравнений (15) примет вид:

;

;

;

,

где G4, C4 - вес и удельная теплоёмкость металла внешней стенки, изолированной от окружающей среды;

F4 - поверхность теплообмена между потоком теплоносителя с температурой Т1 и внешней стенкой;

4 - коэффициент теплоотдачи.

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

,

получаем Т1 = Т4, т. е. температура наружной стенки равна температуре среды, обменивающейся теплом со стенкой. Наличие толстых стенок сказывается на характере переходных процессов в аппаратах, но не на характере стационарного режима их работы.



1.4 Моделирование теплообменных процессов

В теплообменнике типа «труба в трубе» охлаждается жидкость. Хладоагент и охлаждающаяся жидкость движутся прямотоком (рис. 4).

Определить температуру теплоносителей на выходе из аппарата. Построить температурные профили по длине аппарата.

Рис. 4 Теплообменник типа «труба в трубе»

	Горячий теплоноситель	Холодный теплоноситель
Температура, оС	200	35
Объёмная скорость, м3/с	2,310-4	5,110-4
Плотность, кг/м3	900	1000
Теплоёмкость, Дж/кгоС	3,35103	4,19103
Диаметр трубы, м	0,01	0,03

В теплообменнике реализуется режим «вытеснение-вытеснение». Поэтому математическое описание будет иметь вид:

;

.



В стационарном режиме работы теплообменника уравнения теплового баланса примут следующий вид:

(16)

где d - диаметр трубы, м.

Для удобства вычисления введём обозначения:

;

.

Систему полученных дифференциальных уравнений (16) решаем с помощью численного метода Эйлера:

где i - номер шага по длине теплообменника;

h - шаг интегрирования по длине теплообменного аппарата.

Блок-схема алгоритма расчёта теплообменника приведена на рис. 5. Программа расчета теплообменника представлена в Приложении Д.



Рис. 5 Блок-схема программы расчёта теплообменника



Литература

1. Панченков Г. М., Лебедев В. П. Химическая кинетика и катализ. - М.: Химия, 1985. - 589 с.

2. Яблонский Г. С., Быков В. И., Горбань А. И. Кинетические модели каталитических реакций. - Новосибирск: Наука, 198 - 254 с.

Кафаров В. В. Методы кибернетики в химии и химической технологии. - М.: Химия, 1988. - 489 с.

4. Кравцов А. В., Новиков А. А., Коваль П. И. Методы анализа химико- технологических процессов. - Томск: изд-во ТПУ, 1994. - 76 с.

5. Кафаров В. В., Глебов М. В. Математическое моделирование основных процессов химических производств. - М.: Высш. шк., 1991. - 400 с.

6. Мойзес О. Е., Коваль П. И., Баженов Д. А., Кузьменко Е. А. Информатика: учеб. пособие. В 2-х ч. - Томск, 1999. - 150 с.

7. Турчак Л. И. Основы численных методов. - М.: Наука, 1987. - 320 с.

8. Офицеров Д. В., Старых В. А. Программирование в интегрированной среде Турбо-Паскаль. - Минск: Беларусь, 1992. - 240 с.

9. Бесков В. С., Флор К. В. Моделирование каталитических процессов и реакторов. - М.: Химия, 1991. - 252 с.

10. Руд Р., Праустниц Дж., Шервуд Т. Свойства газов и жидкостей

/ под ред. Б. И. Соколова. - Л.: Химия, 1982. - 591 с.

11. Танатаров М. А. и др. Технологические расчеты установок переработки нефти. - М.: Химия, 1987. - 350 с.

12. Жоров Ю. М. Термодинамика химических процессов. - М.: Химия, 1985. - 458 с.

1 Расчеты основных процессов и аппаратов нефтепереработки: справочник / под ред. Е. Н. Судакова. - М.: Химия, 1979. - 568 с.

14. Кафаров В. В. Разделение многокомпонентных систем в химической технологии. Методы расчета. - М.: Московский химико-технологический институт, 1987. - 84 с

Размещено на Allbest.ru