УДК 662.767/94
Математическая модель вепольного горения жидкого
топлива
Громцев С.А.
Санкт-Петербургский Государственный университет низкотемпературных и пищевых технологий
В диапазоне напряженности поля от 20 до 40 кВ/м, что наиболее приемлемо для тепловых аппаратов малой теплопроизводительности, повышение массовой скорости горения составляет от 25 до 50 % и более, что позволяет в зависимости от необходимости или значительно уменьшить размеры топочных камер, или увеличить расход топлива и теплопроизводительность установки с сохранением ее размеров. Оба варианта использования вепольного эффекта имеют большое значение для стационарных и мобильных тепловых установок.
Ключевые слова: математическая модель, жидкое топливо, горение.
Как показал анализ состояния исследований процесса горения жидкого топлива, наиболее приемлемой математической моделью процесса является следующее выражение:
m = (41 / Cp )х 1п{1 +1/L[Cp (Т - Тв) + q0Х01 / W0v0 ]}, (1)
так как оно полностью согласуется с постулативными утверждениями, характеризующими количество вещества, подаваемого в зону реакции для поддержания горения, и скорость его убывания.
Одновременно выражение (1) показывает зависимость массовой скорости горения ш от свойств топлива ш = /(1Ср), теплоты реакции его окисления д0, теплоты испарения Ь и слабую зависимость от температуры горения ш = /(Т- Тв) при независимости от давления.
Анализируя (1) представим зависимость массовой скорости горения ш от ряда переменных величин для конкретного вида топлива:
ш = /МГ,0. (2)
С учетом утверждения [1], что скорость горения топлива является функцией скорости химических реакций, т.е. ш = /(ин), где ин = /(у, а), а скорость химических реакций — через скорость образования активных центров у0 по зависимости (2) и с учетом теоретических предпосылок — зависит от параметров электрического поля, получим:
ш = /(Е). (3)
С учетом изложенного можно записать, что:
т = /(г, АТ, д, Е), (4)
а в общем случае:
т = /(г, АТ, д, Е, Ь, Ср, 1), (5)
при коэффициенте расхода воздуха а = 1.
Выражение (5) отражает влияние на массовую скорость горения жидкого топлива в электрическом поле всех воздействующих факторов и может быть представлено как ее математическая модель.
Из анализа выражения (5), с учетом работ [2, 3], видно, что зависимость т = /(Е) определяется зависимостями радиуса г (а через него общей поверхности капель), температуры горения и теплоты химической реакции горения от напряженности электрического поля, т.е.:
т = Е[(г=/7(Е), АТ=ЯЕ), д=/з(Е), Ь, Ср , 1]. (6)
Из вышеизложенного следует, что математическую модель вепольного горения жидкого топлива в атмосфере окислителя более корректно представить системой уравнений:
т = Е (г, АТ ,во, Ь, Ср ,1)
г = /1( Е), АТ = /2 (Е), во = /з(Е),
С учетом (1), систему (7) можно записать в следующем виде: т = (41-1 /Ср )х 1п{1 +1/ь[Ср (Т - Тъ) + gо¥о /}
г1 = С1 + С2 [С3 + С4 ехр(^1 Е + В1)] ,
АТ = С5 + С6 [С7 + С8 ехр(А2Е + В2 )]-1,
g0 = С9 + С10 [С11 + С12 еХР (А3 Е + В3 )]-1.
Эмпирические коэффициенты А]-, В) и С]- определялись экспериментально с использованием при обработке опытных данных, теории подобия, теоремы Гаусса и вычислительной техники на основе моделирования процессов горения в тепловых установках характерных пищевым производствам.
Так как на входе в канальный газоход неизвестной является не только температура потока газов, но и его скорость, то определение реальной температуры осуществлялось согласно решению нелинейного уравнения:
Тт
guх - g- g Р.с - и2 /2 - | СрёТ = о, (8)
(7)
где gи.с.х — тепло, вносимое в топку; gп.т.p — потери тепла в камере сгорания; gp.с — энергия разрыва связей топлива; и/2 — кинетическая энергия потока;
т
| СрёТ — энтальпия продуктов сгорания.
о
о
Уравнение (8) решалось на ЭВМ при условии, что состав топлива и воздуха, их расход и температура на протяжении одного эксперимента были неизменны. Исходя из формулы (8) и на основании выдвинутых ранее теоретических предпосылок можно утверждать, что прирост температуры сгорания топлива при наложении электрического поля произошел за счет уменьшения энергии разрыва межмолекулярных связей топлива под воздействием энергии наложенного электрического поля. Из рисунка 1 следует, что при увеличении напряженности электрического поля, наложенного на факел пламени, температура горения возрастает и для мазута составляет 180—190°С при увеличении Е до 50 кВ/м. Путем аппроксимации данных эксперимента были получены эмпирические зависимости прироста температуры горения от напряженности наложенного электрического поля: для горения керосина:
Мк = 1,6 х102 -1,5 х102 /{1 + ехр[0,22(Е - 29,2)]} (9)
для горения дизельного топлива:
ЛТд = 1,8х 102 -1,7х 102/{1 + ехр[0,22(Е - 27,5)]} (10)
для горения мазута:
ЛТм = 1,95х 102 -1,91 х 102/{1 + ехр[0,30(Е - 25,1)]} (11)
Удельная теплота химических реакций или теплопроизводительность топлива определялась калориметрированием топок по стандартным методикам. Результаты калориметрирования показали, что при горении керосина прирост теплопроизводительности составляет 7—8, дизтоплива 8—9, а мазута 10—12 кДж/кг (рис.2).
При аппроксимации данных калориметрирования получены зависимости прироста удельной теплоты химических реакций от напряженности электрического поля.
АТ, оС
Дgo, кДж/кг
" Г7 ' I 1 I ' * • I ■ I I I I \ . I .
0.0 5.0 10.0 15.0 20.0 25.0 30.0 35.0 40.0 45.0 50.0
Е, кВ/м
Рис. 1. Прирост температуры горения мазута.
Е, кВ/м
Рис. 2. Прирост удельной теплоты сгорания мазута.
В случае сжигания керосина зависимость имеет вид: Аg^ = 7,7 - 7,5 /{1 + ехр[о,2о(Е - 24)]}; в случае сжигания дизельного топлива:
АgoД = 8,9 - 8,8/{1 + ехр [о,2о(Е - 23)]};
(12)
(13)
(14)
в случае сжигания мазута:
АgoM = 1о,2 - 1о,3 /{1 + ехр [о,21(Е - 2о)]}.
Изменение радиуса капли топлива замеряли в соответствие с [3], в результате апроксимации получили следующие уравнения: для керосина:
-8 , п £ ^ 1 А-5
гк = 5,4 х 1о-8 + 9,6 х 1 о-5 ехр(- о,о3 8Е); для дизельного топлива:
гд = 1,7 х 1 о-5 +1,8 х 1 о-4 ехр- о,о28Е);
для мазута:
гм = -1,6 х 1 о-5 + 5,2 х 1 о-4 ехр(- о,о 15Е). Результаты расчетов представлены в виде графика на рис. 3.
(15)
(16) (17)
гх10,м
Ашх10"6, кг/с
0.0 5.0 10.0 15.0 20.0 25.0 30.0 35.0
Е, кВ/м
Е, кВ/м
Рис. 3. Диспергирование капли мазута. Рис. 4. Прирост массовой скорости
горения: ▲ — мазут, ■ — диз. топливо, • — керосин.
Из анализа состояния исследований процесса горения жидкого топлива, следует, что одной из важных характеристик этого процесса является массовая скорость горения топлива. В отношении форсуночных устройств данная характеристика в значительной степени соответствует понятию «расход топлива» и имеет размерность [кг/с]. При контрольном разогреве теплового устройства с помощью запорной арматуры и приборов контроля за расходом топлива устанавливался стандартный (паспортный) расход топлива, при этом фиксировались размеры факела (главным образом длина) и сопоставлялись с размерами топки. Для этого топка по длине имела несколько смотровых люков.
При работе теплового устройства по вепольной технологии, когда размеры факела изменились, их доводили до стандартных значений путем увеличения или уменьшения величины подачи топлива. Показания расходомеров фиксировались в журнале наблюдений. Таким образом, осуществлялся замер расхода жидкого топлива, в условиях эксперимента соответствующий массовой скорости горения.
В экспериментах производился замер массовой скорости горения жидких топлив, а именно: керосина и дизельного топлива на ПХП-4, мазута на ФТЛ-2, в диапазоне напряженности поля от 5 до 50 кВ/м. Результаты экспериментов показывают, что в указанном диапазоне массовая скорость горения керосина возрастает на 9 кг/ч, дизельного топлива — на 10 кг/ч и мазута — на 13 кг/ч, что в пересчете составит 8х10-6, 24х10-6 и 76х10-6 кг/ч соответственно, с учетом диспергирования топлива (рис. 4).
Результаты экспериментов свидетельствуют о том, что при вепольном сжигании жидкого топлива с повышением напряженности электрического поля массовая скорость горения факела и капель топлива повышается, причем интенсивность указанного процесса возрастает с увеличением вязкости топлива и понижением других качественных характеристик.
Графики подтверждают тот факт, что, во-первых, оптимальным является диапазон напряженности электрического поля от 2о до 4о кВ/м, а во-вторых, полученный эффект обратно пропорционален качеству топлива.
Математическая модель горения мазута в атмосфере окислителя при нормальном атмосферном давлении после определения коэффициентов апроксимации будет выглядеть следующим образом (18) (рис.2):
т = (41 / Ср )х 1п{1 +1/ Ь[Ср (Т - Ть) + gоУо / Ж^ }
гм =-1,6 х1о-0,5 + 5,2 х1о -о,4ехр (- о,о15 Е),
АТм = 1,9х 1о2 -1,9х 1о2/{1 + ехр[о,3о(Е-25,1)]}, ' (18)
g0м = 1о,2х1оо - 1о,3х1оо/{1 + ехр[о,21(Е -2о,о)]}.
Зона отрицательного эффекта при напряженности от о до 2о кВ/м (рис. 4) объясняется затратой энергии на превращение поступательного движения ионов в колебательное с частотой поля и затратами энергии на переориентацию дипольных молекул в соответствие с силовыми линиями поля. Массовая скорость вепольного горения жидких топлив в зависимости от напряженности электрического поля при уменьшении радиуса капли определялась по формуле (1), причем Ти д0 определялись экспериментально.
В диапазоне напряженности поля от 2о до 4о кВ/м, что наиболее приемлемо для тепловых аппаратов малой теплопроизводительности, повышение массовой скорости горения составляет от 25 до 5о% и более, что позволяет в зависимости от необходимости или значительно уменьшить размеры топочных камер, или увеличить расход топлива и теплопроизводительность установки с сохранением ее размеров. Оба варианта использования вепольного эффекта имеют большое значение для стационарных и мобильных тепловых установок.
Список литературы
1. Вулис Л. А. Перекисная теория горения./М.: Госэнергоиздат, 1954 г.
2. Вильямс Ф.А. Теория горения./ М.: Наука, 1971 г.
3. Громцев С. А. Вепольные технологии повышения эффективности энергетических источников техники тыла./Л.: ВАТТ, 1995 г.