Теплоэнергетика, 2021, № 11, стр. 62-70

Анализ применения FGM- и EDM-моделей горения в горелках с частичным предварительным смешиванием для систем отопления, вентиляции и кондиционирования воздуха

N. Jaggu a*, S. R. Alluru a, A. Balakrishna a, V. Kamasetty a

a Trane Technologies Engineering and Technology Centre
560029 Bangalore, Karnataka, Bannerghatta Road, IBC Knowledge Park, 4/1, Tower D, 8th Floor, India

* E-mail: NJaggu@tranetechnologies.com

Поступила в редакцию 26.05.2020
После доработки 18.06.2020
Принята к публикации 26.08.2020

Полный текст (PDF)

Аннотация

Топочные устройства на природном газе являются основным источником получения тепла для систем отопления, вентиляции и кондиционирования воздуха. Использование в них горелок с частичным предварительным смешиванием приобретает в последнее время все большее значение. Использование методов вычислительной гидродинамики (CFD-методов ‒ от англ. computational fluid dynamics) для исследования характеристик пламени и теплопередачи в этих устройствах ограничено. Достоверное моделирование турбулентности в разных камерах сгорания всегда было большой проблемой для инженеров-газодинамиков. В частности, в топочных устройствах, сжигающих газ, необходимо учитывать значительное количество компонентов турбулентного потока, что затрудняет моделирование. В этой связи важно проанализировать различные методы моделирования горения, чтобы предсказать поведение пламени при горении, его структуру, длину факела, температуру и прочие факторы в области быстрых химических реакций. В данной работе предпринята попытка сравнить две технологии моделирования: с использованием простого универсального метода вихревой диссипации (EDM ‒ от англ. eddy dissipation method) и детальную химическую модель коллектора газовых горелок (FGM ‒ от англ. flamelet generated manifold). В большинстве отраслей эти методы моделирования применяются в зависимости от поставленных задач. Влияние соотношения содержания воздуха и топлива на температуру пламени, массовую долю CO2 и КПД топочного устройства было проанализировано при использовании EDM- и FGM-моделей. Исследование показало лучшую корреляцию результатов испытаний при расчетах температуры воздуха и дымовых газов и КПД с применением FGM‑модели по сравнению с результатами расчетов с использованием EDM-модели.

Ключевые слова: моделирование, FGM, EDM, топочное устройство, топливовоздушная смесь, процессы горения, CFD, теплоперенос, температура дымовых газов

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

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

В модели коллектора газовых горелок (FGM‑модель) применяются методы химической редукции пламени для точного моделирования состава газа. В данной статье модели EDM и FGM рассмотрены в контексте сжигания предварительно смешанных газов.

МОДЕЛЬ ВИХРЕВОГО РАССЕЯНИЯ (EDM-MODEL)

При реализации модели вихревого рассеяния предполагается, что скорость реакции определяется турбулентностью без учета времени химического взаимодействия. Для использования EDM-модели вводятся следующие предположения. Скорость горения довольно велика. Она определяется числом Дамкёлера Da [1, 2] и моделью турбулентного горения в газовой фазе при больших значениях Da. В турбулентном диффузионном пламени расход топлива определяется исключительно скоростью смешения топлива и окислителя [3].

В EDM-модели также полагается, что реакция происходит только в том случае, когда турбулентные структуры разрушаются до самых малых размеров [4, 5]. Исходя из этих предположений, авторы [1] представили обобщенное выражение для скорости реакции горения топлива:

(1)
$f + SO = \,\,~\left( {{\text{1}} + S} \right){\kern 1pt} p{\text{;}}$
(2)
$R = \rho \frac{\varepsilon }{k}A\,{\text{min}}\left( {{{Y}_{f}}{\text{,}}\,\,~\frac{{{{Y}_{O}}}}{S}{\text{,}}\,\,B\frac{{{{Y}_{р}}}}{{{\text{1}} + S}}} \right){\text{,}}$
где f ‒ топливо; O ‒ окислитель; p ‒ продукты сгорания; S ‒ стехиометрический коэффициент; R ‒ скорость реакции горения; A = 0.4, B = 0.5 ‒ константы; ρ ‒ плотность; ε ‒ скорость турбулентной диссипации; k ‒ кинетическая энергия турбулентности; Yp, Yo, Yf ‒ массовые доли продуктов, окислителя, топлива; индексы обозначают следующее: f ‒ топливо; O ‒ окислитель; p ‒ продукты сгорания.

