Авторефераты по всем темам  >>  Авторефераты по техническим специальностям  

На правах рукописи

Барочкин Евгений Витальевич

анализ и оптимальный синтез теплообменных

систем со сложной конфигурацией потоков

в энергетических и химических комплексах

05.13.01 Ц Системный анализ, управление и обработка информации

(в химической технологии, нефтехимии и нефтепереработке, биотехнологии)

05.14.14 Ц Тепловые электрические станции, их энергетические

системы и агрегаты

Автореферат

на соискание ученой степени

доктора технических наук

ИВАНОВО

2008

Работа выполнена в ГОУ ВПО Ивановский государственный энергетический университет

Н а у ч н ы й  к о н с у л ь т а н т Ц

доктор технических наук, профессор  Жуков Владимир Павлович

О ф и ц и а л ь н ы е  о п п о н е н т ы:

доктор технических наук, профессор Дорохов Игорь Николаевич

доктор технических наук, профессор Шарапов Владимир Иванович

доктор технических наук, профессор Елин Николай Николаевич

В е д у щ а я  о р г а н и з а ц и я  - Московский государственный университет

инженерной экологии (МГУ ИЭ)

Защита состоится ЕЕЕЕЕЕ. 2008 г. в ЕЕ. часов на заседании Совета по защите докторских и кандидатских диссертаций Д 212.063.05 при ГОУ ВПО Ивановский государственной химико-технологический университет по адресу:

153000, г. Иваново, пр. Энгельса, 7.

С диссертацией можно ознакомиться в библиотеке ГОУ ВПО Ивановский государственный химико-технологический университет.

Автореферат разослан                                2008 г.

Ученый секретарь

совета Д 212.063.05, д.ф.-м.н.                              Зуева Г.А.

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

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

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

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

Актуальность темы работы подтверждается также ее выполнением в рамках ФЦП Интеграция (2.1-А118 Математическое моделирование ресурсосберегающих и экологически безопасных технологий) и международных договоров о научно-техническом сотрудничестве с Ченстоховским политехническим университетом (Польша) и Горным институтом г. Алби (Франция).

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

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

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

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

Научная новизна

1. С позиции системного подхода развиты теоретические основы синтеза математических моделей тепломассообменных аппаратов и их систем со сложной конфигурацией потоков, позволяющих вводить различный уровень декомпозиции моделируемых систем и процессов.
2. Предложена принципиально новая матричная модель для расчета характеристик систем поверхностных теплообменников со сложной конфигурацией потоков теплоносителей, учитывающая возможность фазовых переходов в теплоносителях, в которой модели подсистем (элементов) представлены конечным числом предложенных дифференциальных описаний.
3. Предложенная модель обобщена на случай расчета характеристик многоступенчатых смешивающих тепломассообменных аппаратов струйного и барботажного типов, на основе которой разработана модель деаэрации воды в многоступенчатых деаэраторах струйного и струйно-барботажного типов с учетом внутренних источников массы, обусловленных химическими реакциями.
4. Разработана система кодификации структуры систем тепломассообменных аппаратов, которая однозначно ставит в соответствие данной структуре ее код, а по коду синтезирует разработанную матричную модель системы, позволяющую рассчитать ее характеристики и оценить эффективность. С использованием предложенной системы кодификации разработан генетический (эволюционный) алгоритм структурной оптимизации сложных систем тепломассообменного оборудования.
5. Разработан алгоритмический подход к моделированию переходных процессов в тепломассообменных аппаратах и их системах со сложной конфигурацией потоков, использующий математический аппарат теории цепей Маркова, позволяющий описывать их характеристики в нестационарных режимах работы и позволяющий достаточно просто использовать эти модели в системах автоматического регулирования.
6. Выполнены экспериментальные исследования тепломассообменных процессов при струйной и барботажной деаэрации, восполняющие пробелы эмпирического обеспечения разработанной модели, и их критериальная обработка. На  их основе выполнена параметрическая идентификация и проведена проверка адекватности предложенных моделей процессов тепло- и массопереноса в ряде промышленных аппаратов и расширены возможности эмпирического обеспечения методов их расчета.
7. По ряду выбранных целевых функций сформулированы и решены задачи многопараметрической и многокритериальной оптимизации систем тепломассообменных аппаратов со сложной конфигурацией потоков теплоносителей и показана нетривиальность расчетных рекомендаций по совершенствованию режимно-конструктивных параметров тепломассобменного оборудования.
8. Поставлена и решена с использованием алгоритма Беллмана нелинейная задача об оптимальном по расходу топлива распределении выработки тепловой и электрической энергии на совокупности турбоагрегатов с различными характеристиками при меняющейся потребности в тепловой и электрической энергии с учетом сложных технологических ограничений на область допустимых значений управляющих параметров.

Практическая ценность

1. Проведенные экспериментальные исследования позволили предложить расчетные зависимости для коэффициентов тепло- и массопередачи при струйной и барботажной деаэрации, которые положены в основу эмпирического обеспечения предложенного метода расчета для ряда промышленных аппаратов.
2. На основе разработанной модели и проведенной идентификации предложен инженерный метод расчета систем тепломассообменных аппаратов поверхностного и смешивающего типов со сложной конфигурацией потоков, базирующийся на единой алгоритмической основе. Разработанный метод расчёта даёт более широкую и достоверную информацию о характеристиках сложных технических систем, что позволяет вырабатывать надежные рекомендации по повышению эффективности их работы и снижению риска экологического ущерба от теплового загрязнения окружающей среды.
3. На основе предложенных методов расчета разработаны алгоритмы и пакеты прикладных программ для расчета тепломассообменных систем со сложной конфигурацией потоков. Разработанные подходы матричной формализации положены в основу расчета энергетических характеристик оборудования и использованы при создании программного комплекса УТЭС ЭКСПЕРТФ, внедрение которого только на Владимирской ТЭЦ позволило сэкономить 14 тыс. Гкал/год тепловой энергии за счет выбора оптимального состава и режима работы оборудования.
4. Разработанные математические модели, методы расчета, их программная реализация, решения на их основе проектных и конструкторских задач внедрены в Ченстоховском политехническом университете (Польша), на Ковровском заводе котельного и сушильного оборудования СОЮЗ, на Костромской ГРЭС (ОАО ОГК-3), Владимирской и Дзержинской ТЭЦ (ОАО ТГК-6), ОАО Северсталь, МУП Теплоэнергия г. Череповца. Автор защищает

1. Развитие теоретических основ и принципов системного подхода к синтезу математических моделей массообменных аппаратов и их систем со сложной конфигурацией потоков, позволяющих вводить различный уровень декомпозиции моделируемых подсистем и процессов.
2. Новые матричные модели для расчета характеристик систем поверхностных и смешивающих теплообменников струйного и барботажного типа, учитывающие возможности фазовых переходов в теплоносителях и универсальный алгоритм компьютерного синтеза этих моделей из моделей отдельных подсистем и процессов.
3. Матричную модель деаэрации воды в многоступенчатых деаэраторах струйного и струйно-барботажного типов с учетом внутренних источников массы, обусловленных химическими реакциями.
4. Единый алгоритмический подход к моделированию переходных процессов в тепломассообменных аппаратах и их системах, основанный на математическом аппарате теории цепей Маркова и позволяющий описывать их характеристики в нестационарных режимах работы и инкорпорировать эти модели в системы автоматического регулирования.
5. Результаты экспериментальных исследований и выполненной на их основе параметрической идентификации предложенных моделей для ряда промышленных тепломассообменных установок.
6. Систему кодификации, которая по заданной структуре системы однозначно формирует ее код, а по заданному коду осуществляет компьютерный синтез математической модели системы, и которая позволяет эффективно осуществлять выбор оптимальных структур с использованием эволюционных (генетических) алгоритмов.
7. Результаты решения расчетных и оптимизационных задач по ряду выбранных целевых функций режимно-структурных параметров тепломассобменного оборудования, в том числе, решения с использованием алгоритма Беллмана нелинейной задачи об оптимальном по расходу топлива распределении выработки тепловой и электрической энергии на совокупности турбоагрегатов с различными характеристиками при меняющейся потребности в тепловой и электрической энергии с учетом сложных технологических ограничений на область допустимых значений управляемых параметров.

