MINISTERUL ADMINISTRAŢIEI ŞI INTERNELOR - IGSU
MINISTERUL ADMINISTRAŢIEI ŞI INTERNELOR - IGSU
MINISTERUL ADMINISTRAŢIEI ŞI INTERNELOR - IGSU
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
În care<br />
*<br />
s<br />
T este temperatura medie de suprafaţă calculată în intervalul de timp<br />
*<br />
* *<br />
considerat. Dacă se defineşte numărul lui Damkohler ca fiind ( / ) T E A D − =<br />
atunci avem:<br />
*<br />
∂ρ<br />
= −D<br />
*<br />
∂t<br />
* *<br />
* *<br />
( ρ − δ ) ⎢(<br />
− E / T )<br />
⎡<br />
⎢⎣<br />
⎛ T<br />
⎜<br />
⎝ T<br />
exp ay<br />
*<br />
ay<br />
*<br />
⎞⎤<br />
−1⎟⎥<br />
⎟<br />
⎠⎥⎦<br />
exp p<br />
(10)<br />
*<br />
∂ρ<br />
E , valoarea derivatei va fi foarte mică când<br />
*<br />
*<br />
Pentru valori ridicate ale lui<br />
∂t<br />
*<br />
*<br />
T este apropiat de temperatura de ardere. În aceste condiţii, se poate utiliza ρ = 1 ca<br />
soluţie aproximativă a densităţii în faza iniţială a arderii (faza de încălzire statică a<br />
suprafeţei).<br />
În aceste condiţii, temperatura de încălzire statică a suprafeţei se obţine din<br />
prima şi a cincia ecuaţie din (9) :<br />
⎪⎧<br />
*<br />
*<br />
⎛ ⎞<br />
⎛<br />
⎞⎪⎫<br />
* 1 y<br />
* 2 *<br />
Γy<br />
+ Γ t<br />
* y<br />
T = + ⎨ ⎜ ⎟ − ⋅ ⎜<br />
⎟<br />
static 1 erfc<br />
⎬<br />
⎪⎩<br />
⎜ ⎟<br />
e erfc<br />
⎜<br />
Γ t +<br />
Γ<br />
*<br />
* ⎟<br />
(11)<br />
⎝ 2 t ⎠<br />
⎝ 2 t ⎠⎪⎭<br />
Pentru stagiul următor (arderea iniţială), sistemul de ecuaţii se scrie din nou,<br />
introducând = + θ<br />
* *<br />
T T :<br />
static<br />
⎧<br />
2<br />
* *<br />
∂θ<br />
∂ θ 1 ∂ρ<br />
⎛ ∂T<br />
⎞<br />
static ∂θ<br />
⎪ = +<br />
⎟<br />
⎜ +<br />
* * 2 * * * *<br />
⎪∂t<br />
∂y<br />
ρ ∂y<br />
⎝ ∂y<br />
∂y<br />
⎠<br />
⎪ *<br />
⎪<br />
⎡⎛<br />
*<br />
*<br />
∂ρ<br />
⎞ ⎛<br />
⎞⎤<br />
* * ( )<br />
⎪<br />
⎢⎜<br />
E Tay<br />
= −D<br />
ρ − δ exp − ⎟ ⋅ ⎜ −1⎟⎥<br />
*<br />
∂<br />
⎪<br />
⎢<br />
⎜ * ⎟ ⎜ *<br />
t<br />
⎟<br />
⎣⎝<br />
Tay<br />
⎠ ⎝ Tstatic<br />
+ θ ⎠⎥⎦<br />
⎪<br />
*<br />
*<br />
∂θ<br />
* ( ) ( 1−<br />
ρ ) ⎛ ∂T<br />
⎞ static Γθ<br />
⎨−<br />
0,<br />
t =<br />
⎟<br />
⎜<br />
⎜−<br />
−<br />
*<br />
*<br />
*<br />
*<br />
⎪ ∂y<br />
ρ ⎝ ∂y<br />
⎠ ρ<br />
⎪ * *<br />
ρ ( y , 0)<br />
= 1<br />
⎪<br />
*<br />
*<br />
⎪θ<br />
( y , 0)<br />
= θ ( ∞,<br />
t ) = 0<br />
⎪<br />
⎪<br />
⎪<br />
⎪⎩<br />
În mod similar ecuaţiei (12) se procedează şi pentru celelalte faze ale arderii.<br />
Folosind forma adimensională pentru câteva valori de intrare impuse:<br />
* 11 *<br />
*<br />
A = 10 ; E = 40;<br />
δ = 0,<br />
3;<br />
Γ = 0,<br />
4<br />
(13)<br />
rezultă valorile numerice pentru temperatura, respectiv timpul de aprindere adimensionale:<br />
*<br />
*<br />
T ay = 1,<br />
5;<br />
tay<br />
= 0,<br />
27<br />
(14)<br />
În mod similar se pot obţine şi celelalte mărimi de calcul din sistemul de<br />
ecuaţii diferenţiale (12), având în vedere că acestea se rezolvă cuplat.<br />
În concluzie, se poate spune că modelul de ardere tip Atreya, reprezintă<br />
precursorul şi stă la baza tuturor modelelor de ardere mai avansate, care de altfel au<br />
fost şi implementate în programele de calcul comerciale de referinţă (Exemplu: FDS<br />
– Fire Dynamic Simulator dezvoltat de NIST, şi nu numai). Prezentarea detaliată a<br />
modelului fizic şi a sistemului de ecuaţii aferent a reprezentat un imbold pentru<br />
crearea propriului program de calcul numeric la Facultatea de Pompieri.<br />
(12)<br />
169