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

16

Размещено на http://www.allbest.ru/

Размещено на http://www.allbest.ru/

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

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

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

Специальность 01.04.14 - Теплофизика и теоретическая теплотехника

Василюк Ольга Владимировна

Рыбинск - 2008

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

Научный руководитель

заслуженный деятель науки и техники РФ, доктор технических наук, профессор Пиралишвили Шота Александрович

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

доктор технических наук, профессор Сапожников Сергей Захарович

кандидат технических наук, доцент Христофоров Игорь Леонидович

Ведущая организация

ОАО «НПО «Сатурн», г. Рыбинск

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

Актуальность работы.

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

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

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

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

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

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

Цель диссертационной работы.

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

Научная новизна работы заключается в том, что:

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

- оценить работу ВГУ при различных режимных параметрах на входе;

- выявить их влияние на местонахождение и величину локальных зон повышенной температуры;

- оценить характеристики процесса смесеобразования;

- улучшить характеристики распыла в закрученном потоке.

Достоверность и обоснованность достигается:

- корректным применением основных термогазодинамических законов;

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

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

- совпадением расчетных и экспериментальных зависимостей дисперсности распыла и полей полной температуры.

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

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

Личный вклад автора заключается:

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

- в создании экспериментального стенда для исследования температурных полей в тракте ВГУ, в проведении опытов по определению дисперсности распыла, температуры самовоспламенения топливовоздушной смеси;

- в проведении расчётных и экспериментальных исследований, обработке и анализе опытных данных;

- в разработке математической модели процесса смесеобразования.

На защиту выносятся следующие основные положения

1. Математическая модель, позволяющая рассчитать газодинамические параметры по объему горелочного устройства.

2. Математическая модель процесса смесеобразования в вихревой камере горелочного устройства.

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

Апробация работы

Основные результаты выполненных исследований докладывались и обсуждались на следующих конференциях:

- IV Всероссийская научно-техническая конференция «Теплофизика процессов горения и охрана окружающей среды», Рыбинск, 1999 г.;

- Х Всероссийская научно-техническая конференция «Теплофизика технологических процессов», Рыбинск, 2000 г.;

- ХIII школа-семинар п/р ак. Леонтьева А.И. «Физические основы экспериментального и теоретического моделирования процессов газодинамики и теплообмена», Санкт-Петербург, 2001г. (присужден диплом II степени);

- Всероссийская конференция «Закрутка потока для повышения эффективности теплообменников», ИВТ РАН, Москва, 2002 г.;

- Всероссийская научно-техническая конференция «Процессы горения, теплообмена и экология тепловых двигателей», Самара, 2002 г.;

- XIV школа-семинар п/р ак..Леонтьева А.И. «Проблемы газодинамики и тепломассообмена в энергетических установках», Рыбинск, 2003 г.;

- VI Всероссийская научно-техническая конференция «Теплофизика процессов горения и охрана окружающей среды», Рыбинск, 2004 г.;

- XX юбилейный международный семинар по струйным, отрывным и нестационарным течениям, Санкт-Петербург, 2004 г.

Публикации.

По материалам диссертации опубликовано 16 работ, из них 3 статьи в изданиях, утверждённых ВАК, 2 доклада в трудах конференций.

Структура и объём диссертации.

Диссертация состоит из введения, пяти глав, заключения, списка используемых источников. Работа изложена на 127 страницах машинописного текста, включая 1 таблицу, 61 рисунок. Список используемых источников включает 97 наименований.

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

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

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

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

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

Для обеспечения критического истечения продуктов сгорания из выходного сопла воспламенителя величина степени расширения воздуха в камере энергетического разделения будет определяться соотношением

, (1)

где - давление в запускаемой камере сгорания.

Среднемассовая температура воздуха в перфорированной камере воспламенителя

, (2)

где - относительный подогрев газа в перфорированной камере.

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

, (3)

где - газодинамическая функция расхода.

Диаметр камеры энергетического разделения в сечении соплового ввода

, (4)

противоточный вихревой горелка теплофизический

где - относительная площадь соплового ввода.

Тогда . Высота соплового ввода , где - число каналов соплового ввода. Ширина соплового ввода . Относительную длину вихревой камеры рекомендуется выбирать в пределах 3 9.

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

