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

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

12.4. SIMPLIFIED EQUATIONS 293<br />

As for Y 2 , its value is given by expressions<br />

Interior region: Y 2 = 1 − Y 1 . (12.18)<br />

Exterior region: Y 2 = 1 − Y 3 . (12.19)<br />

Diffusion velocities ¯v d1 and ¯v d3 are given by Fick’s law 7<br />

Y 1¯v d1 = −D 12 ∇Y 1 , (12.20)<br />

Y 3¯v d3 = −D 23 ∇Y 3 . (12.21)<br />

On surface Σ f<br />

Y 1 = Y 3 = 0; Y 2 = 1. (12.22)<br />

Furthermore, since Y 1 and Y 3 must diffuse towards the flame in the stoichiometric<br />

ratio from (12.11), (12.20) and (12.21) we have on Σ f<br />

νD 12<br />

∂Y 1<br />

∂¯n i<br />

= D 23<br />

∂Y 3<br />

∂¯n e<br />

. (12.23)<br />

Here ν is the ratio between the masses of oxidizer and fuel nee<strong>de</strong>d for complete combustion.<br />

¯n e and ¯n i are normals to Σ f towards the exterior and interior regions respectively.<br />

System of equations (12.16) and (12.17) can be substituted by a single equation<br />

for a new variable Y <strong>de</strong>fined as follows<br />

⎧<br />

⎪⎨ Y 1 in the interior region,<br />

Y =<br />

⎪⎩ − 1 ν Y 3 in the exterior region.<br />

Y must satisfy the following differential equation<br />

(12.24)<br />

ρ(¯v · ∇)Y − ∇ · (ρD∇Ȳ ) = 0, (12.25)<br />

where D takes the value<br />

⎧<br />

⎨<br />

D =<br />

⎩<br />

D 12<br />

D 23<br />

in the interior region,<br />

in the exterior region.<br />

(12.26)<br />

Furthermore, Y is zero on the flame and takes values of opposite sign at both si<strong>de</strong>s of<br />

it. The <strong>de</strong>rivatives of Y at both si<strong>de</strong>s of the flame in normal direction to it must satisfy<br />

condition<br />

7 See Chap. 2.<br />

D 12<br />

∂Y<br />

∂¯n i<br />

= −D 23<br />

∂Y<br />

∂¯n e<br />

. (12.27)

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

Saved successfully!

Ooh no, something went wrong!