В настоящей работе используется двухэтапный обобщенный метановоздушный механизм Вестбрука ‒ Драйера [6]

(3)
${\text{C}}{{{\text{H}}}_{{\text{4}}}}\,\,{\text{ + }}\,\,{\text{1}}{\text{.5}}{{{\text{O}}}_{{\text{2}}}} \to {\text{CO}}\,\,{\text{ + }}\,\,{\text{2}}{{{\text{H}}}_{{\text{2}}}}{\text{O;}}$
(4)
${\text{CO}}\,\,{\text{ + }}\,\,{\text{0}}{\text{.5}}{{{\text{O}}}_{{\text{2}}}} \to {\text{C}}{{{\text{O}}}_{{\text{2}}}}{\text{.}}$

МОДЕЛЬ КОЛЛЕКТОРА ГАЗОВЫХ ГОРЕЛОК (FGM-MODEL)

В методе FGM рассматривается многомерное пламя как набор нескольких одномерных единиц пламени [7]. Таким образом, структура одномерного пламени также репрезентативна для структуры многомерного пламени. В ламинарных моделях пламени предполагается, что непрерывная зона реакции состоит из небольших ламинарных участков, встроенных в турбулентный факел с помощью стохастических методов [8‒10].

В областях пламени, где Da >> 1, химические процессы протекают намного быстрее, чем процессы при турбулентном перемешивании, поэтому здесь предполагается быстрое протекание реакций. Химическое равновесие достигается в каждой точке поля потока, и требуется лишь знание локального состава топливовоздушной смеси или ее концентрации. FGM-метод позволяет определять значения всех параметров с помощью нескольких переменных, таких как концентрация смеси, скорость реакции, и решить уравнения переноса для этих параметров с помощью трехмерного CFD-моделирования (computational fluid dynamics). С помощью FGM-метода можно моделировать погасание пламени путем снижения концентрации смеси газов. Уравнения для этой модели представляют собой комбинацию моделей с предварительным смешиванием газов и без него [11, 12]. По модели без предварительного смешивания можно найти концентрацию смеси, рассчитать концентрации компонентов и функцию плотности вероятности (ФПВ), определяющую температуру, энтальпию, удельную теплоемкость и другие параметры для каждого компонента на основе массовых долей [11]. Функция плотности вероятности предназначена для расчета около 20 компонентов [13].

Горелку топочного устройства можно смоделировать с помощью модели горения предварительно смешанного топлива (premixed flamelet) или модели диффузионного пламени (diffusion flamelet). Поскольку воздух подмешивается преимущественно на выходе из горелки, используется модель диффузионного пламени. В FGM-модели средние скорости кинетических реакций получаются путем интегрирования с помощью общей функции плотности вероятности. Общая ФПВ строится на основе β-ФПВ, форма которой моделируется отклонениями как концентрации смеси, так и переменной течения процесса в соответствии со следующими уравнениями:

(5)
$\frac{{\text{d}}}{{{\text{d}}s}}\left( {\sigma m} \right) = \,\,~ - \rho \sigma K{\text{;}}$
(6)
$\begin{gathered} \frac{{\text{d}}}{{{\text{d}}s}}\left( {\sigma m{{Y}_{i}}} \right) - \,\,~\frac{{\text{d}}}{{{\text{d}}s}}~\left( {\sigma \frac{\lambda }{{{\text{L}}{{{\text{e}}}_{i}}{{c}_{p}}}}\frac{{{\text{d}}{{Y}_{i}}}}{{{\text{d}}s}}} \right) = \,\,~\sigma {{\omega }_{i}} - \rho \sigma K{{Y}_{i}} + {{Q}_{i}}{\text{,}} \\ {\text{\;}}~i~\,\, \in \,\,~\left( {{\text{1,}}{{N}_{s}}} \right){\text{;}} \\ \end{gathered} $
(7)
$\begin{gathered} \frac{{\text{d}}}{{{\text{d}}s}}\left( {\sigma mh} \right) - \,\,~\frac{{\text{d}}}{{{\text{d}}s}}~\left( {\sigma \frac{\lambda }{{{{c}_{p}}}}\frac{{{\text{d}}h}}{{{\text{d}}s}}} \right) = \\ = \,\,\frac{{\text{d}}}{{{\text{d}}s}}\sum\limits_{i = {\text{1}}}^{{{N}_{s}}} {\sigma {{h}_{i}}} \frac{\lambda }{{{{c}_{p}}}}\left( {\frac{{\text{1}}}{{{\text{L}}{{{\text{e}}}_{i}}}} - {\text{1}}} \right)\frac{{{\text{d}}{{Y}_{i}}}}{{{\text{d}}s}} - \sigma \rho KH + {{Q}_{h}}{\text{,}} \\ \end{gathered} $
где многомерные свойства отражены с помощью параметров σ, K, Qi и Qh [14]; σ ‒ поле кривизны; K ‒ скорость растяжения пламени; Qi и Qh ‒ члены возмущения; λ, ср, h ‒ теплопроводность, удельная теплоемкость при постоянном давлении, энтальпия дымовых газов; m ‒ массовая скорость горения; Lei ‒ число Леви; ωi ‒ член химического уравнения; Ns ‒ количество компонентов.

С введением переменной течения процесса горения для частично предварительно подготовленной топливовоздушной смеси можно записать трехмерный процесс в адаптированной для пламени системе координат в стационарных уравнениях сохранения массы, видов и энтальпии, как показано в уравнениях (5)(7). Предполагается, что весь процесс происходит в направлении, перпендикулярном переменной течения процесса. С использованием этой переменной можно создать одномерную FGM-модель путем параметризации решений уравнений (5)‒(7). Решение одномерной FGM-модели заносится в базу данных. При решении задачи используется только уравнение сохранения переменной течения процесса. Реакции не просто связаны со смешением топлива с воздухом, но и подвержены изменениям в самой области пламени. Таким образом, FGM-модель теоретически более точна, но по-прежнему подвержена ошибкам как в оценке отклонения, так и в предполагаемой форме β-ФПВ, что может приводить к неточному определению координат пламени.

МЕТОДИКА ИССЛЕДОВАНИЙ

На рис. 1 показана схема испытательного стенда, который содержит четыре газовые горелки с предварительным частичным перемешиванием и теплообменник, состоящий из четырех труб (I‒IV) (на рисунке показана только одна труба). Вытяжной вентилятор 1 подает в горелки 3, расположенные на входе в трубы теплообменника, смесь газов и выводит из установки отработанные газы. Воздух подается на трубы теплообменника нагнетательным вентилятором 2 и после нагревания выводится из установки. Концентрации CO2 и O2 в отработанных газах измеряли газоанализатором. Температура в различных точках на поверхности труб теплообменника определялась с помощью термопар (рис. 2).

Рис. 1.

Схема экспериментального стенда. 1, 2 ‒ вытяжной и нагнетательный вентилятор; 3 ‒ теплообменник; 4 ‒ горелка

Рис. 2.

Труба теплообменника. 1‒8 ‒ точки измерения температуры стенки

Топливо, используемое в топочном устройстве, представляет собой природный газ, который в основном состоит из метана с концентрацией примерно 94%, а также бутана, пропана, CO2 и азота. Испытания проводили при соотношении топлива и воздуха φ = 0.6 и 0.8.

Для создания конечно-элементной сетки модели использовали программный продукт ANSYS, реализующий CFD-моделирование. Некоторые параметры были установлены в качестве целевых для достижения высокого качества и точности. Параметры конечно-элементной сетки для расчетной области подбирались таким образом, чтобы достигалось значение Yplus (y+)11, близкое к 1.