Несмотря на то, что в настоящее время накоплен обширный экспериментальный материал по исследованию газодинамики закрученных потоков, не существует единого мнения относительно механизма энергоразделения и температурной стратификации в таких течениях. Наиболее обоснованной считается гипотеза взаимодействия вихрей, предложенная А.П. Меркуловым и дополненная Ш.А. Пиралишвили. Последователями этой гипотезы являлись W. Frohlingsdorf, H. Unger, осуществившие численное моделирование феномена энергоразделения в вихревой трубе Ранка - Хилша с помощью кодов системы CFX.

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

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

С использованием программного пакета UG построены трехмерные модели исследуемой схемы вихревого энергоразделителя, на базе которых осуществлялась дискретизация (генерация сетки, рис.2) с помощью универсального сеткопостроителя HEXA.

Постановка задачи завершается заданием граничных и начальных условий в виде условий адиабатности и прилипания. На входе в вихревой энергоразделитель известными считались: полные давление и температура, на выходе - статическое давление. Шаг по времени составил 10^-4 - 10^-5 с.

Рис.1. Расчетная схема противоточного вихревого энергоразделителя

Рис. 2. Сеточная модель вихревого энергоразделителя

Проведенные расчеты подтвердили формирование в камере энергоразделения вихревой трубы двух вращающихся в одном направлении вихрей, перемещающихся в противоположных осевых направлениях, которые впервые экспериментально были обнаружены в опытах Shults-Grunow F. Проведенная компьютерная визуализация позволяет сделать вывод о том, что окружная компонента скорости заметно превосходит осевую. Причем, если осевая составляющая скорости периферийного вихря по мере удаления вдоль камеры энергоразделения от соплового ввода к дросселю изменяется незначительно, то затухание вращательного движения более существенно, что приводит к росту шага винтовых линий. Этот факт также использован в гипотезе взаимодействия вихрей. Именно поэтому радиальный градиент давления от соплового сечения к дросселю падает, а статическое давление на оси растет, вызывая появление осевого градиента давления, направленного от дросселя к оси диафрагмы. В то же время закрутка приосевых масс газа, покидающих камеру энергоразделения через отверстие диафрагмы заметно меньше, что может быть объяснено возрастанием осевой компоненты скорости вынужденного приосевого вихря в сопловом сечении вихревой трубы на выходе из отверстия диафрагмы. Один из самых важных результатов визуализации состоит в том, что впервые расчетным путем были получены крупномасштабные вихревые структуры сдвигового характера, предсказанные на основе косвенных опытов в работах Лукачева С.В. и экспериментально обнаруженные при двухфазной визуализации течения керосино-воздушной смесью. Эпюры распределения осевой и тангенциальной составляющих скоростей по характеру практически не отличаются от полученных экспериментально (рис. 3, 4). Сравнить профиль распределения радиальной составляющей скорости, полученный численным методом, с опытными данными не представляется возможным, так как экспериментальное определение этой составляющей скорости связано с трудностями, вызванными существенными радиальными градиентами термогазодинамических параметров, что усложняет само измерение и заметно снижает его точность. Незначительное отклонение осевого положения датчика может существенно изменить его показание на противоположное, если измерения осуществлять вблизи дислокации крупномасштабной вихревой структуры.

16

Размещено на http://www.allbest.ru/

Размещено на http://www.allbest.ru/

Рис. 3. Профили осевой скорости в сопловом сечении

расчет 1 -; расчет 2 -; расчет 3 -; эксперимент [57]

Рис. 4. Профили тангенциальной скорости в сечении

расчет 1 -; расчет 2 -; расчет 3 -; эксперимент [57]

Проведенное сравнение показало возможность использования составленной модели для расчета закрученных течений с погрешностью, не превышающей 15 %. Поэтому составленная математическая модель применена для расчета газодинамики течения в вихревом горелочном устройстве (рис. 5).

