
[Sponsors] 
turbulentDigitalFilterInlet  how is the integral length scale tensor defined 

LinkBack  Thread Tools  Search this Thread  Display Modes 
January 31, 2020, 04:37 
turbulentDigitalFilterInlet  how is the integral length scale tensor defined

#1 
New Member
surya kaundinya
Join Date: Mar 2017
Posts: 15
Rep Power: 7 
Hello everyone,
I have a question about using the turbulentDigitalFilterInlet boundary condition for defining a turbulent velocity field in LES simulations. This is available in v1906. Do anyone know how to define the integral length scale tensor in the model. In the code I have seen that it is define in the local patch coordinate system. If that is the case, when my patch is normal to either x,y or zaxis then the tensor should be something like (L1 0 0 0 L2 0 0 0 L3). But according to the code the tensor cannot have 0 values for any of the components. Can someone please elaborate on how to define this parameter? 

January 31, 2020, 17:35 

#2 
Senior Member
Herpes Free Engineer
Join Date: Sep 2019
Location: The Home Under The Ground with the Lost Boys
Posts: 932
Rep Power: 11 
I think L definition is as follows:
(Luux Luuy Luuz Lvvx Lvvy Lvvz Lwwx Lwwy Lwwz) That is, for example, for the first Luux: uu correlation in the xdirection; for the second item Luuy: uu correlation in the ydirection and so on. To my knowledge, for turbulent fields, this set of correlations cannot be zero. The code then seems to be guarded from such input to reduce potential user errors. Zero correlation means complete randomness, and turbulence is not random, yet stochastic. Hope this helps.
__________________
The OpenFOAM community is the biggest contributor to OpenFOAM: User guide/Wiki1/Wiki2/Code guide/Code Wiki/Journal Nilsson/Guerrero/Holzinger/Holzmann/Nagy/Santos/Nozaki/Jasak/Primer Governance Bugs/Features: OpenFOAM (ESIOpenCFDTrademark) Bugs/Features: FOAMExtend (WikkiFSB) Bugs: OpenFOAM.org How to create a MWE New: Forkable OpenFOAM mirror 

February 1, 2020, 05:48 

#3 
New Member
surya kaundinya
Join Date: Mar 2017
Posts: 15
Rep Power: 7 
Hi,
I think it does not represent the velocity correlation. The code description says that it represents the "integral length scale" as sshown below: L is the Integral lengthscale set (9comp):{e1,e2,e3}{u,v,w} [m] L = (<Lxu> <Lxv> <Lxw> <Lyu> <Lyv> <Lyw> <Lzu> <Lzv> <Lzw>) It is also mentioned that the "integral length scale set input is in the local patch coordinate system". Also it is specified that "L_e1u, Le1v, Le2w, should always correspond to the length scales that are in association with the convective mean flow direction." From these I think that Lx, Ly, Lz are integral length scales in different directions. And e1, e2, e3 are the local patch coordinates. So my question is how to define the these components where patch is parallel to a specific axis. With regards, Surya. 

February 1, 2020, 07:10 

#4 
Senior Member
Herpes Free Engineer
Join Date: Sep 2019
Location: The Home Under The Ground with the Lost Boys
Posts: 932
Rep Power: 11 
Integral length scale there is only the pesudointegration of the twopoint correlation function of velocity in a given direction. They contain the same information.
__________________
The OpenFOAM community is the biggest contributor to OpenFOAM: User guide/Wiki1/Wiki2/Code guide/Code Wiki/Journal Nilsson/Guerrero/Holzinger/Holzmann/Nagy/Santos/Nozaki/Jasak/Primer Governance Bugs/Features: OpenFOAM (ESIOpenCFDTrademark) Bugs/Features: FOAMExtend (WikkiFSB) Bugs: OpenFOAM.org How to create a MWE New: Forkable OpenFOAM mirror 

February 1, 2020, 07:19 

#5 
Senior Member
Herpes Free Engineer
Join Date: Sep 2019
Location: The Home Under The Ground with the Lost Boys
Posts: 932
Rep Power: 11 
"From these I think that Lx, Ly, Lz are integral length scales in different directions. And e1, e2, e3 are the local patch coordinates. So my question is how to define the these components where patch is parallel to a specific axis."
e1 is your mean convection direction, say x. e2 and e3 are the orthogonal directions to e1, say y and z. Now say, you clockwise rotate your mesh 60 degrees around the global xaxis by keeping the L input the same. In the local coordinate system, which does not feel the rotation, if you measure Ls, you will get the input back. But if you do the integral scale measurement in the rotated mesh on global coordinate system, you will get different set of Ls in the global directions. That's it.
__________________
The OpenFOAM community is the biggest contributor to OpenFOAM: User guide/Wiki1/Wiki2/Code guide/Code Wiki/Journal Nilsson/Guerrero/Holzinger/Holzmann/Nagy/Santos/Nozaki/Jasak/Primer Governance Bugs/Features: OpenFOAM (ESIOpenCFDTrademark) Bugs/Features: FOAMExtend (WikkiFSB) Bugs: OpenFOAM.org How to create a MWE New: Forkable OpenFOAM mirror 

February 2, 2020, 11:14 

#6 
New Member
surya kaundinya
Join Date: Mar 2017
Posts: 15
Rep Power: 7 
Hi HPE,
Thanks a lot for the information. I might sound stupid but can you please answer one last question? For a case where I do not have DNS data , how do we calculate the twopoint correlation? Or if I know the integral length scale how to calculate the twopoint correlation? With regards, Surya. 

September 24, 2022, 10:15 

#7 
New Member
Justin
Join Date: Mar 2022
Posts: 6
Rep Power: 2 
If you do not have DNS data, then you cannot compute the twopoint correlation. The integral length scale is defined as a length scale that corresponds to the correlation time of a turbulent cascade. In principle, if you know the spectral scaling, you can get the structure function and estimate how the correlation function scales, but you cannot get exact values. You may try simulating without turbulent inflow BC and see how that goes.


Thread Tools  Search this Thread 
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Calculation of turbulence length scale  Jinjolee  OpenFOAM Running, Solving & CFD  15  April 21, 2022 11:44 
Problem with large time scale and minute length scale  mech_engg_swag  OpenFOAM  0  August 23, 2016 08:39 
Estimation of turbulence length scale  J.Y.Won  OpenFOAM PreProcessing  0  April 28, 2015 08:16 
length scale option  rashmi  CFX  2  July 25, 2005 00:39 
Characteristic Length Scale  Thomas P. Abraham  Main CFD Forum  2  April 20, 1999 18:16 