[File]

Hi to everyone,

I have a problem with running a icoFoam on my own mesh. It is a simple model with one inlet and one outlet. The intention is to make a pressure driven flow as I use a very small mesh 1.4x1.4x1.4 mm pressure difference is 5 Pa

Most of the files are copied from icoFoam tutorial case (cavity) of course I adopted them to my case.

Mesh is generated in Avizo and then imported to enGrid and from there to openFoam. I must point that the same mesh works in Fluent without any problem.

Case loads all the input files and starts the caluculation but after the first step it brakes. Here is the log.

PHP Code:

/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.2.0                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : 2.2.0-5be49240882f
Exec   : icoFoam
Date   : Jun 18 2013
Time   : 16:43:15
Host   : "zeus"
PID    : 7484
Case   : /home/filip/OpenFOAM/proba/moje_2
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Disallowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0

Reading transportProperties

Reading field p

Reading field U

Reading/calculating face flux field phi


Starting time loop

Time = 0.0005

Courant Number mean: 0 max: 0
DILUPBiCG:  Solving for Ux, Initial residual = 1, Final residual = 9.00268e-06, No Iterations 606
DILUPBiCG:  Solving for Uy, Initial residual = 1, Final residual = 9.63275e-06, No Iterations 600
DILUPBiCG:  Solving for Uz, Initial residual = 1, Final residual = 9.85118e-06, No Iterations 609
#0  Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam220/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1  Foam::sigFpe::sigHandler(int) in "/opt/openfoam220/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2   in "/lib/x86_64-linux-gnu/libc.so.6"
#3  Foam::DICPreconditioner::calcReciprocalD(Foam::Field<double>&, Foam::lduMatrix const&) in "/opt/openfoam220/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#4  Foam::DICPreconditioner::DICPreconditioner(Foam::lduMatrix::solver const&, Foam::dictionary const&) in "/opt/openfoam220/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#5  Foam::lduMatrix::preconditioner::addsymMatrixConstructorToTable<Foam::DICPreconditioner>::New(Foam::lduMatrix::solver const&, Foam::dictionary const&) in "/opt/openfoam220/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#6  Foam::lduMatrix::preconditioner::New(Foam::lduMatrix::solver const&, Foam::dictionary const&) in "/opt/openfoam220/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#7  Foam::PCG::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in "/opt/openfoam220/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#8  Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&) in "/opt/openfoam220/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#9  Foam::fvMatrix<double>::solve(Foam::dictionary const&) in "/opt/openfoam220/platforms/linux64GccDPOpt/bin/icoFoam"
#10  
 at icoFoam.C:0
#11  
 in "/opt/openfoam220/platforms/linux64GccDPOpt/bin/icoFoam"
#12  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#13  
 in "/opt/openfoam220/platforms/linux64GccDPOpt/bin/icoFoam"
Floating point exception (core dumped) 


Here are my input files if somebody wants to take a look
0:
https://dl.dropboxusercontent.com/u/122319307/0.tar.gz
constant:
https://dl.dropboxusercontent.com/u/...onstant.tar.gz
system:
https://dl.dropboxusercontent.com/u/.../system.tar.gz

Also I runned checkMesh and I had some warning which I don't know how to fix so here is the output:

PHP Code:

/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.2.0                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : 2.2.0-5be49240882f
Exec   : checkMesh
Date   : Jun 18 2013
Time   : 16:57:33
Host   : "zeus"
PID    : 7664
Case   : /home/filip/OpenFOAM/proba/moje_2
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Disallowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create polyMesh for time = 0

Time = 0

Mesh stats
    points:           95290
    faces:            947476
    internal faces:   875052
    cells:            455632
    faces per cell:   4
    boundary patches: 3
    point zones:      0
    face zones:       0
    cell zones:       0

Overall number of cells of each type:
    hexahedra:     0
    prisms:        0
    wedges:        0
    pyramids:      0
    tet wedges:    0
    tetrahedra:    455632
    polyhedra:     0

Checking topology...
    Boundary definition OK.
    Cell to face addressing OK.
    Point usage OK.
    Upper triangular ordering OK.
    Face vertices OK.
   *Number of regions: 4
    The mesh has multiple regions which are not connected by any face.
  <<Writing region information to "0/cellToRegion"

Checking patch topology for multiply connected surfaces...
    Patch               Faces    Points   Surface topology                  
    Inlet               1966     1143     ok (non-closed singly connected)  
    Outlet              1962     1132     ok (non-closed singly connected)  
    wall                68496    34516    multiply connected (shared edge)  
  <<Writing 8 conflicting points to set nonManifoldPoints

Checking geometry...
    Overall domain bounding box (-5.2703e-07 -5.57384e-07 -5.82047e-07) (0.00141052 0.0014106 0.00141059)
    Mesh (non-empty, non-wedge) directions (1 1 1)
    Mesh (non-empty) directions (1 1 1)
    Boundary openness (1.4271e-17 4.92582e-17 -2.06799e-17) OK.
    Max cell openness = 2.7799e-16 OK.
    Max aspect ratio = 15.936 OK.
    Minimum face area = 1.72258e-11. Maximum face area = 1.07184e-09.  Face area magnitudes OK.
    Min volume = 5.78296e-17. Max volume = 1.01736e-14.  Total volume = 6.68254e-10.  Cell volumes OK.
    Mesh non-orthogonality Max: 78.3513 average: 21.5589
   *Number of severely non-orthogonal faces: 6.
    Non-orthogonality check OK.
  <<Writing 6 non-orthogonal faces to set nonOrthoFaces
    Face pyramids OK.
    Max skewness = 3.02936 OK.
    Coupled point location match (average 0) OK.

Mesh OK.

End 


tnx