THESE de DOCTORAT - cerfacs

THESE de DOCTORAT - cerfacs THESE de DOCTORAT - cerfacs

22.10.2014 Views

28 Chapter 2: Computation of noise generated by combustion ρ Dh [ 1 Dt = ρc Dp p ρr Dt − T Dρ ρ Dt − T D ] Dt (ln r) ( ) + ρ ∑ h s,k + ∆h 0 DYk f ,k Dt k (2.23) The last term of Eq. (2.23) can be developed so that heat release due to combustion ˙ω T comes up. This term is defined as ˙ω T = − ∑ k (h s,k + ∆h 0 f ,k ) ˙ω k [74]. Using the balance of species (Eq. 2.6) results in ( ) ( ∑ h s,k + ∆h 0 f ,k ˙ω k − ∂J ) k ∂x k i = ∑ k h s,k ˙ω k + ∑ k ∆h 0 f ,k ˙ω ∂J k − ∑ h k s,k k ∂x i − ∑ k ∆h 0 ∂J k f ,k (2.24) ∂x i ∂J = − ˙ω T − ∑ h k k (2.25) ∂x k i Equation (2.8) is combined with Eqs. (2.23),(2.25) so that c p r Dp Dt − c pT Dρ Dt − ρc pT D Dt (ln r) − ˙ω T − ∑ k h k ∂J k ∂x i = Dp Dt − ∂q i ∂x i + τ ij ∂u i ∂x j + ˙Q + ρ f i u i (2.26) Dividing everywhere by (c p T) and neglecting volume forces ( f i = 0) yields ( Dp 1 Dt [ + 1 c p T rT − 1 c p T ) − Dρ Dt − ρ D (ln r) Dt ∂J − ˙ω T − ∑ h k k + ∂q i ∂u − i τ ij − ˙Q ∂x k i ∂x i ∂x j ] = 0 (2.27) A similar derivation for the material derivative of the density ρ can be found in [4]. The resulting expression is really useful to study acoustics, as both Dρ Dt and 1 Dp c 2 Dt appear explicitly. Reorganizing Eq. 2.27 and knowing that c p = rγ/(γ − 1) [ Dρ Dt = 1 Dp (γ − 1) c 2 + Dt c 2 ] ∂J − ˙ω T − ∑ h k k + ∂q i ∂u − i τ ij − ˙Q ∂x k i ∂x i ∂x j − ρ D (ln r) (2.28) Dt where the velocity of sound c = √ γrT has been introduced.

2.3 Hybrid computation of noise: Acoustic Analogies 29 2.3.2 Lighthill’s Analogy When considering frozen vortical structures or a frozen flame, clearly all the forces implied in each of these physical mechanisms are totally balanced and no sound is produced. But once unsteadiness occurs, these forces are not balanced anymore. In the case of turbulence, for example, the reciprocal motion of the surrounding gas tries to compensate for this ‘unexpected’ dynamics. Nevertheless, this reciprocal motion fails in completely absorbing the vortical energy and therefore local compressions and dilatations due to compressibility are transmitted to the adjacent particles: sound has been generated. In the source region (vortex/flame dynamics) the flow can be considered incompressible, forasmuch as the compressible part of the pressure fluctuation p ′ a is much smaller than the pressure fluctuations generated by both unsteady vorticity p ′ h (if sound is generated by turbulence) or unsteady temperature p′ f (if sound is generated by combustion). In this region, pressure fluctuations can be then computed by solving a Poisson equation ∇ 2 p ′ = ψ where ψ is a function of the implied sources. In the far-field, on the contrary, the fluctuations due to acoustics p ′ a are much more important than the significatively much smaller p ′ h ,p′ f . In the case of turbulence, for instance, it is known that hydrodynamic pressure fluctuations are negligible in the far field since they typically decay at least as the inverse third power of the distance to the sources [80]. Pressure fluctuations p ′ a may be assessed then by solving an acoustic wave equation. Lighthill’s analogy is considered as the starting points of aeroacoustics as a research field. The Navier-Stokes equations are rearranged so that the resulting formulation has the ‘shape’ of a inhomogeneous wave equation: in the left hand side (LHS) the acoustic operator and in the right hand side (RHS) the sources responsible for noise generation. Neglecting volume forces, the mass and momentum equations (Eqs. 2.4 and 2.5) can be expressed in their conservative forms: ∂ρ ∂t + ∂ρu j ∂x j = 0 (2.29) ∂ρu i ∂t + ∂ρu iu j ∂x j = − ∂p ∂x i + ∂τ ij ∂x j (2.30) to Eq. (2.30) and sub- Applying the time derivative ∂ ∂t stracting them results in to Eq. (2.29), the spatial derivative ∂ ∂x i ∂ 2 ρ ∂t 2 − ∂2 p ∂xi 2 = ∂2 ( ) ρui u j − τ ij ∂x i ∂x j (2.31) In order to construct the D’Alembertian operator, the term 1 c 2 ∞ ∂ 2 p ∂t 2 is added to both RHS and LHS of Eq. (2.31). c ∞ represents the sound velocity in the quiescent propagation media (far field). It

2.3 Hybrid computation of noise: Acoustic Analogies 29<br />

2.3.2 Lighthill’s Analogy<br />

When consi<strong>de</strong>ring frozen vortical structures or a frozen flame, clearly all the forces implied in<br />

each of these physical mechanisms are totally balanced and no sound is produced. But once<br />

unsteadiness occurs, these forces are not balanced anymore. In the case of turbulence, for example,<br />

the reciprocal motion of the surrounding gas tries to compensate for this ‘unexpected’<br />

dynamics. Nevertheless, this reciprocal motion fails in completely absorbing the vortical energy<br />

and therefore local compressions and dilatations due to compressibility are transmitted<br />

to the adjacent particles: sound has been generated. In the source region (vortex/flame dynamics)<br />

the flow can be consi<strong>de</strong>red incompressible, forasmuch as the compressible part of the<br />

pressure fluctuation p ′ a is much smaller than the pressure fluctuations generated by both unsteady<br />

vorticity p ′ h (if sound is generated by turbulence) or unsteady temperature p′ f<br />

(if sound<br />

is generated by combustion). In this region, pressure fluctuations can be then computed by<br />

solving a Poisson equation ∇ 2 p ′ = ψ where ψ is a function of the implied sources. In the<br />

far-field, on the contrary, the fluctuations due to acoustics p ′ a are much more important than<br />

the significatively much smaller p ′ h ,p′ f<br />

. In the case of turbulence, for instance, it is known that<br />

hydrodynamic pressure fluctuations are negligible in the far field since they typically <strong>de</strong>cay at<br />

least as the inverse third power of the distance to the sources [80]. Pressure fluctuations p ′ a may<br />

be assessed then by solving an acoustic wave equation.<br />

Lighthill’s analogy is consi<strong>de</strong>red as the starting points of aeroacoustics as a research field. The<br />

Navier-Stokes equations are rearranged so that the resulting formulation has the ‘shape’ of a<br />

inhomogeneous wave equation: in the left hand si<strong>de</strong> (LHS) the acoustic operator and in the<br />

right hand si<strong>de</strong> (RHS) the sources responsible for noise generation. Neglecting volume forces,<br />

the mass and momentum equations (Eqs. 2.4 and 2.5) can be expressed in their conservative<br />

forms:<br />

∂ρ<br />

∂t + ∂ρu j<br />

∂x j<br />

= 0 (2.29)<br />

∂ρu i<br />

∂t<br />

+ ∂ρu iu j<br />

∂x j<br />

= − ∂p<br />

∂x i<br />

+ ∂τ ij<br />

∂x j<br />

(2.30)<br />

to Eq. (2.30) and sub-<br />

Applying the time <strong>de</strong>rivative ∂ ∂t<br />

stracting them results in<br />

to Eq. (2.29), the spatial <strong>de</strong>rivative<br />

∂<br />

∂x i<br />

∂ 2 ρ<br />

∂t 2 − ∂2 p<br />

∂xi<br />

2 = ∂2 ( )<br />

ρui u j − τ ij<br />

∂x i ∂x j<br />

(2.31)<br />

In or<strong>de</strong>r to construct the D’Alembertian operator, the term 1<br />

c 2 ∞<br />

∂ 2 p<br />

∂t 2 is ad<strong>de</strong>d to both RHS and LHS<br />

of Eq. (2.31). c ∞ represents the sound velocity in the quiescent propagation media (far field). It

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

Saved successfully!

Ooh no, something went wrong!