Numerical research of interconnected heat and mass transfer processes in the “two hot particles—polymeric material—air” system was executed. The joint effect of several local heat sources on the main integrated characteristic of ignition process (ignition delay time) was established. Two ignition models characterized by the relative positioning of hot particles on a polymeric material surface were revealed. Besides, there were established characteristics of local heat sources and the distance between them (700 K<Tp<1150 K and L>1.5 or Tp>1150 K and 0.25<L<1.5) when regularities of heat and mass transfer processes in the “two hot particles—polymeric material—air” system are similar to regularities of heat and mass transfer processes in the “single hot particle—polymeric material—air” system.
1. Introduction
In recent years, polymeric materials (polymethyl methacrylate, polystyrene, polyethylene, etc.) have been widely adopted in various industries as decorative and constructive elements. Polymeric material products are very susceptible to thermal effects [1–4] even at rather low outside temperature (T≈400–600K). Under conditions of some technological processes (at increased ambient temperature) strength characteristics are changed, melting occurs, dangerous carcinogens and are emitted. Possible temperatures of technological processes for power production can reach more than 1000 K. Under such conditions, the probability of local power sources (metal and nonmetallic particles warmed to high temperature with sizes about several millimeters) formation is high [5–9].
The numerical research results [10, 11] were obtained for thermal conduction and thermal convection processes during a polymeric material ignition by a single metal particle heated to high temperature. Established theoretical consequences can be used for developing guidelines and methods to reduce the flammability, ignition preventing, and subsequent stages of polymeric material combustion processes. However, in practice, several (two, three, etc.) small-size particles heated to high temperature can cause fires. Ignition conditions and heat transfer characteristics may be different for the “single hot particle—polymeric material—air” system and the “two hot particles—polymeric material—air” system. For example, it is known that ignition delay time of solid condensed substance (composite propellant) at high concentration (large number per unit area surfaces) of hot particles in the gas stream equals the values of ignition delay time at condensed substance heating by a massive plate with constant (during ignition period) high temperature [5, 6]. Therefore, more detailed information about characteristics of physical and chemical processes at polymeric material heating by several hot particles is necessary for the development of relevant precautionary activities.
The purpose of the present study was to develop the mathematical model and analyze the characteristics of interconnected heat and mass transfer processes during interaction of two small-size metal particles heated to high temperatures with a polymeric material at the accounting of thermal conduction in condensed substance and thermal decomposition of polymeric material, thermal convection, and diffusion and oxidation reaction in outside gas area.
2. Problem Statement
It was determined [12] that ignition conditions for liquid condensed substances are defined by the distance between two neighboring particles at the various quantities of hot particles (local heat sources) falling to the flammable material surface. Therefore, previous study results [12] were taken into account at the problem statement.
A scheme with two steel particles heated to high temperature in the parallelepiped shape (with the same sizes) situated on the surface of typical polymeric material, polymethyl methacrylate (PMMA), was chosen to simulate the conditions of hot particle flow interaction with polymeric material. We reasonably [12] used the “two hot particles—polymeric material—air” system (Figure 1) instead of the “several hot particles—polymeric material—air” system for heat and mass transfer process investigation. Areas with sizes xL and yL much larger than the sizes of hot particles xp and yp were allocated in polymeric material and air (Figure 1).
A scheme of the solution domain at t=0: 1: oxidizer (gas mixture at 0<t≤td); 2: hot particle; 3: polymeric material.
The rate of PMMA thermal decomposition process accelerated at polymeric material near-surface layer heating at0<t<td by the heat of hot particles. Gaseous products of polymeric material pyrolysis mixed with an oxidizer (air) at the diffusion in outside gas area. The gas mixture was warmed by the thermal convection at its movement along lateral sides (x=x1, x=x2, x=x3, x=x4, and y1<y<y2) of heat sources (Figure 1). The ignition occurred when the temperature and concentration of combustible component (PMMA gas) reached the critical values.
The assumptions for the problem statement of heat and mass transfer process include the following.
A gas substance with known kinetic parameters is formed as a result of PMMA thermal decomposition. The realization of only one “effective” oxidation reaction where one substance reacts was assumed.
A possible burning out of the polymeric material is not considered. It was found by the authors [12] that the burning out of substance near-surface layer has an insignificant effect on the ignition characteristics at local heating during short time period (less than 0.5 s).
Surfaces of hot particles and polymeric material (x1<x<x2, x3<x<x4, y=y1) have ideal thermal contact. The possibility of gas gap formation is not considered.
Ignition conditions were taken into account [13].
The heat release from the oxidation reaction of PMMA gas is more than the heat consumed by both hot particles to the heating of polymeric material and gas mixture.
The temperature of gas area in a zone of intense exothermic reaction exceeds the initial temperature of local heat sources.
3. Mathematical Model and Solution Method
The mathematical model describes the interconnected processes of thermal conduction and thermal decomposition in a condensed substance and thermal convection and diffusion processes and oxidation reaction in the gas area are represented by the system of nonstationary partial differential equations (0<t<td).
For gas mixture (0<x<x1, x2<x<x3, x4<x<xL, y1<y<y2; 0<x<xL, y2<y<yL), the following equations and conditions were implemented.
Poisson’s equation:(1)∂2ψ∂x2+∂2ψ∂y2=-ω.
The equation of gas mixture movement:(2)∂ω∂t+u∂ω∂x+ν∂ω∂y=υ1∂2ω∂x2+∂2ω∂y2+βg∂T1∂x.
The thermal convection equation:(3)ρ1C1∂T1∂t+u∂T1∂x+ν∂T1∂y=λ1∂2T1∂x2+∂2T1∂y2+Q1W1.
The diffusion equation:(4)ρ1∂Cf∂t+u∂Cf∂x+ν∂Cf∂y=ρ1D1∂2T1∂x2+∂2T1∂y2-W1.
The balance equation:(5)Cf+Co=1.
The heat balance equation for hot particles (x1<x<x2, x3<x<x4, y1<y<y2):(6)ρ2C2∂T2∂t=λ2∂2T2∂x2+∂2T2∂y2.
The heat balance equation for the polymeric material (0<x<xL, 0<y<y1):(7)ρ3C3∂T3∂t=λ3∂2T3∂x2+∂2T3∂y2-Q3W3.
The initial conditions (t=0):(8)T=T0at0<x<xL,0<y<y1;T=Tpatx1<x<x2,x3<x<x4,y1<y<y2;T=T0,Cf=0,ψ=0,ω=0at0<x<x1,x2<x<x3,x4<x<xL,y1<y<y2;0<x<xL,y2<y<yL.
The boundary conditions (0≤t≤td):(9)x=0,x=xL,0<y<y1:∂T3∂x=0;x=0,x=xL,y1<y<yL:∂T1∂x=0,∂Cf∂x=0,∂ψ∂x=0;x=x1,x=x3,y1<y<y2:-λ1∂T1∂x=-λ2∂T2∂x,T1=T2,∂Cf∂x=0,∂ψ∂x=0,ψ=0;x=x2,x=x4,y1<y<y2:-λ2∂T2∂x=-λ1∂T1∂x,T2=T1,∂Cf∂x=0,∂ψ∂x=0,ψ=0;y=0,0<x<xL:∂T3∂y=0;y=y1,0<x<x1,x2<x<x3,x4<x<xL:-λ3∂T3∂y=-λ1∂T1∂y+Q1W1,T3=T1,ρ1D1∂Cf∂y=-W3,∂ψ∂y=0;y=y1,x1<x<x2,x3<x<x4:-λ3∂T3∂y=-λ2∂T2∂y,T3=T2;y=y2,x1<x<x2,x3<x<x4:-λ2∂T2∂y=-λ1∂T1∂y,T2=T1,∂Cf∂y=0,∂ψ∂y=0,ψ=0;y=yL,0<x<xL:∂T1∂y=0,∂Cf∂y=0,∂ψ∂y=0.
Stream function ψ and vortex velocity vector ω:(10)u=∂ψ∂y,v=-∂ψ∂x,ω=rotyv→=∂v∂x-∂u∂y.
The mass rate of combustible gas mixture oxidation [14]:(11)W1=ρ1k10СonСfmexp-E1RT1,where n and m are constants (n=m=1).
The mass rate of polymeric material pyrolysis:(12)W3=∫0y1ρ3k30exp-E3RTsdy.
The diffusion coefficient of polymeric material thermal decomposition products in a gas mixture:(13)D1=D0T12731.7.
Volume fractions of gas mixture components were determined by its dimensionless mass concentrations:(14)φf=Cf/ρfCf/ρf+Co/ρo,φo+φf=1.
Thermophysical properties of gas mixture were defined as(15)λ1=λfφf+λoφo,C1=Cfφf+Coφo,ρ1=ρfφf+ρoφo.
The system of (1)–(7) with the corresponding initial and boundary conditions was solved by the finite difference method. The equations of elliptic type (Poisson’s and gas mixture movement) were solved by the alternating direction method. Difference analogues of heat balance and diffusion equations were solved by the locally one-dimensional method. A system of differential equations was solved by the iteration method and the sweep method at each iteration (for nonlinear equations), using the implicit four-point difference scheme.
In developing the algorithm for solving numerically the ignition problem, we used the elements of the algorithm developed for the numerical simulation of conjugate heat transfer processes at the local heating of an area with limited sizes [15]. To improve the accuracy of integrated characteristics computation, we used not less than 500 knots of difference grid for each coordinate and chose a short time step Δt=10-6 s. The verification of numerical research results was executed similar to [10, 11] by the test of conservation for the used difference scheme. The error of the energy conservation law in the solution area does not exceed 2.7% (Figure 1).
4. Results and Discussion
The numerical investigations were carried out for the following values of parameters [16–21]: the heat effect of oxidation reaction Q1=2.5×107 J/kg; the heat effect of PMMA thermal decomposition Q3=106 J/kg; activation energy E1=0.125×106 J/mol, E3=0.13×106 J/mol; preexponential factor k10=1010s-1, k30=2.82×109s-1; the thermal expansion coefficient β=0.0009K-1; the kinematic viscosity coefficient υ1=1.4×10-5 m2/s; the diffusion coefficient D0=8.12×10-6 m2/s; the initial temperature of air and polymeric material T0=300 K; temperature of PMMA pyrolysis beginning Ts=500 K; hot particles sizes xp=4×10-3 m, yp=2×10-3 m; solution area sizes xL=20×10-3 m, yL=20×10-3 m. Thermophysical properties of the substances are presented in Table 1 (Figure 1) [16–19].
Thermophysical properties of substances.
Substance
ρ, kg/m3
C, J/(kg⋅K)
λ, W/(m⋅K)
Air
1.161
1190
0.026
Steel
7831
470
49
PMMA (solid)
1200
1466.5
0.19
PMMA (gas)
1.29
1005.6
0.025
The purpose of interconnected heat and mass transfer process research at the ignition of polymeric material consisted of the investigation of the joint effect of several local heat sources on the main integrated characteristic: ignition delay time td (Figure 1). In previous studies [10, 11], for a single particle, it was established that the heat content of a local heat source (characterized by its initial temperature Tp) is the main factor determining td. It is also reasonable for two particles to carry out the analysis of the influence of the distance between heat sources Δx on ignition delay time (Figure 1). Therefore, the numerical investigation was executed by varying the initial temperature of heat sources between 700 K < Tp < 1500 K and varying the parameter of L (where L=Δx/xP) that characterizes the distance between two particles in the range of 0.25<L<2.
The dependence of polymeric material ignition delay time on the value of L at the initial temperature of heat sources Tp = 900 K is shown in Figure 2. It was found that in the “two hot particles—polymeric material—air” system the value of td rises for increasing L from 0.25 to 1.5. In case of L>1.5, the ignition delay time does not change (td=0.0936 s). It corresponds to the ignition period of PMMA at its interaction with a single hot particle (curve 1 in Figure 3). Thus, the value of L=1.5 is the limit (for Tp = 900 K) when the particles still have a joint effect on the intensity of heat and mass transfer processes in the system (Figure 1).
PMMA ignition delay time td versus L at Tp = 900 K.
PMMA ignition delay time td versus initial temperature of particle Tp: 1: one particle; 2: two particles at L=0.25.
It is seen from Figure 2 that the maximum change of ignition delay time in the “two hot particles—polymeric material—air” system is 6.7% at the parameters 0.25<L<1.5 and Tp = 900 K. The joint effect of several particles at td decreases for increasing the heat sources initial temperature. The dependence of PMMA ignition delay time at the initial temperature of a single particle (curve 1) and two particles (curve 2) at L=0.25 is shown in Figure 3. The values of ignition delay time at Tp > 1150 K do not differ from the cases of a single hot particle and two hot particles. It can be explained by the fact that for increasing Tp the heat content of particles rises. The influence of neighboring particles on the warming of polymeric material near-surface layer by thermal conduction and gas mixture by thermal convection decreases. In this case, the temperature of PMMA pyrolysis products increases and less time is required for warming the mixture of combustible gases and oxidizer. Ignition zones as in case of a single particle [10, 11] are formed near the lateral sides of the local heat sources and move to the “polymeric material—hot particle” border (x=x1, x=x2, x=x3, x=x4, and y→y1). Besides, it was determined that at the conditions of 700 K < Tp < 1150 K and L>1.5 or Tp>1150 K and 0.25<L<1.5 the characteristics of heat and mass transfer processes at the PMMA ignition by several hot particles are identical to values of td [10, 11] calculated at the PMMA ignition by a single hot particle.
The isotherms in the “two hot particles—polymeric material—air” system are shown in Figure 4 at the ignition moment for three different values of L. Opposite to the “single hot particle—polymeric material—air” system [10, 11] at variation of 700 K < Tp < 1150 K and 0.25<L<1.5, two ignition models are realized. It is characterized by localization of an oxidation reaction in the gas area relative to the surfaces of polymeric material and local heat sources. It is seen in Figure 4(a) that one local ignition zone is formed in case of the small (L=0.25) distance between particles. The ignition zone is located on the symmetry axis x=xL/2 of the solution area (Figure 1). The maximum temperature gradients and the largest concentration of combustible products (PMMA gas) occurred in comparison with areas x<x1 and x>x4 as a result of hot particles joint effect. For increasing L and other equal conditions, two ignition areas were formed near the internal lateral sides (x=x2, x=x3, and y1<y<y2) of hot particles (Figures 4(b) and 4(c)).
Isotherms (T,K) at ignition moment td=0.087 s, L=0.25 (a), td=0.092 s, L=0.5 (b), td=0.096 s, and L=1 (c) at Tp=900 K: 1: gas mixture; 2: hot steel particle; 3: PMMA.
It is seen from Figure 5 that two main convective whirlwinds intensifying heat and mass transfer processes formed in the gas area at various distances between hot particles.
Stream function (ψ, m2/s) at ignition moment td=0.087 s, L=0.25 (a), td=0.092 s, L=0.5 (b), td=0.096 s, and L=1 (c) at Tp=900 K.
It leads to the decrease of temperature and concentration of combustible gases near the external lateral sides (x=x1, x=x4, and y1<y<y2) of hot particles. Therefore, in these zones, gases-phase ignition does not occur. At the small distance between heat sources diffusion (Figure 5(a)), convection and heat sink from zones x2<x<x3 and y1<y<y2 is absent (Figure 1). Ignition occurs near the symmetry axis x=xL/2 in the zone with maximum concentration and temperature of gas mixture (Figure 1). For increasing the distance between hot particles (Figures 5(b) and 5(c)) in zones x2<x<x3 and y1<y<y2, the secondary convective whirlwinds forms. It leads to the decrease of temperature and concentration of gas mixture near the symmetry axis x=xL/2. In such conditions, ignition occurs near the internal lateral sides (x=x2, x=x3, and y1<y<y2) of heat sources.
5. Conclusions
The predictive mathematical model of interconnected heat and mass transfer processes at the ignition of polymeric material by several small-size hot metallic particles was developed. It considers thermal conduction and thermal decomposition in a condensed substance, thermal convection, and diffusion and oxidation reaction in the gas area.
Characteristics of a local heat sources and its relative location on a polymeric material surface (700 K < Tp < 1150 K and L>1.5 or Tp>1150 K and 0.25<L<1.5) were established when the regularities of heat and mass transfer processes in the “two hot particles—polymeric material—air” system are similar to regularities of heat and mass transfer processes in the “single hot particle—polymeric material—air” system.
The developed mathematical model can be used in mechanical engineering for a definition of the most fire danger parts of polymeric material constructional products during the process of its interaction with several local heat sources which are formed as a result of technological processes. Besides, the model can be used in chemistry for the definition of effective kinetic characteristics of polymeric material thermal decomposition and oxidation.
Nomenclatures and UnitsC:
Constant specific heat, J/(kg·K)
Cf:
Dimensionless concentration of combustible gases in the gas mixture
Co:
Dimensionless concentration of oxidant in the gas mixture
D0:
Coefficient of diffusion (at T=293 K), m2/s
D1:
Coefficient of diffusion, m2/s
E:
Activation energy, J/mol
g:
Gravitational acceleration, m/s2
k0:
Preexponential factor, s−1
L:
Dimensionless parameter, characterized the distance between two neighboring particles (L=Δx/xP)
Q1:
Heat effect of combustible gas mixture oxidation reaction, J/kg
Q3:
Heat effect of polymeric material thermal decomposition reaction, J/kg
R:
Ideal gas constant, J/(mol·K)
t:
Time, s
Δt:
Time step, s
td:
Ignition delay time, s
T:
Temperature, K
Tp:
Initial temperature of hot particles, K
Ts:
Temperature of the thermal decomposition beginning for polymeric material, K
T0:
Initial temperature of air and polymeric material, K
u, v:
Components of combustible gas velocity at the projection onto the axes x, y, m/s
W1:
Mass rate of combustible gas mixture oxidation reaction, kg/(m3·s)
W3:
Mass rate of polymeric material thermal decomposition reaction, kg/(m3·s)
x, y:
Cartesian coordinates, m
Δx:
Distance between two particles, m
xL, yL:
Solution domain sizes, m
xp, yp:
Hot particle sizes (xp=x2-x1,yp=y2-y1), m.
Greek Symbolsβ:
Coefficient of thermal expansion, K−1
λ:
Thermal conductivity, W/(m·K)
ρ:
Density, kg/m3
υ:
Coefficient of kinematic viscosity, m2/s
φf:
Volume fractions of PMMA gasification products
φo:
Volume fractions of air
ψ:
Stream function, m2/s
ω:
Vortex velocity vector, 1/s.
Subscripts1:
Gas mixture (air and PMMA gas)
2:
Hot particles
3:
Polymeric material.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgment
The reported study was partially supported by the Russian Science Foundation (no. 14-39-00003).
MeissnerK.PoltersdorfB.Thermomechanical load and property changes in processing of polymeric materials199010426527610.1002/adv.1990.0601004032-s2.0-0025577169TelnovA. V.ZavyalovN. V.KhokhlovY. A.SitnikovN. P.SmetaninM. L.TarantasovV. P.ShadrinD. N.ShorikovI. V.LiakumovichA. L.MiryasovaF. K.Radiation degradation of spent butyl rubbers2002633–62452482-s2.0-003612248710.1016/S0969-806X(01)00645-4BonatiA.MerusiF.PolaccoG.FilippiS.GiulianiF.Ignitability and thermal stability of asphalt binders and mastics for flexible pavements in highway tunnels20123766066810.1016/j.conbuildmat.2012.07.0962-s2.0-84865968011WangY.ZhangJ.Influences of specimen size and heating mode on the ignitability of polymeric materials in typical small-scale fire test conditions201236323124010.1002/fam.11032-s2.0-84858081867SummerfieldM.1963661Princeton UniversityGol'dshlegerU. I.BarzykinV. V.MerzhanovA. G.Mechanism and laws of ignition of condensed systems by a two-phase flow19747327728610.1007/BF007428092-s2.0-1642345366BurkinaR. S.MikovaE. A.High-temperature ignition of a reactive material by a hot inert particle with a finite heat reserve200945214415010.1007/s10573-009-0019-32-s2.0-64949112076KuznetsovG. V.StrizhakP. A.Transient heat and mass transfer at the ignition of vapor and gas mixture by a moving hot particle2010535-692393010.1016/j.ijheatmasstransfer.2009.11.0312-s2.0-72649101867GlushkovD. O.KuznetsovG. V.StrizhakP. A.Research of macroscopic regularities of heat and mass transfer at the ignition condition of a liquid high-energy material by an immersed source with a limited energy capacity20142014976453710.1155/2014/7645372-s2.0-84901251372GlushkovD. O.StrizhakP. A.Heat and mass transfer at ignition of solid condensed substance with relatively low calorific power by a local energy source2012211697710.1134/S18102328120100792-s2.0-84857568877GlushkovD. O.KuznetsovG. V.StrizhakP. A.Theoretical assessment of ignition stability of a typical polymeric material by a local power source20142321019KuznetsovG. V.StrizhakP. A.On the scale of “simultaneous” influence of several “hot” particles on the conditions of heat and mass transfer at ignition of liquid condensed substance200918426327010.1134/S18102328090400182-s2.0-70450239868VilyunovV. N.ZarkoV. E.1989Amsterdam, The NetherlandsElsevierFrank-KamenetskyD. A.1969New York, NY, USAPlenum PressKuznetsovG. V.SheremetM. A.Conjugate natural convection in an enclosure with local heat sources20091334136010.1615/ComputThermalScien.v1.i3.602-s2.0-78449273478BhattacharjeeS.KingM. D.PaoliniC.Structure of downward spreading flames: a comparison of numerical simulation, experimental results and a simplified parabolic theory200481233910.1088/1364-7830/8/1/0022-s2.0-2342419112WuK. K.FanW. F.ChenC. H.LiouT. M.PanI. J.Downward flame spread over a thick PMMA slab in an opposed flow environment: experiment and modeling2003132469770710.1016/S0010-2180(02)00520-52-s2.0-0037350452AyaniM. B.EsfahaniJ. A.SousaA. C. M.The effect of surface regression on the downward flame spread over a solid fuel in a quiescent ambient200711267862-s2.0-3444764787010.2298/TSCI0702067AVargaftikN. B.Moscow, RussiaStars2006 (Russian)DakkaS. M.TG/MS of poly(methyl methacrylate) the effect of heating rate on the rate of production of evolved gases20047537657722-s2.0-184245921810.1023/B:JTAN.0000027172.41943.55PetersonJ. D.VyazovkinS.WightC. A.Kinetic study of stabilizing effect of oxygen on thermal degradation of poly(methyl methacrylate)1999103388087809210.1021/jp991582d2-s2.0-0033365781