1. Введение
Актуальной задачей для многих отраслей современной промышленности является совершенствование теплообменных аппаратов, среди которых распространены конструкции типа «труба в трубе». Теплообменные аппараты данного типа отличаются простотой конструкции и надежностью, а возможность их эксплуатации в условиях ледообразования позволяет применять их в качестве кристаллизаторов в криоконцентраторах и ледогенераторах.
При предложении новых конструкций аппаратов необходимо представить обоснование, по каким именно показателям эффективности их работы новые конструкции будут более совершенными в сравнении с известными. Различие в подходах к оценке эффективности работы теплообменных аппаратов делает актуальной разработку комплексной сравнительной оценки, включающей ряд наиболее важных характеристик аппаратов, таких как материалоемкость, компактность, тепловая мощность, потери мощности при эксплуатации аппарата.
Прогнозирование работы новых конструкций аппаратов может быть основано на уточненных уравнениях Навье-Стокса, определенных из массива экспериментальных данных.
Рассмотрим вертикальный кристаллизатор ёмкостного типа (концентрирование продукта на внутренней поверхности теплообменной трубы), приняв конструкцию
типа «труба в трубе» (рисунок 1) [1], [2], [3].
Вертикальный кристаллизатор конструкции типа «труба в трубе»
Такая конструкция может использоваться для концентрирования исходного жидкого продукта путем намораживания льда на внутренней поверхности теплообменной трубы (периодическое действие), а также для предварительного охлаждения среды, движущейся внутри теплообменной трубы (непрерывное действие).
Теплообмен в данном кристаллизаторе осуществляется за счет контакта исходного жидкого продукта с хладоносителем через поверхность теплообменной трубы 2. При этом исходный жидкий продукт движется внутри теплообменной трубы 2, а хладоноситель по межтрубному пространству, образованному теплообменной трубой 2 и корпусом 1. Корпус 1, в свою очередь, оснащен патрубками для входа 3 и выхода 4 хладоносителя.
Ключевое отличие кристаллизаторов от других рекуперативных теплообменников — их циклический принцип действия, состоящий из двух основных этапов. Первый этап — это охлаждение исходной жидкости ниже ее точки замерзания с помощью хладоносителя для кристаллизации. Второй этап — это отделение полученного льда и концентрированного продукта. Для отделения льда могут применяться механические средства (например, вибрация или резка) или тепловой метод, при котором используется теплоноситель с температурой, равной или выше температуры плавления льда.
Для оценки эффективности работы теплообменных аппаратов (включая кристаллизаторы) используются как натурные эксперименты (включая испытания на действующих предприятиях), так и компьютерное моделирование
[4][5][6]
Таким образом, можно провести сравнительную оценку работы типовой конструкции кристаллизатора
[7] и ряда новых конструкций при режимах работы без фазового перехода хладоносителя [9][10][11]
В настоящее время компьютерное моделирование широко используется для изучения процесса теплообмена, включая исследования с сопоставлением результатов моделирования с результатами натурных экспериментов.
В отечественных
[12][13][14][15][16][17][18][19][20][21][22][23]
2. Методы и принципы исследования
Выполним кристаллизатор согласно
нормативному документу [7] со следующими размерами:
Проведем компьютерное моделирование работы кристаллизатора в программном комплексе "SolidWorks", приняв
исходные данные хладоноситель этиленгликоль с массовой концентрацией 20% и температурой на входе в кристаллизатор +1 внутри теплообменной трубы движется вода с температурой на входе +10 ℃. Скорость воды на входе примем равной 0,3 м/с, а скорость хладоносителя будет меняться от 1,5 м/с до 3,9 м/с (что соответствует диапазону значений критерия Рейнольдса 4000–10000), с шагом 0,3 м/с.
При расчете и проектировании рекуперативных теплообменных аппаратов и, в частности, кристаллизаторов одним из определяющих параметров для оценки их эффективности, вне зависимости от области применения, является теплообменная мощность
[8]
где Q теплообменная мощность кристаллизатора, Вт; K коэффициент теплопередачи, Вт/(м2‧К); F площадь поверхности теплообмена, м2; Δt средняя разность температур между исходным жидким продуктом и хладоносителем, К.
Особое внимание уделяется увеличению данной характеристики при сохранении габаритных размеров кристаллизатора
[2][3]
В патентах
[9][10][11]
Эпициклоиды с различным количеством каспов
[24]
Изменение площади поверхности теплообмена трубы-эпициклоиды в зависимости от количества ее каспов, относительно трубы круглого сечения
| Количество каспов эпициклоиды | 3 | 4 | 5 | 6 |
| Увеличение площади поверхности теплообмена, % | +1,9 | +6,2 | +9,2 | +11,5 |
Трубы такой формы могут быть изготовлены путем
применения валковой формовки (профилирования) [25], используя в качестве заготовки трубу круглого сечения.
К недостаткам данного технического решения можно отнести сложность изготовления труб-эпициклоид, что неизбежно приведет к увеличению стоимости кристаллизатора.
Проанализируем конструкции кристаллизаторов с теплообменными трубами-эпициклоидами с количеством каспов от 1 до 10. Примем для их моделирования те же исходные данные и габаритные размеры, что и для типовой конструкции.
Поставленная в работе цель — разработка новых конструкций оборудования холодильной техники, кристаллизаторов, и оценка рациональности их последующего внедрения для проведения опытно-промышленных испытаний на основании сравнения как между собой, так и с типовой конструкцией по показателям энерго- и ресурсосбережения — является актуальной для современной промышленности. Выбранные методы и принципы исследования соответствуют современной мировой и отечественной практике.
3. Основные результаты
По итогам моделирования получены данные, позволяющие определить эмпирические коэффициенты в уравнении определения коэффициента теплоотдачи межтрубного пространства каждой из рассматриваемых конструкций кристаллизаторов в диапазоне значений критерия Рейнольдса 4000-10000. Уточненные уравнения для кристаллизаторов с разным сечением теплообменной трубы перечислены ниже.
Круглое сечение (типовая конструкция):
Эпициклоида с 1 каспом:
Эпициклоида с 2 каспами:
Эпициклоида с 3 каспами:
Эпициклоида с 4 каспами:
Эпициклоида с 5 каспами:
Эпициклоида с 6 каспами:
Эпициклоида с 7 каспами:
Эпициклоида с 8 каспами:
Эпициклоида с 9 каспами:
Эпициклоида с 10 каспами:
где
Nu критерий Нуссельта; Re критерий Рейнольдса; Pr критерий Прандтля.
Представленные уточненные уравнения определения коэффициента теплоотдачи межтрубного пространства позволяют прогнозировать работу кристаллизаторов рассматриваемых конструкций, с учетом физико-химических свойств хладоносителя в межтрубном пространстве, его температуры и скорости на входе в кристаллизатор.
Для оценки эффективности кристаллизаторов предлагается использование ряда удельных показателей.
Отношение теплообменной мощности кристаллизатора к его массе:
где
mу удельная масса кристаллизатора Вт/кг; m – масса кристаллизатора кг.
Отношение теплообменной мощности кристаллизатора к его габаритным размерам:
где Оу компактность кристаллизатора, Вт/м3; О – габаритные размеры кристаллизатора, м3.
Отношение теплообменной мощности кристаллизатора к суммарным потерям мощности на перекачивание его рабочих сред:
где Еу удельный показатель суммарных потерь мощности к количеству передаваемого тепла при эксплуатации кристаллизатора, Вт/Вт; Е суммарные потери мощности на перекачивание рабочих сред кристаллизатора, Вт.
Отношение теплообменной мощности кристаллизатора к объему его внутритрубного пространства (максимальному объему исходного жидкого продукта):
где ОQ отношение теплообменной мощности кристаллизатора к объему его внутритрубного пространства, Вт/м3; ОТР объем внутритрубного пространства кристаллизатора, м3.
Отношение суммарных потерь мощности на перекачивание рабочих сред кристаллизатора к объему его внутритрубного пространства:
где ОЕ отношение суммарных потерь мощности при эксплуатации кристаллизатора к объему его внутритрубного пространства, Вт/м3.
Для комплексной оценки кристаллизаторов различных конструкций предлагается разработанная методика сравнительной оценки, в основе которой лежат сумма отношений перечисленных удельных показателей.
Суть методики заключается в том, что в типовой конструкции величине каждого удельного показателя присваивается оценка 1,000. Оценки других конструкций кристаллизаторов получаем путем соотношения величин их удельных показателей с величинами соответствующих удельных показателей типовой конструкции. Так оценка будет показывать во сколько раз оцениваемая конструкция рациональнее к использованию в сравнении с типовой.
Оценки по удельным показателям (13)-(16) получаем путем деления величины показателя оцениваемой конструкции на величину соответствующего показателя типовой конструкции. Например, предположим, что у оцениваемой конструкции кристаллизатора
отношение теплообменной мощности к его габаритным размерам равно 1500 Вт/м3, а у типовой конструкции оно равно 1234 Вт/м3. Тогда оценка оцениваемого кристаллизатора составит 1,216, а типового 1,000.
Оценки по удельному показателю (17) получаем путем деления величины данного показателя типовой конструкции на величину этого же показателя оцениваемой конструкции. Для примера: у оцениваемой конструкции
ОЕ ОЕ оцениваемого кристаллизатора составит 0,797, а типового 1,000.
Чем выше оценки кристаллизатора по перечисленным удельным показателям, тем более рационально его использование в данном значении критерия Рейнольдса межтрубного пространства, относительно типовой конструкции.
Сложением всех полученных оценок каждой конструкции кристаллизаторов при каждом конкретном значении
критерия Рейнольдса получаем суммарную сравнительную оценку рациональности применения новых конструкций кристаллизаторов для диапазона критерия Рейнольдса межтрубного пространства.
Величины удельных показателей получены по данным компьютерного моделирования всех рассматриваемых конструкций кристаллизаторов. Суммарные сравнительные оценки представлены на рисунке 3, где n — количество каспов эпициклоиды, образующей форму поперечного сечения теплообменной трубы.
Суммарная сравнительная оценка кристаллизаторов
По итогам суммарной сравнительной оценки можно рекомендовать в качестве замены типового кристаллизатора конструкции с теплообменной трубой-эпициклоидой с 9 и 10 каспами. Можно также выделить диапазон критерия Рейнольдса 7500–8500, в котором может быть рекомендована конструкция кристаллизатора с теплообменной трубой-эпициклоидой с 3 каспами.
В отличие от известных подходов к сравнительной оценке работы холодильного и теплообменного оборудования (например, представленные в работах
[26][27][28]
4. Обсуждение
Впервые представленные уравнения определения коэффициента теплоотдачи межтрубного пространства позволяют прогнозировать работу 10 новых конструкций кристаллизаторов с профильными теплообменными трубами. Представленные конструкции имеют потенциал широкого использования в линиях холодильной техники, например, концентрирования жидких пищевых продуктов и низкотемпературного опреснения.
Представленные результаты компьютерного моделирования позволяют рекомендовать к использованию конструкции кристаллизаторов с теплообменными трубами в форме эпициклоиды с 9 и 10 каспами в диапазоне значений критерия Рейнольдса в межтрубном пространстве 4000–10000. Для диапазона критерия Рейнольдса 7500-8500 может быть использована конструкция кристаллизатора с теплообменной трубой-эпициклоидой с 3 каспами. Данные конструкции обладают повышенной эффективностью, в сравнении с типовой конструкцией кристаллизатора по комплексу удельных показателей энерго- и ресурсосбережения: материалоемкость, компактность, потери мощности на перекачивание хладоносителя, теплообменная мощность и производительность (выраженная через максимальный объем исходного жидкого продукта в кристаллизаторе).
5. Заключение
Представленные в статье результаты анализа ряда перспективных и запатентованных конструкций кристаллизаторов актуальны для разработки нового оборудования холодильной техники химической, пищевой и биотехнологической отраслей промышленности.
При совершенствовании кристаллизаторов необходимо стремиться не к интенсификации процесса теплопередачи или снижению совокупных энергетических затрат при их эксплуатации, а рассматривать совокупность удельных показателей энерго- и ресурсосбережения процесса.
С использованием компьютерного моделирования проведен термодинамический анализ разработанных и запатентованных конструкций кристаллизаторов с теплообменными трубами в форме эпициклоид с количеством каспов от 1 до 10. Впервые представлены данные эффективности их работы, а также проведено сопоставление с типовой конструкцией кристаллизатора.
С применением метода вычислительного моделирования получены уточненные уравнения коэффициента теплоотдачи межтрубного пространства для кристаллизаторов типовой и рассмотренных запатентованных конструкций. Впервые представленные уравнения с уточненными эмпирическими коэффициентами позволяют прогнозировать работу кристаллизаторов и проводить углубленную оценку особенностей проявляющихся физических закономерностей хладоносителей при их эксплуатации.
Предложена методика сравнительной оценки кристаллизаторов, позволяющая определять их наиболее рациональную конструкцию с точки зрения экономичности. В отличие от традиционного похода к оценке эффективности работы теплообменного и холодильного оборудования, предложенная методика позволяет учитывать комплекс наиболее важных показателей энерго- и ресурсосбережения, что необходимо при решении актуальной задачи современной промышленности создании энерго- и ресурсоэффективных предприятий всех отраслей промышленности.