CC BY
17
УДК 621.431.74.016
Б.И. Руднев
Дальневосточный государственный технический рыбохозяйственный университет,
690087, г. Владивосток, ул. Луговая, 52б
АНАЛИЗ ВЗАИМОДЕЙСТВИЯ ПОТОКОВ ИЗЛУЧЕНИЯ В КАМЕРЕ СГОРАНИЯ ДИЗЕЛЬНОГО ДВИГАТЕЛЯ
Представлена теоретическая оценка взаимодействия потоков излучения и проанализированы основные особенности системы интегральных уравнений радиационного теплообмена применительно к камере сгорания дизельного двигателя.
Ключевые слова: дизельный двигатель, камера сгорания, потоки излучения.
B.I. Rudnev
THE ANALYSIS OF FLUXES INTERACTION RADIATION IN DIESEL ENGINE COMBUSTION CHAMBER
Theoretical estimation of fluxes interaction radiation is given in this article. Main peculiarities of system of integral acquisitions of radiation heat transfer concerning diesel engine combustion chamber. Key words: diesel engine, combustion chamber, fluxes of radiation heat transfer.
Излучающую систему для условий камеры сгорания (КС) дизельного двигателя можно представить как состоящую из четырех зон: трех поверхностных - поверхности крышки цилиндра F1, втулки F2 и поршня F3, а также одной объемной зоны V, заполненной сажегазовой средой. Следуя [1], общую картину взаимодействия потоков излучения в произвольной точке М]с граничной поверхности F можно представить схемой, показанной на рисунке.
Рассмотрим более подробно воздействие основных потоков излучения для точки М^ принадлежащей граничной поверхности F. В каждую точку граничной поверхности F приходит собственное излучение от объемной зоныMi(V), т.е. Есоб(М). Причем для всех точек граничной поверхности Мk, Мп, Мр и ^^этот поток является падающим потоком.
В то же время каждая из поверхностей, на которых расположены точки Мk, Мп, Мр и M имеет свой эффективный (или исходящий [2]) поток излучения:
Еэф (Мп ) = ЕсобЕ (Mn )+ ЕотрЕ (Mn \ (1)
в свою очередь, отраженный поток для точки Мп может быть представлен следующим образом:
EornPF(Mn) = Епад^, (2)
где RF - коэффициент отражения граничной поверхности F.
Учитывая взаимосвязь между коэффициентом отражения RF и степенью черноты поверхности sF, имеем
EompF (Mn ) = EnadRF = (l — SF KadF = (l — SF КобУ (M ) (3)
Схема взаимодействия потоков излучения в произвольной точке Мк граничной поверхности F Scheme of fluxes interaction radiation in diesel engine combustion chamber at random point Mk
adjacent surface F
Тогда эффективный поток излучения от точки Мп граничной поверхности Г можно представить следующим образом:
ЕэфГ (Мп ) = ЕсобГ (Мп ) + (1 ^ КобУ (М) (4)
Аналогично для точек Мр и М ^:
ЕэфГ (МР )= ЕсобГ МР )+ (1 - ^ КобУ (Мг) (5)
ЕэфГ (Мг)= ЕсобГ (МГ)+ (1 - ^ )ЕсобУ (Мг ) (6)
Таким образом, для произвольной точки Мк граничной поверхности Г (см. рисунок) падающими потоками излучения будут ЕсобУ(Мг), Еэфр(МП), Еэфр(Мр) и ЕэфГ(Му) .
Результирующий поток излучения для точки Мк в общем виде запишется так:
Ерез (Мк )= Епад (Мк )-ЕэфГ (Мк ) (7)
Очевидно, что (7) эффективный поток излучения от точки Мк можно представить аналогично эффективным потокам от точекМп, Мр и М^ (4-6):
ЕэфГ (Мк ) = ЕсобГ (Мк )+ (1 - % КобУ (Мг) (8)
115
Подставляя в (7) соответствующие потоки излучения, определяемые по зависимостям (3-5, 7), получим
Ерез (М к ) = ЕсобУ (Мг ) + Е эфЕ (Мп ) + ЕэфЕ (м р ) + ЕэфЕ (м у)- Е эфЕ {Мк ) = = ЕсобУ (МI) + ЕсобЕ (Мп ) + (1 - СЕ )ЕсобУ (Мг ) + ЕсобЕ (Мр ) + + ( - СЕ )ЕсобУ (Мг ) + ЕсобЕ ((М у ) ++ ( - СЕ )ЕсобУ (Мг ) -- ЕсобЕ (Мк ) - (1 - СЕ )ЕсобУ (М1).
Группируя в (9) члены, содержащие собственный поток излучения объемной зоны ЕсобУМ,), имеем
Ерез (Мк ) = ЕсобУ (М1 ) + (1 - Ср )ЕсобУ (М1 ) + (1 - Ср )ЕсобУ (м1) +
+ (1 - Ср )ЕсобУ ((Ц- ) + (1 - Ср )ЕсобУ ((Ц ) - (1 - Се )ЕсобУ (м," ) + (10)
+ ЕсобЕ (Мп ) + ЕсобЕ (Мр ) + ЕсобЕ (Му ) - ЕсобЕ (Мк )
Проанализируем полученную зависимость (10) для результирующего потока излучения произвольной точки Мк граничной поверхности Е, выражающую функциональную связь основных потоков излучения. Из (10) следует, что результирующий поток излучения для точки Мк формируется при взаимодействии собственного потока излучения объемной зоны ЕсобУ(М;) и собственных потоков излучения поверхностей, образующих граничную поверхность Е ЕсобЕ (мп ), ЕсобЕ (мр Е, ЕсобЕ (Мг) и ЕсобЕ (мк Е .
Причем собственный поток излучения поверхности, которой принадлежит точка Мк, имеет противоположный знак по отношению к собственным потокам излучения других поверхностей. Члены уравнения (10) вида (1 -Ср )Есобу (м)представляют собой отраженный поток излучения от соответствующих элементов граничной поверхности Е. Заметим, что при построении математической модели локального радиационного теплообмена в КС дизеля необходимо величину степени черноты граничной поверхности Ср задавать дифференцированно по отношению к каждому ее элементу, т.е. отдельно
для крышки цилиндра, поршня и втулки.
На основании экспериментальных данных [3], величины степени черноты для отдельных поверхностей, составляющих КС, имеют следующие значения: для днища поршня из алюминиевого сплава, покрытого нагаром толщиной 0,1-0,2мм, с = 0,91-0,95, а при толщине нагара 0,2-0,6 мм, с = 0,94-0,97; для крышки цилиндра из чугуна, покрытой нагаром, с = 0,90-0,95; для втулки цилиндра с тонким слоем масла с = 0,3-0,4. Принимая во внимание величины степени черноты поверхностей КС, следует отметить, что члены вида (1 -Ср )Есо(^у (м) в уравнении (10) важны в определенной мере лишь для поверхности втулки цилиндра. Для поверхности крышки цилиндра и поршня эти члены укладываются в погрешность определения температуры излучателя в КС дизеля и концентрации частиц сажи в ней. Напомним, что указанная погрешность составляет 15-20 %. Однако даже и для втулки цилиндра величина (1 -Ср)Есобу(м)не может быть значительной, так как в период активного процесса радиационного теплообмена в КС (от ВМТ до 80-100° поворота коленчатого вала после ВМТ) большая часть втулки цилиндра перекрыта поршнем. Особенно это относится к дизелям с полуразделенными КС типа ЦНИДИ, ЯМЗ, КАМАЗ и др. Что же касается членов уравнения (10), представляю-
116
щих собой собственные потоки излучения элементов граничной поверхности Г:
ЕсобГ (Мп ), ЕсобГ (Мр ), ЕсобГ Му) и ЕсобГ (Мк ) ,- то з«есь необходимо отметить следующее. Собственный поток излучения поверхности однозначно определяется ее температурой и излучательной способностью (степенью черноты) по зависимости
ЕсобГ = а08ЕТГ.
(И)
Аналогично можно определить и собственный поток излучения объемной зоны
ЕсобУ = а08УТУ-
(12)
В результате сложного взаимодействия рассмотренных выше потоков излучения между объемной зоной и соответствующими точками граничной поверхности в замкнутой системе (см. рисунок) формируются результирующие потоки, определяющие радиационную тепловую нагрузку деталей КС. Этот процесс взаимодействия основных потоков излучения описывается системой неоднородных интегральных уравнений излучения (уравнения Фредгольма второго рода). Следуя [1, 3, 4], приведем указанную систему уравнений для условий в точке Мк:
т
Е пад (Мк )=Е! Е эф (Мг ) Ь(М,, Мк )У
1=1 У
11Еэф(Мг ) К(МГ,Мк
+
+
у=1у
Еэф (Мк ) = Есоб (Мк ) + Кк
11Е эф (м Г)К (м г, Мк )
+
Г=1у
+
т
Еэф (Мг) Ь(М г,Мк )У
1=1у Ерез (Мк ) = гк
11 Еэф (Мг, Мк ) К (Мг, Мк )
Г=1Г
+ Х /Еэф (Мг) ь( ,Мк )у
1=1 У
- Есоб (Мк ),
(13)
где Кк, гк - коэффициенты отражения и поглощения (степени черноты) поверхности Г в
точке Мк; М^- произвольная точка поверхности Г, лежащая вне окрестности точки Мк;
Мг - произвольная точка объема У.
Интегральные операторы К (М у; Мк ) и Ь(М1; М к ) называются ядрами уравнения,
являются геометрическими функциями точек объема и поверхности и характеризуют систему как единый оптико-геометрический объект [1]. Проинтегрировав функции К и Ь по всей граничной поверхности, получим следующие уравнения связи геометрических и разрешающих локальных угловых коэффициентов излучения (ЛУКИ) [2]:
Ф(Мг, Мк ) = | К(мг,Мк № = ехр (- т}р(М7,Мк )
F
Ф (М,, Мк ) = 11 (м Мк ) = ехр (-г>(М|, Мк ),
(14)
Ь ( ,, Мк = ехр (-г)>(М|., Мк ),
F
где т = КЬ - число Бугера.
Следуя [1], под <р(Мг> Мк) в (14) имеется в виду геометрический локальный угловой коэффициент излучения от элемента поверхности dFy в окрестностях точки М^ на элемент поверхности dFk в окрестностях точки Мк, показывающий, какая доля эффективного излучения поверхности dFy достигает поверхности dFk с учетом многократных отражений на границе системы. Коэффициент (р[М у;Мк) имеет тот же смысл для излучения элемента объема dVi в окрестностях точки М. на поверхности dFk . В свою очередь, ф(Мг;Мк ) и ф(м .;Мк ) есть разрешающие коэффициенты, отличающиеся от
соответствующих геометрических ЛУКИ тем, что они учитывают ослабление потока излучения вследствие поглощения и рассеяния энергии в объеме.
Приведенная система интегральных уравнений (13) вытекает из интегродифферен-циального уравнения переноса излучения, соответствующих граничных условий и соотношений, связывающих интенсивность излучения с различными характеристиками полусферического (поверхностного) и сферического (объемного) излучений [2]. По своей физической сущности интегральные уравнения (13) являются уравнениями состояния излучающей системы, поскольку описывают ее состояние как на границах, так и в объеме среды [2].
Для построения математической модели локального радиационного теплообмена в КС дизеля воспользуемся итерационно-зональным методом, предложенным Ю.А. Су-риновым [5], а также работами, в которых раскрыты и проанализированы основные особенности применения зональных методов к расчету радиационного теплообмена в различных энергетических установках [2, 6].
В зональных методах расчета теплообмена излучением для систем, образованных излучающей, поглощающей и рассеивающей средами, ограниченными излучающими и отражающими границами, исходя из предложения, излучающая система может быть приближенно заменена эквивалентной, составленной из оптически однородных изотермических подсистем. При этом поверхность границы излучающей системы F подразделяется на п изотермических оптических однородных непрозрачных зон:
п
F = ^ Fk , а объем V, заполненный ослабляющей и излучающей средами, - на т объ-
к=1
емных изотермических оптических однородных зон с полупрозрачными границами:
т
V = X V . В [1] показано, что итерационно-зональный метод обобщает все основные
к=1
постановки задач о радиационном теплообмене, имеющих единственное решение и различающихся способом задания граничных условий и характеристик поля излучения. Математическую основу этого метода составляет конечная система линейных алгебраических уравнений, аппроксимирующих соответствующую систему интегральных уравнений излучения типа (13).
Как известно [1, 5], для реальных излучающих систем Я Ф 1; е Ф 1 и все интегральные уравнения излучения (13) имеют единственное решение. Для его получения воспользуемся итерационно-зональным методом алгебраической аппроксимации [5]. Сущность этого метода заключается в разделении реальной излучающей системы на конечное число объемных и граничных зон, в пределах которых оптические, энергетические и температурные характеристики могут быть осреднены или заданы аналитически [1, 5]. В этом случае исходное интегральное уравнение излучения может быть заменено соответствующей системой линейных алгебраических уравнений, число которых равно числу зон в системе. Алгебраическая аппроксимация интегрального оператора возможна как в исходном интегральном уравнении, так и в уравнении резольвенты. В связи с этим точное решение исходного интегрального уравнения может быть выражено через резольвенту интегрального оператора или непосредственно через разрешающие угловые коэффициенты излучения. В соответствии с этим алгебраическая интерпретация решения может быть представлена в общем виде или через резольвенту [1, 5]. Алгебраическая аппроксимация точного решения исходных интегральных уравнений излучения (13) согласно [1] будет иметь вид
п ! \ ! \ т
Епад(Мк) = XЕсоб(мг) ф(мг,¥к)+ XЕсоб(М) ф(му,¥к);
у=1 У=1
ЕэфМк)= XЕсобМк) Ф((к,Ег)-
у=1
п / \ / \ т
+ (1 -е)ХЕсобМг) Ф(Му,Ек)+(1 -ек)ХЕсоб(Му) Ф(Му,Ек);
г=1
У=1
Ерез (Мк )=ек
±ЕсобМг) ф(мг,Еу)
7=1
+
тп
X Есоб (Му, Ек )-Х Есоб (Мк ) Ф(М к, Еу)
у=1
Г=1
(15)
Уравнения (15) получены в предложении, что излучающая система состоит из п граничных и т объемных зон, а в качестве первой итерации принято условие равенства эффективных потоков излучения в пределах зоны соответствующим потокам собственного излучения этих зон.
Путем замены в интегральном операторе переменного ЛУКИ постоянным, осред-ненным в пределах рассматриваемой зоны, можно систему интегральных уравнений резольвент поверхностного и объемного излучений привести к соответствующим системам алгебраических уравнений. Принимая, что в пределах каждой граничной зоны приближенно выполняется условие (р{мЕк )= (р.^ , в [1, 5] получены следующие системы аппроксимирующих уравнений резольвент поверхностного излучения:
Ф(му,Ек)- ¿Я^Ф(г,Ек )= (р(ыг,Ек).
у=1
и объемного излучения:
Ф(М,,F)- Е V*ф(мv,,F)l
и
(17)
Y=1
где - средний угловой коэффициент излучения зоны у на зону к, определяемый по зависимости
В результате этих действий система интегральных уравнений излучения (13) приводится к эквивалентной конечной системе линейных алгебраических уравнений. Решение последней позволяет определить мгновенное локальное распределение по элементам граничной поверхности всех видов потоков излучения.
Проведенный анализ системы интегральных уравнений и теоретические оценки взаимодействия потоков излучения позволили сделать ряд существенных уточнений в зональной модели радиационного теплообмена в КС дизельного двигателя [7] и других энергетических установок [8].
1. Петриченко Р.М. Элементы автоматизированного проектирования ДВС: Алгоритмы прикладных программ [Текст] / Р.М. Петриченко, С.А. Батурин, Ю.Н. Исаков [и др.] - Л.: Машиностроение, 1990. - 328 с.
2. Блох А.Г. Теплообмен излучением: справ. [Текст] / А.Г. Блох, Ю.А. Журавлев, Л.Н. Рыжков. - М.: Энергоатомиздат, 1991. - 432 с.
3. Modest M.F. Radiative Heat Transfer. - New York: Academic Press, 2003. - 822 p.
4. Dombrovsky L.A. Radiation Heat Transfer in Disperse Systems. - New York: Begell House, 1996. - 256 p.
5. Суринов Ю.А. Об итерационно-зональном методе исследования и расчета лучистого теплообмена в поглощающей и рассеивающей среде [Текст] / Ю.А. Суринов // Изв. СО АН СССР. - 1978. - № 8. - Вып. 2. - С. 106-125.
6. Суринов Ю.А. Современное состояние и перспективы развития теории переноса излучения и радиационного теплообмена [Текст] / Ю.А. Суринов // Тр. первой Российской национальной конференции по теплообмену: в 9 т. Т. 9. Радиационный и сложный теплообмен. - М.: Изд-во МЭИ, 1994. - С. 238-243.
7. Руднев Б. И. Анализ использования зональной модели радиационного теплообмена в камере сгорания дизельного двигателя [Текст] / Б.И. Руднев, О.В. Повалихина // Тр. пятой Российской национальной конференции по теплообмену: в 8 т. Т. 6. Интенсификация теплообмена. Радиационный и сложный теплообмен. - М.: Издательский дом МЭИ, 2010. - С. 261-263.
8. Кулешов О.Ю. Повышение вычислительной эффективности зонального метода расчета сложного теплообмена на основе методики коррекции базовых оптико-геометрических характеристик излучения [Текст] / О.Ю. Кулешов, В.М. Седелкин // Тр. пятой Российской национальной конференции по теплообмену: в 8 т. Т. 6. Интенсификация теплообмена. Радиационный и сложный теплообмен. - М.: Издательский дом МЭИ, 2010. - С. 227-230.
Сведения об авторе: Руднев Борис Иванович, доктор технических наук, профессор.
Y Y
Список литературы
