Turbulent eddies in the RANS/LES transition region

Transcription

1 Turbulent eddies in the RANS/LES transition region Ugo Piomelli Senthil Radhakrishnan Giuseppe De Prisco University of Maryland College Park, MD, USA Research sponsored by the ONR and AFOSR

2 Outline Motivation The problem: eddy generation at the RANS/LES interface Effects and possible solutions WMLES Zonal RANS Conclusions and directions for improvement

3 Motivation Computational approaches for the simulation of an aircraft (from Spalart, 2000) Accurate methods are infeasible. Feasible methods are (often) inaccurate. Hybrid RANS/LES: Use (U)RANS in regions in which models are accurate. Use LES in non-equilibrium regions (separation, 3D mean flow, high pressure gradients) or where structural information is required (noise emission).

4 DES Attached boundary layer URANS, everything else LES. Detached-eddy simulation (DES)

5 WMLES Contours of u'v' ν T du / dy LES URANS Wall layer URANS, everything else LES. Wall-Modeled LES (WMLES) Oldest hybrid application (logarithmic law)

6 Zonal RANS/LES Attached boundary layer URANS, LES includes attached & separated flows.

7 RANS/LES interface Critical issue: RANS/LES interface. RANS: Reynolds stress supported by the model. ν T du dy? u'v'. Flow in a compressor and prediffuser. From Schlüter et al., AIAA Paper

8 RANS/LES interface Critical issue: RANS/LES interface. RANS: Reynolds stress supported by the model LES: Reynolds stress supported by the eddies. ν T du dy? u'v'. ν T du dy = u 'v'. Flow in a compressor and prediffuser. From Schlüter et al., AIAA Paper

9 RANS/LES interface Critical issue: RANS/LES interface. RANS: Reynolds stress supported by the model ν T du dy? u'v'. LES: Reynolds stress supported by the eddies ν T du dy = u 'v'. Turbulent eddies must be generated at the interface. How? Flow in a compressor and prediffuser. From Schlüter et al., AIAA Paper

10 RANS/LES interface Critical issue: RANS/LES interface. Rapid generation of eddies as the model switches from RANS to LES behavior can be achieved by: Natural amplification of instabilities. o Shear layers: OK. Flow in a compressor and prediffuser. From Schlüter et al., AIAA Paper

11 RANS/LES interface Critical issue: RANS/LES interface. Rapid generation of eddies as the model switches from RANS to LES behavior can be achieved by: Natural amplification of instabilities. o Shear layers: OK. o Attached b.l.: less effective. IDDES.

12 RANS/LES interface Critical issue: RANS/LES interface. Rapid generation of eddies as the model switches from RANS to LES behavior can be achieved by: Natural amplification of instabilities. Artificial forcing. o Synthetic turbulence. o Disturbances from similar calculation. o Controlled forcing. RANS below LES RANS into LES

13 Outline Motivation The problem: eddy generation at the RANS/LES interface Effects and possible solutions WMLES Zonal RANS Conclusions and directions for improvement

14 WMLES using hybrid RANS/LES Two main methodologies: Blending function: Compute RANS and SGS eddy viscosity using different models. Blend them using a specified ad hoc function. (Tokyo), Leschziner (Imperial College), Davidson (Chalmers), Edwards (NCSU)... Detached eddy simulation: Use a single model in the RANS and LES regions. Modify the model (length scale) to account for different physics. Nikitin et al. (2000), Piomelli et al. (2003), Pasinato et al. (2005), Keating and Piomelli (2006), Radhakrishnan et al. (2006). Main effect of the absence of turbulent eddies at the RANS/LES interface: logarithmic law mismatch (LLM).

15 WMLES using hybrid RANS/LES Logarithmic law mismatch RANS log layer LES log layer Plane channel flow, Re τ =5,000

16 WMLES using hybrid RANS/LES Logarithmic law mismatch Resolved stress Modeled stress Plane channel flow, Re τ =5,000

17 WMLES using hybrid RANS/LES Logarithmic law mismatch Nominal LES region y > C DES Δ Resolved stress Actual LES region Resolved > Modeled Modeled stress Transition region (DES buffer layer) Plane channel flow, Re τ =5,000

18 WMLES of the flow over a ramp Experiment: Song & Eaton (2003) Calculations Re θ = 21,000 at reference location x = 2 Co-located curvilinear FD code (2 nd order in space and time) LES with DES-based wall-layer model ( ), RANS. Challenging physics: Shallow, pressure-driven separation. Prediction of the flow after separation depends critically on the accuracy of the mean-velocity prediction.

19 WMLES of the flow over a ramp RANS WMLES Experiment

20 WMLES of the flow over a ramp Isosurfaces of Q = 1 ( 2 S2 Ω 2 ) Contours of u in a near-wall plane

21 WMLES of the flow over a ramp Experiment WMLES

