|August 25, 2011, 06:56||
Use of face centered differences in gradient computation
Join Date: Aug 2011
Location: Bergen, Norway
Posts: 2Rep Power: 0
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?
|August 25, 2011, 10:34||
Steven van Haren
Join Date: Aug 2010
Location: The Netherlands
Posts: 149Rep Power: 5
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:
volScalarField G("RASModel::G", nut_*2*magSqr(symm(fvc::grad(U_))));
Is this what you were referring to?
|Thread||Thread Starter||Forum||Replies||Last Post|
|Import netgen mesh to OpenFOAM||hsieh||Open Source Meshers: Gmsh, Netgen, CGNS, ...||32||September 13, 2011 06:50|
|BlockMeshmergePatchPairs||hjasak||OpenFOAM Native Meshers: blockMesh||11||August 15, 2008 08:36|
|fluent add additional zones for the mesh file||SSL||FLUENT||2||January 26, 2008 12:55|
|Axisymmetrical mesh||Rasmus Gjesing (Gjesing)||OpenFOAM Native Meshers: blockMesh||10||April 2, 2007 15:00|
|Trimmed cell and embedded refinement mesh conversion issues||michele||OpenFOAM Other Meshers: ICEM, Star, Ansys, Pointwise, GridPro, Ansa, ...||2||July 15, 2005 05:15|