Use of face centered differences in gradient computation
According to the recent work of Richrds and Norris (J. Wind Eng Ind.Aerodyn 99 (2011) 257-266), the near ground spike in turbulent kinetic enery in the k-epsilon
model originates from an inconsistency in the discretization. This error can be removed by computing all gradients in terms of face centered differences.
Can this shift be achieved by any of the options in the fvSchemes file?
Are you sure this is not already being done?
I think discritization (53) from the paper you refer to is not possible in an unstructered code.
Furthermore, in kEpsilon.C:
Is this what you were referring to?
I have asked this question before
and I think there is no such scheme available. Please correct me if I am wrong.
|All times are GMT -4. The time now is 20:34.|