Физико-химические основы горения топлива
1. Исходные данные для проектирования.
Курсовой проект выполняется на основании задания на проектирования.
Задание №13
Состав природного газа:
СН4-84,41%
С2Н6-6,72%
С3Н8-1.77,%
С4Н10-0,43%
С5Н12-0,27%
Н2S-0,36%
N2-2,59%
СО2-0,28%
Влагосодержание топлива dT=0,01 мг/м3
Коэффициент избытка воздуха αB=1,1
Влагосодержание воздуха dB=0,015 мг/м3
Температура смеси То=450 К
Давление смеси 101,3 кПа
Высота топки 4000 мм
Диаметр топки 3000 мм
2. Определение состава продуктов полного сгорания топлива.
Теоретическое количество кислорода, необходимого для сгорания 1м3 топлива, определяется по формуле:
,
где: H2S, - процентное содержание сероводорода в 1м3 топлива;
CnHm- процентное соотношение i-го углеводорода в 1м3 топлива;
n, m – число атомов соответственно углерода и водорода в химической формуле 1-го углеводорода.
Так как в атмосфере воздуха содержится 21% кислорода, количество воздуха, теоретически необходимого для полного сгорания 1 м3 топлива, будет в 1/0,21=4,76 раза больше чем кислорода:
Количество двуокиси углерода, образующейся при полном сгорании 1 м3 газового топлива, определяется по следующей зависимости:
где: СО2Т, - процентное содержание двуокиси углерода в балансе топлива;
Количество сернистого газа, образующемся при сгорании 1 м3 топлива:
Количество водяных паров, образующихся при сгорании 1м3 топлива:
Количество кислорода, входящих в состав продуктов сгорания при αB>1 .
Количество азота, образующийся при сгорании 1м3 топлива:
где: N2Т - процентное содержание азота в балансе топлива;
Полный объём продуктов сгорания 1м3 топлива:
где: Vсух- объём сухих продуктов сгорания.
Результаты расчёта сводятся в таблицу 1
Таблица 1
aв |
Выход продуктов сгорания, м3 |
||||||
|
H2О |
СО2 |
SO2 |
N2 |
О2 |
Vсух |
VГ | |
|
aв=1 |
2,20 |
1,06 |
0,004 |
7,77 |
0 |
8,84 |
11,04 |
aв=aв зад |
2,2 |
1,06 |
0,0036 |
8,85 |
0,214 |
9.83 |
12,03 |
Для проверки результатов вычислений выполняется расчёт состав продуктов полного сгорания природных газов. Результаты расчёта сводятся в таблицу 2
Таблица 2
Компоненты газо-воздушной смеси и уравнения реакций горения |
Расход воздуха, м3 |
Выход продуктов горения, м3 | |||||||||
|
O2 |
N2 |
CO2 |
H2O |
N2 |
O2 |
SO2 | |||||
|
СН4 |
0,8441 |
1,688 |
6,348 |
0,844 |
1,688 |
6,348 |
- |
- | |||
С2H6 |
0,0672 |
0,235 |
0,884 |
0,134 |
0,202 |
0,884 |
- |
- | |||
С3H8 |
0,0177 |
0,089 |
0,333 |
0,053 |
0,071 |
0,333 |
- |
- | |||
С4H10 |
0,0043 |
0,028 |
0,105 |
0,017 |
0,022 |
0,105 |
- |
- | |||
С5Н12 |
0,0027 |
0,022 |
0,081 |
0,014 |
0,016 |
0,081 |
- |
- | |||
H2S |
0,0036 |
0,005 |
0,020 |
- |
0,004 |
0,020 |
- |
0,004 |
|||
Итого для горючих компонентов топлива при αB=1,0 |
1 |
2,07 |
1,06 |
2,00 |
7,77 |
0,004 |
|||||
9,84 | |||||||||||
СО2 N2 Водяной пар в 1м3 топлива (dт, кг/м3) |
0,0028 |
- |
- |
- |
- | ||||||
- |
- |
0,0259 |
- |
- | |||||||
- |
- |
0,0124 |
- |
- |
- | ||||||
Избыточный воздух с содержанием пара dв, кг/м3 |
0.984 |
- |
0,201 |
0,777 |
0,21 |
- |
- |
- | |||
Итого для газо-воздушной смеси |
10,822 |
1,0628 |
2,0124 |
8,5745 |
0,2066 |
0,004 |
|||||
3. Определение
адиабатной температуры
3.1 Определение
адиабатной температуры
при постоянном объёме.
,
где: ri-доля компонента
Qнi-теплота сгорания i-го компонента.
ТТ=ТВ=То=450 К=176,85 0С
,
где: Vg. CVg- количество и средняя в диапазоне 0…Та изохорная объёмная теплоёмкость g-го компонента продуктов сгорания.
Тапр=800 0С
Тапр=1000 0С
Тапр=1200 0С
Тапр=1400 0С
Тапр=1600 0С
Тапр=1800 0С
Тапр=2000 0С
Тапр=2200 0С
Тапр=2400 0С
Тапр=2600 0С
По данным значениям строится график 1, по которому вычисляется истинное значение Та.
Рис 1. График зависимости адиабатной температуры горения газовой смеси при постоянном объёме.
По графику Та = 2400 °С
3.2 Определение адиабатной температуры горения газовой смеси,
при постоянном давлении.
,
где: Vg. CРg- количество и средняя в диапазоне 0…Та изобарная объёмная теплоёмкость g-го компонента продуктов сгорания.
Тапр=800 0С
Тапр=1000 0С
Тапр=1200 0С
Тапр=1400 0С
Тапр=1600 0С
Тапр=1800 0С
Тапр=2000 0С
Тапр=2200 0С
По данным значениям строится график 2, по которому вычисляется истинное значение Та.
Рис 1. График зависимости адиабатной температуры горения газовой смеси при постоянном давлений.
По графику Та= 2000 °С
4. Определение основных параметров адиабатного воспламенения
Анализ параметров
адиабатного воспламенения
Безразмерная температура:
Безразмерная концентрация:
Критерий Аррениуса (мера реакционной способности смеси):
Чем больше Аrr, тем труднее протекает реакция горения смеси.
Безразмерное время:
4.1 Определение
основных параметров
горения газовой смеси, при постоянном объеме.
Безразмерная температура для изохорного процесса, определяется по формуле:
,
где: То - начальная температура газовой смеси Т0=450К
ТА - адиабатная температура горения газовой смеси, при постоянном объеме.
ТА=2400+273=2673 К
Критерий Аррениуса, определяется по формуле:
где: Е-энергия активации, принимается Е=124024, КДж/моль
R- универсальная газовая постоянная, КДж/мольК
Определяется безразмерное время индукции адиабатного теплового воспламенения.
где: - температура воспламенения в момент
Определяется теплота в точке максимума тепловыделения:
Безразмерные величины переводятся в размерные:
Определяем период горения (время резкого роста температуры периода индукции):
а) безразмерный период:
Dxг = (1-q0)/(1-q*)*ехр(Arr/q*) = (1- 0,168)/(1-0,865)*ехр(5,57/0,
б) реальный период:
Dtг = Dxг/k0 = 3727,59/5,58*1032 = 7,44*10-30, с;
4.2 Определение
основных параметров
горения газовой смеси, при постоянном давлении.
Безразмерная температура для изобарного процесса, определяется по формуле:
,
где: То - начальная температура газовой смеси Т0=450К
ТА - адиабатная температура горения газовой смеси, при постоянном давлении
ТА=2000+273=2273 К
Критерий Аррениуса, определяется по формуле:
где: Е-энергия активации, принимается Е=124024, КДж/моль
R- универсальная газовая постоянная, КДж/мольК
Определяется безразмерное время индукции адиабатного теплового воспламенения.
где: - температура воспламенения в момент
Определяется теплота в точке максимума тепловыделения:
Безразмерные величины переводятся в размерные:
Определяем период горения (время резкого роста температуры периода индукции):
а) безразмерный период:
Dxг = (1-q0)/(1-q*)∙ехр(Arr/q*) = (1- 0,197)/(1-0,8815)∙ехр(6,56/0,
б) реальный период:
Dtг = Dxг/k0 = 11402,45/6,91*1032 = 1,901*10-29, сек.
Предэкспонента К0, определяется по формуле:
где: R- универсальная газовая постоянная, КДж/мольК
Т0-начальная температура смеси, К
Е- энергия активации реакций метана в смеси с воздухом перед воспламенением, Е=124024, КДж/моль,
ТА-аддиабатная температура горения газовой смеси, при постоянном объёме, или давлении,0К
Определяется значение предэкспонента К0 для изохорного процесса.
Определяется значение предэкспонента К0 для изобарного процесса.
5. Расчёт
коэффициентов
многокомпонентной газовой смеси.
5.1 Определение коэффициентов теплопроводности.
Зависимость коэффициента теплопроводности от температуры смеси выражается формулой Сазерленда, однако при достижении достаточно высоких температур, условно можно принять:
l=lо∙(То/TА)0,5
Для определения lо используется формула Масона-Саксена:
где: lm-коэффициент теплопроводности
m-го компонента смеси;
rm, rj - объемные доли m-го и j-го компонентов смеси;
где : mrj, mrm - молекулярные массы j-го и m-го компонентов смеси;
hj, hm - коэффициенты динамической вязкости j-го и m-го компонентов смеси.
Значения m, h и l приводятся в справочной литературе.
Расчёт lо ведётся по программе Microsoft Exel
lо=0,026 Вт/(м2К)
Так как коэффициент теплопроводности зависит от температуры, то l для процессов с постоянными объемом и давлением определяется по формуле:
l= lо*(Tо/ТА)0,5,
При постоянном давлении Р=const:
При постоянном объёме V=const:
5.2 Определение
фиктивного коэффициента
Фиктивный коэффициент теплоотдачи a определяется по формуле:
a=4∙l/R0,
где: l - коэффициент теплопроводности газовой смеси
R0 - радиус сосуда в котором находится данная смесь
при постоянном давлении
при постоянном объёме
6. Расчёт
пределов распространения
учётом химической кинетики и реакции горения.
Верхним и нижним концентрационным пределом распространения пламени называется максимальное и минимальное содержание горючего в его смеси с окислителем, при котором возможно распространение пламени в смеси на любое расстояние от источника зажигания.
Минимальное флегматизирующее содержание индивидуальных инертных компонентов, определяется по формуле:
где: стандартная теплота образования горючей части газа, при 25 0С, кДж/моль
тс, тн, то-число атомов углерода, водорода и кислорода в условной формуле горючей части газа.
Количество С в горючей смеси: 1+2+3+4+5=15
Количество Н в горючей смеси: 4+5+8+10+12+2=42
Число атомов углерода в условной формуле горючей части газа, определяется по формуле:
где: Спр,, спр, Снепр, снепр Нпр,, нпр, Ннепр, ннепр -число атомов углерода, водорода и количество предельных и непредельных углеводородов.
Сгор- сумма горючих компонентов
Стандартная теплота
образования горючей части
где: -теплота образования и количество предельных и непредельных углеводородов.
Минимальное флегматизирующее содержание всех инертных компонентов, находящихся в газе, определяется по формуле:
где: минимальное флегматизирующее содержание индивидуальных инертных компонентов, %
Кф-коэффициент флегматизации
Коэффициент флегматизации определяется по формуле:
где: , - коэффициенты флегматизации;
, – минимальное флегматизирующее содержание индивидуальных инертных компонентов, %
Коэффициенты
флегматизации минимальное
Значения коэффициентов Аи В находятся по формулам:
где: Н2, СО, Спр, Снепр, N2, СО2, Н2О и О2 – объёмные доли газов, предельных и непредельных углеводородов в смеси, %; Сгор – сумма горючих компонентов, %;
, , , , , , , – соответственно, нижний и верхний пределы распространения пламени газов;
Кф и fф - коэффициент флегматизации и минимальное флегматизирующее содержание всех инертных компонентов, находящихся в газе; Δ - поправка, учитывающая содержание кислорода в газе, %.
Δ – поправка, учитывающая содержание кислорода в газе, %.
Δ=0, т.к О2=0
Пределы распространения пламени определяются следующими соотношениями:
т.к. температура газовоздушной среды отлична от 250С, то нижний и верхний пределы распространения пламени рассчитывают по формулам:
7. Расчёт
процессов распространения
Интенсивность процесса распространения пламени в неподвижных или ламинарно движущихся горючих смесях характеризуется нормальной скоростью распространения пламени и массовой скоростью горения.
1. Нормальная
скорость распространения
Нормальная скорость распространения пламени, определяется по формуле:
,
где: А-безразмерный коэффициент, который зависит от начальной температуры смеси, температуры горения и кинетики реакции горения;
асм - температуропроводность смеси
-время химической реакции во фронте пламени.
Расчёт нормальной скорости по данному соотношению затруднён, поэтому значения ир определяют эксперементально.
2. Массовая скорость горения им- количество горючей смеси сгорающую в единицу времени на единице поверхности фронта пламени.
Массовая скорость горения им , определяется по формуле:
,
где: плотность горючей смеси.
Плотность смеси определяется по формуле:
n
ρсм=0,01∑ Viρi ,
i=1
где: Vi – объем газа, %
n- количество компонентов
-плотность отдельных
ρсм=0,01∙ (84,41∙0,717+6,72∙1,356+1,77∙
+0,36∙1,251)=0,75 кг/м3
Для газообразных горючих смесей, состоящих из родственных компонентов, нормальная скорость подчиняется закону аддитивности. Её максимальное значение для конкретной горючей смеси определяется формулой:
где: -максимальные скорости для нормального горения для газовоздушных смесей каждого из компонентов горючего.
Zi-объёмные процентные содержания в сложной газовоздушной смеси простых (двойных) смесей с составом, соответствующим для каждой смеси, %
где: ri-объёмное содержание горючих газов в сложном несбалансированном газе;
объёмные процентные содержания горючих компонентов в простых газовоздушных смесях, соответствующие ,%
объёмное процентное содержание сложного газа в газовоздушной смеси, %
Объёмное процентное содержание сложного газа в газовоздушной смеси, определяется по формуле:
Если в газе содержится балласта более 5%, то вводится поправка:
где: N2, СО2-процентные содержания азота и диаксида углерода в горючем газе, %
Балласт: N2+CО2=2,59+0,28=2,87%<5%, поправка не вводится.
8. Список используемой литературы.
1. Методическое пособие “Последовательность расчётного определения параметров самовоспламенения природного газа” В. Н. Диденко, Ижевск, 1996 г.
2. Методическое пособие “Воспламенение природных газов (основное положения)” В. Н. Диденко, Ижевск, 1996 г.
3. Методическое пособие ”Определение состава продуктов полного сгорания топлива” В. Н. Диденко.
4. Ионин В. А. “Газоснабжение”, М., Стройиздат, 1989 г.
5. Методическое пособие “Расчет состава и температуры продуктов сгорания природного газа с учетом диссоциаций”. Ижевск ИжГТУ 1998 г.