
[Sponsors] 
November 13, 2012, 08:13 
LES SpalartAllmaras Simulation questions

#1 
New Member
Matthias K.
Join Date: Oct 2012
Location: Germany Erlangen
Posts: 20
Rep Power: 13 
Hi,
I am doing a LES simulation of a ball, which is flowed by air. I am calculating the draft coefficients of the ball with different Reynolds' numbers. My SpalartAllmaras LES model is running, but I have a few questions regarding to the initial values. In my 0directory there is a nuSgs and a B file. What are the functions of the files. How can I calculate these values? Are there any formulas available? Should the nutilda be calculated like the RAS SpalartAllmaras case with the well known formulas? The nut is not used in the LES SpalartAllmaras case, am I right? Regards, Matthias 

November 13, 2012, 09:18 

#2 
New Member
Matthias K.
Join Date: Oct 2012
Location: Germany Erlangen
Posts: 20
Rep Power: 13 
I have only found the following formulas for B and nuSgs in the forum.
B = 2/3 rho*k*I2*mu_t*D nuSgs_ = ck_*delta()*sqrt(k(gradU)); But I am not sure if these are applicable and what the abbreviations mean. Regards Matthias 

November 23, 2012, 09:28 

#3 
New Member
Matthias K.
Join Date: Oct 2012
Location: Germany Erlangen
Posts: 20
Rep Power: 13 
Hi,
My Simulation runs without B. Now then I only need the definitions for the nuSGS and the nuTilda. If the nuTilda is calculated as in the RAS Spalard Allmaras case, I only need the definition of nuSGS. On http://foam.sourceforge.net/docs/cpp/a01896.html I found the following formula. nuSgs = ck*sqrt(k)*delta But ck is only described as a dimensioned skalar and delta isn't described anyway. My Problem is, that the coefficients of drag of the ball are negative, so I think I have to set a value for nuSgs. At the moment nuSGS is set to 1e10. Hopefully somebody is able to help me Regards Matthias 

November 23, 2012, 09:55 

#4 
Member
Meindert de Groot
Join Date: Jun 2012
Location: Netherlands
Posts: 34
Rep Power: 13 
Hi Matthias,
Have a look at http://foam.sourceforge.net/docs/cpp/a08492_source.html. nuTilda is given by the solution to the equation on lines 297311. The function on lines 4549 is called to get nuSgs. A discussion about the values of nut (or nuSgs in this case) and nuTilda came up not too long ago. You might want to have a look at that discussion: http://www.cfdonline.com/Forums/ope...tnutilda.html. 

November 27, 2012, 05:51 

#5 
New Member
Matthias K.
Join Date: Oct 2012
Location: Germany Erlangen
Posts: 20
Rep Power: 13 
Hi,
Thank you for your help. I have read the post of your link. The following formula is given for fv1 x^3/(x^3 + cv1^3) with x=nut/nu cv1=7.1 According to nuSgs=nuTilda*fv1 I should get my initial value for nuSgs. Am I right? Regards Matthias 

November 27, 2012, 06:55 

#6 
Member
Meindert de Groot
Join Date: Jun 2012
Location: Netherlands
Posts: 34
Rep Power: 13 
Be aware that nut and nuTilda are not the same. Actually, nut is called nuSgs in this case. X is defined as nuTilda/nu (see lines 5458 in foam.sourceforge.net/docs/cpp/a08492_source.html).
I do not fully understand your question about nuSgs. If you are asking whether you should initialise nuSgs, then the answer is: Yes, you should. Of course, you could test this very easily by removing nuSgs from your time directory. I expect your simulation will not run. 

November 27, 2012, 16:04 

#7 
New Member
Matthias K.
Join Date: Oct 2012
Location: Germany Erlangen
Posts: 20
Rep Power: 13 
Hi,
Thank you very much for your help. I hadn't known, that nut is called nuSgs in this case. Regards Matthias 

December 26, 2012, 17:03 

#8 
New Member
Matthias K.
Join Date: Oct 2012
Location: Germany Erlangen
Posts: 20
Rep Power: 13 
Hi,
I had a few problems with the 3Dcase and the calculation time was very high. So I set up a 2D case (extrudeMesh) of a long cylinder to get rid of the problems. I chose a low Reynold number of 180. According to the book Large Eddy Simulation from J. Fröhlich, it is suitable to simulate 2D in case of low Reynold numbers. My problem is, that I should determine the drag coefficent of the cylinder. To get my drag coefficients, I inserted forceCoeffs in my control.dict file. forceCoeffs { type forceCoeffs; functionObjectLibs ("libforces.so"); patches (zylinder_region0); //change to your patch name pName p; UName U; rhoName rhoInf; rhoInf 1.225; CofR (0 0 0); liftDir (0 0 1); dragDir (1 0 0); pitchAxis (0 1 0); magUInf 0.0243243243; lRef 0.0555; Aref 1.11e3; outputControl timeStep; outputInterval 10; } To get the dragcoefficient cd of this instable case I made the average of the last 500 values. My problem is, that the figure of cd is in the 100Mio. range, where it should not be. Is this the right way, to determine drag Coefficients of an LES case? My yPlus Values of the cylinder are the following: y+ : min: 0.652849 max: 59.6055 average: 22.989 Thank you for your time Regards, Matthias 

December 26, 2012, 17:08 

#9 
New Member
Matthias K.
Join Date: Oct 2012
Location: Germany Erlangen
Posts: 20
Rep Power: 13 
Here are 2 pics of the case
Regards Matthias 

January 3, 2013, 15:17 

#10 
New Member
Matthias K.
Join Date: Oct 2012
Location: Germany Erlangen
Posts: 20
Rep Power: 13 
Hi,
The problem with the high cd values is solved. There were problems with the boundary conditions. Thanks, Matthias 

May 24, 2015, 07:26 
«les

#11 
New Member
javid hoseinpour
Join Date: Nov 2014
Posts: 1
Rep Power: 0 
Hi
How do I set coefficients k and nusgs for LESSpalart 

June 15, 2015, 06:08 

#12 
New Member
Somsankar Mondal
Join Date: Dec 2014
Posts: 4
Rep Power: 11 

Thread Tools  Search this Thread 
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
cyclone LES simulation prob: ERROR # 004100018  sakurabogoda  CFX  15  November 6, 2012 09:07 
Could you do me a favor: regarding questions in PEMFC simulation  bamiao0823  FLUENT  0  March 20, 2012 02:59 
Rotor stator interaction LES simulation  Thomas_Lloyd  OpenFOAM  0  May 4, 2011 07:05 
How to determinate turbulence scale in LES (laminar simulation)?  Franciswu21  ANSYS  0  October 22, 2009 12:48 
LES simulation of bluff body at Re(d)=1200  arash_nl  FLUENT  0  June 1, 2009 13:16 