Авторефераты по всем темам >>
Авторефераты по физике
На правах рукописи
Иванов Радион Игоревич
ПОВЫШЕНИЕ ЭФФЕКТИВНОСТИ ПРОЦЕССА СМЕСЕОБРАЗОВАНИЯ В ГОРЕЛОЧНЫХ УСТРОЙСТВАХ С ИСПОЛЬЗОВАНИЕМ ОСОБЕННОСТЕЙ ТЕЧЕНИЯ В ВИХРЕВОМ ПРЯМОТОЧНОМ ЭЖЕКТОРЕ
Специальность 01.04.14 - Теплофизика и теоретическая теплотехника
АВТОРЕФЕРАТ
диссертации на соискание ученой степени кандидата технических наук
Рыбинск - 2012
Работа выполнена в федеральном государственном бюджетном образовательном учреждении высшего профессионального образования Рыбинский государственный авиационный технический университет имени П. А. Соловьева
Научный консультант: доктор технических наук, профессор Пиралишвили Шота Александрович
Официальные оппоненты: доктор технических наук, профессор кафедры теплотехники и тепловых двигателей Самарского государственного аэрокосмического университета имени академика С.П. Королева Бирюк Владимир Васильевич;
кандидат технических наук, профессор, заведующий кафедрой промышленной теплоэнергетики Липецкого государственного технического университета Губарев Василий Яковлевич.
Ведущее предприятие: открытое акционерное общество Научнопроизводственное объединение Сатурн, г. Рыбинск
Защита диссертации состоится 23 мая 2012 г. в 1200 часов на заседании диссертационного совета в федеральном государственном бюджетном образовательном учреждении высшего профессионального образования Рыбинский государственный авиационный технический университет имени П.А. Соловьева по адресу: 152934, г. Рыбинск, Ярославская область, ул. Пушкина, 53, ауд. Г-237.
С диссертацией можно ознакомиться в библиотеке федерального государственного бюджетного образовательного учреждения высшего профессионального образования Рыбинский государственный авиационный технический университет имени П.А. Соловьева
Автореферат разослан 23 апреля 2012 г.
Ученый секретарь диссертационного совета Каляева Надежда Анатольевна
ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ
Актуальность работы. Проблема повышения эффективности процессов преобразования энергии требует снижения затрат на их реализацию при соблюдении достигнутых условий по экономичности, компактности и надежности, что обуславливает необходимость поиска оптимальных условий эффективного и низкоэмиссионного сжигания углеводородного топлива.
Характерные особенности ограниченных закрученных течений могут быть полезны при разработке такого рода устройств. Повышение качества рабочего процесса возможно за счет высокой интенсивности турбулентности, наличия рециркуляционных зон, радиального перераспределения полной энтальпии и других свойств закрученного потока.
Вихревой эжектор при малых габаритах и простоте конструкции получил распространение в энергетике, машиностроении и других отраслях промышленности. Несмотря на накопленные экспериментальные данные, остаются вопросы, требующие дополнительного исследования: недостаточно изучено распределение давления вдоль приосевой области эжектора; отсутствует адекватная расчетная модель таких устройств, построенная на критериальных зависимостях; не исследованы акустические характеристики вихревых эжекционных устройств; не изучены процессы смесеобразования в проточной части вихревых эжекторов.
Решение этих задач позволит получить дополнительную информацию о закрученных течениях в ограниченном объеме и расширить область применения вихревых эжекторов.
Цель диссертационной работы состоит в исследовании особенностей течения закрученного потока газов в вихревом прямоточном эжекторе для организации процесса смешения топлива с окислителем и последующего горения с целью создания эффективных горелочных устройств бытового и технологического назначения.
Для достижения поставленной цели необходимо решить ряд задач:
Ц численно и экспериментально исследовать рабочий процесс прямоточного вихревого эжектора с оценкой его функциональной пригодности для указанных технических направлений;
Ц провести анализ распределения концентрации компонентов топлива в выходном сечении горелочного устройства и измерить эмиссионные характеристики;
Ц разработать критериальные уравнения расчета проектируемых конструкций и методику математической обработки результатов эксперимента с использованием теории подобия и метода анализа размерностей;
Ц с использованием полученных данных и разработанной методики расчета спроектировать горелочное устройство и инфракрасную горелку, работающие по механизму обедненной предварительно перемешанной топливовоздушной смеси.
Методы исследований. Для решения поставленных задач использованы: аналитические методы на базе основополагающих закономерностей газовой динамики, термодинамики, теории подобия и методов анализа размерностей, методы постановки и обработки теплофизического эксперимента, численное моделирование турбулентных течений.
Научная новизна работы Выявлены режимы работы вихревого прямоточного эжекторасмесителя, позволяющие организовать подачу сжигаемого газа при пониженных давлениях на входе; повысить качество процесса смесеобразования;
обеспечить приемлемую равномерность поля концентрации топлива на выходе из смесителя ( 0,24) и выбросы оксидов азота NOx не более 10 ppm.
Определены акустические характеристики вихревого эжектора, которые позволяют осуществить необходимые защитные мероприятия при его включении в технологические цепочки различного назначения.
Достоверность и обоснованность научных положений обеспечивается применением основополагающих уравнений термогазодинамики, положений теории подобия и метода анализа размерностей, обработкой опытных данных с использованием методов математической статистики. Экспериментальные результаты получены на сертифицированном метрологическом оборудовании, прошедшим необходимую проверку. Достоверность подтверждается соответствием расчетных и опытных данных, а также совпадением с результатами исследований других авторов.
Практическая значимость Разработанная методика позволяет по данным технического задания спроектировать модули смесеподготовки для различных топливосжигающих устройств бытового и технологического назначения с низким значением эмиссионных выбросов.
Экспериментальные данные по исследованию акустической неустойчивости в эжекторе позволят осуществить необходимые защитные мероприятия при его включении в технологические цепочки различного назначения.
Разработана вихревая инфракрасная горелка, предназначенная для отопления бытовых и производственных помещений, термической обработки материалов в авиационной, машиностроительной, химической и других отраслях. Эмиссионные выбросы горелки по СО и NOx не превышают 10 ppm.
Новизна инфракрасной горелки защищена Патентом РФ.
Апробация работы Основные результаты выполненных исследований докладывались:
Ц I Международная научно-техническая конференция, посвященная летию основателя Рыбинской школы теплофизиков доктора технических наук, профессора Пиралишвили Ш.А.;
Ц Всероссийская молодежная научная конференция с международным участием л Х Королевские чтения 2010;
Ц Ежегодный конкурс научно- исследовательских работ студентов вузов Ярославской области 2010;
Ц Пятая Российская Национальная Конференция по Теплообмену Международная научная школа Проблемы газодинамики и теплообмена в энергетических технологиях 2010;
Ц XVI и XVII Школа-семинар молодых ученых и специалистов под руководством академика РАН А.И. Леонтьева л Проблемы газодинамики и тепломассообмена в аэрокосмических технологиях 2009, 2011;
Ц Международная научная школа Проблемы газодинамики и тепломассообмена в энергетических технологиях, Москва, 2011;
Публикации: ключевые положения опубликованы в 7 научных работах, 3 из них в рецензируемых изданиях по списку ВАК.
Структура и объём диссертации. Диссертация содержит 126 стр. машинописного текста, 123 рис. и состоит из введения, 5-ти глав, заключения, списка литературы из 116 наименований.
ОСНОВНОЕ СОДЕРЖАНИЕ РАБОТЫ
Во введении обоснована актуальность темы диссертации, выполнена краткая характеристика работы, сформулированы цели и задачи исследований, отмечена научная новизна, практическая значимость полученных результатов и положения, выносимые на защиту.
В первой главе диссертации представлена история развития и обзор работ по исследованию вихревых и струйных эжекторов.
Рассмотрены труды М.Г.Дубинского, А.П.Меркулова, В.Т. Волова, В.И. Метенина, В.И. Епифановой, А.Д. Суслова, В.В. Бирюка и др. - ведущих ученых в области изучения процессов вакуумирования и эжектирования, реализуемых за счет формирования приосевой области пониженного давления.
В результате анализа существующих экспериментальных исследований вихревых эжекторов различных конструкций определен диапазон эффективной работы по геометрическим и режимным параметрам.
Аналитический обзор позволил сформировать цель диссертационной работы и выполнить постановку решаемых задач.
Вторая глава работы посвящена численному моделированию газовой динамики воздушного потока проточной части вихревого эжектора.
Для описания процессов течения газа использована система дифференциальных уравнений Навье-Стокса, которая замыкалась моделью турбулентности и уравнением состояния. Согласно проведенным поверочным газодинамическим расчетам оптимальной для численного решения поставленной задачи является SST k - модель турбулентности. Расчетная сетка области течения (рис.1) выполнена гибридной во внутренней полости (тетрагональные элементы в основном потоке, гексагональные элементы в пристеночной области). На входе в сопловой ввод эжектора подавался возРис. 1. Расчетная сетка. 1 - сопловой ввод, 2 - вихревая камера, 3 - камера дух с температурой 293 К при давленисмешения, 4 - осевой диффузор ях 0,15 Е 0,4 МПа, в качестве эжектируемого потока выступал воздух при атмосферном давлении и температуре окружающей среды или метан.
При исследовании распределения статического давления (рис.2) по Рис. 2. Распределение статического давления по длине эжектора при значении коэффициента эжекции n=0: 1) = 2; 2) = 3; 3) = длине эжектора установлено, что зона пониженного давления сосредоточена в пределах небольшой области с относительным диаметром dо 0,3 и растягивается до окончания камеры смешения lо 1,38 (рис. 2). При изменении режимных параметров она не деформируется.
Распределение статического давления в сечении соплового ввода описывается квадратичной функцией, характер которой сохраняется для различных режимных и геометрических параметрах.
инии тока в проточной части эжектора (рис. 3Ц4) имеют достаточно сложный трехмерный характер, обусловленный интенсивными радиальными Рис. 3. Распределение линий тока и зона рециркуляции в Рис. 4. Распределение линий эжекторе тока пассивного потока пульсациями закрученного потока. На периферии диффузорной части эжектора формируется тороидальный вихрь, способствующий интенсификации массообменных процессов. В результате исследований обнаружено возвратное течение воздуха, граница которого располагается в начале камеры смешения и смещается вверх или вниз по потоку в зависимости от режимных параметров.
Установлено, что независимо от относительной длины и диаметра сопла пассивного потока максимальному значению коэффициента эжекции соответствует степень повышения давления = 2 (рис. 5Ц6).
Рис. 5. Зависимость коэффициента эжекции Рис. 6. Зависимость коэффициента эжекции от относительной длины сопла пассивного от относительного диаметра сопла пассивного 1) = 2; 2) = 3; 3) = 1,5; 4) = lс 0,9 :
потока dс 0,1: 1) = 2; 2) = 4; 3) = Определена относительная длина сопла пассивного потока lс 0,8 0,9, при которой достигается максимальное значение коэффициента эжекции (рис. 5). Причем длина сопла, соответствующая минимуму давления, не совпадает с длиной соответствующей максимальному значению коэффициента эжекции. Таким образом, пониженное давление играет не доминирующую роль в процессе вихревой эжекции. Повысить коэффициент эжекции позволяет периферийный вихрь, обладающий высокой окружной компонентой скорости.
При исследовании влияния диаметра сопла пассивного потока на значение коэффициента эжекции установлено, что рост коэффициента эжекции на любом режиме работы продолжается до определенного значения dс 0,(рис. 6), при котором выполняется условие равенства статического давления в вихре давлению среды, в которую происходит истечение. При дальнейшем увеличении диаметра dс 0,1 эжекция снижается.
Предварительные исследования рабочего процесса вихревого эжектора позволили определить диапазон работы по геометрическим и режимным параметрам, обеспечивающим организацию процесса смешения.
При подаче в качестве пассивного потока метана исследовались процессы смешения (рис. 7). В результате взаимодействия периферийного вихря, обладающего высокой окружной компонентой скорости, с приосевой струей метана на границе раздела потоков возникают касательные сдвиговые напряжения, которые придают эжектируемой струе окружной импульс. Вследствие этого активная струя увлекает за собой пассивный поток топлива, который по мере продвижения приобретает окружную компоненту скорости.
Таким образом, периферийный воздушный вихрь формирует в приосевой области сонаправленный ему метановый вихрь, который при удалении от сопла подвода пассивного потока перемешивается с периферийным воздушным потоком.
На выходе из смесителя концентрация метана несколько выше на периферии. Это связано тем, что массовый расход окислителя в выходном сечении распределен неравномерно: основная масса воздуха, как и метана, отбрасывается в рециркуляционную зону на периферию, в тоже время часть воздуха подсасывается в центральную область из атмосферы. За счет образования в дифРис. 7. Распределение массовой доли метана фузорном участке тороидальной зоны обратных токов и генерации интенсивной турбулентности в радиальном направлении происходит интенсификация процессов массообмена. Это обеспечивает смешение компонентов топлива и окислителя по объему отмеченной области (рис.8-9).
Рис. 8. Распределение коэффициен- Рис. 9. Распределение коэффициента избытта избытка воздуха в сечении 3Ц3 ка воздуха в сечении 3Ц3.
1) = 2,7; 2) = 3,7; 3) = 1,В главе 3 для реализации необходимых исследований разработан стенд (рис.10) и опытный образец (рис.11). Опытный образец прямоточного вихревого эжектора выполнен с возможностью изменения его геометрии: сопла пассивного потока в диапазоне dс 0,025...0,25 и lс 0...1,9; длины вихревой камеры LВК 11,5; площади проходного сечения соплового ввода в диапазоне Fc 0,03 0,1; конфигурации диффузора и геометрических размеров зазора щели . Все линейные размеры, как принято в теории вихревых труб, отнесены к минимальному диаметру вихревой камеры в сечении е примыкания к сопловому вводу.
Рис. 10. Схема экспериментального стенда: 1 - поршневой компрессор; 2 - ротаметр активного потока; 3 - манометр сжатого воздуха; 4 - термопара; 5 - ротаметр пассивного потока; 6 - газовый баллон; 7 - редуктор давления; 8 - низконапорный резервуар; 9 - опытный образец; 10 - измерительное оборудование При анализе погрешности измерений установлено, что применение измерительных датчиков в комплексе с современным метрологическим обеспечением, позволит провести экспериментальное исследование с допустимой степенью точности.
В четвертой главе приведены результаты опытного исследования вихревого эжектора, направленные на поиски рабочего режима, позволяющего организовать процесс смешеРис. 11. Опытный образец. 1 - сопло пассивного ния при пониженных давлениях потока, 2 - торцевая стенка, 3 - сопловой ввод, газообразного топлива.
4 - патрубок подвода сжатого воздуха, 5 - Распределение статиче вихревая камера, 6 - камера смешения, 7 - конический осевой диффузор, 8 - кольцевая ского давления вдоль оси эжекщель, 9 - стопорная гайка, 10 - регулируемый тора в зависимости от степени конус повышения давления и относительной площади сечения соплового ввода изображено на рис.12. Площадь сечения соплового ввода, в исследуемом диапазоне давлений на входе, слабо влияет на величину вакуумирования и давления на оси эжектора.
Минимальные значения давления зафиксированы в районе относительной длины lo 0...0,6 и относительного радиуса 0,2 r 0,2. С увеличением относительной длины эжектора происходит постепенный рост давления до атмосферного, при приближении к периферии 0,7 r 1 темп роста давления замедляется. Максимальное и устойчивое значение коэффициента эжекции при постоянном значении достигается при относительной длине сопла пассивного потока lс 0,5...0,75 (рис.
13). В случае уменьшения относительной длины сопла пассивного потока lc 0,5 возникают пульсации расхода эжектируемого потока, и эжектор работает в области неустойчивых режимов.
Увеличение относительной длины сопла пассивного потока lc 0,75 способствует равномерному падению коэффициента Рис.12. Распределение статического давления эжекции. В большинстве случаев вдоль оси эжектора в работе эжектора наблюдается эксперимент :1) 2; 3) 3; 5) 4;
вторая область устойчивой рабочисленный расчет : 2) 2; 4) 3; 6) ты, по размерам значительно уступающая первой, которая сосредоточена в диапазоне значений lс 0,95...1,05.
Рис. 13. Зависимость коэффициента эжекции Рис. 14. Зависимость коэффициента эжекции от относительной длины сопла пассивного от относительного диаметра сопла пассивнопотока. го потока эжектора.
LВК 1; dc 0,1; Fc 0,1 LВК 1; lc 0,65; Fc 0,числ. расчет :1) 2; 3) 3; 5) 4; числ. расчет :1) 2; 3) 3; 5) 4;
эксперимент : 2) 2; 4) 3; 6) 4 эксперимент : 2) 2; 4) 3; 6) При увеличении относительной длины вихревой камеры LВК расширяется диапазон неустойчивой работы 0 lc 0,75, область устойчивых режимов при этом смещается 0,9 lc 1,2.
Выполнено исследование зависимости коэффициента эжекции от диаметра сопла пассивного потока, которая имеет экстремальный характер (рис.
14). Максимум функции для сопла dc 0,2 соответствует радиусу, определяющему границу зоны пониженного давления, на которой выполняется условие равенства статического давления в вихре давлению среды, в которую происходит истечение.
При исследовании акустических характеристик вихревого эжектора обнаружено двукратное скачкообразное изменение частоты акустических колебаний при значении = 2,6, что равносильно переходу на вторую гармонику (рис.15). В диапазоне изменения степени расширения от 0 до 2,6 и после скачка частота акустических колебаний как функция перепада давления изменяется монотонно без разрывов.
По результатам экспериментальных исследований получены соответствующие зависимости:
безразмерного давления вдоль оси эжектора 0,0 0,9 е0,3l ;
коэффициента эжекции от длины сопла пассивного потока 0,n 0,18lC (sin( 0,6) log( 0,85));
коэффициента эжекции от площади сечения закручивающего устройства 0.n 1,39 FC0.69;
коэффициента эжекции от относительного диаметра эжектора n 0,83 0,89 (ln(dc ) 0,46) ;
частоты колебаний основного тона от степени повышения давления ( 1)0,.
Sh 0,Рис. 15. Зависимость частоты от степени повышеM ния давления при различных значениях Fc Экспериментальные исследования процесса смесеобразования в вихревом эжекторе подтвердили результаты численного расчета по характеру распределения концентрации метана на выходе (рис. 16Ц17).
Метан распределяется в воздухе достаточно равномерно: значения концентраций вдоль горизонтального и вертикального диаметров отличаются в предельном случае на 6 %; наибольшее расхождение в результатах измерения концентрации в центре и на периферии составило 25 %. При проведении экспериментального исследования установлено, что изменение площади сечения соплового ввода в диапазоне Fc 0,03...0,1 незначительно влияет на эффективность процесса перемешивания, таким образом, завершающая стадия смесеподготовки происходит в рециркуляционной зоне диффузорного участка. Оптимальной для подачи газообразного топлива является относительная длина сопла пассивного потока. Степень равномерности l 0,c распределения метана в смеси оценивалась с помощью параметра :
сmax сmin , сcр сmax, сmin, сср - максимальная, минимальная и средняя концентрации топгде ливного газа в смеси. Минимальное значение параметра , полученное в результате экспериментального исследования, составило 0,12 ;
максмальное - 0,43.
Для описания процесса смешения в выходном сечении эжектора получена соответствующая зависимость ССН 0,13 r 1,1.
Рис. 16. Зависимость концентрации топлива от Рис. 17. Зависимость относительной концентраотносительного радиуса диффузора ции топлива от относительного радиуса - вертикальные измерения;
1 эксперимент, dc 0,1; lc 0,55;
Ц горизонтальные измерения.
2 эксперимент, dc 0,1; lc 0,68;
dc 0,1; lc 0,55; Fc 0,3 эксперимент, данные В.Е.Костюка и др.
Выполненные исследования вихревого эжектора в качестве горелочного устройства показали, что на выходе из диффузора формируется развитая рециркуляционная зона с нарушением симметрии в верхней части пламени, что объясняется наличием естественной конвекции (рис. 18).
При изменении геометрических или режимных параметров форма факела не менялась, происходило лишь его небольшое смещение вверх или вниз по течению.
Пламя распределено не равномерно по занимаемому объему:
между кромкой факела и стенкой диффузора формируется воздушная прослойка, защищающая стенки от прогарания; в Рис. 18. Снимок факела эжектора центральной области также видно наличие воздушного слоя. Таким образом, факел имеет тороидальную форму. При проведении эксперимента изучались эмиссионные характеристики на выходе из эжектора. Выбросы по NOx лежат в пределах 7 ррm, по СО 10 ppm.
Максимальное расхождение между результатами численного и экспериментального исследований вихревого эжектора составило 35 %. В среднем различие расчетных и экспериментальных данных не превышает 20 %. При этом все результаты численного исследования при различных значениях геометрических и режимных параметров качественно верно описывают экспериментально установленные зависимости. Лучшее совпадение результатов численного исследования с экспериментом дает SST - k - модель турбулентности. Максимальное количественное расхождение в результатах наблюдается при исследовании распределения статического давления вдоль оси эжектора: область пониженного давления описывается достаточно точно, но при рассмотрении области относительной длины эжектора lc 1 проявляются существенные противоречия в результатах численного моделирования и эксперимента.
В пятой главе представлена методика расчета и спроектрованные с е использованием эжекционные горелочные устройства.
Инфракрасная газовая горелка предназначена для отопления бытовых и производственных помещений, термической обработки материалов в авиационной, машиностроительной, химической и других отраслях и позволяет сэкономить до 50-60% тепла.
Для предварительных испытаний изготовлен упрощенный опытный образец без рубашки охлаждения и керамической камеры. По результатам экспериментального исследования (рис. 19) установлено, что эмиссионные выбросы по СО и NOx составляют соответственно 10 и 7 ppm, что не превышает существующих норм. Уровень шума находится в диапазоне по частоте 10Ц 2700 Гц и интенсивРис. 19. Снимок инфракрасной эжекционной горелки ности колебаний не более 70 дБ, что также соответствует существующим нормам.
Устройство для сжигания топлива (рис. 20) предназначено для подготовки и сжигания топливовоздушной смеси и может найти применение в авиационных двигателях, теплоэнергетических установках и т.д.
Рис. 20. Снимки факелов богатого и бедного режимов горения На выходе из горелочного устройства четко видна зона обратных токов достаточно симметричной формы (рис. 20). Эмиссионные характеристики на выходе из низконапорной вихревой горелки находятся в пределах нормы:
выбросы по NOx 7 ррm, по СО 10 ppm.
ОСНОВНЫЕ ВЫВОДЫ Определены режимы работы прямоточного вихревого эжектора, удовлетворяющие требованиям качественной смесеподготовки и характеризующиеся отсутствием пульсаций расходов окислителя и топлива.
В результате анализа экспериментальных и расчетных данных выполнена количественная оценка неравномерности состава топливовоздушной смеси в выходном сечении эжектора-смесителя, которая не превышает 20%.
Эмиссионные выбросы, полученные по результатам экспериментального исследования, соответствуют СО 10 ppm; NО 7 ppm.
х Полученная система критериальных уравнений позволяет на начальной стадии оценить влияние геометрических и режимных параметров на эффективность рабочего процесса эжекционного устройства.
Экспериментальные данные по исследованию акустических характеристик работы вихревого эжектора позволили выявить наиболее эффективные геометрические и режимные параметры, удовлетворяющие требованиям по шуму и не снижающие эффективности процессов эжектирования, что обеспечит безопасное включение его в технологические цепочки различного назначения.
На основе развитых в работе представлений предложена методика расчета основных характеристик эжекционных устройств, с использованием которой рассчитаны и спроектированы конструкции инфракрасной газовой горелки и горелочного модуля, работающие по механизму обедненной предварительно перемешанной топливовоздушной смеси.
Основные публикации 1. Иванов, Р.И. Экспериментальное исследование акустических колебаний в вихревом эжекторе [Текст] / Р.И. Иванов, Ш.А. Пиралишвили, М.Н. Сергеев // Вестник Рыбинского государственного технического университета имени П.А. Соловьева. - 2011.Ц №2 (20).Ц С. 201-204.
2. Пиралишвили, Ш.А. Разработка инфракрасного газового горелочного устройства на базе вихревого эжектора [Текст] / Ш.А. Пиралишвили, А.И.Гурьянов, Р.И.Иванов // Вестник Самарского государственного аэрокосмического университета им. С. П. Королева. - 2007. - № 2 (13). - С. 151-154.
3. Пиралишвили, Ш.А. Исследование вихревого эжектора - смесителя [Текст] / Ш.А. Пиралишвили, Р.И.Иванов // Научно-технический и информационно-аналитический журнал Тепловые процессы в технике. - 2010. - №2, Т.2. - С. 57-62.
4. Пиралишвили, Ш.А. Смесеобразование в низкоперепадном вихревом эжекторе [Текст] / Ш.А. Пиралишвили, Р.И.Иванов // Материалы I Международной научно-технической конференции, посвященной 70 летию основателя Рыбинской школы теплофизиков доктора технических наук, профессора Пиралишвили Ш.А. Энергетические установки: тепломассообмен и процессы горения. - Рыбинск. - 2009. - С. 93-5. Пиралишвили, Ш.А. Численное и экспериментальное исследование рабочего процесса вихревого эжектора-смесителя [Текст]/ Ш.А. Пиралишвили, Р.И.Иванов // Труды XVII Школы-семинара молодых ученых и специалистов под руководством академика РАН А.И. Леонтьева Проблемы газодинамики и тепломассообмена в аэрокосмических технологиях. - М.:
Изд-во МЭИ, 2009 г. - Т. 2. - С. 60-64.
6. Пиралишвили, Ш.А. Исследование процесса смесеобразования в прямоточном вихревом эжекторе [Текст] / Ш.А. Пиралишвили, Р.И.Иванов // Труды пятой Российской национальной конференции по теплообмену. - М.: Изд-во МЭИ, 2010 г. - Т.3. - С. 214-27. Пат. 2413131 Российская Федерация МПК F 23 D 14/12 Инфракрасная газовая горелка [Текст] / Пиралишвили Ш.А., Гурьянов А.И., Иванов Р.И.; заявитель и патентообладатель Государственное образовательное учреждение высшего профессионального образования Рыбинская государственная авиационная технологическая академия имени П.А. Соловьева - заявка №2009137759/06, 12.10.2009; опубл. 27.02.11, Бюл. № 6. - 5 с.