Помощь в учёбе, очень быстро...
Работаем вместе до победы

Разработка математических моделей и методов решения комплекса задач по горению природного газа

Диссертация: Разработка математических моделей и методов решения комплекса задач по горению природного газа
ДиссертацияПомощь в написанииУзнать стоимостьмоей работы

Разработана математическая модель светящегося факела, увеличение светимости которого достигается за счет термического разложения части природного газа (2−10%) в камере предварительного горения и сопровождаемая увеличением радиационного теплового потока. Термический крекинг газа уменьшает температуру и снижает концентрацию оксидов азота в продуктах сгорания. Данная модель способствовала созданию… Читать ещё >

Содержание

  • Глава 1. Методы расчета параметров химического равновесия много компонентных, гетерогенных термодинамических систем
    • 1. 1. Комбинированный метод расчета параметров равновесия многокомпонентных термодинамических систем
      • 1. 1. 1. Система уравнений
      • 1. 1. 2. Алгоритм и программа расчетов
      • 1. 1. 3. Расчет термодинамических свойств реагирующей смеси
      • 1. 1. 4. Точность расчета равновесного состава
      • 1. 1. 5. Результаты расчетов
      • 1. 1. 6. Анализ термодинамических характеристик продуктов сгорания
    • 1. 2. Универсальный термодинамический метод расчета параметров химического равновесия многокомпонентных, гетерогенных систем
    • 1. 3. Алгоритм и программа расчета равновесных параметров термодинамических систем
      • 1. 3. 1. Алгоритм расчета
      • 1. 3. 2. Описание программы и краткая инструкция
    • 1. 4. Примеры расчетов и сопоставление с экспериментальными данными

Разработка математических моделей и методов решения комплекса задач по горению природного газа (реферат, курсовая, диплом, контрольная)

В ближайшие годы в России может сложиться крайне напряженная ситуация с обеспечением страны топливом и энергией. Резко упала добыча нефти в связи с истощением ее запасов, проводилась ошибочная политика по свертыванию угледобычи, закрываются шахты даже в тех регионах, где запасы каменного угля огромны и добыча его может продолжаться в течение многих десятков и даже сотен лет. Газовая промышленность, на которую ориентировался весь топливно-энергетический комплекс, также будет не в состоянии покрыть ожидаемый дефицит углеводородного топлива.

Удельный вес природного газа в топливно-энергетическом балансе страны в настоящее время превысил 50% при снижении доли нефти до 31%, угля до 13%), гидрои атомной электроэнергетики до 6%. Ориентация только на один вид энергоносителя — природный газ, недопустима, т.к. это создает угрозу энергетической безопасности России.

Газосбережение является одним из основных путей преодоления дефицита топливно-энергетических ресурсов, наряду с изменением структуры топливно-энергетического баланса и источником газификации коммунально-бытового сектора экономики.

Экономия топливно-энергетических ресурсов определяется в первую очередь снижением удельных расходов ТЭР на единицу продукции. Несмотря на спад производства в нашей стране, который привел к повышению удельных расходов топлива на единицу продукции, потенциальная экономия природного газа за счет снижения удельных расходов топлива до уровня уже достигнутого в развитых странах, может о составить более 22,6 млрд. м (Таблица 1).

ПОТЕНЦИАЛЬНАЯ ЭКОНОМИЯ ПРИРОДНОГО ГАЗА ЗА СЧЕТ СНИЖЕНИЯ УДЕЛЬНЫХ РАСХОДОВ ТОПЛИВА ДО УРОВНЯ В.

РАЗВИТЫХ СТРАНАХ.

Таблица 1.

Сектор РФ, млрд. м3%.

Энергетика 2,6 11,5.

Промышленность и другие отрасли 14,6 64,6.

Коммунально-бытовой сектор 5,4 23,9.

Всего 22,6 100,00.

Этого объема газа хватило бы для газификации всего коммунально-бытового сектора в соответствии с Программой газификации России.

Расчеты показали, что на реализацию мероприятий по экономии потенциально возможных объемов газа (22,6 млрд. м) потребуется на 30−50 млрд руб. меньше капитальных вложений, чем на обеспечение их добычи и транспорта.

Для достижения удельных расходов топлива, приближающихся к мировым показателям, необходима коренная модернизация газоиспользующих агрегатов: котлов, промышленных печей, отопительных системразработка и выпуск современных газосберегающих горелочных устройствприостановление перевода угольных ТЭС на природный газвнедрение новых приборов учета расхода газасовершенствование технологических процессов и агрегатов.

В стране нужна новая энергетическая стратегия. Необходимо ликвидировать сложившуюся диспропорцию между удельным весом различных энергоносителей: газ, уголь, нефть, гидрои атомная энергетика. Природный газ надо использовать квалифицированно — в промышленности на технологию, в химическом производстве, в коммунально-бытовом секторе и т. д.

Гарантированное обеспечение народного хозяйства энергоресурсами может быть решено за счет: экономии энергоресурсовреструктуризации топливно-энергетического баланса в сторону угляподачи газа наиболее экономически эффективным производствамразработке нормативных и регламентирующих документов по эффективному использованию газа.

Совершенствование процессов сжигания топлива и создание современных газогорелочных устройств невозможно без тщательного изучения всех закономерностей открытого и ограниченного факела, сопровождаемого физико-химическими явлениями горения, теплопередачи от факела в окружающее пространство, процессами термодинамики и его аэродинамикой.

Рассматриваются турбулентные, диффузионно-кинетические факела, горящие как на открытом воздухе, так и в камерах сгорания различного назначения. В этих факелах процесс горения определяется диффузией окислителя и горючего в зону, где осуществляется воспламенение и горение этих компонентов. Горение высокоскоростных струй имеет свои особенности, в частности, данные факела оторваны от газовой трубы.

Распределение интегрального химического для таких факелов определяется по методике ВНИПИгаздобыча, однозначно связывающее величину недожога топлива с продольной координатой.

Конвективная теплоотдача определяется градиентом температур, критериями Нуссельта и Грассгофа (для свободного факела) — для ограниченных факелов — критериями Прандтля, Рейнольдса, Нуссельта. Расчет указанных критериев предусматривает определение теплофизических свойств реагирующей смеси (теплоемкости, теплопроводности, вязкости).

Излучение факела определяется его температурой и степенью черноты, зависящей от состава продуктов сгорания, величины оптического луча и гидродинамики струи. Осевая составляющая радиационного потока исключается из рассмотрения из-за ее незначительности.

Математическое моделирование открытых и ограниченных промышленных факелов и создание инженерных методов их расчета является одним из инструментов решения выше указанных задач и весьма актуальна в настоящее время, т.к. наличие универсальных и достаточно гибких компьютерных программ позволит существенно экономить средства на разработку нового газоиспользующего оборудования.

Существующие проблемы экологического и теплового воздействия факела на окружающее пространство не могут быть решены только техническими средствами, а требуют детального изучения теплофизического и химического механизмов его образования. При решении указанной проблемы одним из важных этапов является определение теплофизических параметров факела, в частности, определение полного состава образующихся продуктов сгорания, длины факела, выгорания топлива, температуры, теплотехнических и теплофизических свойств.

Создание инженерных методов расчета факелов, моделирующих высокотемпературные процессы горения, теплообмена существенным образом повысит уровень создаваемых новых газогорелочных устройств с заданными характеристиками, сократит затраты и сроки на их разработку.

В работе была создана универсальная термодинамическая программа, реализующая теплообмен с окружающей средой. Использование данной программы, моделирующей высокотемпературные процессы химического реагирования гетерогенных систем, обеспечивает экономию материальных средств при оценках возможностей технологических процессов.

Создан оригинальный метод послойного расчета открытого и ограниченного турбулентного диффузионно-кинетического пламени, используя который выполнены расчеты выгорания открытого факела для всех крупных газовых и газоконденсатных месторождений России, послужившие основой при создании отраслевой методики расчета факельных установок.

Разработан метод расчета горелок общего и специального назначения, позволивший создать оптимальную конструкцию вихревой горелки для камер сгорания ГПА с пониженным содержанием оксидов азота (25−32 мг/м3 при 15% Ог). Данный метод рекомендуется для расчета вихревых горелок различной мощности.

Создана математическая модель газового факелов, подверженного ветровой нагрузке, для определения длины факела и его теплового воздействия на окружающее пространство с учетом выгорания топлива. Разработанный метод позволяет определить пространственную зону теплового отчуждения вокруг факела в пределах которой человек может быть подвергнут сильному термическому воздействию, связанному с перегревом и последующим биохимическими изменениями верхних слоев кожного покрова.

Данный метод расчета может быть использован при проектировании и усовершенствовании факельных установок. В частности, сформулированы следующие технологические рекомендации для работы факельных горелок:

— при скоростях топливного газа менее 0,2 М на факельной установке должна быть обеспечена подача горячего водяного пара, что исключит наличие в продуктах сгорания несгоревших углеводородов, включая сажу.

— для безопасного сжигания сбросных газов на факельных установках при скоростях газа от 0,2 до 0,8 М горелка должна содержать стабилизационное устройства типа сопла Коандо или Вентури, наличие которого исключает срыв факела, обеспечивает экономию топливного газа на 60% и в 1,6 раза уменьшает содержание токсичных компонентов с продуктами сгорания и плату за эти выбросы от постоянно сжигаемого продувочного газа.

Создана оригинальная математическая модель, алгоритм и программа расчета геометрических размеров факелов жидких углеводородов, сжигаемых на открытом воздухе, и расчета облучения окружающих предметов такими факелами, используя который для различных значений скоростей ветра были рассчитаны значения «критического» расстояния, являющегося зоной отчуждения, необходимые для безопасной работы обслуживающего персонала.

Разработана математическая модель «самокарбюрации» факела, реализуемая в камере предварительного горения при 2−10% природного газа и сопровождаемая увеличением радиационного теплового потока по сравнению с прозрачным факелом. Уменьшение температуры факела снижает концен.