Апробация работы. Основные положения диссертации доложены и обсуждены на следующих конференциях: VI Miedzynarodowa konferencja naukowa ДNowe technologie i osiagniecie w Metalurgii i inzynierii MaterialowejФ, Politechnika Czestochowska (Польша 2005), 34th International Conferense of Slovak of Chemical Engineering (Словакия, 2007), Международной научной конференции Математические методы в технике и технологиях (ММТТ-2001, 2003, 2004, 2005, 2006, 2007), Международной научно-технической конференции Энерго- ресурсосберегающие технологии и оборудование, экологически безопасные производства (Иваново 2004, 2007), Международной научно-технической конференции Состояние и перспективы развития электротехнологии (Бенардосовские чтения) (Иваново 2003, 2006, 2007), Всесоюзной научно-технической конференции. "Химические реакторы "(Чимкент, 1983). Всесоюзной научно-технической конференции. "Физико-химические проблемы материаловедения и новые технологии" (Белгород, 1991).

Публикации. Основное содержание диссертации опубликовано в 51 печатной работе, в том числе в 15 изданиях, предусмотренных перечнем ВАК.

Объем и структура диссертации. Диссертационная работа состоит из введения, пяти глав, основных выводов, списка использованных источников (275 наименований работ) и приложений. Работа содержит 308 страниц, в том числе 294 страниц основного текста, 2 приложения на 14 страницах.

СОДЕРЖАНИЕ РАБОТЫ

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

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

Основы математического моделирования собственно процессов тепломассопереноса заложены в  основополагающих  работах  А.В. Лыкова, С.С Кутателадзе, А.И. Леонтьева и ряда других авторов. В дальнейшем они развивались В.И. Коноваловым, С.В. Федосовым, П.Г. Романковым и другими исследователями. Однако, объектом исследования в этих работах были или собственно процесс, или процесс в отдельном аппарате, а не в сложной их совокупности, включающей, например, и обратные связи.

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

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

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

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

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

Очевидно, что с ростом числа ячеек резко возрастает и число потоков теплоты и (или) массы между ними, что, как правило, требует большого объема экспериментального обеспечения расчета этих потоков. Вместе с тем, переход от локального описания к более интегральному (укрупнение ячеек и уменьшение их числа) всегда возможен и единственен; обратный же переход без привлечения дополнительной информации (или гипотез) невозможен.

Показано, что работа оборудования в переходных режимах часто связана с выходом параметров процесса из допустимых диапазонов, кроме этого, моделирование переходного процесса необходимо при разработке многофункциональных тренажеров и систем автоматического регулирования. Анализу нестационарных и переходных процессов в химических технологиях и энергетике, в том числе и с использованием математического аппарата теории цепей Маркова, посвящены работы И.О. Протодьяконова, В.И. Шарапова, В.Е. Мизонова, А. Тамира, А. Бертье, В.М. Рущинского и др. Однако, вопросы описания переходных процессов в системах теплообменных аппаратов с учетом массообмена, например, при деаэрации остались малоизученными.

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

Показана актуальность в изменившихся условиях функционирования энергетики задачи оптимального распределения нагрузки между энергетическими установками в энергетических комплексах.

В заключении главы формулируются задачи исследования.

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

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

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

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

Модель 1. Фазовый переход теплоносителей отсутствует (ПЦП, ПЦВ, ВЦП, ВЦВ).

Модель 2. Фазовый переход только холодного теплоносителя (НЦВ, НЦП).

Модель 3. Фазовый переход только горячего теплоносителя (ВЦН, ПЦН).

Модель 4. Фазовый переход холодного и горячего теплоносителей (НЦН).

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

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

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

, (1)

где BЦматрица процесса в ступени, ТЦвектор параметров ступени, То - вектор свободных членов. Вид перечисленных матриц для четырех рассматриваемых моделей показан в табл.1.

Для установки из n ступеней с произвольной структурой связей между ними предложена методика синтеза модели системы из моделей подсистем. Для задания направления движения потоков горячего (G), холодного (X), совместно горячего и холодного (GX) теплоносителей предложены матрицы, автоматически приравнивающие значения параметра выходного потока ячейки значению параметра входного потока той ячейки, куда этот поток направляется:

, , . (2)

Система из четырех алгебраических уравнений (1) записывается для n ступеней установки, а полученная система из 4n линейных уравнений приводится к матричному виду

или , (3)

где АЦблочная матрица, определяющая состав и структуру системы, каждый блок которой Aij размера 4×4 соответствует одной из матриц процесса табл. 1 или матрице связи (2). Вектор искомых параметров определяется из решения матричного уравнения (3) в виде

. (4)

Таблица 1

Фазовое состояние теплоносителей	Графики изменения температур горячего и холодного теплоносителей вдоль поверхности нагрева	Математическая модель	Решение модели для теплообменника (или ступени) в матричном виде B*T=To
				T	To
мод.1В-В, П-П П-В В-П
мод.2Н-П Н-В
мод.3П-Н В-Н
мод.4Н-Н

Матрица A для расчета произвольной установки из n ступеней составляется в соответствии со следующими правилами:

1. Каждому теплообменнику (ступени, подсистеме) установки присваивается номер .
2. Теплообменнику c номером i соответствует в блочной матрице A i-я cтрока. В каждой строке располагается не более трех ненулевых блоков.

3. На главной диагонали матрицы размещаются в зависимости от фазового состояния теплоносителей или реализуемой модели матрицы процессов из табл. 1: Aii =Bi.
4. Если горячий теплоноситель ступени подается в j-ый элемент, то Aij =G.
5. Если холодный теплоноситель ступени подается в j-ый элемент, то Aij =Х.
6. Если горячий и холодный теплоносители из i-ой ступени подаются в j-ый элемент, то Aij =GХ.
7. Столбец неизвестных T содержит n блоков или 4n элементов, соответствующих параметрам теплоносителя внутри установки.
8. Ненулевые элементы столбца свободных членов TO равны параметрам известных внешних потоков. При подаче внешнего потока в i-ый элемент в i-ой строке матрицы TO содержится ненулевой элемент.

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

Рис.а3. Схема потоков в трехступенчатом теплообменнике, вид матрицы А, определяющей структуру и состав системы, и пример ее кода

Нумерация элементов выполнена по ходу движения горячего теплоносителя, холодный теплоноситель последовательно проходит ступени 3Ц1Ц2.

Матричная формализация расчета позволила разработать систему кодификации тепломассообменных систем, которая по заданной структуре однозначно формирует код, а по заданному коду восстанавливает структуру системы. Каждой ступени соответствуют четыре позиции кода, разделенные точками и записанные на рис. 3 в круглых скобках. Первая позиция показывает номер ступени, в которую направляется горячий теплоноситель. Вторая позиция соответствует типу канала, в который подается горячий теплоноситель: 1-горячий, 2-холодный. Значения третьей и четвертой позиций выбираются аналогично для холодного теплоносителя. Предложенная система кодификации позволяет определить число возможных вариантов структур системы, которое для системы из N элементов составляет (2N)2N. Для рассматриваемого трехступенчатого теплообменника число схем равно 66.

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

