tomislav_maric

Snippet:

In case you are dealing with sharp fields, there are multiple options to choose from. Usually, the literature recommends the "least squares approach", or linear regression. In order to get somewhere with this, the least squares points should pick up the skew cells , not just the normal own-nei surrounding cells in OpenFOAM. Still, there is a workaround that works fine:

#1) IDW interpolation of cell values to mesh points
#2) create a surface*field
#3) average the...

The primitiveMesh::cellCells returns with a cross cell stencil of 6 cells (the cell in question is excluded). This computation comes from the basis of OpenFOAM: face addressing. Although it is very efficient, this kind of computation introduces errros e.g. in the gradient reconstruction of sharp fields.

Try to generate a field that has a jump in a single layer of cells, whose iso surface is a plane, and try to compute the normal field of this plane from the fraction gradient using...

If you are dealing with some delicate computations that require you to investigate the values of fields in some cells (I use this for debugging), there is a way to do this in Paraview:

View->Selection Inspector

Create Selection

Selection Type: Thresholds

Check the box: Invert selection

Scroll down to the Display Style subwindow and select Cell Label

Check the box: Visible

Choose: "p" for pressure...

If you are visualising a bunch of .vtk files that were generated separately from the standard fields in OpenFOAM, and if you are a lucky person like me at this point (your time steps are constant), there exists a simple way to visualize the files alongside OpenFOAM data:

1) Open the OpenFOAM case as usual.
2) Open the .vtk files. The files have to be named appropriately to be read as time-dependent data: someName-001.vtk, someName-002.vtk, guess what is the next name. ...