Расчет оптимального теплообменника по параметрам эффективности теплопередачи

Содержание

Задание

Описание процесса теплообмена

Обсуждение результатов

Вывод

Список используемых источников

Приложение А

Приложение Б

Задание

Выбрать из ряда типовых теплообменников оптимальный с точки зрения эффективности теплопередачи теплообменник, в котором 100000 кг/ч ацетона меняет свою температуру с t>1н> до 40 ОС за счет теплообмена с 90000 кг/ч дивинила, имеющей начальную температуру 10 ОС и конечную температуру 50 ОС. Определить стоимость теплообменника с учетом факторов удорожания и инсталляции.

Схема процесса теплообмена представлена на рисунке 1.

Рисунок 1 — Схема процесса теплообмена

Данные для расчета (значения теплоемкостей и коэффициенты уравнения Антуана [1]) приведены в таблице 1.

Таблица 1. Данные для расчета

Первый поток

Второй поток

Соединение

ацетон

дивинил

Теплоемкость, кДж/(кг∙К)

2,28

2,65

Коэффициенты уравнения

Антуана

А

16,6513

15,7727

В

2940,46

2142,66

С

–35,93

–34,30

Описание процесса теплообмена

Основную группу теплообменных аппаратов, применяемых в промышленности, составляют поверхностные теплообменники, в которых теплота от горячего теплоносителя передается холодному теплоносителю через разделяющую их стенку. Другую группу составляют теплообменники смешения, в которых теплота передается при непосредственном соприкосновении горячего и холодного теплоносителей.

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

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

Обсуждение результатов

В данной работе будем выполнять расчет теплообменника согласно известной методики [2].

Сначала определим недостающее исходное значение для расчета — это расход второго потока t>1н>.

Количество тепла для первого потока (Q>1>) можно определить по формуле:

Q>1> = ∙ C>1> ∙ Дt>1>, кВт, (1)

G>1> — расход первого потока, кг/ч,

С>1> — теплоемкость первого потока, кДж/(кг∙К),

Дt>1> — разность температур первого потока, ОС.

Аналогично определяется количество тепла для второго потока (Q>2>) по формуле:

Q>2> = ∙ C>2> ∙ Дt>2>, кВт, (2)

G>2> — расход первого потока, кг/ч,

С>2> — теплоемкость первого потока, кДж/(кг∙К).

Дt>2> — разность температур первого потока, ОС.

Определим количество тепла для второго потока (Q>2>):

Q>2> = ∙ 2,65 ∙ (50 – 10) = 2650 кВт.

Из условия равенства количества тепла первого и второго потока (Q>1> = Q>2>) определим недостающее исходное значение для дальнейшего расчета (начальную температуру первого потока t>1н>). Выразим t>2к> из формулы и определим его значение:

t>1н> = + t>1к> , ОС,

t>1н> = + 40 = 81,84 ОС.

Далее рассмотрим вариант теплообмена с прямоточным движением потоков. Температурная схема потоков представлена на рисунке 2.

Рисунок 2 — Температурная схема движения потоков при прямотоке

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

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

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

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

Конструктивно кожухо-трубчатые теплообменники внутри имеют трубное и межтрубное пространство, куда соответственно подается один и другой теплоноситель. Таким образом, возникает необходимость рассмотреть два варианта: первый вариант – когда горячий поток подается в трубное пространство, а холодный поток в межтрубное пространство, и второй вариант – когда наоборот, холодный поток подается в трубное пространство, а горячий поток в межтрубное.

Рисунок 3 — Температурная схема движения потоков при противотоке

Рассмотрим вариант, когда в трубное пространство подается горячий поток (G>1>), а в межтрубное пространство холодный поток (G>2>).

Для расчета стоимости теплообменника необходимо знать величину площади поверхности теплообмена. Однако для этого необходимо знать среднюю разность температур в теплообменнике Дt>ср>, которая определяется в зависимости от значения отношения , где Дt> – большая разность температур в теплообменнике, Дt> – меньшая разность температур.

Если отношение ≤ 2, то среднюю разность определяют как среднее арифметическое:

Дt>ср> = . (3)

Если отношение > 2, то среднюю разность определяют как среднее логарифмическое:

Дt>ср> = . (4)

Определим среднюю разность температур.

Дt> = t>1н> – t>2к> = 81,84 – 50 = 31,84 ОС,