12 трацию оксидов азота в продуктах сгорания. Используя данную модель созданы новые конструкции газогорелочных устройств (ГПГ, ГПГ1), последняя из которых работает на газокислородном топливе. Применение новой технологии отопления высокотемпературных печей, основанной на газокислородном отоплении, обеспечивает экономию природного газа на 20%, снижает содержание оксидов азота на 60% по сравнению с газо-воздушным отоплением с регенератором, увеличивает удельную производительность теплового агрегата не менее чем на 30% за счет роста радиационной составляющей теплового потока, позволяет уменьшить площадь, занимаемую тепловым агрегатом в 2−3 раза.

Выполнен анализ радиационно-конвективной теплопередачи в агрегатах, оснащенных скоростными горелками. Показано наличие экстремальных значений суммарного радиационно-конвективного теплового потока при изменении степени недожога и определенных давлениях в камере сгорания. Проанализирован выход равновесных и замороженных компонентов продуктов сгорания из скоростной горелки с различными условиями истечения. Разработанный метод позволил создать новую скоростную горелку, обеспечивающую экономию топлива на 20%. Удалось снизить содержание оксидов азота в продуктах сгорания до 90 мг/м3. Моделирование снизило затраты на разработку подобной горелки более чем на 30%.

ОБЩЕЕ ЗАКЛЮЧЕНИЕ.

В результате выполненной работы можно сделать следующие выводы:

1. Разработана универсальная термодинамическая программа, в которой реализован теплообмен с окружающей средой. Наличие данной программы, моделирующей высокотемпературные процессы химического реагирования гетерогенных систем, обеспечивает существенную экономию средств при оценках технических возможностей различных технологических процессов.

2. Создан оригинальный метод послойного расчета открытого и ограниченного турбулентного диффузионно-кинетического пламени, используя который выполнены расчеты выгорания открытого факела для всех крупных газовых и газоконденсатных месторождений России. Результаты расчетов и указанный метод расчета были использованы при создании отраслевой методики расчета факельных установок, которая используется при проектировании факельных устройств, включая факельные горелки.

3. Создана математическая модель газового факелов, подверженного ветровой нагрузке, для определения длины факела и его теплового воздействия на окружающее пространство с учетом выгорания топлива. Разработанный метод позволяет определить пространственную зону теплового отчуждения вокруг факела, в пределах которой объект может быть подвергнут сильному термическому воздействию.

4. Сформулированы следующие технические требования к факельным установкам:

— при скоростях топливного газа менее 0,2 М на факельной установке должно быть установлено устройство подачи горячего водяного пара, что исключает появление в продуктах сгорания несгоревших углеводородов, включая сажу;

— для безопасного сжигания сбросных газов на факельных установках при скоростях газа от 0,2 до 0,8 М факельная горелка должна содержать стабилизационное устройство типа сопла Коандо или Вентури, что исключает срыв факела, обеспечивает экономию топливного газа на 60%, уменьшает содержание токсичных компонентов с продуктами и плату за выбросы от постоянно сжигаемого продувочного газа в 1,6 раза.

5. Создана оригинальная математическая модель, алгоритм и программа расчета геометрических размеров факела жидких углеводородов и расчета облучения таким факелом окружающих предметов, что позволяет определить зону теплового отчуждения вокруг факела при различной скорости ветра.

6. Моделирование горелок общего и специального назначения позволило создать оптимальную конструкцию вихревой горелки для камер сгорания ГПА с пониженным содержанием оксидов азота (25−32 мг/м3 при 15% Ог).

7. Разработана математическая модель светящегося факела, увеличение светимости которого достигается за счет термического разложения части природного газа (2−10%) в камере предварительного горения и сопровождаемая увеличением радиационного теплового потока. Термический крекинг газа уменьшает температуру и снижает концентрацию оксидов азота в продуктах сгорания. Данная модель способствовала созданию новых конструкций газого-релочных устройств (ГПГ, ГПГ1), последняя из которых, использующая газокислородное топливо, обеспечивает экономию природного газа на 20%, снижает содержание оксидов азота по сравнению с газо-воздушным агрегатом с регенератором на 60%, увеличивает его удельную производительность на 30%, позволяет уменьшить площадь печи в 2−3 раза.

8. Выполнен анализ радиационно-конвективной теплопередачи в агрегатах, оснащенных скоростными горелками. Метод расчета позволил создать новую конструкцию скоростной горелки, обеспечивающую экономию топлива на 20% по сравнению с горелкой ГНП, снизить содержание оксидов азота в продуктах сгорания 90 мг/м3) за счет многостадийного сжигания топлива.

Показать весь текст

