CFD Online Logo CFD Online URL
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

Need some help debugging OpenFOAM error

Register Blogs Members List Search Today's Posts Mark Forums Read

LinkBack Thread Tools Search this Thread Display Modes
Old   June 26, 2023, 18:54
Default Need some help debugging OpenFOAM error
New Member
Join Date: Jul 2021
Posts: 11
Rep Power: 4
Legacy is on a distinguished road
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/
[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/
[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
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.

Attached Files
File Type: zip (3.5 KB, 1 views)
File Type: zip (4.1 KB, 1 views)
Legacy is offline   Reply With Quote

Old   August 21, 2023, 00:04
yijin Mao
Join Date: May 2010
Location: Columbia, MO
Posts: 59
Rep Power: 15
alundilong is on a distinguished road
please upload the constant folder as well, I can run it on my machine.
alundilong is offline   Reply With Quote


debugging, fpe, openfoam

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On

Similar Threads
Thread Thread Starter Forum Replies Last Post
[OpenFOAM] ParaView command in Foam-extend-4.1 mitu_94 ParaView 0 March 4, 2021 14:46
long error when using make-install SU2_AD. tomp1993 SU2 Installation 3 March 17, 2018 07:25
Undeclared Identifier Errof UDF SteveGoat Fluent UDF and Scheme Programming 7 October 15, 2014 08:11
Ansys Fluent 13.0 UDF compilation problem in Window XP (32 bit) Yogini Fluent UDF and Scheme Programming 7 October 3, 2012 08:24
[swak4Foam] groovyBC: problems compiling: "flex: not found" and "undefined reference to ..." sega OpenFOAM Community Contributions 12 February 17, 2010 10:30

All times are GMT -4. The time now is 08:53.