Дt> = t>1к> – t>2н> = 40 – 10 = 30 ОС,

= = 1,06 < 2, тогда Дt>ср> определяется как среднее арифметическое,

Дt>ср> = = = 30,92 ОС.

В зависимости от процесса протекающего в аппарате, принимаем значения коэффициента теплопередачи K [1]. В теплообменнике тепло передается от одной органической жидкости к другой, поэтому принимаем значение коэффициента теплообмена равное 200 Вт/(м2ОС).

Теперь, зная среднюю разность температур в теплообменнике, коэффициент теплопередачи, определяем значение площади поверхности теплообмена (F):

F = , м2, (5)

F = = 429 м2.

Определенное значение площади поверхности является теоретическим. На практике имеется ряд стандартных теплообменников с определенными площадями теплообмена и другими параметрами. Выбираем из каталога стандартных теплообменников [1] ближний больший по площади теплообменник. Это кожухотрубчатый теплообменник, имеющий следующие параметры и их значения:

– площадь поверхности теплообмена – F = 444 м2;

– длина труб – L = 6 м;

– диаметр кожуха – D> = 1 м;

– наружный диаметр труб – d> = 0,02 м;

– число ходов – z = 1;

– внутренний диаметр труб – d>вн> = 0,016 м;

– количество труб – n = 1178.

Теперь необходимо рассчитать уточненное значение коэффициента теплопередачи (К>). Оно необходимо для того, что бы в дальнейшем проверить, подходит ли нам выбранный из каталога теплообменник. Уточненное значение коэффициента теплопередачи можно определить по формуле:

К> = . (6)

Здесь:

б>тр> и б>мтр> — коэффициенты теплоотдачи для трубного и межтрубного пространства теплообменника соответственно;

У — сумма коэффициентов загрязнения, зависящая от среды (справочное значение [1]).

Коэффициенты теплопроводности б>тр> и б>мтр> определяются в зависимости от режима движения потока. Режим движения определяется критерием Рейнольдса (Re>i>):

Re>i> = , (7)

где: w>i> — скорость движения i-го потока;

d>i> — диаметр i-го сечения;

с>i> — плотность i-го потока (из справочных материалов [1]);

м>i> — вязкость i-го потока.

Скорость движения потока определяется по формуле:

w>i> = , м/с, (8)

где: V>i> — объемный расход i-го потока,

S>сеч_>>i> — площадь сечения трубного или межтрубного пространства для соответствующего потока.

Для трубного и межтрубного пространства площадь сечения определяется по формулам (9) и (10):

S>сеч тр> = ∙ n , м2, (9)

S>сеч тр> = – ∙ n , м2, (10)

S>сеч тр> = ∙ 1178 = 0,2367 м2,

S>сеч тр> = – ∙ 1178 = 0,4151 м2.

Объемный расход i-го потока определяется следующей формулой:

V>i> = , м3/с. (11)

Для горячего потока V>1> = = 0,0352 м3/с,

а для холодного потока V>2> = = 0,0403 м3/с.

Вязкость потока (м>i>) можно определить по выражению:

м>i> = , Па∙с (12)

Здесь:

T — средняя температура потока, К,

VISB и VISTO — параметры для расчета вязкости (справочные данные [3]).

Для горячего потока (ацетон): Т=333,92К; VISB=367,25; VISTO=209,68. Для горячего потока (дивинил): Т=303К; VISB=300,59; VISTO = 163,12.

Определим вязкости горячего и холодного потоков соответственно:

м>1> = = 0,000223 Па∙с,

м>2> = = 0,000141 Па∙с.

Определим скорости движения потоков для горячего и холодного соответственно:

w>1> = = 0,1485 м/с,

w>2> = = 0,0970 м/с.

Определим критерий Рейнольдса для горячего и холодного потоков соответственно:

Re>1> = = 8418,41 (переходный режим),

Re>2> = = 8541,15 (переходный режим).

При расчете коэффициентов теплопроводности потоков имеет место значение критерия Прандтля (Pr>i>), определяемый выражением:

Pr>i> = (13)

л>i>—коэффициент теплопроводности (из справочных материалов [1]), Вт/м∙К.

Определим этот критерий для обоих потоков:

Pr>1> = = 3,63;

Pr>2> = = 3,40.

Коэффициент теплопроводности потока для трубного пространства (б>тр>) при переходном режиме движения потока (1000<Re>i><10000) определяется формулой:

б>тр> = ∙ 0,008 ∙ Re>i>0,9 ∙ Pr>i>0,430,25 . (14)

Выражение 0,25 как правило равно 1, поэтому обычно оно не учитывается.

Найдем значение коэффициента теплопроводности для трубного пространства:

б>тр> = ∙ 0,008 ∙ 8418,410,8 ∙ 3,630,43 = 415,59.

Коэффициент теплопроводности потока для межтрубного пространства (б>мтр>) при Re>i> > 1000 определяется формулой:

б>мтр> = ∙ 0,4 ∙ Re>i>0,6 ∙ Pr>i>0,360,25 ∙ 0,85. (15)

Определяем значение коэффициента (выражение 0,25 также принимаем равным 1):

б>мтр> = ∙ 0,4 ∙ 8541,150,6 ∙ 3,400,36 ∙ 0,85 = 663,69.

Определяем сумму коэффициентов загрязнения на стенках трубопроводов У [1], для данного процесса У = 0,0008 м2ОС/Вт.

Теперь можно определить уточненное значение коэффициента теплопередачи (K>):

К> = = 212,18.

Зная уточненное значение коэффициента теплопередачи, определяем уточненное значение площади поверхности теплообмена (F>):

F> = , м2, (16)

F> = = 404 м2.

Теперь можем определить запас поверхности (Zp) теплообмена, тем самым узнать, подходит ли нам выбранный теплообменник:

Zp = ∙ 100, %, (17)

Zp = ∙ 100 = 9,03 %.

Подходящим принято считать теплообменник, имеющий запас поверхности 15—30 % [2]. Полученный запас поверхности 9,03 % недостаточен. Поэтому следует выбрать из каталога стандартных теплообменников другой теплообменник, больший по площади, и провести аналогичный расчет.

В качестве принимаемого значения коэффициента теплопередачи берем значение 212,18, т.е. уточненное значение из предыдущего расчета. Определено расчетное значение площади теплообмена:

F = , м2,

F> = = 404 м2.

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

– площадь поверхности теплообмена – F = 416 м2;

– длина труб – L = 9 м;

– диаметр кожуха – D> = 0,8 м;

– наружный диаметр труб – d> = 0,020 м;

– число ходов – z = 1;

– внутренний диаметр труб – d>вн> = 0,016 м;

– количество труб – n = 736.

Рассчитан уточненный коэффициент теплопередачи, К>=284,02; запас поверхности составил 27,46 %. Это удовлетворяющий запас поверхности. Поэтому, для варианта подачи теплоносителей "горячий поток в трубное пространство, холодный поток в межтрубное пространство", наиболее удовлетворителен последний теплообменник. Расчет представлен в приложении А.

Для другого варианта подачи теплоносителей "холодный поток в трубное пространство, горячий поток в межтрубное пространство" проведен аналогичный расчет. В качестве принимаемого значения коэффициента теплопередачи взято значение 212,18, т.е. уточненное значение из первого расчета. Определено расчетное значение площади теплообмена:

F = , м2,

F> = = 404 м2.

Из каталога теплообменников выбран стандартный кожухотрубчатый теплообменник, удовлетворяющий нашим условиям, имеющий следующие параметры и их значения:

– площадь поверхности теплообмена – F = 416 м2;

– длина труб – L = 9 м;

– диаметр кожуха – D> = 0,8 м;

– наружный диаметр труб – d> = 0,020 м;

– число ходов – z = 1;

– внутренний диаметр труб – d>вн> = 0,016 м;

– количество труб – n = 736.

Рассчитан уточненный коэффициент теплопередачи, К>=292,92; запас поверхности составил 29 %. Этот теплообменник удовлетворяет рекомендуемому запасу поверхности. Расчет представлен в приложении Б.

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

Определим стоимость теплообменника (Ст>то>):

Ст>то> = (6 + 0,075 ∙ F) ∙ 103 , Ј, (18)

Ст>то> = (6 + 0,075 ∙ 416) ∙ 103 = 37200 Ј.

Следует заметить, что теплообменник работает под определенным давлением, зависящим от состава и температуры потоков. Поэтому рассчитанную стоимость необходимо домножить на так называемый "Фактор удорожания" (F>удор>), учитывающий давление в аппарате.

Определим максимально возможное давления в теплообменнике. Для этого используем уравнение Антуана:

P>i> = , МПа, (19)

где:

А>i>, В>i>, С>i> — значения коэффициентов уравнения Антуана для i-го вещества [1],

T>i> — температура i-го вещества, К.

Через теплообменник проходят два потока, две органические жидкости: ацетон и дивинил. Давление для первого потока, ацетона:

P>ацетон> = = 0,222 МПа.

Давление для второго потока, дивинила:

P>дивинил> = =0,557 МПа.

Из справочного материала для большего давления определяем значение фактора удорожания F>удор>, оно составит 2,0.

Далее определяется стоимость аппарата с учетом фактора удорожания (Ст>то >>F>>удор>):

Ст>то >>F>>удор> = Ст>то> ∙ F>удор А> , Ј, (20)

Ст>то >>F>>удор> = 37200 ∙ 2,0 = 74400 Ј.

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

Общий фактор инсталляции F>инст> определяется по формуле:

F>инст> = 1 + f>уст> + f>труб> + f>инстр> + f>электр> + f>строит> + f>зданий> + f>изол> , (21)

где:

f>уст> – подфактор, учитывающий дополнительные затраты при установке оборудования;

f>труб> – подфактор, учитывающий дополнительные затраты при прокладке трубопроводов;

f>инстр> – подфактор, учитывающий дополнительные затраты на средства контроля и автоматики;

f>электр> – подфактор, учитывающий дополнительные затраты на проводку электроэнергии;

f>строит>, f>зданий> – подфакторы, учитывающий дополнительные затраты на гражданское и специальное строительство;

f>изол> – подфактор, учитывающий дополнительные затраты на изоляцию оборудования;

Инсталляционные подфакторы определяются по справочным данным в зависимости от стоимости оборудования, рассчитанной с учетом фактора удорожания. Определим инсталляционные подфакторы.

Принимаем к установке оборудование, требующее ряд технологических работ по месту установки. Этому соответствует инсталляционный подфактор f>уст> = 0,10.

Для нашего оборудования нам необходимы средние технические и сервисные трубопроводы, чему соответствует инсталляционный подфактор f>тр> = 0,26.

В теплообменнике протекает процесс, в котором необходимо контролировать и регулировать заданные значения конкретных параметров. Для этого необходима автоматизация процесса. Предполагается установка приборов контроля и одного микроконтроллера для управления процессом. Имеет место подфактор f>инст> = 0,13.

Подфактор, учитывающий электроэнергию, инсталляционный подфактор f>электр> = 0,03., т.к. предполагается только освещение и малое потребление электроэнергии системой управления, что ничтожно мало по сравнению с освещением.

Необходимы средние строительные штатные работы (фундамент, 1-й этаж, этажерки и службы), инсталляционный подфактор f>строит> = 0,10. Оборудование предполагается разместить под крышей типа ангара, инсталляционный подфактор f>зданий> = 0,29.

Поскольку теплообменник является термооборудованием, поэтому для снижения потерь предполагается защитная изоляция уровня выше среднего, инсталляционный подфактор f>изол> = 0,06.

F>инст> = 1 + 0,1 + 0,26 + 0,13 + 0,03 + 0,1 + 0,29 + 0,06 = 1,97.

Теперь можно определить полную (окончательную) стоимость теплообменника, т.е. стоимость, рассчитанную с учетом фактора инсталляции (Ст>то полная>):

Ст>то полная> = Ст>то >>F>>удор> ∙ F>инст> , Ј, (22)

Ст>то полная> = 74400 ∙ 1,97 = 146568 Ј.

Вывод

В результате проведенной расчетно–аналитической работы из ряда типовых теплообменников выбран оптимальный с точки зрения эффективности теплопередачи теплообменник, удовлетворяющий поставленным условиям, имеющий следующие параметры и их значения:

– площадь поверхности теплообмена – F = 416 м2;

– длина труб – L = 9 м;

– диаметр кожуха – D> = 0,8 м;

– наружный диаметр труб – d> = 0,020 м;

– число ходов – z = 1;

– внутренний диаметр труб – d>вн> = 0,016 м;

– количество труб – n = 736.

Коэффициент теплопередачи равен 284,02; запас поверхности составил 27,46 %. Стоимость теплообменника с учетом факторов удорожания и инсталляции составила 146568 Ј.

Список используемых источников

1.Холодов «Математическое моделирование и оптимизация химико-технологических процессов» Издательство: ООО "Издательский дом "Профессионал", 2003 г. 480 стр.

