1 / 32

The use of dynamical RG in the development of spectral subgrid models of turbulence

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. Overview of this talk. Brief phenomenology of the statistical theory of turbulence

dionne
Télécharger la présentation

The use of dynamical RG in the development of spectral subgrid models of turbulence

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  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, University of Edinburgh

  2. 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

  3. Phenomenology

  4. relatively arbitrary Incompressible spectral Navier-Stokes eqn We will be working with the divergence-free Fourier transformed Navier-Stokes equation with no mean velocity

  5. 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)

  6. Dimensionless NSE • Move to the dimensionless form of the NSE • Where the local Reynolds number is • work in shorthand notation • where

  7. ew e Statistics (we’ll need this later) ed Characteristic dissipation length scale Richardson energy cascade

  8. From dimension arguments, Kolmogorov showed that for very large Re there exists an intermediate inertial range with scaling McComb (1990) independent of viscosity and forcing. Turbulence ‘forgets it’s roots’ Scaling, self-similarity & K41 Generalized homogeneity N=b1.585 f(x,y)=bf(baxx,bayy) Sierpinski gasket Animation from: http://classes.yale.edu/fractals/IntroToFrac/InitGen/InitGenGasket.html

  9. Here be dragons Coherent structures , etc. log E(k) Inertial Range -5/3 gradient Dissipation Range log k End of known NSE world kL kd

  10. Large Eddy Simulations (LES)- subgrid modeling problem - Approx DNS limitations go as N~Re9/4 5123 -> ~4000 Re Pipe flow transition~2x103 =k0 • 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)

  11. Dynamical RG analysis

  12. 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.

  13. D I S C L A I M E R • 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. 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)

  14. Coarse-graining

  15. Conditional average with asymptotic freedom u(k) - conditional field; w(k) - ensemble realisations

  16. ~ small ~ small

  17. Partitioned equations & the eddy viscosity Partition RG parameter space n NSE1 NSE0 Coarse-grain NSE2 NSEFP Rescale NSE3 Iterate Re Quantities being renormalized: n & local Reynolds # *M. E. Fisher, Rev. Mod. Phys. 70 2 , (1998) [Nice picture of whats happening in RG]

  18. k1=(1-h)k0 k2=(1-h)k1 RG iteration Use for LES kc k3 k2 k1 • Slightly deceptive picture/map of the RG flow -- but good to show validity of our approximations E(k) Where 0 < h < 1 k0 k

  19. RG recursive eqns and approximations ‘Assymptotic freedom’

  20. Results

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

  22. Eddy viscosity (unscaled)

  23. What eddy viscosities should look like from Direct Numerical Simulations * A. Young, PhD Thesis, University of Edinburgh (1999)

  24. E(k)=ae2/3k-5/3 * * K. Sreenivasan, Phys. Fluids 7 11, (1995); P. Yeung, Y. Zhou, PRE 56 2, (1997) Variation of the Kolmogorov constant a with shell width h

  25. TFM 2563 DNS RG 323 LES using the RG subgrid model -- comparisons K41 comparison 2563 comparison • Results from the work of C. Johnston, PhD Thesis, Edinburgh Uni (2000)

  26. Passive scalar convection * H. A. Rose, J. Fluid Mech. 81 4, (1977)

  27. RG fixed point eddy diffusivity (scaled)

  28. Prandtl number Independence Pr*=n*/c*

  29. Variation of the Kolmogorov (a) and Obukhov-Corrsin (b) constants with shell width h Ef(k)=befe-1/3k-5/3

  30. Slaved modes & Near-grid interactions u+ = u~+ + u++ Problems -- pathological divergence over here, have to introduce cut-off -> not desired

  31. ??? Questions ??? • 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)

  32. Parting thoughts • 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. Thank you! End

More Related