Для получения оптимального количества ячеек сетки было проведено исследование независимости решения от разбиения расчетной области. Полагалось, что трубы теплообменника изготовлены из стали. Подход, основанный на подвижной системе отсчета (movable reference frame ‒ MRF), был использован для моделирования вентилятора. Для минимизации высоких турбулентных пульсаций потока от вентилятора был применен метод псевдоперехода. Основные уравнения турбулентного течения решались с помощью SSTk‒ω-модели турбулентности с двумя уравнениями и сочетанием схем первого и второго порядка. Для решения задачи “давление ‒ скорость” использовали связанную программу с псевдопереходным процессом. Модель дискретных ординат (discrete ordinate ‒ DO) применяли для анализа эффектов теплопереноса. Метод дискретных ординат ‒ это наиболее полная модель теплопереноса в приложении Fluent [9]. Этот метод учитывает рассеяние, полупрозрачность среды, зеркальность поверхности и передачу тепла в зависимости от длины волны с использованием слоисто-серой фильтрации. Модель DO применима для всех диапазонов оптической толщины, но требует больше времени на вычисления по сравнению с другими моделями. Граничные условия массового расхода на входе в установку использовали как для потока приточного воздуха, так и для топливовоздушной смеси. В CFD-модели были смоделированы два разных условия для различных соотношений компонентов. Предполагали, что излучательная способность должна быть постоянной и не зависеть от температуры. Потери тепла в окружающую среду не учитывали, а стенки корпуса установки считали адиабатическими.

ОЦЕНКА НЕЗАВИСИМОСТИ КОНЕЧНО-ЭЛЕМЕНТНОЙ СЕТКИ

Чтобы установить лучшую точность результатов CFD-моделирования, исследовали сходимость сетки путем создания сеток с крупными (М1), средними (М2), мелкими (М3) и очень мелкими (М4) ячейками. Анализ процесса горения был выполнен при соотношении компонентов φ = 0.6 с применением всех четырех сеток для прогнозирования температур дымовых газов, температуры стенок труб теплообменника и содержания CO2.

Ключевые параметры всех четырех сеток представлены в табл. 1, изображения сеток ‒ на рис. 3. Из данных табл. 1 видно, что точность прогноза параметров горения и время CFD-моделирования существенно зависят от числа элементов сетки, определяющих ее оптимальный тип по точности и минимальному времени вычислений.

Таблица 1.  

Параметры сеток моделирования

Параметр Ячейка сетки
М1 М2 М3 М4
Количество элементов, млн 25 34 42 78
Время вычисления, ч 10 18 25 56
Температура дымовых газов 0.94 А 0.96 А А А
Содержание СО2 0.82 В 0.87 В 0.90 В 0.91 В
Рис. 3.

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

В настоящей работе все данные представлены в виде нормированных величин. Измеренные значения температуры дымовых газов и содержания CO2 обозначены как A и B (см. табл. 1). Значения параметров, полученные с помощью CFD-моделирования с различными типами сеток, сравниваются с A и B. Сетки M3 и M4 дают лучший прогноз температуры дымовых газов и содержания CO2 по сравнению с сетками M1 и M2. Однако время вычисления для сетки M4 слишком велико.

На рис. 4 для сравнения показаны температуры стенки труб теплообменников, рассчитанные CFD-методом и измеренные для всех типов сеток. Сетки M3 и M4 дали лучший прогноз температуры стенки трубы относительно измеренных значений.

Рис. 4.

Относительная температура стенки трубы теплообменника tст, измеренная (1) и рассчитанная CFD-методом (25) при φ = 0.6. Сетка: 2 ‒ М1; 3 ‒ М2; 4 ‒ М3; 5 ‒ М4

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

РЕЗУЛЬТАТЫ ИССЛЕДОВАНИЙ

На рис. 5 показан контур пламени на начальном участке труб теплообменника. Данные рис. 5 свидетельствуют о том, что длина пламени зависит от соотношения содержания компонентов смеси, а температуры пламени, вычисленные с помощью EDM-модели, выше температур, рассчитанных с помощью FGM-модели. На рис. 6 показана схема сечений трубы теплообменника, на рис. 7 ‒ распределение температуры пламени в этих сечениях, на рис. 8 ‒ то же по длине трубы теплообменника.

