View Issue Details

IDProjectCategoryView StatusLast Update
0000085OpenFOAMBugpublic2010-11-23 14:06
ReporterJuho Assigned Touser2 
PrioritynormalSeveritymajorReproducibilityN/A
Status resolvedResolutionfixed 
Summary0000085: kappatJayatillekeWallFunction: kappat calculation
DescriptionIn kappatJayatillekeWallFunction, when yPlus > yPlusTherm, the kappat is calculated from:

scalar kt = nu*(yPlus/(Prt_/kappa_*log(E_*yPlusTherm) + P) - 1/Pr)

I think the calculation should be:

scalar kt = nu*(yPlus/(Prt_*(log(E_*yPlus)/kappa_ + P)) - 1/Pr)

The modified calculation gives mesh independent results and matches results obtained with the formulation used in the alphaSgsJayatillekeWallFunction. Do you know a good reference for the formulations?
TagsNo tags attached.

Activities

user2

2010-11-23 14:06

  ~0000134

Thanks for the bug report - the code was indeed incorrect - fixed by commit:

    05bdccba11142a2796b25a480943f2341abe89bd

Issue History

Date Modified Username Field Change
2010-11-23 09:33 Juho New Issue
2010-11-23 09:52 user2 Assigned To => user2
2010-11-23 09:52 user2 Status new => assigned
2010-11-23 14:06 user2 Note Added: 0000134
2010-11-23 14:06 user2 Status assigned => resolved
2010-11-23 14:06 user2 Fixed in Version => 1.7.x
2010-11-23 14:06 user2 Resolution open => fixed