# How to calculate Reynolds Stress using LES data?

 Register Blogs Members List Search Today's Posts Mark Forums Read February 22, 2001, 18:04 How to calculate Reynolds Stress using LES data? #1 Jongdae Kim Guest   Posts: n/a I have u,v,w velocities on staggered grid. The velocities are from Large Eddy Simulation. I'd like to compare LES results with experimental data (from LDV). Please let me know how to calculate Reynolds stress with LES data. ************************************************** ****** LES data : space averaged values Reynolds Stress (from Experimental data) : time averaged values  March 1, 2001, 08:18 Re: How to calculate Reynolds Stress using LES dat #2 Lionel Larcheveque Guest   Posts: n/a LES is by nature an unsteady computation. So to calculate components of the Reynolds stress tensor, you need to time average u, v, w, u^2, v^2, w^2, uv, uw, vw as soon as the flow seems establish and for a time period long enough. Then use the formula : =-, and so on, ..., where <> denotes time averaging, with ' defined by u'=u- (Reynolds decomposition, assuming ergodicity, i.e. time average and average on multiple identical experiences are equivalent). I suppose that space averaged values you mentioned correspond to homogeneous direction(s) of the flow. You can use averaged values along this (these) direction(s) at each time step to perform time averaging, in order to accelerate the convergence of statistical values and greatly reduce data size. saidc. likes this.  March 1, 2001, 08:24 Re: How to calculate Reynolds Stress using LES dat #3 Lionel Larcheveque Guest   Posts: n/a It seems that there is a little HTML problem in the formula I've just posted (use of <>). The correct one is : [u'v']=[uv]-[u][v], where [] denotes time averaginf, with ' defined by u'=u-[u]...  March 1, 2001, 08:43 Re: How to calculate Reynolds Stress using LES dat #4 andy Guest   Posts: n/a But what about the SGS contribution? (particularly the normal components).  March 1, 2001, 17:44 Re: How to calculate Reynolds Stress using LES dat #5 Jongdae Kim Guest   Posts: n/a Notations : { * } space (or volume) averaging [ * ] time averaging Procedure of calculation of Reynolds stresses  DNS & Experiments (1) measure u (= [u]+u') for longtime (2) u'=u-[u] (3) u'v' = u' * v' (4) [u'v'] : time averaging of (3) or [uv]-[u][v] : this is not space averaged  RANS equations with Finite volume method : tau = [u' * v'] : calculated from model = {[uv]}-{[u][v]}  LES : implicit box-filter Governing equation : space averaged-NS Numerical scheme : Finite volume method Output of numerical scheme : space averaged velocity R12={uv}-{u}{v} [R12]=[{uv}]-[{u}{v}] Questions : 1. [{*}] = {[*]} .ie. tau=[R12] ? 2. How to compare Reynolds stress of LES output with DNS(Experiment) output ? hityangsir and ahmed banna like this.  March 2, 2001, 08:45 Re: How to calculate Reynolds Stress using LES dat #6 Lionel Larcheveque Guest   Posts: n/a First my notations (slightly different from yours): u=[u]+u' Reynolds decomposition, with[*] as an average on multiple realizations which can become, using ergodicity, space, time or space and time averaging, according to the problem you want to solve (2D, 3D, statistically steady or unsteady). u={u}+u" scales separation using space filtering, with {*} a space filter. LES outputs space filtered velocity {u}. Three remarks : - I haven't understood what you was meaning by space averaging. LES is based on a space filtering : there is from a mathematical point of view a great difference between to average and to filter : the relation {{u}}={u} is false, and then {u"} differs from 0.) - I don't know much about RANS, but, for a flow with one homogeneous direction, it seems to me that : * if you consider steady 2D RANS, the average is both on time and space (you will iterate until you have a converged solution and the statistical quantities are given by the final results without "explicit" time averaging during the convergence process) :[*] means time and space averaging. * if you consider steady 3D RANS, the average is on time, and once you have converged you can average on space (and then you have the relation from your last message : tau={[uv]}-{[u]}[v]} (using your notations)) * if you consider 2D URANS, at each time step you obtain space (and maybe phase time) averaged values that you will time average during the computation. * if you consider 3D URANS, the quantities at each time steep can be regarded as phase average. There is a kind of link between this king of computations and VLES (Very large eddy simulation). - At a rough estimate, I will say that time averaging and space filtering commute on discrete level (they're both integration, but the problem is the infinite bounds of the integrals), but even then the relation you mentioned become : [R12]=[{uv}]-[{u}{v}]={[uv]}-[{u}{v}], which is different from {[uv]}-{[u][v]} (product doesn't commute). Now from the LES point of view :[*] means time averaging, and maybe space averaging if you have homogeneous flow direction(s). Using both the decompositions mentioned above, formally : u'={u}+u"-[{u}]-[u"] => [u'v']=[{u]{v}]-[{u}][{v}] : can be computed +[{u}v"]-[{u}][v"]+[u"{v}]-[u"][{v}] : cross terms +[u"v"]-[u"][v"] : subgrid terms As far as i know, as general rules the comparisons between DNS or experimental data and LES computations are made using the first line of the formula. The idea is that if you have a mesh finer enough (cut length of the implicit filter far in the inertial range of the kinetic turbulent energy spectrum), subgrid scales energy has a minor impact on the global fluctuations. If you really want to go further, here are some possibilities (not an exhaustive list): 1-Filtering your DNS data : the main problem is you don't know precisely the characteristics of the (implicit) LES filter, which is by the way far different from an academical one (gaussian, cutoff) due to the adding filtering effects of the numerical scheme. 2-Using deconvolution techniques to recalculate u". Problems similar to point 1 occur : what are the characteristics of the filter you want to inverse. To reach the same objective, maybe is it possible to use NLDE (Non Linear Disturbance Equation I think) to regenerate an estimation of the fluctuations of the small scales? 3-using models for the third (and the second?) term(s) of the formula. Models based on the Boussinesq approximation (subgrid viscosity) seems less pertinent than those based on a Bardina's approach (scales similarity) who reproduces more accurately the mathematical structure of the subgrid tensor (Is there a reason to consider turbulence only from a dissipative point of view while modeling the effects of the subgrid terms on the Reynolds stresses : the stability of the scheme is not a problem in an a posteriori computation of the Reynolds stresses). I heard about some attempts to use such approach, but it seems that few were conclusive. I hope this little text (I apologize for my English) is helping to clarify your problem. Good luck hityangsir and adhiraj like this.  February 12, 2020, 18:41 #7
