[Legacy]

Hello Foamers! I am trying to simulate flow past tandem cylinders in OpenFOAM. I'm using pisoFoam (tired pimpleFoam too) and I am getting same error. I generated the mesh via Fidelity Pointwise and exported it to OpenFOAM.

My simulation was crashing so I decided to create a simple rectangular block instead of two cylinder in blockMesh and the simulation finished without any errors. I assume it is something mesh related. I created the same mesh (block instead of cylinder in Pointwise and my simulation crashes).

I installed the OpenFOAM in debug mode but couldn't figure out what is going. Here is the output of the simulation. I also attached the 0 and system directory.

I get floating point exemption error so I know it's division by zero or something but could not figure it out.

 at /work/myCluster/OpenFOAM-v2212/src/OSspecific/POSIX/printStack/printStack.C:280
[15] #1  Foam::sigFpe::sigHandler(int)- at /work/myCluster/OpenFOAM-v2212/src/OSspecific/POSIX/signals/sigFpe.C:126
[15] #2  ? in /lib64/libpthread.so.0
[15] #3  Foam::GaussSeidelSmoother::smooth(Foam::word const&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::Field<double> const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, unsigned char, int)f at /work/myCluster/OpenFOAM-v2212/src/OpenFOAM/matrices/lduMatrix/smoothers/GaussSeidel/GaussSeidelSmoother.C:161
[15] #4  Foam::GaussSeidelSmoother::smooth(Foam::Field<double>&, Foam::Field<double> const&, unsigned char, int) const  at /work/myCluster/OpenFOAM-v2212/src/OpenFOAM/matrices/lduMatrix/smoothers/GaussSeidel/GaussSeidelSmoother.C:184 (discriminator 1)
[15] #5  Foam::smoothSolvertail: ::solve(Fwarning: following standard input indefinitely is ineffectiveoam:
:Field<double>&, Foam::Field<double> const&, unsigned char) const at /work/myCluster/OpenFOAM-v2212/src/OpenFOAM/matrices/lduMatrix/solvers/smoothSolver/smoothSolver.C:190
[15] #6  Foam::fvMatrix<Foam::Vector<double> >::solveSegregated(Foam::dictionary const&) at /work/myCluster/OpenFOAM-v2212/src/finiteVolume/lnInclude/fvMatrixSolve.C:220 (discriminator 9)
[15] #7  Foam::fvMatrix<Foam::Vector<double> >::solveSegregatedOrCoupled(Foam::dictionary const&) at /work/myCluster/OpenFOAM-v2212/src/finiteVolume/lnInclude/fvMatrixSolve.C:92
[15] #8  Foam::fvMesh::solve(Foam::fvMatrix<Foam::Vector<double> >&, Foam::dictionary const&) const at /work/myCluster/OpenFOAM-v2212/src/finiteVolume/fvMesh/fvMesh.C:595
[15] #9  ? at /work/myCluster/OpenFOAM-v2212/src/finiteVolume/lnInclude/fvMatrixSolve.C:324
[15] #10  ? at /work/myCluster/OpenFOAM-v2212/src/finiteVolume/lnInclude/fvMatrixSolve.C:346
[15] #11  ? at /work/myCluster/OpenFOAM-v2212/src/finiteVolume/lnInclude/fvMatrix.C:2009
[15] #12  ? at /work/myCluster/OpenFOAM-v2212/applications/solvers/incompressible/pisoFoam/UEqn.H:20 (discriminator 14)
[15] #13  __libc_start_main in /lib64/libc.so.6
[15] #14  ? in /work/myCluster/OpenFOAM-v2212/platforms/linux64GccDPInt32Debug/bin/pisoFoam
[smic002:14413:0:14413] Caught signal 8 (Floating point exception: tkill(2) or tgkill(2))
==== backtrace (tid:  14413) ====
 0 0x000000000004d455 ucs_debug_print_backtrace()  ???:0
 1 0x000000000000f4fb raise()  ???:0
 2 0x0000000000c4901c Foam::sigFpe::sigHandler()  /work/myCluster/OpenFOAM-v2212/src/OSspecific/POSIX/signals/sigFpe.C:126
 3 0x00000000009e2205 Foam::GaussSeidelSmoother::smooth()  /work/myCluster/OpenFOAM-v2212/src/OpenFOAM/matrices/lduMatrix/smoothers/GaussSeidel/GaussSeidelSmoother.C:161
 4 0x00000000009e2386 Foam::GaussSeidelSmoother::smooth()  /work/myCluster/OpenFOAM-v2212/src/OpenFOAM/matrices/lduMatrix/smoothers/GaussSeidel/GaussSeidelSmoother.C:184
 5 0x00000000009dad45 Foam::smoothSolver::solve()  /work/myCluster/OpenFOAM-v2212/src/OpenFOAM/matrices/lduMatrix/solvers/smoothSolver/smoothSolver.C:182
 6 0x0000000000d6802e Foam::fvMatrix<Foam::Vector<double> >::solveSegregated()  /work/myCluster/OpenFOAM-v2212/src/finiteVolume/lnInclude/fvMatrixSolve.C:220
 7 0x0000000000d60080 Foam::fvMatrix<Foam::Vector<double> >::solveSegregatedOrCoupled()  /work/myCluster/OpenFOAM-v2212/src/finiteVolume/lnInclude/fvMatrixSolve.C:92
 8 0x0000000000d59331 Foam::fvMesh::solve()  /work/myCluster/OpenFOAM-v2212/src/finiteVolume/fvMesh/fvMesh.C:594
 9 0x0000000000471d65 Foam::fvMatrix<Foam::Vector<double> >::solve()  /work/myCluster/OpenFOAM-v2212/src/finiteVolume/lnInclude/fvMatrixSolve.C:323
