CC BY
17
УДК 621.184
СТРУКТУРНО-ПАРАМЕТРИЧЕСКАЯ ОПТИМИЗАЦИЯ МНОГОПОТОЧНЫХ СИСТЕМ ПЛАСТИНЧАТЫХ ТЕПЛООБМЕННЫХ АППАРАТОВ
ЖУКОВ В.П., д-р техн. наук, БАРОЧКИН Е.В., канд. техн. наук
Разработана математическая модель многопоточных пластинчатых теплообменных аппаратов, сформулирована задача выбора оптимальной структуры потоков и параметров установки.
Ключевые слова: пластинчатые теплообменные аппараты, математическая модель, теплоноситель.
THE STRUCTURAL-PARAMETRIC MULTITHREAD SYSTEM OPTIMIZATION OF PLATE HEAT-EXCHANGE UNITS
ZHUKOV V.P., Ph.D., BAROCHKIN E.V., Ph.D.
The article contains a developed mathematical model of multithread plate heat-exchange units, a formulated problem of selection of optimum flow structure and installation parameters.
Key words: plate heat-exchange units, mathematical model, heat carrying agent.
В качестве объекта моделирования выбрана многопоточная многоступенчатая система пластинчатых теплообменных аппаратов. Термин «многопоточная» в данном случае понимается как возможность распределенной подачи произвольного числа исходных теплоносителей в систему, а также возможность деления на выходе и объединения на входе ступени произвольного числа потоков. Термин «многоступенчатая» означает наличие в системе произвольного числа ступеней с произвольным числом связей между ними.
Матричная формализация расчета теплообменных аппаратов со сложной конфигурацией потоков [1] предполагает представление каждой ступени или аппарата четырехполюсником с двумя входными и двумя выходными потоками для холодного и горячего теплоносителя, соответственно. В пластинчатом теплообменнике каждый теплоноситель в ступени взаимодействует с двумя соседними теплоносителями (рис. 1,а), поэтому метод матричной формализации [1] для пластинчатых аппаратов нуждается в корректировке. Для адаптации метода принимается следующее допущение: каждая ступень вдоль направления движения теплоносителя условно разбивается на две теплоизолированные между собой части, в каждой из которых теплоноситель контактирует только с одним из соседних теплоносителей. Данный расчетный прием позволяет представить ступень в виде двух четырехполюсников. Выделенные потоки перемешиваются на выходе из ступени и направляются в следующую ступень.
Порядок составления расчетной схемы и матричной модели демонстрируется на примере пластинчатого теплообменника (рис. 1). В крайних секциях теплообмен осуществляется только через одну стенку, что позволяет при составлении расчетной схемы ограничиться тремя четырехполюсниками. Однако для удобства выполнения расчетного анализа и наглядности представления результатов численных экспериментов вводится дополнительный четвертый элемент, входные параметры которого соответствуют выходным параметрам рассматриваемой установки. Таким образом, для составления расчетной схемы и модели предлагается схема из четырех элементов (рис. 1,б).
Рис. 1. Схема потоков в четырехступенчатом пластинчатом теплообменнике: а - технологическая; б - расчетная (сплошными и штриховыми линиями показаны соответственно потоки горячего и холодного теплоносителя)
Расчет пластинчатого теплообменника выполняется в два этапа: сначала определяются расходы теплоносителей через каждую ступень, затем вычисляются значения температуры теплоносителя в произвольной точке аппарата.
При выполнении первого этапа в качестве параметров потоков Х выбираются расходы горяче-
го и холодного теплоносителей: X =
. Потоки в
поверхностном подогревателе разделены стенкой, поэтому их расходы не меняются, а матрица процесса в ступени соответствует единичной матрице 1 0",
Матрица коммутации состоит из четы-
B ,
,0 1
рех элементов Kjj =
КХц
«21
а12
а22
каждый из кото-
рых показывает долю потока из ]-го в 1-й элемент: а11 - горячего теплоносителя в горячий; а12 -
холодного в горячий; а21 - горячего в холодный;
а22 - холодного в холодный. Согласно расчетной
схеме (рис. 1,б), матрицы коммутации записываются следующим образом:
= к12 = к13 = к14 = к22 = к23 = к24 = к
34
12
к
41
43
к
44
к
31
0.5 0
, к
14 22
0 0 , 0 0 1 0
к
, к
21
24 32
0.5 0 0 0 10
0 1/' 42 ” ^0 0^г'43_ ^0 1
На первом этапе расчета матричная модель принимает вид
л
33
к
43
кц -1 к12 *13 к14 л 1х л -[Хо ]і
к21 кц -1 3 сч к 4 2 к Х2 -[Хо ]2
*31 3 к кц -1 4 3 к Х3 -[Хо ]3
к41 2 4 к 3 4 к кц - I / 1Х4 / 1-[Хо ]4
(1)
Из решения (1) находятся расходы теплоносителей через ступени установки. Следует отметить, что предложенный подход справедлив и при наличии внутренних рециркуляций, рассчитать которые посредством последовательных арифметических вычислений не всегда возможно.
На втором этапе расчета вектор признаков составляется из температуры горячего и холодного (I, ^
теплоносителя X = I ^ I. Матрица процесса в ступени В для каждого четырехполюсника, согласно [2], записывается в следующем виде:
В
Ь11 -
Ь11 Ь12
Ь21 Ь22,
( а2 — |
1 + е 1
1 а1
(а1+а2)^
Ь21 =
^2 ^2
— —2е
^1 ^1
а2
(а1+а2)Р
Ь22 =| 1 + 1^е-(а1+а2)р
0 02
- диагональная матрица расходов
теплоносителей через ступень; [0Хо ] =
01
01
- век-
тор-столбец из произведений температуры на расход внешних потоков.
На втором этапе решения задачи вектор внешних входных потоков записывается в виде
Х0 = -[С10110 ^20^20 ^2 0 ^20^20 ^2 ООО 0]'.
Решение системы (2) позволяет определить значение температуры теплоносителей в произвольной точке установки.
Рассматривается пример решения задачи структурно-параметрической оптимизации многопоточной системы пластинчатых теплообменных аппаратов. Установка состоит из трех пластинчатых теплообменных аппаратов, включенных параллельно по горячему и холодному теплоносителю. В качестве параметров оптимизации выбираются доли горячего х1 и холодного х2 теплоносителя, подаваемые на вход в первый подогреватель. При этом доли внешних потоков, направляемые во второй и третий подогреватели, принимаются равными. В качестве целевой функции оптимизации выбирается температура холодного теплоносителя на выходе установки. Необходимо определить распределение внешних потоков по ступеням, которое обеспечивает максимальную температуру холодного теплоносителя на выходе, то есть максимальную эффективность теплообмена. Температура и суммарный расход входных потоков в установку при этом не изменяются.
Результаты решения сформулированной задачи представлены в виде зависимостей температуры холодного теплоносителя на выходе от доли горячего теплоносителя х1 (рис. 2), анализ которых показывает, что максимальное значение целевой функции соответствует равномерному распределению потоков холодного и горячего теплоносителя по ступеням.
а1 =-
а2 = ■
С-|0-| О202
где к - коэффициент теплопередачи; с - удельная теплоемкость; О - расход; 1 - температура теплоносителя; индекс «1» относится к горячему, «2» - к холодному теплоносителю.
Температура не является аддитивным параметром, поэтому для составления балансовых уравнений температура умножается на расход теплоносителя и теплоемкость. Полученное произведение есть энергия (или энтальпия) потока и при смешении потоков их энергии могут складываться. С учетом сделанных замечаний матричная модель принимает вид
Э11
к12022В2 -022 к32022В2
к21011В1 к31011В1 к 41011В1 к42022В2
к13033В3 к23033В3 -033 к43033В3
к14044В4 ^
к24044В4
к34044В4
44
(Хі л [0Хо ]і
Х2 -[0Х о ]2
Х3 -[0Хо ]3
1Х 4 / 1-[0Х о ]4
(2)
Рис. 2. Зависимость температуры холодного теплоносителя на выходе от доли потока горячего теплоносителя в первый теплообменник при различных значениях доли холодного теплоносителя, подаваемого в первый теплообменник: 1 - х2 = 0,01; 2 - 0,33; 3 - 0,99
Результаты исследований могут использоваться для постановки и решения оптимизационных задач при наладке и проектировании пластинчатых теплообменных установок.
X
Список литературы
1. Моделирование тепломассообмена в струйных деаэраторах со сложной конфигурацией потоков / Е.В. Ба-рочкин, В.П. Жуков, Г.В. Ледуховский и др. // Известия вузов. Химия и химическая технология. - 2004. - Т. 47. -Вып. 9.
2. Барочкин Е.В., Жуков В.П., Ледуховский Г.В.
Обобщенная модель каскадных теплообменных аппаратов с учетом фазовых переходов // Известия вузов. Химия и химическая технология. - 2004. - Т. 47. - Вып. 3.
Жуков Владимир Павлович,
ГОУВПО «Ивановский государственный энергетический университет имени В.И. Ленина», доктор технических наук, профессор кафедры прикладной математики, телефон (4932) 26-97-45, e-mail: [email protected]
Барочкин Евгений Витальевич,
ГОУВПО «Ивановский государственный энергетический университет имени В.И. Ленина», кандидат технических наук, доцент кафедры тепловых электрических станций, телефон (4932) 41-60-56, e-mail: [email protected]
