Description: On certain cases
 

Description:
On certain cases, simpleFoam solves for U and p but for some reason does not print out the solution for velocity (i.e. BICCG solver residuals) in stdout or the redirected log file. The case in question happens to be a non-conformal 2D case. Mesh generated using blockMesh. Not sure whether this problem applies to 3D cases too?

Solver/Application: Out of the box simpleFoam -> No modifications.

Source file:
simpleFoam.C ?

Testcase:
Can be downloaded from here -> http://www.ualberta.ca/~madhavan/sample.tar.gz

Platform:
x86_64 GNU/Linux Distro: Scientific Linux 4.1 (basically an RHEL clone).

Version:
OpenFOAM version 1.3

Notes:
Related discussion thread: http://www.cfd-online.com/OpenFOAM_D...tml?1176418785

Thanks for the bug report. Th
 

Thanks for the bug report. The issue is associated with the tolerance used in checking the dimensionality of the case. Apparently the complexity of your mesh is causing an accumulation of round-off error which means the empty patches are not exactly flat. To fix this problem edit

src/OpenFOAM/meshes/polyMesh/polyMesh.C

and in polyMesh::calcDirections() change SMALL to 1e-6 and recompile the OpenFOAM library.

alternatively download the 1.4 pack I uploaded to sourceforge today which already includes this fix.

Henry

Brilliant! Will upgrade to 1.4
 

Brilliant! Will upgrade to 1.4 as soon as possible. Thanks!!