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

Calculation of yPlus in MRF simulations

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

LinkBack Thread Tools Search this Thread Display Modes
Old   February 17, 2015, 04:27
Default Calculation of yPlus in MRF simulations
New Member
Matej Andrejasic
Join Date: Sep 2010
Posts: 2
Rep Power: 0
MatejAndrej is on a distinguished road
Dear All,

I'm simulating airflow passing a propeller using MRF. I'm using incompressible solver simpleFoam and kOmegaSST turbulence model with nutkWallFunction as a wall function at nut. I don't have problems with simulations, it converges beautifully to nice results, but I would like to ask how do I calculate yPlus in my situation?

Since it is High-Re type of simulation I ran yPlusRAS utility, but resulting yPlus is not ok even though I recalculated U by adding \omega x R (rotating term).

After a little bit of research in doxygen I found that since I'm using nutkWallFunction the utility yPlusRAS consequently uses equations: u_{\tau}=C_{\mu}^{1/4}\sqrt{k} together with y+ definition y+=y u_{tau}/\nu. Since turbulent kinetic energy k is being calculated in rotating frame of reference it apparently doesn't result in correct y+ calculation?!?

Then I changed kOmegaSST in constant/RASProperties to SpalartAllmaras and I changed all nutkWallFunction in nut file to nutUWallFunction in the last simulation time folder. Then I ran yPlusRAS again and this time it calculated "logical" y+ distribution.

nutUWallFunction uses an iterative calculation in order to calculate y+ described in: (standard wall function).

Since I ran simulation with nutkWallFunction and calculated yPlus with nutUWallFunction I don't know what to think about these results?!?

How do you calculate yPlus in MRF?

Thank you very much for your help.

MatejAndrej is offline   Reply With Quote

Old   April 28, 2015, 12:01
Senior Member
M. Montero
Join Date: Mar 2009
Location: Madrid
Posts: 134
Rep Power: 14
be_inspired is on a distinguished road
Dear MatejAndrej,

Excuse for using this same topic but I think you can help me the little I need.
I am trying to use the same approach that yours to simulate a wind turbine using MRF approach. All goes ok until I activate turbulence ( I am trying to detect where the problem is).
For the walls inside the rotating region I have put nutkWallFunction, KqRWallFunction and epsilonWallFunction.
I have also tried nutUSpaldingWallFunction with no luck

Is it needed to change in some way the wall function for MRF? The k and epsilon values at wall increase a lot and then the simulation crash.

I am using OF2.3.x and kEpsilon model ( the mesh is ready for y+ around 100)
be_inspired is offline   Reply With Quote


mrf, nutkwallfunction, nutuwallfunction, yplus

Thread Tools Search this Thread
Search this Thread:

Advanced Search
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 Off
Pingbacks are On
Refbacks are On

Similar Threads
Thread Thread Starter Forum Replies Last Post
[Other] Contribution a new utility: refine wall layer mesh based on yPlus field lakeat OpenFOAM Community Contributions 57 February 1, 2015 08:25
Solver Yplus of Automatic wall treatment justjhy CFX 2 May 4, 2013 07:19
Wrong calculation of yPlus for low-Re models bastil OpenFOAM Bugs 4 February 9, 2010 07:39
which yplus (SST) sanchezz CFX 17 January 11, 2010 04:45
DIVERGENCE with MRF calculation in FLUENT FLUENT USER Main CFD Forum 3 June 28, 1999 15:44

All times are GMT -4. The time now is 13:43.