Для иллюстрации возможностей и порядка построения модели сложной системы из моделей подсистем рассмотрен пример расчета регенеративного подогревателя ПН-400-26-2 (рис. 4), используемого для подогрева основного конденсата паровой турбины К-300-240. Подогреватель имеет три хода греющего пара и четыре хода воды.

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

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

Нормальный режим эксплуатации установки соответствует практически полной конденсации пара (кривая 2). На рис. 5.б приведена зависимость расхода воды от расхода пара при условии полной его конденсации на выходе из установки (x1=0). Сопоставление расчетных и экспериментальных результатов для промышленных установок позволяет диагностировать состояние теплообменников. В частности, подбирая значение коэффициента теплопередачи для совпадения расчета с экспериментом и сравнивая его с нормативными данными, можно диагностировать состояние поверхностей нагрева. Разработанный метод позволяет также учитывать подачу в аппарат дополнительных теплоносителей или изменение поверхностей теплообмена.

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

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

Из баланса энергии и массы в элементах с поверхностью обмена dF (рис.6) записана система дифференциальных уравнений, описывающая изменение температурного напора τ и массового расхода теплоносителей G вдоль определяющей координаты F

(5)

где kЦкоэффициент теплопередачи, сЦудельная теплоемкость, GЦрасход, tЦтемпература теплоносителя, rЦудельная теплота парообразования, индекс л1 относится к горячему, л2Цк холодному теплоносителю, n - к состоянию насыщения. Решение (5) при начальных условиях получено в виде

  (6)

Результаты расчета смешивающего теплообменника (F=200 м2, G10=5 кг/с, G20=200 кг/с,  c2 =  4 кДж/(кгК), k  =  5000 Вт/(м2K), τ0=15 оС, r =2258 кДж/кг) приведены на рис. 7 в виде зависимости температурного напора и расхода воды от площади поверхности нагрева.

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

  (7)

Аналитическое решение (7) при начальных условиях может быть записано в матричном виде

или

,  (8)

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

Сопоставление результатов расчета по упрощенной (7) и точной (5) модели на рис. 7 показывает их практическое совпадение.

Матричное представление (8) позволяет интегрировать модель в расчет многоступенчатых установок, матричная формализация которых рассмотрена ранее. Приводится пример расчета тарельчатого смешивающего подогревателя, показанного на рис. 8. В подогревателе выделяется 4 ступени, нумерация которых выполнена по ходу воды. В соответствии с предложенными правилами блочная матрица А в уравнении (3) для системы из четырех смешивающих подогревателей имеет вид

.  (9)

Решение системы (3) с учетом (8), (9) позволяет определить температуру холодного и степень сухости горячего теплоносителя в каждой ступени.

Системы смешивающих тепломассообменных аппаратов - деаэраторов

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

• теплообмен между паром и водой,
• массообмен между паром и водой,
• массообмен между газом, растворенным в воде, и газом паровой фазы.

Температура пара в состоянии насыщения остается постоянной, поэтому температурный напор τ меняется только вследствие изменения температуры воды t2. Теплообмен между паром и водой осуществляется через поверхность раздела фаз - поверхность струй. При постоянной температуре греющего пара теплообмен мало зависит от направления движения теплоносителей, поэтому при составлении расчетной схемы рассматривается прямоточный теплообменник (рис. 10). В качестве определяющей координаты выбирается площадь поверхности теплообмена F. На элементарном участке dF вода нагревается за счет теплопередачи от пара и смешения с его конденсатом. За счет разности концентраций происходит удаление газа из воды.

Из баланса теплоты и массы (рис. 10) получена система дифференциальных уравнений, описывающая изменение температурного напора τ, массового расхода теплоносителей G и концентрации газа в воде (cg2) и паре (cg1) вдоль определяющей координаты F

(10)

где Цтемпературный напор, с - удельная  теплоемкость, сg - концентрация газа, kmЦкоэффициент массопередачи, kgЦкоэффициент, определяющий связь между концентрацией газа в воде и равновесной концентрацией газа в паровой фазе.

Первое уравнение системы (10) определяет изменение температурного напора между паром и водой, второе - уменьшение расхода пара за счет его конденсации, третье - увеличение расхода воды за счет смешения с конденсатом пара. Движущей силой дегазации является разность между равновесной концентрацией газа в паровой фазе и концентрацией газа в паровой фазе. Равновесная концентрация газа определяется в соответствии с законом Генри. Система из пяти дифференциальных уравнений (10) описывает изменение пяти неизвестных функций (τ,G1,G2,cg1,cg2) вдоль координаты F. Процессы тепло- и массообмена между паром и водой сопряжены с большими массо- и энергопотоками и от процесса деаэрации практически не зависят, что позволяет решать систему (10) в два этапа. На первом этапе решаются первые три уравнения системы (10) аналогично решению системы (4). Найденные из решения (3) распределения температур и расходов теплоносителей по ступеням позволяют перейти ко второму этапу моделирования - расчету деаэрации. На втором этапе решаются уравнения газообмена (четвертое и пятое уравнение системы (10)). С учетом допущения о постоянном расходе пара и воды внутри ступени и начальных условий решение записывается в виде

,  (11)

где

,

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

Описание ступени деаэрации в матричном виде (11) позволяет по аналогии с теплообменом для нахождения значений концентрации газа в произвольной точке установки воспользоваться правилами матричной формализации и построить матричную модель в виде

,  (12)

где А - блочная матрица размером n×n блоков, определяющая состав и структуру установки, С, СОЦматрицы-столбцы неизвестных концентраций и свободных членов.

Решение системы (8) позволяет определить значение температуры и расхода, а решение (12)Цконцентрацию газа в произвольной точке многоступенчатой установки.

Уравнения (3), (12) представляют замкнутое описание струйной деаэрации при произвольной структуре потоков. Разработанная модель деаэрации построена на допущении прямоточного движения теплоносителей. Для исследования влияния выбора расчетной схемы системы на результаты расчета процесса тепломассопереноса ступень - зона между двумя тарелками - разбивается различными способами на 1, 3 и 9 подсистем (рис. 9). Результаты расчета, приведенные в табл. 2, показывают, что число элементов и их расположение не влияют на результаты расчета процесса тепломассообмена для воды и пара и влияют незначительно на процесс деаэрации. Рассмотренные модели пригодны в тех случаях, когда фазовые состояния в каждой ступени теплообменного аппарата остаются постоянными и известны заранее. В этом случае для описания каждого элемента (ступени) выбирается одна из моделей табл.1.

Таблица 2

Номер модели (см. рис. 9)	Число элементов	G1, кг/с	G2, кг/с	t2, оC	cg1, мкг/дм3	cg2, мкг/ дм3
1	1	3,57	201,43	89,03	271,5	15,2
2	3	3,57	201,43	89,03	278,0	15,1
3	3	3,57	201,43	89,03	279,9	15,0
4	9	3,57	201,43	89,03	283,0	15,0

Если границы фазовых переходов теплоносителей заранее неизвестны, то для расчета ступени предлагается следующий алгоритм:

1. Поверхность теплообмена разбивается на участки площадью ΔF каждый.
2. На каждом шаге определяется состояние теплоносителей и из табл. 1 выбирается соответствующая расчетная модель. По заданным входным параметрам определяются выходные параметры ступени. Значение площади поверхности теплообмена увеличивается на величину шага.
3. Расчет согласно п. 2 выполняется до достижения конечного значения площади поверхности нагрева.

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

Обобщенный метод матричной формализации. Рассмотренные подходы применимы для описания технологических систем при подаче всего потока из одной ступени на вход другой и при отсутствии рециркуляции или байпасирования потоков теплоносителей. Для устранения указанных ограничений разработан новый обобщенный метод матричной формализации. В основу математического описания положены также уравнения материального и теплового баланса. Пусть установка включает n элементов (ступеней). Ступень теплообмена представляется четырехполюсником с входными и выходными потоками для холодного и горячего теплоносителей. Матрица-столбец (вектор признаков) X составляется из аддитивных параметров теплоносителей, значения которых можно складывать при смешении потоков. Такими параметрами являются, например, теплота или масса.

Наиболее общий случай соединения ступеней предполагает возможность подачи на вход в i-ый элемент потоков из всех остальных элементов. На рис. 11 схематично показано формирование потока на входе в i-ый элемент установки. На вход в первый элемент подаются потоки холодного и горячего теплоносителя, которые характеризуются набором признаков [X0]1, индекс л0 внутри квадратных скобок указывает на входные параметры, индекс л1 за квадратными скобками относится к номеру ступени. Процесс теплообмена, массообмена или их совместное протекание описывается матрицей процесса B, произведение которой на вектор входных параметров позволяет определить выходные параметры ступени [X]1=B1[X0]1. Для указания направления движения потока после ступени сформируется матрица коммутации Ki1, элементы которой показывают доли потока из первой ступени в i-ю. Матричное произведение Ki1B1[X0]1 определяет параметры потока, подаваемого из первого элемента в i-й. Очевидно, что на вход в i-ю ступень могут подаваться потоки из остальных элементов схемы и внешние для установки потоки. Входной вектор признаков для i-го элемента определится как сумма аддитивных характеристик смешиваемых на входе в него потоков

(13)

где индекс вх указывает на внешний поток, подаваемый на вход ступени. Уравнения аналогичные (13) записываются для каждой из n ступеней установки. Полученная система уравнений после преобразований представляется в матричном виде

,  (14)

где IЦединичная матрица.

При известных матрицах B и K решение системы (14) позволяет определить значения параметров теплоносителей в любой точке установки. Система (14) для целого ряда практически важных случаев является системой линейных уравнений. Рассмотрим порядок составления системы уравнений (14) и ее решение на ряде конкретных примеров применительно к широко используемым в промышленности пароводяным и водоводяным теплообменникам.

Пример 1. Каскад поверхностных теплообменников без учета теплообмена (теплоносители с одинаковыми температурами). Вектор признаков составляется из массовых расходов горячего G1 и холодного G2 теплоносителей . Теплоносители в поверхностном подогревателе разделены стенкой, поэтому их расходы не меняются, а матрица процесса равна единичной матрице B=I. Матрица коммутации состоит из четырех элементов , каждый из которых показывает долю потока из j-го в i-ый элемент: -горячего теплоносителя в горячий, -холодного в горячий, -горячего в холодный, - холодного в холодный. Наиболее простой и практически значимый вид матрицы коммутации представлен в табл. 3. Задавая дробные значения элементов матриц коммутаций на отрезке [0,1], можно описывать деление теплоносителей на произвольное число потоков.

Таблица 3

Вид матрицы коммутации при подаче из j-го в i-ый элемент
горячего теплоносителя	холодного теплоносителя	горячего и холодного теплоносителя	теплоносители не подаются

Решение системы (14) при известных матрицах и позволяет найти расходы теплоносителей в любой точке установки. Следует отметить, что предложенный подход справедлив и при наличии внутренних рециркуляций.

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

, .

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

  (15)

где Цдиагональная матрица расходов теплоносителей через ступень, Цвектор-столбец из произведений температуры на расход внешних потоков. Из решения (15) находится температура теплоносителя в любой точке установки.

Пример 3. Каскад смешивающих теплообменников. В смешивающих подогревателях процессы тепло- и массообмена происходят одновременно. Температура горячего теплоносителя считается постоянной и равной температуре насыщения. В качестве параметров, определяющих процесс, выбраны степень сухости пара и температура холодного теплоносителя . Уравнение связи входных и выходных параметров и вспомогательные матрицы А, B записываются виде

С учетом сделанных замечаний система (14) приводится к виду

.  (16)

Расчет смешивающего подогревателя проводится в два этапа: определяются аналогично примеру 1 начальные расходы теплоносителей, затем согласно (16) находятся температуры и уточненные расходы теплоносителей во всех точках установки.

Пример 4. Расчет многоступенчатых деаэраторов. В деаэраторах наряду с подогревом воды происходит удаление из неё растворенных газов. Процесс деаэрации связан с незначительными массо- и энергопотоками, что позволяет рассматривать процессы теплообмена и деаэрации последовательно. На первом этапе решается задача определения температур и расходов теплоносителей аналогично примеру 3. Движущей силой деаэрации является разность концентраций газа. В качестве вектора параметров на втором этапе решения задачи выбирается концентрация газа в паровой фазе и воде . Матрица процесса деаэрации, полученная ранее, записывается в виде

, ,

где kmЦкоэффициент массопередачи, kgЦкоэффициент, определяющий связь между концентрацией газа в воде и равновесной концентрацией газа в паровой фазе.

Составление материальных балансов по газу приводит систему (14) к виду

(17)

где матрица определяет изменение концентрации газа за счет смешения конденсата пара с водой, - изменение степени сухости пара внутри ступени. Решение системы (17) позволяет определить концентрации газа внутри установки.

Рассматривается пример расчета деаэратора  (G10 = 3 кг/с, G20 = 100 кг/с, c2=4 кДж/кг/К, k = 5000 Вт/м2*K, τ0 = 15 оС, r = 2258 кДж/кг, сg10 = 0 мкг/дм3, сg20=21 мкг/дм3), имеющего четыре струйных отсека в деаэрационной колонке и один барботажный отсек в водяном объеме деаэраторного бака. Технологическая и расчетная схемы этого деаэратора приведены на рис. 12.

ю

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

Алгоритм расчета включает следующие операции:

• расчет согласно (16) температуры воды и расхода пара и воды по ступеням,
• расчет согласно (17) концентрации газа в паровой и водяной фазах.

Изменение температуры воды и расхода пара вдоль ступеней теплообмена показано на рис. 13. Изменение концентрации газа в паровой и водяной фазах внутри деаэратора приведены на рис. 14. Экстремальное значение концентрации газа в воде перед второй тарелкой (рис. 14,а) свидетельствует о неэффективной работе первого струйного отсека.

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

Построение такой модели рассмотрено на примере прямоточного кожухотрубного теплообменника поверхностного типа (рис. 15). Теплообменник разбивается по длине на n элементов. Рассматривается процесс теплообмена без учета фазового перехода в теплоносителях. Вводится вектор состояния системы, составленный из температур горячего (t1) и холодного (t2 ) теплоносителей , где первый индекс характеризует теплоноситель: л1 - горячий, л2 - холодный, второй индекс соответствует номеру ячейки. Параметры среды внутри каждой ячейки равномерно распределены по ее объему. Состояние процесса рассматривается через конечные, малые промежутки времени Δ, в течение которых перенос аддитивного свойства возможен только в соседние ячейки, но не далее. Текущее дискретное время определяется как , где j=1,2Е-номер перехода или дискретный аналог времени. Параметры (j+1) и j-го состояния связаны матричным равенством

, (18)

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

Для построения переходной матрицы, соответствующей матрице переходных вероятностей в теории цепей Маркова, рассматриваются тепловые балансы для ячейки горячего и холодного теплоносителей (рис. 15.). Изменение энергии в ячейке происходит за счет внутренних и внешних перетоков теплоносителей и теплопередачи от горячего теплоносителя холодному. Тепловая энергия среды в ячейке представляется как произведение массы G, теплоемкости c и температуры теплоносителя в виде , где верхний индекс показывает номер шага по времени. Считается, что режим движения теплоносителя соответствует идеальному вытеснению со скоростью , где QЦобъемный расход теплоносителя, SЦплощадь сечения канала. Доля теплоносителя, покидающего ячейку за время Δ, или вероятность перехода в соседнюю ячейку находится как отношение перемещения среды к длине ячейки:,

