CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

LES SpalartAllmaras Simulation questions

Register Blogs Members List Search Today's Posts Mark Forums Read

Like Tree3Likes
  • 1 Post By meindert
  • 1 Post By meindert
  • 1 Post By Matt1986

Reply
 
LinkBack Thread Tools Display Modes
Old   November 13, 2012, 09:13
Default LES SpalartAllmaras Simulation questions
  #1
New Member
 
Matt1986's Avatar
 
Matthias K.
Join Date: Oct 2012
Location: Germany Erlangen
Posts: 20
Rep Power: 4
Matt1986 is on a distinguished road
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 0-directory 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
Matt1986 is offline   Reply With Quote

Old   November 13, 2012, 10:18
Default
  #2
New Member
 
Matt1986's Avatar
 
Matthias K.
Join Date: Oct 2012
Location: Germany Erlangen
Posts: 20
Rep Power: 4
Matt1986 is on a distinguished road
I have only found the following formulas for B and nuSgs in the forum.

B = 2/3 rho*k*I-2*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
Matt1986 is offline   Reply With Quote

Old   November 23, 2012, 10:28
Default
  #3
New Member
 
Matt1986's Avatar
 
Matthias K.
Join Date: Oct 2012
Location: Germany Erlangen
Posts: 20
Rep Power: 4
Matt1986 is on a distinguished road
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 1e-10.

Hopefully somebody is able to help me

Regards Matthias
Attached Images
File Type: jpg spalardAllmaras.jpg (40.2 KB, 26 views)
Matt1986 is offline   Reply With Quote

Old   November 23, 2012, 10:55
Default
  #4
Member
 
Meindert de Groot
Join Date: Jun 2012
Location: Netherlands
Posts: 34
Rep Power: 5
meindert is on a distinguished road
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 297-311. The function on lines 45-49 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: nut and nuTilda.
Matt1986 likes this.
meindert is offline   Reply With Quote

Old   November 27, 2012, 06:51
Default
  #5
New Member
 
Matt1986's Avatar
 
Matthias K.
Join Date: Oct 2012
Location: Germany Erlangen
Posts: 20
Rep Power: 4
Matt1986 is on a distinguished road
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
Matt1986 is offline   Reply With Quote

Old   November 27, 2012, 07:55
Default
  #6
Member
 
Meindert de Groot
Join Date: Jun 2012
Location: Netherlands
Posts: 34
Rep Power: 5
meindert is on a distinguished road
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 54-58 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.
Matt1986 likes this.
meindert is offline   Reply With Quote

Old   November 27, 2012, 17:04
Default
  #7
New Member
 
Matt1986's Avatar
 
Matthias K.
Join Date: Oct 2012
Location: Germany Erlangen
Posts: 20
Rep Power: 4
Matt1986 is on a distinguished road
Hi,

Thank you very much for your help.
I hadn't known, that nut is called nuSgs in this case.


Regards Matthias
Matt1986 is offline   Reply With Quote

Old   December 26, 2012, 18:03
Default
  #8
New Member
 
Matt1986's Avatar
 
Matthias K.
Join Date: Oct 2012
Location: Germany Erlangen
Posts: 20
Rep Power: 4
Matt1986 is on a distinguished road
Hi,

I had a few problems with the 3D-case 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.11e-3;
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
Matt1986 is offline   Reply With Quote

Old   December 26, 2012, 18:08
Default
  #9
New Member
 
Matt1986's Avatar
 
Matthias K.
Join Date: Oct 2012
Location: Germany Erlangen
Posts: 20
Rep Power: 4
Matt1986 is on a distinguished road
Here are 2 pics of the case

Regards Matthias
Attached Images
File Type: jpg 90518bec.l.jpg (14.3 KB, 43 views)
File Type: jpg e794953b.l.jpg (76.1 KB, 45 views)
Matt1986 is offline   Reply With Quote

Old   January 3, 2013, 16:17
Default
  #10
New Member
 
Matt1986's Avatar
 
Matthias K.
Join Date: Oct 2012
Location: Germany Erlangen
Posts: 20
Rep Power: 4
Matt1986 is on a distinguished road
Hi,

The problem with the high cd values is solved.
There were problems with the boundary conditions.

Thanks,
Matthias
somFoam likes this.
Matt1986 is offline   Reply With Quote

Old   May 24, 2015, 07:26
Default «les
  #11
New Member
 
javid hoseinpour
Join Date: Nov 2014
Posts: 1
Rep Power: 0
javidFoam is on a distinguished road
Hi
How do I set coefficients k and nusgs for LES-Spalart
javidFoam is offline   Reply With Quote

Old   June 15, 2015, 06:08
Default
  #12
New Member
 
Somsankar Mondal
Join Date: Dec 2014
Posts: 6
Rep Power: 2
somFoam is on a distinguished road
Quote:
Originally Posted by Matt1986 View Post
Hi,

The problem with the high cd values is solved.
There were problems with the boundary conditions.

Thanks,
Matthias
Can you please give some reference how to calculate the nuSgs for simple smagorinsky model.
somFoam is offline   Reply With Quote

Reply

Thread Tools
Display Modes

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 Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
cyclone LES simulation prob: ERROR # 004100018 sakurabogoda CFX 15 November 6, 2012 10:07
Could you do me a favor: regarding questions in PEMFC simulation bamiao0823 FLUENT 0 March 20, 2012 03: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


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