Рис. 5.

Относительная температура пламени tпл на начальном участке труб теплообменника по FGM-модели (а, в) и EDM-модели (б, г) при φ = 0.6 (а, б) и φ = 0.8 (в, г)

Рис. 6.

Схема сечений трубы теплообменника

Рис. 7.

Относительная температура пламени в сечениях труб теплообменника, рассчитанная по FGM-модели (а, в) и EDM-модели (б, г) при φ = 0.6 (а, б) и φ = 0.8 (в, г)

Рис. 8.

Относительная температура пламени по длине трубы теплообменника l, рассчитанная с помощью FGM-модели (1) и EDM-модели (2) при φ = 0.6 (а) и φ = 0.8 (б)

В табл. 2 представлен КПД теплового процесса, полученный при проведении эксперимента и расчета с использованием EDM- и FGM-моделей. Данные испытаний были приняты за X% и Y%. Как следует из данных табл. 2, FGM-модель дает более точный прогноз КПД (в пределах 1%), тогда как EDM-модель дает разницу около 5%. В основном это обусловлено распределением температуры и изменением содержания CO2 по длине трубы теплообменника.

Таблица 2.  

Коэффициент полезного действия теплового процесса

Источник φ
0.6 0.8
Данные эксперимента Х% Y %
Расчет по модели:    
FGM 0.99 Х% 0.99 Y%
EDM 0.95 Х% 0.95 Y%

В табл. 3 приведено содержание CO2 на выходе из дымовой трубы, которое повышается с увеличением соотношения компонентов, но различие содержания минимально для обоих значений φ. Содержание CO2, рассчитанное по EDM-модели, немного больше по сравнению с результатами FGM-моделирования в начальной зоне пламени (рис. 9). Более детальное сравнение содержания СО2 в дымовых газах по длине трубы теплообменника показано на рис. 10.

Таблица 3.  

Содержание СО2, молярные доли, в дымовых газах

Источник φ
0.6 0.8
Данные эксперимента 0.059 0.076
Расчет по модели:    
FGM 0.053 0.068
EDM 0.055 0.072
Рис. 9.

Содержание СО2 в дымовых газах по FGM-модели (а, в) и EDM-модели (б, г) при φ = 0.6 (а, б) и φ = 0.8 (в, г)

Рис. 10.

Содержание СО2 в дымовых газах по длине трубы теплообменника, рассчитанное по FGM-модели (1) и ЕDM-модели (2) при φ = 0.6 (а) и φ = 0.8 (б)

На рис. 11 приведены температуры стенок труб теплообменника. Результаты испытаний и CFD-моделирования хорошо согласуются для вариантов с различными значениями φ, хотя наблюдаются и некоторые различия. Более высокие температуры на стенках трубы наблюдаются в EDM-модели, особенно после первого изгиба трубы.

Рис. 11.

Относительная температура стенки труб теплообменника, измеренная (1) и рассчитанная по EDM-модели (2) и FGM-модели (3) при φ = 0.6 (а) и φ = 0.8 (б)

На рис. 11 точки 1‒4 на поверхности трубы теплообменника представляют собой область, где начинается распространение пламени. EDM-модель хуже прогнозирует температуру в этой области, так как развитие пламени происходит после определенного расстояния от точки воспламенения. Точки 5‒7 представляют собой среднюю часть трубы теплообменника, в которой пламя развивается максимально. На выходе из теплообменника EDM-модель прогнозирует более высокие температуры дымовых газов, нежели измеренные, для обоих значений φ. Точка 8 находится рядом с концом трубы, где от нее отводится большая часть тепла. Видно, что FGM-модель лучше прогнозирует заниженные значения температуры в силу более высокого значения коэффициента теплопередачи, заложенного в эту модель Общие результаты расчетов показывают, что FGM-модель обеспечивает более точный прогноз по сравнению с EDM-моделью.

ВЫВОДЫ

