Исследование аэродинамики в реверсивной топочной камере

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

Full description

Saved in:
Bibliographic Details
Published in:Промышленная теплотехника
Date:2012
Main Authors: Демченко, В.Г., Дуняк, О.В.
Format: Article
Language:Russian
Published: Інститут технічної теплофізики НАН України 2012
Subjects:
Online Access:https://nasplib.isofts.kiev.ua/handle/123456789/60266
Tags: Add Tag
No Tags, Be the first to tag this record!
Journal Title:Digital Library of Periodicals of National Academy of Sciences of Ukraine
Cite this:Исследование аэродинамики в реверсивной топочной камере / В.Г. Демченко, О.В. Дуняк // Промышленная теплотехника. — 2012. — Т. 34, № 6. — С. 46-52. — Бібліогр.: 6 назв. — рос.

Institution

Digital Library of Periodicals of National Academy of Sciences of Ukraine
_version_ 1859830734804484096
author Демченко, В.Г.
Дуняк, О.В.
author_facet Демченко, В.Г.
Дуняк, О.В.
citation_txt Исследование аэродинамики в реверсивной топочной камере / В.Г. Демченко, О.В. Дуняк // Промышленная теплотехника. — 2012. — Т. 34, № 6. — С. 46-52. — Бібліогр.: 6 назв. — рос.
collection DSpace DC
container_title Промышленная теплотехника
description Представлены результаты экспериментального исследования скорости и температуры топочных газов при наличии внутренней рециркуляции. Предложена математическая модель для определения осевых и радиальных скоростей, позволяющая определить значения коэффициента внутренней рециркуляции по геометрическим параметрам топочной камеры. Проведено сравнение полученных экспериментальных данных с результатами компьютерного и численного моделирования. Представлено результати експериментального дослідження швидкості і температури топкових газів за наявності внутрішньої рециркуляції. Запропонована математична модель для визначення осьових і радіальних швидкостей, що дозволяє визначити значення коефіцієнта внутрішньої рециркуляції по геометричних параметрах топкової камери. Проведено порівняння отриманих експериментальних даних з результатами комп'ютерного і чисельного моделювання. Results of an experimental research of speed and temperature of top internal gases in the presence of internal re-circulation are presented. The mathematical model for definition of the axial and radial speeds is offered, allowing defining values of factor of internal re-circulation on geometrical parameters of a fire chamber. Comparison of the received experimental data with results of computer and numerical modeling is spent.
first_indexed 2025-12-07T15:32:05Z
format Article
fulltext ISSN 0204-3602. Пром. теплотехника, 2012, т. 34, №646 УДК. 536.24:697.326 Демченко В. Г., Дуняк О. В. Институт технической теплофизики НАН Украины ИССЛЕДОВАНИЕ АЭРОДИНАМИКИ В РЕВЕРСИВНОЙ ТОПОЧНОЙ КАМЕРЕ Представлено результати ек- спериментального дослідження швидкості і температури топко- вих газів за наявності внутрішньої рециркуляції. Запропонована ма- тематична модель для визначен- ня осьових і радіальних швид- костей, що дозволяє визначити значення коефіцієнта внутрішньої рециркуляції по геометричних па- раметрах топкової камери. Про- ведено порівняння отриманих експериментальних даних з ре- зультатами комп'ютерного і чи- сельного моделювання. Представлены результаты экс- периментального исследования ско- рости и температуры топочных га- зов при наличии внутренней ре- циркуляции. Предложена математи- ческая модель для определения осе- вых и радиальных скоростей, по- зволяющая определить значения коэффициента внутренней рецирку- ляции по геометрическим параме- трам топочной камеры. Проведено сравнение полученных эксперимен- тальных данных с результатами ком- пьютерного и численного моделиро- вания. Results of an experimental research of speed and temperature of top internal gases in the presence of internal re-circulation are presented. The mathematical model for definition of the axial and radial speeds is offered, allowing defining values of factor of internal re-circulation on geometrical parameters of a fire chamber. Comparison of the received experimental data with results of computer and numerical modeling is spent. ИСПОЛЬЗОВАНИЕ И СЖИГАНИЕ ТОПЛИВА Создание и модернизация топочных камер котлов и новых способов сжигания органиче- ских топлив имеют важное научное и практи- ческие значение. Передача теплоты от факела к поверхностям нагрева топочной камеры яв- ляется одним из наиболее сложных видов теп- лообмена. Теплообмен усложнён тем, что он протекает параллельно с процессом горения, который создаёт в излучающей топочной ат- мосфере внутренние источники теплоты. Соот- ношение между мощностью источника тепло- выделения и интенсивностью теплоотдачи к поверхностям топочной камеры определяет как уровень температуры топочных газов, так и ха- рактер их изменений. Условия горения, теплообмена и теплового напряжения топочного объёма зависят от спо- соба сжигания топлива. Современные моно- блочные горелочные устройства обеспечивают высокую степень предварительного смешения топлива с воздухом. Это приводит к увеличе- нию температуры, при уменьшении длины и степени черноты факела. В промышленных ус- тановках интенсивность горения в топочной камере оценивают по величине θv – удельного тепловыделения, отнесённого к единице объ- ёма системы, кВт/м3: θv = B·Qн/V, (1) где B – расход топлива, кг/с; Qн – теплотворная низшая способность топлива, кДж/кг; V – объ- ём топочной камеры, м3. Из этого следует, что модернизации котлов должна быть направлена на оптимизацию то- почного объёма, что обеспечит локализацию и устойчивость реакций горения, создание опти- а, в – безразмерные параметры; B – расход топлива; Q – теплотворная способность топлива; r – радиус; V – объём топочной камеры; Р – давление; Φ – безразмерный корректирующий коэффи- циент; α, β, w – безразмерные параметры топочной камеры; ε – коэффициент рециркуляции; θ – удельное тепловыделение; υ – скорость потока; ρ – плотность газа. Индексы нижние: н – низшая; r – радиальная; дин – динамическое. ISSN 0204-3602. Пром. теплотехника, 2012, т. 34, №6 47 ИСПОЛЬЗОВАНИЕ И СЖИГАНИЕ ТОПЛИВА мальных условий их протекания и поддержа- ние наиболее выгодных режимов эксплуатации позволяющих получить максимальный КПД топочной камеры и снизить вредные выбросы в атмосферу. Данному вопросу посвящено много работ [1], но до сего времени по существующим дан- ным нельзя разработать соответствующие ре- комендации и нормативы. Известно, что интенсивность горения, опре- деляемая по скорости расхода горючих веществ, зависит не только от скорости протекания хи- мической реакции, но и от скорости процесса смесеобразования, определяющим фактором которой является интенсивность протекания турбулентной и молекулярной диффузии. Поэ- тому, одним из методов интенсификации топоч- ного теплообмена является достижение макси- мальной степени черноты топочной камеры [2], которая зависит от степени черноты факела и степени загрязнения экранов. Данное условие обеспечивает организация внутренней рециркуляции дымовых газов в то- почной камере котла и изменение её конструк- ции. Характер температурного поля в попереч- ном разрезе топочной камеры сильно влияет на условия теплообмена между факелом и экран- ными поверхностями нагрева. При высокой степени неизотермичности более холодные пристенные слои топочных газов, которые име- ют повышенную поглощающую способность, могут экранировать излучение центральных более горячих зон факела, что приводит к сни- жению плотности потока излучения, которое падает на поверхность нагрева, поэтому увели- чение степени черноты факела может привести к снижению интенсивности теплообмена в то- почной камере [3]. При этом исследованию аэродинамических процессов протекающих в реверсивных топоч- ных камерах при внутренней рециркуляции и их влиянию на интенсификацию топочного те- плообмена уделено по нашему мнению не до- статочно внимания. При сжигании газа и дизельного топли- ва в топочных камерах котлов большую роль играют размещенные в топочных камерах ог- неупорные излучающие насадки, горки или перегородки. Эти так называемые вторичные излучатели значительно повышают тепловос- приятие радиационных поверхностей нагре- ва. В нормативном методе теплового расчета котельных агрегатов расчет теплообмена в топочных камерах, имеющих вторичные из- лучатели, не рассматривается и рекомендации по их применению отсутствуют, недооценка роли вторичных излучателей может привести к большим ошибкам при эксплуатации, расчете и конструировании топочных камер. Известно, что прямая отдача теплоты в то- почных камерах зависит не только от темпера- туры, но и от степени светимости излучателя. Любое нагретое тело излучает тепловую энер- гию в инфракрасном диапазоне спектра элек- тромагнитных волн и передает её лучистым теплообменом другим телам. Излучающее тело отдает свою энергию в виде теплоты, вслед- ствие чего его температура уменьшается и на- оборот, если поверхность какого-нибудь тела подвергается тепловому излучению, то тело, поглощая энергию, увеличивает свою темпера- туру. Наибольший эффект в этом случае дают вторичные излучатели – раскалившись, они из- лучают теплоту на поверхности нагрева. На- пример, в печах теплоотдача экранным по- верхностям от вторичных излучателей весьма значительна и соизмерима с теплоотдачей из- лучением от факела и газовой среды [4]. Постановка задачи Внутренняя рециркуляция топочных газов может быть обеспечена путём изменения гео- метрии топочной камеры при установке в неё вторичных излучателей, разработанных ИТТФ НАНУ [5]. В этом случае меняется аэродина- мика топочной камеры, появляются зоны пря- мых и обратных токов, зоны разделения пото- ков и пр. Особый интерес представляет определение скоростей и температур топочных газов в зоне бифуркации, в кольцевом канале и щелевом за- зоре образованном стенками жаровой трубы и вторичного излучателя. ISSN 0204-3602. Пром. теплотехника, 2012, т. 34, №648 ИСПОЛЬЗОВАНИЕ И СЖИГАНИЕ ТОПЛИВА Для проведения исследования нами был выбран серийно выпускаемый водогрейный жаротрубный котёл «Виктор-80», мощностью 80 кВт (рис. 1). Котёл имеет реверсивную топочную ка- меру образованную жаровой трубой. Размеры топочной камеры позволяют установить вто- ричный излучатель в форме продольно ореб- ренной трубы. Котёл оборудован моноблочной газовой либо дизельной горелкой расположен- ной на фронтальной стенке по центральной оси топочной камеры. Исследования проводились при работе котла на дизельном топливе. При стационарной работе котла, как по- казывает компьютерное моделирование [6] в нижней и верхней области факела образуются устойчивые вихревые потоки, которые направ- ляют факел вверх, образуя тем самым зону ло- кального перегрева в верхней части жаровой трубы. Топливовоздушная смесь попадает в топоч- ную камеру котла через сопло вентиляторной горелки. Горение протекает в объёме топочной камеры с распространением факела по её дли- не. Топочные газы, ударяясь о заднюю стенку топки, теряют свою скорость, поворачиваются на 180° и возвращаются на фронт котла. Часть топочных газов попадает в корень факела, а больший массовый объём в сборную коробку дымовых газов и оттуда через конвективный пучок образованный дымогарными трубами с турбулизаторами потока в сборную коробку дымовых газов и далее в дымовую трубу. После установки вторичного излучателя в топочную камеру движение топочных газов в котле коренным образом изменяется. После сгорания в объёме вторичного излучателя то- почные газы попадают в поворотную камеру, образованную торцевой стенкой топочной ка- меры и вторичным излучателем. Оттуда, повер- нувшись на 180°, они направляются на фронт котла по кольцевому каналу образованному стенками жаровой трубы и вторичного излуча- теля. На выходе из кольцевого канала продук- ты сгорания разделяются на два потока – один поступает в конвективный пучок и далее, как описано ранее, а второй возвращается в каме- ру сгорания, образованной вторичным излуча- телем. Поле скоростей в этой области можно принять ассиметричным, а форму средних тра- екторий частиц продуктов сгорания достаточно ограничить кривыми второго порядка. Уравне- ние неразрывности потока решается при вве- дении следующих упрощений: распределение плотности газа ρ1, ρ2, ρ3 усредняется по трём областям: - на выходе из кольцевого канала; - в области между жаровой трубой и конвектив- ным пучком; - в середине трубного пучка. Граничные условия отображают отсутствие радиальной скорости на стенках и линии раз- деления потоков. Зависимость для определения компонент скоростей (r1 < r < r2; 0 < z < z0) полу- чены в виде (2): ( ) ( ) ( ) ( )( ) 2 0 22 2 2 2 0 0 0 ln − ϑ = ⋅ π − − − − r b z zq r r r r a z z r r ; Рис. 1. Принципиальная схема организации движения топочных газов в котле «ВИКТОР-80», при установке вторичного излучателя: A – горелка; B – вторичный излучатель, образующий камеру сгорания; C – поворотная камера; D – поворотная камера между фронтальной трубной доской котла и изолированной огнеупором дверкой; E – конвективный пучок (дымогарные трубы); F – сборная коробка дымовых газов; T1 – подающий трубопровод; T2 – обратный трубопровод; 1, 2, 3, 4 – точки замеров динамических давлений и температур. ISSN 0204-3602. Пром. теплотехника, 2012, т. 34, №6 49 ИСПОЛЬЗОВАНИЕ И СЖИГАНИЕ ТОПЛИВА ( ) ( ) ( ) ( )( ) 2 0 2 22 2 2 0 02 0 0 ln − ϑ = ⋅ − − −π − z b r rq r r a z z r rr r r . (2) где a, b – безразмерные параметры, зависящие от геометрии топочного пространства и коор- динат размещения в ней вторичного излучателя. С помощью мероморфной функции, у кото- рой количество полюсов равно количеству ды- могарных труб, можно найти зависимость для определения компонент скоростей газового по- тока на входе в дымогарную трубу в полярных координатах, которая имеет вид (3): ( ) ( )( ) T 2 20 2 0 1, , 2 cos   ϑ ϕ = ϑ − + − ϕ−ϕ    z k k k k zzr z r r r r z z r (3) Она была получена при следующих гранич- ных условиях: - все компоненты скоростей на стенках равны нулю; - поток радиальной скорости через цилиндри- ческую поверхность r = r2 (внутренний радиус жаровой трубы) равен потоку осевой скорости через торцевое сечение всех дымогарных труб; - поток радиальной скорости υr через боковую поверхность конического сечения осевых ско- ростей равен потоку осевой скорости υz через его основание. Аналитический метод расчета компонентов скоростей позволяют определить температур- ное распределение и должны рассматриваться совместно при решении задач связанных с те- пломассопереносом численными методами. Радиальная координата r0 полностью опре- деляет коэффициент рециркуляции ε топочных газов идущих на повторное сжигание, по фор- муле (4). Коэффициент рециркуляции показы- вает массовый объём топочного газа в кольце- вом канале, возвращаемый в камеру сгорания и определяется по формуле: 0 1 1 2 / .ε = ρϑ ρϑ∫ ∫ r r z z r r rdr rdr (4) Данная зависимость получена при условии невлияния изменений плотности на скорость газового потока. Введя гипотетический радиус разделения газового потока r0 (рис. 2), в зоне бифуркации можно получить зависимость для определения коэффициента рециркуляции в виде показательно-степенной функции и пу- тём решения обратной задачи определить не- известные параметры (5): ( ) α −βε = Φ ⋅ = Φ hf h wh e , (5) где: α, β, w – безразмерные параметры топоч- ной камеры, которые определяются произволь- ными опорными точками (εi, hi); Φ – безразмер- ный корректирующий коэффициент, который зависит от соотношения радиусов r1 и r2. Экспериментальное исследование Для определения реальных теплофизичес- ких процессов происходящих в реверсивной топочной камере при внутренней рециркуля- ции топочных газов проведены замеры скоро- сти и температуры в топочных газов в четырёх реперных точках (рис. 3) численные значения которых, характеризуют изменения аэродина- мики при наличии внутренней рециркуляции. Измерения динамических давлений прово- дились пневмометрическими трубками, спе- циально разработанными отделом ПТТ ИТТФ НАНУ для высокотемпературных потоков. Трубки подключались к многофункционально- му измерительному прибору Testo 435-4. Поправочный коэффициент пневмометри- ческих трубок К = 0,73 (отношение истинного значения к значению замеряемого динамичес- кого давления). Торрировка пневмометричес- ких трубок проводилась на специально создан- ном стенде. Измерение температуры проводи- лось термопарами, подключенными к многока- нальному аппаратно-программному комплексу «Triton 6004 ТС». В процессе проведения замеров были от- мечены резкие колебания давлений, вызванные динамическими нагрузками, возникающими при работе котла. Полученные данные усред- нены. Температура дымовых газов в точках за- мера практически постоянна. Данные замеров представлены в таблице 1. При сопоставлении результатов получен- ных при моделировании с результатами, полу- . ISSN 0204-3602. Пром. теплотехника, 2012, т. 34, №650 ИСПОЛЬЗОВАНИЕ И СЖИГАНИЕ ТОПЛИВА а) б) Рис. 3. Установка пневмометрических трубок в топочную камеру котла: а) Рдин.1 и Рдин.2; б) Рдин.3 и Рдин.4. а) б) Рис. 2. Схема движения газовых потоков в топочной камере с вторичным излучателем: а) аксонометрия расчётной области; б) вертикальное сечение топочной камеры по центральной оси горелочного устройства; 1 – жаровая труба; 2 – устье горелки; 3 – топливная дюза; 4 – щелевой зазор; 5 – вторичный излучатель; 6 – торцевая поворотная камера; Lт – длина топочной камеры, м; Lм – длина торцевой пово-ротной камеры, м; Lи – длина вторичного излучателя, м; h – длина фронтального щелевого зазора, м; r1 – радиус вторичного излучателя, м; r2 – радиус жаровой трубы, м; r0 – радиус бифуркации, м. ISSN 0204-3602. Пром. теплотехника, 2012, т. 34, №6 51 ИСПОЛЬЗОВАНИЕ И СЖИГАНИЕ ТОПЛИВА ченными экспериментальным путем установ- лено, что относительная погрешность в точках замеров даёт сходимость по скорости – до 15 % и по температуре до 3 %. После установки вторичного излучателя температура на выходе из топочной камеры снижается на 18…28 %, что свидетельствует о повышении её КПД. Скорость потока топочных газов с вторич- ным излучателем на 39…47 % ниже, что сви- детельствует об увеличении времени пребы- вания топочных газов в камере сгорания и тор- можении потока щелевым зазором образован- ном стенками жаровой трубы и излучателя. Выводы 1. Анализ состояния теории и практики интенсификации теплопередачи в топочных камерах котлов выявил целесообразность при- менения для этой цели вторичных излучателей. Разработанная конструкция вторичного излу- чателя имеет цилиндрическую форму, изготав- ливается из нержавеющей стали с приварива- емыми к ее внешним поверхностям ребрами, для симметричного ее размещения по оси го- релочного устройства. Такая конструкция обе- спечивает появление дополнительного хода то- почных газов и возникновение рециркуляции газов в корень факела. 2. Разработанная математическая модель аэродинамики топочных газов внутри котла позволяет определять расход циркулирующих и транзитных топочных газовых потоков. Пред- ложенная методика расчёта позволяет рассчи- тать распределение расходов топочных газов по пучку дымогарных труб в конвективной ча- сти котла. 3. Скорости в зоне бифуркации потоков продуктов сгорания практически имеют оди- наковые значения. Это указывает на наличие рециркуляции продуктов сгорания в корень фа- кела. 4. Выполненные расчеты по разработан- ной аналитической модели свидетельствуют о том, что не менее 50 % топочных газов про- ходят повторный дожег. 5. Согласие результатов проведенных рас- четов по двум различным методикам свиде- тельствуют об адекватности разработанной аналитической модели и CFD расчета, что под- тверждается полученными эксперименталь- ными данными. ЛИТЕРАТУРА 1. Акилов В.А., Бридун Е.В., Ватачин М.Ю. и др. Актуальные проблемы устойчивого разви- тия. Применение новых технологий сжигания топлива. – Киев: О-во «Знание», 2003. – 430 с. 2. Сполдинг Д.Б. Основы теории горения. – М-Л.: Госэнергоиздат, 1959. – 320 с. 3. Блох А.Г. Теплообмен в топках паровых котлов. – Л.: Энергоатомиздат. Ленинградское отделение, 1984. – 240 с. Табл. 1. Результаты замеров значений динамических давлений и температур То чк а за ме ра Количество замеров Температура, °С Динамическое давление, мбар Скорость, м/с Бе з и зл уч ат ел я С и зл уч ат ел ем Бе з и зл уч ат ел я С и зл уч ат ел ем Бе з и зл уч ат ел я С и зл уч ат ел ем Бе з и зл уч ат ел я С и зл уч ат ел ем № 1 - 45 - 900 - 0,08 - 6,2 № 2 - 45 - 800 - 0,15 - 8,1 № 3 182 210 820 680 0,33 0,14 12,2 7,4 № 4 186 240 820 590 0,46 0,16 14,4 7,6 ISSN 0204-3602. Пром. теплотехника, 2012, т. 34, №652 ИСПОЛЬЗОВАНИЕ И СЖИГАНИЕ ТОПЛИВА 4. Кнорре Г.Ф. Теория топочных процес- сов. – М.: Энергия, 1966. – 134 с. 5. Демченко В.Г., Долинский А.А., Сигал А.И. Организация внутренней рециркуляции ды- мовых газов в реверсивных водоохлаждаемых топках с целью модернизации отопительных котлов. // NATO workshop conference «Advanced combustion and aerothermal technologies» Kiev 2006. 6. Басок Б.И., Демченко В.Г., Мартыненко М.П. Численное моделирование процессов аэ- родинамики в топке водогрейного котла с вто- ричным излучателем // Промышленная тепло- техника. – 2006. – №1. – С. 17 – 22. Получено 28.09.2012 г.
id nasplib_isofts_kiev_ua-123456789-60266
institution Digital Library of Periodicals of National Academy of Sciences of Ukraine
issn 0204-3602
language Russian
last_indexed 2025-12-07T15:32:05Z
publishDate 2012
publisher Інститут технічної теплофізики НАН України
record_format dspace
spelling Демченко, В.Г.
Дуняк, О.В.
2014-04-13T16:38:17Z
2014-04-13T16:38:17Z
2012
Исследование аэродинамики в реверсивной топочной камере / В.Г. Демченко, О.В. Дуняк // Промышленная теплотехника. — 2012. — Т. 34, № 6. — С. 46-52. — Бібліогр.: 6 назв. — рос.
0204-3602
https://nasplib.isofts.kiev.ua/handle/123456789/60266
536.24:697.326
Представлены результаты экспериментального исследования скорости и температуры топочных газов при наличии внутренней рециркуляции. Предложена математическая модель для определения осевых и радиальных скоростей, позволяющая определить значения коэффициента внутренней рециркуляции по геометрическим параметрам топочной камеры. Проведено сравнение полученных экспериментальных данных с результатами компьютерного и численного моделирования.
Представлено результати експериментального дослідження швидкості і температури топкових газів за наявності внутрішньої рециркуляції. Запропонована математична модель для визначення осьових і радіальних швидкостей, що дозволяє визначити значення коефіцієнта внутрішньої рециркуляції по геометричних параметрах топкової камери. Проведено порівняння отриманих експериментальних даних з результатами комп'ютерного і чисельного моделювання.
Results of an experimental research of speed and temperature of top internal gases in the presence of internal re-circulation are presented. The mathematical model for definition of the axial and radial speeds is offered, allowing defining values of factor of internal re-circulation on geometrical parameters of a fire chamber. Comparison of the received experimental data with results of computer and numerical modeling is spent.
ru
Інститут технічної теплофізики НАН України
Промышленная теплотехника
Использование и сжигание топлива
Исследование аэродинамики в реверсивной топочной камере
Research of aerodynamics of reversing combustion chamber
Article
published earlier
spellingShingle Исследование аэродинамики в реверсивной топочной камере
Демченко, В.Г.
Дуняк, О.В.
Использование и сжигание топлива
title Исследование аэродинамики в реверсивной топочной камере
title_alt Research of aerodynamics of reversing combustion chamber
title_full Исследование аэродинамики в реверсивной топочной камере
title_fullStr Исследование аэродинамики в реверсивной топочной камере
title_full_unstemmed Исследование аэродинамики в реверсивной топочной камере
title_short Исследование аэродинамики в реверсивной топочной камере
title_sort исследование аэродинамики в реверсивной топочной камере
topic Использование и сжигание топлива
topic_facet Использование и сжигание топлива
url https://nasplib.isofts.kiev.ua/handle/123456789/60266
work_keys_str_mv AT demčenkovg issledovanieaérodinamikivreversivnoitopočnoikamere
AT dunâkov issledovanieaérodinamikivreversivnoitopočnoikamere
AT demčenkovg researchofaerodynamicsofreversingcombustionchamber
AT dunâkov researchofaerodynamicsofreversingcombustionchamber