Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги / Метанол как топливо для транспортных двигателей

..pdf
Скачиваний:
14
Добавлен:
12.11.2023
Размер:
7.82 Mб
Скачать

 

 

1

 

 

Т

= 2500 К

 

 

т = 2750 К

 

T = 3000 К

 

Компо­

Т = 2250 К

 

 

 

 

 

Si-ù Wc

 

 

 

*

 

 

 

 

 

 

 

 

ненты

 

figue

а Д:

^ и с

Zi дИс

6i

Ôguc.

Нс

 

<5дас.

 

 

нДж/кмоль

%

% кДж/кмоль

%

% кДк/кмоль

%

% кДж/кмоль

%

%

 

и

16,44

1,02

-

86,405

2,03

-

337,39

3,66

-

1052,38

6,16

-

 

0

6,436

0,39

-

39,84

0,938

-

172,94

1,88

-

568,81

3,33

-

U 1

})

0,0011

-

-

0,0143 -

 

 

0,117

0,0013

-

0,673

0,0039

-

ГТ1

Н2

355,78

22,08

«

833,94

20,1

-

1754,31

19,04

 

3194,16

18,7

-

~

-

 

он

189,096

11,73

 

558,82

13,15

-

1323,16

14,36

2614,34

15,3

-

 

с о

1043,78

64,77 -

2710,2

63,78 -

5628,02

61,07

-

9654,08

56,51

-

00%

1611,53

100

1,96

4249,22

100

5,Î6

9215,94

100

11,2

17084,44

100

20,76

1399,56

86,85

 

3564,14

83,88 4,33 7382,33

'80,11I

8,97

12848,24

75,2

16,61

СО%*он1588,66

98,58

1,93

4222,96

-97,03 5,01 8705,49

94,46

10,58

15463,58

90,51

18,79

 

2,4.

Расчет процесса

сгорания

и образования Л/О

 

 

 

 

в цилиндре двигателя при работе на метаноле

 

 

 

 

Разработанная методика основьтается на двухзонной модели

[17] и позволяет производить расчет давления

и

температур

(све­

жей

смеси и продуктов

сгорания),

индикаторных

показателей

цикла,

а также текущей

концентрации

/0Л

по

камере

сгорания в

функции

угла

поворота кривошипа. С ее помощью

можно

оценить влияние

на

индикаторные

показатели и выход

л Ю

таких

параметров,

как нача­

ло и

продолжительность

сгорания,

состав смеси,

потери теплоты

от

газа

в стенки,

степень

сжатия

и др.

 

 

 

 

 

 

 

В данной модели заряд цилиндра в процессе сгорания условно

разделен на две зоны: зону несгоревшей смеси и зону продуктов

до

сгорания. Принимаем, что топливовоэдушная- (несгоревшая)

смесь

и в ходе процесса сгорания является гомогенной. Зона продуктов сгорания в калщый момент цикла имеет однородный состав и темпера­ туру. Расчет процесса сгорания производится по заданному закону

выгорания топлива с учетом диссоциации продуктов

сгорания. Равно­

весные концентрации отдельных продуктов сгорания,

необходимые для

учета диссоциации

(СО и W2

) и образования

А/О

(02 ,

, О ),

определяются по аппроксимирующим полиномам (см. табл.2 .2 ).

 

Для решения поставленных задач разработана программа расче­

та процесса

сгорания и образования tJO применительно к ЭВМ

типа

ЕЕ на языке Фортран-ГУ.

 

 

 

заклю­

Методика расчета процесса сгорания и образования

чается в следующем.

топлива задаем уравнением [17]

 

Закон

сгорания

 

 

 

 

х*о,ьО-cos£$>% ,

 

 

 

где %■ >

%

1

Ч>

- угол

начала сгорания,

угол продолжитель­

ности видимого сгорания и текущее значение угла поворота криво­

шипа соответственно. Для расчета давления и температур использу­

ем систему из

