CFD Online URL
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

lowReOneEqEddy sgsModel

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

Reply
 
LinkBack Thread Tools Display Modes
Old   September 5, 2009, 14:49
Default lowReOneEqEddy sgsModel
  #1
Member
 
cosimo bianchini
Join Date: Mar 2009
Location: Florence, Tuscany, Italy
Posts: 88
Rep Power: 7
cosimobianchini is on a distinguished road
Send a message via Skype™ to cosimobianchini
Hi foamers,
has any of you ever worked with lowReOneEqEddy?

I could see from the comments in the header file that this model implements
nuSgs this way:

nuSgsHiRe = ck*sqrt(k)*delta;
nuSgs = (nu/beta)*(1 - exp(-beta*nuSgsHiRe/nu));

What is this correction supposed to be for?
Is it supposed to act as a damping function for near wall treatment or the lowRe simply stands for low free stream velocity, and the exponential function is just used to limit the maximum nuSgs to 100 times the laminar viscosity?
Does anybody have a reference, to share or to suggest, for this model?

Anyhow looking at the source code, there appears to be an incosistency between the comments and the code itself:

void lowReOneEqEddy::updateSubGridScaleFields()
{
// High Re eddy viscosity
muSgs_ = ck_*rho()*sqrt(k_)*delta();

// low Re no corrected eddy viscosity
muSgs_ -= (mu()/beta_)*(scalar(1) - exp(-beta_*muSgs_/mu()));
muSgs_.correctBoundaryConditions();

alphaSgs_ = muSgs_/Prt_;
alphaSgs_.correctBoundaryConditions();
}

In particular shouldn't the -= operator be just = ?

Thanks a lot for sharing any comment,
Cosimo
__________________
Cosimo Bianchini

Energy Engineering Department "S. Stecco"
University of Florence
Via di S.Marta, 3
50139 Florence - ITALY
Tel: +39 055 4796575
Fax: +39 055 4796342
Mob: +39 320 9460153
e-mail: cosimo.bianchini@htc.de.unifi.it
URL: www.htc.de.unifi.it
cosimobianchini 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



All times are GMT -4. The time now is 10:27.