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.

3 Experiments on Droplet Array <strong>Combustion</strong><br />

stream <strong>of</strong> <strong>the</strong> probe orifices in <strong>the</strong> plane that is flush with <strong>the</strong> combustion<br />

chamber ceiling). It is driven by <strong>the</strong> pressure gradient between <strong>the</strong> evacuated<br />

sample cylinder and <strong>the</strong> isobaric combustion chamber. To provide <strong>the</strong><br />

necessary input for this boundary condition, <strong>the</strong> total time <strong>of</strong> <strong>the</strong> gas sampling<br />

process ∆t sampling was empirically determined by hardware tests, but <strong>the</strong><br />

remaining parameters could be derived from gas dynamics. Choked flow is<br />

assumed downstream <strong>of</strong> <strong>the</strong> physical probe orifices, where <strong>the</strong> four probes<br />

are fitted to one piping joint. The fluid properties in <strong>the</strong> combustion chamber<br />

are expected to be homogeneous everywhere and are equal to <strong>the</strong> stagnation<br />

properties, denoted by <strong>the</strong> subscript zero [224, 396]. Thus, sonic (critical)<br />

conditions are reached at <strong>the</strong> “throat” <strong>of</strong> <strong>the</strong> piping. They are denoted by an<br />

asterisk. Employing <strong>the</strong> ratio <strong>of</strong> specific heats <strong>of</strong> a gas,<br />

κ≡ c p<br />

c v<br />

, (3.14)<br />

where<br />

c p = κ<br />

κ−1 R, c v = 1<br />

κ−1 R, and R = c p− c v , (3.15)<br />

critical temperature T ∗ and critical pressure p ∗ can be calculated via <strong>the</strong> isentropic<br />

relations <strong>of</strong> Equations (3.16) and (3.17), respectively. As a first estimate,<br />

κ= 1.4 is applied here, as generally valid for diatomic gases.<br />

T ∗<br />

= 2<br />

T 0 κ+1<br />

p ∗ ( 2<br />

=<br />

p 0 κ+1<br />

) κ<br />

κ−1<br />

(3.16)<br />

(3.17)<br />

Introducing <strong>the</strong> definition <strong>of</strong> <strong>the</strong> speed <strong>of</strong> sound for a perfect gas, <strong>the</strong> onedimensional<br />

continuity equation for steady flows gives:<br />

ṁ ∗ = p ∗ A ∗( κ<br />

) 1<br />

2<br />

. (3.18)<br />

R T ∗<br />

Everything but <strong>the</strong> throat area A ∗ is given, which is derivable from <strong>the</strong> total<br />

sampling time ∆t sampling and <strong>the</strong> sampling time under critical conditions ∆t ∗<br />

using a goal seek function. The time ∆t ∗ , in turn, can be obtained from <strong>the</strong> total<br />

mass collected during critical flow (Eq. (3.19)) and <strong>the</strong> sample temperature<br />

112

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

Saved successfully!

Ooh no, something went wrong!