New Member

ike
Join Date: Feb 2020
Posts: 4
Rep Power: 5 Hello Lionel,
I am taking analysis of turbulence and I am very new to this subject. I wanted to use matlab for analysis as well. I was wondering if you would be able to help me out with an a quesion. I have data acquired by Particle image velocimetry. I am also willing to pay for your time. please let me know

Quote:   February 13, 2020, 03:09 #8
Senior Member

Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,599
Rep Power: 70   Quote:
 Originally Posted by khattak Hello Lionel, I am taking analysis of turbulence and I am very new to this subject. I wanted to use matlab for analysis as well. I was wondering if you would be able to help me out with an a quesion. I have data acquired by Particle image velocimetry. I am also willing to pay for your time. please let me know
Have you seen the post is 20 years old?
Open your own post and write all the details of your question  Thread Tools Search this Thread Show Printable Version Email this Page Search this Thread: Advanced Search Display Modes Linear Mode Switch to Hybrid Mode Switch to Threaded Mode Posting Rules You may not post new threads You may not post replies You may not post attachments You may not edit your posts BB code is On Smilies are On [IMG] code is On HTML code is OffTrackbacks are Off Pingbacks are On Refbacks are On Forum Rules Similar Threads Thread Thread Starter Forum Replies Last Post quattro123 ANSYS 1 January 11, 2011 12:22 tstorm FLUENT 0 July 27, 2009 14:44 Roberto CFX 0 December 3, 2008 03:29 Francois FLUENT 0 July 7, 2005 11:15 Kim Bindesboll Andersen Main CFD Forum 1 September 26, 2002 13:54

All times are GMT -4. The time now is 22:24.