close

Вход

Забыли?

вход по аккаунту

?

bd000103342

код для вставкиСкачать
На правах рукописи
Cbejuxju^^o
Демидов Максим Игоревич
ОБРАЗОВАНИЕ ОКСИДОВ УГЛЕРОДА И АЗОТА П Р И Г О Р Е Н И И И
ДОГОРАНИИ В Д В И Г А Т Е Л Е В Н У Т Р Е Н Н Е Г О С Г О Р А Н И Я
Специальность 05.04.02 - Тепловые двигатели
АВТОРЕФЕРАТ
диссертации на соискание ученой степени
кандидата технических наук
ТУЛА 2005
\
Работа выполнена на кафедре «Автомобили и авгомобильнос хозяйство»
ГОУ ВПО Тульский государственного университета
Научный руководитель
кандидат технических наук, доцент
Чесноков Сергей Александрович
Официальные оппоненты:
доктор технических наук, профессор
Поляков Евгений Павлович
кандидат технических наук, доцент
Котюхов Федор Алексеевич
Ведущее предприятие: f^l^filC , ffjfimfi^^f^^'^^' ^^' ^■^'
^^^^^'^
Защита состоится « l b » p l g K f t g p J t 2005 года в * £ часов на заседании
диссертационного совета КР 212.271.44~ при ГОУ ВПО Тульский государственный
университет по адресу: 300600, г.Тула, пр-т Ленина 92 (3 уч. корпус ауд.316)
С диссертацией можно ознакомиться в библиотеке Тульского государственного
университета
Автореферат разослан « /^» UoJtolpX
2005 г.
Ученый секретарь
диссертационного совета
л f ^
f Lf?-^'¥~—
И.Ь.лгуреев
\
p*/^s^
3
ЛЛ^^^
ОБЫЧАЯ ХАРАКТЕРИСТИКА РАБОТЫ
Актуальность. Двигатели внутреннего сгорания (ДВС) ифают существенную
роль в загрязнении окружающей среды. В крупных городах они являются одним из
главных источников токсичных веществ, выбрасываемых в атмосферу.
Объемное содержание токсичных веществ в отработавших газах (ОГ) сравни­
тельно невелико и составляет 0,2 - 2 % . При этом около 90 % от общей массы ток­
сичных компонентов ОГ приходится на долю пяти основных веществ: СО, NO,, СН,,
альдегидов RCHO, диоксида серы SO2. Нормируемыми токсичными компонентами
ОГ в соответствии с современными нормативными документами (EUR0-4) являются
СО, N0,, СИх и твердые частицы (для дизелей). Нормирование выбросов оксидов се­
ры осуществляется косвенно, через ограничение содержания серы в топливе. Наи­
большее содержание альдегидов в ОГ отмечается при работе пспрофетого двигателя
на режимах пуска и холосто! о хода, их содержание в ОГ, как правило, не нормирует­
ся, поскольку составляет сравнительно небольшую часть в общей массе токсичных
веществ. Эти оценки, а также последующие, взяты из наиболее весомых монографий
последнего десятилетия Маркова В.А. и др.; Иващенко Н.А. и др.; Звонова В.А.; Гор­
бунова В.В. и Патрахальцева Н.Н.: Кульчицкого А.Р.; Морозова К.А.
Концентрация монооксида углерода СО в ОГ двигателя внутреннего сгорания
может достигать значительных величин. Так, на режимах холодного холостого хода
при значительном избытке топлива объемная доля СО в ОГ достигает 6 - 12 % . Мо­
нооксид углерода СО окисляется в атмосфере до нетоксичного диоксида углерода
СОг в течение 2 - 4 0 месяцев. Из оксидов азота, содержащихся в ОГ, 99 % объема
приходится на монооксид NO. В атмосфере он окисляется до менее токсичного NO2 в
течение от 1 до 100 часов в зависимости от условий окисления.
Следует отметить последовательное ужесточение норм токсичных выбросов в
Европе в течение последних десятилетий, причем нормы RURO-4, соответствуют
очень низко.му содержанию токсичных веществ в ОГ - для легковых автомобилей с
расходом топлива ~ 7 л /100 км предельные объемные доли гсо = 0,11%; г^о =" 0,0083
% ; Усн = 0,010 % . Такое содержание токсических веществ в ОГ можно обеспечить
лишь при высоком уровне организации горения и рабочего процесса и, конечно, при
наличии в выпускной магистрали специальных нейтрализаторов.
В течение последних десятилетий общепринятое математическое описание ра­
бочего процесса и горения в ДВС является нульмерным (термодинамическим), и по­
луэмпирическим (модель Вибе). Оно предполагает наличие однородных, или, в слу­
чае многозонных моделей, ступенчато-однородных полей температуры, концентра­
ций компонентов смеси и, тем самым, коэффициента избытка воздуха горючей смеси
в камере сгорания (КС). Такая идеализация процесса суживает диапазон возмож1Юстей для снижения токсичности ОГ.
Чтобы избежать этого, необходимо использовать более гибкие системы смесе­
образования, такие как программируемый (неоднократный) впрыск топлива непо­
средственно в цилиндр, позволяющий получать в КС неоднородные и наиболее вы­
годные поля коэффициента избытка воздуха а перед зажиганием. Такие поля должны
обеспечивать надежное воспламенение в окрестности свечи зажигания (где смесь
должна быть достаточно богатой), отсутствие детонации к концу горения и малый
общий расход топлива (то есть в среднем по объему КС высокие значения аср=
1...2,5). Известны работы академика РАН В.Е. Алемасова и его учеников А.Л. Абдуллина, А.В. Демина, В.Г. Крюкова, В.И. Наумова, учитывающие гепломассообмен
и неравновесие химического состава в HeolH|igwi4i.j^lJU;iM "Имоточных камер сго-
I
SNMMOTEKA
I
i'"33W/'5 J
рания; работы Д.Д. Матиевского, П.К. Сеначина, М.Ю. Свердлова, С В . Пешкова по
изучению процессов самовоспламенения топлив в ДВС на основе химической кине­
тики, а также работы многих других ученых.
В то же время, до сих пор в отечественной литературе отсутствуют решения
задачи для процесса горения в ДВС, как для химически реагирующей турбулентной
смеси газов в камере сгорания и полости цилиндра изменяющейся сложной геомет­
рии. В связи с этим сформулируем следующее.
Целью работы является повышение экологической чистоты рабочих процес­
сов в ДВС с искровым зажиганием по данным теоретического прогнозирования со­
держания оксидов углерода и азота в отработавших газах, проводимого на основе
уравнений тепломассообмена и химической кинетики для процессов горения и турбу­
лентного догорания в цилиндре двигателя.
Для достижения поставленной цели необходимо решение следующих задач:
1) моделирование нульмерной химической кинетики горения для потока, пересекаю­
щего фронт пламени в ДВС; определение температуры и состава продуктов горе­
ния на выходе фронта (продуктов оттока);
2) разработка, с использованием методов химической кинетики, математического
описания процесса тепломассообмена для описания догорания продуктов оттока
за фронтом горения и в процессе расширения смеси; определение кинетического
механизма догорания продуктов; оценка влияния турбулентности на константы
скорости реакций;
3) моделирование турбулентного догорания продуктов оттока в объеме камеры сго­
рания (и цилиндра), анализ расчетного содержания оксидов углерода и азота в от­
работавших газах и сравнение с опубликованными экспериментальными данными.
Решение перечисленных задач получено для перспективного ДВС Mitsubishi
Galant (1996 г.) с непосредственным впрыском топлива (Gasoline Direct Injection - в
дальнейшем GDI) и двигателя ВАЗ-21011, для которого имеется обширный банк
опытных данных. В работе использованы известные методики и программы для рас­
чета: турбулентных полей скорости в полости цилиндра в процессах горения и рас­
ширения (программа GAS-2, проф. Дунаев В.А., ТулГУ); полей коэффициента избыт­
ка воздуха в камере сгорания и движения фронта горения в неоднородной горючей
смеси (доцент Чесноков С.А.).
Научная новизна результатов работы заключается в следующем:
- предложен сокращенный летальный кинетический механизм реакций догорания
продуктов оттока за фронтом горения в условиях ДВС (52 реакции для 14 частиц);
- определены поправки к константам скорости химических реакций, учитывающие
турбулентные пульсации температур;
- на основе уравнений тепломассообмена и химической кинетики горения разрабо­
тана математическая модель для описания процесса турбулентного догорания за
фронтом пламени и при расширении продуктов сгорания в ДВС; по результатам рас­
четов проведено сравнение с известными экспериментальными данными.
Методы исследования базируются на основных положениях химической ки­
нетики и механики гомогенных смесей газов; в работе используются приемы матема­
тического анализа, методы решений уравнений математической физики, а также ма­
тематическое моделирование на основе численных решений систем дифференциаль­
ных уравнений.
Научная значимость работы заключается в решении задачи турбулентного
тепломассообмена на основе химической кинетики при догорании продуктов оттока
за фронтом горения, как основного процесса, определяющего образование токсичных
оксидов в д в е .
Практическая значимость работы заключается в разработке эффективных
алгоритмов и профамм расчета на Э В М , позволяющих произвести расчеты содержа­
ния оксидов углерода и азота при горении, догорании и выхлопе из двигателя.
Реализация результатов работы. Результаты работы внедрены
- в виде программ расчета содержания оксидов углерода и азота в отработав­
ших газах двигателей внутреннего сгорания (акт внедрения в практику ОАО А К "Туламашзавод" от 16.08.05);
- в учебный процесс ТулГУ и используется в лекционных курсах «Теплотехни­
ка», а так же при выполнении курсового проектирования по специальности 101200
(акт внедрения от 12.09.2005).
Апробация работы. Результаты исследования докладывались на X X X I X Меж­
дународной научно-технической конференции. Секция "Поршневые и газотурбинные
двигатели" (МГТУ "МАМИ", 2002 г.), НТК "Луканинские чтения. Проблемы и пер­
спективы развития автотранспортного комплекса" (МАДИ (ГТУ), 2003 г.), 3-й Все­
российской научно-технической конференции "Современные тенденции развития авгомобилестроения в России" (Тольятти; Тольят. гос. ун-т, 2004 г.). Всероссийской на­
учной конференции студентов и аспирантов «Молодые исследователи-региону» (Во­
логда: ВоГТУ, 2004 г.), 1 -й всероссийской научно-технической конференции студен­
тов и аспирантов «Идеи молодых - новой России» (Тула: ТулГУ, 2004 г.).
Структура и объем работы. Диссертация состоит из введения, четырех глав,
заключения, списка литературы и пршюжения. Общий объем работы составляет
страницы, из них 133 страниц основного текста, 35 рисунка, 41 таблиц и список лите­
ратуры, включающий 94 источника.
Основные положения, выносимые на защиту:
- детальный кинетический механизм химических реакций в ДВС при догорании
продуктов оттока за фронтом горения;
- методика оценки влияния турбулентных пульсаций температур на величину констан i химических реакций;
- математическая модель турбулентного догорания, полученная для описания про­
цессов догорания и расширения; модель догорания при выхлопе;
- результаты расчетов этих процессов и сравнение с известными экспериментальны­
ми данными.
ОСНОВНОЕ СОДЕРЖАНИЕ Р А Б О Т Ы
Сокращения:
ГС - горючая смесь.
ОКМ - общий кинетический механизм
ДКМ-детальный кинетический механизм
(266 реакций, 33 частицы)
ДКМД-ДКМ догорания (26реа<ш1Й, Мчаспщ) ОК-основные компоненты (доля иоторых>0,1%)
КС - камера сгорания
ПС - продукты сгорания
КТО - коэффициент i урбулентного обмена ТД - турбулентное догорание
ОГ - отработавшие газы двигателя.
GDI - Gasoline Direct Iniection - ДВС с впрыском топлива в Г1илиндр (Mitsubishi Galant,1996)
Основные обозначения:
а - коэффициент температуропроводности, м^/с, Q - выделение тепла при горении, Дж/кг,
D - коэффициент диффузии, м/с;
г, - мольная (объемная) доля компонента,
Е - энергия активации. Дж/моль
Sj - скорость химической реакции,
gi - массовая доля компонента;
Тс - температура горения. К,
HI - тепловой эффект реакции, Дж/моль,
U^ - скорость выделения тепла, Вт/м',
kj - константа скорости реакции,
u,w,v- проекции скорости, м/с;
Л/ - кратность коэфф-та турбулентного обмена, «г- турбулентная скорость горения, м/с,
р-давление, Па;
а - коэффициент избытка воздуха топливной смеси,
fi - молярная масса, кг/моль,
V'- кинематический коэффициент вязкости, м^/с;
р - плотность смеси, кг/м"',
PI - распределенная (парциальная) плотность компонента, кг/м ,
Нижние индексы
F - фронт горения, G - горение, (- номер компонента, у - номер реакции, w - граница, н - на­
чальное, см смесь, Г - турбулентное
Верхние индексы: р - прямая реакция; о - обратная
Нумерация разделов, формул, рисунков и таблиц сохранена такой тке как в диссертаiiuu Ссылки на публикации доцента Чеснокова С А. обозначены как [Ч]
Во введении обоснована актуальность работы, сформулированы цели и задачи ис­
следования, рассмотрены научное и практическое значение диссертационной работы.
1. Обзор литературных данных.
Среди рассмотренных работ ведущее положение занимают публикации Инсти­
тута химической физики Р А Н , в одной из которых полностью опубликован Д К М го­
рения метана в условиях Д В С проф. Басевича В.Я. (256 реакций для 29 частиц). Этот
Д К М положен в основу используемых ниже кинетических механизмов. Данные по
кинетике тяжелых частиц, образующихся при распаде бензина в начале горения в у с ­
ловиях низких температур, малочисленны - полный детальный механизм распада ос­
новных фракций неизвестен. Это обстоятельство делает невозможным, в частности,
определение толщины слоя затухания пламени вблизи холодной стенки камеры сго­
рания и, тем самым, оценку содержания несгоревших углеводородов в П С .
По известным литературным данным в перспективных двигателях G D I
Mitsubishi Galant (1996 г.), Pajero Pinin и Carisma (1997 г) применяются два варианта
впрыска - основной для работы на стехиометрической и форсированный для сверх­
бедной смеси. В первом случае впрыск производится широким коническим факелом
в течение такта впуска (рис. 1.4). Длительное, в течение полутора тактов, испарение
способствует образованию гомогенной (но неоднородной) смеси (Оср* 1). Поршень
двигателя имеет сферическую выемку, что обеспечивает разворот струи воздуха, на­
сыщенного парами бензина, в сторону свечи зажигания и, тем самым, обогащенную
смесь в ее окрестности.
впуск, клапан
свецд
Рис. 1.4. Схема конвектив­
ных процессов при всасыва­
нии - впрыске.
Рис.3.1. Зона догорания про­
дуктов оттока от фронта дви­
гателя G D I
рис. 1.4
оис.3.1
Во втором случае впрыск производится узким факелом в конце сжатия, непо­
средственно перед зажиганием. Для этого варианта впрыска характерно одновремен­
ное протекание процессов испарения частиц, горения смеси и перемешивания про­
дуктов сгорания, что па сегодня практически делает невозможным моделирование
процесса в целом. В диссертации рассмотрен первый вариант впрыска (рис. 1.4), для
которого разделены во времени процессы испарения частиц, горения смеси и пере­
мешивания ПС.
2.Нульмерная химическая кинетика горения. Преаварительное тестирование и оценки
В качестве модельного горючего, заменяющего бензин, в данной работе ис­
пользован н-октан CgHis, что представляется оправданным в связи с близкими значе­
ниями характеристик oicrana и бензина, используемых в кинетических расчетах. Ме­
ханизм термического распада октана соответствует опубликованным для различных
тяжелых углеводородов (рис. 2.3):
С8Н18-*СзНб+2С2Н4+СН4,
СзНб^С2Н4+0,5С2Н2+Н.
Выбор скорости распада н-октана, как для тетрадекана обусловлен отсутствием
литературных данных по распаду н-октана и требует оценки влияния этого процесса
на дальнейший кинетический процесс.
В рамках проведенного тестирования по образованию оксида азота NO показа­
на необходимость совместного учета трех реакций механизма Зельдовича Я.Б. и тер­
мической реакции окисления азота Nj, предложепной для условий ДВС Звоновым
В.Л. (Токсичность двигателей внутреннего сгорания. М.: Машиностроение. 1981):
N + N 0 о N2 + О, N + 02 <^ N 0 + О; N + ОН <^ NO + Н; N2 + Oj <^ NO +N0.
На основе Д К М горения легких углеводородов в условиях ДВС Басевича В.Я. с
учетом реакций образования NO и эмпирических реакций распада октана, как мо­
дельного топлива заменяющего бензин, образован общий кинетический механизм
(ОКМ) горения топлива, содержащий 266 реакций для 33 частиц. Механизм уравно­
вешен для условий воспламенения и горения в ДВС: при температурах Т = 1000,
2000, 3000 К и, соответственно, давлениях р = 2; 3,5; 5 МПа.
Применение метода Гира для решения кинетических уравнений нульмерной
задачи позволяет уменьшить маишнное время в 10' раз, сократив его от нескольких
суток до нескольких секунд (см. Чесноков С.А. Моделирование высокотемператур­
ных реакций горения. Тула: Изд-во ТулГУ, 2002).
По многочисленным оценкам в ряде работ турбулентный фронт горения ДВС
имеет микроламинарную структуру. Для потока, пересекающего ламинарный участок
фронта горения на основе ОКМ, решается система уравнений химической кинетики и
уравнение энергии
1=?^.'
^Ar)f.i2:«^„
S.j = ±к^г„Гр ;
Со =
=
-X''^''"^■' '" = S''''"'.
(2.17)
(2.19)
Здесь а, Р- различные значения i; ц - текущая средняя молярная масса смеси; d"p,(T) —
мольные теплоемкости основных компонентов, аппроксимация которых проведена по
Результаты решения
1. Расчет процессагорения,описываемый нульмерной кинетикой ОКМ удается
"запустить" при начальной температуре горючей смеси не ниже 1100 - 1200К (рис.
2.1). Причиной является то, что подогрев смеси с 600-700 К (температура в цилиндре
д в е в конце сжатия) до 1200 К происходит за счет неучтенных низкотемпературных
реакций тяжелых молекул. Для того, чтобы конечная температура горения соответст­
вовала экспериментальным данным (2800...3000 К ) эта неточность компенсируется
завышением удельной теплоемкости смеси температурной поправкой. Ниже (п. 3) ус­
тановлено, что конечный состав смеси соответствует равновесному состоянию при
температуре горения (кроме N0, доля которого мала), поэтому повышение тешюемкости влияет лишь на промежуточную кинетику реакций, не искажая конечный со­
став продуктов горения. График изменения во времени скорости тепловыделения за
счет химических реакций Uy имеет два характерных максимума - в диссертации об­
суждается возможный механизм их образования.
т,к
2800
Uv
К"^ ^°
а=1
lO'^i
2000
Imin
1200
1000
^
' V
- l /
31
/"
/
;
Uv
*
'
Ю т
ми(^
Рис. 2.1 - Температура и скорость
тепловыделения для потока, пересекаюшего фронт горения (р=4 мПа)
Г)
'10-2
С2Н4ч
^ ^ ' ^
■
^
10^
'
C2H2
10-Ч
^
Ч
10"»
a=l
y. H
'.^^f^^^^y'""^
A
^
^^ \
cfMv^--/-"
10-»
Гн
^■^r
C3H6 \
10-i
1(Г*
\ \
,
'--'
la'
lin-»
T, с
Рис. 2.3 - Мольные доли компонентов
при распаде октана в начале процесса
(р = 4МПа)
В начале процесса (при т < 1 мкс) расчетная скорость тепловыделения за
счс! химических реакций является существенно отрицательной из-за сильного по­
глощения тепла реакциями распада октана. В значителыюй степени отрицательное
тепловыделение реакций распада компенсируется при сгорании атомарного водо­
рода Н (рис. 2.3). В итоге временное снижение температуры смеси составляет 100200 К. Изменение скорости распада октана в широком интервале в большую или
меньшую стороны (рис. 2.2) совершенно не влияет на температуру и кинетику про­
цесса в течение основного периода (при т > 1 мкс), что оправдывает достаточно
произвольный выбор констант скорости распада.
Метин СН является типичным представителем углеводородной группы радикшюв,
концентрация которых достигает максимума уже в начале процесса, после чего на­
чинается "развитое" горение смеси. Наибольшее содержание гидроксила ОН на­
блюдается в конце горения, его доля определяет ведущий процесс этого периода выгорание оксида СО. В настоящее время неизвестно соотношение концентраций
возбужденных радикалов СН*, ОН* и "спокойных" радикалов СН, ОН, доля кото­
рых определяется в кинетических расчетах, и это соотношение может быть суще­
ственно нелинейным. Тем более показательным и подтверждающим достоверность
кинетики Басевича В.Я является качественное соответствие экспериментальных и
расчетных (кинетических) данных рис. 2.4 для режима горячего холостого хода
ВАЗ-21011.
3. П о данным кинетических расчетов с различными а = 0,8 . . . 1,5 при горении
смеси конечное содержание основных компонентов смеси соответствует равно­
весному состоянию при температуре горения, с отклонениями не более 12 % . С о ­
держание оксида азота N 0 при горении является неравновесным. Ниже показано,
что конечное состояние при расширении смеси в двигателе является существен­
но неравновесным. Полученная при расчетах температура горения, тепловыделе­
ние и экспериментально оцененная [ Ч ] турбулентная скорость горения (рис. 2.5)
имеют максимальные значения при а ~ 0,9, что соответствует известным опытным
данным. Температура горения и состав продуктов на выходе фронта горения будут
использованы в гл. 3 и 4 как фапичные условия задачи Т Д .
4. На основе Д К М Басевича В . Я . для температуры горения при сравнительной оценке
скоростей реакций в уравнениях кинетики получен сокращенный д е т а л ь н ы й к и ­
нетический механизм д л я догорания смеси ( Д К М Д ) при ее расширении в ци­
линдре, включающий 52 реакции (из 266) для 14 компонентов (из 33). О н будет
использован при решении многомерных задач Т Д в гл. 3 и 4. Изменение констант
отдельных реакций Д К М Д в пределах, определяемых по данным различных лите­
ратурных источников, приводит к отклонениям в содержании основных компо­
нентов смеси не более, чем на 25 % .
Т,К f\2
1100
1000
а=1
tj
Г
а=1
10-^
С8Н18
И
10-^
3
iF
io^ 10-* х,с
\
^
10^ \
10"'
.
\
Iff"
б
'
Т,с
Рис. 2.2 - Температура (а), содержание октана (б) и атомарного водорода ( в )
в начале процесса при различной скорости распада октана:
1 - номинальной: 2 - уменьшенной на два порядка; 3 - увеличенной на порядок
Q, кДж/кг
2500
2000
1.5 а
Рис. 2. 4-Расчетные графики относительной
доли радикалов и экспериментальные точки
[Ч] относительной интенсивности их излу­
чения в возбужденном состоянии:
Д-ОН*; П - С Н *
Рис. 2.5 - Расчетные графики
температуры горения и тепловы­
деления, а также эксперимен­
тальной [Ч] турбулентной скоро­
сти горения для ВЛЗ-21011
10
3. Турбулентное догорание в ДВС. Общая постановка и решение тестовой
(двухмерной) задачи
По результатам расчета с помощью программы GAS-2 получены плоские тур­
булентные поля скорости нереагирующей смеси в периоды тактов наполнения и сжа­
тия, а также при горении и расширении продуктов сгорания:
для двигателя GDI - в вертикальном диаметральном сечении цилиндра (и КС);
для двигателя ВАЗ-21011 - в характерных горизонтальных сечениях КС и циливдра
Результаты расчета для обоих двигателей показывают, что к концу горения, ря­
дом со струей оттока продуктов от фронта горения, со стороны наиболее динамично­
го сокращения длины фронта образуется вихревая зона, определяющая затем общую
циркуляцию газов в процессе расширюния.
Задача турбулентного догорания. В гл. 2 на основе ДКМ Басевича В Я . опре­
делен сокращенный детальный кинетический механизм для догорания смеси
(ДКМД) за фронтом пламени (табл. 3.3), реакции которого уравновешены при пара­
метрах, характерных для процесса расширения: Т = 3000 / 2000 /1000 К и, соответст­
венно, р = 5 / 3 / 0,5 МПа. Для учета влияния турбулентных пульсаций температур на
константы химических реакций к введен турбулентный коэффициент к > 1. Его
значения определены для гармонических колебаний температуры, амплитуда которых
по данным спектральных измерений в КС двигателя ВАЗ-21011 составляет примерно
20 %. Результаты представлены полиномом, что позволяет вычислять эффективные
значения констант скорости реакций к(Т)= кАТехр(-Е/КТ) в процессе решения задачи.
Таблица 3.3 - Ведущие обратимые реакции догорания (ДКМД)
№
1
2
3
4
5
6
7
8
9
Реакции
ОН+Н2=Н+Н20
0Н+0=Н+02
0Н+Н=0+Н2
ОН+ОНЮ+НгО
0Н+Н=Н20
Н+Н02=Н2+02
Н+02=Н20
н+но2=он ю н
0+Н02=02+ОН
№
10
11
12
13
14
15
16
17
18
Реакции
ОН+Н02=02+Н20
ОН+Н202=Н02+Н20
ОН+ОН=Н202
0+Н202=02+Н20
Н02+Н02=Н202+02
CO+0H=C02+H
СО+Н02=С02+ОН
С0+ОС02
СО+02=С02+0
№
19
20
21
22
23
24
25
26
Реакции
Н+ССИНСО
НСО+02=Н02+СО
0Н+НСО=СО+Н20
Н+НСО=Н2+СО
N+N0=N2+0
N+02=N0+0
N+OH=NO+H
N2+02=NO+NO
Основные допущения задачи турбулентного догорания.
1. В условиях высоких давлений (~ 50 атм) температуры компонентов непре­
рывно выравниваются за счет эффективного обмена энергией при многочисленных
соударениях молекул, и внутри каждого компонента устанавливаются общие для всей
смеси поля температур. Тем самым считаем, что, несмотря на интенсивные химиче­
ские реакции, равновесное распределение молекул по энергиям не нарушается.
2. В связи с небольшой скоростью течения (менее 50 м/с) в зоне продуктов
сгорания пренебрегаем кинетической энергией газа по сравнению с внутренней энер­
гией Кроме того, учитывая сравнительно небольшие различия температур в пределах
поля решения теплоемкость смеси считаем постоянной Излучением пренебрегаем.
3. Эффекты термо- и бародиффузии, а также диффузионную теплопровод­
ность и перенос энергии диффузионным потоком не учитываем. Работой диссипативных сил пренебрегаем. Обратимая работа сил давления вычисляется только в период
интенсивного расширения продуктов сгорания.
11
4. Поля турбулентных коэффициентов переноса определяем по данным реше­
ния транспортных уравнений к-е модели турбулентности [Ч] на основе полей скоро­
сти, полученных с помощью программы GAS-2. В трехмерных задачах, для сокраще­
ния машинного времени, поля коэффициентов переноса задаем как однородные и из­
меняющиеся во времени (рис. 3.6). Считаем, что изменение турбулентного коэффи­
циента температуропроводности смеси а г определяется только турбулентной тепло­
проводностью Хъ поскольку объемная теплоемкость ср по сравнению с 1г меняется
незначительно,
5. Догорание горючих компонентов смеси за фронтом пламени характеризу­
ется 52 ведущими реакциями для 14 компонентов (табл. 3.3).
6. Компоненты и смесь в целом являются идеальны.ми газами.
При этих допущениях уравнения сохранения энергии смеси и массы компонен­
тов, полученные на основе математического описания для гомогенной смеси из моно­
графии Нигматулина Р.И. (Динамика многофазных сред. М.: Наука, 1987), имеют вид:
дТ
дТ
VU —
дТ
+W
дТ
-^-V
дт
дх
ду
8gi
dz
dgi
дт
дх
ду
д:
дТЛ
д_
Jx
1
p
^ дх)
д
ду
1
I^;(S;-S;)H
ду
дГ
dz
^1
Hi^'
Og,
ЦрОг
дхК:
дх
0=1,2, ...,N).
Источники теплоты и масс компонентов:
дТ_
ду
dz
+ F„
^Д^^
(3.29)
+ F„
(3.30)
Р dp
F^,=./i,|;(±s;¥5;).
(3.31)
_;=1
P дт
"« J
J=l
Турбулентные характеристики переноса определяются среднеобъемной кратностью
КТО N-1 (как отношении КТО и ламинарной вязкости), поэтому: a-r=a„a.,i Nf,
Di=Dia.M Nj. Значение Л^у определяется по данным рис. 3.6.
Для скорости реакции имеем:
Ря =
S^'^k'/'g&liMA,
k'^'^KArexpi-E/RT).
(3 34)
Система (3.29)-(3.34) включает 15 уравнений и позволяет произвести расчет полей
температуры и массовых долей 14-ти компонентов: Н, Н^, ОН, НОг, Н2О2, НгО, НСО,
0,02,CO,C02,NO,N,N2.
При решении использованы: метод переменного порядка точности Гира для
решения уравнений химической кинетики, метод расщепления шага по времени для
много.мерных задач тепломассообмена, метод прогонки для одномерных решений
(алгоритмический язык - Фортран-95).
Тестирование двухмерных задач турбулентного догорания проведено на
примере первого впрыска рабочего режима двигателя GDI (п = 3000 об/мин, p„aicc=
60 10' Па). Камера сгорания этого двигателя, образованная выемкой в поршне, и оба
клапана расположены в диаметральной плоскости цилиндра (рис. 1.4), поэтому двух­
мерная постановка задачи ТД является достаточно обоснованной. Приняты парамет­
ры двигателя, характерные для легкового автомобиля: диаметр цилиндра 0,08 м, ра­
бочий объем 0,37310'' м', степень сжатия 12, опережение зажигания 18 °пкв. В каче­
стве исходных данных используются:
- радиус начального сферического объема продуктов сгорания вокруг точки зажига­
ния 6 мм; начальное давление в этом объеме р„ =2 МПа, начальная температура
Т.= 2900 К;
- ламинарный коэффициент температуропроводности смеси «,„,=0,1910''' MVC И
средний ламинарный коэффициент диффузии дтякомгюненговсмеси Ди„=0,25-10''' м^/с;
- геометрия поля решения (форма и перемещение фронта горения и поршня); поля
скорости и коэффициента избьп'ка воздуха а в КС (рис. 3.9); а также давление в
функции времени по результатам расчета рабочего процесса;
- параметры констант скорости и тепловые эффекты реакций ДКМД (табл 3.3);
- температуры фронта горения То и состав продуктов оттока для различных а (гл.2).
За.мена полученных в работах [Ч] полей кратности КТО средним объемным
значением Л'/, изменяющимся во времени, приводит к погрешностям содержания ос­
новных ко.мпонентов к концу расширения в пределах 5 %. Поэтому изменение сред­
ней кратности КТО во времени процесса, полученное для обоих двигателей (рис. 3.6),
будет использовано в трехмерных задачах ТД. Это позволяет уменьшить время реше­
ния на Э В М без заметного снижения точности.
Поля содержания компонентов ПС анализируются в следующей главе, тем не
менее отметим следующее. По сравнению с реакциями во фронте горения для реак­
ций догорания водорода и углеводородов характерна относительная уравновешен­
ность - близость скоростей прямых и обратньгх реакций. Это определяет сильное
влияние содержания компонентов на результирующую скорость реакций (как раз­
ности скоростей прямой и обратной реакций) и, тем самым, на интенсивность источ­
ников массы компонентов и теплоты. Тем не менее, все рассмотренные реакции при
догорании в камере сгорания ДВС не являются полностью уравновешенными. Мед­
ленные реакции окисления азота, скорость которых ~10°...10 моль! / (моль с), и бы­
стрые реакции окисления водорода и углерода (~10^...10*) определяют результирую­
щие скорости близких порядков ~10' ...10*. Результирующие скорости образования
радикалов в значительной степени компенсируют друг друга вследствие их акгивносги.
S
А-А.
Рис. 3.9 - Поля коэффициента избыт­ Рис. 3.10 - Поля коэффициента избытка воздуха
ка воздуха в КС перед зажиганием в КС перед зажиганием (п = 3000 об/мин) для
(п = 3000 об/мин) для двигателя GDI
двигателя ВАЗ 21011 [Ч]
Проверка на точность численного решения задачи ТД показала, что уменьше­
ние шага по координатам в 10 раз изменяет содержание основных компонентов не
более, чем на .5 %, а для остальных компонентов- не более, чем на 20 %. Уменьшение
на порядок шага по времени приводит к меньшим погрешностям, ~ I . 2 % . Учиты­
вая, что время решения двухмерной задачи на ПЭВМ Pentium-4 (1,9 ГГц) возрастает
при этом с 0,5 часа до 60 часов, исходная величина шагов может считаться опти­
мальной.
поле решения
496 ф."пкв
Рис. 3 6 - Изменение средней по объему кратности КТО
в период горения (G) и расширения смеси для двигателей.
1 - GDI, 2 - ВАЗ, рабочий ход; 3 - ВАЗ, холостой ход [Ч]
Рис. 4.1 - Схема тепло­
обмена на границе поля
решения
4 Турбулентное догорание в ДВС. Трехмерные задачи и анализ
Двигатель G D I
Горение и расширение. Система уравнений ТД (3.29)-(3,34) решена в трехмер­
ной постановке. Параметры двигателя GDI (рис. 1.4) и необходимые исходные дан­
ные приведены в гл. 3. Термическое сопротивление ламинарного подслоя толщиной
5л было учтено в виде простого соотношения стационарной теплопередачи в форме
граничного условия 4-го рода (рис. 4.1). Из равенства тепловых потоков на границе
поля решения (qj = Чл) температура на границе
Tw=(To+ATcp)/(l+A), где А = 8л N T / L T .
(4.15)
Примем в первом приближении, что относительная ширина поля LT / 5л = const = Cw,
тогда А = N T / Cw- Кратность КТО N j уменьшается во времени процесса в соответ­
ствии с рис. 3.6. Значение коэффициента граничных условий Cw подбиралось при
тестовых решениях задачи ТД таким образом, чтобы температура смеси в начале вы­
пускной магистрали соответствовала экспериментальной. Оценим толщину ламинар­
ного подслоя 5л. Считая, что L^ = Н/4, где Н ~ 0,025 м - высота камеры сгорания дви­
гателя GDI, получим 5л = L T / Cw = 0,006/100 = 0,6-10"* м, что хорошо соответствует
оценкам в монографии Петриченко P.M. (1983), где 5л ~ Ю"^ м.
Выхлоп и выпуск. При открытии выпускного клапана большая часть продуктов
сгорания быстро выбрасывается в выпускную магистраль (выхлоп) и затем, вместе с
принудительно вытесняемыми газами, уходит через нее в атмосферу (вьшуск). Учи­
тывая, что основное (> 90 % ) изменение содержания СО (и N 0 ) наблюдается в про­
цессах горения и расширения, рассмотрена упрощенная постановка задачи выхлопа и
выпуска.
Предположим, что быстрый процесс выхлопа (продолжительностью < 2 мс) яв­
ляется внешне адиабатическим процессом, происходящим без теплообмена с поверх­
ностью цилиндра, поршня и клапана. Примем, что при выхлопе смесь мгновенно пе­
ремешивается и становится идеально однородной. При этом в ней ускоряются
химические реакции догорания, в частности, окисление оксида углерода с выделени­
ем тепла. В процессе выхлопа и после его окончания, в течение -100 мс, смесь пере­
мещается по выпускной магистрали, где реакции продолжаются. Процессы выхлопа и
выпуска можно описать в рамках нульмерной кинетики реакций догоршшя для пото­
ка, вычисляя по известной из эксперимента кривой давления на каждом шаге по вре-
мени адиабатную температуру и выделение геплоты химических реакций в течение
шага, повышающее эту температуру.
Анализ результатов. К концу горения турбулентное поле скорости, получен­
ное с помощью программы GAS-2 (рис. 4.3,а), имеет вихревую зону с центром С].
Ломаной линией показан фронт горения F Пунктиром выделены зоны источников
теплоты или массы компонента, положительные + или отрицательные -. Температура
2230 К (б) является температурой границы поля решения Т „ (рис. 4.1). Вблизи стенок
камеры сгорания формируется турбулентный пограничный слой с температурами
2230-2500 К.
Для нескольких контрольных точек поля решения проведен анализ скоростей
химических реакций Низкие температуры слоя определяют меньшие скорости реак­
ций, но из-за их неуравновешенности суммарные скорости образования большинства
компонентов в несколько раз больше, чем в объеме КС. Наибольпгая концентрация
окиси углерода СО (г) сохраняется в верхней части КС, где горела богатая топливная
смесь. В основном поле камеры доля СО снижается во времени под влиянием веду­
щих реакций догорания: СО+ОН-^СОа+Н, СО+НОг-уСОг+ОН, СО+О-^СОз и
Н+СО->НСО. В пограничном слое для первых двух реакций результируюпше скоро­
сти возрастают в несколько раз, в результате чего выгорание СО ускоряется.
Максимальное содержание оксида азота N 0 (д) наблюдается в зоне сгоревшей
стсхиометрической горючей смеси. В области высоких температур все реакции
N2+0->NO+N, N+Oj^NO+O, N+0H->N0+H и N2+02->NO+NO интенсивно генери­
руют оксид азота. Вдали от фронта горения концентрация N 0 уменьшается под воз­
действием ведущих реакций раскисления азота. N+NO-->N2+0, H+NO—>OH+N. Это
уменьшение значительно усиливается в пограничном слое за счет реакции
N0+N0->N2+O2.
В процессе расширения смеси (рис. 4.5) помимо основного вихря с центром
С] (а) образуется вихревая зона Сг с обратным вращением. Максимум температуры Т
(б) наблюдается в центре основного вихря. В вихревых зонах образуются зоны отри­
цательного (-) тепловыделения за счет преобладания эндотермических реакций.
Рис. 4.3 - Поля скорости (а), гемпературы (б) и массовых долей компонентов (в-д)
в КС двигателя GDI для конца горения (т =1 мс; R =20%; п =3000об/мин; ф=370°пкв)
IS
Поля радикалов ОН, О, НОг и, в какой-то степени, оксида углерода СО форми­
руются быстрыми химическими реакциями с участием этих компонентов. В цешре
вихрей конвективное перемешивание отсутствует, реакции расходования радикалов
замедляются, что способствует сохранению их повышенной концентрации. В отличие
от радикалов поля молекулярных компонентов; кислорода Ог (и N0), СОг, НгО и др.
определяются в основном конвекцией и турбулентной диффузией. Для всех ком­
понентов результирующие скорости ведущих реакций максимальны в зоне высо­
ких температур. В турбулентном пограничном слое, в отличие от периода горения,
значения результирующих скоростей реакций на порядок и более ниже, тем не менее
суммарная скорость расходования каждого компонента (кроме N0) по абсолютной
величине на порядок выше. Это определяет, в частности, интенсивное выгорание ок­
сида углерода СО в пограничном слое. Заметное раскисление оксида азота N 0 на­
блюдается в зоне высоких температур, в течение периода, когда они превышают 2000
К. В пограничном холодном слое суммарная скорость реакций раскисления в десятки
раз меньше.
Методом пробных решений задачи ТД получен оптимальный режим работы
двигателя GDI (обеспечивающий низкое содержание СО, N 0 ) при среднеобъемном
значении коэффициента избытка воздуха а^ = 1,00, изменении а в процессе выгора­
ния неоднородной топливной смеси от 0,90 до 1,06 и применении рециркуляции ОГ в
объеме 20 % . Мольное содержание компонентов в ОГ составляет: оксид углерода СО
- 0,299 % , оксид азота N 0 - 0,252 % , кислород Oj - 0,288 % .
Рис. 4.5 - Поля скорости (а), температуры (б) и массовых долей компонентов (г-з)
в диаметральном сечении цилиндра двигателя GDI для середины процесса
расширения (т = 4 мс; R = 20%; п = 3000 об/мин; ф = 424°пкв)
Расчеты при постоянных в полости КС зггачениях коэффициентах избытка воз­
духа (однородное поле а) в диапазоне а = 0,8... 1,2 для рециркуляции R -= 20 % (рис.
4.6,а) показали удовлетворительное соответствие остаточного содержания кислорода
Оз, оксида углерода СО и оксида азота N 0 известным экспериментальным данным
Полные данные по содержанию СО и NO в ОГ" (рис. 4.6,6), позволяют произве­
сти для данного ДВС предварительный прогноз уровня токсичности О Г по извест­
ным характеристики поля а перед зажиганием- среднему в объеме КС коэффициенту
избытка воздуха асн и неоднородности поля Да = ОМАХ - «мш- Возможен обратный
подход - при заданной токсичности ОГ определить необходимые Оср и А«- Зги па­
раметры можно варьировать, регулируя систему впрыска топлива. В частное ги, неод­
нородность поля Да можно из.менять за счет продолжительности впрыска топлива
[Ч].
Г]
t
% Vco
Ог//
V/
/ //
//
/ /
/ /
/ /
о
0,8
0,9
1
N0
0
1,1
а
Рис 4 6,а - Объемное содержание компонентов
в отработавших газах при различном избытке
воздуха в однородных полях а.
—»— результаты расчетов,
экспериментальные данные монографии'
Вырубов Я Н, Ивашснко Н.А., Ивин В.И. и др.
«ср
1
1,1 0,32
Рис 4 6,6 - Объемное содержание ком­
понентов в ОГ двигателя GDI при раз­
личных значениях' среднего в объеме КС
коэффициента избытка воздуха аср
и неоднородности поля Да
(R = 20 %, п = 3000 об/мин)
Двигатель ВАЗ-21011
Описан алгоритм решения трехмерной задачи ТД (3.29)-(3 34) для двигателя
ВЛЗ-21011. При решении используется неоднородное поле а = 0,9... 1,1 (рис. 3.10).
Изменение средней кратности КТО N^- в период расширения (рис. 3.6) получено при
использовании к-е модели турбулентности [Ч]
Анализ результатов. Hajinnne остаточных газов в количестве у = 10 % приво­
дит к уменьшению конечного содержания оксида азота N 0 на 45 %, а оксида углеро­
да СО на 4%. Причиной сильного снижения доли N 0 является уменьшение темпера­
туры горения (на 100 К). В начале горения для поля температур характерно наличие
"плато", которое в дальнейидем, по мере развития турбулентного тепломассообмена
значительно уменьшается.
Из графиков рис. 4.8 следует, что основное изменение температуры и содержа­
ния Компонентов наблюдается при горении и в течение такта расширения. При вы­
хлопе из цилиндра доля СО снижается на протяжении третьей части выпускной маги­
страли, пока температура смеси превышает 1000 К. При меньших температурах реак­
ции с участием СО "замораживаются". Этот же эффект определяет постоянное со­
держание оксида азота N 0 в смеси при температурах менее 2200...2500 К. Представ­
лены экспериментальные данные с доверительными интервалами, указанными в виде
скобок. Термопарные измерения температуры газа ( В и р ) проведены на входе в вы­
пускную магистраль, эти результаты использованы для подбора коэффициента гра­
ничных условий Cff.
8 10
а
т- ^с 100
Т, мс
Рис. 4.8 - Изменение средних в объеме цилиндра ВАЗ-210П температуры и мольных
долей основньп< компонентов во времени процесса - от горения до выпуска:
а - рабочий режим двигателя (п = 30CiO об/мин; у = О, I; рмлх - 6 МПа); б -горячий
холостой ход (п = 800 об/мин; у = 0,24; рмА\= 1 МПа);
А, В,.. .,G - данные экспериментов [Ч]
Температуры горения в режимах рабочего хода (А) и горячего холостого хода
(Е) измерялись цветовым методом. Содержание оксидов N 0 и СО в отработавших га­
зах (С, D, G) определялось серийным инфракрасным газоанализатором ИКАФ-057.
Соответствие расчетных данных экспериментальным представляется вполне удовле­
творительным.
ОСНОВНЫЕ РЕЗУЛЬТАТЫ И ВЫВОДЫ РАБОТЫ
В диссертации на основе уравнений тепломассообмена решена задача турбу­
лентного догорания (ТД) в ДВС с искровым зажиганием с целью прогнозирования
содержания оксидов углерода и азота в отработавших газах. Автору не известны ана­
логичные решения и публикации в отечественной литературе.
Химическая кинетика 1орения применялась в рамках кинетики легких углево­
дородов, как наиболее разработанной по литературным данным, и оказывающей оп­
ределяющее влияние на содержание оксидов СО и N 0 в продуктах сгорания; низко-
18
температурная кинетика тяжелых частиц такого влияния не имеет. При решении за­
дачи проведены расчеты турбулентных характеристик переноса и оценка влияния
турбулентных пульсаций гемперагуры на величину констант скорости химических
реакций.
Основные научные и практические результаты состоят в следующем:
1. Разделение общей задачи ТД и предварительное решение задач: турбулент­
ной газовой динамики (проф. Дунаев В.А.), образования горючей смеси в К С и дви­
жения фронта I орения в неоднородных полях коэффициента избытка воздуоса (проф.
Чесноков С.А.) себя полиостью оправдало. Известные поля скорости и коэффициента
избытка воздуха использовались для моделирования собственно процесса ТД. Такой
подход позволил создать эффективные математические модели и программ1юс обес­
печение, с помощью которых решение задачи ТД при использовании П Э В М Pentium4 (1,9 ГГц) требует не более 8 часов машинного времени. Альтернативным можно
считать моделирование общего процесса ТД, проведенное Риглсром У, (Германия,
Штутгартский ун-т, 1999) при затратах машинного времени до 8 недель.
2. С использованием общего кинетического механизма (266 реакций для 33
частиц), созданного на базе детального кинетического механизма (ДКМ) Басевича
В.Я., решена нульмерная задача химической кинетики (ХК). Температура горения и
тепловыделение максимальны при коэффициенте избытка воздуха а ~ 0,9, что соот­
ветствует известным экспериментальным данным. Состав продуктов сюрания на вы­
ходе фронта горения близок к равновесному, за исключением оксида азота NO (и
Ог). Применение метода переменного порядка точности Гира позволило уменьшить
время решения задач нульмерной Х К на П Э В М в 10' раз, доведя его до нескольких
секунд.
3. Разработана математическая модель процесса тепломассообмена при догора­
нии продуктов оттока за фронтом горения и в процессе расширения смеси.
Выбор сокращенного ДКМ догорания (52 реакции для 14 частиц) является дос­
таточно обоснованным - расширение схемы за счет неиспользованных реакций угле­
водородов ДКМ Басевича В.Я. не приводит к существенному изменению результатов.
Изменение констант скоростей отдельных реакций в пределах, указанных различны­
ми литературными источниками, приводит к отклонениям концентраций основных
компонентов ОГ не более, чем на 25 %.
Увеличение констан! скорости химических реакций при турбулентных пульса­
циях температуры в двигателе GDI определяет, как показывают расчеты, снижение ее
среднего по объему значения на 100. .130 К и повышение содержания оксидов угле­
рода и азота в ОГ, соответственно, на 18 и 36 %.
4. На основе гюлученных данных и ДКМ догорания решены трехмерные задачи
ТД ДJIя обоих двигателей при наличии остаточных газов. Состав смеси при догорании
продуктов сгорания в течение такта расширения - существенно не равновесный.
Для двигателя GDI методом пробных решений получен оптимальный режим
(с низким содержанием СО, N0) при среднеобъемном значении коэффициента из­
бытка воздухаttcp= 1,00, изменении а в процессе выгорания неоднородной топлив­
ной смеси от 0.90 до 1,06 и применении рециркуляции ОГ в объеме 20 % . Мольное
содержание компонентов в ОГ составляет: СО - 0.299 %, NO - 0,252 %, О2 - 0,288 % .
Расчеты для однородных гюлей а ~ const в КС показали небольшое отличие содержа­
ния оксидов углерода и азота в ОГ от известных экспериментальных данных (рис. 4.6,а).
Полные данные по содержанию СО и N 0 в ОГ (рис. 4.6,6) позволяют произве­
сти для данного ДВС прогноз уровня токсичности ОГ по известным характеристи­
кам поля коэффициента избытка воздуха а в КС перед зажиганием' среднему в объ-
еме КС коэффициенту избытка воздуха ОГР И неоднородности поля Да. Возможен и
обратный подход - по заданной токсичности ОГ определить необходимые
значения ас? и Аа, которые должна обеспечить система впрыска топлива. Не­
однородность поля Да можно изменять, варьируя, например, продолжительность
впрыска.
Для двигателя ВАЗ-21011 расчетные значения мольных долей оксидов углеро­
да и азота в ОГ (с учетом остаточных газов, рис. 4.8) соответствуют измеренным про­
мышленным датчиком ИКАФ-057 для одного из цилиндров в рабочем режиме и в ре­
жиме горячего холостого хода.
Основное содержание работы отражено в следующих публикациях:
1. Чесноков С.А., Рыбаков Г.П., Демидов М.И. Одномерная модель двухфазных тече­
ний во впускной магистрали ДВС. Тезисы докладов X X X I X Международной научнотехнической конференции. Секция "Поршневые и газотурбинные двигатели".МГТУ"МАМИ",2002.-С.56-60.
2. Чесноков С.А., Демидов М.И., Кузьмина И.В., Рыбаков Г.П. Математическая мо­
дель смесеобразования при впрыске топлива в цилиндр ДВС. Тезисы докладов НТК
"Луканинские чтения. Проблемы и перспективы развития автотранспортного ком­
плекса".- М.: МАДИ (ГТУ), 2003. - С. 57 - 59.
3. Чесноков С.А., Демидов М.И. Моделирование тепломассообмена и химической
кинетики образования окиси азота в ДВС с искровым зажиганием. Известия ТулГУ,
сер. Автомобильный транспорт, вып. 7 - Тула: ТулГУ, 2003. - С. 255-265.
4. Чесноков С.А., Кузьмина И.В., Демидов М.И., Рябых А.Н. Химическая кинетика
горения в двигателях ВАЗ. 3-я Всероссийская научно-техническая конференция "Со­
временные тенденции развития автомобилестроения в России". Тольятти: Тольят. гос.
ун-т, 2004. - С. 159-165.
5. Чесноков С.А., Демидов М.И.. Химический турбулентный тепломассообмен за
фронтом пламени в ДВС. Часть 1. Химическая кинетика. Известия ТулГУ, сер. Авто­
мобильный транспорт, вып. 8 - Тула: ТулГУ, 2004. - С. ПО-124.
6. Чесноков С.А., Демидов М.И. Химическая кинетика при турбулентных пульсациях
температуры в ДВС. Известия ТулГУ, сер. Автомобильный транспорт, вып. 8 - Тула:
ТулГУ, 2004.-С. 124-131.
7. Демидов М.И. Моделирование тепломассообмена и химической кинетики оксида
углерода в ДВС. Всероссийская научная конференция студентов и аспирантов «Мо­
лодые исследователи-региону». Вологда: ВоГТУ, 2004. -с.191-193.
8. Демидов М.И. Моделирование тепломассообмена и химической кинетики образо­
вания оксида азота в ДВС с искровым зажиганием. 1-я всероссийская научнотехническая конференция студентов и аспирантов «Идеи молодых - новой России».
Тула: ТулГУ, 2004. - регистрационное свидетельство Х23928, номер гос. регистрации
0320400297.
t2Z82 0
РНБ Русский фонд
2006-4
24682
Иад лиц. ЛР № 020300 от 12.02.97 . Подписано в печать
Формат бумаги 60x84'/,,. Бумага офсетная.
Усл. пен. л. >*• . Уч.-изд. л Р,5 .
Тираж ■/РСэю. Заказ5*
Тульский госуаарственный университет.
300600, г. Тула, пр. Ленина, 92.
Отпечатано в Издательстве ТулГУ
300600, г. Тула, ул. Боядина, 151.
Документ
Категория
Без категории
Просмотров
0
Размер файла
1 026 Кб
Теги
bd000103342
1/--страниц
Пожаловаться на содержимое документа