21.11.2014 Views

Untitled - Aerobib - Universidad Politécnica de Madrid

Untitled - Aerobib - Universidad Politécnica de Madrid

Untitled - Aerobib - Universidad Politécnica de Madrid

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

5.2. WAVE EQUATIONS 109<br />

2) Burnt gases,<br />

x → +∞ : p → p f , ρ → ρ f , T → T f , v → v f , Y → 1, ε → 1. (5.8)<br />

Such conditions imply, naturally, the following ones,<br />

x → ±∞ :<br />

dp<br />

dx → 0, dρ<br />

dx → 0, dT<br />

dx → 0, dY<br />

dx → 0,<br />

dε<br />

→ 0. (5.9)<br />

dx<br />

Herein, the possibility of satisfying the previous boundary conditions, will not<br />

be discussed since it is subjected to a careful study in the chapter <strong>de</strong>dicated to <strong>de</strong>flagrations.<br />

Therein, it will be seen that the boundary conditions relative to the cold<br />

boundary of the flame give rise to a problem which, so far, has not been satisfactorily<br />

solved.<br />

In the present study it will be assumed that the boundary conditions can be<br />

satisfied and, therefore, that the problem is completely <strong>de</strong>termined.<br />

By taking the boundary conditions (5.7), (5.8) and (5.9) into the system of<br />

equations (5.1.a), (5.4.a)) and (5.4.a), the following laws of conservation are obtained,<br />

which agree with the invariants <strong>de</strong>duced in the preceding chapter<br />

ρ 0 v 0 = ρ f v f , (5.10)<br />

p 0 + ρ 0 v 2 0 = p f + ρ f v 2 f , (5.11)<br />

1<br />

2 v2 0 + c p T 0 + q = 1 2 v2 f + c p T f . (5.12)<br />

The boundary conditions relative, for instance, to the burnt gases, make possible<br />

the elimination of constants i, e and f from equations (5.3.a), (5.4.a) and (5.5.a),<br />

thus obtaining<br />

p + mv − 4 3 µ dv<br />

dx = p f + mv f , (5.13)<br />

1<br />

)<br />

m(<br />

2 v2 + c p T − qε − λ dT<br />

dx − 4 dv<br />

( 1<br />

)<br />

µv<br />

3 dx = m 2 v2 f + c p T f − q , (5.14)<br />

ρD dY<br />

− Y + ε = 0. (5.15)<br />

m dx<br />

The discussion and analysis of the possible solutions of the previous system are<br />

rather complicated. An example of a discussion for a similar system (but in which the<br />

diffusion neglected) can be found in Friedrichs’s work [3]. As a result of his analysis,<br />

Friedrichs conclu<strong>de</strong>s that except in the case where the reaction rate w is exceptionally<br />

high (and takes a well <strong>de</strong>fined value) weak <strong>de</strong>tonations are impossible. Consequently,<br />

the only possible <strong>de</strong>tonations with a normal reaction rate are strong <strong>de</strong>tonations and

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!