где w - скорость теплоносителя,  Δ - шаг по  времени, Δx - длина ячейки, ΔV - ее объем, B - массовая производительность, G - масса теплоносителя в ячейке. Перейдя от вероятности к долям, будем называть матрицу Р просто переходной матрицей.

Считается, что размеры канала, расход и скорость движения не меняются вдоль аппарата, но различаются для горячего и холодного теплоносителей. В следующую i+1 ячейку переходит количество теплоты, определяемое произведением . Теплота, приходящая с теплоносителем из предыдущей i-1 ячейки, равна . За счет внешних потоков энергия изменяется на величину , за счет теплопередачи уменьшение тепловой энергии горячего теплоносителя составит . Баланс энергии для ячейки горячего теплоносителя записывается в виде

.

Считая заполнение ячеек по длине реактора одинаковым, а теплоемкость постоянной, уравнение преобразуется к виду

, (19)

где .

Аналогичные рассуждения для холодного теплоносителя позволяют записать

,  (20)

где .

Выражения (19) и (20), записанные для всех ячеек аппарата, формируют искомую переходную блочную матрицу P

;;;;.

Матрица Pа учитывает переход тепловой энергии от горячего теплоносителя к холодному. Матрица Pt описывает переходы теплоносителей из ячейки в ячейку. Матрица ТР характеризует внешние для аппарата потоки, которые могут подаваться в любые ячейки. Матричное уравнение (18) позволяет определить распределение температуры теплоносителей вдоль канала в произвольный момент времени.

Пример расчета переходного процесса для прямоточного кожухотрубного подогревателя проиллюстрирован на рис. 16. Эскиз подогревателя и направления движения теплоносителей показаны на рис. 16.а. Теплообменник по длине условно разделен на 5 участков, холодному и горячему теплоносителю на каждом элементе длины соответствует своя ячейка. Вектор состояния системы, составленный из температур теплоносителей, приведен на рис. 16.в. Внешние для аппарата потоки холодного и горячего теплоносителей подаются в первую ячейку.

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

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

где SijЦпараметр теплоносителя в ячейке, индекс i=1 соответствует массе горячего, i=2Цмассе холодного, i=3Цэнергии горячего, i=4Цэнергии холодного теплоносителя, индекс j=1,2,..n показывает номер ячейки.

Вектор состояния Sg- описывает содержание газа в паровой и водяной фазе

где SgijЦмасса газа в ячейке, индекс i=1,2 характеризует соответственно горячий и холодный теплоноситель, индекс j=1,2,..n показывает номер ячейки.

На каждом временном шаге k моделирование процесса осуществляется в два этапа: на первом этапе описывается преобразование по времени вектора состояния S,

,

где ,

где p1, p2Цвероятность перехода в соседнюю ячейку горячего и холодного теплоносителя соответственно. Вектор питания системы Q характеризует внешние массо- и энергопотоки теплоносителей, которые могут подаваться в любые ячейки.

Матрица свободных членов С размера 4n1 содержит для каждой ячейки следующие четыре элемента

,

где t- температура теплоносителя,  индекс  n относится к состоянию насыщения.

Найденные распределения массы и энергии теплоносителей по ячейкам позволяют перейти ко второму этапу решения задачи - расчету деаэрации

,

где переходная матрица

 

где Qg Цвектор питания, km - коэффициент массопередачи, kg - коэффициент, определяющий связь между концентрацией газа в воде и равновесной концентрацией газа в паровой фазе.

В четвертой главе изложено описание условий и методики проведения экспериментальных исследований. Целью этих исследований является получение опытных данных для идентификации, верификации модели и экспериментального обеспечения метода расчета процесса тепломассообмена в многоступенчатых деаэраторных установках. Эксперименты проводились на установке подпитки теплосети с деаэратором струйного типа филиала ОАО ОГК-3 Костромская ГРЭС (КГРЭС) и на установке подготовки технологической воды с деаэратором струйно-барботажного типа ОАО Северсталь. На рис. 17 приведена фотография пробоотборных устройств и схема экспериментальной установки промышленного деаэратора ДСА-300 струйно-барботажного типа с указанием замеряемых параметров. Экспериментальные данные условно разделяются на теплотехнические и химические. К первой группе относятся расход и температура теплоносителей, ко второй - концентрации растворенных в воде газов, ее щелочность и др.

	Рис.18.Идентификация расчетной модели по полученным экспериментальным данным: а)Цзависимость рассчитанной (линии) и замеренной (точки) температуры воды от номера ступени; b)Цзависимость рассчитанной (линии) и замеренной (точки) концентрации газа в воде от номера ступени; с), d) сопоставление рассчитанных и найденных из эксперимента коэффициентов теплопередачи и массопередачи соответственно

Расчет деаэратора (глава 2) выполняется также в два этапа: сначала определяются расходы и температуры теплоносителей, затем содержание газов в паровой и водяной фазах. В соответствии с вышесказанным идентификация модели (рис. 18) проводится тоже в два этапа: на первом этапе при выполнении теплового расчета, используя теплотехнические экспериментальные данные, определяются коэффициенты теплопередачи.

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

Для разработанной модели выполнена оценка показателя точности. Характеристика точности модели - величина относительной погрешности расчета коэффициентов тепло- и массопередачи относительно их значений, определенных по результатам идентификации математической модели. Корреляции между относительной погрешностью и уровнем коэффициентов тепло- и массопередачи не выявлено. В  среднем относительная  погрешность расчета  коэффициента  теплопередачи составила 2,4 %, коэффициента  массопередачи 8,5 %.

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

1. Разработка алгоритма и программная реализация методов расчета процесса тепломассопереноса в установках со сложной конфигурацией потоков.
2. Использование разработанных алгоритмов синтеза сложных структур теплообменных аппаратов для решения следующих практических задач: расчет деаэратора, расчет барабанного котла, реконструкция водогрейного котла в паровой котел, расчет систем пластинчатых теплообменников.
3. С использованием принципа оптимальности по Парето решены задачи системного анализа работы водогрейного котла (рис. 19), струйно-барботажного деаэратора, системы пластинчатых теплообменников.
4. Разработка методики построения энергетических характеристик оборудования и решения на их основе с использованием принципа оптимальности по Беллману задач оптимального распределения нагрузок между турбогенераторами ТЭЦ (рис. 20). В качестве критерия оптимизации при распределении нагрузки используется величина удельного расхода топлива на выработку электрической энергии, которая является величиной обратной коэффициенту полезного действия.

Структура программного комплекса ТЭС-ЭКСПЕРТ с указанием информационных

Рис. 21. Состав информационных модулей и структура информационных потоков в программном комплексе УТЭС ЭКСПЕРТФ: БД - база данных; ТБ - модуль сведения тепловых балансов;

ПВБ - модуль сведения пароводяных балансов; ЭХ - расчет и компьютерное представление энергетических характеристик оборудования; ТЭП - технико-экономические показатели

потоков между основными информационными и  расчетными  модулями  приведена на рис. 21.

Комплекс позволяет оптимизировать состав и режим работы оборудования по выбранной целевой функции, как при оперативном ведении режима, так и при перспективном планировании работы энергетического комплекса. Внедрение комплекса только на одном энергетическом предприятии (Владимирском филиале ОАО ТГК-6) позволило снизить суммарный расход тепловой энергии на 14 тыс. Гкал/год, что эквивалентно экономии более 5 миллионов рублей в год.

