bug in dnsFOAM
I discovered a bug in dnsFOAM!
More precisely, the last line in kShellIntegration.C is wrong.
It should be changed from
return graph("E(k)", "k", "E(k)", Ek1D, k1D);
return graph("E(k)", "k", "E(k)", k1D, Ek1D);
otherwise the following manipulation and the resulting
energy spectra do not make any sense!
I would also like to mention that it would be very helpful if a reference to the used forcing scheme
would be given somewhere in the code or maybe a tiny little description of the forcing parameters sigma and alpha.
This would spare the users of dnsFoam a lot headache while trying to fiddle with these parameters!
Thanks for the bug-report, I will push the fix into OpenFOAM-1.6.x shortly.
You might find what you need to know in this paper or in the references:
author = "Fureby, C. and Tabor, G. and Weller, H.G. and Gosman, A.D.",
title = "A Comparative Study of Sub-Grid Scale Models in Homogeneous
journal = "Phys. Fluids",
year = 1997,
volume = 9,
number = 5,
pages = "1416-1429",
month = "May"
|All times are GMT -4. The time now is 20:46.|