По результатам расчета построены поля температуры, давления и скоростей. Выявлено, что область максимальной температуры меняет свое месторасположение от середины камеры энергоразделения до конца перфорированной камеры в зависимости от режимных параметров на входе: давления и температуры. Причем с повышением давления на входе повышается максимальная температура (при = 0,2 МПа = 1,1; = 0,3 МПа = 1,11; = 0,4 МПа = 1,15; = 0,5 МПа = 1,18; = 0,6 МПа = 1,3).

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

Проведенный расчет позволил изучить структуру течения и получить распределение газодинамических параметров потока по объему горелочного устройства. Численный расчет показал, что осредненная скорость потока по тракту горелки сильно снижается (с 472 м/с в сечении завихрителя до 24 м/с внутри перфокамеры), следовательно можно уменьшить длину горелки, сократить ее габариты и при этом избежать образования застойных зон.

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

. (5)

Уравнение замыкается уравнениями тепло- и массообмена для сферической капли и условиями кинематики.

Расчет прогрева и испарения капель топлива показал, что крупные капли (с диаметром > 18,5 мкм) подвергаются вторичному дроблению (критерий Вебера We = 14), капли с диаметром 5 15 мкм вылетают из горелки и являются стабилизаторами пламени, а большинство капель (с диаметром 15 мкм) отбрасывается на стенку, испаряясь с нее, образуя спиральные жгуты, которые обеспечивают надежное пленочное охлаждение, гарантирующее длительный ресурс при высокой теплонапряженности.

Четвёртая глава посвящена экспериментальным исследованиям дисперсности распыла методом малоуглового рассеяния света. Сравнение расчетных и экспериментальных данных показало, что суммарная погрешность определения V_s(d) в закрученном потоке методом малоуглового рассеяния света не превышает 5% (рис. 10). Проведенный эксперимент по исследованию дисперсности распыла топлива показал, что качество распыла внутри вихревого горелочного устройства улучшается (дисперсность распыла по сравнению с незакрученным потоком увеличилась с 80 до 20 мкм, что объясняется более лучшим дроблением и испарением капель топлива в закрученном потоке).

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

Температура движущихся закрученных масс воздуха измерялась в 5 сечениях = 1,8; 3,6; 5,3; 6,7; 8,9 от стенки до оси камеры энергетического разделения. В плоскости выше среза сопла перфорированной камеры (= 6,7) наблюдались максимальные эффекты подогрева ( = 1,51) при =0,86, =0,22. Таким образом, в этом месте возможно самовоспламенение топливовоздушной смеси (рис. 11). Причем при использовании завихрителя с = 0,22 на всем протяжении камеры энергоразделения достигаются максимальные эффекты подогрева, которые имеют минимальное значение в сечении = 1,8. Остальные завихрители дают гораздо меньший эффект подогрева. Самый маленький эффект подогрева достигается при использовании завихрителя с относительной площадью проходного сечения = 0,05 (= 1,23). С увеличением давления на входе эффекты подогрева возрастают (рис. 12). Эта зависимость прослеживается на всех значениях , во всех сечениях.

Рис. 13. Схема экспериментального стенда

1 - компрессор; 2, 3, 18, 19 - вентили; 4 - ССУ; 5 - U - образный манометр; 6, 15, 21 - пружинный манометр; 7, 13, 16 - термопара ХК; 8 - объект исследования; 9 - термопара ХА; 10 - антипомпажная магистраль; 11 - топливный бак; 12 - мерная трубка; 14 - электронагреватель; 17 - секундомер; 20 - топливный насос

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

. (6)

Эксперимент показал более высокое значение температуры теплового возгорания по сравнению с расчетным прогнозом (рис. 14). Это объясняется тем, что часть энергии расходуется на испарение топлива, вводимого в жидкой фазе. Границы самовоспламенения лежат в достаточно узком диапазоне по коэффициенту избытка воздуха. Самовоспламенение наблюдалось при суммарном коэффициенте избытка воздуха =0,9 0,1. Эти результаты позволили рассмотреть возможность практического применения вихревой горелки для создания воспламенителя камеры сгорания безискрового розжига.

16

Размещено на http://www.allbest.ru/

Размещено на http://www.allbest.ru/

Рис. 14. Граница самовоспламенения

- эксперимент; граница самовоспламенения по [96] (в потоке воздуха); граница самовоспламенения по [97] (в неподвижном воздухе)

