CFD Online Logo CFD Online URL
Home > Forums > OpenFOAM Bugs

Possible bug in nutRoughWall

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

LinkBack Thread Tools Display Modes
Old   October 1, 2009, 18:17
Default Possible bug in nutRoughWall
New Member
John O\'Sullivan
Join Date: Mar 2009
Location: Auckland, New Zealand
Posts: 7
Rep Power: 10
jposunz is on a distinguished road
Hi all,

I think there may be a bug in the rough wall function. Its been discussed a little in the thread:

But basically nutw is set to 0 every time step. This means that

scalar limitingNutw = max(nutw[faceI], nuw[faceI]);

always returns limitingNutw = nuw[faceI].

This in turn means that

nutw[faceI] =
*(yPlus*kappa_/log(max(Edash*yPlus, 1+1e-4)) - 1),
), 0.5*limitingNutw

always returns nutw[faceI] = 2*limitingNutw

As a result the simulation is subjected to a Ks much less than specified.

In the test case attached the velocity should end up close to:

u = u_star/kappa*log((z+z0)/z0)

where in this case u_star = 0.625137, kappa = 0.4 and z0 = 0.01.

I've fixed the problem buy adding the following lines which evaluate the limitingNutw using the values nut in the cell next to the wall.

const tmp<volScalarField> toldnut = rasModel.nut();
const volScalarField& oldnut = toldnut();

scalar limitingNutw = max(oldnut[faceCellI], nuw[faceI]);

Attached Files
File Type: zip (17.9 KB, 12 views)

Last edited by jposunz; October 1, 2009 at 18:40.
jposunz is offline   Reply With Quote

Old   October 2, 2009, 02:56
Senior Member
Join Date: Mar 2009
Posts: 854
Rep Power: 15
henry is on a distinguished road
Yes, we fixed this a couple of days ago and are testing the fix at the moment and will push it into 1.6.x shortly.

henry 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

Similar Threads
Thread Thread Starter Forum Replies Last Post
Maybe a bug in CFX12.0 BETA? frank CFX 4 November 26, 2008 19:08
Bug in twoPhaseEulerFoam wallfunctions alberto OpenFOAM Bugs 1 February 9, 2007 15:15
Bug reports Mattijs Janssens (Mattijs) OpenFOAM 0 January 10, 2005 11:05
Win XP hot fix bug!!! Emilio Siemens 1 March 27, 2003 13:06
Forum y2k Bug Jonas Larsson Main CFD Forum 1 January 5, 2000 11:22

All times are GMT -4. The time now is 15:55.