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

Problem calculate Y in compressible turbulence model

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

Reply
 
LinkBack Thread Tools Display Modes
Old   May 31, 2006, 11:00
Default Hi all, May be It's a studp
  #1
Member
 
Luca M.
Join Date: Mar 2009
Location: Luzern, Switzerland
Posts: 59
Rep Power: 8
luca is on a distinguished road
Hi all,

May be It's a studpid question but I've noticed in file "turbulence.C" in compressible turbulence models that OF calucalte Y+ like:

Yp = pow(Cmu.value(), 0.25)
*y_[patchNo]
*sqrt(k()().boundaryField()[patchNo].patchInternalField())
/(
mu().boundaryField()[patchNo].patchInternalField()
*rho_.boundaryField()[patchNo]<---
);

I think that it's wrong because I need the Kinematic viscosity that is defined as


mu/rho. (I think that mu is refereed to Dynamic viscosity)

and so the right way to calculate Y+ is

Yp = pow(Cmu.value(), 0.25)
*y_[patchNo]
*sqrt(k()().boundaryField()[patchNo].patchInternalField())
/(
mu().boundaryField()[patchNo].patchInternalField()
/rho_.boundaryField()[patchNo]<---
);



I'm wrong ?!

Luca
luca is offline   Reply With Quote

Old   May 31, 2006, 11:41
Default Hi, In the calculation of Y
  #2
Member
 
Pierre Le Fur
Join Date: Mar 2009
Location: UK
Posts: 60
Rep Power: 8
pierre is on a distinguished road
Hi,

In the calculation of Yplus in wallViscosity.H and WallFunctions.H for compressible flow codes, it is written:
"
scalar yPlus =
Cmu25*turbulenceModel::y_[patchi][facei]*sqrt(k_[faceCelli])
/(muw[facei]/rhow[facei]);
"
which is as you suggest, so I guess you might be right. Thought that member function defined in turbulence.C is not used anywhere it seems.

Pierre
pierre is offline   Reply With Quote

Old   May 31, 2006, 11:58
Default Hi Pierre, So It's a bug..
  #3
Member
 
Luca M.
Join Date: Mar 2009
Location: Luzern, Switzerland
Posts: 59
Rep Power: 8
luca is on a distinguished road
Hi Pierre,

So It's a bug..

I use the function Yp in turbulence.C to check y+ for compressible turbulence model that I use. In detail I build the utility starting from "checkYPlus" incompressible utility where, I think, it uses Yp function.

So how do you suggest me to build my utility to check y+ without Yp function in turbulence compressible model?!

Thanks in advance

Luca
luca is offline   Reply With Quote

Old   May 31, 2006, 12:19
Default why don't you calculate Yplus
  #4
Member
 
Pierre Le Fur
Join Date: Mar 2009
Location: UK
Posts: 60
Rep Power: 8
pierre is on a distinguished road
why don't you calculate Yplus directly as you wrote it above without using Yp function?

Pierre
pierre is offline   Reply With Quote

Old   June 1, 2006, 05:16
Default Thanks Luca, it was a bug and
  #5
Super Moderator
 
Mattijs Janssens
Join Date: Mar 2009
Posts: 1,416
Rep Power: 16
mattijs is on a distinguished road
Thanks Luca, it was a bug and has been fixed. Any more of these? If so please report them on the bugs list so we can keep track.
mattijs is offline   Reply With Quote

Old   June 1, 2006, 05:53
Default Ok,Mattijs!! Thanks for the
  #6
Member
 
Luca M.
Join Date: Mar 2009
Location: Luzern, Switzerland
Posts: 59
Rep Power: 8
luca is on a distinguished road
Ok,Mattijs!!

Thanks for the suggestion, and support.

Luca
luca 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
Compressible Turbulence Model Implementations shaun OpenFOAM Running, Solving & CFD 4 October 21, 2009 11:55
compressible turbulence model shuo Main CFD Forum 3 December 22, 2007 10:51
Turbulence model problem? PattiMichelle Phoenics 0 July 11, 2007 19:35
Problem with LaunderSharma compressible turbulence model wrong formulation luca OpenFOAM Running, Solving & CFD 4 July 3, 2006 06:24
LES turbulence model problem-HELP James Willie FLUENT 2 August 23, 2005 04:54


All times are GMT -4. The time now is 09:18.