CC BY
47
GASDYNAMIC ANALYSIS OF THE AIRBAG GAS GENERATOR WITH A POROUS CHARGE
*Tomsk State Pedagogical University "Tomsk State University
В статье представлена математическая модель работы airbag-газогенератора с пористым моноблочным зарядом. Выполнен сравнительный анализ результатов вычислений и экспериментальных данных. На основе анализа развития основных физических процессов предложена и обсуждена новая конструктивная схема airbag-газогенератора.
Mathematical model
In the present time car airbags, in which a source of pressure are low temperature solid propellant gas generators, have received wide spread occurrence and successfully develop. An autonomy, compactness, the large potential reserves of gases in small volumes of energetic material make by their indispensable at problem solving of car safety and in a number of extreme situations. In the most widespread design with a granular (peleet) charge the high temperature of gaseous combustion products and high-temperature particles of
a condensed phase in combustion products are main disadvantages causing difficulty of gas generator application. The given defects of the design are overcome by use of special units of cooling and mechanical filters.
Different perspective design largely deprived indicated disadvantages, are the airbaggas generators with porous solid propellant charge, which combustion will be organized in a wake filtration mode [1].
The design concepts of similar gas generators provide such organization of surface combustion of a porous charge, when the combustion products are essentially cooled, diffusing through a porous structure of a charge.
The filtration of the heat-carrier through a charge in a direction of reaction front results, in turn, in local heating of a combustion zone. The redistribution of energy between gaseous and condensed combustion products allows to get a low temperature gases. Besides the output of condensed combustion products is reduced due to their partial settling in a porous structure of a charge.
In the present paper the results of gas dynamic research of processes accompanying airbag gas generators with porous solid propellant charge (fig. 1) operation are submitted. The following scenario of physical processes development was considered. From the moment of igniter start-up high-temperature pyrotechnic combustion products go in the gas generator combustion chamber, causing a gradual warm-up and ignition of a solid charge, which porous structure is represented by a system connected among themselves propellant’s granules of the cylindrical shape.
and energy and for the variation in the free interior volume in the igniter then becomes [2]:
d_
dt
£
dt
(pW)=(l-z)M,~M„
{pRV/) = (1 - z)R!Mi-RM„
_(p GpW) = (]- z)CpjMj - CpM.,
dt
pW
J-l.
= (1 -z)CpJp,Mt
dW
Fig. 1, Gas generator with porous charge
The front «working» charge surface, connected to combustion, results in increase of pressure in the gas generator combustion chamber head volume; after nozzle gag breaking up the quasistationary process of end face burning begins. The combustion products, diffusing through a porous structure of a charge, cause preliminary (before the shave of combustion front) a warm-up of propellant, providing, in consequent, increased burning rate.
The mathematical model describing indicated processes, is founded on equations expressing laws of mass, momentum and energy conservation, which are formulated for igniter in averaged parameters, flow of combustion products through a porous charge - in onedimensional description of internal ballistic characteristics changes.
The mathematical model describing a change of gas dynamics parameters in a igniter case cavity, is represented by a system of ordinary differential equations for averaged on internal free volume W internal ballistic characteristics. It is supposed, that spherical granules of a pyrotechnic composition forming a igniter charge, flames simultaneously on all surface, and the design of the perforated case prevents take out of igniters grains. The system of equations for the conservation of mass
dm _ dm
where t - time; p, p - pressure and density; R,y- gas constant and adiabatic coefficient (ratio of specific heats); C, r - specific heat of gas components or mixture and isobaric combustion temperature; z - portion of the condensed phase in the combustion products; sub i relate to individual properties of pyrotechnic mix.
The gas flow, M„ through the perforation holes is calculated by the familiar quasi - stationary equations into the external (in the gas generator chamber) pressure [2]. The transient volume of the pyrotechnic grains depends on the flame depth (m0 ip0 >©> 0, mn, po -mass and powder density) and is a function of the grain shape. For a wide spectrum of shapes, the relationship between the transient volume of these grains, co, and their sizes is very simple
m,
P,
3 ^ s
ГЬ
j=i
where e, (j = 1, 2, 3) is the initial geometric grain characteristics.
The variation of the typical size of a grain as a result of the burning is described through the burning rate, w, and the depth, e, by the equations de/dt = u .
The mathematical model of the processes in the combustion chamber of a gas generator is based on the integral equations of gas dynamics derived from the laws of mass, momentum, and energy conservation. Additional special equations define the thermal and physical properties of the multicomponent combustion products from the igniter charge and the air that initially fills the chamber. This system of equations is applied to a reference volume, V, bounded by a closed surface comprised ofboth gas-permeable (A) and gas-impermeable (S) materials. Heat and mass are exchanged between gas flow and gas generator elements. The changes of mass, momentum and energy in the control volume are coupled by the volume interactions, the influence of the external
surroundings on each of the mentioned quantities, and their transfer through A. This system of conservation equations is written below in the integral form with a generalized coordinate system [2]:
JpdV + JpNdA = \MpdS,
at v v s
~ fpudv + JndA + Jims + Jpfas = o,
Ct y y t S
\EdV + \{E + p)NdA + fqdS = jCppTppMpdS,
° t V V s s
\pRdV + IpRNdA = \RpMpdS,
Ot y y S
j- jpCpdV + jpCpNdA = \CppMpdS,
V V .5
y-1 2 '
with u - vector of velocity; N - modulus of the velocity normal to the surface determined by the scalar product jo of the vector of the velocity « and the unit vector nor- 5
mal to the outer surface n, N = III, it); H - flow of CL impulse, Pn + piLV; M - mass supply density; H - enthalpy of flow products: /- friction stress; q - heat flux density; sub p relate to porous propellant charge.
Density of mass Income from combustion surface Mn, is determined by a ratio M.0 =up \ 1 + |, where
p - temperature coefficient of burning rate dependence from «initial» temperature of propellant. Local temperature of a charge T (x, t) is determined by average of a local temperature profile obtained from the solution of a heat conduction equation, describing change of a charge temperature owing to a filtration through it of combustion products. The particular values of surfaces sizes A and S are determined by a factor of a porosity of a charge.
The ignition time calculation is made in frameworks of the solid theory received broad development in activities of V. Vilyunov [3]. The mathematical model of an ignition problem of solids includes a heat conduction equation with allowance for of exothermal responses and equation of a chemical kinetics [3, 4].
The mathematical statement of a problem is closed by appropriate initial and boundary conditions. The systems of ordinary differential equations are decided by the predication - correction schemes; the spatial motion of gas is calculated by a modified Godunov method of any gap decay [5] with the increased order of approximating. For implementation of an ignition equation the adapted computational grid varying depending on conditions of a warm-up is used [4].
Calculation Results
The calculations were realized for the gas generator with a charge with the following main characteristics: Diameter of a charge - 90 mm; Length of a charge L = = (120-140) mm; Length of a charge together with the filter - 150 mm;
Characteristic of a material of a charge:
Density - 1920 kg/m3; Coefficients of propellant thermal conductivity - 6.3 x 10~7 m2/'s; Coefficients of propellant heat conduction - 0.86 W/(mxK); Porosity -(0.3-0.6); Temperature coefficient of the burning rate law (0.001-0.004).
On fig. 2 the results of comparison of computational and experimental data [1 ] for relation of pressure to time for area 1 of head volume and 2 - parameters of gas behind the filter (fig. 1) are submitted.
t. s
Fig, 2. Comparisons of computational and experimental data for relation «pressure - time»
The analysis of calculations results shows their good qualitative conformity to experimental data. Relation of piBssuiBp,{t) (head volume of the gas generator) has rather good quantitative conformity with experiment. The curve />,(/) in all time range is situated above than experimentally determined level. Some rise of pressure on the given curve during completion of a charge combustion is connected as to smaller losses of heat into propellant, so by heating of a filter to the given point of time. It is possible to explain more significant rise of the experimentally measured pressure by the factor of lowering of the filter area, free for pass of gas, owing to a settling in it of condensed particles.
On fig. 3 the pressure profiles on the combustion chamber in various points of time are shown. Qualitatively character of submitted relations is stipulated, mainly, motion of combustion front.
го
а.
0.0
0.2
0.4 0.6
X/L
0.8
1.0
Fig. 3. Pressure profile on length of a charge in various points of time
Character of gas distribution (combustion products) temperature is poor change of its profile with regard to coordinate of a currant position x, of combustion front during gas generator operating time. The particular kind of relation is adduced on fig, 4.
1 200
600
300
I \
\ \ \
\
ro
Q_
CL
<1
Fig. 5. Comparison of experimental and computational values of pressure differential along a charge
As a whole concurrence of the experimental and computational data is satisfactory, that allows to use the offered approach for designing of gas generators with a porous charge.
Based on discovered features of a gas dynamics of processes in the gas generator with end face burning charge it is possible to offer the following design of porous charges for airbag gas generators application. In perforated on an external surface (see fig. 6) a case the charge of the tubular form is placed, igniter is located in channel. The igniter provides ignition of a channel surface of a charge and further mode of cocurrent combustion.
Xf Xf + 5 mm Xf + 10 mm Xf + 15 mm
Fig. 4. A profile of gas temperature in a neighbourhood of combustion front
On fig. 5 the outcomes of a pressure differential size (between head volume and area behind the filter) comparison for charges of a various porosity r| are submitted. The painted area represents a field of a series of experimental results; a solid line - averaged (on operating time of the gas generator) computational values. It is visible, that the computational pressure differential practically is always less significant, that it is possible to explain by difference of actual «internal» geometrical parameters of a porous charge and used model and also absence of the registration of motion and settling of a combustion products condensed phase in a charge.
Fig. 8. Scheme of gas generator
For this design the solution was constructed similarly to foregone scheme - the averaged parameters in a channel of a charge and one-dimensional flow and heat conduction along a radius of a porous charge and filter (but for the axisymmetrical description of flows) were calculated. The geometry of a charge was selected so that the initial combustion surface corresponded to similar value for a previous calculations. On fig. 7 the time
history of averaged on a channel of a charge pressure is shown.
It is necessary' to note qualitative and quantitative difference of processes in comparison with the gas generator, where the porous charge bums on an end face surface. In particular, the relation of pressure to time (fig. 7) gr
has increasing character during all active operating time, that is caused both development of geometry of a combustion surface, and higher degree of a warm-up of a much more thin porous charge mass. Mean temperature of combustion products on an output of the device is higher on 74 K (in comparison with calculation for the gas generator with a end face charge) as combustions products which are flowing past through a rather thin porous charge, have no time to be cooled up to its temperature.
Thus, in the given paper the technique of the gas dynamic analysis of airbag-gas generators with a porous charge ensuring calculation of their main parameters is submitted. The legitimacy of the offered approach application is confirmed by the satisfactory matching with experimental data. The new design of airbag-gas
0 20 40 60 80 100 120
t, ms
Fig. 7. A time history of averaged on a charge channel pressure
generators with porous charges of surface combustion high-energy materials in a mode of a cocurrent filtration is offered; the key features of a gas dynamics of processes devices are ascertained.
This work was partly supported by CRDF,
References
1. ShandakovV,, PuzanovV., Komarov V., Borociikin V. The method of low temperature generating in solid gas generators//Combustion, Explosion and Shock Waves, 1999. № 4.
2. Vorozhtsov A., Bondarchuks., Salko A,, Kondratova 0. Mathematical Simulation of Airbag inflation by Low Temperature Gas Generator Products // An intern. J, Propellants, Explosives, Pyrotechnics, 2000. V. 25. issue 5.
3. Vilyunov V.N., Zarko V.E. Ignition of Solids, Elsevier, Amsterdam, 1989,
4. Borovskoy I.G., Bondarchuk 8,8., Kozlov BA, V.N. Vilyunov, Solution of the Problem of Condensed System Ignition Based on a Difference Scheme with Time - Dependent Grid // Combustion, Explosion and Shock Waves. 1986. № 3.
5. Numerical Solution of Multidimensional Problems of Gas Dynamics / Ed, S.K. Godunov. M,, 1982,
