|
[Sponsors] |
March 11, 2011, 13:01 |
"LRR Turbulent Model" Problems
|
#1 |
Senior Member
|
Dear experts,
By changing KEpsilon to LRR in constants/RASProprties OpenFOAM works fine apparently but the problem is in this model: k and epsilon which is set for inlet has a sharp decrease toward internal domain for example in my case which is a backward facing step it reach from k=7.5 to 0.0001, also same for epsilon. I think it is a reason which it is not possible to reach true answer. Settings are same as pitzDaily tutorial. i used pisoFoam. What is the problem? Any suggestion will be appreciated. Regards. |
|
March 12, 2011, 04:43 |
|
#2 |
Senior Member
Claus Meister
Join Date: Aug 2009
Location: Wiesbaden, Germany
Posts: 241
Rep Power: 18 |
Hard to day: Did u set the inlet conditions for R correctly, meaning everywhere zero apart from the main diagonal and on the k=... I haven't the conditions formula on me, but you can loop up in Versteeg.
|
|
March 12, 2011, 05:46 |
|
#4 |
Senior Member
Claus Meister
Join Date: Aug 2009
Location: Wiesbaden, Germany
Posts: 241
Rep Power: 18 |
You must set the BC for R correctly:
Compute k=3/2*(U_ref*Ti)^2 For U_ref take the velocity at the inlet, look in "U", set Ti = 0.05 (for other cases you muss adjusted it). If you have concrete value for k than set for R at the inlet: type fixedValue; value uniform (k/2 0 0 k/2 0 k); as you can see the k-values are divided by 2; one is untouched, this one points in flow direction (here z-axis). Since you integrate a epsilon equation with these, you must set the BC suitable. Somewhere in the user guide or programmer's guide is an entry, you gotta go for it. |
|
March 12, 2011, 07:55 |
|
#5 | |
Senior Member
|
Quote:
Thanks for your suggestion. i have calculated k and epsilon. k=7.5 and epsilon=500. But by seeing incompressible/simpleFoam/pitzdaily tutorial decide to use (0 0 0 0 0 0) for inlet R. By changing inlet R to what you told answers are seemed to be better. Two questions: 1- How can R be calculated to what value you told? i mean if x direction be streamwise or if it be 3D and ... 2- Is T=(a b c d e f) means: a d e d b f e f c or it has another format for symmetric tensors. i searched in user guide and programmers guide but nothing found about it. |
||
March 12, 2011, 11:36 |
|
#6 |
Senior Member
Claus Meister
Join Date: Aug 2009
Location: Wiesbaden, Germany
Posts: 241
Rep Power: 18 |
(1) Basically, when you consider a canonical reference frame where the inlet is pointing into the x-direction than you must set:
R=(k 0 0 k/2 0 k/2), i.e. in flow direction k and into the others k/2. The reason why k appears here is due to the relationship trace(R)=2*k (look Versteeg). R is a symmetric tensor; to reduce memory consumption the symmetry property is exploited, i.e. 6 entries instead of 9, so (2) is right! R is called the Reynolds Stress Tensor and appearce by reynolds averaging procedure. My suggestion: Try to get Versteeg. Many times posted here and by simply entering it in google; you find in google books, definitely. Cheers and good luck |
|
March 12, 2011, 12:31 |
|
#7 |
Senior Member
|
I have Versteeg. thanks for your suggestion.
if trace of reynolds stress tensor should be 2k and if your told R could be true, format of 6 elements tensor instead of 9 elements couldn't be what i told. it should be: T=(a b c d e f) means: a b c b d e c e f Isn't it? |
|
March 12, 2011, 13:58 |
|
#8 |
Senior Member
Claus Meister
Join Date: Aug 2009
Location: Wiesbaden, Germany
Posts: 241
Rep Power: 18 |
You got it, foamer!
|
|
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Laminar vs Turbulent Navier-Stokes | truman | Main CFD Forum | 8 | July 10, 2017 08:20 |
Stability problem due to turbulent dispersion force in a subcooled boiling model | Edy | OpenFOAM | 7 | August 10, 2011 13:00 |
Turbulent Kinetic Energy | Olga | FLUENT | 2 | October 11, 2002 16:05 |
Problem of Turbulent Viscosity Ratio Limited | David Yang | FLUENT | 3 | June 3, 2002 07:13 |
Turbulent viscosity | Christian Holm | Main CFD Forum | 4 | June 23, 2001 23:04 |