THESE de DOCTORAT - cerfacs
THESE de DOCTORAT - cerfacs THESE de DOCTORAT - cerfacs
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
- Page 1: UNIVERSITE MONTPELLIER II SCIENCES
- Page 5: An expert is a man who has made all
- Page 8 and 9: Y gracias familia mía!!, que así
- Page 11: Abstract Today, much of the current
- Page 14 and 15: 4 CONTENTS 3.4.1 Preconditioning .
- Page 16 and 17: List of Figures 1.1 Main sources of
- Page 18 and 19: 8 LIST OF FIGURES 5.19 Acoustic ene
- Page 20: List of Tables 1.1 EU recommended r
- Page 24 and 25: 14 Chapter 1: General Introduction
- Page 26 and 27: 16 Chapter 1: General Introduction
- Page 28 and 29: 18 Chapter 1: General Introduction
- Page 30 and 31: 2 Computation of noise generated by
- Page 32 and 33: 22 Chapter 2: Computation of noise
- Page 34 and 35: 24 Chapter 2: Computation of noise
- Page 36 and 37: 26 Chapter 2: Computation of noise
- Page 40 and 41: 30 Chapter 2: Computation of noise
- Page 42 and 43: 32 Chapter 2: Computation of noise
- Page 44 and 45: 34 Chapter 2: Computation of noise
- Page 46 and 47: 36 Chapter 2: Computation of noise
- Page 48 and 49: 38 Chapter 2: Computation of noise
- Page 50 and 51: 40 Chapter 2: Computation of noise
- Page 52 and 53: 42 Chapter 3: Development of a nume
- Page 54 and 55: 44 Chapter 3: Development of a nume
- Page 56 and 57: 46 Chapter 3: Development of a nume
- Page 58 and 59: 48 Chapter 3: Development of a nume
- Page 60 and 61: 50 Chapter 3: Development of a nume
- Page 62 and 63: 52 Chapter 3: Development of a nume
- Page 64 and 65: 54 Chapter 3: Development of a nume
- Page 66 and 67: 56 Chapter 3: Development of a nume
- Page 68 and 69: 4 Validation of the acoustic code A
- Page 70 and 71: 60 Chapter 4: Validation of the aco
- Page 72 and 73: 62 Chapter 4: Validation of the aco
- Page 74 and 75: 64 Chapter 4: Validation of the aco
- Page 76 and 77: 66 Chapter 4: Validation of the aco
- Page 78 and 79: 68 Chapter 4: Validation of the aco
- Page 80 and 81: 70 Chapter 4: Validation of the aco
- Page 82 and 83: 5 Assessment of combustion noise in
- Page 84 and 85: 74 Chapter 5: Assessment of combust
- Page 86 and 87: 76 Chapter 5: Assessment of combust
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