1. Пространственное распределение фронта горения и взаимодействия потоков для определения формы и местоположения пламени можно получить наиболее точно с помощью FGM‑модели.

2. Основным недостатком использования FGM-моделирования по сравнению с EDM-моделированием является более длительное время вычислений.

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

Список литературы

  1. Magnussen B.F., Hjertager B.H. On the mathematical modeling of turbulent combustion with special emphasis on shoot formation and combustion // Symp. (Int.) Combust. 1977. V. 16. Is. 1. P. 719–729. https://doi.org/10.1016/S0082-0784(77)80366-4

  2. Tunç G., Yılmaz I. Comparison of the numerical models for the temperature distributions of non-premixed swirling methane flame // J. Polytechn. 2019. № 22 (4). P. 819‒826. https://doi.org/10.2339/politeknik.448529

  3. Spalding D.B. Mixing and chemical reaction in steady confined turbulent flames // Symp. (Int.) Combust. 1971. V. 13. Is.1. P. 649–657. https://doi.org/10.1016/S0082-0784(71)80067-X

  4. Stefan M., Sebastian K., Ognijan B. Development of an eddy dissipation model for the use in numerical hybrid rocket engine combustion simulation // 7th European Conf. on Aeronautics and Space Sciences (EUCASS). Milan, Italy, 02‒06 July 2017. https://doi.org/10.13009/EUCASS2017-157

  5. Lewandowski M.T., Ertesvag I.S. Analysis of eddy dissipation concept formulation for MILD combustion modeling // Fuel. 2018. V. 224. P. 687‒700. https://doi.org/10.1016/j.fuel.2018.03.110

  6. Westbrook C., Dryer F. Chemical kinetic modelling of hydrocarbon combustion // Prog. Energy Combust. Sci. 1984. V. 10. Is. 1. P. 1‒57. https://doi.org/10.1016/0360-1285(84)90118-7

  7. Van Oijen J.A., de Goey L.P.H. Modelling of premixed laminar flames using flamelet-generated manifolds // Combust. Sci. Technol. 2000. V. 161. Is. 1. P. 113‒137. https://doi.org/10.1080/00102200008935814

  8. Stat-of-the-art in premixed combustion modelling using flamelet generated manifolds / J.A. Van Oijen, A. Donini, R.J.M. Bastiaans, J.H.M. ten Thije Boonkkamp, L.P.H. de Goey // Prog. Energy Combust. Sci. 2016. V. 57. P. 30‒74. https://doi.org/10.1016/j.pecs.2016.07.001

  9. Van Oijen J.A., de Goey L.P.Y. Modelling of premixed counter flow flames using the flamelet-generated manifold method // Combust. Theory Modell. 2002. V. 6. Is. 3. P. 463‒478. https://doi.org/10.1088/1364-7830/6/3/305

  10. Van Oijen J.A., Lammers F.A., de Goey L.P.H. Modeling of complex premixed burner systems by using flamlet generated manifolds // Combust. Flame. 2001. V. 127. Is. 3. P. 2124‒2134. https://doi.org/10.1016/S0010-2180(01)00316-9

  11. The flamelet generated manifold method applied to steady planar partially premixed counter flow flames / H. Bongers, J.A. Van Oijen, L.P.H. de Goey, L.M.T. Somers // Combust. Sci. Technol. 2005. V. 177. Is. 12. P. 2373‒2393. https://doi.org/10.1080/00102200500241198

  12. Multidimensional flamelet-generated manifolds for partially premixed combustion / Nguyen Phuc-Danh, L. Vervisch, V. Subramanian, P. Domingo // Combust. Flame. 2010. V. 157. Is. 1. P. 43‒61. https://doi.org/10.1016/j.combustflame.2009.07.008

  13. ANSYS Fluent theory Guide V16.2.

  14. De Goey L.P.H., ten Thije Boonkkamp J.H.M. A flamelet description of premixed laminar flames and relation with flame stretch // Combust. Flame. 1999. V. 119. Is. 3. P. 253‒271. https://doi.org/10.1016/S0010-2180(99)00052-8

Дополнительные материалы отсутствуют.