11УДК 536.423.1
ИСПАРЕНИЕ ХИМИЧЕСКИ ОПАСНЫХ ВЕЩЕСТВ С НИЗКОЙ ТЕМПЕРАТУРОЙ КИПЕНИЯ
Попов В.М.*, к.ф.-м. н., Ильюшонок А.В.**, к.ф.-м.н., доцент, Катков В.Л.***, д.т.н., профессор, Лешенюк Н.С.**, д.ф.-м.н., профессор *Институт физики им. Б.И. Степанова НАН Беларуси, **Командно-инженерный институт МЧС Республики Беларусь ***Объединенный институт проблем информатики НАН Беларуси
e-mail: [email protected]
Рассмотрено испарение жидкости, имеющей низкую температуру кипения и пролитой на поверхность в результате аварии. Описаны три режима испарения, реализуемые в перегретой, кипящей и не кипящей жидкости. В уравнение баланса энергии, определяющего процесс испарения, включены перенос тепла с подстилающей поверхности и из воздуха, поглощение жидкостью коротковолнового излучения Солнца и длинноволнового излучения атмосферы, а также испускание теплового излучения поверхностью жидкости. Приведены примеры расчета испарения жидкости при различных атмосферных условиях.
Evaporation of liquid having a low boiling point is considered after its spilling onto a surface as a result of an accident. Three modes of evaporation condition are described, which are realized in superheated, boiling and not boiling liquid. In the energy balance equation, which determines the evaporation process, heat transfer from underlying surface and from air, the absorption of shortwave solar radiation and atmospheric longwave radiation by liquid, as well as the emission of thermal radiation from a surface of liquid, is included. Examples of calculation of evaporation of the liquid are given under various atmospheric conditions.
(Поступила в редакцию 2 сентября 2010 г.) Введение
Современной промышленностью производится, транспортируется и используется большое количество аварийно химически опасных веществ (АХОВ), имеющих при нормальных условиях низкую температуру кипения Tb, например, хлор и аммиак (Tb = -34,03 °C и Tb = -33,43 °C соответственно). Чаще всего эти вещества хранятся в закрытых резервуарах в сжиженном виде при низкой температуре или при высоком давлении. В результате аварии, приводящей к разгерметизации резервуара, жидкость будет пролита на поверхность земли. Можно считать, что выброс жидкости происходит мгновенно. На поверхности земли, имеющей температуру выше температуры кипения АХОВ, пролитое вещество закипает и испаряется. Температура и скорость испарения пролитой жидкости обычно служат граничными условиями в задаче о рассеянии паров АХОВ в атмосфере и определении границы зоны химического заражения. Кроме подвода тепла от подстилающей поверхности на кипение и испарение оказывают влияние перенос тепла из атмосферы, поглощение излучения Солнца и теплового излучения атмосферы, испускание излучения поверхностью жидкости, а также унос ветром пара от открытой поверхности жидкости. Перенос излучения и тепла из атмосферы зависит, в свою очередь, от ее состояния, которое определяется температурой воздуха, скоростью ветра, временем года и суток, широтой и долготой местности, величиной облачности и рядом других параметров. Целью данной работы является комплексный учет всех процессов, определяющих температуру и скорость испарения пролитой жидкости.
1. Уравнение баланса энергии в объеме жидкости
Будем считать, что АХОВ хранилось в резервуаре в сжиженном состоянии при высоком давлении и температуре, равной температуре окружающей среды. В результате ава-
рии оно было пролито на поверхность земли, ограниченную обваловкой. Площадь свободной поверхности жидкости определяется площадью обваловки и считается известной. Полагаем, что пролитая жидкость не поглощается подстилающей поверхностью и ее температура кипения ниже температуры окружающей среды. В этом случае начальная стадия интенсивного испарения, возникающая сразу после пролива жидкости, протекает за счет уменьшения внутренней энергии жидкости при температуре выше температуры кипения. Эта стадия будет отсутствовать в тех случаях, когда вещество хранилось в сжиженном состоянии при низкой температуре и атмосферном давлении. В процессе кипения при постоянной температуре скорость испарения полностью определяется скоростью подвода тепла от подстилающей поверхности, переносом тепла из атмосферы, величинами потоков поглощенной солнечной радиации, теплового излучения атмосферы и испущенного поверхностью жидкости излучения. После прекращения кипения на тепловой баланс жидкости, кроме перечисленных источников тепла, оказывает влияние перенос ветром пара от ее поверхности.
В общем виде закон сохранения энергии для объема испаряющейся жидкости, имеющего глубину И и площадь свободной поверхности £а, можно представить следующим образом:
*(Р ь^СьТ ) = ОА + ОА + 6А - 4ба (ь+срт), (1)
где I - время, прошедшее с момента пролива;
Т - температура пролитой жидкости;
р ь - плотность жидкости;
Сь - удельная теплоемкость жидкости;
Ср - удельная теплоемкость пара при постоянном давлении;
- площадь поверхности жидкости находящейся в контакте с землей;
О, Оа - тепловые потоки, подводимые к единице площади жидкости соответственно
от подстилающей поверхности и окружающего воздуха;
Ог - суммарный поток солнечной радиации, теплового излучения атмосферы и излучения поверхности жидкости;
ь - удельная теплота парообразования;
4 - масса жидкости, испаряющейся в единицу времени с единицы площади.
Из (1) видно, что изменение со временем внутренней энергии жидкости равно разности потоков энергии, подводимой к объему жидкости и теряемой ею при испарении и испускании теплового излучения.
2. Испарение перегретой жидкости
Сразу после разгерметизации резервуара, когда происходит быстрый выброс жидкости в обваловку, ее испарение осуществляется за счет собственного тепла. На данной стадии всеми потоками тепла в (1), подводимыми к жидкости в течение короткого времени, можно пренебречь. Тогда для полной массы пролитой жидкости вместо (1) получим следующее уравнение:
* (шЬСЬТ) _ йш
ь (Ь + СрТ), (2)
где шь = рЬБаИ - масса всего объема пролитой жидкости.
При условии, что Ь >> СрТ можно считать Ср ~ Сь. Тогда, вычисляя производную в (2), получим
т,
ё (СЬТ) _ ёть
Ж
Ж
Ь .
Интегрируя последнее выражение в пределах от температуры окружающего воздуха Та до температуры кипения Ть, найдем, что масса образовавшегося пара та _ тЬ0 - тЬ определяется следующим выражением [1]:
то _ ть,о I1 - ехР
С (Та - Ть)
Ь
(3)
где тЬ0 - вся масса пролитой жидкости.
Одновременно с паром при кипении выбрасывается аэрозоль. Его масса полагается равной массе пара то или массе оставшейся жидкости ть, когда та > ть [1].
Для большинства задач можно считать, что образовавшееся облако однородно по температуре, а ее величина Ту определяется следующим выражением:
1 та
Т _- [ Тёта
та 0
т-Ь-
а С
ехр
ь у
- ^ (Та - Ть )
Т -Ь-
Ть С
Ь
1 - ехр
-Ь(Та - Ть)
Величину Ту можно использовать в качестве начального значения температуры облака в задаче о рассеянии залпового выброса пара.
3. Потоки тепла к объему жидкости
Тепловой поток от подстилающей поверхности определяется из решения уравнения теплопроводности для полубесконечной среды. Оно имеет вид [2]:
2 _0
а (Т_„- Т)
^а п I
(4)
где ось г направлена вертикально вверх;
Т-х - температура почвы на глубине 2 = -ю;
А - коэффициент теплопроводности почвы; а _ А/(р) - коэффициент температуропроводности;
р - плотность почвы;
С - удельная теплоемкость почвы.
Формула (4) используется в предположении, что распределение температуры в почве по глубине 2 устанавливается быстрее изменения температуры жидкости.
Поток тепла из окружающего воздуха на поверхность пролитой жидкости согласно закону Ньютона-Рихмана пропорционален разности температур жидкости и атмосферы и определяется следующим уравнением [3]:
й (<)_ ка (Та - Т) ,
где коэффициент теплоотдачи ка определяется по эмпирическим формулам, полученным путем обработки результатов измерений методом анализа размерностей.
Обычно он выражается в виде соотношения между числами Нуссельта Ып, Рейнольд-са Яе и Прандтля Рга:
Ып = 0,037 х Рг13 х Яе4/5;
здесь для турбулентного потока Яе>5х105 и Ып = ка • 2гр/Ха , Яе = рап^10-(2гр))па, Рга = ^/А « 0,786,
где гр - радиус свободной поверхности жидкости; Ха - коэффициент теплопроводности воздуха;
ра - плотность воздуха;
пк10 - скорость ветра на высоте 10 м;
Па - динамическая вязкость воздуха;
Рга - число Прандтля для воздуха;
Уа - кинематическая вязкость воздуха;
Ба - коэффициент диффузии воздуха.
Радиационный поток состоит из потоков коротковолнового солнечного излучения
, длинноволнового излучения атмосферы1 и собственного излучения поверхности
жидкости ° [4]:
О =(1 - а) + Ош - Оь.
Здесь а - альбедо поверхности жидкости. Значение а вычислялось по формуле для водной поверхности:
0,05
а = ■
sin hQ + 0,15
где hQ - угловая высота Солнца над горизонтом.
Поток коротковолнового излучения Солнца на поверхности Земли можно определить с помощью формулы Кондратьева [5]:
QW - т^Т i -
1 + ет/ sin hQ
где Q0 - солнечная постоянная;
е - эмпирическая функция высоты Солнца, заданная таблично;
т - оптическая толщина атмосферы для интегрального потока, принимаемая равной 0,105; f = 0,5607 - эмпирический коэффициент; n - облачность в долях единицы.
Угловая высота Солнца определяется следующим выражением:
sin hQ - sin ф • cos 5 + cos ф • cos 5 • cos 0 ,
где ф - широта местности;
1 В атмосферной физике данное излучение называется встречным (back radiation). 74
Вестник Командно-инженерного института МЧС Республики Беларусь, №1 (13), 2011 5 - склонение Солнца;
0 _ (л/12)(ё -12) - часовой угол, отсчитываемый от полудня;
¡ё - суточное время, ч.
Длинноволновое излучение атмосферы является функцией температуры и влажности на высоте 2 м, а также облачности. Эта функция определяется соотношениями, полученными Б.Б. Ибо [6]:
_ £а°Т24(1 + СЬП1);
^а _ СеГ еХР
Схъ _ 0,22;
( 350 ^
Т
V 2 у
(0,15, Т2 < 273,15; С _ [0,14, Т2 > 273,15,
где Т2 - температура воздуха на уровне 2 м;
е2 - парциальное давление водяного пара на уровне 2 м, гПа;
а - постоянная Стефана-Больцмана.
Собственное излучение поверхности жидкости задается известной формулой Стефа-на-Больцмана.
Скорость испарения не кипящей жидкости д' равняется произведению коэффициента переноса массы кт и плотности пара р над ее поверхностью:
Ч'_ктр _ к,(5)
К1
где ц - молярная масса пара;
ру (Т) - давление насыщенного пара, которое является известной функцией температуры жидкости Т.
Выражение для коэффициента переноса массы кт было получено на основе экспериментальных исследований. Оно имеет вид [7]:
кт _ Ст&т • <78 (2Гр)-0,11 • ¿с-0,67, (6)
где Ст&т = 0,004786 м0,33/с0,22; 8с _Уу^а - число Шмидта;
Уу - коэффициент кинематической вязкости пара;
Ва - коэффициент диффузии пара в воздухе.
Отметим, что для газов и паров Бс ~ 0,8. Аналогичное соотношение для коэффициента переноса массы было приведено в работе [8].
4. Испарение кипящей жидкости
На стадии кипения температура жидкости остается постоянной, а подвод тепла от внешних источников компенсирует потери энергии на испарение и излучение с поверхности.
Полагая равной нулю производную от температуры по времени в уравнении теплового баланса (1), для скорости испарения жидкости ч' с единицы площади поверхности в момент времени I получим следующее выражение:
О, (')^ + Оа + О, Ч'(* ) =-\-• (7)
В уравнении (7) зависимость от времени Ч (¿) полностью определяется тепловым потоком от подстилающей поверхности к объему жидкости.
Время, в течение которого происходит кипение жидкости %, находится из равенства скоростей испарения кипящей (7) и не кипящей (5) жидкостей на границе кипения. Из этого равенства с учетом (2) получаем:
Ь (Т-.- Т)
Ру (Ть Оа + О,
(8)
к
ять
ь
Масса испарившейся за это время жидкости равна интегралу по времени в пределах от 0 до % от потока пара с единицы поверхности (7), умноженному на площадь свободной поверхности:
ж=^Ч^Ь + о±а,,А. (9)
Ь \ап Ь
Как показывают расчеты, время кипения жидкости 1ь (например, для хлора и аммиака) невелико. Но за это время поток тепла (4) от подстилающей поверхности уменьшается существенно, а масса испарившейся жидкости составляет относительно небольшую величину. Один из примеров расчета указанных характеристик будет приведен ниже.
5. Испарение не кипящей жидкости
В результате охлаждения прилегающих к жидкости слоев грунта поток тепла (4) от подстилающей поверхности к объему жидкости с течением времени уменьшается, и начиная с некоторого момента его станет недостаточно, чтобы поддерживать кипение. Температура жидкости начнет падать. Скорость уменьшения температуры и скорость испарения будут зависеть от уноса ветром пара с поверхности жидкости. При падении парциального давления пара ниже его равновесного значения, жидкость будет быстро испаряться, восстанавливая с высокой скоростью утерянное равновесие. И чем интенсивнее будет унос пара, тем быстрее будет происходить испарение жидкости. Количественно изменение со временем температуры и массы испаряющейся жидкости определяется путем численного интегрирования уравнения (1), учитывающего приток тепла от всех источников и сток энергии.
6. Результаты расчета
В качестве примера рассмотрим испарение хлора, пролитого на поверхность бетона в результате аварии в трех различных случаях, произошедших 10 июля, 25 октября и 10 января при температуре окружающей среды 303,15 К, 283,15 К и 263,15 К соответственно. Первоначально хлор находился в жидком состоянии в закрытом резервуаре при высоком
давлении и температуре окружающей среды. Другие параметры задачи во всех случаях имели следующие значения: объем хлора - 1 м3, начальная масса хлора - 1427 кг, площадь обва-ловки - 20 м2, скорость ветра - 3 м/с, облачность - 0,5 бала, подстилающая поверхность -тяжелый бетон с коэффициентом теплопроводности 1,3 Дж/(с-м-К) и коэффициентом температуропроводности 5,9-10-7 м2/с. После разгерметизации и уменьшения давления жидкость оказывалась в перегретом состоянии и бурно вскипала. В результате этого часть ее практически мгновенно выбрасывалась в атмосферу в виде газа и аэрозоля. В проливе, в соответствии с указанными выше температурами воздуха в разные дни оказывалось 0,853; 1,019; 1,197 т жидкости, а время кипения этих масс определялось по формуле (8) и составляло 140; 48,6; 12,3 с. После окончания кипения жидкости происходило ее испарение в не кипящем состоянии.
На рис. 1, 2 представлены рассчитанные для различных условий зависимости массы и температуры жидкости в проливе от времени. Расчет производился с начала кипения жидкости при постоянной температуре и описывал испарение в течение 1 ч, если жидкость не испарялась полностью раньше этого времени. Из рис. 2 видно, что время кипения жидкости мало и при выбранном масштабе прослеживается только для кривой 1.
т, тонн
0.5
Т, К
240
230
220
Температура
©- 1
$- 2
V- 3
г, час 210
г, час
0.5
1
0.5
1
0
1
Рисунок 1 - Зависимость от времени Рисунок 2 - Зависимость от времени
массы испаряющейся в проливе температуры испаряющейся в проливе
жидкости при различных атмосферных жидкости при различных атмосферных
условиях, наблюдаемых 10 июля условиях, наблюдаемых 10 июля
(кривая 1), 25 октября (кривая 2) (кривая 1), 25 октября (кривая 2)
и 10 января (кривая 3) и 10 января (кривая 3)
На рис. 3 продемонстрирован вклад различных источников тепла в процесс испарения жидкости. Для этого проведено сопоставление зависимости от времени потока тепла QL = ^ Ь, затрачиваемого на испарение пролитой жидкости, с потоком тепла Qg, подводимым к объему жидкости от подстилающей поверхности, а также с суммой потоков излучения Солнца, атмосферы и потока тепла из окружающего воздуха Qs = (1 - а) Qsw +QLw + Qa. Кривой 1 обозначен поток тепла, затраченного на испарение жидкости, кривой 2 - поток тепла от подстилающей поверхности, кривой 3 - сумма потоков излучения Солнца, атмосферы и потока тепла из окружающего воздуха. В качестве исходных параметров определены условия, относящиеся к 10 июля. Начальными точками на кривых 1, 2 выбраны значения функций при г = 0,0056 ч.
Q, кВт/м2
Потоки тепла
е— i е— 2
10
0
т t, час
0
0.5
Рисунок 3 - Зависимость от времени потоков тепла, участвующих в испарении жидкости
На начальной стадии наибольший вклад в испарение жидкостей с низкой температурой кипения вносит поток тепла от подстилающей поверхности. С течением времени этот поток ослабевает, и его начинает превосходить суммарный поток от других источников.
На основе данной методики была разработана программа, которая позволяет редактировать входные данные, отображать на экране компьютера и сохранять в файле результаты расчета. Полученные значения параметров могут использоваться в качестве граничных условий для моделирования распространения паров жидкости в атмосфере.
1. Методические указания по оценке последствий аварийных выбросов опасных веществ (РД 03-26-2007): http://www.complexdoc.ru/ntd/537889
2. Тихонов, А.Н. Уравнения математической физики / А.Н. Тихонов, А.А. Самарский. М., 1966. - 724 с.
3. Methods for the calculation of physical effects (Yellow Book). - Netherlands : Minsterie van Verkeer en Waterstaat, 2005. - 870 p.
4. Степаненко, В.М. Численная модель процессов тепловлагообмена в системе водоем-почва и ее валидация / В.М. Степаненко // Вычислительные технологии. - 2004. - Т. 9, ч. 1. - С. 112-122.
5. Матвеев, Л.Т. Физика атмосферы / Л.Т. Матвеев. - СПб. : Гидрометеоиздат. 2000. - 778 с.
6. Idso, S.B. A set of equations for full spectrum and 8-14 m and 10.5-12.5 m thermal radiation from cloudless skies / S.B. Idso // Wat. Resour. Res.. -1981. - Vol. 17. - P. 295-304.
7. Kawamura, P.I. The Evaporation of volatile liquids / P.I. Kawamura, D. MacKay // J. of Hazardous Materials. - 1987. - Vol. 15. - P. 365-376.
8. Studer, D.W. Vaporization and Dispersion Modeling of Contained Refrigerated Liquid Spills / D.W. Studer, B.A. Cooper, L.C. Doelp // Plant/Operations Progress. - 1998. - Vol. 7, №. 2.
Литература
