the use of dynamical rg in the development of spectral subgrid models of turbulence khurom kiyani,...

32
The use of dynamical RG in the development of spectral subgrid models of turbulence Khurom Kiyani, David McComb Turbulence Theory group, School of Physics, University of Edinburgh

Upload: rachael-tugwell

Post on 31-Mar-2015

214 views

Category:

Documents


0 download

TRANSCRIPT

Page 1: The use of dynamical RG in the development of spectral subgrid models of turbulence Khurom Kiyani, David McComb Turbulence Theory group, School of Physics,

The use of dynamical RG in the development of spectral subgrid

models of turbulence

Khurom Kiyani, David McCombTurbulence Theory group, School of Physics,

University of Edinburgh

Page 2: The use of dynamical RG in the development of spectral subgrid models of turbulence Khurom Kiyani, David McComb Turbulence Theory group, School of Physics,

Overview of this talk

• Brief phenomenology of the statistical theory of turbulence

• Large-eddy simulations (LES) & subgrid modeling

• Dynamical renormalization group (RG) method

• Results

- Homogeneous & isotropic turbulence

- Passive scalar advection (by above)

- Other LES comparisons

• Problems with the current schemes - introduction of slaved modes to handle near-grid terms

Page 3: The use of dynamical RG in the development of spectral subgrid models of turbulence Khurom Kiyani, David McComb Turbulence Theory group, School of Physics,

Phenomenology

Page 4: The use of dynamical RG in the development of spectral subgrid models of turbulence Khurom Kiyani, David McComb Turbulence Theory group, School of Physics,

Incompressible spectral Navier-Stokes eqn

We will be working with the divergence-free Fourier transformed Navier-Stokes equation with no mean velocity

relatively arbitrary

Page 5: The use of dynamical RG in the development of spectral subgrid models of turbulence Khurom Kiyani, David McComb Turbulence Theory group, School of Physics,

Homogeneous, Isotropic & stationary NSE for infinite fluid

The simplest non-trivial case -- shrink the monster to a smaller monster. Makes the maths a bit easier.

• k-space allows us to deal directly with the many strongly coupled degrees of freedom.

• Statistically steady state - the only reason why we have included f in NSE.

• No mean velocity implicitly implies global isotropy

Leaves us with the most quintessential, unadulterated turbulence -- but pretty artificial(ish)

Page 6: The use of dynamical RG in the development of spectral subgrid models of turbulence Khurom Kiyani, David McComb Turbulence Theory group, School of Physics,

Move to the dimensionless form of the NSE

Where the local Reynolds number is

work in shorthand notation

where

Dimensionless NSE

Page 7: The use of dynamical RG in the development of spectral subgrid models of turbulence Khurom Kiyani, David McComb Turbulence Theory group, School of Physics,

Richardson energy cascade

Statistics (we’ll need this later)

w

d

Characteristic dissipation length scale

Page 8: The use of dynamical RG in the development of spectral subgrid models of turbulence Khurom Kiyani, David McComb Turbulence Theory group, School of Physics,

Scaling, self-similarity & K41

Sierpinski gasket

Animation from: http://classes.yale.edu/fractals/IntroToFrac/InitGen/InitGenGasket.html

N=b1.585

f(x,y)=bf(baxx,bayy)

Generalized homogeneityQuickTime™ and aGIF decompressor

are needed to see this picture.

From dimension arguments, Kolmogorov showed that for very large Re there exists an intermediate inertial range with scaling

independent of viscosity and forcing.

Turbulence ‘forgets it’s roots’

McComb (1990)

Page 9: The use of dynamical RG in the development of spectral subgrid models of turbulence Khurom Kiyani, David McComb Turbulence Theory group, School of Physics,

log E(k)

log k

Here be dragons

Coherent structures , etc.

Inertial Range -5/3 gradient

Dissipation Range

kdkL

End of known NSE world

Page 10: The use of dynamical RG in the development of spectral subgrid models of turbulence Khurom Kiyani, David McComb Turbulence Theory group, School of Physics,

Large Eddy Simulations (LES)- subgrid modeling problem -

