Manfred Kaltenbacher in cooperation with A. Hppe, I. Sim - - PowerPoint PPT Presentation

manfred kaltenbacher in cooperation with a h ppe i sim
SMART_READER_LITE
LIVE PREVIEW

Manfred Kaltenbacher in cooperation with A. Hppe, I. Sim - - PowerPoint PPT Presentation

Spectral finite elements for a mixed formulation in computational acoustics taking flow effects into account Manfred Kaltenbacher in cooperation with A. Hppe, I. Sim (University of Klagenfurt), G. Cohen and S. Imperial (INRIA, Paris) and B.


slide-1
SLIDE 1

Wave Propagation and Scattering, Inverse Problems and Applications in Energy and the Environment, RICAM, 2011

Spectral finite elements for a mixed formulation in computational acoustics taking flow effects into account

Manfred Kaltenbacher in cooperation with

  • A. Hüppe, I. Sim (University of Klagenfurt), G. Cohen and S. Imperial

(INRIA, Paris) and B. Wohlmuth (TU Munich)

Alps-Adriatic University of Klagenfurt, Austria

slide-2
SLIDE 2

Wave Propagation and Scattering, Inverse Problems and Applications in Energy and the Environment, RICAM, 2011

Overview

 Physical modeling

 Pierce equation  Acoustic perturbation equation

 FE formulation (no flow)

 Acoustic conservation equations

  • Mixed formulation
  • Spectral elements

 Comparison to wave equation with pFEM

 FE formulation (with flow)

 Acoustic perturbation equations  Occurring instabilities  Stabilization (flux term and dissipative term)  Application to aeroacoustics

 Multi-Model approach

slide-3
SLIDE 3

Wave Propagation and Scattering, Inverse Problems and Applications in Energy and the Environment, RICAM, 2011

Acoustics in Flowing Media

 Euler’s equations  Idea of decomposition

Mean quantities Alternating quantities (disturbances)

slide-4
SLIDE 4

Wave Propagation and Scattering, Inverse Problems and Applications in Energy and the Environment, RICAM, 2011

Acoustics in Flowing Media

 Pierce equation (just for simple flows)  PML in time domain (Imbo Sim, Poster on Monday)

slide-5
SLIDE 5

Wave Propagation and Scattering, Inverse Problems and Applications in Energy and the Environment, RICAM, 2011

Acoustics in Flowing Media

 Acoustic perturbation equations1

 Subset of linearized Euler equations  Support just

  • acoustic modes
  • no entropy and vorticity modes

 Fully considers

  • convection
  • refraction
  • 1R. Ewert and W. Schröder. Acoustic perturbation equations based on flow decomposition via source filtering.

Journal of Computational Physics, 2003

slide-6
SLIDE 6

Wave Propagation and Scattering, Inverse Problems and Applications in Energy and the Environment, RICAM, 2011

Acoustics (no flow)

 Conservation equations  Linear acoustic wave equation Investigated Methods

  • h-FEM → mesh refinement
  • p-FEM & s-FEM → increase order of approximation

Wave Equation Acoustic Conservation Equations

slide-7
SLIDE 7

Wave Propagation and Scattering, Inverse Problems and Applications in Energy and the Environment, RICAM, 2011

Acoustics (no flow)

 Mixed formulation for conservation equations

 Discrete spaces1

  • 1G. Cohen & S. Fauqueux, Mixed Finite Elements with Mass-Lumping for the Transient Wave Equation

Journal of Computational Acoustics, 2000

Piola transform Lagrange polynomial space of order N Mapping:

slide-8
SLIDE 8

Wave Propagation and Scattering, Inverse Problems and Applications in Energy and the Environment, RICAM, 2011

Acoustics (no flow)

 Properties of Piola transform

 Preserves the normal component!  Term with gradient  Term with divergence

slide-9
SLIDE 9

Wave Propagation and Scattering, Inverse Problems and Applications in Energy and the Environment, RICAM, 2011

Acoustics (no flow)

 Spectral finite elements

slide-10
SLIDE 10

Wave Propagation and Scattering, Inverse Problems and Applications in Energy and the Environment, RICAM, 2011

Acoustics (no flow)

 Consequences of the Choice of Spaces

 Elements  Semidiscrete Galerkin formulation

slide-11
SLIDE 11

Wave Propagation and Scattering, Inverse Problems and Applications in Energy and the Environment, RICAM, 2011

Acoustics (no flow)

 Example

 Excitation with a sine pulse of main wavelength λ  Reference solution obtained with h= λ/120 and Δt= 1/(f 200)  Computational mesh with mean element size of λ/5 and Δt<=λ /(2*c)  Time Stepping:

  • h-FEM & p-FEM: Implicit Newmark scheme
  • s-FEM: Explicit leapfrog time stepping

Setup Defomed mesh

slide-12
SLIDE 12

Wave Propagation and Scattering, Inverse Problems and Applications in Energy and the Environment, RICAM, 2011

Acoustics (no flow)

 Comparison

Pressure Field

slide-13
SLIDE 13

Wave Propagation and Scattering, Inverse Problems and Applications in Energy and the Environment, RICAM, 2011

