CC BY
53
Вычислительные технологии
Том 7, № 1, 2002
ЧИСЛЕННОЕ МОДЕЛИРОВАНИЕ ГАЗОДИНАМИЧЕСКИХ ПРОЦЕССОВ В КАМЕРЕ СГОРАНИЯ АВТОМОБИЛЬНОГО УСТРОЙСТВА БЕЗОПАСНОСТИ (AIRBAG)
А. Д. Рычков, Н.Ю. ШокинА Институт вычислительных технологий СО РАН, Новосибирск, Россия e-mail: rych@ict.nsc.ru, nina.shokina@ict.nsc.ru
The paper is devoted to the numerical modeling of gas-dynamic processes in airbag combustion chamber. The physical model of the processes and the numerical realization of the model are considered. Some calculation results for the particular devices are compared with full-scale tests.
Введение
В настоящее время в системах индивидуальной безопасности на автомобильном транспорте широкое распространение получило устройство эйрбэг (airbag), представляющее собой камеру сгорания, наполненную гранулами твердого унитарного топлива с относительно невысокой температурой горения и соединенную со специальной оболочкой из газонепроницаемой ткани. В исходном состоянии оболочка свернута в компактный рулон. При соударении автомобиля с препятствием срабатывает система воспламенения твердого топлива, продукты сгорания которого за 50 - 80 мс наполняют оболочку, превращая ее в эластичную подушку безопасности, предохраняющую водителя или пассажира от серьезных травм. Такими системами безопасности в настоящее время снабжены практически все автомобили зарубежного производства среднего класса и выше. Производители автомобилей уделяют самое серьезное внимание совершенствованию этих систем. Кратковременность и сложность физико-химических процессов, протекающих в эйрбэгах, заставляют обращаться к численному моделированию, которое на сегодняшний день является эффективным инструментом исследования подобных процессов. Численное моделирование позволяет совершенствовать конструкции эйрбэгов и повышать степень безопасности на автомобильном транспорте.
В работе рассматриваются физическая модель процессов, протекающих в камере сгорания эйрбэга, и ее численная реализация. Приводятся некоторые результаты расчетов для конкретных устройств и их сравнение с результатами натурных испытаний.
1. Физическая модель
Камера сгорания эйрбэга (рис. 1) представляет собой цилиндр, в центральной части которого располагается воспламенитель (бустер), состоящий из гранул твердого топлива
© А.Д. Рычков, Н.Ю. Шокина, 2002.
сферической формы. Продукты сгорания этих гранул содержат твердые частицы с достаточно высокой концентрацией, что способствует интенсификации процесса воспламенения основной массы топливных гранул.
Топливные гранулы, заполняющие остальной объем камеры сгорания, имеют цилиндрическую форму. Состав гранул подбирается таким образом, чтобы, с одной стороны, обеспечить высокую массовую скорость горения, а с другой — иметь низкую температуру продуктов сгорания и обеспечить отсутствие в них высоких концентраций опасных для здоровья веществ. Продукты сгорания через улавливатель твердых частиц (лабиринт) покидают камеру сгорания, проходя через ряд отверстий (сопел), которые с достаточной степенью точности можно моделировать кольцевым зазором.
Сформулируем основные допущения, позволяющие провести численное моделирование процесса работы такого устройства с учетом сложности протекающих процессов и невозможности их детального описания:
1. Течение является нестационарным и осесимметричным. Моделирование проводится в рамках континуальной модели: все основные компоненты (бустер, топливные гранулы, продукты сгорания) рассматриваются как три сплошные взаимопроникающие среды со своими скоростями и температурами. Между средами осуществляется взаимный обмен массой, импульсом и энергией.
2. Форма топливных гранул предполагается сферической, отклонения ее реальной цилиндрической формы от сферы учитываются коэффициентом формы в законах сопротивления. В процессе работы устройства гранулы топлива и бустера считаются неподвижными. Число гранул в единице объема (счетная концентрация) всегда остается постоянным и определяется из начальных условий их загрузки в камеру сгорания. Предполагается, что температура внутри гранул при их прогреве вплоть до воспламенения распределена по некоторому закону, полученному из приближенного решения задачи о прогреве сферы [1].
3. Используется модель горения с постоянной температурой поверхности [2]. Предполагается, что воспламенение гранулы происходит, когда температура ее поверхности достигает значения некоторой эмпирически задаваемой величины, называемой температурой воспламенения. Скорость горения также определяется эмпирической формулой, учитывающей ее зависимость от величины относительного давления в камере сгорания. Теплота, выделяемая при горении гранулы, расходуется на нагревание продуктов сгорания, и часть ее через механизм теплообмена возвращается грануле.
4. Продукты сгорания (несущий газ) являются совершенным газом с постоянным отношением удельных теплоемкостей. Течение несущего газа предполагается ламинарным
Рис. 1. Схема камеры сгорания эйрбэга.
и описывается уравнениями Навье — Стокса. В уравнении энергии пренебрегают работой сил трения и сил давления ввиду невысоких значений скоростей течения. Пренебрегают также кондуктивной теплопроводностью между горящими гранулами.
5. Инициирование воспламенения гранул бустера происходит вследствие поступления высокотемпературных продуктов сгорания воспламенителя (черный порох) через верхнюю границу области, занятой бустером. Продукты сгорания гранул бустера — газ и твердые мелкодисперсные частицы. Размеры частиц достаточно малы, поэтому движение такой двухфазной среды можно считать равновесным.
2. Математическая модель
Система уравнений, описывающая движение двухфазной среды с высокой плотностью загрузки частиц, имеет следующий вид:
дрг дрхщ 1 д
-ж + — + -гдГ (гр1щ)_ 3 + (1 - а) 3'
дргд , дргд Щ 1 д
-Ж + — + гот (гргд_ аЗь - Лосргд'
(р2 _ т .
(Ъ 3Р + Лосргд ,
((рз _ (Ъ _
др1и1 др\Щ, 1 д др д
п. + —^--+ "ТТ" (трхПхVI) + £1— _ —
дъ дх г дг дх дх
ди1 дх
1д
+ - — ( Г^1£1
г дг
ди1 дь1 дг дх
-р2в12Щ - рзв13Щ,
др1ь1 др1и1у1 1 д
дъ
+
дх
1 д др д ( + - — (Гр1ь1) + £1— _ — ^1£1 г дг дг дх \
дь1 ди1 дх дг
-р2@12У1 - рз@13У1,
др1Н1 + др1Н1и1 + 1 д ( ь ) д дъ дх г дг дх
£1
Хд дН1
СР дх
1 д_
г дг
Хд дЬ,1
Ср дг
+р2112 (Тз2 - Т1) + рзЪз (ТзЗ - Т1) + Зд<^д +
р _ рд ВТ1, р1 _ £1рд ,р2 _ £2рр, рз _ £3рь, £1 + £2 + £3 _1,Ы _ Срд Т .
(1)
(2)
(3)
(4)
(5)
1 д ( ду1
+ - — 2г^1£1— -
г дг \ дг
(6)
(7)
(Т (Т
р2С2-(2 _ р2112 (Т1 - Т32) + ЛосргдСгд (Т1 - Т3), рзСз-^з _ рзЪ2 (Т1 - Т3з) , (8)
(9)
Здесь индексы 1, д относятся к несущему газу, 2 — к топливным гранулам, 3 — к бустеру; г _ 2, 3; Ъ — время; и, V — компоненты вектора скорости по осям Ох и Ог цилиндрической системы координат; рд, рР, рь, Срд, С2, Сз — истинные (физические) плотности несущего газа, топливных и бустерных гранул и их удельные теплоемкости соответственно; ргд, а, Сгд — счетная плотность твердой фазы в продуктах сгорания бустера (масса в единице объема камеры), ее доля в этих продуктах и удельная теплоемкость; Х1 _ £1Хд — коэффициент теплопроводности; в^, — межфазные коэффициенты сопротивления и теплообмена; Тзр, Т3ь — температуры поверхности гранул топлива и бустера. Для их
определения используются результаты приближенного решения задачи о прогреве сферы
7 - Т? 1 + 0.2Б1 1 ;'
а? ¿? ^ ^ 4а? г _ 0.1Б12
— — критерий Био; Го = —2 ? ¿2
где Б1 = ? ? — критерий Био; Го = —--критерий Фурье; ^ — -;-^,
Л? а,- 3 + 1.2Б1 + 0.13Б1
а? = 0 ; индексы 3 = 2, 3 относятся к гранулам топлива и бустера соответственно. Р/Ч?'
Коэффициент сопротивления в? вычисляется при £д < 0.8 по формуле Эргана:
в? = 15»Г?^ + 1-75, 3 = 2,3,
а при £д > 0.8 как коэффициент сопротивления сферы, которая эквивалентна реальной грануле по массе:
в =3 с Ра |ид - и? I £—2.65 • = 23
в? = 4с°? а?р?р0 £1 ' 3 = 2,3' 24
С™ = — (1 + 0.15Кв0.68^ при Ке? < 103 Ке? 4 ? 7
сд? = 0.44 при Ке? > 103,
Ке?
|ид - и?1 Рд £д
Мд
где и = (и, V) — вектор скорости; р? — коэффициент формы; 0? — эквивалентный диаметр сферы.
Коэффициенты теплообмена определяются по формулам
Ас = аос|и!| 4 0^2,
71. = 6^? Ке? рг1/3ф(£ рг= С^ . = 23
714 = а2р0 ' ^ = 4 ' Рг= Ла ' 3 = 2,3'
Ф (£) = 0.508 - 0.56 (1 - £д) при £д < 0.4, ф(£) = 1 - 1.164(1 - £д)2/3 при £д > 0.4,
где аос — эмпирический коэффициент, учитывающий отклонение траекторий частиц твердой фазы при обтекании гранул несущим газом.
Изменение размеров горящих гранул 0?, а также массоприходы 3? с их поверхностей и скорости горения г? находятся по формулам
= -2г?, 3? = р0£?г?п?, г? = а? + Ь? (р/ро)ъ' , 3 = 2, 3, (10)
где Б? = — площадь поверхности гранулы; а?, Ь?, 7? — эмпирические константы; р0 — атмосферное давление.
Для системы уравнений (1)-(9) задаются следующие граничные условия: 1) Р1 и1 = Огд, = 0, 7 = Тгд при г < ггд, х = хк, 0 < г < г1, где Огд, Тгд — массовый приход продуктов сгорания воспламенителя и их температура, ггд — время работы воспламенителя;
dTi
2) u1 = v1 = 0, —— = 0 — условия прилипания и отсутствия теплообмена на непрони-
dx
цаемых стенках;
3) "мягкие" граничные условия в выходном сечении r = rk, 0 < x < x1.
3. Разностная схема
Для численного решения системы уравнений (1) - (9) с граничными условиями 1-3 используется разностная противопотоковая ЬИ-схема второго порядка точности по пространству и времени, обладающая ТУБ-свойствами. Эта схема близка к схеме, реализованной в работе [3]. Особое внимание уделяется методике расчета течения несущего газа в окрестности разрыва плотности р2 (х, г) на границе х2 < х < х3, г = г3. Здесь используются условия на скачке, отделяющем область течения двухфазной смеси от области течения чистого газа, полученные в работе [4]:
[Р1 Мп1] = 0, [Срд Т + ^/2] = 0, [р1Иг 1] = 0.
Эти условия отражают сохранение массы и энергии на скачке, где ип1, ит 1 — нормальная и касательная проекции вектора скорости газа на поверхности разрыва. Тогда остальные параметры течения при дополнительном задании одного из них определяются из условий на разрыве:
[Р = j2
1 , Ч 1 1
_(1 - е2)-1 __
1а- pg+_
■ 1 ■ 1 1
.P1. pg+ pg-
(1 _ £2) 1 , j = PlUnl-
Величины со знаком "+" берутся в области двухфазной среды, со знаком "—" — в области чистого газа. Так, при ип1 > 0 за скачком задается значение давления полученное из разностной схемы, а при ип1 < 0 — соответственно р_. После этого значения всех параметров течения по обе стороны скачка в каждом узле разностной сетки находятся из решения системы нелинейных уравнений с помощью метода Ньютона. Такой подход позволяет полностью избежать появления паразитных нефизичных осцилляций численного решения, возникающих при сквозном счете окрестности разрыва величины р1.
4. Некоторые результаты расчетов
Были проведены расчеты процессов в камере сгорания реального эйрбэга при условиях, соответствующих натурным стендовым испытаниям. Общая масса топливных и бустерных гранул составляла 60 и 2 г соответственно. Топливные гранулы имели форму цилиндров диаметром 6 мм и высотой 4 мм. Диаметр сферических гранул бустера составлял 0.1 мм. Тепло, выделяемое при их горении, равнялось Qp = 3591 Дж/г и Q¿ = 3223 Дж/г соответственно.
На рис. 2 приведены кривые изменения давления в зависимости от времени в верхней части камеры сгорания, полученные в расчете (кривая pcaic) и в эксперименте (кривая Рехр). Видно, что даже с учетом всех упрощающих предположений результаты оказываются достаточно близкими. Более высокий "заброс" давления вследствие воспламенения и горения бустера связан, по-видимому, с использованием упрощенной модели горения и с неточными значениями коэффициентов в формулах (10) и в дальнейшем может быть скорректирован.
Р
50
40
30
20
10
20
40
60
ВО
Рис. 2. Кривые изменения давления.
Структура поля течения в камере сгорания в момент времени I =10 мс приведена на рис. 3-6. На рис. 3. показаны траектории движения "частиц" несущего газа, на рис. 4 изображены изолинии чисел Маха, на рис. 5 — изолинии температуры, а на рис. 6 — изолинии степени выгорания топливных гранул.
Как и ожидалось, структура течения является наиболее сложной в области лабиринта, где имеют место отрывы потока и застойные зоны. Интенсивное искривление траекторий "частиц" позволяет судить об эффективности осаждения в лабиринте несгоревших мелкодисперсных остатков топливных гранул, уносимых потоком несущего газа из камеры сгорания. Анализ изолиний чисел Маха показывает, что в выходном сечении сопла имеет место сверхзвуковая скорость истечения при наличии достаточно "толстого" пограничного слоя.
Распределение температуры в камере сгорания оказывается практически однородным. Некоторое возрастание температуры на периферии связано с продолжающимся горением топливных гранул после того, как в центральной части камеры оно закончилось. При этом, как следует из расположения изолиний степени выгорания топлива, горение гранул происходит параллельными слоями, что является несомненным достоинством данной газодинамической схемы эйрбэга.
Заключение
1. Разработана математическая модель газодинамических процессов в камере сгорания автомобильного устройства безопасности или эйрбэга, учитывающая процессы воспламенения и горения гранул твердого топлива.
2. Сравнение результатов расчетов и экспериментов свидетельствует о вполне удовлетворительном качестве предложенной математической модели, способной достаточно адекватно описывать основные процессы, протекающие в подобных устройствах.
3. Разработанная модель может быть легко обобщена на случай пространственного течения в эйрбэге.
Рис. 3. Траектории движения "частиц" несущего газа.
Рис. 4. Изолинии чисел Маха.
Рис. 5. Изолинии температуры.
Рис. 6. Изолинии степени выгорания топливных гранул.
Список литературы
[1] Жуков М.Ф., СолонЕнко О. П. Высокотемпературные запыленные струи в процессах обработки порошковых материалов / Отв. ред. В. Е. Накоряков. Новосибирск: ИТ СО АН СССР, 1990. 516 с.
[2] Теория горения порохов и взрывчатых веществ / Отв. ред. Я. Б. Зельдович, А. Ф. Беляев. М.: Наука, 1982. 244 с.
[3] Chen K.H., Shuen J.S. Three-dimensional coupled implicit methods for spray combustion flows at all speeds // AIAA 94-3047. P. 1-19.
[4] СтЕРНИН Л. Е. Основы газодинамики двухфазных течений в соплах. М.: Машиностроение, 1974. 212 с.
Поступила в редакцию 7 декабря 2001 г.