10 0x00000000004654f2 Foam::fvMatrix<Foam::Vector<double> >::solve()  /work/myCluster/OpenFOAM-v2212/src/finiteVolume/lnInclude/fvMatrixSolve.C:345
11 0x0000000000459796 Foam::solve<Foam::Vector<double> >()  /work/myCluster/OpenFOAM-v2212/src/finiteVolume/lnInclude/fvMatrix.C:2007
12 0x0000000000450d63 main()  /work/myCluster/OpenFOAM-v2212/applications/solvers/incompressible/pisoFoam/UEqn.H:20
13 0x0000000000022555 __libc_start_main()  ???:0
14 0x000000000044dd89 _start()  ???:0

Checkmesh

Mesh stats
    points:           443072
    faces:            1275388
    internal faces:   1222580
    cells:            416328
    faces per cell:   6
    boundary patches: 7
    point zones:      0
    face zones:       0
    cell zones:       0

Overall number of cells of each type:
    hexahedra:     416328
    prisms:        0
    wedges:        0
    pyramids:      0
    tet wedges:    0
    tetrahedra:    0
    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: 1 (OK).

Checking patch topology for multiply connected surfaces...
    Patch               Faces    Points   Surface topology
    bottom              18924    19264    ok (non-closed singly connected)
    cylinder            1936     2024     ok (non-closed singly connected)
    inlet               2156     2277     ok (non-closed singly connected)
    outlet              2156     2277     ok (non-closed singly connected)
    sideLeft            4356     4577     ok (non-closed singly connected)
    sideRight           4356     4577     ok (non-closed singly connected)
    top                 18924    19264    ok (non-closed singly connected)

Checking faceZone topology for multiply connected surfaces...
    No faceZones found.

Checking basic cellZone addressing...
    No cellZones found.

Checking basic pointZone addressing...
    No pointZones found.

Checking geometry...
    Overall domain bounding box (-0.5 -0.15 0) (0.5 0.15 0.1)
    Mesh has 3 geometric (non-empty/wedge) directions (1 1 1)
    Mesh has 3 solution (non-empty) directions (1 1 1)
    Boundary openness (-2.42231e-16 -6.91153e-16 -1.18655e-13) OK.
    Max cell openness = 3.20423e-16 OK.
    Max aspect ratio = 5.28697 OK.
    Minimum face area = 2.70684e-06. Maximum face area = 4.5476e-05.  Face area magnitudes OK.
    Min volume = 1.11426e-08. Max volume = 1.56e-07.  Total volume = 0.0293222.  Cell volumes OK.
    Mesh non-orthogonality Max: 0 average: 0
    Non-orthogonality check OK.
d    Face pyramids OK.                                                                                                                                                                                           Max skewness = 2.19941e-13 OK.
    Coupled point location match (average 0) OK.

Mesh OK.

End

Thanks!

[alundilong]

please upload the constant folder as well, I can run it on my machine.