CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (https://www.cfd-online.com/Forums/openfoam-solving/)
-   -   Simulation of scrubber using simpleFOAM (https://www.cfd-online.com/Forums/openfoam-solving/137853-simulation-scrubber-using-simplefoam.html)

harsha_kulkarni June 24, 2014 05:02

Simulation of scrubber using simpleFOAM
 
Hello friends

I am trying simulation single phase scrubber using simpleFOAM. I did geometry and

meshing using gambit. During the solving I found the error of sigFPE

[0] #0 Foam::error::printStack(Foam::Ostream&) at ??:?
[0] #1 Foam::sigFpe::sigHandler(int) at ??:?
[0] #2 Uninterpreted:
[0] #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) at ??:?
[0] #4 Foam::GaussSeidelSmoother::smooth(Foam::Field<doub le>&, Foam::Field<double> const&, unsigned char, int) const at ??:?
[0] #5 Foam::GAMGSolver::Vcycle(Foam::PtrList<Foam::lduMa trix::smoother> const&, Foam::Field<double>&, Foam::Field<double> const&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::PtrList<Foam::Field<double> >&, Foam::PtrList<Foam::Field<double> >&, unsigned char) const at ??:?
[0] #6 Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
[0] #7 Foam::fvMatrix<double>::solveSegregated(Foam::dict ionary const&) at ??:?
[0] #8 Foam::fvMatrix<double>::solve(Foam::dictionary const&) at ??:?
[0] #9
[0] at ??:?
[0] #10
[0] at ??:?
[0] #11 __libc_start_main in "/lib/i386-linux-gnu/libc.so.6"
[0] #12
[0] at ??:?
[ict-desktop:05913] *** Process received signal ***
[ict-desktop:05913] Signal: Floating point exception (8)
[ict-desktop:05913] Signal code: (-6)
[ict-desktop:05913] Failing at address: 0x1719
[ict-desktop:05913] [ 0] [0xb76f740c]
[ict-desktop:05913] [ 1] [0xb76f7424]
[ict-desktop:05913] [ 2] /lib/i386-linux-gnu/libc.so.6(gsignal+0x4f) [0xb50cdaff]
[ict-desktop:05913] [ 3] /home/ict/OpenFOAM/OpenFOAM-2.3.0/platforms/linuxGccDPOpt/lib/libOpenFOAM.so(_ZN4Foam6sigFpe10sigHandlerEi+0x5a) [0xb593a2aa]
[ict-desktop:05913] [ 4] [0xb76f7400]
[ict-desktop:05913] [ 5] /home/ict/OpenFOAM/OpenFOAM-2.3.0/platforms/linuxGccDPOpt/lib/libOpenFOAM.so(_ZN4Foam19GaussSeidelSmoother6smoot hERKNS_4wordERNS_5FieldIdEERKNS_9lduMatrixERKS5_RK NS_10FieldFieldIS4_dEERKNS_8UPtrListIKNS_17lduInte rfaceFieldEEEhi+0x2e8) [0xb56dcee8]
[ict-desktop:05913] [ 6] /home/ict/OpenFOAM/OpenFOAM-2.3.0/platforms/linuxGccDPOpt/lib/libOpenFOAM.so(_ZNK4Foam19GaussSeidelSmoother6smoo thERNS_5FieldIdEERKS2_hi+0x54) [0xb56dd254]
[ict-desktop:05913] [ 7] /home/ict/OpenFOAM/OpenFOAM-2.3.0/platforms/linuxGccDPOpt/lib/libOpenFOAM.so(_ZNK4Foam10GAMGSolver6VcycleERKNS_7 PtrListINS_9lduMatrix8smootherEEERNS_5FieldIdEERKS 8_S9_S9_S9_S9_S9_RNS1_IS8_EESD_h+0x1d4) [0xb56f7c54]
[ict-desktop:05913] [ 8] /home/ict/OpenFOAM/OpenFOAM-2.3.0/platforms/linuxGccDPOpt/lib/libOpenFOAM.so(_ZNK4Foam10GAMGSolver5solveERNS_5Fi eldIdEERKS2_h+0x3f4) [0xb56fadc4]
[ict-desktop:05913] [ 9] /home/ict/OpenFOAM/OpenFOAM-2.3.0/platforms/linuxGccDPOpt/lib/libfiniteVolume.so(_ZN4Foam8fvMatrixIdE15solveSegr egatedERKNS_10dictionaryE+0x14d) [0xb6ab494d]
[ict-desktop:05913] [10] simpleFoam(_ZN4Foam8fvMatrixIdE5solveERKNS_10dicti onaryE+0x146) [0x808cda6]
[ict-desktop:05913] [11] simpleFoam() [0x808d0cf]
[ict-desktop:05913] [12] simpleFoam() [0x80605f3]
[ict-desktop:05913] [13] /lib/i386-linux-gnu/libc.so.6(__libc_start_main+0xf5) [0xb50b8905]
[ict-desktop:05913] [14] simpleFoam() [0x8061327]
[ict-desktop:05913] *** End of error message ***
--------------------------------------------------------------------------
mpirun noticed that process rank 0 with PID 5913 on node ict-desktop exited on signal 8 (Floating point exception)



can you help me for that

nithishgupta June 24, 2014 05:58

Hello Harsha,

can you post your case ?!

harsha_kulkarni June 24, 2014 06:06

scruuber_case
 
1 Attachment(s)
I'm not able to attach mesh file because it's file size exceeds

nithishgupta June 24, 2014 06:30

upload it in dropbox

harsha_kulkarni June 24, 2014 06:54

I have upload it in drop box

Tushar@cfd June 25, 2014 05:45

Quote:

Hello friends

I am trying simulation single phase scrubber using simpleFOAM. I did geometry and

meshing using gambit. During the solving I found the error of sigFPE

[0] #0 Foam::error::printStack(Foam::Ostream&) at ??:?
[0] #1 Foam::sigFpe::sigHandler(int) at ??:?
[0] #2 Uninterpreted:
[0] #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) at ??:?
[0] #4 Foam::GaussSeidelSmoother::smooth(Foam::Field<doub le>&, Foam::Field<double> const&, unsigned char, int) const at ??:?
[0] #5 Foam::GAMGSolver::Vcycle(Foam::PtrList<Foam::lduMa trix::smoother> const&, Foam::Field<double>&, Foam::Field<double> const&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::PtrList<Foam::Field<double> >&, Foam::PtrList<Foam::Field<double> >&, unsigned char) const at ??:?
[0] #6 Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
[0] #7 Foam::fvMatrix<double>::solveSegregated(Foam::dict ionary const&) at ??:?
[0] #8 Foam::fvMatrix<double>::solve(Foam::dictionary const&) at ??:?
[0] #9
[0] at ??:?
[0] #10
[0] at ??:?
[0] #11 __libc_start_main in "/lib/i386-linux-gnu/libc.so.6"
[0] #12
[0] at ??:?
[ict-desktop:05913] *** Process received signal ***
[ict-desktop:05913] Signal: Floating point exception (8)
[ict-desktop:05913] Signal code: (-6)
[ict-desktop:05913] Failing at address: 0x1719
[ict-desktop:05913] [ 0] [0xb76f740c]
[ict-desktop:05913] [ 1] [0xb76f7424]
[ict-desktop:05913] [ 2] /lib/i386-linux-gnu/libc.so.6(gsignal+0x4f) [0xb50cdaff]
[ict-desktop:05913] [ 3] /home/ict/OpenFOAM/OpenFOAM-2.3.0/platforms/linuxGccDPOpt/lib/libOpenFOAM.so(_ZN4Foam6sigFpe10sigHandlerEi+0x5a) [0xb593a2aa]
[ict-desktop:05913] [ 4] [0xb76f7400]
[ict-desktop:05913] [ 5] /home/ict/OpenFOAM/OpenFOAM-2.3.0/platforms/linuxGccDPOpt/lib/libOpenFOAM.so(_ZN4Foam19GaussSeidelSmoother6smoot hERKNS_4wordERNS_5FieldIdEERKNS_9lduMatrixERKS5_RK NS_10FieldFieldIS4_dEERKNS_8UPtrListIKNS_17lduInte rfaceFieldEEEhi+0x2e8) [0xb56dcee8]
[ict-desktop:05913] [ 6] /home/ict/OpenFOAM/OpenFOAM-2.3.0/platforms/linuxGccDPOpt/lib/libOpenFOAM.so(_ZNK4Foam19GaussSeidelSmoother6smoo thERNS_5FieldIdEERKS2_hi+0x54) [0xb56dd254]
[ict-desktop:05913] [ 7] /home/ict/OpenFOAM/OpenFOAM-2.3.0/platforms/linuxGccDPOpt/lib/libOpenFOAM.so(_ZNK4Foam10GAMGSolver6VcycleERKNS_7 PtrListINS_9lduMatrix8smootherEEERNS_5FieldIdEERKS 8_S9_S9_S9_S9_S9_RNS1_IS8_EESD_h+0x1d4) [0xb56f7c54]
[ict-desktop:05913] [ 8] /home/ict/OpenFOAM/OpenFOAM-2.3.0/platforms/linuxGccDPOpt/lib/libOpenFOAM.so(_ZNK4Foam10GAMGSolver5solveERNS_5Fi eldIdEERKS2_h+0x3f4) [0xb56fadc4]
[ict-desktop:05913] [ 9] /home/ict/OpenFOAM/OpenFOAM-2.3.0/platforms/linuxGccDPOpt/lib/libfiniteVolume.so(_ZN4Foam8fvMatrixIdE15solveSegr egatedERKNS_10dictionaryE+0x14d) [0xb6ab494d]
[ict-desktop:05913] [10] simpleFoam(_ZN4Foam8fvMatrixIdE5solveERKNS_10dicti onaryE+0x146) [0x808cda6]
[ict-desktop:05913] [11] simpleFoam() [0x808d0cf]
[ict-desktop:05913] [12] simpleFoam() [0x80605f3]
[ict-desktop:05913] [13] /lib/i386-linux-gnu/libc.so.6(__libc_start_main+0xf5) [0xb50b8905]
[ict-desktop:05913] [14] simpleFoam() [0x8061327]
[ict-desktop:05913] *** End of error message ***
--------------------------------------------------------------------------
mpirun noticed that process rank 0 with PID 5913 on node ict-desktop exited on signal 8 (Floating point exception)

Check your mesh, it seems the solver failed due to failure in mesh.

Check mesh quality by running the following command (in terminal):
checkMesh

Try correcting skewed cells

-
Best Luck!


All times are GMT -4. The time now is 01:58.