27
F
|
2
|
Ясная ночь
|
Инверсия (отрицательный градиент)
|
Решением системы вышеуказанных уравнений являются зависимости
Ma = Ma(t), Т= Т(t), r= r(t).
Для решения системы уравнений необходимы дополнительные соотношения
ρg.a = (Ma + Mg ) / (Ma / ρa + Mg / ρg ) ( Ta / T ). (Г.4)
В качестве критерия окончания фазы падения принимают выполнение условия:
(ρg.a - ρa ) / ρg.a < 10-3 . (Г.5)
Зависимость h = h(t) находим из соотношения
h(t)=(Ma / ρa + Mg /ρg ) (T /Ta )( 1/(π r(t)2). (Г.6)
Г. 1.1.3 Когда плотность паровоздушного облака незначительно отличается от плотности воздуха (т. е. после окончания фазы падения), его движение определяется как фаза пассивной дисперсии и описывается процессами турбулентной диффузии.
Концентрацию газа в точке с координатами (х, у, z) в фазе пассивной дисперсии определяют из формулы
(Г.7)
где σy, σz — среднеквадратичные отклонения, зависящие от величины xc - x0 ;
хc — координата центра облака в направлении ветра, м
x0 — координата точки окончания фазы падения, м;
σy (xc - x0); σz (xc - x0) зависят от класса устойчивости по Паскуиллу.
При xc = x0 принимается σy0 = r / 2,14, σz0 = h / 2,14;
при xc > x0 
Г.1.2 Непрерывное истечение СУГ
Для описания непрерывного истечения СУГ из резервуаров или иных аппаратов предполагается, что результирующая концентрация газа в паровоздушном облаке является суммой концентраций от отдельных элементарных газовых объемов и рассчитывается по формуле
, (Г.8)
где Q = т· τ j,— масса СУГ в j-м элементарном объеме, кг;
т — массовая скорость истечения СУГ, кг/с;
xj— координата центра j-го элементарного объема, м;
— среднеквадратичные отклонения распределения концентраций в j-м элементарном объеме, м.
- определяют аналогично 
в Г. 1.1.3.
Пример — Расчет динамики паровоздушного облака в открытом пространстве
Для расчета динамики паровоздушного облака (движения в пространстве границы облака, определяемой НКПВ) принимается, что в некоторый момент времени t0 начинается истечение пропана с массовой скоростью 1,3 кг/с, скорость ветра составляет 1 м/с, градиент температуры составляет 0,667 К/м.
Процедура расчета, реализованная на ПЭВМ, представлена на блок-схеме (рисунок Г.1). Результаты расчета границы облака для двух значений времени t0 + 10 с и t0 + 300 с представлены на рисунке Г.2.