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.

7.2. TURBULENT COMBUSTION THEORIES 213<br />

BURNED<br />

σ<br />

l<br />

UNBURNED<br />

σ<br />

GASES<br />

GASES<br />

Figure 7.4: Structure of the turbulent flame in the case l ≫ d l and v ′ u l , according to<br />

Shelkin.<br />

the laminar flame front thus increasing its surface as shown in Fig. 7.4. Let σ l be the<br />

effective surface of combustion and σ the surface of access to the flame of the unburnt<br />

gases. Velocity u t is given in this case by expression<br />

u t<br />

u l<br />

= σ l<br />

σ . (7.13)<br />

The problem lies then in estimating the value for σ l , for which several mo<strong>de</strong>ls have<br />

been proposed.<br />

Shelkin assumes that σ l is formed by a system of cones whose base is proportional<br />

to the square of the scale of turbulence l 2 and whose height is proportional to<br />

distance v ′ l/u l travelled by a mass of gas with a diameter l due to turbulent oscillations<br />

during the time l/u l nee<strong>de</strong>d by the laminar flame to cross this mass. Here, σ l /σ<br />

is the ratio of the lateral to the base area of the cones, and from (7.13), it results<br />

√<br />

( )<br />

u t<br />

v<br />

′ 2<br />

= 1 + k , (7.14)<br />

u l u l<br />

where k is a numerical coefficient of the or<strong>de</strong>r of magnitu<strong>de</strong> unity. Therefore, u t is<br />

also in<strong>de</strong>pen<strong>de</strong>nt from the scale of turbulence. When turbulence is weak, Eq. (7.14)<br />

shows that its effect is of the second or<strong>de</strong>r, whilst for very intense turbulence Eq. (7.14)<br />

reduces to (7.11). Through a different analysis M. Tucker [15] obtains the following<br />

expression for the case of weak turbulence<br />

( )<br />

u t<br />

v<br />

′ 2<br />

= 1 + k(θ) , (7.15)<br />

u l u l<br />

which is valid if v ′ /u l ≪ l. In this formula, k(θ) is a coefficient <strong>de</strong>pending on ratio<br />

θ of temperature T f of the burnt gases to temperature T 0 of unburnt gases. It is seen<br />

that when v ′ /u l ≪ l. Eqs. (7.14) and (7.15) agree.<br />

Karlovitz [16] reasons as follows. Let ¯X =<br />

√ ¯X2 be the root mean square<br />

displacement of a particle due to turbulent oscillations. ¯X is a increasing function of

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

Saved successfully!

Ooh no, something went wrong!