ВведениеЭффективность извлечения ионов тяжелых металлов из водных растворов предложенным серосодержащим сорбентом доказана в [1, 2], аппаратурное оформление технологического процесса его синтеза приводится в [3]. Для синтеза сорбента используется химический емкостной реактор периодического действия с непрерывным перемешиванием ингредиентов. Масштаб производственной установки не отражается на пропорциях ингредиентов, используемых при синтезе сорбента, т. к. их масса повышается согласно масштабу [3]. Изменяется только время перемешивания, нагрева, термостабилизации и охлаждения реакционной смеси.При расчете реактора использовались инженерные методы расчета, поскольку еще в XX в. было доказано, что невозможно выполнить масштабный переход, используя физическое подобие, из-за несовместности одновременного моделирования массообменных и химических процессов [4, 5]. Актуальность работы обосновывается важностьюи необходимостью автоматизированного расчета реактора, обеспечивающего сокращение трудоемкости, времени разработки проекта установки, производственных затрат и себестоимости продукта [6, 7].  Алгоритм теплового расчета стадии термостабилизации реакционной смесиНа рис. 1 приведена блок-схема алгоритма расчета реактора, исполненного из стали 20, со стандартным эллиптическим днищем. Для перемешивания реакционной смеси используется пропеллерная мешалка [3]. Рис. 1. Алгоритм расчета реактора с пропеллерной мешалкой Fig. 1. Calculation algorithm of a propeller stirrer reactor  Конструктивные размеры реактора приведены в [10] и табл. 1. Таблица 1Table 1Конструктивные размеры реактора с рубашкойStructural dimensions of a jacket reactorПоказатель ОбозначениеЗначениеВнутренний диаметр реактора, мD1,20 Толщина стенки реактора, мs0,006 Наружный диаметр реактора, мD11,212 Длина цилиндрической части реактора, мLц1,214 Объем цилиндрической части реактора, м3Vц1,374 Объем реактора, м3Vн1,6 Площадь поверхности эллиптического днища реактора, м2Fд1,786Объем эллиптического днища реактора, м3Vд0,226Внутренний диаметр цилиндрической обечайки рубашки, мD21,300Внешний диаметр рубашки, м D31,308Объем донной части рубашки, м3Vруб.д0,054Полный объем рубашки, м3 Vруб0,265Площадь теплообмена рубашки с реактором при охлаждении смеси, м2Fохл6,444Площадь рубашки при термостабилизации смеси, м2Fруб7,054  Этап формирования полимерной молекулы сорбента характеризуется изотермичностью реакционной зоны при t = 45 oC (318 К). Для обеспечения изотермичности вода в рубашке должна иметь среднюю температуру, равную этому значению. При этом на внешней стенке рубашки будут иметь место теплопотери в окружающую среду (ОС). Исходя из сказанного, вода изменяет температуру от (318 + х) до (318 – х) К, где x – отклонение температуры относительно среднего значения, являющегося целевым при термостабилизации реакционной смеси. Таким образом, технология термостабилизации заключается в регулировке температуры воды, поступающей в рубашку, на уровне (318 + х) К при заданном ее расходе.Примем температуру ОС 20 оС (293 К). Оценим значение х при скорости теплоносителя W = 0,001 м/с. Средняя разность температур ΔTср определяется по формуле среднелогарифмической разности температур [11]                                                       (1)где ΔT1 – разность температур теплоносителя на входе в рубашку (318 + х) и ОС (293), К; ΔT2 – разность температур теплоносителя на выходе из рубашки (318 – х) и температуры ОС:При этом тепловые потери через стенку рубашки в окружающую среду Q1, Дж, составят  где k – коэффициент теплопередачи через стенку рубашки, Вт/(м2∙град).Источником теплопотерь является охлаждение воды в рубашке: где Q2 – объем тепла, отдаваемого водой в рубашке, Дж; ср – удельная теплоемкость воды в рубашке, Дж/(кг·град); Мв – массовый расход воды на термостабилизацию, кг/с; Δt – изменение температуры воды от входа в рубашку до выхода из нее, К.Тепловой баланс в данной ситуации выражается равенством Q1 = Q2, которое приводит к уравнению относительно неизвестного х: .Решение данного уравнения: .Для оценки х необходимо определить входящие в формулу параметры. При средней температуре t = 45 oC (318 К) вода в рубашке имеет свойства, приведенные в табл. 2 . Таблица 2Table 2Свойства воды и воздуха при заданной средней температуреWater and air properties at a given average temperatureСвойстваЗначение свойств водыЗначение свойств воздуха при средней температуре 318 Кпри среднейтемпературе 288,4 Кпри среднейтемпературе 318 КПлотность ρ, кг/м3999,325990,11,1105Динамическая вязкость µ, Па∙с1 228,75 · 10–3601,35 · 10–319,35 · 10–6Кинематическая вязкость v, м2/с1,231 · 10–60,6075 · 10–617,455 · 10–6Теплопроводность λ, Вт/(м·град) 0,57550,64150,02795Удельная теплоемкость ср, Дж/(кг·град) 4 1894 177,51 005Скорость воды в рубашке W, м/с0,020,001– Поскольку задача не имеет однозначного решения, каким-то параметром необходимо задаться априори. Для этого удобно задаться скоростью потока теплоносителя. Примем скорость воды в рубашке W = 0,001 м/с. Тогда критерий Рейнольдса в рубашке составит где Dэкв – эквивалентный диаметр рубашки в цилиндрической части, м, равный  где S, П – площадь и смоченный периметр поперечного сечения потока жидкости.Вычислим Мв, кг/c:  . Число Прандтля Pr для воды в рубашке [12]: .Так как значение критерия Рейнольдса существенно меньше 2 000, то режим движения воды в рубашке ламинарный. Для ламинарного режима критерий Нуссельта Nu определяется по формуле [12]   Коэффициент теплоотдачи в рубашке от воды, Вт/(м2·К):Стенка рубашки выполнена из стали 20 толщиной  4  мм с внутренним диаметром 1,3 м и диаметром наружной поверхности 1,308 м. При средней температуре t = 45 oC (318 К) сталь 20 имеет теплопроводность λст = 84,08 Вт/(м·град) На наружной поверхности происходит охлаждение воздухом. Воздух при средней температуре t = 45 oC (318 К) имеет свойства, значения которых приведены в табл. 2. Высота реактора с двумя стандартными эллиптическими днищами, м, равна   где   – высота стандартного эллиптического днища, м. Для стенки такой высоты критерийГрасгофа равен [12]   , где параметр β = 1 / 293, К–1; g – ускорение свободного падения, м/с2.Охлаждение рубашки происходит путем свободной конвекции атмосферного воздуха вдоль вертикальной стенки, для которой в ламинарном режиме критерий Нуссельта равен [12]:Тогда коэффициент теплоотдачи от рубашки к атмосферному воздуху, Вт/(м2·К), равенгде λв – теплопроводность воды, Вт/(м·град).Суммарный коэффициент теплопередачи при теплопотерях от рубашки в окружающий атмосферный воздух,  Вт/(м2·К), равенВозвращаясь к определению неизвестного значения х, ºC, получим:Тогда   между рубашкой и ОС, ºC, составитМассовый расход воды для термостабилизации смеси Mв = 0,081 кг/с. Температура воды на входе в рубашку должна быть 46,81 oC, а на выходе 43,19 oC.Таким образом, применительно к тепловому расчету химического реактора для синтеза серосодержащего сорбента получен новый алгоритм, описывающий процесс термостабилизации реакционной смеси при подаче теплоносителя в рубашку с заранее заданной скоростью 0,001 м/си вычисленной входной и выходной температурами, такими, что теплопотери от рубашки в окружающую среду сбалансированы таким образом, что средняя температура в рубашке была равна стабилизируемой температуре реакционной смеси. Алгоритм теплового расчета стадии охлаждения реакционной смесиОхлаждению подвергается реакционная масса с температуры 45 до 20 °С.При условии постоянства потока теплоносителя (вода 10 °С = 283 К, нагревающаяся до 15 °С = 288 К) и неизменности его термодинамических свойств ΔТср определяется по формуле (1), где ΔT1 – разность начальных температур рабочей среды (45 °С = = 318 К) и охлаждающей воды (10 °С = 283 К); ΔT2 – разность конечных температур рабочей среды (20 °С == 293 К) и охлаждающей воды (15 °С = 288 К).Для режима охлаждения значение ΔТср, К, вычисляется по формуле   Критерий Рейнольдса в рубашке в режиме охлаждения реакционной смеси вычисляется по формуле .Так как критерий Рейнольдса менее 2 000, режим движения ламинарный. Значение критерия Прандтля Pr определяется по формуле [12]  Критерий Нуссельта для ламинарного режима [12]:  Коэффициент теплоотдачи в рубашке от воды, Вт/(м2·К), рассчитывается как отношение [12]Коэффициент теплопередачи k, Вт/(м2·К), определяется по формуле [12]  где αс – коэффициент теплоотдачи от реакционной смеси к стенке реактора; αТ – коэффициент теплоотдачи от стенки реактора к охлаждающей воде; коэффициент теплопроводности стали 20λм = 79,64 Вт/(м·К) при 353 К.Время охлаждения реакционной смеси τохл, ч, рассчитывается по формуле   , где Q – тепловая нагрузка, кДж.Массовый расход воды при охлаждении смеси, кг/ч, вычисляется по формуле   Затраты Z холодной воды, кг, при учете 5 % тепловых потерь составятТогда можно определить значения объемного расхода воды V в рубашке, м3/с, ее скорость W, мм/с:Здесь Sж = 0,174 м2 – площадь кольцевого сечения рубашки с диаметрами 1,212 и 1,3 м.Таким образом, применительно к тепловому расчету химического реактора для синтеза серосодержащего сорбента разработан алгоритм для расчета стадии охлаждения реакционной смеси, реализующий задачу теплопередачи через стенку реактора. Для оценки сопутствующих теплопотерь в окружающую среду на внешней границе рассматриваемой системы теплоносителей взята приемлемая величина в 5 % от общего теплового потока в процессе охлаждения реакционной смеси.   Время полного цикла работы реактораЦелесообразность  использования  химического реактора периодического действия обусловливается его КПД [13]:                  (2)где τр = 9 ч – продолжительность химических реакций [1]; τц – время полного цикла работы химического реактора; τв – время проведения вспомогательных операций:                                              (3)где τп.зп – время подготовки реактора, заполнения его суспензией и выгрузки продукта; τнагр = 0,88 ч – время нагрева суспензии [10]; τохл = 1,254 ч – время охлаждения суспензии.Примем КПД = 0,7. Тогда из (2) получаем: 0,7 = 9 / (9 – τв);τв = 3,857 ч.Тогда из формулы (3) следует:   Время полного цикла работы реактора, ч, составит:В данном разделе, исходя из взятого КПД = 0,7 для химического реактора периодического действия и вычисленных продолжительностей основных стадий синтеза серосодержащего сорбента, получена оценка времени полного цикла работы реактора на одну загрузку. Программа автоматизированного расчетаНа рис. 2, 3 представлен интерфейс программы, реализованной на языке C#, с результатами автоматизированного теплового расчета реактора по описанным алгоритмам.   Рис. 2. Интерфейс программы расчета реактора при термостабилизации реакционной смеси Fig. 2. Interface of the reactor calculation program at reaction mixture thermal stabilization    Рис. 3. Интерфейс программы расчета реактора при охлаждении реакционной смеси Fig. 3. Interface of the reactor calculation program at reaction mixture cooling   Работоспособность программы расчета реактора для синтеза серосодержащего сорбента подтверждается совпадением ее результатов расчета с результатами ручного расчета. ЗаключениеПредложены новые алгоритмы расчета теплообмена при термостабилизации и охлаждении реакционной массы, а также времени полного цикла работы реактора синтеза серосодержащего сорбента периодического действия. Создана компьютерная программа, выполняющая автоматизированный расчет реактора по заданным и вычисленным параметрам, а также расчет времени полного цикла работы реактора.