22 Resolved-eddy enhancement A transition problem? Smooth, laminar-like flow in the inner layer. Turbulent flow in the outer layer. How to accelerate the transition to turbulence in the LES region? Diffusion dominated advection dominated regime

23 Resolved-eddy enhancement A transition problem? Smooth, laminar-like flow in the inner layer. Turbulent flow in the outer layer How to accelerate the transition to turbulence in the LES region? Diffusion dominated advection dominated regime Possible solution: add perturbations to stir the flow. Piomelli et al. (2003) Random forcing to generate small-scale fluctuations in the RANS/LES transition region. The random fluctuations are massaged by the strain field and become eddies. Forcing amplitude set to match resolved and modelled Reynolds stresses over the transition region:

24 WMLES of the flow over a ramp Isosurfaces of Q = 1 ( 2 S2 Ω 2 ) Contours of u in a near-wall plane

25 WMLES of the flow over a ramp

26 WMLES of the flow over a ramp RANS WMLES Experiment WMLES, stochastic force

27 WMLES of the flow over a ramp Experiment WMLES stochastic force WMLES no force RANS

28 WMLES of the flow over a ramp Experiment WMLES, no force WMLES, stochastic force

29 Outline Motivation The problem: eddy generation at the RANS/LES interface Effects and possible solutions WMLES Zonal RANS Conclusions and directions for improvement

30 Zonal Hybrid RANS/LES strategies Two approaches: Integrated simulation (DES, Menon, ) Single grid, model changes. Separate simulation (CTR, Sagaut, ) RANS data used to assign boundary conditions for LES. Equivalent to inflow assignment for DNS/LES. Generation of eddies by: Growth of natural disturbances Synthetic turbulence Synthetic turbulence + controlled forcing

31 Information transfer between RANS & LES RANS gives: Mean flow Reynolds stresses Always u v Sometimes TKE Sometimes u u, v v and w w LES requires: Instantaneous u, v and w. Spectra and phase relations. Synthetic turbulence can be constructed to give Assigned mean flow and Reynolds stresses Assigned spectra No phase relations

32 Channel flow. Synthetic turbulence at the RANS/LES interface Controlled

33 Channel flow. Synthetic turbulence at the RANS/LES interface The flow rapidly loses turbulent kinetic energy and begins to relaminarize. Reference Eventually, the flow transitions and reaches acceptable turbulence levels 20δ downstream of the inflow. Synthetic Shear stress Mean velocity x/δ = 10 x/δ = 15 x/δ = 20

34 Controlled forcing at the RANS/LES interface Philosophy: Generate reasonably realistic turbulence through inflow conditions or forcing. Spectra Stresses Selectively amplify bursts to establish the correct shear stress profile. Ingredients: Synthetic turbulence Controlled forcing

35 Synthetic turbulence Batten, Goldberg and Chakravarthy AIAA J. 42, 485 (2004) Three-dimensional, unsteady velocity field Mean flow from RANS data Fluctuations with TKE and u v from RANS data. Length and time scales from the RANS data. E(k) ~ k 2 exp(- k 4 ) Possiblyanisotropic

36 Controlled forcing Spille-Kohoff and Kaltenbach. In DNS/LES Progress and Challenges (Liu, Sakell & Beutner eds.) 319 (2001) Add forcing term to the v momentum equation at a number of control planes downstream of the interface. Use a controller to drive the Reynolds shear stress towards a target Reynolds shear stress.

37 Channel flow. Controlled forcing at the RANS/LES interface The flow adjusts within 10-15δ Reference Controlled forcing Synthetic Shear stress Mean velocity x/δ = 10 x/δ = 15 x/δ = 20

38 Channel flow. Controlled forcing at the RANS/LES interface Synthetic Controlled

39 Decelerating boundary layer Calculations of the flow on a flat plate with variable freestream velocity. Cartesian staggered code, 2 nd order in space and time. Freestream velocity points (reference calculation) points (hybrid calculation) at the inlet

40 Decelerating boundary layer Freestream velocity Skin-friction coefficient SA-RANS Controlled Synthetic

41 Decelerating boundary layer SA-RANS Synthetic Controlled Synthetic Controlled SA-RANS

42 Decelerating boundary layer Reference Synthetic turbulence + controlled forcing

43 Conclusions The interface between RANS and LES zones may affect critically the accuracy of the flow predictions. Separation. Turbulent kinetic energy levels The need for turbulent eddies in the LES region is recognized. Several solutions have been proposed. Synthetic turbulence Forcing (DNS databases, controlled,.) Decreased eddy viscosity Partial success so far. Phase information is crucial. Some flows are more forgiving.

44 Directions for future work Improved integration between turbulent physics and model. Better understanding of the stability characteristics of the system: Smooth, laminar-like flow in the inner layer. Diffusion dominated. Turbulent flow in the outer layer. Advection dominated. Identification of optimal disturbances.