четырех уравнений [17 1

 

 

 

Q = A Ü + L ;

(2.23

 

M m C p p i T * Мсмmcpm Ta ,*MwmcPfnTllc ;

 

2 .(25)

 

Тсн-ТиС— ) 1^

 

(2.26)

Количество теплоты,

сообщенное рабочему телу к данному

моменту

процесса

сгорания,

определяем по уравнению

 

 

 

Q - ( Q HBU, X - Qgüc)ypt

 

(2.27)

где В ц

- количество метанола в цилиндре за цикл;

у

- коэф­

фициент

потерь теплоты в результате теплообмена. Лотери

теплоты

Qguc учитывают

невыделйвшеееяее количество вследствие не­

полноты сгорания (при cL <. I) и диссоциации и находятся по содержанию иСО в продуктах сгорания, равновесные концен­

трации которых определяются по аппроксимирующим полиномам второй степени.

Изменение внутренней энергии определяем по уравнению

 

 

 

A Ü - U - UH

,

 

где

il

 

t/ц —

внутренние энергии

газа в цилиндре *в любой мо-

момент

и

в начале процесса сгорания

соответственно,

 

 

 

U«M m C w n T «

I V М аmи С% тм .

Здесь

М , М е и

- зардц цилиндра в процессе и в начале сгорания;

m

CVmf m C% -

средние мольные теплоемкостй заряда в процессе

сгорания и

свежей

смеси; Т

Тн

- температуры в цилиндре в

процессе и в начале сгорания.

 

 

Теплоемкость свежей смеси в цилиндре находим с учетом со­

держания в

заряде

воздуха,

остаточных

газов и паров топлива:

m c Z - 4 m c }m t гг пи$т + гттс?т .

Здесь

Zg

Sÿ.

Z7

относительное объемное

содержание воз­

духа,

остаточных

газов

и паров топлива соответственно, определя­

емое

по формулам

 

 

 

 

 

 

 

*

9

 

оСМо

 

 

 

 

 

 

(оСМа* i/mrXj tf)

 

 

 

г

ь

г

 

 

 

 

Ч

 

T T

ÿ - »

 

 

 

 

*т'

 

1

 

 

 

 

mrXeCMo+Vm,#-#)

 

т с £ ,

_ 1 Г ас

ГГ1cjm

- средние мольные

теплоемкости

т t-vm

 

воздуха, остаточных газов и паров топлива. Средние мольные тепло­

емкости воздуха и топлива

(метанола) находим из уравнений [42J:

т с £ * Î9.8

+ 209,3-« Г 5 Т

m Cjm= 21,15

+ 7 ,09-ГО"2 Г

среднюю мольную теплоемкость остаточных газов определяем по ап­

проксимирующему

полиному.

 

 

Выразим

тСут Ти

через энтальпию остаточных

газов:

m C L V

m

 

где ni GpmТн » W ?

- энтальпия

остаточных газов,

устанавли­

ваемая по аппроксимирующему .полиному в зависимости отТц , Л

Тогда внутреннюю энергию газа в начале процесса сгорания находим из выражения

U M . Мсм[Сгё т с£ ,+ггтС£,)Тн * гг Н?-ЬЫЬTj.i2.26)

Температура

Тм

и давление

рн смеси в начале видимого

сгорания

могут

быть

наЭДены по уравнению политропы сжатия,

а ко­

личество

заряда

в цилиндре к началу

сгорания - из уравнёния

сос­

тояния при известных

Тц

р н

и

VM

Количество газа в цилин­

дре в процессе сгорания определяем из

выражения

М с М

см£*1 +

 

)],

 

 

теплоемкость зардца в процессе сгорания - по формуле

 

 

т

C J « Q-x)т C/.Ï Т см + JCnnCvmC Tn.c.

 

Внутреннюю энергию в процессе сгорания U

находим аналогично*

U H с использованием выражений для энтальпий остаточных газов

и продуктов

сгорания

ИТ

 

Окончательное выражение для

а

имеет вид

 

 

 

 

 

 

 

 

u* MÔ-X)f a m

 

ZftncJmÿr*

 

 

 

)

Значение

Н т

 

 

устанавливаем по аппроксимирующему уравнению в

зависимости от

 

T

,ci.

, Р

к данному моменту цикла,

опреде­

 

Работу,

совершаемую

газом

ляем по

