Заглавная страница Избранные статьи Случайная статья Познавательные статьи Новые добавления Обратная связь FAQ Написать работу КАТЕГОРИИ: АрхеологияБиология Генетика География Информатика История Логика Маркетинг Математика Менеджмент Механика Педагогика Религия Социология Технологии Физика Философия Финансы Химия Экология ТОП 10 на сайте Приготовление дезинфицирующих растворов различной концентрацииТехника нижней прямой подачи мяча. Франко-прусская война (причины и последствия) Организация работы процедурного кабинета Смысловое и механическое запоминание, их место и роль в усвоении знаний Коммуникативные барьеры и пути их преодоления Обработка изделий медицинского назначения многократного применения Образцы текста публицистического стиля Четыре типа изменения баланса Задачи с ответами для Всероссийской олимпиады по праву Мы поможем в написании ваших работ! ЗНАЕТЕ ЛИ ВЫ?
Влияние общества на человека
Приготовление дезинфицирующих растворов различной концентрации Практические работы по географии для 6 класса Организация работы процедурного кабинета Изменения в неживой природе осенью Уборка процедурного кабинета Сольфеджио. Все правила по сольфеджио Балочные системы. Определение реакций опор и моментов защемления |
Метод расчета интенсивности теплового излучения при пожарах проливов лвж и ГЖСодержание книги
Похожие статьи вашей тематики
Поиск на нашем сайте
В. 1 Интенсивность теплового излучения q, кВт/м2, рассчитывают по формуле q = Ef · Fq · t, (B.1) где Ef — среднеповерхностная плотность теплового излучения пламени, кВт/м2; Fq — угловой коэффициент облученности; t — коэффициент пропускания атмосферы. В.2 Ef принимают на основе имеющихся экспериментальных данных. Для некоторых жидких углеводородных топлив указанные данные приведены в таблице В. 1.
Таблица B.1— Среднеповерхностная плотность теплового излучения пламени в зависимости от диаметра очага и удельная массовая скорость выгорания для некоторых жидких углеводородных топлив
При отсутствии данных допускается Ef принимать равной 100 кВт/м2 для СУГ, 40 кВт/м2 для нефтепродуктов. 8.3 Рассчитывают эффективный диаметр пролива d, м, по формуле , (В.2) где S — площадь пролива, м2. 8.4 Рассчитывают высоту пламени Н, м, по формуле , (В.3) где т — удельная массовая скорость выгорания топлива, кг/(м · с); rв — плотность окружающего воздуха, кг/м3; g— ускорение свободного падения, равное 9,81 м/с2. 8.5 Определяют угловой коэффициент облученности Fq по формуле , (В.4) где ,(В.5) где А = (h2 + + 1) / 2S1 , (В.6) Sl = 2 r/d (r— расстояние от геометрического центра пролива до облучаемого объекта), (В. 7) h = 2H/d; (B.8) , (В.9) B = (1+S2 ) / (2S), (B.10) B.6 Определяют коэффициент пропускания атмосферы t по формуле t = exp[ -7,0 · 10 -4 (r - 0,5 d)] (B.11) Пример — Расчет теплового излучения от пожара пролива бензина площадью 300 м2 на расстоянии 40 м от центра пролива. Расчет Определяем эффективный диаметр пролива d по формуле (В. 2) м. Находим высоту пламени по формуле (В.3), принимая т = 0,06 кг / (м2 · с), g = 9,81 м/с2 и rв = 1,2 кг/м3: Находим угловой коэффициент облученности Fq по формулам (В.4) — (В. 10), принимая r = 40 м: h = 2 · 26,5 / 19,5 = 2,72, S1 =2 · 40 / 19,5= 4,10, А = (2,722 + 4,102 + 1) / (2 · 4,1) = 3,08, B = (1 + 4,12) / (2 · 4,1) =2,17, Определяем коэффициент пропускания атмосферы т по формуле (В. 11) t = exp [ - 7,0 · 10 -4 (40 - 0,5 · 19,5)] = 0,979. Находим интенсивность теплового излучения q по формуле (В.1), принимая Еf= 47 кВт/м2 в соответствии с таблицей В. 1: q = 47 · 0,0324 · 0,979 = 1,5 кВт/м2.
ПРИЛОЖЕНИЕ Г (рекомендуемое)
МЕТОД РАСЧЕТА РАЗМЕРОВ ЗОН РАСПРОСТРАНЕНИЯ ОБЛАКА ГОРЮЧИХ ГАЗОВ И ПАРОВ ПРИ АВАРИИ Г.1 Сущность метода В настоящем приложении установлен порядок расчета изменения во времени концентрации газа в облаке при мгновенном выбросе и непрерывном истечении сжиженного углеводородного газа (СУГ), плотность которого больше плотности воздуха. Г.1.1 Мгновенный выброс СУГ Г. 1.1.1 Мгновенный выброс СУГ может происходить при повреждении резервуара или иного аппарата, в котором СУГ находится под давлением. За счет внутренней энергии СУГ его массовая доля d мгновенно испаряется, образуя с капельками жидкости облако аэрозоля. За счет больших скоростей вихревых потоков происходит быстрое вовлечение в облако воздуха и быстрое испарение оставшейся части СУГ. Массу воздуха Ма0, кг, мгновенно вовлекающуюся в облако для такого испарения, рассчитывают по формуле Ма0 = (l - d) MgLg / (Cp.a (Ta - Tg) + XwLw), (Г.1) где Мg — масса выброшенного СУГ, кг; Ср.a — удельная теплоемкость воздуха, Дж/(кг·К); Lg — удельная теплота парообразования СУГ, Дж/кг; Ta — температура окружающего воздуха, К; Тg — температура кипения СУГ при атмосферном давлении, К; Хw — массовая доля водяных паров в воздухе; Lw — удельная теплота парообразования воды, Дж/кг. d определяют из соотношения d = 1 - ехр (- Сp.g (Ta - Tg) / Lg), (Г.2) где Cp.g — удельная теплоемкость СУГ, Дж/(кг · К). Г. 1.1.2 Принимают, что образовавшееся облако дрейфует по ветру со скоростью vd = 0,6nв (nв — скорость ветра) и имеет в начальный момент форму цилиндра, высота которого равна его радиусу. С течением времени высота облака уменьшается, а радиус растет. Изменение во времени радиуса, высоты облака и концентрации газа в нем в этой фазе (называемой фазой падения) определяется путем решения методом Рунге-Кутта (реализованным в виде стандартной программы на ЭВМ) системы обыкновенных дифференциальных уравнений: dMa / dt = ra p r2 a2 a3 nв Ri-1 + 2 ra а1 (dr / dt) p r h, dT / dt =((dMa / dt) Cp.a (Ta - T) + p r2 (Tgr - T)1,333) / (Ma Cp.a + Mg Cp.g), (Г.З) dr / dt = a4 (gh (rg.a - ra) / rg.a) 0,5, где Ma — масса воздуха в облаке, кг; ra — плотность воздуха, кг/м3; r — радиус облака, м; а1, a2, a3, a4 — коэффициенты (а1 = 0,7, а2 = 0,5, a4 = 1,07, a3 = 0,3 для классов устойчивости А—В (классы устойчивости даны по Паскуиллу, таблица Г.1); 0,24 — для С—В; 0,16 — для E—F); Ri — число Ричардсона, определяемое из соотношения Ri = (5,88 h0,48 g / (a32 nв2)) (rg.a - ra) / ra; h — высота облака, м; Т— температура облака. К; Тgr — температура земной поверхности. К; rg.a — плотность паровоздушного облака, кг/м3.
Таблица Г.1— Классы устойчивости атмосферы по Паскуиллу
Решением системы вышеуказанных уравнений являются зависимости Ma = Ma (t), Т= Т(t), r= r(t). Для решения системы уравнений необходимы дополнительные соотношения rg.a = (Ma + Mg) / (Ma / ra + Mg / rg) (Ta / T). (Г.4) В качестве критерия окончания фазы падения принимают выполнение условия: (rg.a - ra) / rg.a < 10-3. (Г.5) Зависимость h = h(t) находим из соотношения h(t)=(Ma / ra + Mg / rg ) (T /Ta)(1/(p r(t)2). (Г.6) Г. 1.1.3 Когда плотность паровоздушного облака незначительно отличается от плотности воздуха (т. е. после окончания фазы падения), его движение определяется как фаза пассивной дисперсии и описывается процессами турбулентной диффузии. Концентрацию газа в точке с координатами ( х, у, z) в фазе пассивной дисперсии определяют из формулы (Г.7) где sy, sz — среднеквадратичные отклонения, зависящие от величины xc - x0; хc — координата центра облака в направлении ветра, м x0 — координата точки окончания фазы падения, м; sy (xc - x0); sz (xc - x0) зависят от класса устойчивости по Паскуиллу. При xc = x0 принимается sy0 = r / 2,14, sz0 = h / 2,14; при xc > x0 Г.1.2 Непрерывное истечение СУГ Для описания непрерывного истечения СУГ из резервуаров или иных аппаратов предполагается, что результирующая концентрация газа в паровоздушном облаке является суммой концентраций от отдельных элементарных газовых объемов и рассчитывается по формуле , (Г.8) где Q = т· t j,— масса СУГ в j -м элементарном объеме, кг; т — массовая скорость истечения СУГ, кг/с; x j— координата центра j -го элементарного объема, м; — среднеквадратичные отклонения распределения концентраций в j -м элементарном объеме, м. - определяют аналогично в Г. 1.1.3. Пример — Расчет динамики паровоздушного облака в открытом пространстве Для расчета динамики паровоздушного облака (движения в пространстве границы облака, определяемой НКПВ) принимается, что в некоторый момент времени t0 начинается истечение пропана с массовой скоростью 1,3 кг/с, скорость ветра составляет 1 м/с, градиент температуры составляет 0,667 К/м. Процедура расчета, реализованная на ПЭВМ, представлена на блок-схеме (рисунок Г.1). Результаты расчета границы облака для двух значений времени t0 + 10 с и t0 + 300 с представлены на рисунке Г.2.
Рисунок Г. 1 — Алгоритм расчета параметров паровоздушного облака t0 — время начала истечения Рисунок Г. 2 — Границы паровоздушного облака по НКПВ на различные моменты времени от начала истечения
ПРИЛОЖЕНИЕ Д (рекомендуемое)
|
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Последнее изменение этой страницы: 2016-04-20; просмотров: 511; Нарушение авторского права страницы; Мы поможем в написании вашей работы! infopedia.su Все материалы представленные на сайте исключительно с целью ознакомления читателями и не преследуют коммерческих целей или нарушение авторских прав. Обратная связь - 18.191.165.192 (0.011 с.) |