Разработанные методики расчета, их программное обеспечение, а также результаты решения практических задач с их использованием реализованы на Костромской ГРЭС ОАО ОГК-3, Владимирской и Дзержинской ТЭЦ ОАО ТГК-6, ОАО Северсталь, МУП Теплоэнергия (г. Череповец), на Ковровском заводе котельного и сушильного оборудования Союз (г. Ковров), на кафедре Котлов и термодинамики Ченстоховского политехнического университета (Польша).

Основные результаты работы

1. С позиции системного подхода развиты теоретические основы синтеза математических моделей тепломассообменных аппаратов и их систем со сложной конфигурацией потоков, позволяющих вводить различный уровень декомпозиции моделируемых систем и процессов.
2. Предложена принципиально новая матричная модель для расчета характеристик систем поверхностных теплообменников со сложной конфигурацией потоков теплоносителей, учитывающая возможность фазовых переходов в теплоносителях, в которой модели подсистем (элементов) представлены конечным числом предложенных дифференциальных описаний.
3. Предложенная модель обобщена на случай расчета характеристик многоступенчатых смешивающих тепломассообменных аппаратов струйного и барботажного типов, на основе которой разработана модель деаэрации воды в многоступенчатых деаэраторах струйного и барботажного типов с учетом внутренних источников массы, обусловленных химическими реакциями.
4. Разработана система кодификации структуры систем тепломассообменных аппаратов, которая однозначно ставит в соответствие данной структуре ее код, а по коду синтезирует разработанную матричную модель системы, позволяющую рассчитать ее характеристики и оценить эффективность. С использованием предложенной системы кодификации разработан генетический (эволюционный) алгоритм структурной оптимизации сложных систем тепломассообменного оборудования.
5. Разработан алгоритмический подход к моделированию переходных процессов в тепломассообменных аппаратах и их системах со сложной конфигурацией потоков, использующий математический аппарат теории цепей Маркова, позволяющий описывать их характеристики в нестационарных режимах работы и позволяющий достаточно просто использовать эти модели в системах автоматического регулирования.
6. Выполнены экспериментальные исследования тепломассообменных процессов при струйной и барботажной деаэрации, восполняющие пробелы эмпирического обеспечения разработанной модели, и их критериальная обработка. На их основе выполнена параметрическая идентификация и проведена проверка адекватности предложенных моделей процессов тепло- и массопереноса в ряде промышленных аппаратов и расширены возможности эмпирического обеспечения методов их расчета.
7. По ряду выбранных целевых функций сформулированы и решены задачи многопараметрической и многокритериальной оптимизации систем тепломассообменных аппаратов со сложной конфигурацией потоков теплоносителей и показана нетривиальность расчетных рекомендаций по совершенствованию режимно-конструктивных параметров тепломассобменного оборудования.
8. Поставлена и решена с использованием алгоритма Беллмана нелинейная задача об оптимальном по расходу топлива распределении выработки тепловой и электрической энергии на совокупности турбоагрегатов с различными характеристиками при меняющемся потребности в тепловой и электрической энергии с учетом сложных технологических ограничений на область допустимых значений управляющих параметров.
9. Разработан алгоритм и пакет прикладных программ для реализации метода расчета теплообменных установок со сложной конфигурацией потоков. Результаты работы реализованы при разработке компьютерного модуля ТЭС-ЭКСПЕРТ для оптимального распределения нагрузки между турбоагрегатам.
10. Разработанные методики расчета и их программное обеспечение переданы Костромской ГРЭС (КГРЭС) ОАО ОГК-3, Владимирской и Дзержинской ТЭЦ ОАО ТГК-6, ОАО Северсталь, МУП Теплоэнергия г. Череповца, на Ковровский завод котельного и сушильного оборудования Союз (г. Ковров), на кафедру Котлов и термодинамики Ченстоховского политехнического университета (Польша), которые используется при выполнении исследовательских и проектных работ.

Основные положения диссертации опубликованы:

В ведущих рецензируемых журналах и изданиях (по списку ВАК):

1. Барочкин, Е. В. Системный анализ многоступенчатых теплообменных установок [Текст] / Е. В. Барочкин // Изв. ВУЗов, Химия и химическая технология. - 2005. - Т. 48. - Вып. 11. - С. 123Ц126.
2. Барочкин, Е. В., Структурно-параметрический синтез модели и системный анализ многоступенчатых  деаэраторов [Текст] / Е. В. Барочкин, В. П. Жуков, В. Е. Мизонов, Е.В. Магдиев // Химическая промышленность сегодня. - 2005. - Вып. 3. - С. 28Ц32.
3. Барочкин, Е. В. Математическое моделирование многоступенчатых теплообменников сложной конфигурации [Текст] / Е. В. Барочкин, В. П. Жуков, Г. В. Ледуховский // Изв. ВУЗов, Химия и химическая технология. - 2004. - Т. 47. - Вып. 2. - С. 45Ц47.
4. Барочкин, Е. В. Опыт работы автоматизированной установки ионообменного умягчения  природной воды [Текст] / Е. В. Барочкин, М. Ю. Опарин, А. А. Ильичев, А. Б. Ларин // Теплоэнергетика. - 2005, - № 10, - С. 18Ц23.
5. Барочкин, Е. В. Метод расчета многоступенчатых теплообменных аппаратов с учетом фазового перехода [Текст] / Е. В. Барочкин, В. П. Жуков, Г. В. Ледуховский, Х. Отвиновский // Изв. ВУЗов, Химия и химическая технология. - 2004. - Т. 47. - Вып. 2. - С. 170Ц173.
6. Барочкин, Е. В. Обобщенная модель каскадных теплообменных аппаратов с учетом фазовых переходов [Текст] / Е. В. Барочкин, В. П. Жуков, Г. В. Ледуховский // Изв.ВУЗов, Химия и химическая технология, - 2004. - Т. 47. - Вып. 3. - С. 67Ц69.
7. Барочкин, Е.В. Моделирование тепломассообмена в смешивающих подогревателях со сложной конфигурацией потоков [Текст] / Е. В. Барочкин, В. П. Жуков, Г. В. Ледуховский // Изв. ВУЗов, Химия и химическая технология. - 2004. - Т. 47. - Вып. 4. - С. 164Ц166.
8. Барочкин, Е.В. Моделирование тепломассообмена в смешивающих подогревателях со сложной конфигурацией потоков [Текст] / Е. В. Барочкин, В. П. Жуков, Г. В. Ледуховский // Изв. ВУЗов, Химия и химическая технология. - 2004. - Т. 47. Вып. 4. - С. 164Ц166.
9. Барочкин, Е.В. Матричная формализация расчета барабанных паровых котлов [Текст] / Е.В Барочкин, В. П. Жуков // Изв. ВУЗов, Химия и химическая технология. - 2006. - Т. 49. - № 10. - С. 98Ц100.
10. Барочкин, Е.В. Моделирование тепломассообмена в струйных деаэраторах со сложной конфигурацией потоков [Текст] / Е. В. Барочкин, В. П. Жуков, Г. В. Ледуховский, А. А. Борисов // Изв. ВУЗов, Химия и химическая технология. - 2004. - Т. 47. - Вып. 9. - С.76Ц79.
11. Барочкин, Е.В. Матричная модель пластинчатых теплообменников [Текст] / Е. В. Барочкин, В. П. Жуков, И. В. Степин, А. А. Борисов // Изв. ВУЗов, Химия и химическая технология. - 2007. - Т. 50. - № 9. - С. 121Ц123.
12. Барочкин, Е.В.  Применение  теории цепей  Маркова к  моделированию  тепломассообмена  в  струйных  деаэраторах  [Текст] / Е. В. Барочкин,  В. П. Жуков, Е. В. Магдиев, В. Е Мизонов // Изв. ВУЗов, Химия и химическая технология. - 2007. - Т. 50. - № 10. - С. 99Ц101.
13. Барочкин, Е.В. Обобщенный метод расчета многоступенчатых деаэраторов [Текст] / Е. В. Барочкин, В. П. Жуков, Г. В. Ледуховский, А. А. Борисов // Изв.ВУЗов, Химия и химическая технология. - 2004. - Т. 47. - вып. 9. - С. 100Ц103.
14. Барочкин, Е.В. Рециркуляция теплоносителей в теплообменных установках со сложной конфигурацией потоков [Текст] / Е. В. Барочкин, В. П. Жуков, В. Е. Мизонов, Н. Otwinovski // Изв.ВУЗов, Химия и химическая технология. - 2005. - Т. 48. - Вып. 1. - С. 124Ц128.
15. Жуков, В.П. Применение теории цепей Маркова к динамическому моделированию теплообменных аппаратов [Текст] / Е. В. Барочкин, В. П. Жуков, В. Е. Мизонов, Г. В. Ледуховский // Изв. ВУЗов, Химия и химическая технология. - 2005. - Т.48. - Вып. 4. - С. 87Ц89.

В монографии и в других журналах и материалах конференций:

16. Энергоаудит  и  энергосбережение  на ТЭС.  [Текст] / Под ред.  А.В. Мошкарина,  Е.В. Барочкина,  Ю.М. Пашковского. ИГЭУ.  Иваново.  1999. - 208 с.  ISBN-5-89482-076-6.
17. Барочкин, Е. В. Алгоритм структурно-параметрического синтеза систем тепломассобменных аппаратов со сложной конфигурацией потоков [Текст] / Е. В. Барочкин // Вестник ИГЭУ. - 2006. - Вып. 4. - С. 66-68.
18. Барочкин, Е. В. О моделировании газообмена в пароводяном тракте ТЭС [Текст] / Е. В. Барочкин, В. П. Жуков, И. А Шатова, В. Н. Виноградов, Г. В. Ледуховский // Вестник ИГЭУ. - 2006. - Вып. 2. - С. 28Ц31.
19. Барочкин, Е.В. Программный комплекс  ТЭС-ЭКСПЕРТ: опыт оптимизации режимов работы оборудования ТЭЦ [Текст] / Е. В. Барочкин, А. А. Поспелов, В. П. Жуков А. А. Андреев, Г. В. Ледуховский, А.А. Борисов // Вестник ИГЭУ. Ц2006. - Вып. 4. - С.3Ц6.
20. Барочкин, Е.В. Структурно-параметрическая оптимизация многопоточных систем пластинчатых теплообменных аппаратов [Текст] / Е. В. Барочкин, В. П. Жуков // Вестник ИГЭУ. - 2007. - Вып. 3. - С. 138Ц139.
21. Барочкин, Е.В. Расчет декарбонизации воды в многоступенчатом деаэраторе [Текст] / Е. В. Барочкин, В. Н. Виноградов, В. П. Жуков, Г. В. Ледуховский, И. В. Степин // Вестник ИГЭУ. - 2004. - Вып. 5. - С.3Ц6.
22. Ледуховский, Г.В. Метод расчета многоступенчатых теплообменников сложной конфигурации с учетом фазового перехода теплоносителей [Текст] / Г. В. Ледуховский, В. П. Жуков, Е. В. Барочкин // Вестник ИГЭУ. - 2004. - Вып. 3. - С. 138Ц139.
23. едуховский, Г.В. Оптимизация режимов работы ТЭС [Текст] / Г. В. Ледуховский, А. А. Борисов, А. А. Поспелов, Е. В. Барочкин, В. П. Жуков // Вестник ИГЭУ. - 2005. - Вып.4. - С. 170.
24. Барочкин, Е. В. Обобщенный метод расчета каскада подогревателей поверхностного и смешивающего типа [Текст] / Труды XVIII Международной конференции Математические методы в технике и технологиях ММТТ-18, - Казань. - 2005. - С. 209.
25. Барочкин, Е. В. Расчет тепломассообмена в контуре циркуляции  барабанного  котла [Текст] / Е. В. Барочкин, В. П. Жуков, В. Е. Мизонов // Региональное  приложение к журналу Современные наукоемкие технологии. - Иваново. - 2005. - С. 34Ц38.
26. Барочкин, Е.В. Расчет многоступенчатых подогревателей с произвольной структурой связей между ступенями [Текст] / Е. В. Барочкин, В. П. Жуков, Г. В Ледуховский // Труды XVII Международной конференции Математические методы в технике и технологиях ММТТ-17. - Кострома. - 2004. - Т. 9. - С. 120Ц121.
27. Барочкин, Е. В. Моделирование тепломассообмена в каскаде аппаратов сложной конфигурации [Текст] / Е. В. Барочкин, Г. В. Ледуховский, В. П. Жуков // Труды Междунар. науч. конф. Энерго- ресурсосберегающие технологии и оборудование, экологически безопасные производства. - Иваново. Ц  2004. - Т. 2. - С. 36.
28. Барочкин, Е. В. Обобщенный  метод  расчета  сложных  теплообменных деаэраторных установок [Текст] / Е. В. Барочкин, В. П. Жуков, Г. В. Ледуховский, Д. В. Русин // Материалы XII Межд. науч. техн. конф. Бенардосовские чтения Состояние и перспективы развития электротехнологии. - Иваново. - 2005. - С. 131.
29. Barochkin, E. V. Modelowanie wielostopniowych wymiennikov ciepla. Modelling of Multistage heat exchangers [Text] / E. V. Barochkin, V.P. Zhukov, H. Otwinovski, D. Urbaniak // VI Miedzynarodowa konferencja naukowa Nowe technologie i osiagniecie w Metalurgii i inzynierii Materialowej. - Politechnika Czestochowska. - 2005. - Р. 112-116.
30. Барочкин, Е. В. Выбор схемы включения прямоточных деаэрационных устройств [Текст] / Е. В. Барочкин, И. А. Шатова, Г. В. Ледуховский. // Труды IV Российской научно-практическая конф. Повышение эффективности теплоэнергетического оборудования. - Иваново. - 2005. - С. 63 - 65.
31. Барочкин, Е. В. Задача системного анализа пластинчатых теплообменных установок [Текст] / Е. В. Барочкин, В. П. Жуков., М. Ю. Зорин // Материалы XIV Межд. науч. техн. конф. Бенардосовские чтения Состояние и перспективы развития электротехнологии. - Иваново. - 2007. - Т. 2. - С. 183.
32. Барочкин, Е. В. Оценка влияния октадециламина на деаэрацию химочищенной воды [Текст] / Е. В. Барочкин, И. А. Шатова, Г. В. Ледуховский // Труды IV Российской научно-практической конференции Повышение эффективности теплоэнергетического оборудования. - Иваново. Ц2005. - С. 66Ц69.
33. Барочкин, Е. В. Опыт наладки деаэрационно-питательных установки при резко переменных нагрузках [Текст] / Е. В. Барочкин, В.Н. Виноградов, Г. В. Ледуховский // Материалы XII Межд. науч. техн. конф. Бенардосовские чтения Состояние и перспективы развития электротехнологии. Ц  2005. - Т. 1. - С. 137.
34. Барочкин Е. В. Повышение эффективности атмосферной деаэрации с использованием охладителей выпара [Текст] / Е. В. Барочкин., В. Н. Виноградов, Г. В. Ледуховский // Материалы XII Межд. науч. техн. конф. Бенардосовские чтения Состояние и перспективы развития электротехнологии. - Иваново. - 2005. - Т. 1. - С. 138.
35. Барочкин Е. В. Моделирование, оптимизация и управление энерго- и ресурсосберегающими технологиями в  сложных  теплообменных системах [Текст] / Е. В. Барочкин, В. П Жуков // Труды Междунар. научн. конф. Теоретические основы создания оптимизации и управления энерго- ресурсосберегающие процесссами и оборудованием. - Иваново. - 2007. - Т 1. - С. 219Ц227.
36. Барочкин, Е. В. Метод расчета многоступенчатого деаэратора с учетом декарбонизации воды [Текст] / Е. В. Барочкин, В. П. Жуков, И. В. Степин, В. Н. Виноградов // Сборник трудов XX межд. конф. Математические методы в технике и технологии ММТТ-20. - Ярославль. - 2007. - Т. 5. - С. 158Ц159 .
37. Жуков, В. П. Структурная и  режимная  оптимизация сложных технологических схем измельчения [Текст] / В. П. Жуков, А. Е. Иванов, В. Е. Мизонов, Е. В. Барочкин. // Тез. докл. XI Межд. научн. конф. Математические методы в химии и технологиях ММХТ-XI. - Владимир. - 1998. - Т. 3. - С. 106.
38. Жуков, В. П. Динамическое моделирование теплообменных аппаратов [Текст] / В. П. Жуков, Е. В. Барочкин, В. Е. Мизонов, Г. В. Ледуховский. // Материалы XII Межд. науч. техн. конф. Бенардосовские чтения Состояние и перспективы развития электротехнологии. - Иваново. - 2005. - С. 123.
39. Ледуховский, Г.В. Моделирование тепломассообмена в многоступенчатых струйных деаэраторах [Текст] / Г.В. Ледуховский, В. П. Жуков, Е. В. Барочкин // Труды XVIII Международной конференции Математические методы в технике и технологиях ММТТ-18. - Казань. - 2005. - С. 210.
40. Межеумов, Г.Г. Математическое моделирование процессов измельчения с переменным кинетическим потенциалом [Текст] / Г.Г. Межеумов, В.П. Жуков, Е.В. Барочкин, H. Otwinowski // Сборник трудов XIV межд. конф. Математические методы в технике и технологии ММТТ-14. - Смоленск. - 2001. - Т. 3. - С. 99 .
41. Жуков, В.П. Матричная формализация расчета многоступенчатой испарительной установки [Текст] / Жуков В.П., Мошкарин А.В., Барочкин Е.В. // Материалы XI Межд. науч. техн. конф. Бенардосовские чтения Состояние и перспективы развития электротехнологии. - Иваново. - 2003. - Т. 1. - С. 111.
42. Козлов, М. А. Математическая модель агломерации сыпучих материалов в кипящем слое [Текст] / М. А. Козлов, В. П. Жуков, Л. Н. Овчинников, Е. В. Барочкин // Сборник трудов XVI межд. конф. Математические методы в технике и технологии ММТТ-16. - Санкт-Петербург. - 2003. - С. 114.
43. Ледуховский, Г. В. Расчетно-экспериментальные исследования процессов тепломассообмена в условиях струйно-барботажной деаэрации воды [Текст] / Г. В. Ледуховский, Е. В. Барочкин, В. П. Жуков, В. Н. Виноградов // Материалы XIII Межд. науч. техн. конф. Бенардосовские чтения Состояние и перспективы развития электротехнологии. - Иваново. - 2006. - С. 56Ц57.
44. Магдиев, Е. В. Применение теории цепей Маркова к моделированию процессов в теплообменных аппаратах [Текст] / Е. В. Магдиев, В. П. Жуков, Е. В. Барочкин // Материалы XIII Межд. науч. техн. конф. Бенардосовские чтения Состояние и перспективы развития электротехнологии. Ц  Иваново. - 2006. - С. 53 - 55.
45. Жуков, В. П. Системный анализ многоступенчатых деаэраторов струйного типа [Текст] / В. П. Жуков, В. Е. Мизонов, Е. В. Барочкин // Вестник Международной академии системных исследований. - Москва. - 2006. - Т. 9. - Ч. 1. - С. 45Ц50.
46. Магдиев, Е. В. Моделирование переходных процессов в теплообменных аппаратах  смешивающего типа на основе  цепей теории Маркова [Текст] / Е. В. Магдиев, В. П. Жуков, Е. В. Барочкин // Труды Междунар. науч. конф. Теоретические основы создания оптимизации и управления энерго- ресурсосберегающие процесссами и оборудованием. - Иваново. - 2007. - Т. 2. - С. 80Ц82.
47. Магдиев, Е. В. Моделирование  переходных  процессов в теплообменных аппаратах с учетом  фазового перехода [Текст] / Е. В. Магдиев, В. П. Жуков, Е. В. Барочкин. // Сборник трудов XX межд. конф. Математические методы в технике и технологии ММТТ-20. - Ярославль. - 2007. - Т. 5. - С. 115Ц116.
48. Жуков, В. П. Метод расчета пластинчатых теплообменных аппаратов [Текст] / В. П. Жуков, Е. В. Барочкин, Д. Г. Денисов, И. В. Степин // Сборник трудов XX межд. конф. Математические методы в технике и технологии ММТТ-20. - Ярославль. - 2007. - Т. 5. - С. 114Ц115.
49. Pastucha, L. Heat exchanger modelling based on Markov chains [Text] / L. Pastucha, H. Otwinowski, E. V. Barokin, V.P. ukov [Text] // Proceedings 34th International Conferense of Slovak of Chemicfl Engineering. - Slovakia. - 2007. - Po-We-4. - S 91. - pdf.
50. Мошкарин, А. В. Выбор основного и вспомогательного оборудования ТЭС. Методическое пособие. [Текст] / А. В. Мошкарин, Е. В. Барочкин, М. Ю. Зорин, Г. В. Ледуховский. - Иваново. - ИГЭУ. - 2004. - 56 с.
51. Барочкин, Е. В. Тепловые схемы и показатели работы теплофикационных паротурбинных установок ТЭС и АЭС. Методическое  пособие. [Текст] / Е.В. Барочкин, А. В. Мошкарин, А. Л. Виноградов. - Иваново. - ИГЭУ. - 2005. - 68 с.