ЗАКЛЮЧЕНИЕ

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

На основании выполненных расчетных и экспериментальных исследований можно сделать следующие выводы:

1. Создана методика расчета газодинамических характеристик ВГУ (давления, температуры, скоростей, чисел Маха и т. д.) с помощью численного метода, позволяющая рассчитать все значения отмеченных параметров на различных режимах работы по всему объему ВГУ с точностью до 10 %.

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

3. Разработана математическая модель процесса смесеобразования в ВГУ, позволяющая улучшить характеристики распыла в закрученном потоке. Дисперсность распыла уменьшилась с 80 до 20 мкм. Большинство капель (с диаметром 18,5 мкм) подвергаются вторичному дроблению, а капли с диаметром 15 мкм отбрасывается на стенку и испаряются с нее, образуя спиральные жгуты, которые обеспечивают надежное пленочное охлаждение, гарантирующее длительный ресурс при высокой теплонапряженности.

4. Экспериментальное исследование поля температуры внутри вихревой горелки выявило, что наибольшие эффекты подогрева (более 50 %) достигаются при относительной площади соплового ввода = 0,22 и относительном радиусе отверстия сопла-диафрагмы = 0,86.

ОСНОВНЫЕ ПУБЛИКАЦИИ ПО ТЕМЕ ДИССЕРТАЦИОННОЙ РАБОТЫ

1 Казанцева, О. В. Исследование процесса смесеобразования в вихревом воспламенителе [Текст] / О. В. Казанцева, Ш. А. Пиралишвили // Труды XIII Школы-семинара молодых ученых и специалистов под руководством академика РАН А. И. Леонтьева «Физические основы экспериментального и теоретического моделирования процессов газодинамики и теплообмена». - Москва: изд-во МЭИ, 2001. - т.1. - с. 232 - 234.

2 Казанцева, О. В. Исследование смесеобразования в вихревом воспламенителе [Текст] / О. В. Казанцева, Ш. А. Пиралишвили, Н. П. Лякина // Известия Академии наук «Энергетика». -2002. - №5. - с.162 - 166. (перечень ВАК).

3 Казанцева, О. В. Численное моделирование и расчет вихревого термотрансформатора [Текст] / О. В. Казанцева, Ш. А. Пиралишвили, Д. К. Василюк, М. В. Медведева // Труды XIV школы семинара молодых ученых и специалистов п/р акад. РАИ А. И. Леонтьева «Проблемы газодинамики и тепломассообмена в энергетических установках». - Рыбинск, 2003. - т.1. - с. 73 - 76.

4 Василюк, Д. К. Численное моделирование закрученных течений в вихревых трубах [Текст] / Д. К. Василюк, О. В. Казанцева, Ш. А. Пиралишвили // Материалы Международной научно-технической конференции «Проблемы и перспективы развития двигателестроения». - Самара: СГАУ, 2003. - с. 268 - 271.

5 Казанцева, О. В. Влияние испаряемости топлива на характеристики распыла [Текст] / О. В. Казанцева // Материалы IV Всероссийской научно-технической конференции «Теплофизика процессов горения и охрана окружающей среды». - Рыбинск: РГАТА, 2004. - с. 17 - 23.

6 Гурьянов, А. И. Вихревые горелочные устройства [Текст] / А. И. Гурьянов, О. В. Казанцева, М. В. Медведева, Ш. А. Пиралишвили // Справочник. Инженерный журнал. - 2005. - 43, № 5. - С. 8 - 15. (перечень ВАК).

7 Казанцева, О. В. Численное моделирование закрученных течений в вихревых трубах [Текст] / О. В. Казанцева, Ш. А. Пиралишвили, А. А. Фузеева //Теплофизика высоких температур. - 2005. - 43, № 4. - с. 606 - 611. (перечень ВАК).

8 Василюк, Д. К. Математическое моделирование структуры потока в трубах Ранка [Текст] / Д. К. Василюк, О. В. Казанцева, Ш. А. Пиралишвили // Материалы Всероссийской научно-технической конференции «Теплофизика технологических процессов». - Рыбинск: РГАТА, 2005. - С. 220 - 222.

Размещено на Allbest.ru