• Aim: To model the large scales of a turbulent flow whilst accounting for the missing scales in an appropriate way.

• Using a sharp spectral filter (Heaviside unit-step fn)

=k0

Approx DNS limitations go as

N~Re9/4

5123 -> ~4000 Re

Pipe flow transition~2x103

Page 11: The use of dynamical RG in the development of spectral subgrid models of turbulence Khurom Kiyani, David McComb Turbulence Theory group, School of Physics,

Dynamical RG analysis

Page 12: The use of dynamical RG in the development of spectral subgrid models of turbulence Khurom Kiyani, David McComb Turbulence Theory group, School of Physics,

Renormalization Group (RG)We can find what kc is and the form of the eddy viscosity using Renormalization Group (RG) techniques.

What is RG?

• RG is an iterative method for reducing the number of degrees of freedom (DoF’s) in a problem involving many DoF’s.

•In our context of fluid turbulence, this can be interpreted as the elimination of Fourier velocity fluctuation modes.

RG in k-space

• Coarse-grain or average out the effect of the high-k modes and add it onto the kinematic viscosity.

•Rescale the variables so that the new renormalized NSE look like the original one.

•Repeat until you get to a fixed point - picture does not change.

Page 13: The use of dynamical RG in the development of spectral subgrid models of turbulence Khurom Kiyani, David McComb Turbulence Theory group, School of Physics,

• Non-equillibrium phenomena different (nastier, richer) monster from equillibrium physics -- analogies to ferromagnetism etc. quite hard; Don’t quite know what the order parameter is* (ask me about this at the end).

• Confining ourselves to LES - so no critical exponents etc. calculated -- don’t think anyone has obtained K41 from NSE using RG.

• RG has to be formulated appropriately/delicately -- not a magic black box -> exponents, renormalized quantities etc. You really have to have an inclination of the ‘physics’ before you start RG’ing.• Involves approximations (often) and blatant abuses.

However…Very deep and profound ideas of the perceived physics of the system and explanation of universality in physically distinct systems

D I S C L A I M E R

D. Forster et al., Large-distance and long-time properties of a randomly stirred fluid, PRA 16 2, (1977)

* M. Nelkin, PRA 9,1 (1974); Zhou, McComb, Vahala -- icase 36 (1997)

Page 14: The use of dynamical RG in the development of spectral subgrid models of turbulence Khurom Kiyani, David McComb Turbulence Theory group, School of Physics,

Coarse-graining

Page 15: The use of dynamical RG in the development of spectral subgrid models of turbulence Khurom Kiyani, David McComb Turbulence Theory group, School of Physics,

Conditional average with asymptotic freedom

u(k) - conditional field; w(k) - ensemble realisations

Page 16: The use of dynamical RG in the development of spectral subgrid models of turbulence Khurom Kiyani, David McComb Turbulence Theory group, School of Physics,

~ small

~ small

Page 17: The use of dynamical RG in the development of spectral subgrid models of turbulence Khurom Kiyani, David McComb Turbulence Theory group, School of Physics,

Partitioned equations & the eddy viscosity

Partition

Coarse-grain

Rescale

Quantities being renormalized: & local Reynolds #

Iterate

NSE0NSE1

NSE2

NSE3

NSEFP

Re

RG parameter

space

*M. E. Fisher, Rev. Mod. Phys. 70 2 , (1998) [Nice picture of whats happening in RG]

Page 18: The use of dynamical RG in the development of spectral subgrid models of turbulence Khurom Kiyani, David McComb Turbulence Theory group, School of Physics,

k

E(k)

k0k1

k1=(1-)k0

Where

0 < < 1

k3kc

RG iteration

Use for LES

k2

k2=(1-)k1

• Slightly deceptive picture/map of the RG flow -- but good to show validity of our approximations

Page 19: The use of dynamical RG in the development of spectral subgrid models of turbulence Khurom Kiyani, David McComb Turbulence Theory group, School of Physics,

RG recursive eqns and approximations

‘Assymptotic freedom’

Page 20: The use of dynamical RG in the development of spectral subgrid models of turbulence Khurom Kiyani, David McComb Turbulence Theory group, School of Physics,