Список условных обозначений:

t, oC - температура теплоносителя; , oC - температурный напор; F, м2 - поверхности нагрева; k, Вт/ oC м2 Цкоэффициент теплопередачи; km, кг/с м2 - коэффициент массопередачи; с, Дж/(кг oC) - удельная теплоемкость; G, кг/с - массовый расход теплоносителя; Q,м3/с - объемный расход теплоносителя; V, м3 - объем; хЦстепень сухости; r, Дж/кг - удельная теплота парообразования;w, м/с - скорость движения теплоносителя; p - вероятность.

Матрицы:

ВЦматрица процесса в ступени (подсистеме); А Ц матрица процесса в системе; G Ц матрица связи для горячего теплоносителя; X Ц матрица связи для холодного теплоносителя; GX Ц матрица связи для горячего и холодного теплоносителя; K Ц матрица коммутации; T, Х Ц вектор состояния; Р Ц переходная матрица.

Индексы: 0 - начальные условия, известные параметры; 1 - горячий теплоноситель; 2 - холодный теплоноситель; g - газ, nЦнасыщение.

БАРОЧКИН ЕВГЕНИЙ ВИТАЛЬЕВИЧ

анализ и оптимальный синтез теплообменных

систем со сложной конфигурацией потоков

в энергетических и химических комплексах

Автореферат

диссертации на соискание ученой степени

доктора технических наук

ицензия ИД № 05285 от 4 июля 2001 г.

Подписано в печать 28.01.2008. Формат 5084 1/16

Печать плоская. Усл. печ. л. 1,86.

Тираж 100 экз. Заказ № 151.

     Авторефераты по всем темам  >>  Авторефераты по техническим специальностям