2.«Математическое моделирование химико-технологических систем» Д.В. Саулин

3.www.XuMuK.ru

4.http://www.chemport.ru/

5.http://www.chmm.spb.ru/lectures.php?type=sysanalisys

6.Закгейм А.Ю. «Введение в моделирование химико-технологических процессов» Издание 2 

Приложение А

Первая жидкость

Вторая жидкость

Ацетон

Дивинил

Расход

100 000,00

90 000,00

кг/ч

Темпер нач

81,84

10,00

оС

Темпер кон

40,00

50,00

оС

C

2,28

2,65

кДж/(кг К)

K

212,18

Вт/(м2 °С)

Плотность с

790,00

621,00

кг/м3

VISB

367,25

300,59

VISTO

209,68

163,12

Коэфф т/провод л

0,14

0,11

Вт/(м К)

Количество тепла 1 потока Q1

2 650,00

кВт

Количество тепла 2 потока Q2

2 650,00

кВт

Противоток

Дtб

30,000

Дtм

31,842

Дtб/Дtм

0,942

Дtср

30,921

F

403,912

выбираем теплообменник

F (принятая)

416,000

м2

Поверх теплообмена

L

9,000

м

Длина труб

0,800

м

Диаметр кожуха

0,020

м

Наруж диам труб

z

1,000

Число ходов

dВН

0,016

м

Внутр диам труб

n

736,000

Число труб

Горячий поток в трубном пространстве, а холодный в межтрубном

Об расх V1

0,0352

м3/сек

Об расх V2

0,0403

м3/сек

Площ сеч труб пр-ва

0,1479

м2

Площ сеч м/тр пр-ва

0,2713

м2

Скорость W1 (тр)

0,2377

м/сек

Скорость W2 (м/тр)

0,1484

м/сек

Сред темп потока 1

60,921

°С

333,921

Сред темп потока 2

30,000

оС

303,000

Вязкость м1

0,000223

Па/с

Вязкость м2

0,000141

Па/с

Re1

13474,034

турбулентный

Re2

13068,755

турбулентный

Pr1

3,631945

Pr2

3,397391

б тр

643,703

б м/тр

856,632

Сумма коэфф загряз

0,0008

Коэфф К

284,0203895

Уточненная поверхность теплообмена

F уточненная

301,746392

Запас поверхности

27,46480972

%

Приложение Б

Первая жидкость

Вторая жидкость

Ацетон

Дивинил

Расход

100 000,00

90 000,00

кг/ч

Темпер нач

81,84

10,00

оС

Темпер кон

40,00

50,00

оС

C

2,28

2,65

кДж/(кг К)

K

212,18

Вт/(м2 °С)

Плотность с

790,00

621,00

кг/м3

VISB

367,25

300,59

VISTO

209,68

163,12

Коэфф т/провод л

0,14

0,11

Вт/(м К)

Количество тепла 1 потока Q1

2 650,00

кВт

Количество тепла 2 потока Q2

2 650,00

кВт

Противоток

Дtб

30,000

Дtм

31,842

Дtб/Дtм

0,942

Дtср

30,921

F

403,912

выбираем теплообменник

F (принятая)

416,000

м2

Поверх теплообмена

L

9,000

м

Длина труб

0,800

м

Диаметр кожуха

0,020

м

Наруж диам труб

z

1,000

Число ходов

dВН

0,016

м

Внутр диам труб

n

736,000

Число труб

Холодный поток в трубном пространстве, а горячий в межтрубном

Об расх V1

0,0403

м3/сек

Об расх V2

0,0352

м3/сек

Площ сеч труб пр-ва

0,1479

м2

Площ сеч м/тр пр-ва

0,2713

м2

Скорость W1 (тр)

0,2722

м/сек

Скорость W2 (м/тр)

0,1296

м/сек

Сред темп потока 1

60,921

°С

333,921

Сред темп потока 2

30,000

оС

303,000

Вязкость м1

0,000223

Па/с

Вязкость м2

0,000141

Па/с

Re1 труб

19176,977

турбулентный

Re2 м/тр

9182,304

переходный

Pr1

3,397391

Pr2

3,631945

б тр

829,553

б м/тр

710,007

Сумма коэфф загряз

0,0008

Коэфф К

292,9196809

Уточненная поверхность теплообмена

F уточненная

292,578933

Запас поверхности

29,66852574

%