May 20, 2018, 06:55
|
Foam fatal IO error
|
#1
|
New Member
Albert
Join Date: May 2018
Posts: 4
Rep Power: 7
|
Hello everyone, I am new in Open foam and I'm just facing some problems that I am unable to solve. Once I created my mesh and I have done gmshTofoam and set all the values in the folders 0 constant and system, when I run checkMesh, an error appears but I don't know why, I leave you here the code:
Quote:
Create time
Create mesh for time = 0
Time = 0
Mesh stats
points: 167260
internal points: 0
faces: 334670
internal faces: 166123
cells: 83680
faces per cell: 5.98462
boundary patches: 8
point zones: 0
face zones: 0
cell zones: 1
Overall number of cells of each type:
hexahedra: 82393
prisms: 1287
wedges: 0
pyramids: 0
tet wedges: 0
tetrahedra: 0
polyhedra: 0
Checking topology...
Boundary definition OK.
***Total number of faces on empty patches is not divisible by the number of cells in the mesh. Hence this mesh is not 1D or 2D.
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
back 83680 83630 ok (non-closed singly connected)
front 82958 82982 ok (non-closed singly connected)
inlet 27 56 ok (non-closed singly connected)
up 47 96 ok (non-closed singly connected)
outlet 87 176 ok (non-closed singly connected)
down 57 116 ok (non-closed singly connected)
wall 969 1938 ok (non-closed singly connected)
defaultFaces 722 800 ok (non-closed singly connected)
Checking faceZone topology for multiply connected surfaces...
No faceZones found.
Checking basic cellZone addressing...
CellZone Cells Points BoundingBox
fluid 83680 167260 (-2 -1.5 0) (3 1.5 0.1)
Checking geometry...
Overall domain bounding box (-2 -1.5 0) (3 1.5 0.1)
Mesh has 2 geometric (non-empty/wedge) directions (1 1 0)
Mesh has 2 solution (non-empty) directions (1 1 0)
All edges aligned with or perpendicular to non-empty directions.
Boundary openness (-2.6261144e-18 -3.6988763e-18 -1.0586381e-14) OK.
Max cell openness = 6.2447076e-16 OK.
Max aspect ratio = 24.088342 OK.
Minimum face area = 3.3701894e-08. Maximum face area = 0.050357489. Face area magnitudes OK.
Min volume = 3.3701894e-09. Max volume = 0.0050357489. Total volume = 1.4917998. Cell volumes OK.
Mesh non-orthogonality Max: 80.129444 average: 12.836373
*Number of severely non-orthogonal (> 70 degrees) faces: 8.
Non-orthogonality check OK.
<<Writing 8 non-orthogonal faces to set nonOrthoFaces
Face pyramids OK.
***Max skewness = 5.1180369, 15 highly skew faces detected which may impair the quality of the results
<<Writing 15 skew faces to set skewFaces
Coupled point location match (average 0) OK.
Failed 1 mesh checks.
End
|
On the other hand, when I run icoFoam, another error appears and I can't contiune with my simulation:
Quote:
Create time
Create mesh for time = 0
PISO: Operating solver in PISO mode
Reading transportProperties
Reading field p
--> FOAM FATAL IO ERROR:
Cannot find patchField entry for defaultFaces
file: /home/albert/Documents/upc/Mecanicadefluids/Projecte/Laminar/0/p.boundaryField
From function void Foam::GeometricField<Type, PatchField, GeoMesh>::Boundary::readField(const Foam:imensionedField<TypeR, GeoMesh>&, const Foam::dictionary&) [with Type = double; PatchField = Foam::fvPatchField; GeoMesh = Foam::volMesh]
in file /home/albert/OpenFOAM/OpenFOAM-v1712/src/OpenFOAM/lnInclude/GeometricBoundaryField.C at line 191.
FOAM exiting
|
Thank you very much,
Albert
|
|
|