автореферат диссертации по энергетическому, металлургическому и химическому машиностроению, 05.04.02, диссертация на тему:Улучшение показателей двигателей внутреннего сгорания на основе совершенствования локальных характеристик процессов газообмена
Автореферат диссертации по теме "Улучшение показателей двигателей внутреннего сгорания на основе совершенствования локальных характеристик процессов газообмена"
ХАРЬКОВСКАЯ ГОСУДАРСТВЕННАЯ АКАДЕМИЯ ЖЕЛЕЗНОДОРОЖНОГО ТРАНСПОРТА
„ р НА ПРАВАХ рукописи
V <
И
ПЕЛЕПЕЙЧЕНКО Владимир Игоревич
УЛУЧШЕНИЕ ПОКАЗАТЕЛЕЙ ДВИГАТЕЛЕЙ ВНУТРЕННЕГО СГОРАНИЯ НА ОСНОВЕ СОВЕРШЕНСТВОВАНИЯ ЛОКАЛЬНЫХ ХАРАКТЕРИСТИК ПРОЦЕССОВ ГАЗООБМЕНА
05.04.02 - тепловые двигатели
АВТОРЕФЕРАТ
ДИССЕРТАЦИИ НА СОИСКАНИЕ УЧЕНОЙ СТЕПЕНИ ДОКТОРА ТЕХНИЧЕСКИХ НАУК
Харьков - 1995
- г -
Работа выполнена на кафедре "ТЕПЛОТЕХНИКА И ТЕПЛОВЫЕ ДВИГАТЕЛИ" Харьковской государственной академии железнодорожного транспорта.
Научный консультант - доктор технических наук,
профессор САХАРЕВИЧ В.Д.
Официальные оппоненты - доктор технических наук,
профессор ДРАГАНОВ Б.Х.
- доктор технических наук, профессор ЕАГМУТ Г.А.
- доктор технических наук, профессор СТЕФАНОВСКИИ B.C.
Ведущая организация -Институт проблем машиностроения HAH Украины, г.Харьков
Защита состоится "Мг ".^ЙС.. 1995г. в. И часов в аудЖЗ на заседании специализированного совета Д 02.15.02 по тепловым двигателям при Харьковской государственной академии железнодорожного транспорта, 310050, г. Харьков-50,пл. Фейербаха 7.
С диссертацией и авторефератом можно ознакомиться в библиотеке Харьковской государственной академии железнодорожного транспорта.
Автреферат разослан "fff. "./ff? . 1995г.
Ученый секретарь специализированного совета Д 02.15.02, к.т.н. ,доц. ЛЯЛЮК В.М.
- 3 - ■
ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ, АКТУАЛЬНОСТИ И СТЕПЕНИ ИССЛЕ-ДОВАННОСТИ ТЕМАТИКИ
Необходимш условием получения высоких технико-экономических и экологических показателей дизеля является обеспечение требуемых показателей газообмена. Для оценки качества газообмена используют показатели, основными из которых являются коэффициент избытка воздуха при сгорании, коэффициент остаточных газов, коэффициент наполнения, отнесенный к условиям во впускном ресивере, коэффициент продувки.
Перечисленные показатели являются интегральными, так как характеризуют итог сложных нестационарных процессов, протекающих в цилиндре и коллекторах в период газообмена. Однако они не несут никакой информации о характере движения заряда в цилиндре в разные моменты времени и в камере сгорания (КС) в период топливоподачи.
В настоящее время общепризнано исключительно важное значение организованного движения заряда для качественного смесеобразования и достижения высокого индикаторного КПД.
При проектировании и доводке двигателей возникает проблема количественной оценки влияния конструктивных и режимных факторов на показатели газобмена и локальные значения параметров состояния заряда в цилиндре и КС. Кроме информации о псшях скоростей в некоторых случаях необходимы данные о локальных значениях температуры и концентрации остаточных газов или топлива.
Применение эффективных методик расчета указанных' пара-
метров дало бы возможность повысить уровень проектирования, целенаправленно формировать характеристики двигателей.
Течение заряда в цилиндре и КС в общем случае представляет собой нестационарное трехмерное турбулетное движение газовой смеси, сопровождающееся трением о стенки и теплообменом при изменяющихся во времени граничных условиях.
В настоящее время существует ряд нерешенных прин^ципи-альных проблем, связанных как с построением математической модели такого течения, так и с выбором эффективного метода решения системы нелинейных дифференциальных уравнений в частных производных, образующих модель.
Известные методы расчета носят либо частный характер, либо имеют недостаточную точность, что ограничивает возможность их эффективного применения.
В представленной диссертационной работе содержится развитие и обобщение методов расчета внутрицилиндровых течений, что позволило повысить их точность и расширить область применения на новый класс задач. Внедрение б практику проектирования созданного на основе теоретических разработок программного комплекса позволяет квалифицировать работу как актуальную, направленную на решение крупной научно-технической проблемы, имеющей важное прикладное значение - снижение затрат на создание и доводку дизелей и повышение их технического уровня за счет рационального выбора параметров, определяющих характер процессов газообмена.
ЦЕЛЬ И ОСНОВНЫЕ ЗАДАЧИ ИССЛЕДОВАНИЯ Целями исследования являются:
1. Разработка программного комплекса, позволяющего на стадии проектирования определять интегральные показатели газообмена и локальные значения параметров состояния заряда в цилиндре и камере сгорания в периоды газообмена и сжатия.
2. Изучение основных особенностей и закономерностей внутрицилиндровых течений заряда в двигателях разных типов и разработка практических рекомендаций по повышению их технико - экономических показателей.
Для достижения указанных целей необходимо было решить следующие основные задачи.
1. Разработать метод решения системы уравнении,входящих -в математическую модель, описывающую нестационарное трехмерное турбулентное течение двухкомпонентного сжимаемого газа.
2. Разработать метод учета турбулентного обмена в модели нестационарного течения.
3. Обосновать выбор модели / гипотезы / турбулентности, учитывающей особенности течения в цилиндре ЛВС.
4. Уточнить, а в необходимых случаях разработать способы задания граничных и начальных условий.
5. Разработать и отладить комплекс прикладных программ, реализующих математическую модель течения заряда в цилиндре.
6. Разработать методику проведения исследования по выбору параметров системы воздухоснабжения конкретного двигателя.
7. Разработать методику количественной оценки степени равномерности распределения топлива в объеме камеры сгора-
ния.
ХАРАКТЕРИСТИКА МЕТОДОЛОГИИ,МЕТОДА ИССЛЕДОВАНИЯ ПРЕДМЕТА И ОБЪЕКТА
Предмет исследования - вопрос о связях между конструктивными особенностями двигателя,режимом его работы и мгновенными локальными параметрами заряда в цилиндре в периоды газообмена и сжатия.
Объект исследования - цилиндры четырехтактных и двухтактных двигателей.
Методология исследования - моделирование газодинамических процессов.
Метод исследования - расчет на ЭВМ мгновенных локальных параметров потока в цилиндре, использование экспериментальных данных для идентификации математических моделей и подтверждения теоретических выводов.
ОБОСНОВАНИЕ ТЕОРЕТИЧЕСКОЙ И ПРАКТИЧЕСКОЙ ЦЕННОСТИ ИССЛЕДОВАНИЯ, ЕГО НАУЧНОЙ НОВИЗНЫ Теоретическую ценность имеет методика математического моделирования нестационарного пространственного турбулентного течения двухкомпонентного сжимаемого газа в замкнутой области сложной и меняющейся во времени формы при нестационарных граничных условиях.
Практическую ценность имеют;
- программный комплекс,позволяющий рассчитывать локальные параметры потока в цилиндре и КС (скорость, концентрацию, давление, температуру, плотность, интенсивность турбулентности), а также интегральные показатели газообмена;
- методика оценки степени равномерности распределения компонентов смеси (топливо-воздух) в объеме КС;
- данные о закономерностях формирования полей параметров потока в цилиндрических КС дизелей автотракторного типа;
- уравнения регрессии, устанавливающее связь мевду показателями газообмена дизеля 6ДН12/2x12, режимом его работы и конструктивными параметрами окон;
- рекомендации по выбору параметров системы воздухо-снабжения и конструкции окон модификации транспортного дизеля 6ДН12/2х12, предназначенной для дизель-поезда.
- данные о характере распределения газообразного топлива в объеме камеры сгорания дизеля 6ГЧН31.8/33.
Научную новизну имеют:
1. Алгоритм решения системы уравнений газовой динамики, объединяющий подходы методов крупных частиц (МКЧ) и распада произвольного разрыва (РПР).
2. Метод учета турбулентного переноса в математических моделях течения, ориентированных на использование алгоритмов методов расщепления по физическим факторам.
3. Метод коррекции масштаба турбулентности в модели Н.И.Булеева в зависимости от числа Ричардсона и формула для вычисления числа Ричардсона в условиях течения в цилиндре. Применение коррекции обеспечило учет влияния градиентов плотности, возникающих во вращающемся потоке, на пульсаци-онную составляющую скорости.
4. Опытные данные о значении коэффициента трения движущегося воздушного заряда о стенку цилиндра.
5. Метод оценки степени равномерности распределения
компонентов газовой смеси в объеме КС и рабочие формулы для вычисления интегрального показателя равномерности.
6. Уравнения в форме полиномов второй степени, устанавливающие связь между конструктивными и режимными параметрами и показателями газообмена дизеля 6ТД (6ДН12/2х12).
7. Результаты исследования основных закономерностей эволюции локальных параметров внутрицилиндровых потоков в двигателях разных типов.
УРОВЕНЬ РЕАЛИЗАЦИИ, ВНЕДРЕНИЯ НАУЧНЫХ РАЗРАБОТОК
Диссертация выполнена в соответствии с программой ГКНТ СССР N 535 от 16.05.89 " Создание и освоение производства новых типов двигателеи внутреннего сгорания и агрегатов на их основе", постановлением СМ СССР N 1041 от 28.08.86 "О мерах по увеличению использования в XII пятилетке природного газа как моторного топлива", планами НИР и ОКР Харьковского конструкторского бюро по двигателестроению (ХКБД).
Результаты научного исследования в виде программ расчета локальных параметров заряда в цилиндре и интегральных показателей газообмена переданы в ХКБД , на Харьковскии завод тракторных двигателеи, в Головное специализированное конструкторское бюро по двигателям средней мощности, НТП КБ среднеоборотных двигателеи, используются в практике НИР и учебном процессе в ХДРГАЖТ.
Рекомендации по выбору параметров системы воздухо-снабжения и конструкции окон тепловозного варианта дизеля 6ДН12/2х12 переданы в ХКБД.
- у -
АПРОБАЦИЯ И ПУБЛИКАЦИЯ РЕЗУЛЬТАТОВ НАУЧНОГО ИССЛЕДОВАНИЯ, СТРУКТУРА И ОБЪЕМ ДИССЕРТАЦИОННОЙ РАБОТЫ Апробация работы, основные положения и результаты диссертационной работы доложены на Всесоюзных научно-технических конференциях:
"Проблемы создания и использования двигателей с высоким наддувом", г. Харьков, 1979г.
"Создание и совершенствование быстроходных дизелей", г.Токмак, Юждизельмаш, 1982г.
"Современный уровень и пути совершенствования экономических и экологических показателей двигателей внутреннего сгорания", г.Луганск, 1983г.
"Теория и расчет мобильных машин и двигателей внутреннего сгорания", г. Тбилиси, 1985г.
"Актуальные проблемы развития двигателей внутреннего сгорания и дизельных установок", г. Ленинград, 1990г. на Международных научно-технических конференциях: "Компьютер-.наука, техника, технология, здоровье", Харьковский политехнический институт - Мишкольцкий университет, Венгрия, 1993г.,1995г.
на научно-технических конференциях Харьковской государственной академии железнодорожного транспорта, 1979-1994гг.
Публикации: по результатам диссертационного исследования опубликована 21 научная работа.
Объем и структура работы. Диссертация состоит из введения, пяти глав, заключения, списка использованных литературных источников, включающих 226 наименований.
- lü - 1
Работа содержит 254 страницы основного текста, 7 таблиц и 67 рисунков.
ЛИЧНЫЙ ВКЛАД В РАЗРАБОТКУ НАУЧНЫХ РЕЗУЛЬТАТОВ, ВЫНОСИМЫХ НА ЗАЩИТУ
В ходе выполнения работы лично автором получены следующие научные результаты, которые выносятся на защиту.
1. Метод решения системы уравнений газовой динамики, объединяющий подходы МКЧ и PHP.
2. Метод учета турбулентного переноса в математических моделях течения, ориентированных на использование алгоритмов методов расщепления по физическим факторам.
3. Метод коррекции масштаба турбулентности в модели Н.И.Булеева в зависимости от числа Ричардсона и формула для вычисления числа Ричардсона в условиях течения в цилиндре, позволяющие учесть влияние градиентов плотности газа, находящегося в поле массовых сил, на интенсивность турбулентности.
4. Методика обработки опытных данных об изменении поля скоростей заряда в процессе сжатия, позволяющая определить значение коэффициента трения Еоздуха о стенки цилиндра .
5. Метод оценки равномерности распределения топлива в объеме КС по суммарному числу столкновений молекул топлива и воздуха в единицу времени.
6. Методика проведения исследований по выбору параметров системы воэдухоснабжения и конструктивных параметров окон тепловозной модификации дизеля 6ДН12/2Х12.
- 11 -СОДЕРЖАНИЕ РАБОТЫ Во введении обоснована актуальность темы диссертационной работы. Отмечен постоянный и устоичивыи интерес ведущих исследовательских центров и фирм, таких как МГТУ, ХПУ, С-ПбПУ, ЧелПУ, С-ПбКИ, МАДИ, МАМИ, МИИТ, РИИЖТ, ХабИИЖТ, ХабПИ.ВгПИ, ВПИ, ЛИВТ, НИВТ, ВПИ, УфАИ, ХАИ, БелПИ, ЗалИИ, ХИМЭСХ, КАДИ, ХАДИ, АзИНГ, ПО "3-д им. Малышева", "Коломенскии з-д", БМЗ, "Звезда", "Ювдизельмаш", "Первомаисвдизельмаш", "Пенздизель-маш", "Дальдизель", КАМЮ, ГАЗ, ВгМЗ, ЦКИДИ, ХКБД, ГСКБД, ИШаш, ИММ, НАТИ, НАШ, AYL, Ricardo, Ley land , B&W, Fiat, Sulzer, GMT, GMC, Pielstick, Cummins, Perkins, и других к проблеме совершенствования процессов газообмена в двигателях разных типов, развития теории и методов исследования процессов движения заряда в цилиндре.
Обоснованы и сформулированы цели и основные задачи исследования, изложенные выше. Указаны возможные области расширенного применения результатов работы.
ПЕРВАЯ ГЛАВА содержит анализ известных методов расчета интегральных показателей газообмена и локальных параметров газового потока в цилиндре и камере сгорания.
Эти методы можно условно разделить на три основные группы:
- термодинамические;
- основанные на использовании уравнения закона сохранения момента количества движения заряда относительно оси цилиндра;
- газодинамические.
Показаны ограниченные возможности методов первой и второй групп. Сделан вывод, что только газодинамические методы позво
-ляют определять локальные параметры потока при течении в областях сложной формы.
Рассмотрена общая система уравнений газовой динамики, описывающая нестационарное трехмерное турбулентное течение вязкого сжимаемого газа.
Приведен• список возможных допущений, которые используются для упрощения исходной системы.
Показано, как при принятии той или иной комбинации допущений из общей системы уравнений получаются частные модели, описанные в литературе и применявшиеся для решения конкретных задач. Выполнен анализ преимуществ и недостатков частных моделей, обусловленных уровнем соответствия или несоответствия принятых допущений реальной картине развития течения и на этой основе установлены границу применимости частных моделей.
Показано, что на сегодняшний день наиболее совершенными по I остановке являются модели турбулентного движения заряда в цилиндре, предложенные в работах М.Г.Кругдова, С.Р.Березина, О.Н.Агапитова, О.Н.Лебедева, А.В.Бастерса, В.П.Хлынина, М.А.Гороховского, С.Н.Чиркова, В.П.Закржевского, М.В.Дульге-ра, Е.М.Иткиса.
Для замыкания системы уравнений турбулентного движения должны быть использованы дополнительные соотношения, которые дают модели (гипотезы) турбулентности.
Проведен анализ преимуществ и недостатков известных моделей турбулентности и рассмотрен опыт применения этих моделей для расчетов течений в цилиндрах ДВС. Сделан вывод, что для решения прикладных задач целесообразно применять модели
- 13 - '
локальной турбулентности, которые могут обеспечить хорошее совпадение результатов расчета и эксперимента при минимальном усложнении полной системы уравнений. Особо выделена модель локальной турбулентности Н.И.Булеева, поскольку в отличие от других она учитывает особенности пространственного течения в каналах сложной формы. Однако эта модель еще не применялась для полностью замкнутых объемов, в том числе и для цилиндров ДВС. Эмпирические коэффициенты, входящие в модель, определены по данным о течении в трубах. Поэтому правомерность использования модели Н.И.Булеева для условий движения в цилиндре должна быть доказана.
Приведена классификация условий однозначности и сделан анализ практики их задания в известных работах. Показано, что задание условия "прилипания" потока на стенке при размерах сетки на 1...2 порядка больше толщины пограничного слоя приводит к искусственному, "схемному", увеличению толщины пограничного слоя, а следовательно, к искажению картины течения вблизи стенки. Выходом могло бы быть сгущение сетки в пристеночной области, однако это приводит к резкому увеличению общего числа узлов сетки и времени решения задачи.
Рассмотрен метод объединения аналитических решений для распределения скорости в пограничном слое с численным решением для основной области течения. Показано, что относительная погрешность в определении этим методом напряжений трения на стенке в условиях течения в цилиндре ДВС может достигать 250Х.
На основании анализа литературных данных определены пути совершенствования модели течения заряда в цилиндре ДВС.
- 14 - '
ВО ВТОРОЙ главе приведена классификация методов решения системы уравнении газовой динамики и рассмотрена практика применения этих методов для решения системы уравнений,входящих в известные математические модели, описывающие эволюцию параметров потока в цилиндре и КС. Показано, что чем сложнее постановка задачи, тем ухе круг возможных методов.
Для нелинейной систеш дифференциальных уравнений в частных производных, описывающих нестационарное трехмерное турбулентное течение сжимаемого двухкомпонентного газа б области сложной и изменяющейся во времени формы при нестационарных граничных условиях необходимо применение численных методов. В настоящее время наибольшее распространение получили метод "крупных частиц" (МКЧ) О.М.Белоцерковского и Ю.М.Давыдова и метод "распада произвольного разрыва" (PIP) С.К.Годунова. Выполнено сопоставление этих методов. Отмечено, что каждый из этих методов имеет свои преимущества и недостатки. МКЧ, относясь к группе методов расщепления, позволяет использовать многопроцессорные ЭВМ. Это качество особенно важно при решении многовариантных или оптимизационных задач. Но в МКЧ нет строгого обоснования формул для определения параметров на границах расчетных ячеек,что приводит к в случае расчетов пространственных течении к накоплению ошибок и потере устойчивости разностной схемы.
В методе РПР формулы на разрывах / на гранях ячеек / получены строго, но в предположении, что состав среды в смежных ячейках одинаков. Однако в случае расчета продувки цилиндра ДВС среду следует рассматривать как состоящую минимум из двух компонентов - воздуха и остаточных газов, причем нез-
- 15 - '
бежно существование градиентов их концентрации. Метод РПР не ориентирован на применение многопроцессорных ЭВМ. .
Разностные схемы обоих методов консервативны и обладают "схемной вязкостью", однако этой вязкости во многих случаях недостаточно для адекватного воспроизведения реальных дисси-пативных процессов, обусловленных переносом турбулентными молями массы, импульса, энергии и концентрации из одной области течения в другую. Практика применения метода РПР и МКЧ для расчета турбулентных течений многокомпонентных газов чрезвычайно ограничена.
Исходя из изложенного, в настоящей работе была поставлена задача разработки метода и алгоритма решения системы уравнении газовой динамики, в котором были бы соединены положительные свойства МКЧ и РПР и учитывался турбулентный обмен.
ГЛАВА ТРЕТЬЯ содержит описание разработанной математической модели и программного комплекса для расчета показателей газообмена и характеристик движения заряда в цилиндре и КС. Приведен анализ причин худших результатов применения в алгоритме МКЧ разностных формул второго порядка точности по сравнению с формулами первого порядка. Использование формул второго порядка для расчета параметров на гранях ячеек эквивалентно принятию допущения о том, что малые возмущения в потоке распространяются со скоростями, большими скорости звука. Противоречие реальной физической картине и приводит к худшим результатам. В МКЧ скорость и давление на гранях ячейки определяется как среднеарифметическое значение по параметрам в двух смежных ячейках. Эти формулы введены без фи-
- lü - 1
зического обоснования. При этом оказывается, что скорость на гранях не зависит явно от распределения давления по смежным ячейкам, а давление на грани не зависит явно от распределения скорости. Разумеется, существует неявная зависимость, однако она проявляется через один или более временных шагов. Временной сдвиг между изменением параметров в ячейках и реакцией на это изменение параметров на гранях ячейки может в условиях расчетов трехмерных течений приводить к неустойчивости разностной схемы.
В методе РПР параметры на разрывах, то есть фактически на гранях ячейки, определяются решением уравнений законов сохранения. Поэтому предлагается оставить алгоритм МКЧ, предусматривающий расщепление по физическим факторам на "эйлеров", "лагранжев" и заключительный этап, но параметры на гранях определять по формулам метода РПР. Могут использоваться либо линеаризованные формулы для "больших" величин, либо "акустическое" приближение, имеющее вид:
а?„.= ^ -Г- - ■ cd
2-рЛо
Р." =
рл ♦ р.:,.
р.-с.-
П ь
c*'>j ы
(2)
где U и Р - скорость и давление; j>e и Сс~ плотность и скорость звука.
- 17 - '
Целые индексы означают, что данная величина относится к параметрам в ячейке, дробные - что данная величина относится к грани, например с номером .
Выражения в квадратных скобках полностью совпадают с формулами МНИ для параметров на гранях. Выражения в фигурных скобках учитывают свойства газа и состояние потока в окрестности рассматриваемой точки.
Оценка относительных различий в определении параметров на гранях по методам РПР и МКЧ в условиях течения в цилиндре показала, что различие в скорости может достигать 16%, в то время как различие в давлении и плотности не превосходит 3,1%. Поэтому целесообразно корректировать только формулы цля определения скорости на гранях ячейки.
Разработана методика учета турбулентного обмена. Известно, что в турбулентных течениях вне пограничного слоя коэффициенты турбулентной вязкости и диффузии на несколько юрядков превосходят соответствующие молекулярные коэффици-гнты, поэтому эффектами молекулярного переноса можно пренеб-эечь. Тогда диссипативные процессы будут определяться переюсом движущимися тубулентными молями массы, импульса и »нергии из одной области течения в другую и последующим »ассеиванием этих молей. Предлагается в алгоритм МКЧ вклю-[ить механизм учета такого переноса.
Рассматривается система интегральных уравнений законов ¡охранения, имеющая вид:
- 1В - '
+ =0, (3)
V У
+ (4)
V 5 $
^р-Е)-^^ ^РьУ^-с^гО, (5)
5 5
IX2
Р-- > Е = £ + • (б,?)
Принято, что через каждую из поверхностей, ограничивающих ячейку, во встречных направлениях движутся турбулентные моли, причем объем моля, покидающего ячейку, равен объему моля, входящего в ячейку через данную поверхность. Допущение о равенстве объемов молей, пересекающих границу ячейки, позволяет исключить на "эйлеровом" этапе учет работы их проталкивания. Турбулентность мелкомасштабная , то есть длина свободного пробега моля и его характерный размер меньше размеров ячейки. Моль, вошедший в ячейку, теряет свою индивидуальность. Из принятых допущений следует, что на поверхности каждой грани ячейки существует три участка с разными значениями нормальной составляющей скорости потока йп. Схема распределения скоростей приведена на рис.1.
На участке А скорость равна усредненной скорости на поверхности . На участках В и С скорости определяются из условий
М'иЬ = М/П + И/и 5 Пне : Й'и ~ М'И з
где 11^- пульсационная составляющая скорости. Поверхности участков связаны соотношением
s&= Sc - 4"S0 ,
где ^ - эмпирический коэффициент, характеризующий относительное сечение моля. При принятых допущениях Ч^О.Б.
На практике возможны два варианта соотношения скоростей Йи и и/„. Вариант "а": и вариант "б" |и/„|>|й„|.
При варианте "а" моли из ячейки, находящиеся "ниже" по течению, не могут проникать в ячейку, находящуюся "выше" по течению, так как сносятся потоком.
Учет турбулентного переноса осуществляется в случае, когда JH'hJ>|Hh|.
На "лагранхевом" этапе вычисляются интегралы по поверхности ячейки, выражающие в формулах (3), (4), (5) потоки массы, импульса и энергии. Их можно представить в виде:
s l-.t
+ E^if^g-p^xj-Cl^i-lii.i), (в)
где £ - число поверхностей, ограничивающих ячейку; дБ- площадь грани ячейки, ( X =1 для формулы (3), *-Ц, для (4), Х=Е для (5) ).
Свободный подстрочный индекс [*] означает,что данная величина соответствует параметрам ячейки, из которой происходит течение через поверхность с номером 1 , индекс "а" -что данная величина относится к параметрам в рассматриваемой ячейке, индекс "в" - что данная величина относится к параметрам в соседней ячейке.
Значение ty следует принимать равным 0,5. Физически это означает, что вся область течения охвачена турбулентным дви-
жением. Первая суша в правой части (8) представляет поток, обусловленный усредненной скоростью, вторая - поток, вызвая-' ньш турбулентными пульсациями.
Для практического использования формулы (8) необходимо каким-либо образом определить значение пульсационной скорости И.
Анализ известных моделей турбулентности показал, что "модель локальной турбулентности Н.И.Вулеева позволяет учесть влияние формы области течения на величину пульсационной скорости . Согласно этой модели модуль пульсационной скорости определяется как
oli .
где JU, - эмпирический коэффициент (уЦ =1,25),| модуль деформации поля скоростей усредненного движения; L - масштаб турбулентности, вычисляемый по формуле
ь «г J е > сю)
где С- расстояние от рассматриваемой точки до стенки канала в направлении, задаваемом углом СЭ . В случае замкнутой области интегрирование необходимо выполнять в пределах телесного угла 4 . Выполнено сопоставление расчетных значений пульсационных составляющих скорости, найденных по формуле (9), с опытными данными о распределении компонентов вектора скорости усредненного и пульсационного движений в цилиндрах и КС двигателей разных типов. Отмечено, что относительная разница расчетных и экспериментальных значений и' в некоторых случаях не превосходит 10...20%, а в некоторых
достигает 200% и более.
Сделан вывод, что причиной нестабильных результатов непосредственного применения модели (9) к течениям в цилиндре является отсутствия в ней учета влияния градиента плотности среды, находящейся в поле массовых сил, на характеристики турбулентности.
Градиенты плотности и массовые силы порождаются закруткой заряда относительно оси цилиндра. Для учета этих эффектов предложено корректировать масштаб турбулентности Ь с использованием формулы, аналогичной по своей структуре формуле Монша-Обухова, применяющейся в физике атмосферы:
Ь=Ц,*(1+р-1?1). (И)
где масштаб, вычисляемый по формуле (10); р- опытный коэффициент; И- безразмерное число Ричардсона, вычисляемое по формуле „ - -2.
В диссертации приведен вывод формулы для определения локального значения I?! в условиях цилиндра и КС:
Л /Эй,чЛ~а
Ш= —-—--(13)
■С \9н / '
где Ц ~ локальное, значение тангенциальной составляющей
скорости; *2. - расстояние до оси цилиндра; С - скорость звука
при параметрах газа на оси цилиндра; II- деформация поля
' V и 1
скоростей.
Установлено, что наилучшие результаты обеспечиваются при
р =0,4. Предложенная корректировка масштаба турбулентности позволила получить относительную разницу расчетных и экспериментальных данных о турбулентных пульсациях скорости, не превышающую во всех случаях 201.
В разработанной модели принято, что на гранях пристеночных ячеек действуют напряжения трения, вызванные движением газа относительно стенки. Эти напряжения можно рассчитать по формуле
- С^ • "¿Г" . (14)
где коэффициент трения, Ц,- скорость газа относительно стенки, определяемая по сеточным параметрам для пристеночной ячейки.
В литературе отсутствуют данные о значении в условиях движения в цилиндре ДВС. Поэтому была разработана методика определения Cj.no опытным данным об изменении параметров движущегося заряда при сжатии. Обработка опытных данных по двигателям разных типов позволила сделать вывод, что в диапазо-5" 6
не 10 $ 5?е 4 10 С^ автомоделей по числу Рейнольдса и равен 0.012.
На основе изложенной методики разработан программный комплекс, ориентированный на ресурсы 1ВМ РС АТ. Он состоит из двух основных модулей. Первый представляет собой программу расчета замкнутого рабочего цикла комбинированного двигателя. второй - программу расчета параметров движущегося в цилиндре заряда. В первом модуле реализована термодинамическая модель процессов, которая позволяет определять изменение во времени термодинамических параметров состояния рабочего тела в каждом из цилиндров, впускных и выпускных кол-
- 23 - '
лекторах, трубопроводах в зависимости от набора большого числа конструктивных и режимных факторов. Модуль может использоваться и автономно для решения задач оптимизации параметров системы газотурбинного наддува. Но он не позволяет рассчитывать локальные параметры заряда в цилиндре. Роль первого модуля в программном комплексе заключается в определении начальных и граничных условий для расчета процессов течения заряда в цилиндре при газообмене.
Связь между модулями осуществляется через общую базу данных, из которой второй модуль получает информацию о давлении и температуре газа в коллекторах . Результат работы второго модуля - интегральные показатели газообмена и информация о полях давлений, температур, концентраций остаточных газов, компонент векторов скорости для любого момента времени процессов газообмена и сжатия. Данные о текущих локальных параметрах потока выводятся в процессе счета в графической форме на монитор в виде векторных полей скорости или в виде графиков, показывающих изменение по радиусу составляющих вектора скорости и доли остаточных газов. Это позволяет изучать развитие течения в цилиндре в динамике. Предусмотрена возможность получения в любой момент времени копии графического экрана на принтере.
Адекватность модели и программы расчета параметров движущегося заряда проверялась сопоставлением результатов расчета полей скоростей с экспериментальными данными.
Для примера на рис. 2 показаны опытные и расчетные графики изменения осевой и тангенциальной составляющих скорости в цилиндре дизеля 6ДН12/2х12 при его статической продувке. В
диссертации приведены результаты проверки на адекватность и по данным для двигателей других типов.
В ЧЕТВЕРТОЙ ГЛАВЕ приведены примеры решения прикладных .задач с использованием разработанного программного комплекса Выполнен анализ структуры поля скоростей в цилиндрической КС дизеля 6ЧН 10,5/12. Рассмотрены варианты несмещенной и смещенной относительно оси цилиндра оси КС. Моделирование течения выполнялось на трехмерных сетках. На рис.3 показаны выведенные на принтер векторные поля скоростей в разные моменты сжатия. В нижней части рис.3 показана схема течения. На рис.4 приведены диаграммы изменения компонентов вектора скорости по радиусу на трех уровнях по глубине КС. Поле скоростей в смещерной КС не обладает осевой симметрией. В большей мере это относится к осевой составляющей, в меньшей - к тангенциальной, которая меняется не только по радиусу, но и по глубине камеры. Осевая и радиальная составляющие скорости вблизи стенок камеры малы. Поэтому основную роль в смесеобразовании играет тангенциальная составляющая, возникающая из-за предварительной закрутки заряда. Разрыв графика радиальной составляющей скорости в центре камеры на уровне А-А связан исключительно с принятым правилом знаков - скорость положительна, если ее вектор направлен от центра к периферии.
Информация о состоянии поля скоростей позволяет обоснованно выбирать направление топливных факелов.
Выполнено исследование по выбору параметров системы воздухоснабжения и конструкции окон транспортного двигателя
- 25 - :
6ДН12/2х12, предназначенного для дизель-поезда.
В цилиндре двухтактного двигателя с прямоточной продувкой и закруткой заряда на впуске имеет место сложное пространственное турбулентное движение воздуха и остаточных газов. Существуют условия для образования плохо продуваемой зоны у оси цилиндра. В этой зоне может возникать вредное явление - течение, направленное навстречу основному потоку продувочного воздуха.
Для примера на рис. 5 показано векторное поле скоростей, на рис. 6 - диаграммы изменения компонентов вектора скорости и доли остаточных газов по радиусу цилиндра в сечении форсуночного пояса в период продувки при максимальном открытии продувочного окна.
Большой набор параметров, влияющих на газообмен, делает необходимым разработку методов, позволяющих на этапах проектирования и доводки определять их рациональное сочетание.
Ставилась задача установить, какими должны быть значения фаз газораспределения, давлений во впускном и выпускном коллекторах наибольшего угла закрутки продувочного окна по высоте*^ чтобы на режимах тепловозной характеристики обеспечить заданное сочетание массы свежего воздуха в цилиндре М$, коэффициента остаточных газов ^ , коэффициента продувки,^, а также интенсивности движения заряда в КС в конце сжатия. За меру интенсивности движения заряда принималась Путях" максимальная величина тангенциальной составляющей скорости.
Рассмотрен вариант, когда закрытие выпускных и продувочных окон происходит одновременно.Тогда при постоянном угле заклинки коленчатых валов значение всех фаз газораспреде-
- 25 - '
ления определяется какой-либо одной из них, например, углом начала открытия выходного окна Неодновременная реализация заданных интегральных показателей газообмена и интенсивности движения заряда является необходимым условием получения высокого индикаторного КОД.
Задача решалась в три этапа. На первом выполнена серия расчетов процессов газообмена и сжатия при различных сочетаниях факторов, влияющих на газообмен. Сочетания задавались в соответствии с G-опгимальным планом пятифакторного эксперимента с варьированием на трех уровнях. Определены коэффициенты в уравнениях регрессии, устанавливающие зависимости
) Ü* = УьЛ.Рг,**,*») I У =
I =
где п - частота вращения коленчатого вала. Уравнение регрессии имеет вид полиномов второй степени с парным взаимодействием и каждое содержит 21 коэффициент. На втором этапе для режима номинальной мощности (п=2200мин~', N - 552 кВт) определяли значения Чс, ^з обеспечивающие заданные значения Mj.ll^ V». Для этого численными методами решали систему четырех нелинейных уравнений регрессии вида (15) . Получено, что заданные значения М =4,2.10 кг, Y=0,06, Фа. =1,2, =42м/с при п=2200 мин. могут быть реализованы при Р5=0,252Ша,Рг=0,22МПа, % -108? п.к.в., =40° . Há третьем этапе проверялись показатели газообмена, которые будут получены на режимах тепловозной характеристики.
При правильно подобранных параметрах системы воздухо-
снабжения не возникает обратных течений в осевой зоне и "перезавихрения" заряда в период топливоподачи.
На рис.7...рис.9 показаны параметры дизеля 6ДН12/2Х12 по тепловозной характеристике в серийном и опытном вариантах. Реализация рекомендованных параметров системы позволила снизить удельный эффективный расход топлива на режимах тепловозной характеристики на 10..12 г/(кВт.ч).
Повышение экономичности достигнуто благодаря перераспределению воздуха, поступающего в двигатель- увеличению коэффициента избытка воздуха для сгорания и уменьшению коэффициента продувки.
Выполнен анализ влияния различия давлений воздуха на входе в продувочные окна, расположенные с противоположных сторон гильзы цилиндра дизеля 6ДН12/2х12, на интегральные показатели газообмена и состояние поля скоростей заряда.
Расчет течения выполнялся на трехмерной сетке. Приведены диаграммы распределения компонентов векторов скорости и доли остаточных газов по радиусу в различных сечениях цилиндра. Показано, что несмотря на возникающую асимметрию полей скоростей в цилиндре в период продувки, интегральные показатели газообмена, вычисленные с учетом неравномерности распределения и вычисленные по среднему давлению, практически одинаковы (относительная разница не превосходит 52.)
Приведены результаты расчета локальных концентраций газообразного топлива (метана), подаваемого в цилиндр газодизеля Д50 (6ГЧН31.8/33) на такте впуска через клапан, установленный на месте клапана пневматического пуска. На рис.10 показано распределение локальных значений коэффициента из-
бытка воздуха в двух взаимно-перпендикулярных осевых сечениях цилиндра в момент начала подачи запальной порции жидкого топлива. В дизеле Д50 отсутствует предварительная закрутка заряда при впуске,а канал для подачи газа смещен от оси цилиндра к периферии. Это приводит к образованию вблизи выходного отверстия канала зоны, где локальные коэффициенты избытка воздуха меньше единицы. В процессе сжатия заряда равномерность распределения топлива в объеме КС несколько улучшается, однако из-за малой интенсивности движения заряда в цилиндре переобагащенная зона существует вплоть до момента начала подачи топлива.
Информация о локальных концентрациях газообразного топ-дива в КС позволяет оптимизировать процесс смесеобразования, верно выбирать ориентацию сопел распылителя форсунки, через которые подается жидкое запальное топливо, исключить его подачу в зоны, переобагащениые газообразным топливом.
Разработанный программный комплекс позволяет исследовать влияние конструктивных и режимных факторов на пространственное распределение концентраций газообразного топлива и проводить оптимизацию процесса смесеобразования.
ПЯТАЯ ГЛАВА посвящена разработке метода оценки равномерности распределения топлива в объеме КС.
В работах Ю.Б.Свиридова, А.С.Лышевского. Я.М.Майера, Н.Ф.Разлейцева, А.П.Кудряша, И.И.Тимченко, М.Л.Копылова, С.В.Сахаревича.Д.А.Москаленко и других предложены методики расчета локальных концентраций топлива в объеме КС в период топливоподачи. Характер макрораспределения топлива зависит от геометрии камеры, конструкции распылителей, характеристик
топливолодачи, интенсивности движения воздушного заряда. При проведении расчетов смесеобразования получают, как правило, серию различных макрораспределений и тогда возникает проблема их сравнения. Структура макрораспределений столь сложна, что без применения специальных методов невозможно дать им даже чисто качественную оценку, по принципу "более равномерно" или "менее равномерно" распределено топливо. Известные методы оценки равномерности исходят из двухзонной модели, которая предполагает, что какая-то часть воздуха доступна, а какая-то недоступна для сгорания. Локальный характер распределения внутри "доступной" зоны при этом не рассматривается.
Нами в качестве показателя равномерности распределения топлива предлагается использовать отношение числа столкновений в единицу времени молекул топлива в паровой фазе с молекулами воздуха при реальном распределении топлива в объеме
КС к числу столкновений, которое было бы при равномерном распределении того же количества топлива и воздуха в камере. При расчете числа столкновений следует учитывать только те из них, которые происходят в зонах, где концентрация топлива не выходит за пределы воспламенения.
В молекулярно-кинетической теории газов показано, что число столкновений молекул топлива и воздуха в единице объема в единицу времени пропорционально произведению плотности и 62 числа разнородных частиц
6г . (16) .
С использованием соотношения для газовых смесей получена формула для вычисления показателя равномерности
е с -30'
- = ~--- , (17)
Zog hogn0r
где t- число ячеек, в которых концентрация находится в пределах воспламенения;относительный объем I -ой ячейки; число молей воздуха и топлива в элементарной ячейке с номером I \ tf<Tобщее число молей воздуха и топлива в объёме КС.
Для случая, когда макрораспределение представлено в виде локальных коэффициентов избытка воздуха, получена формула
е
, Г 1 г l_
(18)
где и среднее по цилиндру и локальное значение
ср
коэффициента избытка воздуха; мт и /ц - молекулярные массы топлива и воздуха; х - доля топлива, выгоревшего к данному моменту времени.
При равномерном распределении ^ =1, при предельно неравномерном ^ =0. Приведены примеры расчетов показателя равномерности. Для газодизеля 6ЧН31.8/33 с подачей газа на такте впуска через отверстие клапана пневмопуска ^ =0,66, для камеры дизеля 10Д100 ^ =0,6, для открытой цилиндрической камеры дизеля 6ЧН10.5/12 ^ =0,38. Использование показателя равномерности позволяет выявить резервы совершенствования процесса смесеобразования.
ОСНОВНЫЕ РЕЗУЛЬТАТЫ РАБОТЫ И ВЫВОДЫ 1. В цилиндрах ДВС, представляющих замкнутую область
сложной и меняющейся во времени формы в периоды газообмена и сжатия имеет место пространственное турбулентное течение многокомпонентной сжимаемой газовой смеси при нестационарных условиях взаимодействия со стенками /трение, теплообмен / и со смежными областями /впускным и выпускным коллекторами /.
При построении математической модели такого течения возникает проблема обоснованного выбора гипотезы турбулентности, способов задания начальных и граничных условий, констант, входящих в эмпирические формулы. Отсутствует единое мнение о целесообразности применения различных методов решения системы уравнений, входящих в математическую модель. Это сдерживает создание эффективных прикладных программ, позволяющих рассчитывать локальные параметры потока и интегральные показатели газообмена.
Создание методик расчета внутрицилиндровых течений и программных комплексов, позволяющих решать прикладные зада-чи-актуальная проблема теории ДВС.
2. Разработан метод решения системы уравнении газовой динамики, объединяющий подходы метода крупных частиц и метода распада произвольного разрыва. Суть этого метода в том, что сохраняется схема расщепления исходных уравнений по физическим факторам, свойственная МКЧ, но параметры на гранях ячеек определяются с использованием соотношений метода РПР.
Это позволяет получить устойчивую разностную схему, аппроксимирующую уравнения трехмерного течения газа.
3. Предложен метод учета турбулентного обмена в рамках общей модели, объединяющей подходы РПР и МКЧ. Метод предпо-
латает выделение на "лагранжевом" этапе расчета из общих потоков через грань ячейки потока, вызванного усредненной скоростью, и потока, вызванного турбулентными пульсациями скорости. Пульсации скорости должны определяться с использованием какой-либо модели турбулентности. Из уравнений предложенного метода как частный случай получаются известные уравнения метода крупных частиц.
4. Обоснован выбор модели турбулентности для условии течения заряда в цилиндре ДВС. Показано, что модель локальной турбулентности Н.И.Вулеева позволяет учесть специфику течения в замкнутом объеме сложной формы, каким является цилиндр и камера сгорания. Однако эта модель не применялась ранее для расчета закрученных течений, в которых существуют градиенты плотности потока. Это приводит в некоторых случаях к относительным погрешностям определения пульсационной скорости, превышающим 250 7..
Предложена формула для корректировки масштаба турбулентности в модели Н.И.Вулеева в зависимости от числа Ричардсона, и формула для вычисления числа Ричардсона в условиях течения в цилиндре. Применение этих формул обеспечило погрешность определения пульсационной скорости не более 20 %.
5. На основании анализа опытных данных об изменении поля скоростей заряда в процессе сжатия в двигателях разных типов показано, что коэффициент трения газа о стенку автомоделей по числу Рейнольдса и равен 0,012. Сделан вывод о неправомерности использования аналогий с обтеканием пластин для вычисления коэффициента трения в условиях течения газа в ци-
линдре.
6. Разработана и отлажена программа расчета параметров нестационарного трехмерного турбулентного течения заряда в цилиндре в процессах газообмена и сжатия. Программа позволяет определять поля всех параметров газового потока, а также поля концентраций для любого момента времени. Рассчитываются также интегральные показатели газообмена.
7. Разработан программный комплекс, состоящий из модуля расчета замкнутого рабочего цикла многоцилиндрового дизеля с газотурбинным наддувом, охлаждением наддувочного воздуха, многоколлекторной выпускной системой и модуля расчета пространственного течения в цилиндре и камере сгорания в периоды газообмена и сжатия.
Модуль расчета замкнутого рабочего цикла дизеля обеспечивает определение начальных и граничных условий для расчета течения заряда в цилиндре. Связь между модулями осуществляется через общую базу данных.
Сопоставление расчетных и экспериментальных данных позволило сделать вывод, что относительная погрешность определения компонентов локальных усредненных скоростей 1 не превосходит 15%,что приемлемо для решения прикладных задач.
8. Количественную оценку степени равномерности распределения топлива в объеме КС предлагается выполнять с помощью показателя равномерности, вычисляемого как отношение числа столкновений молекул топлива и воздуха при реальном распределении топлива к числу столкновений, которое было бы при равномерном распределении в объеме КС того же количества топлива и воздуха. Показатель равномерности позволяет срав-
нивать макрораспределения топлива при различной организации процесса смесеобразования.
9. Выполнен теоретический анализ закономерностей развития течения в цилиндрических камерах сгорания двигателеи автотракторного типа. Исследовано течение в несмещенной и смещенной относительно оси цилиндра камере. Показано, что существуют градиенты тангенциальной составляющей скорости по глубине камеры. Смещение камеры относительно оси цилиндра приводит к асимметрии распределения осевых и радиальных составляющих вектора скорости вдоль радиуса КС, но тангенциальная составляющая остается осесимметричной.
10. В соответствии с планом пятифакторного эксперимента выполнена серия расчетов, позволивших получить уравнения множественной регрессии, связывающие интегральные показатели газообмена и характеристики движения заряда в камере дизеля 6ДН12/2х12 с конструктивными параметрами продувочных окон, параметрами газа во впускном и выпускном коллекторах, скоростным режимом работы двигателя. Уравнения имеют вид полиномов второй степени и содержат 21 коэффициент. Они могут быть использованы при создании различных модификации данного дизеля.
11. Разработаны рекомендации по выбору конструктивных параметров окон и параметров системы воздухоснабжения тепловозной модификации дизеля 6ДН12/2х12. Требуемый уровень показателей газообмена и интенсивности вращения заряда в КС может быть получен при углах открытия выходного окна
= 108°, продувочного окна = 118°, угле закрутки продувочного окна по высоте = 40°, давлении наддува
- ЗЬ - 1
{5 = 0,252 МПа, среднем давлении в выпускном коллекторе
= 0,220 МПа. Реализация указанных параметров позволила снизить удельный эффективный расход топлива по тепловозной характеристике на 10...12 г/(кВт.ч) по сравнению с базовым вариантом двигателя.
12. Исследовано влияние неравномерности распределения давления наддувочного воздуха у продувочных окон дизеля 6ДН12/2х12 на показатели газообмена. Показано что неравномерность приводит к асимметрии поля скоростей в цилиндре в период продувки. Однако интегральные показатели газообмена, вычисленные с учетом неравномерности распределения давления и по условному среднему давлен™ воздуха у окон практически одинаковы. Это позволяет выполнять расчет газообмена двигателя данного типа по среднему давлению у окон.
13. Выполнено исследование распределения метана в КС газодизеля Д50 при подаче газа на такте впуска через отверстие клапана лневмопуска. Показано наличие переобогащенных зон, при попадании в которые запальное жидкое топливо будет находиться вне пределов воспламенения. Определены значения показателя равномерности "распределения газообразного топлива в КС.
14. Результаты работы в виде программного комплекса внедрены в Харьковском конструкторском бюро по двигателе-строению, на Харьковском заводе тракторных двигателеи, ГСКБ по двигателям средней мощности, НТП КБ среднеоборотных двигателеи, используются при выполнеии НИР и в учебном процессе в ХарГАЖТ. Рекомендаций по выбору параметров тепловозной модификации дизеля 6ДН12/2х12 внедрены ХКБД.
ОСНОВНЫЕ ПОЛОЖЕНИЯ ДИССЕРТАЦИИ ОПУБЛИКОВАНЫ В РАБОТАХ
1. Пелепейченко В.И. Моделирование турбулентного движения заряда в цилиндре двигателя внутреннего сгорания на основе метода крупных частиц.- Харьк. гос. . академия ж.-д. транспорта, Харьков, 1994.13с. Деп. в ГНТБ Украины, N 127-Ук 95, 16.01.95.
2. Пелепейченко В.И. Определение мгновенных параметров газового потока и характеристик турбин при исследовании дизелей с импульсными системами наддува.// Тезисы докладов Всесоюзной научно-технич. конф."Проблемы создания и использования двигателей с высоким наддувом".- Харьков, 1979. -С. 193.
3. Пелепейченко В.И. Модификация метода крупных частиц для расчета трехмерного нестационарного движения заряда в цилиндре ДВС.- Харьк.гос. академия ж.-д. транспорта.- Харьков, 1994.-13с. .-Деп. ГНТБ Украины. N 126-Ук 95, 16.01.95.
4. Пелепейченко В.И. Идентификация модели турбулентного движения воздушного заряда в камере сгорания дизеля.- Харьк. гос. академия ж.-д. транспорта.-Харьков, 1994.-13с.-Деп. ГНТБ Украины, N 129-Ук 95, 16.01.95.
5. Пелепейченко В.И., Назаренко Д. А., Крушедольский А. Г. Метод экспериментального определения основных параметров газового потока в коллекторах ДВС. // Совершенствование сельскохозяйственных, тракторов и автомобилей: Сб. научных трудов. МИИСП, -г XV, ВЫП.4.-М. 1978.-с.69...74.
6. Каминский В.Н., Мамин Б.В., Пелепейченко В.И. Выбор системы газотурбинного наддува для дизеля КамАЗ-740.//Тезисы
докладов Всесоюзной научно-техн. конф. "Проблемы создания и использования двигателей с высоким наддувом".- Харьков, 1979. - с.191-192.
7. Симеон А.Э., Сахаревич В.Д., Пелепейченко В.И.,Пет-росянц В.А. Оптимизация конструктивных параметров дизеля и турбокомпрессора. //Тезисы докладов Всесоюзной научно-техн. конф. "Актуальные проблемы развития двигателей внутреннего сгорания".- Л.: Судостроение, 1990.-с.16-17.
8. Крушедольский А.Г., Пелепейченко В.И., Симеон А.Э., Шаройко H.A. Математическое моделирование рабочего цикла развернутого дизеля 10Д100 при его работе в эксплуатации.-Харьк. институт инж. ж.-д. транспорта.- Харьков, 1985.-10с. Деп. ЦНШТЭИТЯЖМАШ 26.03.1985. N 3020.
9. Пелепейченко В.И., Сахаревич В.Д. Математическая модель для расчета параметров нестационарного движения заряда
• в цилиндре в период наполнения-сжатия. "Компьютер: наука, техника, технология, здоровье": //Тезисы докладов международной научно-техн. конф,- Харьков, 1993.-е.
10. Линник A.B., Курдюков А.П., Гапон В.В., Пелепейченко В.И. Четырехтактный дизель с газотурбинным наддувом. A.C. СССР N 1870274, зарегестр. 1.10.1987.
11. Пелепейченко В.И. Симеон А.Э., Леховицер М.А., Кравченко A.B., Гапон В.В. Двигатель внутреннего сгорания A.C. СССР N 1257262, зарегестр. 15.05.1986.
12. Симеон А.Э., Сахаревич В.Д., Пелепейченко В.И. Метод определения-расходной характеристики турбины, работающей в импульсной системе наддува ДВС. //Вопросы теплообмена и гидродинамики транспортных и промышленных установок.-Вып.
634.-М. :МЖГ, 1979.- с.69-79.
13. Симеон А.Э.-, Счастный Е.Е., Лесовицкий И.В., Пеле-пейченко В. И. Оценка путей реализации характеристики двигателя постоянной мощности дизеля 6ЧН12/14. Деп. ЦНИИТЭТЯЖ-МАШ.4.09.89, N 445, ТМ-89.
. 14. Пелепейченко В.И., Карнаухов Ю.И. Оценка резервов повышения экономичности дизеля 6ЧН13/11.5 путем выравнивания показателей газообмена цилиндров правого и левого рядов. Деп. ЦНШГЭИТЯЖМАШ, 24.06, 1985, N 1491, ТМ-85.
15. Пелепейченко В.И., Сахаревич В.Д. Показник для виз-начення р1вном!рност1 розпод1лу палива у камер1 згорян-ня.//Тези допов1дей 55-i науково1 конференцИ спещал1ст1в зал1зничного транспорту.- XapKiB. iHCT. 1нж. зал1зничн. транспорту, 1993. -с. 17.
16. Курдюков А.П., Пелепейченко В.И., Сахаревич В.Д., Таратушка В.А. Выбор системы ГТН перспективного дизеля 12ЧН14/14 для маневрового тепловоза.//Повышение топливной экономичности тепловозных дизелей и теплоэнергетических установок. Межвуз. сборник научн. трудов.- Харьков: ХИИТ, 1986.-е.51-63.
17. Сахаревич В.Д., Пелепейченко В.И. Способ оценки равномерности распределения топлива в камере сгорания дизеля при расчетах смесеобразования. - Харьк. гос. академия ж.-д. транспорта.-Харьков, 1994.-8с. - Деп. в ГНГБ Украины. N 123 - Ук.95, 16.01.95.
18. Алехин С.А., Пелепейченко В.И. Определение коэффициента трения вращающегося воздушного заряда о стенки цилиндра двигателя внутреннего сгорания.- Харьк. гос. академия
ж.-д. транспорта.- Харьков, 1994, -13с. Деп. в ГНТБ Украины. N 124 - УК 95. 16.01.95.
19. Алехин С.А., Пелепейченко В.И. Многофакторные зависимости для определения показателей газообмена цилиндра дизеля ДН12/2х12.- Харьк. гос. академия ж.-д. транспорта. Харьков, 1994.-12с. Деп. в ГНТБ Украины. N 128 - Ук 95. 16.01.95.
20. Алехин O.A., Пелепейченко В.И. Выбор параметров системы воздухоснабжения тепловозной модификации дизеля 6ДН12/2Х12 из условия достижения заданных показателей газообмена. -Харьк. гос. академия ж.-д. транспорта.- Харьков, 1994. - 15с. - Деп. в ГНТБ Украины. N 125 - Ук 95, 16.01.95.
21. Пелепейченко В.И., Сахаревич В.Д. Расчетно- теоретический анализ состояния поля скоростей в смещенной относительно оси цилиндра камере сгорания дизеля.- Харьк. гос. академия ж.-д. транспорта. - Харьков, 1994. - 8с. -Деп. в ГНТБ Украины. N 1121-Ук 95, 11.05.95.
Варианты соотношения между усредненной и пульсационной скоростями
& 1 1 >1 1 а/ ь
1
А и-* и,» А
С С
г
а) ' и, Рлс.1 в)
Компоненты вектора скорости в сечении форсуночного пояса дизеля ДН12/2х12
й 1.4
1,0 0,8 0,6 о/. 0.2 0
-Ч -2. 0 2 ^ ЯЛ0,м.
— опыт, -расчет
Рис.2
Осевой Тангенциальный
/ / ~ /» \\
/ / ч\
V
л \\ /V / '
\\ у * ' у X /
V г/ //
т и и
1
^ / /* ^ —
_Ю и/с
I > . | \ \ \ I * \ 1 »N11
ч \ \ 1
^=328°
/ / / у- ^ -—
—10 М/с
\ N
- - ч
- — N
ч\\ Ч \ 1 Г —
' • 1 } ^ \ \ \ V
\ ч V \
/ 1 Ч V \ / !
< ^ , X 1 ■ ✓ 1
1 V . ! 1 / ! 1
1 \ ч \ I / / 1 / 1 1
1 \ л \ 1 / / / / ; 1
\ \ \ \ 1 I 1 1 1 Г 1
"р.-ззГ
— 10 п/с
\ \ ч
• » I
/ I I
/ - _ N N
у _ ч ч
/
г i
i i
Схема течения
Рис.3
Компоненты вектооа скорости в различных сечениях смещенной камеры сгорания
Чг,м/с 10
0
-5
ЧД 0
-5 -10 V* /#
15 10
5 о
А/ Л"
Б
.С
-0 -П.Ь -0,2 □ 0 ОА 0,6 0,& +
Рис Л
Поле скоростей в цилиндре дизеля Ш12/2х12
Компоненты векторов скорости, доля остаточных газов £
т
50
г кг
0,5
Рг
ц
1,2
1.0
__-- —
±1 ** у-
у / X > / / ^ - * ■ / — « •
* к/
г > И'тах .
/
% / ■
500 300
тСС
•• <50
■ № ■ ■ 50
г(о
1.0
•40 ..20
1 Ш 1200 «00 2000 П, МИН Рис.7
Рис. 8
-------ОПЫТНАЯ £>е3 "/'ОМ. ОХЛОХ/V.
-*- оиитиа;» с пром. окланСА ■
Рис. 3
Распределение газа в камере сгорания дизеля Д50 в момент начала подачи запальной порции топлива
/ ?=ззбГ ^ =о.еб /
Рис. -10
•ьъ-
Pelepeichenko V.I. Improving of the Indices of Internal Combustion Engines by means of Perfection of Local Gas Exchange Characteristics.
The dissertation for presenting doctor's degree of technical sciences on speciality 05.04.02-heat engines, Kharkov State Academy of Railway Transport, Kharkov, 1995.
19 research works and 2 author certificates, which include the method of the calculation of unstationary multi-composition multi-dimensional turbulent gas flow paramétrés in cylinders and combustion chambers during intake and compression, are presented. The guantitative inflluence of the engine construction and thermo-dynamic paramétrés upon the local flow paramétrés and the gas exchange characteristics has been determined. The results of the study have been used in research practice and design.
Пелепейченко В.И. Улучшение показателей двигателей внутреннего сгорания на основе совершенствования локальных характеристик газообмена.
Диссертация на соискание ученой степени доктора технических наук по специальности 05.04.02-тепловые двигатели, Харьковская государственная академия железнодорожного транспорта, Харьков, 1995.
Защищается 19 научных работ и 2 авторских свидетельства, содержащих методику расчета параметров нестационарного турбулентного многокомпонентного газового потока в цилиндрах ЛВС в периоды газообмена-сжатия. Установлено количественное влияние конструктивных и режимных параметров двигателя на локальные параметры потока, характеристики газообмена.
Результаты исследовании внедрены в практику НИР и ОКР. Ключов1 слова: газообмен, поля параметр!в, економ!чн1сть.
-
Похожие работы
- Методы исследования и пути совершенствования процессов газообмена и сгорания в бензиновых двигателях
- Влияние нестационарных явлений на статические и динамические характеристики двигателей внутреннего сгорания
- Повышение эффективности функционирования силовых агрегатов в приводе транспортных машин
- Исследование процесса сгорания и образования токсичных веществ в двигателе с двухстадийным сгоранием
- Повышение эффективности процесса сгорания в тракторных дизелях совершенствованим элементов систем впуска и управления топливоподачей
-
- Котлы, парогенераторы и камеры сгорания
- Тепловые двигатели
- Машины и аппараты, процессы холодильной и криогенной техники, систем кондиционирования и жизнеобеспечения
- Машины и агрегаты металлургического производства
- Технология и машины сварочного производства
- Вакуумная, компрессорная техника и пневмосистемы
- Машины и агрегаты нефтяной и газовой промышленности
- Машины и агрегаты нефтеперерабатывающих и химических производств
- Атомное реакторостроение, машины, агрегаты и технология материалов атомной промышленности
- Турбомашины и комбинированные турбоустановки
- Гидравлические машины и гидропневмоагрегаты
- Плазменные энергетические и технологические установки