CFD Online Logo CFD Online URL
Home > Forums > OpenFOAM Post-Processing

Problems with YPlusRAS and wallShearStress

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

Like Tree23Likes

LinkBack Thread Tools Display Modes
Old   May 5, 2013, 08:42
Senior Member
immortality's Avatar
Join Date: Oct 2012
Location: Iran
Posts: 2,209
Rep Power: 20
immortality is on a distinguished road
in mutkWallFunctionFvPatchScalarField the equation uses for y-plus is related to y-star not y-plus

y-plus=rho*u_star*y/mu which u*=aqrt(taw_w/rho)
y-star is related to turbulent kinetic energy while y-plus relates to friction velocity(that relates to wall skin friction in its own)

in the code of mut
forAll(mutw, faceI)     {         label faceCellI = patch().faceCells()[faceI];          scalar yPlus =             Cmu25*y[faceI]*sqrt(k[faceCellI])/(muw[faceI]/rhow[faceI]);          if (yPlus > yPlusLam_)         {             mutw[faceI] = muw[faceI]*(yPlus*kappa_/log(E_*yPlus) - 1);         }     }
Injustice Anywhere is a Threat for Justice Everywhere.Martin Luther King.
To Be or Not To Be,Thats the Question!
The Only Stupid Question Is the One that Goes Unasked.
immortality is offline   Reply With Quote

Old   May 4, 2016, 05:15
Default Correction necessary on Niklas's code
New Member
vikramaditya gaonkar
Join Date: Jun 2012
Location: Stuttgart, Germany
Posts: 18
Rep Power: 7
vikramaditya91 is on a distinguished road
I believe I am very late to the party. I am currently using the OpenFOAM 2.4.0 and I have to say that they still have not made any changes to the code, even though I believe it would not be a involve a lot of coding. There was a post on the official openFOAM page, in case you might want to have a look, here it is

Well i was using the code Niklas had made, but I believe it requires a correction. Of course it needed to be edited slightly for the OpenFOAM 2.4 but that is not the correction I am talking about.
When it calculates the sqrt(viscosity*gradient of velocity), it has to consider the sum of kinematic viscosity+turbulent viscosity , but he has only used the turbulent viscosity (mut.boundaryField()[patchi]).
Practically, the turbulent viscoisty can be ignored and we are only left with the kinematic viscosisty. Very close to the wall i.e y+<3, the turbulent viscosity is negligible since it lies in the viscous laminar sub layer where there are negligible turbulent effects. The line I am talking about would be 3rd in the code attached which is from his Niklas's code.
yPlus.boundaryField()[patchi] = 
            y[patchi]           * sqrt              (            
   *  mag(U.boundaryField()[patchi].snGrad())       
   /  rho.boundaryField()[patchi]                 )     
      /  (RASModel->mu().boundaryField()[patchi]
I would like to hear what you guys think about it. have a look at this for reference. Check out the section near equation (21) in this page

Last edited by vikramaditya91; May 4, 2016 at 06:38. Reason: Grammatik
vikramaditya91 is offline   Reply With Quote


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 06:01.