Results

Page 21: The use of dynamical RG in the development of spectral subgrid models of turbulence Khurom Kiyani, David McComb Turbulence Theory group, School of Physics,

RG map - Evolution of (scaled) eddy viscosity with RG iteration

Page 22: The use of dynamical RG in the development of spectral subgrid models of turbulence Khurom Kiyani, David McComb Turbulence Theory group, School of Physics,

Eddy viscosity (unscaled)

Page 23: The use of dynamical RG in the development of spectral subgrid models of turbulence Khurom Kiyani, David McComb Turbulence Theory group, School of Physics,

What eddy viscosities should look like from Direct Numerical Simulations

* A. Young, PhD Thesis, University of Edinburgh (1999)

Page 24: The use of dynamical RG in the development of spectral subgrid models of turbulence Khurom Kiyani, David McComb Turbulence Theory group, School of Physics,

Variation of the Kolmogorov constant with shell width

E(k)=k-5/3

*

* K. Sreenivasan, Phys. Fluids 7 11, (1995); P. Yeung, Y. Zhou, PRE 56 2, (1997)

Page 25: The use of dynamical RG in the development of spectral subgrid models of turbulence Khurom Kiyani, David McComb Turbulence Theory group, School of Physics,

323 LES using the RG subgrid model -- comparisons

RG

TFM

2563 DNS

Model Chi2

RG 48.6

TFM 77.3

DNS 88.21

Model Chi2

RG 205.9

TFM 276.3

DNS 35.9

K41 comparison

2563 comparison

Results from the work of C. Johnston, PhD Thesis, Edinburgh Uni (2000)

Page 26: The use of dynamical RG in the development of spectral subgrid models of turbulence Khurom Kiyani, David McComb Turbulence Theory group, School of Physics,

Passive scalar convection

* H. A. Rose, J. Fluid Mech. 81 4, (1977)

Page 27: The use of dynamical RG in the development of spectral subgrid models of turbulence Khurom Kiyani, David McComb Turbulence Theory group, School of Physics,

RG fixed point eddy diffusivity (scaled)

Page 28: The use of dynamical RG in the development of spectral subgrid models of turbulence Khurom Kiyani, David McComb Turbulence Theory group, School of Physics,

Prandtl number Independence

Pr*=

Page 29: The use of dynamical RG in the development of spectral subgrid models of turbulence Khurom Kiyani, David McComb Turbulence Theory group, School of Physics,

Variation of the Kolmogorov ( and Obukhov-Corrsin () constants with shell width

E(k)=k-5/3

Page 30: The use of dynamical RG in the development of spectral subgrid models of turbulence Khurom Kiyani, David McComb Turbulence Theory group, School of Physics,

Slaved modes & Near-grid interactions

u+ = u~+ + u++Problems -- pathological divergence over here, have to introduce cut-off -> not desired

Page 31: The use of dynamical RG in the development of spectral subgrid models of turbulence Khurom Kiyani, David McComb Turbulence Theory group, School of Physics,

• The reason why we do not introduce extra couplings is due to us not wanting to compute higher order terms like u-u-u- in an LES -- it would be a poor subgrid model.

• Pessimistic - Possible existence of infinite number of marginal scaling fields (your approximations are never good enough)*. • Optimistic - Apart from cusp behaviour results are not doing too bad. Get pretty good values for ‘universal’ constants. Eddy viscosity performs just as well as other leading brands**.

* G. Eyink, Phys. Fluids 6 9, (1994)** McComb et. al -- (see next slide)

??? Questions ???

Page 32: The use of dynamical RG in the development of spectral subgrid models of turbulence Khurom Kiyani, David McComb Turbulence Theory group, School of Physics,

• RG of McComb et al. has been used in actual LES.

- W. D. McComb et al., Phys. Fluids 13 7 (2001)

- C. Johnston PhD Thesis, Edinburgh Uni (2000)

• Need more analysis on including near-grid cross terms. Look at some way of ascertaining fixed point behaviour of different terms/couplings (relevant scaling fields etc.)

• Maybe have a look at non-perturbative variational approaches.

Parting thoughts

Thank you!

End