Linearized Navier-Stokes solver for in-duct aeroacoustics
GUNILLA EFRAIMSSON Department of Aeronautical and Vehicle Engineering KTH, Stockholm, Sweden [email protected]
Axel Kierkegaard, Wei Na, Mats Åbom, Ciarán O’Reilly, Susann Boij, Hans Bodén Outline
• Tango-aspects
• Linearized Navier-Stokes equations solver
• Two tecases - Area expansion - Thick orifice
• Conclusions and future work Tango setting
Q˙ 0 Unsteady heat release
Generation of sound waves Perturbation flow & flame
Pressure Oscillation p0
Study the coupling sound , flow and heat release using linearized Navier-Stokes equations Generation and Propagation of Sound
• Generation Fans, rotor-stator interaction, constrictions, jets, combustion - Acoustic analogies, direct methods - Transient phenomena è transient simulations (LES,, DES, uRANS) or scaling laws
• Propagation Free air, pipes, influence of flow, influence of solid surfaces (shielding), coupling of sound and flow (whistling, damping of sound), duct terminations - Linear equations in time or frequency domain: Helmholtz equation, wave equation, linearized Euler equations etc
Propagation of sound in ducts with constrictions
Straight cylindrical duct with orifice plate Turbulent meanflow Incoming waves are partly reflected and partly transmitted Reflections from downstream obstacles (or boundary conditions)
Zone I Zone III Zone IV
Zone II x 0 1 = x2 = 0 Propagation of sound in ducts with constrictions
Straight cylindrical duct with orifice plate Turbulent meanflow Incoming waves are partly reflected and partly transmitted Reflections from downstream obstacles (or boundary conditions)
Zone I Zone III Zone IV
Zone II x 0 1 = x2 = 0 Propagation of sound in ducts with constrictions
Straight cylindrical duct with orifice plate Turbulent meanflow Incoming waves are partly reflected and partly transmitted Reflections from downstream obstacles (or boundary conditions)
Zone I Zone III Zone IV
Zone II x 0 1 = x2 = 0 Propagation of sound in ducts with constrictions
Straight cylindrical duct with orifice plate Turbulent meanflow Incoming waves are partly reflected and partly transmitted Reflections from downstream obstacles (or boundary conditions)
Zone I Zone III Zone IV
Zone II x 0 1 = x2 = 0 Propagation of sound in ducts with constrictions
Straight cylindrical duct with orifice plate Turbulent meanflow Incoming waves are partly reflected and partly transmitted Reflections from downstream obstacles (or boundary conditions)
Zone I Zone III Zone IV
Zone II x 0 1 = x2 = 0 Propagation of sound in ducts with constrictions
Straight cylindrical duct with orifice plate Turbulent meanflow Incoming waves are partly reflected and partly transmitted Reflections from downstream obstacles (or boundary conditions)
Zone I Zone III Zone IV
Zone II x 0 1 = x2 = 0 Efficient simulation methodology
The methodology should be efficient for
• Confined flows
• Mean flows effects
• Dissipation due to acoustic-flow interaction
Available simulation methods Linearized Navier-Stokes equations in frequency domain
Assumptions • Linear acoustics • Frequency domain e-iωt • Isentropy. Temperature dependence via Sutherland’s formula and ideal gas law.
• Solved using Comsol Multiphysics (FEM, direct solver)
• Previous work on orifices and area expansions in 2D with good results. Now extended to 3D. Kierkegaard et al JASA 2010, JSV 2011, 2011
Inclusion of • Viscous effects e.g. acoustic boundary layers in the vicinity of edges
Linearized Navier-Stokes equations in frequency domain Linearized Navier-Stokes equations in frequency domain
Two stages
1. Calculate mean flow (using RANS)
2. Calculate perturbations around mean flow using linearized Navier-Stokes equations
Frequency domain è One calculation for every frequency in a frequency sweep
Case 1 Area expansion
Experiment by Ronneberger (1987) Theory by Boij & Nilsson (2003).
Step 1: Mean flow
• Incompressible, steady-state, non-linear Reynolds Averaged Navier-Stokes (RANS). • Realizable k-ε turbulence model • Flow speed 14.5 m/s (Mach number = 0.08) Step 2: Linearized N-S eqns
• Perturbations around mean flow. • Single frequency calculations. • Frequency sweep 100 < f < 3000 Hz, ∆f = 100 Hz.
• Finite Element Method (FEM), COMSOL Multiphysics. Density perturbations
• Frequency 1000 Hz • Both sound waves and vorticity
• ρ’(x,t) = Re{ρ(x, ω)exp(-iωt)}
Density perturbations
• Frequency 500 Hz
• Frequency 2500 Hz Side a Side b Scattering matrix
(ω)
• Relates outgoing waves to incident waves.
• Every matrix element is a funtion of frequency.
Scattering matrix, magnitude
Experiments Simulations Analytical solution Scattering matrix, phase
Experiments Simulations Analytical solution Case 2 Thick orifice Mean flow • M = 0.042, RANS k-ε Calulate perturbations
• Density perturpations • f = 1800 Hz (maximum amplification) Scattering matrix (S)
Solid black line: Simulations
Dashed black: Experiments at KTH
Solid grey: Experiments at Le Mans, France (Testud et al 2009) Whistling potentiality Amplicifaction or dissipation of sound Energy power balance è Sign of Eigenvalue of (1-SHS)
Solid black line: Simulations Linearized Navier-Stokes equations Dashed black line: Simulations Linearized Euler equations Dotted black: Experiments at KTH Solid grey: Experiments at Le Mans, France (Testud et al 2009) Conclusions
• Frequency domain linearized Navier-Stokes can predict acoustic scattering in ducts, such as area expansions, orifices, impedance tubes
• Acoustic boundary layers needed in case of sound-vortex interaction
• No special treatment at edges (eg Rutta conditions)
• Fast compared to time domain in 2 space dimensions
• Arbitrary geometries possible
• Influence of mean flow different on different scattering matrix components, different flow field Present and future work
• Inclusion of energy equations done
• Include thermo-acoustic sources in process
• Evaluation for higher order modes in process
• Iterative solver
Thank you for your attention