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 107<br />

of diffusion, this <strong>de</strong>gree of advancement of the combustion differs from the mass<br />

fraction Y (x) of the burnt gases that <strong>de</strong>fines the mixture composition at point x,<br />

that is, the fraction Y (x) of burnt gases that would be obtained by analyzing the<br />

composition of a gas sample taken from the said point.<br />

The simplifying assumptions previously stated do not limit the extent of the<br />

study that follows, which is of a qualitative nature. On the other hand, these assumptions<br />

simplify calculations.<br />

The wave equations are obtained by particularizing the general equations of<br />

continuity, momentum and energy 1 to the case of a one-dimensional stationary motion<br />

(∂/∂y = ∂/∂z = ∂/∂t = 0). When this is done, and in addition to the preceding<br />

assumptions are taken into account, the following system of equations is obtained:<br />

a) Continuity equation.<br />

that is<br />

d (ρv)<br />

dx<br />

ρv = m,<br />

= 0, (5.1)<br />

(5.1.a)<br />

where m is the mass flow normal to the wave, per unit surface.<br />

b) Continuity equation for the burnt gases. Since the mass flow per unit surface is<br />

ρv and the burnt fraction is ε, the mass flow of the burnt gases at section x is<br />

ρvε. Its variation with x is due to chemical reactions. The equation that gives<br />

this variation is<br />

d (ρvε)<br />

= w, (5.2)<br />

dx<br />

where w is the reaction rate. In virtue of (5.1.a), equation (5.2) can also be<br />

written in the form<br />

c) Momentum equation.<br />

m dε<br />

dx = w.<br />

ρv dv<br />

dx = − dp<br />

dx + 4 3<br />

(5.2.a)<br />

(<br />

d<br />

µ dv )<br />

. (5.3)<br />

dx dx<br />

This equation can be immediately integrated by taking into account (5.1.a), giving<br />

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

dx = i,<br />

(5.3.a)<br />

where i is an integration constant that must be <strong>de</strong>termined by the boundary conditions.<br />

2<br />

1 See chapter 3.<br />

2 In equations (5.3) and (5.4) it has been assumed that the second viscosity coefficient of the mixture<br />

(see chapter 2) is zero. When not, it is sufficient to substitute µ for `µ + 3 4 µ′´ in (5.3) and (5.4).

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

Saved successfully!

Ooh no, something went wrong!