CFD Online Discussion Forums

CFD Online Discussion Forums (
-   OpenFOAM Running, Solving & CFD (
-   -   BuoyantBoussinesqSimpleFoam error (

bego November 30, 2011 08:32

BuoyantBoussinesqSimpleFoam error
Hi all,

I solved a problem buoyant air flow inside a cone with a diffuser inside. Without the diffuser solves the problem perfectly, but I get the diffuser insert floating-point exception and ends the iteration. I tried using pressure PCG, and now I have changed to GAMG do not know if I'm wrong in using these methods. As for the other parameters I'm using PBiCG and DILU as preconditioner.
I don't know where is the error.

bego November 30, 2011 08:37

I forget say that I had a mesh mistake, my checkMesh log is the following file

Create time

Create polyMesh for time = 0

Time = 0

Mesh stats
points: 188872
faces: 463519
internal faces: 403803
cells: 139587
boundary patches: 6
point zones: 0
face zones: 0
cell zones: 0

Overall number of cells of each type:
hexahedra: 114624
prisms: 4832
wedges: 8
pyramids: 0
tet wedges: 55
tetrahedra: 0
polyhedra: 20068

Checking topology...
Boundary definition OK.
Cell to face addressing OK.
Point usage OK.
Upper triangular ordering OK.
Face vertices OK.
Number of regions: 1 (OK).

Checking patch topology for multiply connected surfaces ...
Patch Faces Points Surface topology
paredExterior 3232 3427 ok (non-closed singly connected)
salida 8236 8482 ok (non-closed singly connected)
copa115cerradoYarriba.stl_copa115cerradoY-arriba4846 5693 ok (non-closed singly connected)
horno_cono_nuevo.stl_horno_cono_nuevo35497 39761 ok (non-closed singly connected)
pisadorabiertohorno.stl_pisadorabiertohorno7701 8251 ok (non-closed singly connected)
tapa_8mm.stl_tapa_8mm204 229 ok (non-closed singly connected)

Checking geometry...
Overall domain bounding box (-0.1 -0.24 -0.14) (0.1 0.05 0.14)
Mesh (non-empty, non-wedge) directions (1 1 1)
Mesh (non-empty) directions (1 1 1)
Boundary openness (-1.3848e-16 -1.52287e-16 2.37733e-16) OK.
Max cell openness = 3.17735e-16 OK.
Max aspect ratio = 28.6561 OK.
Minumum face area = 2.78817e-07. Maximum face area = 6.65257e-05. Face area magnitudes OK.
Min volume = 4.8752e-10. Max volume = 3.52085e-07. Total volume = 0.0146593. Cell volumes OK.
Mesh non-orthogonality Max: 63.2908 average: 12.3099
Non-orthogonality check OK.
Face pyramids OK.
***Max skewness = 8.02469, 8 highly skew faces detected which may impair the quality of the results
<<Writing 8 skew faces to set skewFaces

Failed 1 mesh checks.



bego December 1, 2011 08:45

Hi foamers,
I have solved my problem checking maxBoundarySkewness and maxInternalSkewness to 0.5.
I wish help somebody with this,
Thanks for reed my post

All times are GMT -4. The time now is 21:59.