Список литературы

  1. Термодинамические и теплофизические свойства продуктов сгорания, / В. Е. Алемасов, А. Ф. Дрегалин, А. П. Тишин, В. А. Худяков. М.: ВИНИТИ, АН СССР, т. 1, 1971.-266 с.
  2. И.Б. Рождественский, В. П. Шевелев, К. К. Олевинский. Расчет состава и термодинамических функций произвольных реагирующих газовых систем. Сб. «Свойства газов при высоких температурах». Наука, М., 1967.
  3. Г. Б. Синярев, Н. А. Ватолин, Б. Г. Трусов, Т. К. Моисеев. Применение ЭВМ для термодинамических расчетов металлургических процессов. М.: Наука, 1982. — 263с.
  4. Е.А. Термодинамический метод расчета равновесных химически реагирующих систем // Газовая промышленность, 1997 г, № 9, с. 56-^57.
  5. Н.А. Ватолин, Г. К Моисеев, Б. Г. Трусов Термодинамическое моделирование в высокотемпературных неорганических системах. М.: Металлургия, 1994. — 352с.
  6. Е.А. Компьютерная программа для создания позонного метода расчета факелов // Газовая промышленность, 1998 г, № 7, с. 59-К31.
  7. В.М., Шибаева Л. И. //Инженерная методика расчета длины факела и выгорания газообразного топлива применительно к горелкам общего назначения, г. Саратов: ВНИПИГаздобыча, 1981. — 24с.
  8. Термодинамические свойства индивидуальных веществ. Справочное издание: в 4-х т./ Л. В. Гурвич, И. В. Вейц, В. А. Медведев и др. М.: Наука, 1982.
  9. JANAF Thermochemical tables: 2-nd edition. NSRDS-NBS 37. -Waschington: US Gov. Print. Office, 1971. 1141p.
  10. Schick H.L. Thermodinamics of certain refractory compounds, v. 1−2. N.-Y. — London Acad. Press, 1966. V. 1. — 632 p.- v. 2. — 775 p.
  11. А.П. Термодинамический расчет горения и истечения в области высоких температур, М., БНТ, 1947.
  12. В.П. Жидкое топливо для реактивных двигателей. М. ВВА им. Н. Е. Жуковского, 1936.
  13. Я.Б., Полярный А. И. Расчеты тепловых процессов при высокой температуре. М., БНТ, 1947.
  14. В.Е.Алемасов, А. Ф. Дрегалин, А. П. Тишин / Теория ракетных двигателей, М., Изд-во «Машиностроение», 1969.
  15. А.В. / Расчет процессов в камере сгорания и сопле жидкостного ракетного двигателя. М., Оборонгиз, 1957.
  16. Я.Б., Ривин М. А., Франк-Каменецкий Д.А. / Импульс реактивной силы пороховых ракет. М., Оборонгиз, 1963.
  17. С. Химия ракетных топлив. М., Изд-во «Мир», 1969.
  18. Р.Е. Газотермодинамика ракетных двигателей на твердом топливе. М., Изд-во «Наука», 1967.
  19. М.Б. Эффективность использования топлива, Изд-во «Наука», М, 1977.-343 с.
  20. Расчеты горения топлива. Нормали. М., Стальпроект, 1963.
  21. И.Н., Сорока Б. С., Дашевский Л. Н., Семернина С. Д. Продукты сгорания природного газа при высоких температурах (состав и термодинамические свойства). Киев. Изд-во «Техника», 1968.
  22. С.Л. Термодинамические свойства воздуха и продуктов сгорания топлив. М., Госэнергоиздат, 1962.
  23. Н.Р. Справочник по углеводородным топливам и их продуктам сгорания. М., Госэнергоиздат, 1962.
  24. В.Н., Борцов В. И., Сапрыкин Г. С. Таблицы термодинамических свойств продуктов сгорания саратовского природного газа при высоких температурах. Саратов, Саратовский университет, 1964.
  25. Л.К., Щеголев Г. М. Энтропийные диаграммы продуктов сгорания (до 4000°К). Киев, изд-во «Наукова думка», 1968.
  26. А.П. Термодинамический расчет горения и истечения в области высоких температур. М., БНТ, 1947.
  27. A.M., Шаулов Ю. Х. Термодинамические исследования методом взрыва и расчеты процессов горения. М., Изд-во МГУ, 1955.
  28. Внутрикамерное сжигание природного газа в кислороде и воздухе / Н. В. Лавров, Е. А. Пацков, Е. А. Плужников, Н. А. Федоров. М.: 1972 г. — 38с.
  29. И., Дефей Р. Химическая термодинамика. Наука, Новосибирск, 1966. — 510 с.
  30. М.И., Иссерлин А. С., Ежов Е. И., Рыжов О. А. Получение высокотемпературного факела с контролируемой газовой фазой // Теплоэнергетика, 1978 г, № 9, с. 38−40.
  31. М.Н., Шурыгин А. П. Огневая переработка и обезвреживание промышленных отходов. М.: Химия, 1990, 302с.
  32. С.Ф., Зайцев С. И. Рос.хим.ж. (Ж.Рос. химического общества им Д.И.Менделеева), 1995, t. XXXIX, № 4, с.24−26.
  33. С.Ф. ЖПХ, 1990, т.63, N 9, с.2088−2090.
  34. Б. Льюис, Г. Эльбе Горение, пламя и взрывы в газах. Второе издание под ред. Щелкина К. И., Изд-во «Мир», М. 1968, с. 592.
  35. Г. С., Стрекалова Л. В., Шарихина Л. В., Прокофьева А. М./ Методика расчета параметров выбросов и валовых выбросов вредных веществ от факельных установок сжигания углеводородных смесей, РАО «Газпром», ВНИИгаз, ИРЦ «Газпром». М. 1996 г., с. 47.
  36. А.В. Арсеев. Сжигание природного газа, Металлургиздат, М. 1963, с. 407.
  37. Л.А. Вулис, Ш. А. Ершин, Л. П. Ярин. Основы теории газового факела, «Энергия», ЛО, 1968, с. 204.
  38. В.Г. Аэродинамические характеристики факела в условиях действия подъемных сил. ИВУЗ «Черная металлургия», 1969, N4, с.143−149
  39. Н.В., Арсеев А. В., Китаев Б. И. Исследование закономерностей строения горящей свободной струи газа. В сб. «Теория и практика сжигания газа», Л., «Недра», 1967, вып. III, с.84-ь103.
  40. Г. Н. Теория турбулентных струй. М., Физматгиз, 1960, 716с.
  41. О. Г., Пацков Е. А. Расчет выбросов загрязняющих веществ на факельных установках // Газовая промышленность. 1997. — N 1. — С. 24−26.
  42. В.Г., Китаев Б. И., Кокарев Н. И. Усовершенствование методов сжигания природного газа в сталеплавильных печах. М., «Металлургия», 1977, 280с.
  43. С.И., Семикин И. Д. Расчет длины турбулентного газового факела. «Известия Вузов, Черная металлургия», 1965, № 4, с. 202−211.
  44. С.И., Семикин И. Д. Влияние различных факторов на длину турбулентного газового факела. «Известия Вузов, Черная металлургия», 1965, № Ю, с. 146−152.
  45. С.Л., Ривин Э. М. Методы расчета вредных выбросов в атмосферу из нефтехимического оборудования. М., ЦНИИТнефтехим, «Нефтеперерабатывающая и нефтехимическая промышленность», Серия: Охрана окружающей среды. Обзорная информация, Вып.1, 1991, 80с.
  46. Методика расчета вредных выбросов в атмосферу из нефтехимического оборудования (РМ62−91−90). Гипрокаучук, Воронежский филиал. Воронеж, 1991, 78с.
  47. Г. Н. О температурном поле жидкости, горящей со свободной поверхности, и о факеле над ней. Известия АН СССР, ОТН, № 7, 1951, с. 1015 -1024.
  48. В. И., Худяков Г. Н. Диффузионное горение жидкости. М., АН СССР, 1961,208 с.
  49. О. М. Пожарная безопасность резервуаров с нефтепродуктами. М., «Недра», 1984, 152 с.
  50. О.Г., Пацков Е. А. Расчет выбросов при горении жидких углеводородов на наземных факельных установках // Газовая промышленность, N 6, 1997 г, с. 28−29.
  51. NAO Inc., General Catalog.
  52. Основы практической теории горения: Учебное пособие для вузов/ В. В. Померанцев, К. М. Арефьев, Д. Б. Ахмедов и др.- Под ред. В. В. Померанцева. 2-е изд. переработанное и дополненное Л.: Энергоатомиздат. Ленингр. Отделение, 1986.-312 е.: ил.
  53. И. Я. Сигал, Защита воздушного бассейна при сжигании топлива. 2-е изд., переработанное и дополненное — Л.: Недра, 1988. — 312с.: ил.
  54. Р. Б. Дутьевые газогорелочные устройства. М., Недра. 1977.
  55. ГОСТ 29 134–91. Горелки газовые промышленные. Методы испытаний = Industrial gas burners. Test methods. Введен 01.01.93 — M.: Изд-во стандартов, 1992 — 22с.: ил. УДК 662. 951. 2. 001. 4:006. 354. Группа Е29(47) СССР.
  56. И.И., ЭльНатанов А.И. // Факельные установки, М.: Химия, 1979, — 184 е., с ил.
  57. W. R., Wedell D.S., Hottel Н.С. // 3-rd Symposium on Combustion and Flame and Explosion Phenomena, Baltimore, 1949, 266−288.
  58. K., Gazley C., Kapp N. // 3-rd Symposium on Combustion and Flame and Explosion Phenomena, Baltimore, 1949, p.288.
  59. K., Kapp N.M., Gazley C. // Meteor Report UAC-26, September, 1948.
  60. K., Kapp N. M., Gazley C. // 3-rd Symposium on Combustion and Flame and Explosion Phenomena, Baltimore, 1949, p. 3.
  61. D.A., Garside J.E. // 3rd Symposium on Combustion and Flame and Explosion Phenomena, Baltimore, 1949, p.102.
  62. W. Forstall. Momentum and Mass Transfer in a Submerged water Jet. // Jour. Of Applied Mechanics AS ofMEAM. 1954.
  63. H. Hottel, Burning in Laminar and Turbulent Fuel Jets // Fourth Symposium Int. On Combustion. Massachusetts Inst. Of Technology, Sep., 1952.
  64. G. Kent. Practical Desing of Flare Stacks // Hydrocarbon Processing. 1964, 43, 8, p. 121−125.
  65. W. Crocker. Assessment of mathematical models for fire and explosion hazards of LPG. // J. Of Hazardous Materials. 1988, 29, p.109−135.
  66. Г. Г., Комов В. Ф. Влияние негорючих примесей на длину присоплового осесимметричного турбулентного диффузионного факела пламени // Физика горения и взрыва, 1977, № 5, с. 667−670.
  67. API RP-521 American Petroleum Institute, Division of Refining. 1969.
  68. R. Mc. Murray. Flare radiation estimated. // Hydrocarbon Processing, 1985, May, 89−91.
  69. D. Cook. Flaring of natural gas from inclined vent stacks. // Inst. Chem. Eng. Symp. Ser., 1987, 102, p.289−300.
  70. T. Brzustowski. Turbulent Combustion (Prog. Astro and Alro) AIAA, 1978, 58.407.
  71. D.M. De Faveri. Estimate flare radiation intensity. // Hydrocarbon Processing, 1985, May, p.89−91.
  72. H. Becker, D.Liang. Visible Length of vertical free turbulent diffusion flames. // Combustion and Flame., 1978, 32, p. 115−137.
  73. G. Kalghatgi. The visible shape and size of a Turbulent Hydrocarbon Jet Diffusion Flame in a Cross-wind // Combustion and Flame, 1983, 52, 91−106.
  74. Ш. А., Ярин Л. П. Исследование диффузионных пламен. В сб.: «Прикладная теплофикация», Алма-Ата, изд-во АН КазССР, 1964, вып.1, с. 101 139.
  75. G. Chamberlin. Developments in design methods for predicting thermal radiation from flares. // Chem. Eng. Res. Des., 1987, 65, pp. 299−309.
  76. Я.Б. К теории горения не перемешанных газов. Ж.Т.Ф., 1949, т.19, вып. 10, с.1199−1121.
  77. Mudan К. Thermal radiation hazard from hydrocarbon pool fires // Property Energy and Combustion Science, 1984, 10, 1, p.59−80.
  78. И. Я. Сигал, A. M. Марковский, H. А. Гуревич, С. С. Нижник. Образование окислов азота в топках котельных агрегатов. «Теплоэнергетика», 1971,4, с. 57-г 60.
  79. И. Я. Сигал. Горение газа в котлах и атмосфера городов. «Газовая промышленность», 1969, № 2, с. 30−35.
  80. R. Е. Jeorge, R. L. Chass. «Journ. of the Air poll. Control ASS», 1967, vol 17,1 6.
  81. Я. Б. Зельдович, П. Я. Садовников, Д. А. Франк-Каменецкий. Окисление азота при горении. Изд. АН СССР, М.-Л., 1947.
  82. Ю. П. Райзер. Образование окислов азота в ударной волне при сильном взрыве воздуха. МФХ, т. XXXIII, № 3, 1959, 700−709.
  83. J. Longwell, М. Weiss. «Jnd. Eng. Chem», 47, 1955, 1634.
  84. Ellison. «Process optimization in control of air pollution». «Pap. Amer. Mech. Eng.», 1970, NWA/APC-2.
  85. Southern California Edison limits Nox with firring modification, dispatching technique. «Eles. World», 1970, 174, 1 9.
  86. Weidemann Horst, Fahring Peter, Andermann Werner. Brenner zur Beseitigung nitroser Gase aus Industrie abgasen. Пат. ГДР, кл. 24 е., 10(F 23d), № 76 567, заявл. 18. 12. 68., опубл. 5. 10. 70.
  87. James D. W. Coping with Nox: growing problem. «Eles World», 1971 175, № 3, 44−47.
  88. А. Г. Тумановский, В. А. Христич, A. M. Шевленко. Влияние типа горелочного устройства на образование окислов азота в камерах сгорания ГТУ при сжигании природного газа «Теплоэнергетика», 1970, № 5, 35−38.
  89. A. W. Bell, Valo N. de Bayard, В. F. Breen. Combustion control for elimination of nitric oxide emissions from fossil field power plant. «13th Sump. Int. Combustion Salt Lake City, Utah, 1970, Abstr. Pap», Pittsburgh, Pd, 60−62.
  90. В. А. Кривандин. Светящееся пламя природного газа. М., «Металлургия», 1973, с. 136.
  91. Е. А. Пацков, Э. И. Розенфельд, В. В. Федоров. К расчету влияния светимости газового факела на теплообмен в печи.//"Газовая промышленность", 1973, № 9
  92. Е.А., Розенфельд Э. И., Федоров Н. А. «Снижение образования оксидов азота в котлах электростанций и промышленных газовых топках.// Теплоэнергетика, 1974, № 8, с.55−60.
  93. О. Л. Магдасалиев. Экспериментальное исследование теплообмена при охлаждении высокотемпературного потока в трубе. Сб. Тепло- и массобмен в одно- и двухфазных средах. Изд. „Наука“, М., 1971, 135−140.
  94. А.С. № 859 762 Победин Ю. В., Теляков В. В., Герасимов К. И.
  95. Е.А., Плужников А. И. Метод расчета газогорелочных устройств в топках тепловых агрегатов. // Газовая промышленность, 1993 г., № 8 с. 30−32.
  96. Е.А., Розенфельд Э. И. Основы расчета горелок с высокоскоростным потоком продуктов сгорания. // Газовая промышленность, 1972 г., № 8 с. 30−34.
  97. Лиз. Конвективный теплообмен при наличии подвода вещества и химических реакций. Сб. Газодинамика и теплообмен при наличии химических реакций, пер. с англ. Изд-во иностр. лит-ры, М., 1962.
  98. Т. А. Агекян. Основы теории ошибок для астрономов и физиков. Изд. „Наука“, М., 1972.
  99. В.А., Акопова Г. С., Синицын Ю. Н., Корнеев В. И., Гладкая Н. Г. Снижение выбросов загрязняющих веществ с отходящими газами газотурбинных ГПА. М.: ВНИИЭгазпром, 1991, 61с.
  100. И .Я. //Газовая промышленность, 1969, № 2, с.30−3 5.
  101. А.В., Арсеева Н. В. Загрязнение атмосферы окислами азота продуктов сгорания топлива. Научно-технический обзор. Серия: Использование газа в народном хозяйстве, ВНИИЭгазпром, М., 1974.
  102. А.Г., Христич В. А., Шевченко A.M. „Теплоэнергетика“, 1970, N 5, с.35−38.
  103. В.Н. Константы скоростей газофазных реакций (Справочник), М., Изд-во „Наука“, 1970.
  104. Е.А., Плужников А. И., Федоров Н. А. Термодинамический анализ горения метана. Газовая промышленность, 12, 1972, с.42−46.
  105. Е. С. Щетинков. Физика горения газов. Изд. „Наука, М., 1965.
  106. Л. М., Звягинцев К. Н., Фалин А. А., Щепнова Е. С. Оптимизация двухстадийного сжигания топлива в скоростных горелках. „Газовая Промышленность“ № 11, 1983, с. 38−39.
  107. Автогенные процессы в цветной металлургии /Мечев В.В., Быстров В. П., Тарасов и др. М. Металлургия, 1991 — 413с.
  108. FK (9,100), NC (20), G (2), FR (100,9), NFR (100),
  109. ESMB (100,3), RS (100), RP (8), Q (80) DOUBLE PRECISION A, T
  110. CALL PP100P (P1, P2, L1, L2, M, EL, NO, KODI, NVAR) IF (M.eq.777) goto 99
  111. NQUE=0 rewind 17 rewind 12call PP100C (M, K, NO, KODI, MSK, FK, NC, FR, NFR, ESMB, NVAR) if (M.eq.777) go to 991. NVAR=NVAR+11. PARS=P2(1)
  112. PARS=PARS+P2(2) IF (P2(2).LT.0.0) GO TO 6 IF (PARS.GT.P2(3)) GO TO 7 GO TO 2
  113. IF (PARS.LT.P2(3)) GO TO 7 GO TO 2
  114. PARS=P2(1) IF (ABS (P1(4)-.77 777).LT.1E-4.AND.ABS (P1(3)77777).GT.1E-4) GO TO 9 DO 8 1=1,9
  115. P1(I)=P1(I+1) P1(10)=.777 771. (ABS (P1(1)-.77 777).LT.1E-4) go to 97 GO TO 2
  116. P1(1)=P1(1)+P1(2) IF (P1(2).LT.0.0) GOTO 10 IF (P1(1).GT.P1(3)) GO TO 97 GO TO 2
  117. IF (P1(1).GE.P1(3)) GOTO 2 97 WRITE (5,101)
  118. READ (*, 100) Q 97 continue read (4,100) Q1.(Q (1).NE.'H,.AND.Q (1).EQ., h') GO TO 27 IF (Q (2).NE.'E, AND. Q (2).NE.'e') GO TO 27
  119. WRITE (5,103) DO 25 1=1,74 READ (4,100) Q IF (I.LT.69) GO TO 25 WRITE (5,100) (Q (J), J=1,79)
  120. CONTINUE REWIND 4 GO TO 9727 DO 11 1=1,101.(Q (l).EQ.'Y'.OR.Q (l).EQ.y) GOTO 110 IF (Q (l).EQ.'fl,.OR.Q (l).EQ., fl') GOTO 110 IF (Q (l).EQ.'B'.OR.Q (l).EQ.'6') GO TO 11 011 CONTINUE STOPcontinue100 FORMAT (80A1)
  121. FORMAT (' Конец исходных данных. Будете про', должать работу? ')
  122. FORMAT ('- На любой вопрос можно ответить словом'
  123. НЕ1РУ5Х,'и Вам будет оказана консультация'////) 103 FORMAT (24(/))1. END
  124. SUBROUTINE РР100Р (P1,P2,L1,L2,M, ELNO, KODI, NVAR) DIMENSION np (6), EL (10), NO (10), P1 (10), P2(10), CN (82),
  125. NS (94), AME (94), COEF (5) DIMENSION KN (30), KS (400), Q (80) CHARACTER CR1 (4)1. EQUIVALENCE (CR1(1), CR4)
  126. DATA COEF/0.101 325,273.15,0.001,4.184.0.101 972/
  127. DATA NP/31,25,37,34,24,36/
  128. WRITE (5,108) READ (*, 103) Q70 continue read (4,103) Q DO 701 1=1,80 CR4=Q (I) J=ICHAR (CR1(1))1. (J.GE.97.AND.J.LE.122) J=J-32 CR1(1)=CHAR (J) Q (I)=CR4 il=0
  129. DO 700 J=1,82 IF (Q (I).EQ.CN (J)) ll=J700 CONTINUE1. (ll.GT.64) 11=11−19 701 Q (I)=CN (II)1. (Q (1).EQ.CN (24).AND.Q (2).EQ.CN (21)) GO TO 134 J=0
  130. READ (*, 103) Q 90 READ (1,103) Q WRITE (*, 112) Q1. (Q (1).EQ.CN (24).AND.Q (2).EQ.CN (21GO TO 135 J=0
  131. DO 76 1=1,80 IF (Q (I).EQ.CN (45)) J=J+176 CONTINUE1. (J.GT.1) GO TO 75 DO 77 1=1,78 JJ=79-I
  132. Q (JJ+2)=Q (JJ) IF (J.EQ.1) J3=J2 Q (1)=CN (J3) Q (2)=CN (45) GO TO 8478 KEY=11. (J3.EQ.J2) GO TO 79 J3=J2 GO TO 75 79 WRITE (5,110)79 continue80 READ (*, 103) Q
  133. READ (2,103) Q WRITE (*, 112) Q1. (Q (1).EQ.CN (24).AND.Q (2).EQ.CN (21)) GO TO 1361. KEY=31. DO 81 1=1,801. J=81-l1. (Q (J).NE.CN (11)) GO TO 8281 CONTINUE GO TO 80
  134. IF (Q (J).EQ.CN (12).OR.Q (J).EQ.CN (13)) GOTO 87 J1=0
  135. DO 83 1=1,10 IF (Q (J).EQ.CN (I)) J1=183 CONTINUE1. (Q (J).EQ.CN (14)) J1=1 IF (J1.EQ.1) Q (J+1)=CN (13) GO TO 8784 DO 85 1=1,80 J=81-l1. (Q (J).NE.CN (11)) GO TO 8685 CONTINUE GO TO 90
  136. OR.KS (l+2).EQ.35)) S=COEF (5)1. (I3.EQ.52.0R.I3.EQ.18) S=COEF (2) IF (I3.EQ.59) S=COEF (3) IF ((I3.EQ.53.0R.I3.EQ.26).AND.(KS (I+2).EQ.53
  137. M=M+1 N0(M)=I1 SS=1.0 N1=3 GOTO 17
  138. IF (N1.NE.2) GO TO 28 P2(N2)=W1. GO TO 2928 EL (M)=W29 CONTINUE GO TO 3050 WRITE (5,100) GO TO 6051 WRITE (5,101) • GO TO 6052 WRITE (5,102) GO TO 6053 WRITE (5,104) GO TO 6054 WRITE (5,105) GO TO 60
  139. WAM=0. DO 121 J=1,M LL=NO (J)
  140. WAM=WAM+EL (J)*AME (LL) WAM=1000.0/WAM1. DO 122 J=1,M
  141. FORMAT (1 OX,'Ошибка в данных: непредусмотренный символ.')
  142. FORMAT (1 Ох,'Ошибка в данных: непредусмотренный ', т.д. параметр.')
  143. FORMAT (10Х,'Ошибочный символ химического элемента.')103 FORMAT (80A1)
  144. FORMAT (1 ОХ/Отсутствуют исходные данные.')
  145. FORMAT (10х,'Ошибка в данных: непредусмотренная ', размерность.')106 FORMAT (23(/))
  146. FORMATC-- Назовите два параметра равновесия ', 1. P, T, V, S, I, U):')109 FORMAT (' -- ', А1,'=')
  147. FORMAT (2x,'Мольное содержание элементов:')
  148. FORMATC -- Нужно учитывать ионы? ')
  149. FORMATC Для продолжения нажмите (У, Д)')
  150. FORMAT (2Х, 78А1) 120 FORMAT (F10.4)1. END
  151. SUBROUTINE РР100С (M, N, NO, KODI, MSK, F, NC, FR, NFR, ESMB, NVAR)
  152. DOUBLE PRECISION SK1, SK2,SK3,SK4,SNN
  153. DIMENSION MSK (1100), F (900), NC (20), FR (100,9), ESMB (100,3),
  154. Q (80), KS (100), SNN (100,2), NFR (100), CN (51), NS (94),
  155. N1(5), SYM (3), N2(5), FF (10), NO (M), Q4(20), SNAME (4,3) CHARACTER SYM1(12), SB1(4), QLL (80)
  156. EQUIVALENCE (SYM (1), SYM1(1)),(SBQ, SB1(1)),(Q4(1), QLL (1)) =======================================================:=========CС
  157. Ввод коэффициентов полиномов, аппроксимирующих С термодинамические свойства индивидуальных веществ СС
  158. Входные параметры: M, NO, KODI, NVAK С
  159. М количество химических элементов- С
  160. NO массив номеров химических элементов Спо периодической таблице- С
  161. KODI константа, устанавливающая необходи- Смость рассмотрения ионов- С
  162. NVAR номер варианта расчета. С
  163. Выходные параметры: N, MSK, F, NC, FR, NFR, ESMB С
  164. N количество компонентов равновесия- С
  165. MSK массив стехиометрических коэффициентов- С F — массив коэффициентов, аппроксимирующих С полиномов- С
  166. NC массив программных констант- Сс FR массив дополнительных коэффициентов, Саппроксимирующих полиномов- С
  167. NFR массив номеров индивидуальных веществ, С к которым относятся соответсвующие Сдополнительные полиномы- С
  168. ESMB массив имен индивидуальных веществ- С=================С
  169. DATA CN /'О', Т,'273','4','576', 778','9',-/A'.'B'.'C'/D',
  170. E'/F'/G'/H'/l'/J'/K'/L'/M'/N', 1. OVPVQVRVS'.T/U'.VVWVX',
  171. Т.Т/Г.Т^'/Ж'.'Э'/Л'.^.'Г/З'/
  172. MSK (l)=0 IF (NVAR.EQ.1) OPEN (UN^S.FILE^COMP.FOR1) IF (NVAR.GT.1) REWIND 3 WRITE (5,232), 1. NR1=0 11. N=0 ION=0
  173. DO 99 1=1,100 99 NFR (l)=0 MA=11 K=0
  174. READ (3,113,END=501) Q DO 5 1=1,801. (Q (I).EQ.CN (11)) GO TO 5 IF (Q (I).EQ.CN (13)) GO TO 6 DO 3 J=1,511. (Q (I).EQ.CN (J)) GO TO 43 CONTINUE GO TO 2014 K=K+1 KS (K)=J-15 CONTINUE GO TO 26 DO 7 1=1,161. (KS (I).EQ.11) GO TO 87 CONTINUE GO TO 202
  175. N1(M1)=N1(M1)*50+KS (J) GO TO 2320 IF (N2(M1)) 21,21,2221 N2(M1)=KS (J) GO TO 23
  176. N2(M1)=N2(M 1)*10+KS (J) M1=M1+123 CONTINUE1. (N2(M1).EQ.O.AND.N1(M1).NE.O) N2(M1)=1
  177. CONTINUE MSK (MA)=0 MSK (MA+M+1)=0
  178. IF (KS (IS-1).EQ.15) W=-W FF (I10)=W1.(KI.EQ.K) GOTO 3736 CONTINUE
  179. IF (I10.NE.10) GO TO 204 FF (6)=FF (6)-FF (10)*1.0E-4 SK3=N1(1)
  180. SK1=SK1+SK3*100 000 000.0 SK3=N2(1)
  181. SK2=SK2+SK3*100 000 000.0 IF (N.EQ.0) GO TO 41 DO 38 1=1,N
  182. SK3=DABS (SK1-SNN (I, 1)) SK4=DABS (SK2-SNN (I, 2)) IF (SK3.LT.0.5.AND.SK4.LT.0.5) GO TO 3938 CONTINUE GO TO 41
  183. NR1=NR1+1 NFR (NR1)=I DO 40 1=1,940 FR (NR1,I)=FF (I) GOTO 1
  184. N=N+1 SNN (N, 1)=SK1 SNN (N, 2)=SK2 J=(N-1)*91. DO 42 1=1,9
  185. SNAME (J, 1)=SYM (1) SNAME (J, 2)=SYM (2) SNAME (J, 3)=SYM (3) J25=J1. (J.LT.4) GO TO 1 33 continue
  186. WRITE (5,231) N-3,(SNAME (1,KA), KA=1,3), N-2,
  187. SNAME (2,KA), KA=1,3), N-1,(SNAME (3,KA), KA=1,3),
  188. N,(SNAME (4,KA), KA=1,3) GO TO 1
  189. IF (J25.EQ.4) GO TO 301 IF (J25.EQ.1) WRITE (5,231) N-J25+1,(St, AME (1,KA), KA=1,3) IF (J25.EQ.2) WRITE (5,231) N-J25+1,(SNAME (1,KA), KA=1,3),
  190. N-J25+2,(SNAME (2,KA), KA=1,3) IF (J25.EQ.3) WRITE (5,231) N-J25+1,(SNAME (1,KA), KA=1,3),
  191. N-J25+2,(SNAME (2,KA), KA=1,3), N-J25+3,(SNAME (3,KA), 1. KA=1,3) GO TO 30 101 WRITE (5,221) Q
  192. GO TO 300 °02 WRITE (5,222) Q GO TO 30 003 WRITE (5,223) Q GO TO 30 004 WRITE (5,224) Q300 M=777 301 NC (10)=777 NC (12)=KOOI1.(KODI.EQ.O.OR.ION.EQ.I) GO TO 3031. M1=M+21. K1=M+11. DO 302 1=1,N1. J=(I-1)*M11. J1=(I-1)*K11. DO 302 11=1, K1
  193. MSK (I1+J1)=MSK (IH-J) NC (12)=0
  194. IF (M.EQ.777) GO TO 304 WRITE (5,232)
  195. DO 140 IS=NAA, 99 DO 137 M1=1,6
  196. ESMB (IS, M1)=ESMB (IS+1,M1) J=(IS-1)*91. J1=IS*91. DO 138 M1=1,9
  197. F (M1+J)=F (M1+J1) J1=M+KODI+1 J=(IS-1)*J11. DO 139 M1 = 1, J1139 MSK (J+M1)=MSK (J+J1+M1)140 CONTINUE DO 141 IS=1,1001. (NFR (IS).EQ.NAA) NFR (IS)=1000 IF (NFR (IS).GT.NAA) NFR (I)=NFR (IS)-1141 CONTINUE
  198. CONTINUE DO 143 1=1,80 J=81-l1. (Q (J).NE.CN (11)) GO TO 144 143 CONTINUE GO TO 304
  199. WRITE (5,234) DO 146 1=1,69 READ (4,230) Q1. (I.LT.62) GOTO 146 WRITE (5,230) (Q (J), J=1,79)
  200. CONTINUE REWIND 4 GOTO 134
  201. IF (Q (J).EQ.CN (12)) GOTO 134 304 N=N-NA1. WRITE (5,232)1. RETURN 13 FORMAT (80A1)
  202. FORMAT (5X,'Непредусмотренный символ: 71x, 80A1)
  203. FORMAT (5X,'Недопустимая длина хим. формулы :71х, 80А1)
  204. FORMAT (5Х,'Несуществующий символ хим. элемента :'/1х, 80А1)
  205. FORMAT (5Х,'Ошибочный полином :'/1х, 80А1)
  206. FORMAT (' Число компонентов превышает 100')
  207. FORMAT (' Число „резервных“ полиномов больше 100') 29 FORMAT (' Какие вещества нужно исключить ', из рассмотрения (по номерам) ?'/' ')
  208. FORMAT (' Представлен полный перечень возможных веществ.'
  209. Для продолжения нажмите ')30 FORMAT (80А1)
  210. FORMAT (1X, 4(I3, ', ЗА4,2Х))
  211. FORMAT (ЗХ/) 34 FORMAT (23(/))1. END
  212. SUBROUTINE PP100U (KOD, P1, P2,M, K, EL, RP, RS, MSK, FK, NC, G,1. FR, NFR, ESMB)
  213. DOUBLE PRECISION A, RZ, T, STR, X, EX, R0, S, H, CP,
  214. W, U, R1, EK, V, Y, EPS, EPS1, W, VW, SS, VXI
  215. DIMENSION NC (20), G (2), FR (100,9), NFR (100), ESMB (100,3), RP (8),
  216. RS (100), FK (9,100), MSK (1100), A (30,31), T (30), RZ (30),
  217. STR (31), X (100), EX (100), S (100), H (100), CP (100),
  218. W (30), VW (100), ELR (11), NJI (27)
  219. DIMENSION EL (M), SS (9), R1(6), PP (2), RPK (6), STIME (2) CHARACTER НРОЦ32) ================================== -=============================CС
  220. Организация итерационного процесса определения Схарактеристик равновесия, фазового и химического Ссостава термодинамической системы С
  221. Входные параметры: KOD, P1, P2,M, K, EL, MSK, FK, FR, NFR, ESMB С KOD константа, определяющая условия расчета- С Р1, Р2 — значения т.д. параметров- С
  222. М количество химических элементов- С
  223. К количество компонентов равновесия- С
  224. ESMB массив имен индивидуальных веществ С
  225. Выходные параметры: RP, RS, NC, G С
  226. RP значения равновесных параметров- С
  227. RS значения концентраций компонентов равновесия С NC — массив программных констант- С
  228. G заданные термодинамические параметры С ================================================================С
  229. DATA R0/1.98719D0/, R1/1.0D0,1000.D0,0.D0,1D+4,1 D+4,1D+4/,
  230. CONTINUE NC (6)=KODI L1=M+NC (5)+6 L3=M1.=71.=7+M1. ME=M-KODI2 L2=L1+1 KOD=31. (L1.GT.30) GO TO 61 KOD=41. (K.GT.100) GO TO 611.=L5−11.=L4−41.=L4−11. (NC (1).GT.1) GO TO 61. DO 3 1=1,К1. X (l)=1.5
  231. EX (I)=4.48 169 DO 4 l=L4,L14 RZ (I)=-15. DO 5 1=1,65 RZ (I)=R1(I)
  232. DO 7 1=3,4 J=NC (I) RZ (J)=PP (l-2)1. (J.EQ.1) RZ (1)=RZ (1)*9.869 231. (J.EQ.3) RZ (3)=DLOG (0.820 5662D-4*RZ (2)/RZ (3))1. (J.GE.4) RZ (J)=RZ (J)*239.0057 CONTINUE1Т=0 ITER=0
  233. CALL PP100S (K, S, H, CP, RZ, NC, FR, NFR, FK, TMIN, ТМАХ)
  234. EPS=0. DO 15 1=1,1.1 DO 15 J=1,L2
  235. A (l, J)=0. DO 16 J=1,2 l=NC (J+2) A (J, L2)=PP (J)1. (I.GE.4) A (J, L2)=PP (J)*239.005 IF (I.EQ.1) A (J, L2)=PP (J)*9.86 923 IF (I.EQ.3) A (J, L2)=DLOG (0.820 5662D-4*RZ (2)/PP (J))-1.0 IF (I.EQ.3) A (J, 2)=-1.0/RZ (2)16 A (J, I)=1.0
  236. A (3,1)=-DEXP (-RZ (3)) DO 17 1=1, ME17 A (I+3,L2)=ELR (I)1. (NC (2).EQ.1) GO TO 191. DO 18 1=4,61.=l+L318 A (I1,I)=-1.0 IK=L5
  237. XI—1D+6 DO 41 1=1,К J1=(M+1)*I-M+1 J3=J1−1 J2=J3+M1. DO 20 J=J1,J2
  238. NJI (J-J3)=MSK (J) J1=MSK (J1−1)1. (J1.NE.3) GO TO 23 W=H (I)*0.5 032 231 D0/RZ (2) STR (1)=W*EX (I)/RZ (2) STR (2)=-EX (I)
  239. STR (3)=(S (I)*0.5 032 231 D0-W-W+1.0-X (l))*EX (l) W=EX (l)*0.5 032 231 DO DO 21 J=1, L3
  240. A (N, L2)=(RZ (IK)-1.)*Y*EK+A (N, L2) GO TO 2624 A (N, 2)=STR (1)*Y+A (N, 2)
  241. A (N, 3)=STR (2)*Y+A (N, 3) A (N, L2)=A (N, L2)-STR (3)*Y DO 25 J=L4,L625 A (N, J)=STR (J-L7)*Y+A (N, J)26 CONTINUE1. (NC (2).EQ.1) GO TO 33 N=L3+41. (J1.EQ.3) GO TO 27 W=CP (I)*3.0
  242. A (N, 2)=W*E K/RZ (2)+A (N, 2) A (N, IK)=S (I)*EK
  243. A (N, L2)=((RZ (IK)-1.0)*S (I)+W)*EK+A (N, L2) GO TO 29
  244. A (N, L2)=((RZ (IK)-1.0)*U+W*RZ (2))*EK+A (N, L2) GO TO 32
  245. A (N, 2)=CP (I)*EX (I)+STR (1)*U+A (N, 2) A (N, 3)=STR (2)*U+A (N, 3) A (N, L2)=CP (I)*EX (I)*RZ (2)-STR (3)*U+A (N, L2) DO 31 J=L4,L6
  246. A (N, J)=STR (J-L7)*U+A (N, J) U=H (I)-R0*RZ (2)32 CONTINUE
  247. IF (J1.EQ.3) GO TO 38 IF (IT.GE.1) GO TO 36
  248. W=H (I)/RZ (2) A (IK, 2)=W/RZ (2) A (IK, IK)=-.001 A (IK, L2)=W+W-S (I)34 DO 35 J=1,M
  249. A (IK, J+L4−1)=NJI (J) GO TO 40
  250. IF (RZ (IK).LT.-IO.) GO TO 330 A (IK, 2)=H (I)/RZ (2)**2 A (IK, L2)=H (I)*2.0/RZ (2)-S (I) GO TO 34
  251. IF (NC (6).EQ.0) GO TO 40 IF (NJI (L3).EQ.O) GO TO 40 N=ME+4 Y=NJI (L3)
  252. U=DEXP (X (I)+VXI)*Y W=H (l)*0.5 032 231 D0/RZ (2) A (N, 2)=W*U/RZ (2)+A (N, 2)1. A (N, 3)=A (N, 3)-U
  253. A (N, L2)=(W+W-1.0+X (l)-S (l)*0.5 032 231 D0)*U+A (N, L2) W=U*0.5 032 231 DO DO 39 J=L4,L6
  254. A (N, J)=NJI (J-L7−3)*W+A (N, J) IF (X (I).GT.XI) XI=X (I)
  255. IF (J1.EQ.1.0R.J1.EQ.2) IK=IK+141 CONTINUE VXI=-XI N=L11. (IT.EQ.O) N=L6 CALL PP100G (L1,N, T, A) IK=L51. (IT.EQ.O) GO TO 45 IF (EPS1.GT.1D-4) GO TO 43 JV=JV+1 DO 42 1=1,L1
  256. W (I)=W (I)+T (I) DO 420 1=1,К
  257. VW (I)=VW (I)+EX (I) GO TO 4543 JV=01. DO 44 1=1,L144 W (l)=0. DO 440 1=1,К440 VW (l)=0.
  258. DO 51 1=1,К J1=(M+1)*I~M+1 J3=J1−1 J2=J3+M1. DO 46 J=J1,J2
  259. NJI (J-J3)=MSK (J) J1=MSK (J1−1)1. (J1.NE.3) GO TO 50 U=H (l)*0.5 032 231 D0/RZ (2) W=U/RZ (2)*T (2)-T (3)+S (I)*0.503 2231D0-U-U DO 47 J=1,M
  260. W=NJI (J)*0.5 032 231 D0*T (J+L8)+W IF (W.LT.X (I)) GO TO 481. U=1.01. (X (l).LT.-7.0) U=5.0 IF (W-X (I).LT.U) GO TO 49 W=X (I)+U GO TO 49
  261. IF (X (I)-W.LT.3.0) GO TO 49 W=X (l)-3.0
  262. IF (W.GT.10.0) W=10.0 IF (W.LT.-50.0) W=-50.0 U=W-X (I)1. (IT.EQ.1) EX (I)=(W+1.0-X (l))*DEXP (X (l))1.(IT.EQ.O) EX (I)=DEXP (W)1. X (I)=W1. (DABS (W).LT.0.5.0R.W.LT.-20.0) GO TO 51 U=DABS (U/W)1. (U.GT.EPS) EPS=U GO TO 51
  263. X (I)=RZ (IK) EX (I)=DEXP (RZ (IK)) IK=IK+151 CONTINUE1. (NC (3).LE.3.AND.NC (4).LE.3) GOTO 52 IF (T (2).GT.TMAX) T (2)=TMAX+0.001 IF (T (2).LT.TMIN) T (2)=TMIN-0.001
  264. CALL PP100E (EPS, L5, RZ, STR, T, A, NC)
  265. ITER=ITER+1 W=-DLOG (EPS)*2.17 145+5.0 J=W1. (J.LE.1) J=1 IF (J.GT.30) J=31 DO 161 1=1,3061 HPOL (l)=7 DO 162 1=1,31,5
  266. HPOL (l)='+' HPOL (32)=RPK (1) HPOL (J+1)='>' HPOL (J)='=' l=ITER/10 J2=ITER-I*10
  267. WRITE (5,101) HPOL, l, J2 EPS1=EPS1. (ITER.NE.100) GOTO 62 KOD=661 l=NC (3) W=RPK (I) l=NC (4) U=RPK (I)
  268. WRITE (6,100) W, P1, U, P2, KOD GOTO 673
  269. IF (EPS.GT.1D-4) GOTO 68 IF (IT.GT.0) GOTO 631.=11. EPS1 = 1.0 GOTO 68
  270. IF (EPS.GT.1D-5.AND.JV.NE.10) GOTO 68 IF (EPS.LE.1D-5) GOTO 661. (EPS.LT.0.5D-4) GOTO 641. KOD=21. NC (11)=JV
  271. DO 65 1=1,L1 RZ (I)=W (I)/JV5 T (I)=RZ (I)
  272. DO 650 1=1, К 50 EX (I)=VW (I)/JV i6 IF (NC (7).NE.1) GOTO 6701. (KOD.NE.2) KOD=1 70 RP (7)=0. RP (8)=0. DO 671 1=1,К
  273. RP (8)=RP (8)+EX (I) RP (7)=RP (7)+CP (l)*EX (l)/239.005 IF (EX (l).LT.1D-20) EX (l)=0.0 71 RS (I)=EX (I) RP (1)=T (1)/9.86 923 RP (2)=T (2)
  274. RP (3)=0.820 5662D-4*T (2)/DEXP (T (3))1. G (1)=P11. G (2)=P21. (NC (2).EQ.2) GOTO 674 DO 672 1=4,7672 RP (l)=0.0673 M=M-KODI RETURN674 DO 675 1=4,6
  275. RP (l)=T (l)/239.005 GOTO 67 368 KOD=01. (NC (3).EQ.2.0R.NC (4).EQ.2) GOTO 14 GOTO 8
  276. FORMAT (///5X, 1H (, A4, E12.5,A4,E12.5,3H), 3X,'Аварий', * 'ноезавершение, код-, И///)
  277. FORMAT ('+ ', 32А1,1Н (, 211,1Н)) END
  278. SUBROUTINE PP100G (N, NCH, T, A) double precision T, A, S, W, Z DIMENSION T (30), A (30,31)
  279. IF (K.EQ.K2) GOTO 7 DO 6 1=1,N S=A (I, K) A (I, K)=A (I, K2)
  280. A (I, K2)=S K3=M (K) M (K)=M (K2) M (K2)=K3
  281. DO 8 l=KK, N Z=A (I, K)/A (K, K) DO 8 J=K, N1
  282. A (j, J)=A (l, J)-Z*A (K, J) DO 11 J=1,N NN=J-1l=N1-J S=0.1. (N.EQ.I) GOTO 10 DO 9 K=1,NN K1=N1-K9 S=A (K1, N1)*A (I, K1)+S10 K3=M (I)1. (DABS (A (l, l)).LT.1.0D-30) GOTO 11
  283. A (I, N1)=(A (I, N1)-S)/A (I, I)1. T (K3)=A (I, N1)11 CONTINUE RETURN END
  284. SUBROUTINE PP100E (EPS, L5, RZ, STR, T, A, NC) DOUBLE PRECISION TASTR, EPS, RZ, W, U
  285. DIMENSION T (30), A (930), NC (20), RZ (30), STR (31) *************************
  286. ВЫЧИСЛЕНИЕ СТРОКИ КОЭФФИЦИЕНТОВ1. STR (1)=RZ (1)*0.21. STR (2)=200.01. STR (3)=0.51. W=DABS (RZ (3))1. (W.GT.1.0) STR (3)=W*0.51. STR (4)=RZ (4)*0.51. STR (5)=1D+91. STR (6)=1D±952 DO 56 J=1,61. (DABS (T (J)-RZ (J))-STR (J)) 53,53,5453 W=T (J) GOTO 55
  287. IF (T (J).LT.RZ (J)) STR (J)=-STR (J) W=RZ (J)+STR (J)
  288. IF (RZ (I)-T (I).LT.3.0) GOTO 58 W=RZ (l)-3.0
  289. IF (W.GT.5.0) W=5.0 IF (W.LT.-50.) W=-50. IF (T (I).GT.5.0) T (l)=5.0 U=W-RZ (I)1. RZ (I)=W1. (DABS (W).LT.0.5) GOTO 59 U=DABS (U/W)*0.1 IF (U.GT.EPS) EPS=U59 CONTINUE60 RETURN END
  290. SUBROUTINE PP100S (K, S, H, CP, RZ, NC, FR, NFR, FK, TMIN, TMAX) DOUBLE PRECISION T, A, RZ, S, H, CP, RO, U, SS, W, V, Y,1. S1, H1, CP1, A1, A2, CP2
  291. DIMENSION NC (20), FR (100,9), NFR (100>, FK (9,100),
  292. S (100), H (100), CP (100), RZ (100), SS (9) **************************************************
  293. Вычисление термодинамических свойствиндивидуальных веществ (компонентов) **************************************************
  294. W=DABS (SS (1)-RZ (2)) V=DABS (SS (2)-RZ (2))1. (V.LT.W) W=V V=DABS (FR (J, 1)-RZ (2)) Y=DABS (FR (J, 2)-RZ (2)) IF (Y.LT.V) V=Y IF (W.LT.V) GOTO 135 DO 134 J1 = 1,9134 SS (J1)=FR (J, J1)135 CONTINUE
  295. IF (RZ (2).LE.SS (2)) GOTO 11 W=SS (2)*1D-41. NC (7)=1
  296. IF (RZ (2).GE.SS (1)) GOTO 12 W=SS (1)*1D-41. NC (7)=1
  297. S (I)=SS (3)+SS (4)*(1.0+DLOG (W))-SS (5)*W**(-2)+4.0*SS (9)*W+3.0*SS (8))*W+2.0*SS (7))*W H (I)=((((3.0*SS (9)*W+2.0*SS (8))*W+SS (7))*W+
  298. SS (4))*W-SS (6)-2.0*SS (5)A/V)*1D+4 CP (I)=SS (4)+2.0*(((2.0*SS (9)*W+SS (8))*3.0*VV+
  299. DABS (RZ (2)-FR (J, 1)).GT. 1.0) GOTO 123 J2=21. DO 122 J1 = 1,9122 SS (J1)=FR (J, J1) GOTO 12 123 CONTINUE GOTO 13 124 A2=SS (2)-1.01. ((RZ (2)-A2).LT.0.0) GOTO 1251. A1=SS (2)+1.01. GOTO 126 125 A1=SS (1)-1.0 A2=SS (1)+1.0
  300. Y=(RZ (2)-A 1)/(A2-A 1) CP2=DABS (H (I)-H1)*0.51. (CP2.GT.CP1) CP1=CP2
  301. S (I)=(S (I)-S1)*Y+S1 H (I)=(H (I)-H1)*Y+H1 CP (I)=CP1 J2=1
  302. GOTO 13 127 H (I)=H (I)+CP (I)*(RZ (2)-W*1 D+4)
  303. S (l)=S (l)+CP (l)*DLOG (RZ (2)*1D-4/W) 13 CONTINUE RETURN END
  304. SUBROUTINE PP1 OOF (K, RP, RS, NO, EL, M, NC, G, ESMB, NQUE, NQUE1, dt) DIMENSION NC (20), G (2), ESMB (100,3), RP (8), RS (100),
  305. Q (13), QQ (4), QD (4), CN (50), R (10,5), NO (M), EL (M),
  306. NS (94), QREP (80), SY (3), RSP (3), BB1(50) CHARACTER ES (12,3), SYM (12) CHARACTERS2 CC, CC1(50) EQUIVALENCE (SY (1), SYM (1))
  307. DOUBLE PRECISION DQ, A, D, FB, AMu, ALi, A1, Qp3, Qp4,Qpot ============================================================0С
  308. Вывод результатов расчетов СС
  309. Входные параметры: K, RP, RS, NO, EL, M, NC, G, ESMB С К количество компонентов равновесия- С
  310. RP значения равновесных параметров- С
  311. RS значения концентраций компонентов равновесия С N0 — массив номеров химических элементов по Спериодической таблице С
  312. EL массив значений мольного содержания Сэлементов моль/кг. С
  313. М количество химических элементов- С
  314. NC массив программных констант- С
  315. G заданные термодинамические параметры- С
  316. ESMB массив имен индивидуальных веществ- С
  317. Выходные параметры: NQUE, NQUE1 С
  318. NQUE.NQUE1 константы, устанавливающие вывод С результатов на печатающее устройство С ============================================================0
  319. DATA Q/' Р’ЛМПА',' V ТК','4Н У[7куб.',
  320. M/Kr'/StKA'.^K'/r.KJ'/ltKfl'.WKr'/UtKfl1/ DATA CN /'717 273 747 576 777 856', 7,7−7-7+7−7A7B7C7D',
  321. E'.'F'/G'/H'.T/J'/K'/L'.'M'.'N', 07P7Q7R7S7T7U7V7W7X',
  322. Y,'Z','','.','-,'Ж','Д',' 7 7 '/
  323. NQUE1=0 999 NQUE=1 DO 7 1=3,4 DO 1 J=1,31 QQ (J)=Q (3) QQ (4)=G (l-2)1. (NC (I).NE.1) GOTO 2 QQ (2)=Q (1) QQ (3)=Q (2) GOTO 5
  324. IF (NC (I).EQ.2) QQ (3)=Q (4) IF (NC (I).NE.3) GOTO 3 QQ (1)=Q (5) QQ (2)=Q (6) QQ (3)=Q (7)1. GOTO 5
  325. IF (NC (I).NE.4) GOTO 4 QQ (1)=Q (8) QQ (2)=Q (9) QQ (3)=Q (10)
  326. IF (NC (I).LE.4) GOTO 5 QQ (2)=Q (11) QQ (3)=Q (12)1. (NC (I).EQ.6) QQ (2)=Q (13)
  327. IF (I.EQ.4) GOTO 7 DO 6 J=1,46 QD (J)=QQ (J)
  328. CONTINUE DO 8 1=1,10 DO 8 J=1,5
  329. WRITE (5,103) ((R (I, J), J=1,5), EL (I), I=1,J1) IF (NQUE1.EQ.1) WRITE (6,103) ((R (I, J), J=1,5), EL (I), I=1,J1) IF (M.LE.J1) GOTO 10 J1=M1. (M.GT.8) J1 =8
  330. WRITE (5,104) ((R (I, J), J=1,5), EL (I), I=5,J1)1. (NQUE1.EQ.1) WRITE (6,104) ((R (I, J), J=1,5), EL (I), I=5,J1)1. (M.LE.J1) GOTO 10
  331. DO 19 1=1,К SY (1)=ESMB (I, 1) SY (2)=ESMB (I, 2) SY (3)=ESMB (I, 3)
  332. IF (SYM (12).NE.“) GOTO 17 DO 16 J=12,2,-1
  333. SYM (J)=SYM (J-1) SYM (1)='' GOTO 1517 KP=KP+1 DO 18 J=1,12
  334. ES (J, KP)=SYM (J) RSP (KP)=RS (I)1. (KP.LT.3) GOTO 19
  335. WRITE (5,110) ((ES (J, J1), J=1,12), RSP (J 1), J1=1,3) IF (NQUE1.EQ.1) WRITE (6,110) ((ES (J, J1), J=1,12), RSP (J1), J1 = 1,3) WRITE (7,117) ((ES (J, J1), J=1,12), J1 = 1,3) WRITE (8,118) (RSP (J1), J1=1,3) KP=0
  336. CONTINUE REWIND 7 REWIND 8 REWIND 9
  337. READ (9,118) T, V, CP, AM DO 300 j=1,k-1 READ (7,120) CC READ (8,118) BB1.(CC.EQ.' K*C') WRITE (9,118) BB IF (CC.EQ.' H20') WRITE (9,118) BB
  338. CONTINUE REWIND 7 REWIND 8 REWIND 9
  339. READ (9,118) Tfi, Vk, Cp, AM, AMc, H20 OPEN (UNIT=10,FILE='RESUL') WRITE (10,102) QD, QQ WRITE (10,103) ((R (I, J), J=1,5), EL (I), I=1,J1) WRITE (10,106)
  340. WRITE (10,107) (RP (I), I=1,3), RP (5), RP (6), RP (8), RP (4), RP (7) WRITE (10,108)
  341. WRITE (5,103) ((R (I, J), J=1,5), EL (I), I=1,J1- WRITE (5,106)
  342. WRITE (5,107) (RP (I), I=1,3), RP (5), RP (6), RP (8), RP (4), RP (7)1. WRITE (5,108)1. JJ=1
  343. DO 301 J=1,K-1 READ (7,120) CC READ (8,118) BB IF (BB.GT.1E-6) THEN JJ=JJ+1 CC1(JJ)=CC
  344. BB1 (JJ)=100*BB/(AM-H20) BBB=100*BB/(AM-H20) < IF (CC.NE.' H20') WRITE (10,121) CC. BBB ELSE GOTO 301 ENDIF
  345. CONTINUE REWIND 7 REWIND 8
  346. DO 302 J=1,K-1 READ (7,120) CC READ (8,118) BB1.(CC.EQ.' NIO'.OR.CC.EQ.1 N102'.0R.CC.EQ. #' N20'.0R.CC.EQ.' CIO'.OR.CC.EQ. #' H20') THEN WRITE (9,118) BB ELSE GOTO 302 ENDIF302 CONTINUE
  347. Расчет энтальпии (Icm) i-ro участка для ограниченного факела OPEN (UNIT=14,FILE='ENTRYLT') READ (14,200) Bg, Tg, Tb, Pbar, Dg READ (14,119) IM
  348. READ (14,200) aLf, Qn, Gg, Rg, aLo, Rb
  349. OPEN (UNIT=12lFILE=, bxodslg') READ (12,119) I READ (12,122) DX, DQ, a, D, FB OPEN (UNIT=18,FILE-bixodfak') REWIND 1 read (12,119) 11 read (12,122) T, CP1, A1,aad
  350. DATA Cs/4.96/, Cb/0.239/, Cg/0.52/, T1/1 ООО/, Е1/0.35/1. TO=Tfi/10 001. AMM=1000/AM1.(a.GE.1.0.AND.a.LE.100.0) THEN
  351. AMu=(1.664+0.913*TO-0.2*TC)**2)*SQRT (AMM*Tfi)/10**8
  352. Ali=Cp*(0.041+0.175*TO-0.015*TO**2)1. ENDIF1.(a.GE.0.1 .AND.a.LT. 1.0) THEN
  353. Eg=1-exp (-pok) Ef=((1-E1)**2)*(1-Eg)**2 Ec=(1-Eg)*E 1 +(1 -E1)*(1 -Eg)**2 Q1=(T1/100)**4 Q2=(T2/100)**4 Qg=(Tfi/100)**4
  354. Qs=E1*(Eg*Qg-(1-Ec)*Q1)/(1-Ef) R1=Rb/Rg
  355. QpO=Qn*dq/(Rg*(1 +a*aLo*R1))
  356. Qp 1 =(1 +a1 *aLo*R1)*Cp1*t/(1 +a*aLo*R1)1. Qp2=Cp*Tfi
  357. Qp3=4.1868*Cs*fb*Qs/(Gg*(1+a*aLo*R1)) ddT=Tfi-T1
  358. Qp4=ANu*Ali*fb*ddT/(d*Gg*(1+a*aLo*R1)) Alev=QpO+Qp1 Prav=Qp2+Qp3+Qp4dd=abs ((Alev-Prav)/Alev) Qpot=0.5*Qp2+Qp3+Qp4+0.9*Qp0 Qpot=0.5*Qp2+Qp3+Qp4 READ (1,122) P, ai ai1=aad-Qpot ai1=(ai+ai1)/2 rewind 1
  359. READ (9,118) Tfi, Vk, Cp, AM, AMc, H20 READ (9,118) Cn1o, Cn1o2, Cn2o, c1o, h2o Vc=(Vk/AM)*(AM-h2o) print *, AM, Vc
  360. Cco=102.6*C10*(Tfi/Vc)*(0.323*a-0.032) Cno2=Cn1o2+1,533*Cn1o+1.045*Cn2o
  361. Cnox=168.52*Cno2*(Tfi/Vc)*(0.323*a-0.032)
  362. WRITE (10,150) Cnox. Cco Cnox=(46 007*Cno2/Vc)*1.11 *(a-0.1) Cnox1=2.054*Cnox WRITE (10,151) Cnox, Cnox1 W=aLo*Rb/Rg1. CM=(1 +A*W)/1000
  363. WRITE (5,110) ((ES (J, J1), J=1,12), RSP (J1), J1 = 1, KP) IF (NQUE1.EQ.1) WRITE (6,110) ((ES (J, J1), J=1,12), RSP (J1), J1=1,KP) 20 IF (NC (7).EQ.1) WRITE (5,111)1. (NC (11).NE.O) WRITE (5,112) NC (11) WRITE (5,113)1. (NQUE1.EQ.1) WRITE (6,114) RETURN
  364. WRITE (5,115) DO 14 1=1,78 READ (4,101) QREP IF (I.LT.74) GOTO 14
  365. WRITE (5,116) (QREP (J), J=1,79)14 CONTINUE REWIND 4 GOTO 11
  366. FORMAT (/' Нужно записывать полный состав ?')101 FORMAT (80А1)
  367. FORMAT (1X, 76(1 H-)/10X,'Заданные значения 4х, ЗА4,'.=',
  368. G12.5/10Х,'т.д. параметров:, 4X, 3A4,'.=', G12.5)
  369. FORMAT (' Исходный состав рабочего тела ', моль/кг):'/5Х, 4(5А1, F9.5,3X))
  370. FORMAT (5X, 4(5A1,F9.5,3X))
  371. FORMAT (1X, 76(1 H-)/19X,'Параметры равновесного ', состояния занесены в файл')
  372. FORMAT (1Х, 76(1 Н-)//19Х,'Параметры равновесного ', состояния:'/)
  373. FORMAT (6X,'P (Mna)=', G12.5,7X, T (K)=', G12.5,
  374. V (Ky6.M/Kr)=', G12.5/2X,' Н (кДж/кг)=', С12.5,иСкДж/кгКС^.б.гх,' М (моль/кг)=', С12.5/11Х, 8(кДж/кг.К))=', С12.5,Зх,'Ср (кДж/(кг.К))=', 1. G12.5)
  375. FORMAT (19Х,'Содержание компонентов (моль/кг):')
  376. FORMAT (19Х,'Содержание компонентов (%, обьемн.(сух)):')
  377. FORMAT (/19Х,'Содержание компонентов (моль/кг):'/)
  378. FORMAT (1X, 3(12A1,1X, G11.5,2X))
  379. FORMAT (' ** Внимание! ** Выполняласьэкстраполяция термодинамических свойств.')
  380. FORMAT (' ** Внимание! ** Решение результат', осреднения по ', 12,' последним итерациям',»)113 FORMAT (1Х, 76(1 Н-))114 FORMAT (/1Х, 76(1 Н-)/)115 FORMAT (23(/>)116 FORMAT (1Х.79А1)117 FORMAT (12А1)118 FORMAT (e11.5)119 F0RMAT (I3)120 FORMAT (A12)
  381. FORMAT (3(A12,1X, E11.5,1X))122 FORMAT (f10.4)
  382. FORMAT (2x,'Характеристики', i3,'-ro участка факела:1)
  383. FORMAT (1x,'alfaK/(M*4ac*K)=', e10.4)
  384. FORMAT (1x,'Pr=', e10.4,1x,'Re-.elOAIx/Ni^'.elO^)
  385. FORMAT (2x,'Температура факела', 1x, i3,'-слоя Tfi-1(K)=', f6.1)
  386. FORMAT (2x, e10.4,1x, e10.4,1x, f8.4,1x, f8.4,1x, f6.3,2x, i3)
  387. FORMAT (2x,'(Tfi-Tfip)=', f6.1)
  388. FO RMAT (1 x, 'Tf-Tfp=', f 10.3,1 x,'Пересчет для 1ст (кДж/кг)=', П0.3)
  389. FORMAT (1x,'Tf-Tfp=', f10.3,1x,'Конец расчета для i=', i3--ro слоя') 136 FORMAT (1x,'Qn=', f10.4,1 x,'dq=', f10.4,1x,'R1 =', f10.4,1 x,'Rg=', f10.4)
  390. FORMAT (1x,'31=^10.4,1x,'Cp1=', f10.4,1x, Tfi-1=', f10.4)134 format (1x, Cg=', f6.4,1x,'Tg=', f6.2,1x,'Cb-, f6.4)135 format (1x,'Cs=', f6.2,1x,'fb=', f8.4,1x,'Qs=', e10.4)
  391. FORMAT (1x,'ANu=', e10.4,1x,'ALi=', e10.4,1x, lTfi-Tb=', f10.4,1x, #'d=', f10.4)145 format (1x,'Alev1-, e10.4,1x,'Alev2=', e10.4)146 format (1x,'Prav1 =', e10.4,1 x,'Prav2=', e10.4,1 x/Prav3=', e10.4)
  392. FORMAT (1X, Qp1=', e10.4,1x,'Qp2=', e1u.4,1x,'Qp3=', e10.4J1x, #'Qp4=', e10.4)
  393. FORMAT (1 x,'NOx (Mr/M3,15%02)=', e10.4,1 х,'СО (мг/мЗ, 15%02)=', e10.4)
  394. FORMAT (1X,'NOx (Mr/M3,a=1,0)=', e10.4,1x,'NOx (ppm, a=1,0)=', e10.4) 200 FORMAT (F12.4)133 STOP END
Заполнить форму текущей работой

Пора сдавать?