уравнению

 

 

 

 

 

 

 

 

 

 

 

 

L < +

Q5(p^+àpj A V Ii

 

 

 

где

p.|

, p

 

 

- давление газа в начале

и в конце расчетного ин­

тервала;

дV

 

- изменение объема в течение рассматриваемого ин­

тервала,

A V B V-V(|

р н ,

 

в начале каждого

шага рас­

 

Параметры

газа

Т4

чета известны.

 

Давление в конце

расчетного

интервала р

может

быть

определено

по уравнению состояния. Тогда

 

Подстановка значений Q , U H il, и L , нашейных по

уравнениям (2.27) - (2.30), в уравнение (2.23) позволяет привести

его,

после

соответствующих

преобразований,

к

квадратному уравне­

нию относительно средней температуры газа

в цилиндре,

которая оп­

ределяется затем

обычным способом.

текущее значение в

 

После' вычисления Т

можно определить

цилиндре по уравнению (2.24).и температуру смеси по уравнению

(2.26). Температуру в зоне

продуктов сгорания

Т П|С|

получаем,

используя

исходное уравнение

(2.25):

 

 

 

 

ЛП Срт Tptс* -^^Мя1Срт Т “Мсм Ш .Срт Тсм).

 

Заменив значение*

р 0«С м

 

 

 

 

11 П-С

ГПLpm |п.с

энтальпией продуктов сгорания Н т

определяемой по аппроксимирующему полиному

в функции

от Т

аС

р

(для диапазона

температур 1700. ..3300 К), получим:

Н “ - М т С р Д * M « m с £ Т с « - о .

МПС

таппроксимирующий полином второго порядка,

после соответствующих преобразований получим квадратное уравнение

относительно Гас , которое решается обычным способом.

Индикаторные показатели цикла вычисляем, используя значение индикаторной работы цикла, равное

где L Q .

Lp

Lc>#c ~ работы сгорания,

расширения и

сжатия

соответственно.

 

при сгорании

исполь­

Для расчета процесса образования hlQ

зуем значения температур в зоне продуктов сгорания и принимаем,

что образование tJO происходит по двум реакциям цепного и одной

реакции бимолекулярного механизма, причем определяющим является

цепной механизм, расчет образования N 0 выполнялся по уравнению

[17J

с/г,JO

.

Р - г л .

/ 1

JZ.HQ

*\ / K m

<£а

fc,f0

ZOL\

d\f>

 

■?49,4-пТп.с_\

 

 

 

 

2

1

где 2 ^

 

 

 

<2но

-

концентрации

компонентов

про-

дуктов сгорания,

определяемые по условиям равновесия;

К*п

,

К 4

К^п

- константы

скоростей реакций.

Расчет проводится

методом

 

численного интегрирования Эйлера - Коши. Концентрация

в зоне

продуктов

сгорания

на

конец расчетного участка

с ите­

рациями на каждом шаге

 

 

 

 

 

 

 

 

 

 

 

o)i

+

 

 

 

 

 

Проверка адекватности двухзонной модели. Двухзонная модель процесса сгорания и образования оксидов азота разработана для

двигателя с

электрическим

зажиганием как при работе на бензине

[ 17 ]» так

и на метаноле.

Проверка ее адекватности опытным дан­

ным производилась сравнением диаграмм давления газа в цилиндре,

кинетики образования оксидов азота и их интегральных выбросов. В

качестве опытной диаграммы давления использовалась средняя диаг­ рамма, полученная по 5Ï2 последовательным циклам. Образование ок­

сидов азота в камере сгорания в функции угла поворота кривошипа исследовалось экспериментально с помощью стробоскопического отбо­ ра проб газа газоотборными клапанами (17].

При использовании в качестве топлива бензина сопоставление опытных и расчетных результатов показало хорошую адекватность мо­ дели. Дополнительно производилась проверка адекватности модели

для случая использования метанола по

интегральным выбросам-^0 в

зависимости от коэффициента избытка

воздуха при сгорании (рис.

2 .6), подтвердившая хорошую адекватность модели и в этом случае, lia рис. 2.7 приведена часть ин^рмации по процессу сгорания с помощью двухзснной модели.. Кроме указанных данных на печать вы­ водится таг-* информация по средней температуре газа в цилиндре, текущим и равновесным концентрациям tiO в зоне продуктов сгора­

ния, концентрации-молекулярного и атомарного кислорода в зонё продуктов сгорания.

Рис. 2.6. Зависимость концентрации IslO от состава

смеси: П = 2500 мин-1; % =3 5 0 ° =

= 45°; £ = 9 , 5 н Э — опытные точки для

метанола

Рис. 2.7. Сравнение расчетных концентраций

М О и

показателей процесса сгорания для метанола

и бензина: Л =

2500

мин-*; t

9,5;

ifo =

350°;

=

45°; _ _ _ _-

метанол;

__ _ _ _

.- бензин

 

 

 

3. Способы применения метанола в *ДВС

ЗЛ . Применение метанола в качестве добавки

кбензинам

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

лировок системы питания, агрессивностью смесей к некоторым мате­

риалам системы питания.

Фазовая стабильность и подбор стабилизаторов. Учитывая вы­ сокую гигроскопичность метанола, при использовании ВМС необходи­ мо решить проблему низкотемпературной ее стабильности, которая оценивается по температуре расслоения согласно ГОСТ 5066-56. Фа­ зовая стабильность BùC исследовалась на бензинах различного угле­ водородного состава с применением в качестве стабилизаторов выс­

ших

спиртов (пропанола, иэобутанола), диметил- и 3,5-метилпиразо-

ла,

тетрасидрофур^урилового спирта

и др.

 

В результате предварительных

экспериментов установлено, что

изобутанол обладает наилучшим депрессаторным эффектом. Из углево­ дородов, входящих в состав бензинов, положительное влияние на

стабильность оказывают ароматические. Исследовались в основном

смеси с содержанием 10.. Л 5 %метанола и 3.>Л2% иэобутанола по массе на бензинах, содержащих 12...40%по массе, ароматических

углеводородов.

Установлено, что влияние иэобутанола на стабильность смесей более значительно, чем ароматических углеводородов. Например, при

увеличении содержания изобутанола в два раза (с 4 до 8%) темпе­

ратура расслоения снижается на 30 К. Увеличение содержания арома­

тических углеводородов в бензине с 12 до 24 %снижает ее на 20 К«

В результате обработки экспериментальных данных была получена за­ висимость Тр от отношения метанол/изобутанол, представленная на рис. 3.1. Анализ показывает, что превалирующее влияние на ста­

бильность

оказывает отношение спиртов метанол/изобутанол м / С ^

в добавке,

а не количество самой добавки. Причем влияние

Рис. 3.Î. Зависимость температуры расслоения

от соотно­

шения спиртов Q/c^

в ВМС. Лр =

12%:

 

• —

- - 1 0 %

=

 

; —

о —

-

12,5;

 

о--- 15. Ар

18 % :—

1- - - -

-

10 %

 

А--- — 12,5;

 

 

■ А—

15• Ар

= 30 % —

-- о--- ; Ар =

40 % --- *--

 

 

 

и

Рис. 3.2. Влияние

сроков

хранения на стабильность

фазовый

состав

ВМС

(^i/c

=

1,67)

при

опыт­

но-эксплуатационных

условиях

 

 

 

 

неодинаково. При

»

< 2 , 5

происходит

наиболее интенсивный

рост стабильности,

а при

 

 

3

кривые .стабильности прохо­

дят полого. Характер влияния спиртовой добавки на стабильность ВМС для базового бензина с различным содержанием ароматики прак­ тически одинаков.

С учетом высокой гигроскопичности метанола необходимо бы­ ло оценить влияние сроков хранения ВМС на стабильность. Исследо­ вание проводилось в лабораторных и опытно-эксплуатационных усло­ виях. Во втором случае смесь, состоящая ($ по массе) из бензина - 76; метанола - 15 и иэобутанола - 9, хранилась в течение трех месяцев (с марта по май) в’стандартной прицеп-цистерне ПЦ-5,6-817 для перевозки нефтепродуктов.

Анализ полученных данных (рис. 3.2) показывает, что в опыт-

Соседние файлы в папке книги