Авторефераты по всем темам >>
Авторефераты по физике
На правах рукописи
СУХИХ Андрей Анатольевич
ИССЛЕДОВАНИЕ ТЕРМОДИНАМИЧЕСКИХ СВОЙСТВ И ТЕПЛОТЕХНИЧЕСКИХ ХАРАКТЕРИСТИК ФТОРОРГАНИЧЕСКИХ РАБОЧИХ ВЕЩЕСТВ
Специальность - 01.04.14 Теплофизика и теоретическая теплотехника
АВТОРЕФЕРАТ
диссертации на соискание ученой степени доктора технических наук
МОСКВА - 2012
Работа выполнена в ФГБОУ ВПО Национальный исследовательский университет МЭИ на кафедре Теоретических основ теплотехники им.
М.П. Вукаловича.
Официальные оппоненты:
доктор технических наук, профессор Архаров Иван Алексеевич, проректор по международным связям МГТУ им Э. Н. Баумана;
доктор технических наук, профессор Герасимов Анатолий Алексеевич, заведующий кафедрой Теплогазоснабжения и вентиляции ФГБОУ ВПО Калининградский государственный технический университет;
доктор технических наук, профессор Цветков Олег Борисович, заведующий кафедрой Теоретических основ тепло- и хладотехники ФГБОУ ВПО Санкт-Петербургский НИУ Информационных технологий, механики и оптики.
Ведущая организация: ОАО ЭНИН им. Кржижановского
Защита состоится 19 октября 2012 г. в 10:00 на заседании диссертационного совета Д 212.157.04 при ФГБОУ ВПО Национальный исследовательский университет Московский энергетический институт по адресу: г. Москва, ул.
Красноказарменная, д. 14, Малый актовый зал.
Отзывы на автореферат с подписями, заверенными печатью учреждения, просим направлять по адресу: 111250, Москва, Красноказарменная ул., д.14, Ученый Совет ФГБОУ ВПО НИУ МЭИ.
С диссертацией можно ознакомиться в научно-технической библиотеке ФГБОУ ВПО Национальный исследовательский университет МЭИ.
Автореферат разослан л 2012 г.
Ученый секретарь диссертационного совета Д 212.157.к.т.н. Ястребов А.К.
ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ
Актуальность темы Выбор в качестве объекта исследований фторорганических рабочих веществ (РВ) обусловлен необходимостью решения как ряда теплофизических, так и теплотехнических проблем при внедрении новых рабочих тел в теплосиловых установках специального назначения в сфере малой энергетики, атомной отрасли, а также в теплонасосной технике. Проведение теплофизических исследований данных веществ является актуальным, обоснованным и необходимым при разработке новых энергоустановок с рабочим телом фторорганического состава.
Особое внимание в работе уделено классу фторуглеродов.
Термодинамические свойства некоторых фторуглеродов хорошо изучены и представлены для низкотемпературного диапазона, а при повышенных температурах в известной литературе такие данные отсутствовали. Поэтому постановка и проведение исследований в широком диапазоне параметров были необходимы и актуальны.
Применение в теплосиловых энергоустановках в качестве рабочих веществ фторуглеродов является новаторским и требует выполнения комплекса научных исследований, опытных и расчетных проработок, в том числе проведения теплотехнических испытаний. Разработка и термодинамический анализ теплосиловых циклов на предлагаемых рабочих веществах актуальны и необходимы для подготовки этапа проектирования.
На современном этапе развития энергетики всё больший вес приобретают проекты в сфере малой, или распределенной, энергетики. На первый план выходят вопросы эффективного использования энергетических ресурсов, снижения энергозатрат в промышленности и коммунальном хозяйстве, ведутся непрерывные разработки лальтернативных, нетрадиционных источников энергии. Одним из таких направлений является разработка электрогенерирующих установок мощностью в диапазоне 100 кВт - 5 МВт на неводных РВ, работающих на основе использования вторичных и возобновляемых источников энергии. Подобные установки в полной мере отвечают современной концепции децентрализованного электроснабжения и могут использоваться как мини-электростанции для небольших посёлков, промышленных объектов, удалённых от электросетей, а также как вспомогательные, резервные, либо аварийные источники электроэнергии автономных объектов различного назначения.
Результаты настоящей работы призваны доказать, что электрогенерирующие установки такого назначения на фторорганическом РВ будут иметь более высокий уровень как термодинамической эффективности, так и технологической безопасности по сравнению с установками на водяном паре либо на углеводородах, работающих по так называемому лорганическому циклу Ренкина.
Существенной частью диссертационной работы является также исследование термодинамических свойств и теплотехнических характеристик фторорганических рабочих веществ для теплонасосной техники. При внедрении нового поколения рабочих веществ парокомпрессионных теплонасосных установок (ТНУ) необходимо решать вопросы, связанные с обоснованием выбора новых рабочих тел, определением энергетической эффективности циклов, интенсификацией тепломассообмена в испарителе/конденсаторе и др. Постановка и решение данных проблем весьма актуальны. Тематика данной работы развивалась при выполнении ряда хоздоговорных и бюджетных работ в течение последних лет. Автор был научным руководителем более десяти хоздоговорных работ (в том числе, международного контракта с Praxair Inc. (USA) и российскими заказчиками:
НПКФ ЭКИП, ОАО НИКИЭТ, ПКФ ОАО Концерн Энергоатом), а также трех грантов РФФИ.
Основные результаты работы вошли в состав комплекса исследований, представленного на соискание премии Правительства РФ коллективом под научным руководством профессора Сычева В.В., и были высоко оценены научным сообществом. Автор настоящей диссертационной работы стал лауреатом премии Правительства РФ в области науки и техники 2008 года за Разработку и внедрение комплекса прецизионных данных о теплофизических свойствах рабочих веществ криогенной и холодильной техники и тепловых насосов.
Целью работы является разработка теплофизических основ и теплотехнических баз данных, необходимых для внедрения новых рабочих тел в теплосиловые установки специального назначения в сфере малой энергетики, атомной отрасли, а также в теплонасосную технику.
Ход исследования и структуру диссертационной работы определили поставленные задачи:
1) проведение измерений рvT соотношений предлагаемых к внедрению новых рабочих веществ;
2) построение на их основе экспериментально-обоснованных расчетных моделей;
3) разработка процессов, циклов, схем теплосиловых и теплонасосных установок на фторорганических рабочих веществах;
4) проведение теплотехнических испытаний прототипов основных конструкционных элементов таких энергоустановок на новых и традиционных рабочих телах.
Научная новизна работы заключается в следующем:
1. Получены опытные данные о плотности октафторпропана (C3F8) в сверхкритической области и области перегретого пара при повышенных температурах (в интервале температур от 373 до 773 К и давлении до 10 МПа).
2. Получены опытные данные о плотности декафторбутана (C4F10) в сверхкритической области и области перегретого пара при повышенных температурах до 773 К и давлениях до 10 МПа; получены опытные данные о плотности в области жидкости, а также - при умеренных температурах вблизи кривой конденсации, по которым рассчитаны ортобарические плотности и .
3. Разработаны уравнения состояния октафторпропана и декафторбутана вириального типа для газовой фазы в диапазонах до 12 МПа и до 823, 15 К.
4. Получены опытные данные для фторэфира HFE347mcc по давлению ps насыщения в диапазоне температур Т = 313,15Е353,15 К и экспериментальные pvT -данные в диапазоне температур от 333,15 до 353,15 К и давлений от 105,93 до 372,56 кПа.
5. Получена аппроксимационная зависимость для давления насыщения фторэфира HFE347mcc для диапазона Т = 260Е437,7 К. Для кривой насыщения построено уравнение состояния кубического типа КарнаханаЦСтарлингаЦде Сантиса с использованием как собственных данных на кривой конденсации, так и данных японских исследователей. Для расчета термодинамических величин в газовой фазе было разработано уравнение состояния эфира HFE-347mcc вириального типа.
6. Экспериментальные pvT -данные для бинарных смесей фторэфира HFE347mcc с хладонами FС218, HFC23 и FC14 в газовой фазе и на кривых конденсации являются на сегодняшний день единственными.
7. Разработаны уравнения состояния смесевых хладагентов на основе фторэфира HFE347mcc в форме вириального разложения. По разработанным уравнениям были рассчитаны таблицы термодинамических свойств смесевых хладагентов в газовой фазе.
8. Получены экспериментальные pvT - данные в паровой фазе и на кривой конденсации системы SF6 - C3F8 диапазоне температур Т = 238,15Е333,15 К при давлениях до 3500 кПа для двух составов.
9. Проанализированы термодинамические особенности циклов на фторуглеродах и сформулированы предпосылки для внедрения таких веществ в качестве рабочих тел теплосилового контура реакторной установки на быстрых нейтронах.
10. Выполнен анализ термодинамической эффективности схем ПГУ при замещении водяного контура ПТУ на фторуглеродный. Сформулирован ряд технологических и конструкционных преимуществ, возможных при внедрении фторуглеродов в качестве рабочих веществ теплосиловых установок.
11. Разработаны тепловые схемы и циклы установок для утилизации теплоты генераторного газа.
12. Сформулированы общие особенности теплосиловых циклов на фторуглеродах.
13. С целью подтверждения расчетно-теоретических разработок был построен циркуляционный стенд на октафтопропане (C3F8) тепловой мощностью котельного агрегата 80 кВт в качестве маломасштабного прототипа реальной теплосиловой установки на фторуглеродном рабочем веществе.
14. Результаты ресурсных и теплотехнических испытаний позволили определить характерные режимы работы основных аппаратов стенда и подтвердили высокую термическую стабильность октафторпропана.
15. Разработана математическая модель процессов термогидродинамики в регенеративных теплообменниках стендовой установки.
16. Исследованы теплотехнические характеристики ТНУ на диоксиде углерода, а также на хладагентах R134a и R22.
17. Разработаны оригинальные схемы и теплообменные аппараты ТНУ, на основные технические решения которых получены авторские свидетельства.
Получены экспериментально обоснованные соотношения для расчета теплоотдачи при кипении и конденсации выбранных хладагентов на рельефных поверхностях теплообмена в элементах ТНУ.
18. На основе собственных теплотехнических измерений выполнен термодинамический анализ и произведено сопоставление эффективности теплонасосных установок на диоксиде углерода и наиболее распространенных хладагентах R134a и R22.
19. Выполнено исследование процесса пузырькового кипения в большом объёме на микроструктурированных поверхностях, изготовленных с помощью новых технологий и рекомендуемых в качестве рабочих элементов основных теплообменных аппаратов ТНУ. Эксперименты проведены при кипении двух рабочих веществ: хладона R134a и теплоносителя фторуглеродного состава FC3284 в большом объеме при различных давлениях и тепловых нагрузках.
Практическая значимость и ценность научной работы Разработанные расчетные модели и полученные данные о калорических свойствах октафторпропана и декафторбутана достаточно надежны во всей области заявленных параметров, и могут быть использованы при проектировании термодинамических циклов и для расчетов процессов и аппаратов энергетических установок.
Проведенный термодинамический анализ показал, что применение фторуглеродов в качестве рабочих веществ теплосиловых циклов АЭС позволит реализовать сверхкритический прямой теплосиловой цикл с давлением до 14 МПа и температурой перед турбиной до 580 оС, с конденсацией пара при избыточных давлениях. Конфигурация такого цикла с расширением РВ в турбине на перегретом паре и сжатием жидкой фазы насосом обеспечивает высокую интегральную температуру подвода теплоты и объединяет преимущества газовых циклов (Брайтона) и пароводяных (Ренкина).
В настоящей работе сформулированы термодинамические особенности циклов на фторуглеродах и разработаны теплотехнические основы для внедрения таких веществ в качестве рабочих тел теплосилового контура реакторной установки. Предполагается, что внедрение фторуглеродов в качестве рабочих тел теплосилового контура АЭС на быстрых нейтронах позволит повысить не только термодинамическую эффективность, но и экологическую и технологическую безопасность энергоблока АЭС.
Анализ термодинамической эффективности схем ПГУ при замещении водяного контура ПТУ на фторуглеродный, проведенный на расчетной базе (экспериментально-обоснованных уравнениях состояния), подготовленной на предыдущих этапах настоящей работы, показал, что практическая реализация одной из этих схем обеспечит существенное увеличение действительного КПД цикла.
Практически всё разрабатываемое оборудование в контуре фторуглеродной ПТУ не имеет промышленных аналогов. Однако высокая термодинамическая эффективность предлагаемых схем, наличие отечественной базы по производству фторуглеродов, возможность существенного удешевления их производства делает подобные проекты вполне реальными.
Тепловые схемы и циклы установок по утилизации теплоты генераторного газа, разработанные автором, могут быть использованы на этапе проектирования опытно-промышленных образцов.
На основные технические решения по предлагаемым схемам и теплообменным аппаратам ТНУ получены авторские свидетельства. Получены экспериментально обоснованные соотношения для расчета теплоотдачи при кипении и конденсации исследованных хладагентов на рельефных поверхностях теплообмена в элементах ТНУ. Экспериментально подтверждена возможность повышения коэффициента теплопередачи при использовании таких поверхностей в 3-5 раз. Полученные результаты могут быть использованы при проектировании нового поколения теплонасосных установок, работающих в составе систем тепло- и холодоснабжения жилых и административных зданий, а также технологических объектов различного назначения.
Результаты исследования проблем интенсификации процессов тепломассообмена в аппаратах ТНУ имеют непосредственное прикладное значение. Интенсификация теплообмена в испарителе/конденсаторе ТНУ является одним из способов повышения компактности и экономичности установки в целом.
В настоящей работе был выполнен комплекс исследований интегральных коэффициентов теплоотдачи в основных аппаратах ТНУ и локальных коэффициентов теплоотдачи на микроструктурированных поверхностях, предлагаемых в качестве рабочих элементов испарителей и конденсаторов ТНУ.
На защиту выносятся:
1. Экспериментальные данные о плотности октафторпропана (C3F8) в сверхкритической области и области перегретого пара при повышенных температурах (в интервале температур от 373 и до 773 К и давлении до 10 МПа).
2. Экспериментальные данные о плотности декафторбутана (C4F10) в сверхкритической области и области перегретого пара при повышенных температурах до 773 К и давлениях до 10 МПа; экспериментальные данные о плотности в области жидкости, а также при умеренных температурах вблизи кривой конденсации.
3. Уравнения состояния октафторпропана и декафторбутана вириального типа для газовой фазы в диапазонах давления до 12 МПа и температуры до 823,15 К.
4. Опытные данные для фторэфира HFE347mcc по давлению насыщения рs в диапазоне температур Т = 313,15Е353,15 К и экспериментальные pvT -данные в диапазоне температур от 333,15 до 353,15 К и давлений от 105,93 до 372,56 кПа.
5. Уравнение состояния кубического типа КарнаханаЦ СтарлингаЦде Сантиса для кривых насыщения и аппроксимационная зависимость логарифмического типа для давления насыщения фторэфира HFE347mcc в диапазоне Т = 260Е437,7 К.
Уравнение состояния эфира HFE-347mcc вириального типа для расчета термодинамических величин в газовой фазе.
6. Экспериментальные pvT -данные для бинарных смесей фторэфира HFE347mcc с хладонами FС218, HFC23 и FC14 в газовой фазе и на кривых конденсации.
7. Уравнения состояния смесевых хладагентов на основе фторэфира HFE347mcc.
8. Экспериментальные pvT - данные в паровой фазе и на кривой конденсации системы SF6 - C3F8 в диапазоне температур Т = 238,15Е333,15 К при давлениях до 3500 кПа для двух составов.
9. Теплосиловые циклы и схемы на фторуглеродных рабочих веществах:
второго контура реакторных установок на быстрых нейтронах, установок для утилизации генераторного газа и выхлопов ГТУ.
10. Циркуляционный стенд на октафторпропане (C3F8) тепловой мощностью котельного агрегата 80 кВт и результаты его теплотехнических испытаний.
11. Математическая модель процессов термогидродинамики в регенеративных теплообменниках стендовой установки на октафторпропане.
12. Схемы и теплообменные аппараты теплонасосных установок (ТНУ) и результаты их теплотехнических испытаний.
13. Экспериментально обоснованные соотношения для расчета теплоотдачи при кипении и конденсации фторорганических хладагентов на рельефных поверхностях теплообмена в элементах ТНУ.
14. Результаты термодинамического анализа и сопоставления эффективности теплонасосных установок на диоксиде углерода и наиболее распространенных хладагентах R134a и R22.
15. Результаты исследования процесса пузырькового кипения хладона R134a и теплоносителя фторуглеродного состава FC-3284 в большом объеме при различных давлениях и тепловых нагрузках в большом объёме на микроструктурированных поверхностях, изготовленных с помощью новых технологий и рекомендуемых в качестве рабочих элементов основных теплообменных аппаратов ТНУ.
Апробация работы Основные положения и результаты работы докладывались и обсуждались:
1. На семинаре международной теплофизической школы (г. Тамбов: ТГТУ.
1995).
2. На международных научно-технических конференциях студентов и аспирантов, проводимых на базе МЭИ (2002-2003 гг.).
3. На Х Российской конференции по теплофизическим свойствам веществ (г.
Казань: КГТУ. 2002 г.).
4. На Всероссийской научно-практической конференции Человеческое измерение в информационном обществе (Москва: ВВЦ, павильон №57.
Минобразования РФ. 2003 г.).
5. На 2-ой международной научно-практической конференции Низкотемпературные и пищевые технологии в XXI веке (г. Санкт-Петербург:
ГУНТиПТ. 2003 г.).
6. На 2-ой международной научно-практической конференции: Природные холодильные агенты - альтернатива глобальному потеплению (г. СанктПетербург: ГУНТиПТ. 2003 г.).
7. На республиканском научно-техническом семинаре Тенденции развития теоретической теплотехники. Создание современных средств и технологий в теплоэнергетике (Ташкент: ТГТУ. 2004 г.).
8. На электронной конференции по подпрограмме Топливо и энергетика научно-технической программы Научные исследования высшей школы по приоритетным направлениям науки и техники (Москва: МЭИ. 2004 г.).
9. На 11-ой Российской конференции по теплофизическим свойствам веществ (г. Санкт-Петербург. 2005 г.).
10. На XI национальной конференции по теплоэнергетике НКТЭ - 2006 (г.
Казань. 2006г.).
11. На XIII Российской конференции по теплофизическим свойствам веществ (г. Новосибирск. 2011 г.).
12. На национальной конференции Повышение эффективности, надежности и безопасности работы энергетического оборудования ТЭС и АЭС. (Москва:
ФБГОУ ВПО НИУ МЭИ, 2012 г.).
ичный вклад автора В основу диссертации легли результаты исследований, выполняемые автором на кафедре Теоретических основ теплотехники им. М.П. Вукаловича ФГБОУ ВПО НИУ МЭИ в течение 13 лет. В диссертации обобщены результаты исследований, выполненных автором как самостоятельно, так и в качестве научного руководителя соискателей и аспирантов.
ичное участие автора заключается в постановке и реализации комплексного экспериментального и расчетно-теоретического исследования:
Х постановка задач и выбор способов их решения;
Х создание экспериментальной установки л pvTх , разработка и реализация метода изохорически связанных расширений;
Х экспериментальное исследование термодинамической поверхности как индивидуальных, так и смесевых рабочих веществ;
Х обработка, обобщение и анализ опытных результатов;
Х разработка экспериментально-обоснованных уравнений состояния на базе, в основном, собственных данных;
Х разработка и анализ термодинамических циклов на фторуглеродных рабочих веществах;
Х создание циркуляционного стенда на октафторпропане и проведение теплотехнических испытаний;
Х создание экспериментальных теплонасосных установок для исследования теплотехнических характеристик фторорганических рабочих веществ;
Х разработка оригинальных конструкций испарителя и конденсатора теплонасосных установок, обеспечивающих интенсификацию процессов теплообмена;
Х проведение теплотехнических измерений на теплонасосных установках;
Х разработка рекомендаций для постановки дальнейших НИР и дальнейшего развития направления повышение энергоэффективности теплонасосных установок;
Х научное руководство комплексом работ по исследованию процессов кипения фторорганических рабочих веществ на микроструктурированных поверхностях.
Степень достоверности результатов проведенных исследований, выводов и рекомендаций, сформулированных в диссертации Достоверность научных положений, выводов и рекомендаций, сформулированных в диссертации, подтверждается экспериментальным материалом, полученным с достаточно высокой точностью при использовании современных средств измерения и на основе апробированных методов исследования. Используемые в работе измерительные элементы тарировались и испытывались на достоверность в условиях, описанных в соответствующих методических инструкциях, и обеспечили высокий метрологический уровень полученных опытных данных.
Сопоставление рассчитанных по экспериментально-обоснованным уравнениям состояния термодинамических свойств и известных из публикаций, например теплоемкостей, энтальпий, скоростей звука октафторпропана, проводилось в диапазонах, имеющих статус стандартных справочных данных, и подтвердило достоверность построенных моделей.
При проведении экспериментальных работ по исследованию термодинамической поверхности рабочих веществ использовались различные установки и независимые методики обработки опытных данных. Использовались также специальные методы проверки внутренней согласованности полученных данных.
Математическая модель термогидравлического расчета регенератора маломасштабной теплосиловой установки на фторуглеродном рабочем веществе была верифицирована на базе опытных данных. Термическая и коррозионная стойкость октафторпропана как рабочего вещества была подтверждена результатами спектрометрического анализа.
Достоверность расчетов схем утилизационных и теплонасосных установок, а также результатов теплотехнических испытаний подтверждалась согласованием тепловых балансов их основных узлов и аппаратов.
Основным методическим подходом, реализованным в диссертационной работе, является экспериментальное обоснование всех теоретических выводов и расчетных результатов.
Публикации По теме диссертации опубликовано 35 научных работ, в том числе, 13 в журналах, включенных в перечень ВАК РФ, 20 - тезисы докладов и статьи, опубликованные в материалах конференций, 2 - изобретения по теме диссертации.
Структура и объем работы Диссертационная работа состоит из введения, 9 разделов, заключения, списка использованной литературы из 163 наименований, приложений с таблицами рассчитанных термодинамических свойств исследованных веществ. Содержание работы изложено на 380 страницах, включает 123 таблицы и 105 рисунков.
СОДЕРЖАНИЕ РАБОТЫ
ВВЕДЕНИЕ Обоснована актуальность темы, определены объекты и сформулирована цель исследования, изложены научные и практические результаты, выносимые на защиту.
Поиск заменителей водяных и газовых рабочих тел теплосиловых установок специального назначения проводился в классе веществ, относящихся к фторуглеродам (C3F8, C4F8, C4F10Е). Для решения указанных задач использовались pT две экспериментальные установки: л (исследование термодинамической pvTх поверхности при повышенных температурах), л (исследование термических параметров на кривой конденсации и в газовой области), а также была создана и использовалась установка ТСС-20 (циркуляционный стенд на октафторпропане) тепловой мощностью до 80 кВт (исследование циклов и аппаратов - прототипов промышленных теплосиловых установок). На основе результатов проведенных термодинамических исследований разрабатывались уравнения состояния, рассчитывались процессы, циклы, тепловые схемы энергоустановок. Результаты теплотехнических и ресурсных испытаний циркуляционного стенда на октафторпропане были получены впервые и подтвердили работоспособность, термическую стабильность, высокую эффективность фторуглеродного РВ. Были также измерены интегральные коэффициенты теплоотдачи, термогидравлические характеристики основных аппаратов стенда, т.е. подготовлена база для внедрения фторуглеродов в качестве РВ теплосиловых установок. Решение подобных задач необходимо и весьма актуально при переходе к этапу проектирования электрогенерирующих установок широкого ряда мощностей.
Поиск заменителей озоноразрушающих и радиационно-активных хладонов (R114, R11, R12,Е) проводился в классах фторуглеродов CyFx, фторпропанов C3FxH8-x, фторэфиров HFE и их смесей, в том числе и с гексафторидом серы SF6, обладающих предельно малыми значениями экологических потенциалов.
Экспериментальное исследование рvT свойств гептафторбутанолового эфира HFE347mcc и его смесей c R218, R14 и R23 проводилось в соответствии с концепцией продвижения на рынок озонобезопасных рабочих тел с низкими потенциалами глобального потепления. Для решения указанных задач pvTх использовались, как вышеупомянутая экспериментальная установка л , так и теплонасосные установки ТН300У на R134a, ТН-14 на R22, ТН-18 на диоксиде углерода (исследование циклов и аппаратов ТНУ) и многофункциональный комплекс программ обработки результатов теплотехнических измерений.
Для исследования проблем интенсификации теплообмена в основных аппаратах ТНУ была использована экспериментальная установка Института Тепловой и Энергетической Техники Университета Падеборна (Германия), разработанная для измерения кривых кипения на трубах с микроструктурированными поверхностями. В качестве исследуемых веществ были выбраны: широко распространенный в холодильной и теплонасосной технике фреон R134a и синтезированный теплоноситель фторуглеродного состава с высокой смачиваемостью FC-3284.
Расчетно-теоретическое исследование в настоящей работе проводилось с помощью уравнений состояния, разработанных, в основном, на базе собственных экспериментальных данных и, в меньшей мере, данных других авторов. Расчетные модели позволили провести разработку процессов и циклов энергоустановок на фторорганических рабочих веществах, оценить их термодинамические и теплотехнические характеристики. Новые модели и опытные данные существенно расширили базу для проектирования энергетических установок на новых рабочих телах фторорганического состава. Результаты исследований предназначены, в первую очередь, для внедрения в системах малой распределенной энергетики.
1 ЭКСПЕРИМЕНТАЛЬНЫЕ МЕТОДЫ ИССЛЕДОВАНИЯ ТЕРМОДИНАМИЧЕСКИХ СВОЙСТВ ФТОРОРГАНИЧЕСКИХ РАБОЧИХ ВЕЩЕСТВ pvTх Экспериментальная установка л и метод изохорически связанных последовательных расширений. Для исследования кривой насыщения рS(T) и получения pvT-зависимости в перегретых и насыщенных парах индивидуальных и смесевых рабочих веществ была создана экспериментальная установка л pvTх , схема которой представлена на рис.1. На установке реализуется метод последовательных расширений опытного образца, как в классическом, так и в специально разработанном изохорически связанном, вариантах.
Основным элементом опытной установки является пьезометрическая система, которая состоит из: двухкамерной ячейки равновесия (VА+VВ), дифференциального мембранного блока (ДМБ), циркуляционного контура с микронасосом и трех вентилей постоянного объема. Дифференциальный мембранный блок находится в малом объеме VА, его электрическая часть - снаружи ячейки; ДМБ отделяет исследуемую среду от буферного газа (Ar) в системе измерения давления; имеет линейную характеристику в диапазоне разности давлений от -1 до +1кПа. Порог чувствительность ДМБ менее 2 Па.
Установка оснащена системами: измерения давления; термостатирования и регулирования температуры; вакуумирования; анализа состава смеси.
Манометрический комплекс оснащен грузопоршневыми манометрами МП2,5 класса 0,02; МП-6, МП-60, МП-600 класса 0,05 (с газомасляными разделителями) и электрическими манометрами с цифровым выходом ИПДЦ-1 и ИПДЦ-2 класса 0,06.
Рис. 1 - Схема экспериментальной установки л pvTх Блок пьезометров вместе с соединительными трубопроводами, вентилями и циркуляционным насосом установлен в стандартном криостате КГ-15,150-1. В качестве термостатирующей жидкости при температуре от 180 до 350 К используется кремнийорганическая жидкость ПМС-1,5Р. Электродвигатель мешалки находится под крышкой криостата и имеет автономное охлаждение корпуса водой. Все выводы механических и электрических линий из криостата герметизированы - это исключает возможность выпаривания термостатирующей жидкости (загрязнение окружающей среды) или конденсации паров из воздуха в криостат. Термометрическая система создана на базе десятиомного образцового термометра сопротивления ПТС-10 и потенциометра Р-348 класса точности 0,001.
Для поддержания температуры в статическом режиме используется специальная система, собранная на базе прецизионного регулятора ПРТ-2М. Состав газовых смесей контролировался с помощью газового хроматографа ЛХМ-80/1 в комплекте с электронным интегратором ИЦ-26.
Двухкамерная пьезометрическая система позволяет реализовать классический метод последовательных изотермических расширений опытного образца. При использовании этого метода основное расчетное соотношение для коэффициента сжимаемости имеет следующий вид:
j j Zj = pj (Z0 / p0)N, (1) i i=где N= (VВ + VА)/VВ - геометрическая константа - величина, характеризующая пьезометрическую систему экспериментальной установки; (p0/Z0) j = RT/v = А Цлкоэффициент заполнения; - произведение барических i i=поправок для всех j - перепусков. Геометрическая константа N определяется специальной тарировкой. Ее получают, используя предельный переход lim0 p / p = (VI + VII )0 /(VI )0 = N, (2) j-1 j Pj экстраполяцией полинома в виде значений pj-1/pj как функции от (pj), коэффициенты которого определяются процедурой по минимизации среднеквадратичных отклонений.
Коэффициент заполнения A = p0/Z0 определяется аналогичным способом при использовании соотношения:
j j (p0 / Z0) = lim pj N, (3) i Pj i=но его величина для каждого опыта имеет свое значение.
Точность определения N и А зависит от точности измерения давлений и от достаточно большого количества точек на заключительной стадии опыта при низких давлениях.
В тарировочных опытах для определения N использовался гелий-4 высокой чистоты (99,985%) и азот особой чистоты (99,997%).
Надежность опытных данных по He и N2 оценивалась по отклонению экспериментальных значений удельных объемов от стандартных справочных данных (максимальные отклонения были до 0,1%), а также - по вторым вириальным коэффициентам, полученным при экстраполяции зависимости (Zj-1)/j = (j) при j 0.
В результате статистической обработки величин, имеющих нормальное распределение, в интервале с доверительной вероятностью Р = 95% значение N=1,407251 ( 0,002%).
( ( ( Состав газовых смесей в зависимости от способа их приготовления определялся: 1) в случае подготовки смеси непосредственно в ячейке с помощью уравнений состояния индивидуальных компонентов; 2) путем поочередного взвешивания заправочных баллончиков, когда уравнения состояния для какоголибо компонента отсутствовали; 3) с помощью газового хроматографа ЛХМ-80/при использовании готовой газовой смеси.
Дополнительным достоинством первого способа является его высокая информативность: измеряя термические эффекты смешения (pЕ и vE), можно сразу (до окончания основной серии перепусков) определить Vсм, Zсм при Т = Топ и p = pсм, коэффициент заполнения А0 = 0RT = p0 / Z0.
Основные измерения в широком диапазоне температур и давлений проводили изохорически связанным методом последовательных расширений, причем фиксировали pТ - параметры равновесных состояний, как в однофазной, так и в двухфазной, областях однокомпонентного вещества или смеси.
Схема проведения изохорически связанного эксперимента поясняется на рис. 2, применительно к бинарной смеси постоянного состава. Переход смеси из однофазного состояния в двухфазное устанавливался по излому изохор.
Смесь заданного состава в количестве М0 молей с максимальным давлением ( p01) готовили в полном объеме пьезометрической ячейки (VI+VII) из компонентов, находящихся первоначально в разных камерах VI и VII. Наивысшая температура Т(1) выбиралась достаточно далеко от критической, где адсорбционные эффекты практически отсутствуют и, следовательно, не могут исказить перераспределение массы газа между двумя камерами в тот момент, когда происх двумя роисходит их разделение (перед очередным д перепуском уском).
Температур термостата ературу понижали соответственно до жали Т(2), Т(3), Т(4), Т(5),ЕT(k) и, измеряли равновесные яли давления ( ( ( ( ( p02), p03), p04), p05),...p0k ) при, изомолярно переходе (М0 = олярном const) на первой квазиизохоре ) на п ( ( v0k ) = v01) (T, p), где (T, p) - полином третьей степени, ном коэффицие которого фициенты учитывают термическую и вают барическую деформации ескую измеритель ячейки, а k - ительной указывает н вает номер изотермы (k = 2, 3, 4, 5, Е) после каждого 4, 5, изохоричес рического перемещения опытного о ного образца.
Рис. 2 - Иллюстрация к измерениям p,Т - параметров смеси постоянного ация тров состава в газовой и двухфазной областях при изомоляр переходах по и д олярных квазиизохорам и послед азовой изотерме Т(1) оследовательных расширениях на базово Затем пьезометри азовой метрическую систему нагревали до базово температуры Т(1), ( восстанавливая первона осле рвоначальное состояние p01), и, после удаления газа из пьезометра VII, производ изотермическое расширение оставшегося количества изводили е оста ( М1 газовой смеси того же состава до давления p11). Далее измеряли равновесные того ж ее изм ( ( (4 ( ( давления p02), p03), p04), p05),...p0k ) при температурах Т(2), Т(3), Т(4), Т(5),ЕT(k) на Т второй квазиизохоре (М1 = const) и т.д.
М и расч изохорам на термическую Поправка при расчете удельных объемов по квазиизохор деформацию пьезометро пределялась из обработки метров для каждого i-ого перепуска опреде представительного масси измерений на гелии как веществ имеющем статус массива ществе, стандартных справочных данных, в изохорически связанном опыте в следующей очных нном форме:
(1) (1) (1) vi(k ) = vi(1)[1+ a1k (Tk - T ) + a2k (Tk - T )2 + a3k (Tk - T )3] T, (4) (1) Tk пературы хорическом перемещении где T, - температ на базовой и при k-том изохорич изотермах; а индекс i соо соответствует номеру перепуска.
и изме солютная погрешность T Погрешности измеряемых величин: предельная абсолют в измерении температур предел ратуры равна 0,03 К; максимальная предельная относительная погрешность измерения давления p = p/p = 0,05%. Относите ения осительные погрешности величин, определяемых при реализации изохорически связанного метода яемых ски последовательных расширений опытного образца: коэффици заполнения расш ффициента вов имаемости и удельных А 0,03% ; составов Y = 0,22% ; коэффициентов сжимаем объемов: Z = 0,15%, v = 0,16%.
v На данной установке исследованы pvT свойства: 1) декафторбутана (R31-10) в диапазоне давлений 0,6Е1 МПа и температур 279,15Е353,15 К; 2) гексафторида серы (R846) и его смесей с октафторпропаном (R218) в диапазоне температур Т = 238,15Е333,15 К при давлениях до 3,5 МПа, в том числе, азеотропной смеси R510;
3) гептафторбутанолового эфира HFE347mcc и его смесей (по три состава каждой композиции) с тетрафторметаном (R14), октафторпропаном (R218) и трифторметаном (R23) в диапазоне температур 213,15Е353,15 К и давлений до 2,90 МПа.
pT Экспериментальная установка л и метод пьезометра постоянного объема для измерения плотности вещества в жидком и газообразном состоянии реализованы в следующей конфигурации (рис. 3). Система создания и измерения давления включает в себя гидравлические прессы с грузопоршневыми манометрами МП-60 класса 0,02 и МП-6 класса 0,02 (поз. 9, 10), заправленными трансформаторным маслом, соединенными через масло-масляный разделитель (11) с грузопоршневым манометром МП-600 класса 0,02. Давление, создаваемое гидропрессами, уравновешивает давление в пьезометре (2) через мембранный разделитель (12). Состояние мембранного разделителя регистрируется с помощью дифференциального трансформаторного датчика (13) со вторичным прибором (14).
Полость разделителя над мембраной, заполненная исследуемым веществом, соединена с пьезометром тонким капилляром (16). Заполнение линий исследуемым веществом и маслом, передающим давление, осуществляется после предварительного вакуумирования форвакуумным насосом ВН-461 (19).
pT Рис. 3 - Схема экспериментальной установки л Система измерения и регулирования температуры включает в себя жидкостный термостат (1), изготовленный по индивидуальному проекту, обеспечивающий в комплекте с дополняющим его оборудованием поддержание температуры с устойчивостью не менее 0,01 С в интервале температур 30 - 550 С при условии правильно подобранного теплоносителя. Теплоноситель перемешивается по кольцевой схеме при помощи электроприводной мешалки (6), установленной на охлаждаемой муфте (5) с ременным приводом от электродвигателя (7). Подвод необходимой постоянной тепловой мощности достигается основным электронагревателем (3), питание которого обеспечивается точно установленным стабилизированным напряжением системы электропитания.
Поддержание температуры в термостате осуществляется при помощи автоматического регулятора температуры, состоящего из регулирующего блока Р111 в комплекте с БУТ-1 (неполный комплект ВРТ-3) и блока, задающего регулируемую температуру, работающего на базе алмазного датчика температуры.
Измерение температуры осуществляется платиновым десятиомным термометром сопротивления ПТС-10, работающим в комплекте с полуавтоматическим потенциометром Р-363/1 класса 0,001, которые обеспечивают измерение температуры в интервале 0-630 С с погрешностью не более 0,01 С.
Система электропитания подключена в трехфазную сеть напряжением 380/220 В с установленной максимальной потребляемой мощностью 6 кВт. В состав системы входит пульт пакетных выключателей и переключателей, обеспечивающий коммутацию нагревателей, агрегатов и приборов.
К системе измерения плотности относятся: пьезометр постоянного объема (2), выполненный из сплава ЭИ-617 в виде толстостенного сосуда шарообразной формы; заправочный баллон с исследуемым веществом; узел для отбора вещества (15); мерный баллон оригинальной конструкции и аналитические весы (17) ВЛА200 с ценой деления шкалы 0,1 мг.
Определение объема пьезометра. Точность экспериментальных данных по плотности пропорционально зависит от точности определения объема пьезометра.
Тарировочный опыт основывался на использовании в качестве эталонного вещества диоксида углерода особой чистоты (99,997 %) и проводился в изотермическом режиме при температуре 32,61 С. Измерения проводились при давлениях 5, 6, 8 и 10 МПа.
Объем пьезометра в каждом состоянии находился, как функция от массы вещества в пьезометре и плотности в зависимости от давления и температуры с учетом погрешностей от внесения поправок и погрешностей отнесения:
[ m ] - m вып бал, (5) V = ( p, t ) CO где - коэффициент приведения объема к начальным условиям (t = 20 C, р = 0,1 МПа); mвып - суммарная масса выпущенного диоксида углерода от данного состояния до последнего выпуска, mбал - масса диоксида углерода в балластном объеме в данном состоянии, CO2(p,T) - плотность диоксида углерода при параметрах данного состояния. Числитель дроби в представленном выражении соответствует массе вещества в пьезометре при данных условиях опыта.
Усредненное значение объема пьезометра с учетом округлений составило 204,530,01 см3.
Методика проведения опытов. В ходе опыта, который проводится по изотермам (это выгодно при изменении состояний) и начинается с максимального давления, производятся частичные выпуски вещества в баллон известной массы с последующим взвешиванием его на аналитических весах. При каждом выпуске давление в пьезометре понижается и измеряется после возвращения системы на стационарный температурный режим. Таким образом, в конце опыта известна вся масса вещества в каждом состоянии:
n k m = m + m - m (6) i 0 i k 1 , где mk - масса вещества в данном состоянии k, считая первым - состояние при максимальном давлении опыта; i - номер выпуска, mi - масса выпущенного вещества при переходе в соответствующее состояние, n - количество всех выпусков, m0 - остаточная масса вещества в пьезометре, рассчитанная по уравнению состояния.
Учет термической деформации производится следующим образом:
V = V20 [1 + 3(t - 20 ) + 3 (t - 20 )2 ], (7) где V - объем пьезометра при температуре опыта, V20 - объем пьезометра при С; - коэффициент линейного расширения при соответствующей температуре.
Известное соотношение для барической деформации сферического пьезометра корректировалось специальной тарировкой и, таким образом, было учтено наличие сварных швов.
Погрешность определения объема пьезометра была оценена в 0,034%.
Погрешность полученных значений плотности в газовой фазе, без учета погрешностей отнесения, оценена в 0,16Ц0,25%. Погрешность измерения плотности в жидкой фазе при измерении по изохорам - не более 0,13%.
pT На данной установке исследованы свойства: октафторпропана (R218) в диапазоне температур 373,15Е773,15 К и давлений 1Е10 МПа в газовой фазе; декафторбутана (R31-10) в диапазоне температур 293,15Е773,15 К и давлений 1Е10 МПа, как в газовой фазе, так и в жидкой, а также на линиях насыщения.
2 РАСЧЕТНЫЕ И ТЕОРЕТИЧЕСКИЕ МЕТОДЫ ИССЛЕДОВАНИЯ Методы построения уравнений состояния для расчета термодинамических свойств индивидуальных веществ, смесей и растворов.
В настоящей работе полиномиальная зависимость коэффициента сжимаемости для индивидуальных веществ формировалась в виде:
Si r - j z = 1+ i (8) bij i=1 j=где z = pv / RT - коэффициент сжимаемости; = / кр - приведенная плотность;
= Т /Ткр - приведенная температура. Для получения коэффициентов bij уравнения использовался метод наименьших квадратов на основе минимизации функционала:
Si n r - j S = z -1- i, (9) Wi bij k =1 i=1 j= Wi = 1/( zэ )2, (10) k где Wi - весовая функция; k- номер экспериментальной точки; к - предельная относительная погрешность экспериментальных значений сжимаемости.
Для описания использовалась компьютерная программа, написанная на языке FORTRAN в пакете компьютерных программ FORTRAN POWER STATION 4.0. Для массива экспериментальных рT-данных величина среднеквадратического отклонения аппроксимации вычислялась по соотношению:
n p э = 100 (1- z / z ) /(n - m) (11) z 1 где zэ, zр - соответственно, экспериментальное и расчетное значение сжимаемости в каждой точке; n - число экспериментальных точек; m - число эмпирических констант уравнения.
Для данного варианта термического уравнения состояния (9) при расчете термодинамических свойств использовалась методика, подробно описанная в работах. Для удобства расчета свойств были использованы комплексы:
Si Si r Si r r i bi, j i i i, j i, j A1 = A0 = A2 = ; ; ;
(i +1)b (i -1)b j j j i=1 j=0 i=1 j=0 i=1 j=Si Si r r Si r bi, j i bi, j i bi, j i A3 = ;. (12) (i -1) A5 =-(i -1) j (i +1) ; A4 = j j j i i i i=1 j=0 i=1 j=0 i=1 j=В данных комплексах символы i, j, , аналогичны символам в формуле (9).
Из соотношений (12) были найдены термодинамические функции в приведенной h h= + Aбезразмерной форме: сжимаемость Z = 1 + A0 ; энтальпия h*= ;
RT RT s s0 с c = - ln + A4 ; изохорная теплоемкость cv*= = + A5 ;
энтропия s* = R R 0 R R сp cp (1+ A2) = +(A5 -1)+ изобарная теплоемкость cp* = ; скорость звука w* = R R (1+ A1) w (A2 - A1) = 1+ A1 ; адиабатный дроссель-эффект * = =.
w0 0 (1+ A1) Для расчета калорических свойств были использованы идеально-газовые функции сp0, h0, S0, представленные в соответствующих литературных источниках.
Для смесевых композиций зависимость коэффициента сжимаемости формировалась в виде квазивириального уравнения состояния в широких рТ - диапазонах газовой фазы:
ni m l i j k Zm = 1+ ( / )y(13) bijk i=1 j=1 k =p T где Zm = ; = ; =.
mRT N TN В качестве N и TN могут быть выбраны критические параметры либо другие численные величины. Индекс m соответствует параметрам смеси.
Уравнение (13) предназначено для полного концентрационного интервала и может быть названо унитарным для компонентов и смесей произвольного состава.
При недостаточном количестве экспериментальной информации в некоторых областях состояний в качестве базовой модели уравнения состояния, предназначенного для оперативных расчетов термодинамических свойств сосуществующих фаз жидкость - газ смесей низкокипящих веществ, использовалось трёхпараметрическое кубическое уравнение состояния ПателаRT a(T ) p = Теджа:. (14) v - b v(v - b) + c(v - b) Оптимальные значения регулируемых параметров в формулах для аi, bi, ci определяли на основании машинной обработки табличных рvT - данных для компонентов, а температурную зависимость коэффициента бинарного взаимодействия 12 = (A + ВТ) - по собственным рvTx(y)-данным для сосуществующих фаз жидкость - газ смесей.
Все рассмотренные выше уравнения состояния являются экспериментально обоснованными, широкодиапазонными, и применимы в полном концентрационном интервале (x1 = 0 - 100% мольн.).
Термодинамические эффекты смешения являются мерой неидеальности растворов, их измерение позволяет методически независимо определять pvTx - поверхность газовой смеси. В настоящей работе они использовались для подтверждения внутренней согласованности данных, полученных, в основном, методом последовательгых расширений.
Для определения объемного эффекта смешения газов vЕ использовали экспериментальную установку с двухкамерной пьезометрической ячейкой с предварительно измеренным отношением объемов пьезометров N = (VI +VII ) VI = N (1+ aP) и интенсивным циркуляционным перемешиванием E p исследуемого вещества. При определении в начальной стадии изотермического опыта чистые компоненты находятся в разных пьезометрах при одинаковых с погрешностью манометрического прибора давлениях pI = pII = p0.
После открытия вентилей постоянного объема, соединяющих пьезометры VI и VII, включали циркуляционный насос и производили перемешивание до полной стабилизации давления pm образующейся смеси. Зафиксированный прирост давления соответствует pE = (pт - pо), а мольная доля компонента, находящегося, например, в меньшем пьезометре VII, вычисляется по формуле yII = ZII /[ZII + ZI (N -1)], (15) где ZI и ZII - коэффициенты сжимаемости компонентов в VI и VII при {pо;Т}.
Значения vЕ и vm вычислялись по соотношениям:
vE ( p,T, y) = yk [vok ( p - pE,T ) - vok ( p,T )], vE ( p,T, y) = yk[vok ( p - pE,T ) - vok ( p,T )] (16) Необходимым условием реализации такого подхода является наличие достоверных уравнений состояния компонентов смеси.
3 ИССЛЕДОВАНИЕ ТЕРМОДИНАМИЧЕСКИХ СВОЙСТВ ОКТАФТОРПРОПАНА (C3F8) И ДЕКАФТОРБУТАНА (C4F10) Обзор литературных данных о теплофизических свойствах октафторпропана и декафторбутана выполнен до 1947г.
На рис. 4 новые данные по C3F8 МЭИ отмечены знаком (); данные Brown J.A., полученные в 1963г. - (); данные, полученные В.П. Барышевым в 1980г. - ().
Данные по C4F10 () получены в 1947 г. (R.D. Fowler, J.M. Hamilton, J.S. Kasper) и в 1954г. (Brown J.A., Mears W.H.). Всего получено 65 опытных точек в интервале температур от 373 до 773 К и давлениях до 10 МПа. Новые данные по C4F10 МЭИ получены при повышенных температурах (до 500 С) и давлениях до 10 МПа на установке, реализующей метод пьезометра постоянного объема; опытные точки представлены на рисунке 4 (б) - с условным обозначением л (всего 65 точек).
а) б) Рис. 4 - Диап ений плотности:
Диапазон известных и новых измерений а) октафто тана (C4F10) тафторпропана (C3F8); б) декафторбутана Также получено 27 точек в области жидкости на этой же установке и ено 2 на эт рассчитано 4 точки с определением на левой пограни кривой. При ки с ограничной умеренных температурах вблизи кривой насыщения на установке, реализующей атурах устан метод последовательных расширений, было получено 55 точек, которые условно льных точе обозначены л.
Погрешность полученных экспериментальных данных оценивается в пол данн пределах 0,2% без учета погрешностей отнесения для данных, полученных з уче ля да методом пьезометра постоянного объема, и 0,15% - методом последовательных а пос етодом расширений.Максимальн погрешности измерений с учетом вероятных ошибок мальные четом отнесения в надкрити 0,критической области не превышают 0,45 %. Однако такие состояния не входят в область предполагаемых расчетов. В диапазоне расчетов ят в о В теплосиловых циклов полная погрешность измерения плотнос по отности октафторпропана и декафторбутана состав оставляет 0,19Е0,27%.
Разработанное уравнение состояния декафторбутана вириального типа (9) е ура на ви содержит 26 коэффиц исх ффициентов и описывает массив исходных данных со среднеквадратическим отклонением sk = 0,24%. С целью повышения устойчивости ким от повыш данного уравнения на первом этапе были разработаны 25 эквивалентных я на таны уравнений и проверен их внутренняя согласованность при сопоставлении верена ость получаемых калорически величин. На втором этапе было определено уравнение ических ло оп состояния с усредненным а процедура оптимизации енными коэффициентами и применена проц числа этих коэффициенто дены циентов, их численные значения приведены в табл. 1.
Табл. 1 - Коэффициенты ояния декафторбутана в иенты усредненного уравнения состояния газовой фазе b10.895517138+0 b22.305204027+01 b51.877806277+138+b11 -.451977419+0 b33 -.210191305+03 b52 -.257695651+419+b12 -.725628973+0 b34.528515254+03 b23.272666551+973+b13.151229654+0 b40 -.242956271+03 b30.262613803+654+b14 -.113766528+0 b42 -.132395234+03 b31.119690376+528+b15.294819457+0 b43.214308648+03 b53.220203588+457+b16 -.727144512+0 b44.235303664+03 b54 -.167702407+512+b20.209679895+0 b45 -.588753006+03 b55 -.333616999+895+b21 -.169746079+0 b50.270925138+079+Разработано уравнение состояния октафторпропана вириального типа (9) с уравн а вир 27 коэффициентами (табл. 2), существенно расширяющее расчетный диапазон по (таб е рас давлениям до 12 МПа и температурам до 823,15 К. Среднее квадратическое отклонение описания опытных данных составило 0,24 %.
Результаты тестирования уравнения состояния октафторпропана при сопоставлении расчетных калорических величин и соответствующих прямых экспериментальных измерений подтвердили его высокую достоверность.
Например, основные расхождения данных по изобарной теплоемкости в интервале 4% температур 320Е420 К находятся в пределах, и в большинстве точек имеют разный знак.
Табл. 2 - Коэффициенты уравнения состояния октофторпропана b10 10,0991666 b23 487,8757060 b41 101,61997b11 -69,1750574 b24 -272,6943450 b42 -238,39646b12 188,5019770 b25 54,6259367 b43 234,83000b13 -252,6102350 b30 28,9916996 b44 -84,14328b14 162,5383900 b31 -172,1530880 b50 -8,88306b15 -40,4872052 b32 376,6919530 b51 29,56638b20 -23,1360654 b33 -359,6544650 b52 -30,92879b21 157,2898870 b34 127,1500710 b53 8,07896b22 -403,9661080 b40 -14,7356347 b54 2,38218Рассчитаны таблицы основных термодинамических свойств октафторпропана и декафторбутана в газовой фазе (коэффициент сжимаемости z, удельные энтальпия h, энтропия s, изобарная теплоемкость cp, изохорная cv и скорость звука W) в диапазоне давлений 0,1Е10 МПа и температур 293,15Е823,15 К.
4 ИССЛЕДОВАНИЕ ТЕРМОДИНАМИЧЕСКИХ СВОЙСТВ ГЕПТПАФТОРБУТАНОЛОВОГО ЭФИРА HFE347mcc И ЕГО СМЕСЕЙ С ОКТАФТОРПРОПАНОМ FC218, ТРИФТОРМЕТАНОМ HFC23, ТЕТРАФТОРМЕТАНОМ FC Сделан подробный обзор литературных данных о термодинамических свойствах компонентов смесей на основе гептафторбутанолового эфира. Данных о каких-либо измерениях для смесевых композиций на основе фторэфира не обнаружено.
Фторэфир HFE347mcc (CF3CF2CF2OCH3) принадлежит к новому поколению озонобезопасных хладагентов с малыми значениями экологических потенциалов. У этого вещества время жизни LT[y] = 5.6, GWP[100y] = 368-430, HGWP меньше 0.15, и его рассматривают в качестве возможного заменителя фтортрихлорметана (СFCl3; R11). При изучении смесей основными компонентами явились трифторметан (СHF3; HFC23), октафторпропан (C3F8; FC218) и тетрафторметан (C1F4; FC14), которые также относятся к озонобезопасным веществам с ODP = 0, позволяют повышать давление насыщения для смесевого хладагента и подбирать конфигурацию кривых кипения и конденсации для каждого варианта предполагаемой замены.
Измерение давления насыщения фторэфира. Для фторэфира было проведено три серии опытов. В первой серии значения рs измеряли в области состояний жидкость-газ с содержанием основного вещества не менее 99,9 масс.
% в диапазоне температур Т = 313,15Е353,15 К (8 точек). На основании совместной обработки трёх групп данных в интервале Т = 260Е437,7 К получена i аппроксимационная зависимость для давления насыщения вида: ln ps = T, ai i=-a где рs в барах, Т - в Кельвинах, а коэффициенты равны:
i а-1 = -1,7470332105; а0 = 2,8749816103; а1 = -1,9897807101; а2 = 7,348913710-2;
а5 = -7,609130010-11; а3 = -1,520712410-4; а4 = 1,670652110-7.
Определение pvT- данных в паровой фазе и на кривой насыщения фторэфира HFE347mcc.
Во второй и третьей сериях были получены экспериментальные pvT-данные для чистого фторэфира HFE347mcc на двух базовых изотермах и пяти квазиизохорах в диапазоне температур от 333,15 до 353,15 К и давлений от 105,до 372,56 кПа. По опытным рТ-данным рассчитаны параметры состояния в точках росы, а также значения мольных объёмов и коэффициентов сжимаемости для перегретого пара и двухфазной области. Всего 45 опытных точек. Новые данные существенно дополнили известные к настоящему моменту измерения зарубежных исследователей, а проведенный анализ показал их согласованность.
Для кривой насыщения построено уравнение состояния кубического типа Карнахана - Старлинга - де Сантиса:
p 1+ y + y2 - y3 A = RT RT( + B) (1- y) (17) с использованием, как собственных данных на кривой конденсации (шесть точек), так и данных японских исследователей (четыре точки на кривой конденсации и точек на кривой кипения). Коэффициенты А(Т) и В(Т) уравнения состояния Карнахана - Старлинга - де Сантиса эфира HFE-347 mcc (18) определялись для каждой экспериментальной точки при условии равенства химических потенциалов на кривой кипения и конденсации. Величины А(Т) и В(Т) аппроксимированы полиномами:
A(T ) = a0 exp(a1T + a2T ) B(T ) = b0 + b1T + b2T, (18) где a0 =.52392795+04; b0 =.24838381+00;
a1 =.48038101-03; b1 = -.11139012-03;
a2 = -.56315361-05; b2 = -33576770-06.
Среднеквадратическая погрешность определения А(Т) составила 0,061%, а В(Т) - 0,053%.
Для расчета термодинамических величин в газовой фазе было разработано уравнение состояния эфира HFE-347mcc вириального типа (9) с коэффициентами (табл. 3), которые были определены на основе собственных данных и данных японских исследователей.
Табл. 3 - Коэффициенты уравнения состояния фторэфира HFE 347mcc b10 = -.597032449D+03 b22 = -.189548938D+02 b44 = -.906546099D+b11 =.269420644D+04 b24 =.507531473D+01 b50 = -.168145672D+b12 = -.486147701D+04 b30 = -.311913653D+02 b51 =.454986408D+b13 =.437274212D+04 b31 =.161725316D+02 b53 =.111076841D+b14 = -.195859323D+04 b33 =.132357599D+02 b60 =.418988762D+b 15 =.348654586D+03 b40 =.271768963D+02 b62 = -.393046881D+b20 =.153403744D+02 b42 = -.161357085D+В программу нахождения коэффициентов bij уравнения состояния была заложена возможность произвольного набора матрицы индексов суммирования по плотности и температуре. Среднеквадратическая погрешность аппроксимации составила 0,32 %.
Данные по идеально-газовой теплоемкости фторэфира в литературе отсутствовали и определялись с помощью методики, учитывающей вклад структурных групп. Такая методика была апробирована на более изученных хладонах. С ее помощью было разработано соотношение для вычисления идеально-газовой теплоемкости эфира HFE-347mcc:
Cp0 / R = -1,0243 -1 +13, 6432 0 +45,23 1 -23,537 2, (19) где = T/1000.
Результаты сравнения значений Cp0/R, рассчитанных в приближении жесткий ротатор гармонический осциллятор с использованием имевшихся спектроскопических данных, показали, что согласие сравниваемых данных удовлетворительно в интервале температур 200Е800 К, где расхождения не превышают 1-2 %. Указанная величина была принята в качестве оценки погрешности расчетных значений теплоемкости Cp0.
Исследование термодинамических свойств смесей фторэфира с октафторпропаном, трифтрометаном и тетрафторметаном.
pvT -данные для бинарных смесей фторэфира Экспериментальные HFE347mcc с хладонами FС218, HFC23 и FC14 в газовой фазе и на кривых конденсации были получены впервые.
Экспериментальные данные для бинарной системы FC218 - HFE347mcc получены для трех составов: 5; 10,04 и 20,64 масс. % второго компонента и охватывают диапазон температур от 299,15 до 353,15 К и давлений от 101 до 20кПа. Всего получено 53 опытные точки в области перегретого пара, 14 - в двухфазной области состояний и 5 точек на кривой конденсации.
Экспериментальные данные для бинарной системы R23 - HFE 347mcc получены:
- на пяти квазиизохорах в диапазоне температур 288,15Е353,15 К и при давлениях до 2,904 МПа для смесей трех составов с содержанием HFE 347mcc 12,22; 16,40 и 30,77 масс. %;
- на базовой изотерме 353,15 К в интервале давлений 0,1302Е 1,291 МПа для смеси, содержащей 30,77 масс. % HFE 347mcc;
- на базовой изотерме 333,15 К в интервале давлений 0,109Е5,099 МПа для смеси, содержащей 9,40 масс. % HFE 347mcc;
- на базовой изотерме 313,15 К в интервале давлений 0,1385Е2,421 МПа для смеси, содержащей 12,22 масс. % HFE 347mcc. Всего 82 опытных точки.
Экспериментальные данные для бинарной системы хладагентов R14 - HFE 347mcc получены:
- на трех квазиизохорах в диапазоне температур 213,15Е353,15 К и при давлениях до 2,18 МПа для смесей трех составов с содержанием HFE 347mcc 35,90; 18,70 и 9,40 масс. %;
- на базовой изотерме 353,15 K при двух составах смеси 37,80 и 35,90 масс.
% HFE 347mcc;
- на базовой изотерме 313,15 К в интервале давлений 0,109 Е5,099 МПа для смеси, содержащей 9,40 масс. % HFE 347mcc. Всего 46 опытных точек.
Предельные погрешности измерения температуры и давления в ячейке равновесия оценивались в пределах: Т = 0,03 K; Р = 0,05 %.
Индивидуальные компоненты исследуемых систем содержали не менее 99,90% основного продукта, и, следовательно, только данный фактор всецело определял погрешность состава в пределах 0,2%. Все данные были получены на двухкамерной установке pvTx методом изохорически связанных последовательных расширений опытного образца. Погрешность полученных значений коэффициента сжимаемости (или мольных объемов) без учета ошибок отнесения оценивалась величиной не более 0,15%.
Коэффициенты для трех систем приведены в табл.4-6.
Табл. 4 - Коэффициенты УС бинарной системы FC218-HFE347mcc b100 =-.919537895671030D+03 b200 =-588183911641896D+03 b232 =-.810544161295714D+b101 =.208999714446658D+03 b201 = 112707957133896D+04 b233 =.353709591019647D+b102 =.173084230231090D+04 b202 =-463550628464317D+03 b300 =-.128445852703237D+b103 =.642229854264370D+03 b203 = 262869330135201D+03 b301 =.398123556949017D+b104 =-.169852567778283D+04 b210 = 659647002735874D+03 b303 =-.123736262186052D+b110 =.270118092775896D+04 b211 =-164485390074229D+04 b304 =-.182614491855613D+b111 =-.193628905987250D+04 b212 = 150395358539937D+03 b310 =-.448192173847199D+b112 =-.331250520794357D+04 b213 =-355407260488766D+03 b311 =.400313450414809D+b113 =-.792485772158339D+03 b220 = 392716339883678D+03 b313 =-.319672038475758D+b114 =.342260363703289D+04 b221 = 656801523415933D+03 b314 =.175850171954202D+b120 =-.166001654227743D+04 b222 =.227048470600594D+03 b330 =.650613456636313D+b121 =.127050236220287D+04 b223 =.771863708457340D+02 b331 =-.967588126218346D+b122 =.219764842555139D+04 b230 =-.696186875936918D+03 b333 =.663276344899703D+b123 =-.194583738839398D+03 b231 =.651611004743811D+03 b334 =-.282532383910413D+b124 =-.166025029881409D+Табл. 5 - Коэффициенты УС бинарной системы HFC23-HFE347mcc b100 =.213744549D+02 b200 =-.252878916D+04 b232 =.185694165D+b101 =-.202330452D+03 b201 =.306575685D+04 b233 = -.172073401D+b102 =.500334154D+03 b202 = -.824399163D+03 b300 =.441385079D+b103 =-.596037055D+03 b203 =.327835914D+03 b303 =.507599266D+b104 =.281096202D+03 b210 =.502046082D+04 b305 =-.917580410D+b110 =.523953678D+03 b211 =-.171843238D+04 b320 = -.618198710D+b111 =-.890115423D+03 b212 =-.354808179D+04 b321 =.118306244D+b112 =-.598835483D+03 b213 =.838175726D+02 b323 =-.638314291D+b113 =.200600807D+04 b220 =-.465919327D+04 b325 =.558264806D+b114 =-.105113696D+04 b221 =.107779256D+04 b340 =.180545070D+b121 =.103465806D+04 b222 =.245573040D+04 b341 =-.341969499D+b122 =.181750364D+03 b223 =.133111877D+04 b343 =.137050193D+b123 =-.146219121D+04 b230 =.215131703D+04 b345 =.336309860D+b124 =.781864867D+03 b231 =-.237112416D+Уравнения состояния (УС) смесевых хладагентов на основе фторэфира HFE347mcc строились на массивах собственных экспериментальных pvTx - данных и расчетных значениях плотности фторуглеродных компонентов по справочным данным при параметрах смеси в форме вириального разложения (13). Уравнения состояния воспроизводят исходные экспериментальные pvTx- данные со среднеквадратической погрешностью 0,109, а гистограммы отклонений близки к нормальному распределению.
Табл. 6 - Коэффициенты УС бинарной системы FC14-HFE347mcc b100 =-.919537895671030D+03 b200 =-588183911641896D+03 b233 =.353709591019647D+b101 =.208999714446658D+03 b201 = 112707957133896D+04 b300 =-.128445852703237D+b102 =.173084230231090D+04 b202 =-463550628464317D+03 b301 =.398123556949017D+b103 =.642229854264370D+03 b203 =.262869330135201D+03 b303 =-.123736262186052D+b104 =-.169852567778283D+04 b210 = 659647002735874D+03 b304 =-.182614491855613D+b110 =.270118092775896D+04 b211 =-164485390074229D+04 b310 =-.448192173847199D+b111 =-.193628905987250D+04 b212 =.150395358539937D+03 b311 =.400313450414809D+b112 =-.331250520794357D+04 b213 =-355407260488766D+03 b313 =-.319672038475758D+b113 =-.792485772158339D+03 b220 = 392716339883678D+03 b314 =.175850171954202D+b114 =.342260363703289D+04 b221 =.656801523415933D+03 b330 =.650613456636313D+b120 =-.166001654227743D+04 b222 =.227048470600594D+03 b331 =-.967588126218346D+b121 =.127050236220287D+04 b223 =.771863708457340D+02 b333 =.663276344899703D+b122 =.219764842555139D+04 b230 =-.696186875936918D+03 b334 =-.282532383910413D+b123 =-.194583738839398D+03 b231 =.651611004743811D+03 b233 =.353709591019647D+b124 =-.166025029881409D+04 b232 =-.810544161295714D+ Уравнения состояния трех бинарных систем описывают термодинамические свойства в диапазоне температур 303,15Е353,15 К и концентраций первого (основного) компонента в диапазоне 0,85Е1,0 мольных долей. По разработанным уравнениям были рассчитаны таблицы термодинамических свойств смесевых хладагентов в газовой фазе.
5 ИССЛЕДОВАНИЕ ТЕРМОДИНАМИЧЕСКИХ СВОЙСТВ СМЕСЕЙ ГЕКСАФТОРИДА СЕРЫ С ОКТАФТОРПРОПАНОМ (SF6/C3F8) В обзоре отмечено, что к началу настоящего исследования имелась обширная информация из литературных источников о теплофизических свойствах (ТФС) компонентов данной смесевой композиции. В Московском энергетическом институте (на кафедрах Теоретических основ теплотехники и Инженерной теплофизики) была выполнена крупномасштабная программа экспериментальных исследований ТФС SF6, в первую очередь, в газовой и жидкой фазах (1970-19гг.), а также его смесей с низкокипящими веществами (Не, N2, CF4, CHF3,Е) в период 1982-1995 гг.
Экспериментальные pvTx- данные для системы SF6-C3F8. Для компонентов системы SF6-C3F8 имеются достаточно надёжные уравнения состояния. В этом случае, исследуемую смесь готовили в полном объёме ячейки (V1+V2) в процессе изотермического смешивания компонентов, находившихся первоначально в разных пьезометрах при Т = Т0 и одинаковых давлениях (р1 = р2).
Перемешивание осуществляли до полной стабилизации давления в ячейке на значении р0 и определяли барический эффект изохорно-изотермического смешения компонентов рЕ = р0 - р1,2; затем вычисляли v0, а также мольную долю компонента, находившегося первоначально в объёме V1, по формулам:
v0( p0,T0, x) = xk v0k ( p0 - pE,T ) (20) k=x1 =, (21) v1 (N -1) 1 + vгде v1, v2, v0k - мольные объёмы компонентов при Р1,2 и Т0.
Данная процедура для системы SF6-C3F8 была осуществлена на базовой изотерме Тб = 333,15 К, при этом давление и объём в начальной точке составили: Р= 3432,3 кПа, v0 = 0,43563 м3/кмоль, а содержание второго компонента x2 = 17,мольн. %.
Основные опыты проводили изохорически связанным методом последовательных расширений в широком температурном интервале и фиксировали рТ- параметры, как в однофазной, так и в двухфазной, областях (рис.
2). Значения pvT-параметров на кривой конденсации были получены пересечением линий полиномов (решением системы двух уравнений), аппроксимирующих экспериментальные точки на квазиизохорах в газовой и двухфазной областях.
Новые экспериментальные pvT- данные в паровой фазе и на кривой конденсации системы SF6 - C3F8 с содержанием второго компонента x2 = 17,мольн. % получены в диапазоне температур Т = 288,15Е333,15 К при давлениях до 3,5 МПа. Погрешность мольного объёма, с учётом ошибок отнесения не превышает 0,16 %.
Новые экспериментальные pvT- данные в паровой фазе системы SF6 - C3F8 с содержанием второго компонента x2 = 93,82 мольн. %, (или 0,05 масс. % SF6 - состав азеотропного смесевого хладагента R510), получены в диапазоне температур Т = 238,15Е353,15 К при давлениях до 3,3 МПа. На базовой изотерме (353,15 К) получено 14 опытных точек, а на 11 квазиизохорах получено опытных значения.
Результаты описания опытных pvTx- данных кубическим уравнением состояния. Для построения многоконстантных уравнений квазивириального типа (13) требуются более полные массивы опытных данных, чем полученные для бинарной системы SF6-C3F8. Поэтому для описания полученных экспериментальных данных по двум составам бинарной системы применено трёхпараметрическое кубическое уравнение состояния Патела-Теджа (14).
Значения регулируемых параметров компонентов в формулах для аi, bi, ci определены по pvT-данным в диапазоне Т = 233Е473 К при давлениях до 3,5 Е МПа и равны: с1 = 0,3163451; F1 = 0,7302131; с2 = 0,3306856; F2=0,8787626.
Коэффициенты температурной зависимости параметра бинарного взаимодействия ij =А+ВТ найдены по полученным в этой работе экспериментальным pvTxданным, а регулируемые константы самого уравнения равны: a(T)=1,229423; b= 0,496210-3.
Среднеквадратические отклонения мольных объёмов, рассчитанных по уравнению состояния (14), от исходных данных для компонентов (1 и 2) и смеси (mix), равны соответственно: v1 = 0,943%; v2 =0,678%; vmix =0,676%. В сокращённом диапазоне температур (например, Т = 253Е333 К) все отклонения заметно уменьшаются.
6 РАЗРАБОТКА ТЕПЛОСИЛОВЫХ ЦИКЛОВ НА ФТОРОРГАНИЧЕСКИХ РАБОЧИХ ТЕЛАХ Термодинамическая эффективность фторуглеродов в качестве рабочих тел теплосиловых циклов АЭС.
Разработаны цикл и схема теплосилового контура реакторной установки (РУ) на быстрых нейтронах (рис. 5).
I - линия конденсации; II - линия кипения; К - критическая точка; ТР - теплообменник реактора тепловой мощностью Q0; Т - турбина; Г - генератор, вырабатывающий электрическую мощность Nэ; РТ - регенеративный теплообменник тепловой мощностью Qр; К - конденсатор тепловой мощностью Qк;
Н - насос.
Рис. 5 - Принципиальная схема и цикл теплосилового контура РУ на быстрых нейтронах Результаты расчета параметров цикла второго контура АЭС на декафторбутане С4F10 с конденсацией при положительных и отрицательных температурах приведены в таблице 7.
Действительный КПД цикла находился по формуле:
Т Н (h1 -h2)oi -(h5 -h4')/oi i = (22) (h1 -h6) Табл. 7 - Зависимость КПД цикла второго контура АЭС на декафторбутане от давления при различных температурах перед турбиной и конденсации цикла Р1 = 4 МПа Р1 = 7 МПа Р1 = 10 МПа Р1 = 14 МПа Р1 = 18 МПа Температура конденсации t2 = -10 C 0,467 0,491 0,502 0,512 0,5t1 = 520 C 0,479 0,503 0,514 0,523 0,5t1 = 550 C Температура конденсации t2 = 20C 0,411 0,445 0,460 0,472 0,4t1 = 520 C 0,423 0,457 0,473 0,484 0,4t1 = 550 C Для приведенных ниже расчетов использовались следующие оценки: потери давления в теплообменнике реактора 4%; потери давления в регенераторе по верхней изобаре 3%, по нижней - 12%; потери давления в конденсаторе 10%;
относительный внутренний КПД турбины: 92%; КПД насоса: 85%;
недоохлаждение газа в регенеративном теплообменнике: 5 С.
На основе расчетов сверхкритических конденсаци расче нсационных циклов на фторуглеродах был сделан вывод о том, что замена водопаро теплоносителя сдел опарового закритических параметр на фторуглероды во втором контуре энергетической аметров конт установки позволит дост низких давлениях.
т достичь тех же самых КПД при более низки Термодинамическ ароводяного контура на ический анализ схемы замещения парово фторуглеродный в паро парогазовых установках.
В данном подразделе представлена оценка термодинамической под ка эффективности схем ПГУ при замещении водяного контура ПТУ на хем ого фторуглеродный. Показа что из трех проработанных вариантов один (рис. 6) оказано, вариа дает существенное увеличение внутреннего КПД цикла. Кроме того, отмечен ряд увели Кром технологических и конструкционных преимуществ, возможн при внедрении конс зможных фторуглеродов в качеств ых ус ачестве рабочего вещества теплосиловых установок. Разработка новых схем ПГУ была проведена на расчетной базе (экспериментальнобыл базе обоснованные уравнения нения состояния), описанной выше.
Парогазовая установка, представленная на рис. 6, работает по уст регенеративному циклу ГТУ с утилизацией остаточной теплоты уходящих из иклу теп рекуператора газов в одноконтурном паротурбинном цикле на декафторбутане в о икле (C4F10). Расчет такой схемы показал, что внутренний КПД бинарного цикла, без сх ПД б учета потерь на прокачку котле-утилизаторе, равен качку рабочего тела в рекуператоре и котле 74%. Это позволяет предполагать, что, при удачном констру ет пр нструировании бинарной установки, будет возмож КПД озможным достижение электрического КПД нетто, существенно превышающего 60%.
ГТ - газовая турбина, КС на, КС - камера сгорания,РП - рекуператор; КУ - котел-утилизатор; Н - насос; К - кон р; КУ конденсатор; ЭГ - электрогенератор Рис. 6 - Схема и цикл ой ГТУ и утилизацией цикл бинарной ПГУ с регенеративной ГТ остаточной теп декафторбутане й теплоты рекуператора в цикле на декаф Разработка теплосиловых циклов утилизационных установок в системе тепло ых у распределенной энергетики. Сформулированы особен энер собенности технологий генерации электроэнерг в циклах на органических рабочих веществах.
оэнергии их р Проведено сопоставле ивности тавление энергетической эффективност ряда широко используемых в утилизационных установках рабочих тел (углеводородов) и утилиз тел предлагаемых к внедрени едрению фторуглеродов (рис.7).
Рис.7 - Зависимость внутреннего КПД цикла (i ) от давления пара на входе в турбину (Р1) для ряда рабочих веществ утилизационных установок На данном температурном уровне утилизации теплоты (температуры: пара на входе в турбину - 290 С, конденсации - 40 С) наибольший внутренний КПД цикла был получен на гексане. Но, очевидно, что такие пожароопасные рабочие вещества как пентан, бутан, гексан, могут работать только в двухконтурных схемах с промежуточным теплообменником, через котел-утилизатор которых циркулирует термостойкий органических теплоноситель, тогда как в схеме на одном веществе (фторуглероде) температура перед турбиной может быть повышена минимум на С при переходе на одноконтурную схему. Кроме того, высокая термическая стойкость октафторциклобутана RC318 и декафторбутана R31-10 позволяет достичь высоких КПД цикла при повышении температуры газа перед турбиной до 550 оС. Результаты такого численного анализа показывают прирост КПД от 30% до 43% при повышении температуры перед турбиной.
Сформулированы неоспоримые преимущества фторуглеродов в качестве РВ утилизационных установок перед существующими аналогами: пожаро- и взрывобезопасность, возможность перевода теплосилового контура на одно РВ с исключением промежуточного теплообменника, высокая термодинамическая эффективность при повышении температурного уровня утилизации (до 580 С), а также преимущество в энергоэффективности над водопаровыми установками в диапазоне мощностей до 15 МВт.
Разработана двухстадийная схема утилизации генераторных газов на октафторциклобутане (C4F8).
Вся энергоустановка по утилизации генераторного газа рассчитана на 1 МВт, выделяемый при его сжигании в камере сгорания ГТУ. На первой стадии (рис.8) электрическая энергия вырабатывается в трехступенчатой турбоустановке мощностью 86 кВт с теплообменником (котлом-утилизатором) мощностью 2кВт, использующем высокий температурный потенциал генераторного газа (500С). На второй стадии (рис. 6) электрическая энергия вырабатывается в бинарной ПГУ с регенеративной ГТУ и утилизацией остаточной теплоты рекуператора в ПТУ на октафторциклобутане при непосредственном сжигании генераторного газа в камере сгорания с выделением тепловой мощности в 1 МВт.
Действительный КПД такого двухблочного комплекса с учетом затрат на привод насосов и компрессоров был оценен в 57%.
Рис. 8 - Принципиальная тепловая схема утилизации теплоты генераторных газов в трехступенчатой турбоустановке на октафторциклобутане на I стадии Общие преимущества теплосиловых циклов на фторуглеродах:
1. При одинаковых температурах перед турбиной и конденсации средняя интегральная температура подвода теплоты и, следовательно, действительный КПД цикла выше, чем КПД цикла на водяном либо газовом рабочем теле.
2. Такое преимущество в термодинамической эффективности наиболее заметно при низких и средних давлениях.
3. Фторуглеродная турбина (как и турбина ГТУ) работает в области перегретого пара, регенерация осуществляется не с помощью отборов из турбины, а передачей теплоты в рекуператоре.
4. Предварительные расчеты конструкции турбины показывают существенное снижение числа ступеней турбины и диаметра рабочего колеса, высоты лопаток и, в целом, габаритов фторуглеродной турбины по сравнению с аналогами водопаровых турбин подобной мощности.
5. Конденсация протекает при избыточном давлении, вследствие чего исключаются подсосы воздуха, упрощается конструкция конденсатора и уменьшается его вес. В зимние периоды возможно использование конденсатора с воздушным обдувом и проведение конденсации при отрицательных температурах (например, и при - 20 С), что дает более высокий рост КПД, чем повышение на такую же величину температуры перед турбиной.
6. Так как расширение происходит в газовой фазе, то отпадают проблемы, связанные с влажностью: пониженный КПД турбин, необходимость сепарации влаги, эрозийный износ и др.
7. Фторуглероды химически инертны, тяжелее воздуха, обладают свойствами пропеллентов и могут обеспечить пожаротушение в экстренных случаях.
8. Возможность применения более дешевых низколегированных сталей для основных аппаратов установки, что, особенно важно, скажется на упрощении конструкции и удешевлении турбоагрегата.
Основные недостатки применения фторуглеродов в качестве рабочих веществ теплосиловых циклов:
1. Повышенные требования к герметичности всей установки.
2. Ограничения по термостойкости (работа возможна до температуры 580С).
3. Поскольку в регенеративном теплообменнике передаются большие тепловые мощности, то предъявляются жесткие требования к эффективности регенеративных аппаратов.
4. Более низкие значения теплопроводности по сравнению с водой.
7 ТЕПЛОТЕХНИЧЕСКИЕ ИСПЫТАНИЯ ТЕПЛОСИЛОВОГО СТЕНДА НА ОКТАФТОРПРОПАНЕ Циркуляционный стенд на октафторпропане (C3F8) тепловой мощностью котельного агрегата 80 кВт был построен в качестве маломасштабного прототипа реальной теплосиловой установки на фторуглеродном рабочем веществе.
Принципиальная гидравлическая схема стенда представлена на рис. 9.
Рис. 9 - Принципиальная гидравлическая схема стенда и схема измерений рабочих параметров циркуляционного стенда Также на ней указаны точки измерений рабочих параметров стенда.
Основные аппараты стенда: котел каталитического горения на природном газе;
дроссельное устройство; регенеративный теплообменник; конденсатор; насосдозатор (НД); бак слива конденсата. Стенд оснащен системой измерения температур, давлений и расхода воды через конденсатор. Стенд оснащен также специальными вспомогательными системами: системой вакуумирования; узлом загрузки образцов и системой отбора рабочего вещества для проведения спектрометрического анализа; системой заправки стенда; системой взвешивания баллонов с отобранными пробами на основе прецизионных весов ВЛА-200;
системой термостатиров топливного бака с целью стабилизации давления тирования ю ста сжиженного газа и его расхода в котле каталитического горения; системой и его ского принудительной вентиля нтиляции отходящих продуктов горения.
В результате ресурсных и теплотехнических испытаний стенда их и установлены: параметры рабочих режимов (давление нагнета до 6,033 МПа, метры агнетания температура перед дросселем в диапазоне до 523,9 С); расход октафторпропана дросс расх (до 0,083 кг/с) и массовы опроводах; интегральные ссовые скорости РВ в аппаратах и трубопров коэффициенты теплоотд (например, в регенераторе сред лоотдачи реднелогарифмический температурный напор tср.лог.= 29,3 C, среднеинтегр апор нтегральный линейный коэффициент теплоперед и гидр передачи klср = 12,62 Вт/(мК)); потери гидравлического напора во всех основных аппаратах циркуляционного стенда; потери теплоты в х ап окружающую среду. Результаты спектрометрического анализа, произведенные Ре анал через 100, 300 и 600 часов работы, показали полную идентичность состава 600 ю ид отобранных образцов на стенде МЭИ и образцов ЭНИН, подвергшихся нагреву в ов на, подв специальных капсулах. После проведения полных ресурсных испытаний в течение лах. П ных и 600 часов состав проб не изменился и был эквивалентен исходному образцу проб нтен рабочего вещества. Подт терми одтверждена полная химическая и термическая стабильность октафторпропана в качестве рабочего тела теплосиловой установки в заданном каче ой уст интервале теплотехниче измерений: изменение температу перед турбиной хнических пературы до 525 C, давления - до до 6 МПа.
Создана математ терм атематическая модель процессов термогидродинамики в регенеративных теплооб в плообменниках и котлах-утилизаторах в среде MATHCAD.
Результаты сопоставле экспериментальных и расчетн характеристик тавления асчетных регенератора приведены в табл. 8. Расхождения составили не бо дены более 3%.
Табл.8 - Сопоставле экспериментальных и расчетных значений ставление рас температурных напоров генераторе поров и гидравлических потерь в регенер Температурный Экспери- Расчет Гидравличес- Экспер Эксп кспери- Расчет ме мент мент кий напор, Па мент напор, C 81 81,6 81,6 1,96105 1,881 tвх2- tвыхP 6,3 7,0 0,7,752105 0,7331tвых2- tвх1 P Анализ результат расчета (рис. 10) показывает, что интенсивность льтатов вает, теплоотдачи в регенерат меняется в диапазоне от 2050 до 2800 Вт/(м2.К) для нераторе 50 до потока нагреваемой сред и от 660 до 1170 Вт/(м2.К) среды сверхкритических параметров и от - для греющей среды.
Разработанная на основе проведенного эксперимен вычислительная я на римента программа расчета термогидродинамики теплообменных аппаратов позволяет а тер ых а рассчитывать локальны тепловые и гидродинамические характеристики для льные разных типов поверхностей ипов теплообменников и различного в направления потоко токов теплоносителей.
Рис. 10 - Расч асчетные значения коэффициентов теплоотдачи вдоль в те регенератора: 1 - для сверхкритического флюида C3Fкого внутри трубы (нагреваемая среда); 2 - (нагр для перегретого пара C3F8 в кольцевом о пар зазоре (греющая сре я среда) В результате анализа произведенных расчетов сделан вывод о том, что для повышения эффективности регенератора данной конструкции необходима реализация дополнительных мероприятий по интенсификации теплоотдачи со стороны греющей среды.
8 ТЕПЛОТЕХНИЧЕСКИЕ ХАРАКТЕРИСТИКИ ТЕПЛОНАСОСНЫХ УСТАНОВОК И ИХ ЭЛЕМЕНТОВ Испытания проведены на следующих установках: ТНУ мощностью 300 Вт на хладагенте R134a; ТНУ мощностью 20 кВт на диоксиде углерода; ТНУ мощностью 8 кВт на хладагенте R22.
Принципиальная схема стенда ТНУ на R22 приведена на рис.11.
На стендах установлены системы контроля параметров установки в процессе проведения экспериментов и автоматического снятия показаний с датчиков, размещённых на стенде. Температура и давление контролировались во всех характерных точках фреонового цикла. Кроме того, непосредственно измерялась ваттметром мощность трехфазного электропривода компрессора. В водяном контуре - источнике низкопотенциальной теплоты - измерялись расход, температуры и давления на входе и выходе из испарителя. В водяном контуре - аналоге системы теплоснабжения - измерялись также расход, температуры и давления на входе и выходе из конденсатора.
Рис. 11 - Принципиальная схема стенда ТНУ на RКонструкции конденсатора и испарителя была разработана автором, и базировались на применении пучков труб, ошипованных по технологии деформирующего резания.
Были определены интегральные и локальные коэффициенты теплоотдачи аппаратов ТНУ, сконструированных на основе интенсифицированных элементов.
В качестве примера приведены результаты обработки опытных данных испытаний ТНУ на R22 (табл. 9).
Проведенное сопоставление энергетической эффективности фреоновых ТНУ с ТНУ на диоксиде углерода при стандартных температурах прямой и обратной воды в системе теплоснабжения показало, что при отсутствии низкотемпературного потребителя, термодинамическая эффективность ТНУ на диоксиде углерода существенно ниже (более чем в два раза). Был сделан вывод о том, для повышения эффективности и внедрения ТНУ на диоксиде углерода необходима разработка специальных схем с включением низкотемпературных потребителей.
Методика сопоставления термодинамических циклов была построена на условиях фиксации температурных уровней внешних тепловых источников и итеративном подборе давлений в испарителе и конденсаторе с целью обеспечения одинаковых температурных напоров в каждом аппарате для сравниваемых циклов.
Проведенный расчетно-теоретический анализ на основе опытных данных позволил сделать вывод о том, что получение значительно большего, чем в базовом цикле на R22, значения коэффициента преобразования теплоты в таких циклах возможно при применении неазеотропных смесевых композиций на основе фторуглеродов и добавок гексафторида серы.
Табл. 9 - Теплотехнические характеристики теплообменных аппаратов ТНУ на R№ реж. Тепл. Привед. Средн. Расход Коэфф. Скор. Коэфф. Коэфф.
испыт. мощн. поверхн. темп. воды, теплопер. воды в теплоотд. теплоотд.
Q, кВт теплообм. напор G, кг/с пучках/ от воды в при конд.
Q k = пучка и центр. пучках/ и кипен.
t,оС, ft центр. тр. трубе, центр. фреона, Вт/м2К Fп+Fц, тр, м2 Wв, трубе, в ф м/с Вт/(м2К) Вт/(м2К) Конденсатор, приведенная рабочая поверхность Fк=0,21031 м1 13,67 0,16640 + 15,3 0,333 4250 1,33 /1,19 7970/6280 960,0432 14,14 19,9 0,250 3480 1,00 /0,89 6530/4820 79Испаритель, приведенная рабочая поверхность Fи=0,36005 м1 8,80 0,28488+ 4,9 0,600 4990 2,36 /2,12 7170/6810 1680,0752 9,22 6,5 0,500 3940 1,97 /1,69 7715/5880 85Результаты теплотехнических испытаний ТНУ на R22, R134а, диоксиде углерода, приведенные в данном подразделе, подготовили базу для постановки экспериментальных работ по внедрению смесевых фторуглеродных хладагентов и развитию направления по повышению энергоэффективности ТНУ на последующих этапах исследований.
9 ИССЛЕДОВАНИЕ ПРОЦЕССА ТЕПЛООБМЕНА НА МИКРОСТРУКТУРИРОВАННЫХ ПОВЕРХНОСТЯХ Целью работы было экспериментальное и теоретическое исследование теплообмена на новой микроструктурированной поверхности, предназначенной для интенсификации процессов кипения в испарителях ТНУ. В качестве объектов исследования были выбраны широко распространенный R134a и недавно синтезированный тепл сост теплоноситель фторуглеродного состава с высокой смачиваемостью FC-328 Для проведения исследования была использована 3284. ния экспериментальная установка Института Тепловой и Энергетической Техники уста Энерг Университета Падеборна (Германия), разработанная для исследования процессов еборна иссл кипения на одиночных трубах с микроструктурированны поверхностями.
чных ванными Установка была модифиц кривых кипения на трубах дифицирована с целью исследования кривы с такой поверхностью, расположенных одна над другой в верт, р в вертикальной плоскости (тандеме). Принципиаль схема установки представлена на рис. 12. Установка пиальная на на для исследования теплообмена при кипении в большом объеме на пучке труб тепло м объ состоит из замкнутого контура с циркуляцией исследуемого рабочего вещества. В ого ко состав контура входит рабочая камера А, в которой располагаются трубы со одит распо специально подготовле раб отовленными поверхностями. Под рабочей камерой А располагается дополнит ый о олнительный нагреватель D, который обеспечивает нагрев жидкости, охлажденной в конденсаторе В, до состояния насыщ нной асыщения.
Рис. 12 - Принц льной установки ринципиальная схема экспериментальной Проведено три типа экспериментов с измерением кривых кипения на три т ием к верхней трубе: 1) на обеих трубках устанавливались одинаковые значения на сь од теплового потока; 2) на н лась п на нижней трубке А425 устанавливалась постоянная тепловая нагрузка в 8 кВт/м2, 50 кВт/м2, или 125 кВт/м2, при этом верхней трубке (либо 50 ве А428, либо А437) устанавливались все тепловые нагрузки в ди устана и в диапазоне от 2 кВт/мдо 125 кВт/м2; 3) третий тип экспериментов был сходен второму, но постоянная ретий втор тепловая нагрузка устанавливалась на верхней трубке, при одновременном уста ке, п измении тепловой нагруз нагрузки на нижней трубке.
Геометрические характеристики поверхности микрострукту харак руктурированных трубок, исследованных в данной работе, приведены в таблице 10.
нной Результаты исследован кипен едования представлены в виде кривых кипения, пример которых приведен на рис. 13. В результате исследования установл что на таких ановлено, поверхностях необходи лишь незначительный перегре теплоотдающей бходим ерегрев поверхности для начала устойчивого кипения жидкос а тандемное нача идкости, расположение трубок даёт незначительное улучшение процесса кипения на бок д е пр верхней трубке по сравне сравнению с одиночной тестовой трубой.
Табл.10 - Геометр икроструктурированной ометрические характеристики микрос поверхности Диаметр Объем бъемная Угол наклона Высота Плотность Структура ористость штырьков, штырьков, m штырьков, см-2 штырьков, пористо m m структур уктуры, % градусы A425 51 7,6105 7,63 34 A428 55 5,0106 3,5 46 A437 58 3,3106 3,7 37 +40 / - - верхн верхняя;
- боков боковая;
- нижн нижняя термопары.
Рис. 13 - Кр нии 5 бар на Кривые кипения R134a при давлении мик А4микроструктурированной трубке А4Особое прикладн значение имеют результаты по устойчивости и кладное ты п эффективности процесса цесса кипения при низких давлениях.
ЗАКЛЮЧЕНИЕ Е В рамках провед на ко роведенного исследования разработана концепция внедрения фторорганических рабочих веществ в качестве рабочих тел теплосиловых рабо очих установок в системы малой распределенной энергетики, атомную отрасль и мы м ки, а теплонасосную технику.
нику.
В результате выполнения комплекса экспериментал е вы ентальных и расчетнотеоретических исследова олучены новые данные о едований в течение 13-него периода получе термодинамических свойствах и теплотехнических характе х сво рактеристиках ряда как индивидуальных, так и с ганического состава.
ак и смесевых рабочих веществ фторорганич Результаты измере РВ змерения плотности (удельного объема) РВ были получены на двух экспериментальны установках методически независ альных езависимыми способами, имели высоких метроло но р етрологический уровень и существенно расширили диапазон известных данных. Разработанные экспериментально-обосн Раз обоснованные уравнения состояния на их основе были использованы для расчета теплосиловых циклов в нове тепл диапазонах рабочих пара устан х параметров реальных энергетических установок.
Были разработан принципиальные тепловые схемы второго контура аботаны схем реакторной установки на быстрых нейтронах, а также схемы утилизации теплоты вки на хемы выхлопов ГТУ и теплот генераторного газа. Термодинамич еплоты намический анализ таких схем показал их высокую гическую эффективность, сокую термодинамическую и технологичес а также их преимущества в данной области применения в сравнении с аналогами:
водопаровыми установками и установками на других РВ органического состава.
Высказано также предположение, что внедрение фторуглеродов в качестве рабочих тел теплосилового контура АЭС позволит повысить не только термодинамическую эффективность, но и экологическую и технологическую безопасность энергоблока АЭС.
Результаты ресурсных и теплотехнических испытаний аппаратов циркуляционного стенда - маломасштабного прототипа теплосиловой установки - на октафторпропане подтвердили высокую эффективность, термостабильность, безопасность фторуглеродов в качестве РВ. Опытные теплотехнические характеристики были проанализированы с помощью разработанной математической модели процессов термогидродинамики для аппаратов стендовой установки.
Были разработаны оригинальные схемы и теплообменные аппараты ТНУ, на основные технические решения которых были получены авторские свидетельства.
В результате проведения теплотехнических испытаний ТНУ были получены экспериментально обоснованные соотношения для расчета теплоотдачи при кипении и конденсации исследованных хладагентов на рельефных поверхностях теплообмена в элементах установок. Экспериментально подтверждена возможность повышения коэффициента теплопередачи при использовании таких поверхностей в 3-5 раз. Полученные результаты могут быть использованы при проектировании нового поколения теплонасосных установок. Проведенный расчетно-теоретический анализ показал перспективу внедрения в качестве рабочего вещества ТНУ неазеотропных смесевых композиций на основе фторуглеродов и гексафторида серы и дальнейшего развития направления по повышению энергоэффективности ТНУ.
На кафедре Теоретических основ теплотехники МЭИ за последние 10 лет выполнен ряд хоздоговорных и госбюджетных работ по теме исследования процессов теплообмена для фреоновых рабочих веществ на микроструктурированных и оребренных поверхностях. Результаты исследования теплообмена на новых микроструктурированных поверхностях существенно пополнили имеющуюся базу экспериментальных и расчетных данных и могут быть использованы при разработке интенсифицированных элементов аппаратов ТНУ широкого диапазона мощностей.
ОСНОВНЫЕ РЕЗУЛЬТАТЫ ДИССЕРТАЦИОННОЙ РАБОТЫ ИЗЛОЖЕНЫ В СЛЕДУЮЩИХ ПУБЛИКАЦИЯХ Тезисы докладов и статьи, опубликованные в материалах конференций 1. Алтунин В.В., Сухих А.А. Исследование термических свойств в газовой фазе и на кривых конденсации смесей бесхлорных фторсоединений // Тез. докл.
межд. теплофиз. школы. - Тамбов: ТГТУ. 1995. с. 77.
2. Соколов Е. Ю., Сиденков Д.В., Сухих А.А. Лабораторный стенд и комплекс программ для исследования энергетических характеристик теплонасосной установки // Тез. докл. 8 межд. научно-техн. конференции студентов и аспирантов 28.02-1.03. 2002г., М.: МЭИ. Т. 3, с.130.
3. Закопырин М.А., Сухих А.А. Экспериментальное исследование теплофизических свойств бинарной смеси фреона HFC23 и фторэфира HFE347mcc в газовой фазе и на кривой конденсации: Тез. докл. 8 межд. научно-техн. конфер.
студентов и аспирантов. - М.: МЭИ. 2002. Т.3. с. 119-120.
4. Закопырин М.А., Сухих А.А., Алтунин В.В. Экспериментальное исследование термодинамических свойств в перегретых и насыщенных парах смесей перфторпропана с многоатомными фторсоединениями: Тез. докл. 8 межд.
научно-техн. конфер. студентов и аспирантов. - М.: МЭИ. 2002. Т.3. с. 120-121.
5. Сухих А.А., Закопырин М.А., Алтунин В.В. Экспериментальное исследование теплофизических свойств гептафторбутанолового эфира и его смесей с перфторпропаном // Тез. докл. Х Российской конфер. по теплофизическим свойствам веществ. - Казань: КГТУ. 2002. с. 46 - 47.
6. Сухих А.А., Сиденков Д.В. Исследование энергетических характеристик цикла теплонасосной установки. Натурный и вычислительный эксперимент.
Программа HOLCON-2. // Тез. докл. Всероссийской научно-практич. конфер.
Человеческое измерение в информационном обществе. 2003г., стр.51-52, М.:
ВВЦ, павильон №57. Минобразования РФ.
7. Алтунин В.В., Сухих А.А. Термодинамические исследования фторэфира HFE347mcc и его бинарных смесей с озонобезопасными фторсоединениями как альтернативных хладагентов. Тез. докл. 2 межд. научно-практич. конфер.:
Природные холодильные агенты - альтернатива глобальному потеплению С.Петербург: ГУНТиПТ. 2003. с.123-127.
8. Сухих А.А., Алтунин В.В., Петров М.В. Термодинамические исследования новых рабочих веществ парокомпрессионных тепловых насосов. Тез. докл. 9 межд.
научно-техн. конфер. студентов и аспирантов. 2003г., М.: МЭИ. Т. 3, стр.18-19.
9. Закопырин М.А., Сухих А.А., Алтунин В.В. Давление насыщенного пара и объемные соотношения фторэфира HFE347mcc и его смесей с многоатомными фторсоединениями. Труды 2 межд. научно-практич. конфер.:
Низкотемпературные и пищевые технологии в XXI веке. С.-Петербург:
ГУНТиПТ. 2003. с.235-238.
10. Сухих А.А., Сиденков Д.В. Расчетно-экспериментальное исследование теплообмена в испарителе теплонасосной установки. Сборник научных статей Республиканского научно-технического семинара Тенденции развития теоретической теплотехники. Ташкент: ТГТУ. 2004. с. 11-14.
11. Сухих А.А., Сиденков Д.В., Закопырин М.А. Комплексное исследование энергетической эффективности теплонасосных установок // Тез. докл. электронной конфер. по подпрограмме Топливо и энергетика научно-технической программы Научные исследования высшей школы по приоритетным направлениям науки и техники. М.: МЭИ, ноябрь 2004, с.26.
12. Сухих А.А., Утенков В.Ф., Закопырин М.А. Экспериментально обоснованное уравнение состояния HFE347mcc для расчета циклов тепловых насосов // Тезисы докладов 11 Российской конфер. по теплофизическим свойствам веществ. С.-Петербург, 2005, с.42-43.
13. Сухих А.А., Закопырин М.А. Экспериментальное исследование рvTx - поверхности смесевого хладагента SF6-C3F8. Тез. докл. 11 Российской конфер. по теплофиз. свойствам веществ. С.-Петербург, 2005, с.163-164.
14. Сухих А.А., Закопырин М.А., Утенков В.Ф. Уравнение состояния Карнахана - Старлинга - де Сантиса эфира HFE-347 mcc // Материалы докладов национальной конфер. по теплоэнергетике НКТЭ - 2006. Казань, 2006. С. 69-72.
15. Кузнецов К.И., Скородумов С.В., Сухих А.А., Утенков В.Ф.
Экспериментальные данные о плотности октафторпропана при повышенных температурах // Труды XIII Российской конфер. по теплофизич. свойствам веществ.
Новосибирск, 2011. - 1 CD-ROM/ -ISBN 978-5-89017-030-9.
16. Кузнецов К.И., Скородумов С.В., Сухих А.А., Утенков В.Ф. Разработка уравнения состояния вириального типа для расчета термодинамических свойств октафторпропана в газовой области // Труды XIII Российской конфер. по теплофизич. свойствам веществ. Новосибирск, 2011. - 1 CD-ROM/ -ISBN 978-589017-030-9.
17. Ежов Е.В., Милютин В.А., Сухих А.А. Теплотехнические испытания теплосилового стенда на октафторпропане// Тез. докл. национальной конфер.:
Повышение эффективности, надежности и безопасности работы энергетического оборудования ТЭС и АЭС. - Москва: ФБГОУ ВПО НИУ МЭИ. 2012. с. 109.
18. Антаненкова И.С., Сухих А.А. Разработка методики сравнения термодинамической эффективности холодильных и теплонасосных установок // Тез. докл. национальной конфер.: Повышение эффективности, надежности и безопасности работы энергетического оборудования ТЭС и АЭС. - Москва:
ФБГОУ ВПО НИУ МЭИ. 2012. с. 99.
19. Антаненкова И.С., Сухих А.А. Повышение энергоэффективности теплонасосных установок на диоксиде углерода // Тез. докл. национальной конфер.: Повышение эффективности, надежности и безопасности работы энергетического оборудования ТЭС и АЭС. - Москва: ФБГОУ ВПО НИУ МЭИ. 2012. с. 126.
20. Устинов Д.А., Устинов В.А., Сухих А.А., Устинов А.К. Интенсификация теплообмена на микроструктурированных поверхностях // Тез. докл. национальной конфер.: Повышение эффективности, надежности и безопасности работы энергетического оборудования ТЭС и АЭС. - Москва: ФБГОУ ВПО НИУ МЭИ. 2012. с. 128.
Работы, опубликованные в научных журналах, включенных в перечень ВАК Министерства образования и науки Российской Федерации 1. Сухих А.А., Генералов К.С., Акимов И.А. Испытания теплового насоса для теплоснабжения индивидуального дома // Труды МГУИЭ: Техника низких температур на службе экологии. М: МГУИЭ. 2000г. с.49 -53.
2. Сухих А.А., Алтунин В.В., Закопырин М.А. Экспериментальное исследование термических свойств бинарной системы гексафторид серыперфторпропан // Вестник МЭИ. 2002. №2. с. 86-903.
3. Сухих А.А., Закопырин М.А., Утенков В.Ф. Уравнение состояния вириального типа и таблицы термодинамических свойств альтернативного хладагента HFE-347mcc // Холодильная техника. 2007. №5.
4. Сухих А.А., Закопырин М.А., Утенков В.Ф. Экспериментальное определение плотности и построение локального уравнения состояния бинарной системы хладагентов R14-HFE347mcc в газовой фазе // Холодильная техника. №6.
2008. с.43-45.
5. Сухих А.А., Закопырин М.А., Утенков В.Ф. Экспериментальное исследование термических свойств бинарных смесей фторэфира HFE347mcc с трифторметаном HFC23 и разработка многоконстантного уравнения состояния вириального типа // Вестник МАХ. 2009. №2. с. 16-21.
6. Кузнецов К.И., Скородумов А.В. Сухих А.А. Экспериментальное исследование рvT-поверхности рабочих тел фторуглеродного состава // Энергосбережение и водоподготовка. №2. 2009. с.28-32.
7. Сухих А.А. Исследование плотности декафторбутана в широком диапазоне температур и давлений / Сухих А.А., Кузнецов К.И., Закопырин М.А., Утенков В.Ф., Скородумов С.В // Известия ВУЗов. 2009. №7,8 стр. 27-36.
8. Сухих А.А., Кузнецов К.И., Утенков В.Ф. Построение уравнения состояния декафторбутана в газовой фазе // Проблемы энергетики. 2010. № 11-12.
с.46-51.
9. Устинов В.А., Сухих А.А. Исследование процессов теплообмена на микроструктурированных поверхностях в испарителе теплонасосной установки // Энергосбережение и водоподготовка. 2010. №2 (64). С.43-46.
10. Сухих А.А., Закопырин М.А, Джураева Е. В. Экспериментальное исследование плотности бинарных смесей фторэфира HFE347mcc с хладоном R218 и разработка многоконстантного уравнения состояния вириального типа // Вестник МАХ. 2010. №1. с. 23-29.
11. Сухих А.А., Милютин В.А., Антаненкова И.С. Термодинамическая эффективность фторуглеродов как рабочих тел в теплосиловых циклах АЭС / Электрические станции. 2010. №10. с. 2-8.
12. Сухих А.А. Теплотехнические испытания теплонасосной установки на диоксиде углерода / Сухих А.А., Антаненкова И.С., Кузнецов В.Н., Сотсков С.А. // Вестник МЭИ. 2011. №3. С. 10 - 16.
13. Сухих А. А. Термодинамический анализ схемы замещения пароводяного контура на фторуглеродный в парогазовых установках / А. А. Сухих, К. И.
Кузнецов, В. А. Милютин // Вестник МЭИ. 2011. № 3. с. 17-22.
Изобретения по теме диссертации 1. Сухих А.А., Антаненкова И.С. Патент РФ на полезную модель №782Теплообменный аппарат. Приоритет полезной модели 26.06.2008.
2. Сухих А.А., Антаненкова И.С. Патент РФ на полезную модель №75879 Теплонасосная установка. Приоритет полезной модели 16.04.2008.