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

ce and ck in ‘oneEqEddy’

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

Like Tree7Likes
  • 2 Post By cfdonline2mohsen
  • 3 Post By zxj160
  • 2 Post By zxj160

Reply
 
LinkBack Thread Tools Display Modes
Old   August 10, 2013, 10:54
Default ce and ck in ‘oneEqEddy’
  #1
Senior Member
 
Jian Zhong
Join Date: Feb 2012
Location: Birmingham
Posts: 109
Rep Power: 6
zxj160 is on a distinguished road
Dear foamers,

I want to modify the value of ce and ck in the LES model (i.e. oneEqEddy) and implement a new turbulent model. I found the declaration of ck_ and its value in the 'oneEqEddy.C'. But where can I find the declaration of ce?

My best regards,
Jian
zxj160 is offline   Reply With Quote

Old   August 11, 2013, 17:10
Default
  #2
Senior Member
 
cfdonline2mohsen's Avatar
 
Mohsen KiaMansouri
Join Date: Jan 2010
Location: CFD Lab
Posts: 118
Rep Power: 8
cfdonline2mohsen is on a distinguished road
Dear Jian
I remember that I changed the value of Cs in the standard smagorinsky model once to Cs=0.12 by using the relationship between Cs and ce and ck. (I can't remember now but you can find this relationship in turbulence books)
I remember that only one of the coefficients (ce or ck) were defined in the .c code (probably ck_ that is defined in the code) and the other coefficient (probably ce) were defined in one of the #.H that exists at the beginning of the code.
because Ce will be used by the other turbulence models; so changing its value will affect the other models.
so it is better to derive the new value for ck for your desired Cs by using the relationship between Cs and ce and ck and change cK in this code rather than changing Ce.
zxj160 and ScarFace like this.
__________________
“If you have an apple and I have an apple and we exchange these apples then you and I will still each have one apple. But if you have an idea and I have an idea and we exchange these ideas, then each of us will have two ideas.”
cfdonline2mohsen is offline   Reply With Quote

Old   August 13, 2013, 07:37
Default
  #3
Senior Member
 
Jian Zhong
Join Date: Feb 2012
Location: Birmingham
Posts: 109
Rep Power: 6
zxj160 is on a distinguished road
Dear Kia,

I find a way to modified these value without changing the source code. But I am not sure whether it is right.

Just add the following lines in the 'LESProperties' file under the constant directory.

....

LESModel oneEqEddy;

delta cubeRootVol;

printCoeffs on;

cubeRootVolCoeffs
{
deltaCoeff 1;
}

oneEqEddyCoeffs
{
ck 0.08;
ce 1.05;
}

....


My best regards,
Jian
zxj160 is offline   Reply With Quote

Old   August 13, 2013, 08:49
Default
  #4
Senior Member
 
cfdonline2mohsen's Avatar
 
Mohsen KiaMansouri
Join Date: Jan 2010
Location: CFD Lab
Posts: 118
Rep Power: 8
cfdonline2mohsen is on a distinguished road
Thanks Jian
I did not notice that!!
that was tricky
at the beginning of the run; OpenFOAM writes the solver that you use and the turbulence models coefficients at your computer's screen.
does it print this new values of ce and ck that is specified in oneEqEddyCoeffs?
__________________
“If you have an apple and I have an apple and we exchange these apples then you and I will still each have one apple. But if you have an idea and I have an idea and we exchange these ideas, then each of us will have two ideas.”
cfdonline2mohsen is offline   Reply With Quote

Old   August 13, 2013, 08:53
Default
  #5
Senior Member
 
Jian Zhong
Join Date: Feb 2012
Location: Birmingham
Posts: 109
Rep Power: 6
zxj160 is on a distinguished road
Dear Kia,

Yes, I think so. The new values are outputted at the beginning of the run. So I guess it should work properly.

My best regards,
Jian
cfdonline2mohsen and ScarFace like this.
zxj160 is offline   Reply With Quote

Reply

Tags
oneeeeddy

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



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