Acoustics (no flow)

 Comparison for time domain computations

Conservation equations

slide-14
SLIDE 14

Wave Propagation and Scattering, Inverse Problems and Applications in Energy and the Environment, RICAM, 2011

Acoustics Perturbation Equation

 Formulation  Spaces

Piola transform

slide-15
SLIDE 15

Wave Propagation and Scattering, Inverse Problems and Applications in Energy and the Environment, RICAM, 2011

Acoustics Perturbation Equation

 Semidiscrete Galerkin formulation  Example

 Initial condition  Flow

slide-16
SLIDE 16

Wave Propagation and Scattering, Inverse Problems and Applications in Energy and the Environment, RICAM, 2011

Acoustics in Flowing Media

 Results (cartesian grid)

slide-17
SLIDE 17

Wave Propagation and Scattering, Inverse Problems and Applications in Energy and the Environment, RICAM, 2011

Acoustics in Flowing Media

 Results: Long time simulation (Cartesian grid)

slide-18
SLIDE 18

Wave Propagation and Scattering, Inverse Problems and Applications in Energy and the Environment, RICAM, 2011

Acoustics Perturbation Equation

 Formulation  Spaces

Piola transform

slide-19
SLIDE 19

Wave Propagation and Scattering, Inverse Problems and Applications in Energy and the Environment, RICAM, 2011

Acoustics in Flowing Media

 Stabilization

 Central flux term  Reverse integration by parts on

slide-20
SLIDE 20

Wave Propagation and Scattering, Inverse Problems and Applications in Energy and the Environment, RICAM, 2011

Acoustics in Flowing Media

 Stabilization

 Averaging leads to  Add penalty (dissipative) term

slide-21
SLIDE 21

Wave Propagation and Scattering, Inverse Problems and Applications in Energy and the Environment, RICAM, 2011

Acoustics in Flowing Media

 Results: Long time simulation (Cartesian grid, penalty term)

slide-22
SLIDE 22

Wave Propagation and Scattering, Inverse Problems and Applications in Energy and the Environment, RICAM, 2011

Acoustics in Flowing Media

 Results: Long time simulation (Cartesian grid, penalty + flux term)

slide-23
SLIDE 23

Wave Propagation and Scattering, Inverse Problems and Applications in Energy and the Environment, RICAM, 2011

Acoustics in Flowing Media

 Results: Long time simulation (deformed grid)

slide-24
SLIDE 24

Wave Propagation and Scattering, Inverse Problems and Applications in Energy and the Environment, RICAM, 2011

Acoustics in Flowing Media

 Results: Long time simulation (deformed grid, penalty term)

slide-25
SLIDE 25

Wave Propagation and Scattering, Inverse Problems and Applications in Energy and the Environment, RICAM, 2011

Acoustics in Flowing Media

 Results: Long time simulation (deformed grid, penalty + flux term)

slide-26
SLIDE 26

Wave Propagation and Scattering, Inverse Problems and Applications in Energy and the Environment, RICAM, 2011

Acoustics in Flowing Media

 Results: Long time simulation (deformed grid)

Spurious waves

slide-27
SLIDE 27

Wave Propagation and Scattering, Inverse Problems and Applications in Energy and the Environment, RICAM, 2011

Acoustics in Flowing Media

 Shear flow

slide-28
SLIDE 28

Wave Propagation and Scattering, Inverse Problems and Applications in Energy and the Environment, RICAM, 2011

CAA (Computational Aeroacoustics)

 Air foil

 URANS CFD computations (Fluent, Michele Degenaro, AIT, Vienna)  Mach number about 0.3

slide-29
SLIDE 29

Wave Propagation and Scattering, Inverse Problems and Applications in Energy and the Environment, RICAM, 2011

CAA (Computational Aeroacoustics)

 Acoustic sources

 Lighthill analogy  Acoustic Perturbation equation (APE) RHS of Lighthill’s equation Test function RHS of APE Lamb vector

slide-30
SLIDE 30

Wave Propagation and Scattering, Inverse Problems and Applications in Energy and the Environment, RICAM, 2011

CAA (Computational Aeroacoustics)

 Arbitrary flow

Without flow With flow

slide-31
SLIDE 31

Wave Propagation and Scattering, Inverse Problems and Applications in Energy and the Environment, RICAM, 2011

Multi-Model Approach

 General idea

Acoustic perturbation equation Pierce equation PML layer Non-matching grid interface (Mortar framework)

slide-32
SLIDE 32

Wave Propagation and Scattering, Inverse Problems and Applications in Energy and the Environment, RICAM, 2011

Multi-Model Approach

 Interface conditions

 Continuity of pressure  Continuity of normal component of particle velocity Lagrange multiplier

slide-33
SLIDE 33

Wave Propagation and Scattering, Inverse Problems and Applications in Energy and the Environment, RICAM, 2011

Multi-Model Approach

 Formulation

 Acoustic perturbation equation  Pierce equation  Continuity of pressure in a weak sense

slide-34
SLIDE 34

Wave Propagation and Scattering, Inverse Problems and Applications in Energy and the Environment, RICAM, 2011

The End

Thank you for your attention!