|
[Sponsors] |
November 30, 2011, 08:32 |
BuoyantBoussinesqSimpleFoam error
|
#1 |
New Member
Begoņa
Join Date: Oct 2009
Posts: 28
Rep Power: 16 |
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. Thanks! |
|
November 30, 2011, 08:37 |
|
#2 |
New Member
Begoņa
Join Date: Oct 2009
Posts: 28
Rep Power: 16 |
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. End Thanks |
|
December 1, 2011, 08:45 |
|
#3 |
New Member
Begoņa
Join Date: Oct 2009
Posts: 28
Rep Power: 16 |
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 Regards |
|
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
[swak4Foam] GroovyBC the dynamic cousin of funkySetFields that lives on the suburb of the mesh | gschaider | OpenFOAM Community Contributions | 300 | October 29, 2014 19:00 |
c++ libraries and solver compiling | vaina74 | OpenFOAM Installation | 13 | February 3, 2012 18:43 |
[OpenFOAM] Saving ParaFoam views and case | sail | ParaView | 9 | November 25, 2011 16:46 |
DecomposePar links against liblamso0 with OpenMPI | jens_klostermann | OpenFOAM Bugs | 11 | June 28, 2007 18:51 |
user defined function | cfduser | CFX | 0 | April 29, 2006 11:58 |