Physico-mathematical modeling of fluoride hydrogen desublimation from gas mixtures onto walls of a condenser.pdf Для очистки потоков гексафторида урана (ГФУ) от легких примесей используются очистительные каскады, состоящие из блоков диффузионных машин. Отфильтрованные очистительными каскадами вещества, удаляемые из технологической цепочки, представляют собой смесь газов, основными компонентами которых является ГФУ, фтористый водород (HF), воздух и другие неконденсируемые примеси в незначительном количестве [1]. Газовая смесь подвергается фракционной разгонке для возврата ГФУ в технологию и выделению фтористого водорода, для дальнейшей переработки. Технологическая линия фракционной разгонки газовых смесей состоит из двух последовательно расположенных установок, работающих по принципу физических процессов конденсации и десублимации на охлажденных поверхностях. Установка улавливания ГФУ обеспечивает десублимацию вещества в промежуточных емкостях (ПЕ), охлаждаемых холодным воздухом с температурой 183 К. Де-сублимированный в ПЕ ГФУ возвращается в технологическую цепочку. Установка очистки от HF и прочих примесей предназначена для улавливания проскоков ГФУ после установки улавливания ГФУ и конденсации фтористого водорода в емкостях-осадителях объемом 24 л, охлаждаемых жидким азотом с температурой 77 К [2, 3]. Фтористый водород, сконденсированный в емкостях-осадителях установки, направляется на переработку, а компоненты воздуха и неконденсируемые примеси откачиваются вакуумными насосами, фильтруются специальными установками и выводятся в атмосферу. Для десублимации фтористого водорода из газовой смеси в производстве ОАО «Сибирского химического комбината» используются емкости-десублиматоры, охлаждаемые жидким азотом до температуры 77 К. Экономически более целесообразно охлаждать десублиматоры холодным воздухом. В настоящей статье проводится математическое моделирование процесса десублимации фтористого водорода в емкостях-десублиматорах при температурах 77 и 113 К. Десублиматор представляет собой цилиндрическую емкость с 6 пластинами-сегментами. Пластины-сегменты расположены поперек проходного сечения потоку газовой смеси, и повернуты под углами 60° по отношению друг к другу. Стенки десублиматора охлаждаются либо жидким азотом при температуре 77 К, либо охлажденным до температуры 113 К воздухом. В верхнюю часть десублиматора подается газовая смесь фтористого водорода и воздуха. Математическая модель строится при следующих допущениях: течение в де-сублиматоре предполагается квазиодномерным; периметр потока полагается равным площади охлаждающей поверхности, деленной на высоту цилиндра десуб-лиматора; температура внутренних поверхностей десублиматора равна температуре охладителя; десублимация фтористого водорода происходит только на охлаждаемых стенках. Десублимация фтористого водорода при наличии неконденсирующегося воздуха сопровождается диффузией HF через пограничный слой к холодной поверхности, на которой HF десублимируется. Принимается, что все сопротивление переносу массы и теплоты сосредоточено в этом слое. Диффузия через слой описывается уравнением [4]: (1) где nHF - мольный поток HF; D - коэффициент диффузии HF по отношению к неконденсирующемуся воздуху; n - общая молярная концентрация; yHF - мольная доля HF; s - расстояние от поверхности десублимации. Для расчетов коэффициента диффузии в воздухе используется формула [4] (2) p (HF+U/3)' где uHF, ub - молярные диффузионные объемы HF и воздуха; |xHF, - молярные массы этих веществ, кг/кмоль; p - давление в Па; T - температура в К. Диффузионный объем воздуха равен 201, диффузионный объем HF вычисляется через объемы составляющих его атомов и равен 7.46. При вычислениях удобно использовать произведение Dn , входящее в формулу (1). Заменяя в (2) давление из уравнения состояния идеального газа p = nRT , запишем (3) r (HF+U/3)' Если принять в формуле (1) в качестве произведения Dn среднее значение Dn поперек пограничного слоя, то можно получить интеграл уравнения (1): nHF = ^ ln f ^ 1, (4) sf К 1 - УН ) где yHF0, yHF - мольные концентрации HF на границе раздела фаз и в потоке соответственно; Sf - толщина пограничного слоя. Для записи квазиодномерного уравнения сохранения числа молей nHF в де-сублиматоре введена координата z . Координата z направлена параллельно оси десублиматора от точки входа первоначального потока к его выходу из устройства. Тогда для числа молей HF можно записать: dSun ж = ППш, (5) dz где S - площадь поперечного сечения десублиматора; u - скорость потока; П -периметр. Подставляя в (5) величину потока из (4) и выражая относительные концентрации через числа молей HF и воздуха в единице объема, получим dSunHF = п Dn dz sf ln | 1 + | - ln | 1 + nb Левую и правую части этого уравнения поделим на известную постоянную величину расхода воздуха (Sunb) = Gb и введем функцию z) = nHF/nb . В результате получим уравнение # = - П -D- [ln (1 +1) - ln(1 + ^)], (6) dz SfGb которое не содержит неизвестную u , скорость потока. Коэффициент Dn в этом уравнении зависит от температуры. Так как температура меняется поперек пограничного слоя при интегрировании (6), то вычислялась средняя величина Dn по формуле т °,75 + г,,",75Y 1 1л1/2 nD = 0,0003204^- HF Ц b п{ 1/3 1/3\2 R (UHF +Ub ) в которой т - температура газа в потоке, т0 - его температура на поверхности охлаждения. Для вычисления температуры использовалось уравнение тепловой энергии dS (CP HFnHF + CPbnb)uT n . „(„ „) -1-d---= ПСР HFnHF - аП (T - T0 ), (7) dz где cP HF - молярная теплоемкость HF; cPb - молярная теплоемкость воздуха; a - коэффициент теплоотдачи. В правой части (7) первое слагаемое описывает унос из потока тепла в результате диффузии, второе - в результате теплоотдачи. С помощью (5) уравнение (7) приводится к виду P HFnHF + CPbnb )) = ~аП (T - T0 ) . Su ( (8) Поделив левую и правую части (8) на расход воздуха Gb , получим )CP HF + CPbnb )) =-аП (T - T0 ). dz Gb Для определения коэффициента теплоотдачи выразим его через число Нус-сельта Nu . Так как число Рейнольдса течения мало, то согласно [4] Nu = 4 . Система уравнений (6) и (9) решается численно. Необходимое для расчетов значение nHF0 на охлаждающей поверхности вычисляется при давлении насыщенного пара на этой поверхности. Граничное значение функции ^ (0) рассчитывалось из заданного значения отношения расходов 4(0) = (nHF/nb) = Ghf/Gb . Температура при z = 0 полагалась равной 298 К. Расчеты проводились для случая охлаждения жидким азотом с температурой равной T0 = 77 K и холодным воздухом с температурой T0 = 113 K . На рис. 1, а показана зависимость десятичного логарифма отношения концентраций lg (nHF/nb) от расстояния z . Видно, что на участке z < 0,18 м концентрация фтористого водорода вследствие вымораживания быстро уменьшается до величины nHF0 , соответствующей давлению насыщенного пара HF при температуре охлаждающего жидкого азота. Согласно расчетам, эта концентрация равна 1,5 -10-10 кмоль/м3. На рис. 1, б для этого случая показана зависимость температуры от z . Видно, что на участке z < 0,18 м температура газа быстро уменьшается до температуры азота T0 = 77 K . (9) Т, К 250 200 150 100 50 2 0 0 ~1-1-I- 0,1 0,2 - 0,1 -I-1-I-1-I-1-I 0,2 0,3 0,4 z, м 0,3 0,4 z, м Рис. 1. Зависимость десятичного логарифма отношения концентраций фтористого водорода и воздуха (а) и температуры газовой смеси (б) от высоты емкости. T0 = 77 К Аналогичные расчеты были проведены для случая охлаждения емкости при температуре T0 = 113 K . На рис. 2, а приведена зависимость lg (nHF/nb) для этого случая. Из рисунка следует, что при охлаждении воздухом вымораживание газообразного HF происходит быстрее, чем вымораживание с помощью азота. Этот факт объясняется увеличением коэффициента диффузии с ростом температуры, которая при охлаждении воздухом выше, чем при охлаждении азотом. Остаточная концентрация паров HF при охлаждении воздухом также выше, чем при охлаждении азотом. Ее величина согласно расчетам составляет 1,04 -10-6 кмоль/м3. 11 0 -1- Т, К -| 280240200 160 -12080 к ад -2 -3-4 "Г "Г "Г "г "Г 0,1 0,2 0,3 0,4 z, м 0 0,1 0,2 0,3 0,4 z, м 0 Рис. 2. Зависимость десятичного логарифма отношения концентраций фтористого водорода и воздуха (а) и температуры газовой смеси (б) от высоты емкости. T0 = 113 К Полученные расчетные концентрации HF не учитывают возможной конденсации перенасыщенного пара HF внутри пограничного слоя в виде капелек тумана, диффузия которых на охлаждающую поверхность будет затруднена вследствие их больших размеров по сравнению с размером молекул. Поэтому полученные результаты по концентрации можно рассматривать как нижний предел вероятных значений nHF. В [5] показано, что при низких температурах вблизи давления насыщенных паров молекулы фтористого водорода образуют ассоциации, содержащие до 6 молекул HF в одном комплексе. Эти ассоциации обладают большим молекулярным весом и в соответствии с теорией большими диффузионными объемами. Это приводит, как следует из формулы (2), к уменьшению коэффициента диффузии и, следовательно, к уменьшению диффузионного потока на охлаждаемые стенки де-сублиматора. На рис. 3, а показано распределение десятичного логарифма отношения концентраций фтористого водорода и воздуха по высоте емкости при охлаждении жидким азотом с учетом наличия комплекса HF, содержащего 6 молекул фтористого водорода. Видно, что вследствие уменьшения потока на стенки концентрация в объеме уменьшается значительно медленнее, чем для неассоциированного HF. На рис. 3, б показано распределение температуры для этого случая. На рис. 4, а показано распределение десятичного логарифма отношения концентраций фтористого водорода и воздуха по высоте емкости при охлаждении воздухом с учетом наличия комплекса HF, содержащего 6 молекул фтористого водорода. На рис. 4, б показано распределение температуры для этого случая. Т, К-250200 15010050 1-1-I-1-I-1-I-1-I-1-I 0 0,1 0,2 0,3 0,4 z, м 2-, 0- J -2- -4- 0 т-1-I-1-I-1-г 0,1 0,2 0,3 0,4 z, м "ад Рис. 3. Зависимость десятичного логарифма отношения концентраций фтористого водорода и воздуха (а) и температуры газовой смеси (б) от высоты емкости. T0 = 77 К Т, К 280240200160120 80 -|-I-|-I-|-I-|-I-| 0 0,1 0,2 0,3 0,4 z, м 0 4 -I-1-I-1-I-1-I-1-I 0,1 0,2 0,3 0,4 z, м Рис. 4. Зависимость десятичного логарифма отношения концентраций фтористого водорода и воздуха (а) и температуры газовой смеси (б) от высоты емкости. T0 = 113 К С увеличением степени ассоциации скорость очистки газовой смеси от HF существенно уменьшается. Тепловое равновесие между температурой стенок емкости и температурой потока смеси газов устанавливается на одной трети пути газа по высоте емкости. Концентрация фтористого водорода стремится к равновесному значению, соответствующему давлению насыщенных паров при температуре стенок емкости.

Скачать электронную версию публикации