14.09.2014 Views

On the Formation of Nitrogen Oxides During the Combustion of ...

On the Formation of Nitrogen Oxides During the Combustion of ...

On the Formation of Nitrogen Oxides During the Combustion of ...

SHOW MORE
SHOW LESS

Create successful ePaper yourself

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

4.5 Simulation <strong>of</strong> Single Droplets<br />

<strong>the</strong> Dufour effect, Annamalai and Puri [19] state that heat flux is dominated<br />

by ordinary heat conduction. Consequently, <strong>the</strong> diffusion velocity in <strong>the</strong> gas<br />

phase g simplifies to<br />

∆u m,r =− D m ∂X g ,m<br />

. (4.47)<br />

X g ,m ∂r<br />

This is <strong>the</strong> best first-order approximation to <strong>the</strong> exact formulation 5 [147, 180,<br />

341] and avoids evaluating a linear system <strong>of</strong> equations [325]. According to<br />

Poinsot and Veynante [341] and Annamalai and Puri [19] <strong>the</strong> diffusion coefficient<br />

<strong>of</strong> a single species m in a multi-component gas can be approximated<br />

by Equation (4.48). In this case, <strong>the</strong> coefficient D m is not a binary diffusion<br />

coefficient but an equivalent diffusion coefficient <strong>of</strong> species m into <strong>the</strong> rest <strong>of</strong><br />

<strong>the</strong> mixture.<br />

D m = 1−Y m<br />

∑n;n≠m X n<br />

D mn<br />

(4.48)<br />

Besides, as <strong>the</strong> mass flow <strong>of</strong> <strong>the</strong> vaporized liquid ṁ(t ) is equal to <strong>the</strong> mass flow<br />

through a spherical shell with radius r , <strong>the</strong> radial velocity u r is proportional to<br />

1/r 2 [297, 298].<br />

The mass balance (Eq. (4.3)) and species conservation equation (Eq. (4.4)) remain<br />

unaffected by those simplifications for <strong>the</strong> gas phase, apart from <strong>the</strong> outlined<br />

changes in <strong>the</strong> diffusion velocity ∆u m,r . The viscous terms in <strong>the</strong> radial<br />

momentum equation (Eq. (4.5)) cancel each o<strong>the</strong>r out [363, 364]. This is a<br />

meaningful fact, as Euler flow (neglecting viscous effects), in <strong>the</strong> strict sense,<br />

is only valid here for constant vaporization mass flow ṁ, constant density ρ,<br />

and constant viscosity η g . Fur<strong>the</strong>rmore, a low Mach-Number is presumed.<br />

This finally reduces <strong>the</strong> momentum equation (Eq. (4.5)) to a single term, still<br />

accounting for Stefan flow [19, 443]. In <strong>the</strong> energy equation (Eq. (4.6)), viscous<br />

terms are also neglected and Euler flow is presumed at a constant pressure<br />

p g in space. The heat flux ˙q g ,r is evaluated neglecting <strong>the</strong> so-called Dufour<br />

term, <strong>the</strong> second term on <strong>the</strong> right hand side <strong>of</strong> Equation (4.13). Altoge<strong>the</strong>r,<br />

<strong>the</strong> governing equations for <strong>the</strong> gas phase as applied in <strong>the</strong> model are given<br />

below.<br />

5 In-depth analyses <strong>of</strong> both second-order effects (Soret and Dufour) can be found in literature [328, 335, 361].<br />

137

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

Saved successfully!

Ooh no, something went wrong!