Shock/boundary layer interactions

Transcription

1 Shock/boundary layer interactions Turbulent compressible channel flows F.S. Godeferd Laboratoire de Mécanique des Fluides et d Acoustique Ecole Centrale de Lyon, France Journée Calcul Intensif en Rhône Alpes 1

2 Contents Presentation of the UFAST project Numerical method Shock in a curved channel Shock reflection Perspectives 2

3 An example of normal shock/boundary layer interaction [Experiment at the University of Cambridge - Bruce & Babinsky] 3

4 An example of normal shock/boundary layer interaction [Experiment at the University of Cambridge - Bruce & Babinsky] 3-a

5 The UFAST project Unsteady effects of shock wave induced SeparaTion European STREP project 18 academic and industrial partners Experiments RANS, URANS LES Control 4

6 Motivation of the study unsteady shock wave boundary layer interaction Aeronautical industry shock waves on wings/profiles, nozzle flows and inlet flows Interaction unsteadiness initiated and/or generated by SWBLI; often destabilized by the outer flow field; response of shock wave and separation to periodic excitations Control methods: synthetic jets, electro-hydrodynamic actuators, stream-wise vortex generators and transpiration flow Can we reproduce the unsteady interaction with URANS? Need of costly LES? 5

7 LMFA: URANS using platform elsa (ONERA) Unsteady Reynolds-Averaged Navier-Stokes Simulations Physical model mean flow Reynolds-averaged Navier-Stokes equations turbulence two-equation model k-l (turbulent kinetic energy, mixing length); or one-equation model Spalart-Allmaras for turbulent viscosity ν t. ( Sutherland viscosity law ν(t) = ν 0 T 0 +C T+C Adiabatic walls ) 3/2 T T 0 6

8 Numerical method Conservative finite volume method Roe fluxes with limiters Implicit Euler timestepping Multi-blocks structured mesh Parallel resolution 7

9 Initialization - Example: Channel flow with bump [expe. in Queens U. Belfast] Uniform initialization Euler laminar Navier-Stokes turbulent Navier-Stokes Subsonic freestream: Ma = Peak: Ma = Normal shock Animation 8

10 Performance Location Architecture # proc. # cores cpu/pt/ite Mem/pt Static max seconds bytes speedup LMFA Opteron 280 2Ghz P2CHPD Opteron Ghz x a IDRIS Nec SX ÉCL α-server EV7 1.15Ghz

11 URANS of curved channel flow at Ma = 1.45 Polish Institute of Mechanics experiment Inlet conditions: P=101kPa; T=290K. Turbulence: Tu=1%; L=1% channel height. Outlet conditions: pressure ratio specified in URANS to match experimental shock location 10

12 Splitted geometry with N = 8 procs; grid points Bottom wall boundary layer resolution: y + 2 Side walls: y + 20 Need to test the dependence of the solution on various elements of the simulation 11

13 Details of flow in shock zone The simulation allows a detailed investigation of the flow Mach contours, showing Ma = 1.45 upstream the shock 12

14 Oil flow visualization URANS streaklines Detachment/re-attachment length 13

15 URANS streamlines access to 3D flow structure 14

16 Adjustment of shock position: dependence on geometry 1 degree opening of channel at the outlet Mach number contours in a transverse wall for choked and unchoked geometries. The curve shows the Mach number at mid-section in this plane. 15

17 Dependence on expression of fluxes and outlet pressure Effect of counterpressure. Effect of fluxes scheme: Roe or Jameson. 16

18 Flow unsteadiness Main shock oscillations in the experiment: power spectrum No unsteadiness in the URANS need to introduce explicit fluctuations in URANS 17

19 Shock reflection case: expe. TU Delft 18

20 19

21 URANS Turbulence model: k L Mach number at two times and fluctuations at point in shock region 20

22 Unphysical solution! Huge detachment and unsteadiness due to poorly performing turbulence model. Reason: turbulence is strongly damped in incoming boundary layer. 21

23 URANS Turbulence model: Spalart-Allmaras More reasonable shock detachment but steady flow 22

24 Perspective: introduce explicit fluctuations at inlet BL fluctuations are sufficient to trigger SWBLI unsteadiness Synthetic fluctuations injection 23

25 Perspectives of UFAST project Control of experimental flows Simulations with control Large eddy simulations and hybrid RANS/LES methods Comparison experiment/simulation results and synthesis, with industry partners Database available on web: