CC BY
52
УДК: 614.8
А. А. Салин, А. Д. Галеев, С. И. Поникаров ОЦЕНКА ПРИМЕНИМОСТИ ДОПУЩЕНИЯ ОБ ИДЕАЛЬНОМ КОНВЕКТИВНОМ ПЕРЕМЕШИВАНИИ РАСТВОРА ПРИ МОДЕЛИРОВАНИИ ИСПАРЕНИЯ
Ключевые слова: термогравитационная конвекция, растворы, испарение, аварийный пролив, математическая модель.
В статье представлена математическая модель испарения бинарного неидеального раствора с поверхности аварийного пролива, позволяющая определить характеристики процесса с учетом тепломассопереноса в жидкой фазе. В работе приводится сравнение результатов численных экспериментов, полученных: а) с учетом термогравитационной конвекции в слоях испаряющейся жидкости, б) с применением допущения об идеальном перемешивании во всем объеме раствора.
Keywords: thermogravitational convection, solutions, evaporation, emergency spill, mathematical model.
The article presents the mathematical model of non-ideal binary solution evaporation from the emergency spill surface, which allows determining the process characteristics with the heat and mass transfer in the liquid phase. This study compares the results of numerical experiments obtained by: a) subject to thermo-gravitational convection in the layers of the evaporating liquid, b) using the assumption of ideal mixing in the whole volume of the solution.
Введение
Неотъемлемой частью процедуры анализа риска на опасных производственных объектах нефтехимической отрасли промышленности является оценка концентрационных полей и размеров зон возможного токсического поражения при аварийных разливах химически- и пожароопасных веществ.
Для решения подобных задач необходимо применять комплексные математические модели, одновременно описывающие парообразование опасного вещества у границы раздела фаз и дальнейшее его распространение в атмосфере [1]. Данные тепло-массобменные процессы взаимосвязаны и значительно зависят как от метеорологических условий, так и от физических свойств жидкой фазы и испаряющегося компонента. В настоящее время рекомендуются к использованию модели, изложенные в разработанных и апробированных нормативных методиках [2,3]. При расчетах по методике [2] не принимается во внимание нестационарный характер испарения опасных веществ, обусловленный изменением состава жидкости. К тому же предлагаемая формула для расчета интенсивности испарения не учитывает влияние сил плавучести, проявляющихся в усилении или подавлении турбулентного перемешивания, что не может не сказываться на точности полученных результатов.
Нередко математическое описание процесса испарения сводится к использованию упрощенных аналитических зависимостей и эмпирических корреляций [4,5]. Использование СРБ-технологий позволяет учесть сложную взаимосвязь процессов поступления опасного вещества из пролива и его рассеяния в атмосфере благодаря численному решению трехмерных нестационарных уравнений гидродинамики и тепло-массообмена.
Ранее была разработана математическая модель испарения бинарных растворов [6], позволяющая определить характеристики процесса с учетом изменения состава и температуры жидкой фазы и трансформации турбулентной структуры потока
над поверхностью пролива. Модель основана на решении дифференциальных уравнений
сохранения переноса импульса, массы, энергии, замыкаемых RNG k-e моделью турбулентности и реализована с использованием программного комплекса FLUENT (Certificate Number - ANS2012-S021; Customer # 657938). Применение в указанной модели допущения о том, что жидкость идеально перемешана по высоте слоя, позволяет значительно упростить математическое описание процесса и экономить машинное время.
Целью данной работы являлось исследование влияния термогравитационной конвекции в слое жидкости на её перемешивание и, как следствие, на динамику интенсивности испарения раствора.
Упрощенная математическая модель испарения
При испарении раствора у границы раздела фаз снижается концентрация легколетучих компонентов и происходит охлаждение поверхностного слоя, что приводит к образованию плотностных и температурных градиентов в объеме жидкости. Предположение о возникновении движущей силы, инициирующей объемное перемешивание раствора, которое послужило основой для применения допущения о равномерном распределении компонентов раствора в объеме пролива. Для подтверждения правомерности применения этого допущения решалась упрощенная двумерная задача испарения раствора аммиака на базе программного комплекса FLUENT.
Массовый поток компонентов с поверхности пролива определялся из уравнения (1), с учетом поправки на стефановский поток [7]. Данное уравнение справедливо в случае отсутствия возмущающего действия испаряющегося компонента на структуру потока над проливом.
Vi^ Ci\Ks
R0 'TS
(1)
где МI - молекулярная масса /-ого компонента; -универсальная газовая постоянная; Кэ - поправка на
Стефановский поток; Р^Т^О ) - парциальное давление /-ого компонента в растворе, зависящее от температуры Т и концентрации компонента О,.
Коэффициент массоотдачи к рассчитывается по формуле [5]:
кг-0,029.и0'78.0-°'11.8сГ0>67 (2)
где и - скорость ветра [м/час]; й - диаметр пролива, или другой характерный геометрический размер [м]; Бй/ - число Шмидта, которое определяется по следующему уравнению:
" (3)
Р-О;
где у - вязкость газа над проливом; р - плотность газовой фазы над проливом; й - коэффициент молекулярной диффузии компонента в газовой фазе:
Экспериментальные значения парциального давления -ого компонента водного раствора Р приведены в работе [8] и в модели определяются интерполяцией в зависимости от температуры Т и содержания аммиака в жидкости О..
Температура тонкого приповерхностного слоя жидкости ТБ определяется из уравнения теплового баланса с учетом испарительного охлаждения и естественного нагрева пролива за счет солнечной радиации. Теплообмен с атмосферным воздухом не учитывался, ввиду малых значений соответствующего теплового потока [9].
Поле температур в слое жидкой фазы получено путем решения системы двухмерных нестационарных уравнений Навье-Стокса, неразрывности и сохранения энергии в приближении Буссинеска [10]:
ди ■ дх ■
= 0
ди: -+-
д дXj
дд
Ро
дх-
д
дt
М
дх ■
д
С я ^ ди-
V-
д
V и
\
7
дТ
дхП
(4)
(5)
(6)
Модель Буссинеска рассматривает плотность как постоянное значение во всех решаемых уравнениях, за исключением составляющей плавучести в уравнении моментов:
(р-р0)д* -р0^-Т0)д (7)
Это уравнение получено с использованием аппроксимации Буссинеска р=ро(1-рлТ), для того, чтобы исключить р из составляющей плавучести.
Отток теплоты при испарении равен:
дИСП=£ г-ил .дн,.) (8)
где АН, - удельная теплота парообразования /'-ого компонента.
При определении теплового потока от подстилающей поверхности для твердого слоя решалось уравнение теплопроводности:
дТ = српРп ~^ = Ап
сд2т
д
Я 2
ду ;
На границе, разделяющей твердое тело и жидкость (у=0), выполняется условие равенства тепловых потоков:
дТ
' ду
= -Ап —
дТ
д
(10)
у=-0
у=+0
При моделировании использовалось предположение, что теплофизические свойства жидкости и подстилающего слоя незначительно зависят от температуры.
Результаты расчетов
Один ряд численных экспериментов проводился с учетом термогравитационной конвекции в слоях испаряющейся жидкости, другой - с применением допущения об идеальном перемешивании во всем объеме раствора, при прочих равных условиях: начальная температура жидкости и подстилающего слоя принималась равной 303 °К, высота слоя жидкости составляла 0,05 м, длина расчетной области 1 м, и глубина грунта 2 м. Каждая граница расчетной области является плоскостью симметрии. Скорость ветра и в уравнении (2) варьировалась и принимала значения 1; 2,5 и 5 м/с. В качестве исследуемой жидкости был выбран 25% (масс.) водный раствор аммиака.
Рис. 1 - Масса испарившегося аммиака при скорости ветра 1 м/с: расчет с учетом термоконвективного перемешивания; .......расчет с
учетом допущения об идеальном перемешивании
Из рисунка 1 видно, что математическая модель с допущением об идеальном перемешивании жидкой фазы по высоте слоя дает на 5% завышенные значения массы испарившегося аммиака, по сравнению с расчетами в которых учитывается термогравитационная конвекция в слое жидкости. В последнем случае испаряющийся компонент переходит в газовую фазу лишь с приповерхностного слоя, в котором на начальном этапе резко снижается концентрация аммиака и температура.
Из рисунка 2 можно сделать вывод, что температура приповерхностной пленки ниже температуры основного объема раствора. Охлаждение приповерхностного слоя жидкости вследствие испа-
800 1000 ^ сек
Рис. 3 - Масса испарившегося аммиака при скорости ветра 2,5 м/с: расчет с учетом термоконвективного перемешивания; .......расчет с учетом допущения об идеальном перемешивании
304 302 300 298 296 294 292 290 288 286
800 1000 ^ сек
1400 1600
, 2,5 Е
0,5
400
1600
284
400
1200
рения приводит к увеличению его плотности. В результате возникают условия для развития конвективных движений, увлекающих тяжелые массы жидкости в глубину и осуществляющих транспорт более легких теплых масс к поверхности.
^ сек
Рис. 2 - Средняя температура раствора в зависимости от времени при 1 м/с: температура
приповерхностного слоя жидкости, 0 температура основного объема - при расчете с учетом
термоконвективного перемешивания; .......
расчет с учетом допущения об идеальном перемешивании
Благодаря температурным градиентам происходит интенсивное перераспределение массы и энергии в рассматриваемой области. Конвективные потоки способствуют выравниванию температуры и концентрации по высоте слоя. Отличие температур на поверхности и в глубине слоя для аммиачной воды в основном не превышало 1,5° при скорости ветра 1 м/с за исключением резкого снижения температуры поверхности на 4 градуса в первые 20 секунд.
Полученные результаты свидетельствуют о том, что испаряющуюся жидкость можно считать перемешанной по высоте слоя. Это позволяет использовать для определения парциального давления компонентов раствора их средние по объему концентрацию и температуру. В случае «идеального перемешивания» раствора средняя температура жидкости в начале эксперимента равна температуре основного объема, а через 15 минут приближается к температуре поверхности.
При скорости 2,5 м/с отличие в значениях массы испарившегося аммиака в случаях идеального перемешивания и термогравитационной конвекции также составляет около 5% по истечении 30 минут, что видно из рисунка 3. При данной скорости в обоих случаях за рассматриваемые промежуток времени аммиака испарилось примерно на 30% больше чем при 1 м/с.
В целом характер охлаждения раствора при 2,5 м/с (рисунок 4) подобен результатам при 1 м/с. После 30 минут численного эксперимента температура раствора в среднем на 5° ниже по сравнению с расчетами при 1 м/с, а первоначальный скачок температуры приповерхностного слоя составил 7°, вместо 4°.
Рис. 4 - Средняя температура раствора в зависимости от времени при 2,5 м/с: температура приповерхностного слоя жидкости, в температура основного объема - при расчете с учетом термоконвективного перемешивания;
.......расчет с учетом допущения об идеальном
перемешивании
При скорости ветра 5 м/с расчетные значения массы испарившегося аммиака, полученные с применением допущения об идеальном перемешивании и с учетом термогравитационной конвекции, различаются на 7%.
Из полученных результатов можно сделать вывод, что разница в значениях массы испарившегося вещества, полученная с использованием разных подходов возрастает с увеличением скорости ветра. Это связано с тем, что с увеличением скорости ветра лимитирующее влияние тепломассопереноса в слое жидкости усиливается, так как подвод компонента из глубинных слоёв жидкости происходит более медленно, чем его отвод от поверхности пролива в газовую фазу.
Температура идеально перемешанного раствора при скорости ветра 5 м/с в конце численного эксперимента не только приблизилась к температуре поверхности раствора (в расчетах без допущения), но и опустилась на 1° ниже как видно из рисунка 6, что обусловлено более высокой интенсивностью испарения в случае идеального перемешивания.
0 200 400 600 800 1000 1200 1400 1600 1800 t, сек
Рис. 5 - Масса испарившегося аммиака при скорости ветра 5 м/с: расчет с учетом термоконвективного перемешивания; .......расчет с
учетом допущения об идеальном перемешивании
При проведении анализа предполагалось, что режим течения в слое жидкости является ламинарным. Дополнительное перемешивание раствора за счет волнообразования будет приближать истинную температуру слоя к средней температуре раствора в случае с идеальным перемешиванием.
0 200 400 600 800 1000 1200 1400 1600 1800 t, се к
Рис. 6 - Средняя температура раствора в зависимости от времени при 5 м/с: температура
в темпе-
ратуры по толщине слоя испаряющейся жидкости. Таким образом, испаряющуюся жидкость можно считать перемешанной по высоте слоя. Это позволяет использовать для определения парциального давления компонентов раствора осредненные по высоте слоя жидкости значения концентрации и температуры.
Кроме того, допущение об идеальном перемешивании в слое пролива позволяет значительно упростить математическое описание процесса испарения, а несколько завышенные результаты интенсивности испарения позволяют увеличить запас надежности при прогнозировании концентрационных полей токсичных и взрывопожароопасных веществ в случае аварийной разгерметизации оборудования.
Литература
3
приповерхностного слоя жидкости, ратура основного объема - при расчете с учетом
термоконвективного перемешивания; .......
расчет с учетом допущения об идеальном перемешивании
Выводы
В данной работе была проведена оценка влияния термогравитационной конвекции в слое испаряющегося неидеального раствора на характеристики испарения.
Полученные результаты численного моделирования показали, что действие термогравитационной конвекции приводит к выравниванию темпе© А. А. Салин - асп.т каф. машин и аппаратов химических производств КНИТУ, [email protected]; А. Д. Галеев - канд. техн. наук, доц. той же кафедры, [email protected]; С. И. Поникаров - д-р техн. наук, проф., зав. каф. машин и аппаратов химических производств КНИТУ, [email protected].
Старовойтова Е.В. Моделирование и оценка последствий аварийного выброса сжиженного аммиака / Е. В. Старовойтова, А. Д. Галеев, С. И. Поникаров // Вестник Казан. технол. ун-та. - 2011. - № 13. - С.175-179. РД-03-26-2007. Методические указания по оценке последствий аварийных выбросов опасных веществ (утв. приказом Федеральной службы по экологическому, технологическому и атомному надзору от 14 декабря 2007 г. №859).
TNO Yellow Book. Methods for the calculation of physical effects due to releases of hazardous materials (liquids and gases) // Committee for the Prevention of Disasters, CRP 14E, The Hague, The Netherlands, - 1997. - P. 870.
4. Mackay D. Evaporation rates of liquid hydrocarbon spills on land and water / D. Mackay, R. S. Matsugu // Canadian Journal of Chemical Engineering. - 1973. - Vol. 51. - № 4. - PP. 434-439.
5. Kawamura P. I., Mackay D. The evaporation of volatile liquids / P. I. Kawamura, D. Mackay // Journal of Hazardous Materials. - 1987. - Vol. 15. - Issue 3. - PP. 343-364.
6. Салин А.А. Сравнение результатов численного моделирования водного раствора аммиака с экспериментальными данными / А. А. Салин, А. Д. Галеев, С. И. Пони-каров // Вестник Казан. технол. ун-та. - 2012. - № 15. -С. 225-228
7. Law, C.K. Combustion physics / C.K. Law - Cambridge University Press, 2006.
8. Wilson, T.A. The total and partial vapor pressures of aqueous ammonia solutions / T.A. Wilson - Engineering Experiment Station, 1925. - Bulletin №146.
9. Старовойтова Е.В. Численный анализ процесса парообразования при кипении аварийного пролива сжиженного газа / Е. В. Старовойтова, А. Д. Галеев, С. И. Пони-каров // Пожаровзрывобезопасность. - 2011. - № 2. -С.24-28.
10. Галеев А.Д. Численный анализ процесса испарения однокомпонентной жидкости с поверхности аварийного пролива в условиях неизотермичности / А.Д. Галеев, С.И. Поникаров, М.А. Долгова // Химическая промышленность сегодня. - 2012. - № 2. - С. 43-48.
5
4
, 3 Е
2
0
305
300
295
290
285
280
275
