CC BY
159
ТЕПЛОВЫЕ ДВИГАТЕЛИ
DOI: 10.18698/0236-3941-2015-6-59-73 УДК 621.434:536
ВЛИЯНИЕ ФОРМЫ ВПУСКНЫХ КАНАЛОВ НА ЭФФЕКТИВНЫЕ И ЭКОЛОГИЧЕСКИЕ ПОКАЗАТЕЛИ СРЕДНЕОБОРОТНОГО ДИЗЕЛЯ
Р.З. Кавтарадзе, А.А. Зеленцов
МГТУ им. Н.Э. Баумана, Москва, Российская Федерация e-mail: kavtaradzerz@mail.ru; zelentsov.aa@gmail.com
Проведено численное исследование организации вихревого движения во впускной системе среднеоборотного дизеля, а также процессов смесеобразования и сгорания топливовоздушной смеси в камере сгорания. Математическая модель основана на осредненных по Рейнольдсу уравнениях переноса количества движения, энергии, массы и концентрации реагирующих веществ, которые дополняются (k — Z — f )-моделью турбулентности и моделью сгорания (CFM-моделью когерентного пламени) для определения скорости химической реакции окисления топлива. Задача решена в два этапа: сначала разработана форма впускных каналов, обеспечивающая интенсификацию вихревого движения потока в горизонтальной плоскости цилиндра, после чего рассчитаны эффективные и экологические показатели двигателя с модифицированной конструкцией впускных каналов в крышке.
Ключевые слова: среднеоборотные дизели, математическое моделирование, вихревое число, эффективные показатели.
INFLUENCE OF INTAKE PORT SHAPE ON EFFECTIVE
AND ECOLOGICAL PARAMETERS OF MIDDLE-SPEED DIESEL ENGINE
R.Z. Kavtaradze, A.A. Zelentsov
Bauman Moscow State Technical University, Moscow, Russian Federation e-mail: kavtaradzerz@mail.ru; zelentsov.aa@gmail.com
The paper presents the numerical analysis of swirl motion generation inside the intake system of a middle-speed diesel engine. It also analyzes the processes of both mixture formation and air-fuel mixture combustion in the combustion chamber. The mathematical model is based on the three-dimensional non-stationary equations of momentum, energy, motion, mass transfer, and consistency of reacting substances in the Reynolds form with the addition of k—Z—f turbulence model and the Coherent Flame Model (CFM), which are used for estimating the combustion chemical reaction rate. The authors solved the problem as follows: at first, an intake port shape with a higher swirl number was developed, and then effective and ecological parameters of the diesel engine with a modified shape of the intake ports in the cylinder head were calculated.
Keywords: middle-speed diesel engine, mathematical modeling, swirl number, effective parameters.
Поршневые двигатели внутреннего сгорания остаются на сегодняшний день основой силовых установок для наземного транспорта. Обеспечение высоких эффективных и экологических показателей
поршневого двигателя неразрывно связано с совершенствованием процессов в камере сгорания. Один из способов улучшения процессов смесеобразования и сгорания топливовоздушной смеси заключается в интенсификации вихревого движения свежего заряда в цилиндре двигателя.
Наиболее распространенным методом оценки интенсивности вихревого движения заряда в цилиндре дизеля является метод статической и изотермической продувки на холодных физических моделях [1]. Еще в 20-е годы прошлого века немецкий исследователь Ф. Засс из фирмы AEG устанавливал пластину с тензо- и индуктивным датчиками в камеру сгорания дизеля и исследовал воздействие потока впускного воздуха на пластину. Таким путем он определял энергетические возможности потока в цилиндре. Наличие пластины в цилиндре, конечно, деформировало поток, и поэтому Г. Рикардо, основатель фирмы Ricardo, впервые заменил пластину крыльчатым анемометром. Регистрируя частоту вращения крыльчатки, расположенной в цилиндре, он ввел понятие "вихревое отношение" (вихревое число или число Рикардо), представляющее собой отношение частот вращения n в крыльчатке и n коленчатого вала двигателя. Ясно, что частота вращения крыльчатки в таком случае отождествляется с частотой вращения воздушного вихря в плоскости, перпендикулярной оси цилиндра. В настоящее время используются и отличные от упомянутого определения вихревого числа, такие как отношение тангенциальной (окружной) Cu и аксиальной (осевой) Ca скоростей потока Dc = Cu/Ca, а также вих-MRp
ревое число DM = ——, определяемое через измеренные значения
m 2
крутящего момента М, созданного вращающимся потоком, и расхода воздуха m через цилиндр (R — радиус цилиндра, p — плотность воздуха при впуске). В [2] показано, что между указанными показателями интенсивности вихревого движения в цилиндре поршневого двигателя существует однозначная зависимость.
Целью работы является повышение эффективных показателей среднеоборотного дизеля с обеспечением заданных экологических показателей путем модернизации конструкции впускных каналов, генерирующих вихревое движение в цилиндре в процессе впуска. Для решения поставленной задачи определяется форма впускных каналов, обеспечивающая интенсификацию закрутки потока в горизонтальной плоскости (в плоскости, перпендикулярной оси цилиндра), после чего рассчитываются эффективные и экологические показатели двигателя с модифицированной конструкцией крышки с предложенными формами впускных каналов.
Математическая модель. Математическая модель нестационарных процессов переноса количества движения, энергии, массы и концентрации реагирующих веществ в расчетном объеме достаточно подробно описана в [1, 3, 4]. В табл. 1 эти уравнения представлены как
Таблица 1
Уравнения законов сохранения в исходной форме и в форме Рейнольдса
Уравнения законов сохранения как частные случаи обобщенного дифференциального уравнения
дх
дх,-
дх,-
dxjJ
Уравнение сохранения количества движения (уравнение Навье Стокса)
Ф II]. Г; LI.
dp дх.
Sa, =G,~ ^r- +11,
Исходная форма
Форма Рейнольдса
Новые неизвестные
DU] ^ dp д р-- = Gi —— + —
Dt дх, дх.
ГII] д1Г. 2, dWk M-l —L + —~—8,—-
дх,- дх; 3 ' дх,.
_DWj _ — dp | д Dt ' дх , дх,-
dir, dir, 2 dirk . _„.,„., H —+ - rll/ll
дх,-
дх.
дх,.
р/Г/ II ". -тензор рейнольдсовых турбулентных
напряжений, определенный по пульсационным
составляющим скорости;
(г/г, д¥, 2 я Г/Г; ^ т„ = ц -+—;---8,у- -тензор вязких
■' [ дх1 дх, 3 ■' дхк )
(турбулентных) напряжений, определенный по усредненным значениям компонентов скорости
Уравнение сохранения энергии
Ф = Н, Та, = —
Исходная форма
DH др д
р— ал . +—+—
Dt ■' ■' дт дх,
дх,-
Форма Рейнольдса
Новая неизвестная
• (1.11'. +w,.Or +■
сррТ' И'1 - турбулентный перенос энтальпии рсрТ' посредством флуктуации скорости \Г]
дх,-
Уравнение сохранения массы (неразрывности)
Ф = 1, Гл =0, =0
Исходная форма
дх дх,- ■"
Форма Рейнольдса
дх дх,- ^ " 3
i = 0
Уравнение диффузии (концентрации)
Ф = —. гф =А, р
= т
Исходная форма
DC 3 ЗС
-=- D-I + т
Dx дх,- дх,-
Форма Рейнольдса
Новая неизвестная
DC д \ дС . ■ -=- D--( II |+ т
Dt дх, дх.
С'1Г- - турбулентный диффузионный перенос массы компонента посредством флуктуации скорости /Г_(С - концентрация компонента)_
Примечание. В уравнениях применяется правило суммирования по индексу (У, у, к = 1, 2, 3), повторяющемуся дважды, а усреднения выполняются по Фавру, т.е. с использованием плотности в качестве весовой функции.
частные случаи обобщенного дифференциального уравнения, там же приведена тензорная форма записи этого уравнения для декартовой системы координат.
В уравнениях, приведенных в табл. 1, приняты следующие обозначения: Б/(Бт) — субстанциональная производная; р — плотность; р — давление; С — проекция вектора плотности объемных сил на ось Охг прямоугольной декартовой системы координат; С — концентрация;
Н — полная удельная энергия; У^ = ^р grad — член, вы-
ражающий объемную деформацию; р — динамическая вязкость; ср — теплоемкость при постоянном давлении; ыг — скорость химической реакции на единицу объема; — количество выделяемой теплоты на единицу массы; Л — теплопроводность; 8^ — символ Кронекера; Б — коэффициент диффузии; т — интенсивность источника массы (скорость изменения массы химической компоненты в единице объема); IV — вектор скорости; Ф — произвольная зависимая переменная; Гф — коэффициент обмена (диффузии); Sф - источниковый член, который в общем случае можно представить как разность генерации Sфg и аннигиляции Sфa потоков, т.е. Sф = Sфg — Sфa. Конкретный вид Гф и Sф, а также Sфg и Sфa зависит от смысла переменной Ф (см. табл. 1) и после подстановки соответствующих значений из обобщенного дифференциального уравнения получаем уравнения Навье - Стокса, Фурье-Кирхгофа, Фика и сохранения массы (неразрывности).
Решение предусматривает определение локальных параметров газа во всей расчетной области, которые представляются в виде суммы осредненной и пульсационной составляющих. Незамкнутая система уравнений в форме Рейнольдса дополняется (к—(—/)-моделью турбулентности [4-6], где к — кинетическая энергия турбулентности, £ — нормированный масштаб скорости, а / — эллиптическая функция релаксации.
Уравнения (к—(—/)-модели имеют вид [4, 5]:
Dk ч д
pDT = Р (Рк — e) + dXj
р +
Pt Ок
дк
dxi
De c*£lPk — c£2€ д p— = p f1 ^-+
Dt
дх^
, Pt р +--
О f
д€ дх.
д
РТТг = P/ — PkPk + дх
, Pt P +--
К
oz) дхj
(1)
Эллиптическая функция релаксации / определяется из выражения
/ — I
2
дxj дxj
I + Рк
= | Ci + С2 —
2/3 — Z
(2)
где Т — масштаб времени турбулентности; I — масштаб длины, при этом _ _
'к а \ /ТУ3\1/2
T = max
min | —, - I , cT I - I
1 ^ vac« |S| e r Ч е )
l = ci max
k3/2\ W 3/4 — ,Cn
(3)
(4)
Значения параметров модели в (3) и (4) ст = 0,36 и сп = 85 определены эмпирическим путем; Б = \j2SijБ^, где Б^ — среднее значение деформации элементарного объема жидкости в единицу времени.
В выражении для скорости диссипации кинетической энергии турбулентности е была изменена одна из констант:
& = с£1 (1 + 0,045^Г7с). (5)
Данное изменение имело целью более точно отразить влияние пристеночных эффектов в предлагаемой модели. Остальные константы, входящие в выражение (1), аналогичны константам стандартной (к-е-/)-модели: се2 = 1,92; = 1; о£ = 1,3.
Влияние химического процесса горения учитывается с помощью интенсивности внутреннего источника теплоты (¡,и и интенсивности источника массы т. Значения этих параметров могут быть вычислены с помощью скорости ыг химической реакции сгорания
= Qr ыг; т = —ыг, (6)
где Qr — выделившееся в результате химической реакции количество теплоты на единицу массы. Скорость процесса сгорания вычисляется по модели когерентного пламени, основанной на концепции распространения ламинарного пламени [1, 4]. Согласно этой концепции, усредненная по всему фронту пламени скорость и толщина фронта 81 зависят только от давления, температуры и состава свежего заряда. При этом считается, что реакция начинается в относительно тонких слоях, которые отделяют свежий несгоревший газ от продуктов сгорания.
В модель вводится дополнительное дифференциальное уравнение переноса плотности фронта пламени, записываемое относительно Е — площади фронта пламени на единицу объема. Данное уравнение имеет вид [1, 4]
«Е + А № Е) = ® (^ п) + (7)
дт дх¡- ' дх3 \ Ргш дх3)
где V — турбулентная кинематическая вязкость; Ргш — турбулентное диффузионное число Прандтля (число Шмидта); Бе — источниковый член, равный разности генерации протяженности поверхности фронта
пламени в результате деформирования из-за турбулентного горения и аннигиляции поверхности фронта пламени в результате расхода реагентов.
Тогда средняя скорость реакции горения топлива определяется как
Wr = -pTp wi £, (8)
где pTp — парциальная плотность топлива; wl — скорость ламинарного пламени, зависящая от локального давления p, температуры T свежего заряда газа и локального коэффициента избытка воздуха ав.
Расчеты проводились с помощью SD-CRFD-кода FIRE, разработанного фирмой AVL List GmbH (Австрия) [4]. Ядро FIRE основано на численном методе контрольных объемов с использованием усовершенствованного алгоритма SIMPLE. Верификация математической модели осуществлялась методом сравнения экспериментальных и расчетных индикаторных диаграмм, позволяющим в случае необходимости уточнить эмпирические коэффициенты в модели сгорания [7-9].
Результаты численного моделирования закрутки потока на впуске. Объект исследования — четырехтактный среднеоборотный дизель (одноцилиндровый отсек) с наддувом, размерностью d/s = 265/310 мм/мм и степенью сжатия 15,5. Мощность двигателя составляет ne = 331 кВт при частоте вращения коленчатого вала n = 1000 мин-1.
Для определения закрутки воздуха в цилиндре двигателя, обеспечиваемой впускными каналами, на основании трехмерной модели крышки цилиндра строится внутренний объем каналов. Полученная модель впускных каналов дополняется объемом цилиндра экспериментальной установки (диаметр соответствует диаметру цилиндра двигателя D, длина равна 2,5D (рис. 1, а), а также некоторым объемом на впуске для моделирования граничных условий на входе в канал (рис.2,а,б). При расчетах крыльчатка анемометра заменяется выделенной областью цилиндра (рис. 1, б), в которой определяются осевые и аксиальные составляющие скорости потока воздуха.
Рис. 1. Расчетные схемы:
а — экспериментальная установка для определения вихревого числа; б — крыльчатка анемометра (Б — диаметр цилиндра)
а б
Рис. 2. Представление расчетной области (базовые впускные каналы) в виде:
а — твердотельной модели; б — разбиения на контрольные объемы
Результат разбиения на контрольные объемы представлен на рис. 2, б. Максимальный размер элементов составляет 4 мм с измельчением сетки на огневом днище крышки цилиндра до 2 мм и седлах клапанов до 1 мм (число элементов составило 1 417 811). В процессе расчета перемещения клапанов не происходит (расчет стационарный).
Анализ трехмерных результатов показал, что в целом базовая геометрия впускных каналов обеспечивает безотрывное течение воздуха в системе впуска. Несмотря на то что каналы в крышке рассматриваемого двигателя — тангенциальные (обеспечивающие формирование вихревого движения в объеме цилиндра, а не в объеме канала), дальний от впускного сечения канал отчасти выполняет роль вихревого, поскольку в нем осуществляется интенсивное изменение направления движения воздуха (рис. 3). При этом скорость истечения воздуха в клапанной щели этого канала оказывается наибольшей и составляет ~ 40 м/с.
Поскольку на первом этапе исследования течения во впускной системе рассматривается как стационарный, процедура расчета является
Рис.3. Расчетные значения скоростей воздуха (а) и кинетической энергии турбулентности (б) для базовых впускных каналов
итеративной. Изменение вихревого числа в горизонтальной плоскости (перпендикулярно оси цилиндра) Вп, равного отношению частот вращения воздушного вихря и коленчатого вала двигателя на заданном режиме работы, в зависимости от номера итерации приведено на рис. 4.
Следует отметить, что отрицательное значение вихревого числа показывает, что вращение происходит в обратном направлении, принятом в качестве положительного при определении Бп. Из результатов расчетов следует, что базовая форма каналов обеспечивает умеренную закрутку свежего заряда в цилиндре (Бп ^ 1). Оптимальная интенсивность вихревого движения заряда в дизелях определяется прежде всего в зависимости от давления впрыскивания и числа факелов впрыскиваемого топлива. При этом необходимо избежать нежелательного эффекта чрезмерной закрутки потока, приводящего к ухудшению процессов смесеобразования и сгорания [1, 10]. Аналогичное ограничение следует иметь в виду и в случае бензиновых двигателей с непосредственным впрыскиванием, в которых, в отличие от двигателей с внешним смесеобразованием [8, 9, 11], также существует ограничение по интенсивности вихря. В рассматриваемом случае для среднеоборотного дизеля было установлено, что интенсификация вихревого движения может привести к существенному росту эффективных показателей двигателя, что требует обеспечения высокого (по сравнению с базовым двигателем) уровня закрутки потока на впуске.
Повышение вихревого числа может быть обеспечено различными вариантами модификации существующих каналов, часть из которых может быть осуществлена с минимальными изменениями исходной геометрии крышки цилиндра, что актуально с точки зрения снижения затрат на ее изготовление, часть реализуется при значительном изменении конструкции (например, замена дальнего от впускного сечения канала на вихревой). Однако второй вариант, хотя и улучшает условия для поступления воздуха в цилиндр двигателя, но конструктивно оказывается гораздо сложнее. Кроме того, сужение сечения канала в горизонтальном направлении при его повороте вокруг оси клапана приводит к необходимости увеличения размера канала по вертикальной оси (для сохранения проходного сечения), а следовательно, габаритных размеров крышки в целом. К тому же растут потери на поворот потока при впуске.
Рис.4. Зависимость вихревого числа в горизонтальной плоскости Оп от номера итерации (итерация № 3500 соответствует окончанию расчета)
Рис. 5. Зависимости вихревого числа в горизонтальной плоскости (Оп) от формы впускных каналов:
а — угол наклона оси канала к плоскости огневого днища крышки 100, б — угол наклона 170
При сохранении конфигурации и проходных сечений базовых каналов изменяется угол наклона оси канала (до их разделения) относительно горизонтальной плоскости (угол наклона оси базовых каналов составляет 13°). При этом рассматриваются варианты с углами наклона 10о и 17° (см. далее).
Уменьшение угла наклона оси канала относительно горизонтальной плоскости привело к росту максимальной скорости истечения газа в клапанном канале (с 40 м/с при угле 13° до 53 м/с при 10°). Данное обстоятельство оказывает положительное влияние на создание вихревого потока в цилиндре двигателя (рис. 5, а): происходит интенсификация вихревых потоков в горизонтальной плоскости (число Бп растет от 0,99 до 1,28). При большем угле наклона оси канала к плоскости огневого днища крышки (рис. 5, б) отмечается снижение максимальных скоростей потока газа (до 38 м/с в клапанной щели). Данное обстоятельство вкупе со стремлением потока на входе в цилиндр к осевому объясняет уменьшение интенсивности горизонтального вихря Бп с одновременным увеличением вертикального вихря Тп.
Обобщенные результаты трехмерных расчетов течения газа во впускной системе в зависимости от формы каналов в крышке рассматриваемого двигателя следующие: при углах наклона оси каналов к плоскости огневого днища крышки 10°, 13° и 17° коэффициенты расхода составляют 0,569, 0,563 и 0,596, а вихревые числа — 1,28, 0,99 и 0,78 соответственно.
Таким образом, для повышения интенсивности вихревого движения заряда в цилиндре исследуемого двигателя рекомендуется использование каналов с уменьшенным углом наклона оси относительно горизонтальной плоскости (совпадающей с плоскостью огневого днища крышки цилиндра), как обеспечивающих рост в 1,29 раза вихревого числа Бп.
Анализ влияния интенсивности закрутки потока на впуске на эффективные и экологические показатели дизеля. Второй, заключительный, этап расчета предусматривал исследование нестационарных турбулентных процессов переноса, горения и образования вредных веществ в цилиндре исследуемого дизеля в целях оценки эффективных и экологических показателей двигателя при использовании модифицированных впускных каналов.
Моделирование рабочего процесса двигателя было проведено в трехмерной постановке, что позволило учесть сложные течения внутри камеры сгорания в поршне (перетекание рабочего тела с периферии при движении поршня вблизи ВМТ и т.п.), а также особенности внутреннего смесеобразования (геометрические параметры распылителя, место его установки и ориентация относительно оси цилиндра). Расчетный объем представлен сектором камеры сгорания с подвижной нижней границей (огневым днищем поршня) и неподвижной верхней (головкой цилиндра).
Разбиение на контрольные объемы (КО) показано на рис. 6. При этом использовались следующие параметры разбиения: максимальный размер элементов 2,5 мм, число слоев вблизи стенок — 2, высота ячеек для описания пограничного слоя — 0,2 мм. Итоговое число элементов при разбиении меняется в зависимости от положения поршня, максимальное число достигается при положении поршня в НМТ — 39 664 контрольных объемов (для одного сектора камеры сгорания, число секторов соответствует числу сопловых отверстий в форсунке, для исследуемого двигателя их 8). Таким образом, максимальное число контрольных объемов для всего объема камеры сгорания (см. рис.6) составляет 317 312.
Контроль достоверности результатов численного эксперимента является одним из основных этапов исследования рабочего процесса
Рис. 6. Разбиение сектора камеры сгорания двигателя Д500 на контрольные объемы (положение поршня в ВМТ, степень сжатия е = 15,5)
в цилиндре поршневого двигателя. В общем случае верификация проводится путем сравнения расчетных результатов с экспериментальными данными, в случае отсутствия таковых допустимо использование данных, полученных из нульмерного моделирования с применением апробированных математических моделей.
Верификация математической модели осуществлялась на основании экспериментальных данных. При верификации проводилось сравнение индикаторных диаграмм и зависимостей осредненных по объему камеры сгорания температур от угла п.к.в. (рис. 7, 8), показавшее удовлетворительное согласование расчетных данных, особенно зависимостей давления от угла п.к.в. (в этом случае погрешность по значению максимального давления рг не превышает 1,6%, а по углу его достижения — менее 1 %). При этом наибольшие значения температур различаются менее чем на 3 % (см. рис. 8).
Сравнение индикаторных диаграмм исследуемого двигателя при различных значениях закрутки заряда на впуске приведено на рис. 9. Некоторое повышение максимальных давлений цикла происходит с ростом интенсивности вихревого движения заряда в цилиндре, что
570 610 650 690 730 770 810 850 Угол п.к.в., град
Рис.7. Сравнение расчетной и экспериментальной индикаторных диаграмм исследуемого двигателя (степень сжатия в = 15,5; 720° п.к.в. соответствуют ВМТ)
2000 а 1600 й 1200 I 800
с
| 400 н
ol-1-1-1-1-1-1-
570 610 650 690 730 770 810 850 Угол п.к.в., град
Рис. 8. Сравнение осредненных по объему КС температур рабочего тела исследуемого двигателя (степень сжатия в = 15,5; 720° п.к.в. соответствуют ВМТ)
— Расчет
Рис.9. Зависимости давления в цилиндре двигателя Д500 от угла поворота коленчатого вала при различных значениях вихревого числа Оп
Рис. 10. Зависимости массовых долей оксидов азота от угла п.к.в. для различных значений Оп
объясняется улучшением условий для смесеобразования и сгорания топливовоздушной смеси.
На рис. 10 также приведено сравнение массовых долей выделившихся в результате сгорания топлива вредных веществ (оксидов азота). Зависимость для КОх оказывается аналогичной — рост Бп ведет также к некоторому увеличению выбросов оксидов азота, однако при этом возрастают и показатели эффективности работы двигателя (табл. 2).
Таблица 2
Сравнение эффективных показателей дизеля при различных значениях
вихревого числа
Параметр Значения вихревого числа Бп
0,78 0,99 1,28
Ne, кВт 320,19 325,64 328,97
ge, г/(КВт-ч) 221,1 217,4 215,2
0,42 0,43 0,436
Следует отметить, что некоторая погрешность в определение индикаторных и эффективных показателей двигателя вносится принятыми при расчетах допущениями о том, что рассматривается только сжатие
и рабочий ход двигателя, впуск и выпуск учитываются только при задании граничных условий расчета (в том числе и при задании закрутки потока на впуске). Таким образом, расчет ведется для углов п.к.в. от 570° до 846° (здесь угол п.к.в. 720° соответствует ВМТ).
Заключение. Разработаны конструкции впускной системы дизеля, позволяющие добиться повышенной (Юп = 1,28) и пониженной (Юп = 0,78) относительно базовых каналов (Юп = 0,99) интенсивности вихревого движения в цилиндре дизеля. Вихревые числа для исследуемых каналов были определены моделированием процесса впуска в трехмерной постановке с имитацией условий испытательного стенда с крыльчатым анемометром для определения расходных характеристик канала при максимальном подъеме клапанов.
Для повышения интенсивности вихревого движения заряда в цилиндре исследуемого двигателя рекомендуется использование каналов с уменьшенным углом наклона оси относительно горизонтальной плоскости (совпадающей с плоскостью огневого днища крышки цилиндра), как обеспечивающих рост вихревого числа Бп в 1,29 раза.
Трехмерный расчет с учетом локальных нестационарных температур рабочего тела, а также сложной формы камеры сгорания позволяет оценить влияние закрутки воздуха на впуске в двигатель на его эффективные и экологические показатели. В результате расчетов было получено повышение эффективных показателей двигателя при интенсификации вихревого движения воздушного потока. В то же время увеличение вихревого числа Бп, характеризующего закрутку заряда в горизонтальной плоскости цилиндра, приводит к некоторому повышению эмиссии оксидов азота. При уменьшении вихревого числа наблюдается обратная зависимость.
Таким образом, изменение геометрии впускных каналов позволяет добиться улучшения либо эффективных, либо экологических показателей двигателя. Одновременного снижения и сажи, и оксидов азота с приемлемым уровнем эффективных показателей также можно добиться, изменяя тип рабочего процесса (частично-гомогенное сгорание топлива), либо используя систему пост-обработки отработавших газов.
Работа выполнена при финансовой поддержке Министерства промышленности и торговли РФ (государственный контракт № 13411.1370399.20.011 от 18.06.2013, шифр "Выброс").
ЛИТЕРАТУРА
1. Кавтарадзе Р.З. Теория поршневых двигателей. Специальные главы. М.: Изд-во МГТУ им. Н.Э. Баумана, 2008. 719 с.
2. Вошни Г., Цайлингер К., Кавтарадзе Р.З. Вихревое движение воздуха в быстроходном дизеле с четырьмя клапанами на цилиндр // Вестник МГТУ им. Н.Э. Баумана. Сер. Машиностроение. 1997. № 1. С. 74-84.
3. Basshuysen R., Schäfer F. (Hrsg.) Handbuch Verbrennungsmotor. 4. Auflage. Vieweg und Sohn Verlag. Wiesbaden, 2007. 1032 s.
4. FIRE. Users Manual Version 2013. AVL List GmbH Graz, Austria, 2014 (License Agreement for Use of the Simulation Software AVL FIRE between Bauman Moscow State Technical Univ. and AVL List GmbH, 2014).
5. Hanjalic K., Popovac M., Hadziabdic M. A Robust Near-Wall Elliptic-Relaxation Eddy-Viscosity Turbulence Model for CFD // Int. J. Heat Fluid Flow. 2004. No. 25. P. 897-901.
6. Advanced Turbulent Heat Transfer Modeling for IC-Engine Applications Using AVL FIRE / R. Tatschl [et al.] // International Multidimensional Engine Modeling User's Group Meeting. Detroit, 2006. P. 1-10.
7. Kavtaradze R.Z., Onishchenko D.O., Zelentsov A.A., Sergeev S.S. The influence of rotational charge motion intensity on nitric oxide formation in gas-engine cylinder // International J. of Heat and Mass Transfer. 2009. Т. 52. No. 19-20. Р. 4308-4316.
8. Кавтарадзе Р.З., Онищенко Д.О., Зеленцов A.A., Финкельберг Л.А., Костюченков А.Н.Моделирование процессов в системе "впускной коллектор-цилиндр" авиационного поршневого двигателя с распределенным впрыскиванием топлива // Вестник МГТУ им. Н.Э. Баумана. Сер. Машиностроение. 2012. № 4. С. 3-16.
9. Ланшин А.И., Финкельберг Л.А., Костюченков А.Н., Зеленцов A.A., Баканов М.А. Исследование влияния предварительной закрутки потока на характеристики авиационного поршневого двигателя // Вестник ВГТУ 2012. Т. 8. № 2. С. 96-99.
10. Кавтарадзе Р.З., Сергеев С.С. Новый альтернативный (частично гомогенный) процесс сгорания как способ снижения концентраций оксидов азота и сажи в продуктах сгорания дизеля // РАН. Теплофизика высоких температур. 2014. Т. 52. № 2. С. 294-309.
11. Kavtaradze R.Z., Zelentsov A.A., Skripnik A.A. 3D Investigation of Features of Realization of Working Processes of Aircraft Piston Engines. Proceedings International User Conference Advanced Simulation Technologies (AST), 18-20 June 2013, Austria, Graz. 2013. P. 1-9.
REFERENCES
[1] Kavtaradze R.Z. Teoriya porshnevykh dvigateley. Spetsial'nye glavy [Theory of reciprocating engines. Specific chapters]. Moscow, MGTU im. N.E. Baumana Publ., 2008. 719 p.
[2] Woschni G., Zeilinger K., Kavtaradze R.Z. Air vortex motion in a high-speed diesel with four valves on a cylinder. Vestn. Mosk. Gos. Tekh. Univ. im. N.E. Baumana, Mashinostr. [Herald of the Bauman Moscow State Tech. Univ., Mech. Eng.], 1997, no. 1, pp. 74-84 (in Russ.).
[3] Basshuysen R., Schafer F. (Hrsg.) Handbuch Verbrennungsmotor. 4. Auflage. Vieweg und Sohn Verlag. Wiesbaden, 2007. 1032 s.
[4] FIRE. Users Manual Version 2013. AVL List GmbH Graz, Austria, 2014 (License Agreement for Use of the Simulation Software AVL FIRE between Bauman Moscow State Technical Univ. and AVL List GmbH, 2014).
[5] Hanjalic K., Popovac M., Hadziabdic M. A Robust Near-Wall Elliptic-Relaxation Eddy-Viscosity Turbulence Model for CFD. Int. J. Heat Fluid Flow, 2004, no. 25, pp. 897-901.
[6] Tatschl R. et al. Advanced Turbulent Heat Transfer Modeling for IC-Engine Applications Using AVL FIRE. International Multidimensional Engine Modeling User's Group Meeting. Detroit, 2006, pp. 1-10.
[7] Kavtaradze R.Z., Onishchenko D.O., Zelentsov A.A., Sergeev S.S. The Influence of Rotational Charge Motion Intensity on Nitric Oxide Formation in Gas-Engine Cylinder. International J. of Heat and Mass Transfer, 2009, vol. 52, no. 19-20, pp. 4308-4316.
[8] Kavtaradze R.Z., Onishchenko D.O., Zelentsov A.A., Finkel'berg L.A., Kostyuchenkov A.N. Simulation of processes in the "inlet collector-cylinder" system of aviation reciprocating engine with distributed fuel injection. Vestn. Mosk. Gos. Tekh. Univ. im. N.E. Baumana, Mashinostr. [Herald of the Bauman Moscow State Tech. Univ., Mech. Eng.], 2012, no. 4, pp. 3-16 (in Russ.).
[9] Lanshin A.I., Finkelberg L.A., Kostyuchenkov A.N., Zelentsov A.A., Bakanov M.A. Investigation of Inlet Flow Swirl Influence on the Gasoline Aviation Piston Engine Characteristics. Vestnik VGTU [The Bulletin of VSTU], 2012, vol. 8, no. 2, pp. 96-99 (in Russ.).
[10] Kavtaradze R.Z., Sergeev S.S. New Alternative (Partially Homogeneous) Combustion Process as a Method for Reduction of Concentrations of Nitric Oxides and Soot in Combustion Products of Diesel. High Temperature (TVT), 2014, vol. 52, no. 2, pp. 294-309.
[11] Kavtaradze R.Z., Zelentsov A.A., Skripnik A.A. 3D Investigation of Features of Realization of Working Processes of Aircraft Piston Engines. Proceedings International User Conference Advanced Simulation Technologies (AST), 18-20 June 2013, Austria, Graz. 2013, pp. 1-9.
Статья поступила в редакцию 27.03.2015
Кавтарадзе Реваз Зурабович — д-р техн. наук, профессор кафедры "Поршневые двигатели" МГТУ им. Н.Э. Баумана.
МГТУ им. Н.Э. Баумана, Российская Федерация, 105005, Москва, 2-я Бауманская ул., д. 5.
Kavtaradze R.Z. — D.Sc. (Eng.), Professor, Department of Piston Engines, Bauman Moscow State Technical University.
Bauman Moscow State Technical University, 2-ya Baumanskaya ul. 5, Moscow, 105005 Russian Federation.
Зеленцов Андрей Александрович — канд. техн. наук, доцент кафедры "Поршневые двигатели" МГТУ им. Н.Э. Баумана.
МГТУ им. Н.Э. Баумана, Российская Федерация, 105005, Москва, 2-я Бауманская ул., д. 5.
Zelentsov A.A. — Ph.D. (Eng.), Associate Professor, Department of Piston Engines, Bauman Moscow State Technical University.
Bauman Moscow State Technical University, 2-ya Baumanskaya ul. 5, Moscow, 105005 Russian Federation.
Просьба ссылаться на эту статью следующим образом:
Кавтарадзе P.3., Зеленцов А.А. Влияние формы впускных каналов на эффективные и экологические показатели среднеоборотного дизеля // Вестник МГТУ им. Н.Э. Баумана. Сер. Машиностроение. 2015. № 6. C. 59-73.
Please cite this article in English as:
Kavtaradze R.Z., Zelentsov A.A. Influence of intake port shape on effective and ecological parameters of middle-speed diesel engine. Vestn. Mosk. Gos. Tekh. Univ. im. N.E. Baumana, Mashinostr. [Herald of the Bauman Moscow State Tech. Univ., Mech. Eng.], 2015, no. 6, pp. 59-73.
