[vw.cfd]

Hi,

Can anybody tell me what this error log means..

[1] #0 Foam::error:rintStack(Foam::Ostream&) in "/appl/openfoam/1.5/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM.so"
[1] #1 Foam::sigFpe::sigFpeHandler(int) in "/appl/openfoam/1.5/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM.so"
[1] #2 __restore_rt in "/lib64/tls/libc.so.6"
[1] #3 void Foam:rocessorLduInterface::compressedSend<double >(Foam::Pstream::commsTypes, Foam::UList<double> const&) const in "/appl/openfoam/1.5/OpenFOAM-1.5/lib/linux64GccDPOpt/libfiniteVolume.so"
[1] #4 Foam:rocessorGAMGInterfaceField::initInterfaceMa trixUpdate(Foam::Field<double> const&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::Field<double> const&, unsigned char, Foam::Pstream::commsTypes) const in "/appl/openfoam/1.5/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM.so"
[1] #5 Foam::lduMatrix::initMatrixInterfaces(Foam::FieldF ield<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, Foam::Field<double> const&, Foam::Field<double>&, unsigned char) const in "/appl/openfoam/1.5/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM.so"
[1] #6 Foam::lduMatrix::Amul(Foam::Field<double>&, Foam::tmp<Foam::Field<double> > const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, unsigned char) const in "/appl/openfoam/1.5/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM.so"
[1] #7 Foam::PCG::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in "/appl/openfoam/1.5/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM.so"
[1] #8 Foam::GAMGSolver::solveCoarsestLevel(Foam::Field<d ouble>&, Foam::Field<double> const&) const in "/appl/openfoam/1.5/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM.so"
[1] #9 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::PtrList<Foam::Field<double> >&, Foam::PtrList<Foam::Field<double> >&, unsigned char) const in "/appl/openfoam/1.5/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM.so"
[1] #10 Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in "/appl/openfoam/1.5/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM.so"
[1] #11 Foam::fvMatrix<double>::solve(Foam::Istream&) in "/appl/openfoam/1.5/OpenFOAM-1.5/lib/linux64GccDPOpt/libfiniteVolume.so"
[1] #12 main in "/appl/openfoam/1.5/OpenFOAM-1.5/applications/bin/linux64GccDPOpt/rasCavitatingFoam"
[1] #13 __libc_start_main in "/lib64/tls/libc.so.6"
[1] #14 Foam::regIOobject::readIfModified() in "/appl/openfoam/1.5/OpenFOAM-1.5/applications/bin/linux64GccDPOpt/rasCavitatingFoam"
[pbplnx11:18402] *** Process received signal ***
[pbplnx11:18402] Signal: Floating point exception (8)
[pbplnx11:18402] Signal code: (-6)
[pbplnx11:18402] Failing at address: 0xb1cb000047e2
[pbplnx11:18402] [ 0] /lib64/tls/libc.so.6 [0x3fd2c2e300]
[pbplnx11:18402] [ 1] /lib64/tls/libc.so.6(gsignal+0x3d) [0x3fd2c2e26d]
[pbplnx11:18402] [ 2] /lib64/tls/libc.so.6 [0x3fd2c2e300]
[pbplnx11:18402] [ 3] /appl/openfoam/1.5/OpenFOAM-1.5/lib/linux64GccDPOpt/libfiniteVolume.so(_ZNK4Foam21processorLduInterfac e14compressedSendIdEEvNS_7Pstream10commsTypesERKNS _5UListIT_EE+0x125) [0x2a95f786f5]
[pbplnx11:18402] [ 4] /appl/openfoam/1.5/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM.so(_ZNK4Foam27processorGAMGInterfaceFi eld25initInterfaceMatrixUpdateERKNS_5FieldIdEERS2_ RKNS_9lduMatrixES4_hNS_7Pstream10commsTypesE+0x47) [0x2a96c7f5e7]
[pbplnx11:18402] [ 5] /appl/openfoam/1.5/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM.so(_ZNK4Foam9lduMatrix20initMatrixInte rfacesERKNS_10FieldFieldINS_5FieldEdEERKNS_8UPtrLi stIKNS_17lduInterfaceFieldEEERKNS2_IdEERSC_h+0x357 ) [0x2a96c52dd7]
[pbplnx11:18402] [ 6] /appl/openfoam/1.5/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM.so(_ZNK4Foam9lduMatrix4AmulERNS_5Field IdEERKNS_3tmpIS2_EERKNS_10FieldFieldIS1_dEERKNS_8U PtrListIKNS_17lduInterfaceFieldEEEh+0xda) [0x2a96c514fa]
[pbplnx11:18402] [ 7] /appl/openfoam/1.5/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM.so(_ZNK4Foam3PCG5solveERNS_5FieldIdEER KS2_h+0xb3f) [0x2a96c5a57f]
[pbplnx11:18402] [ 8] /appl/openfoam/1.5/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM.so(_ZNK4Foam10GAMGSolver18solveCoarses tLevelERNS_5FieldIdEERKS2_+0x495) [0x2a96c715f5]
[pbplnx11:18402] [ 9] /appl/openfoam/1.5/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM.so(_ZNK4Foam10GAMGSolver6VcycleERKNS_7 PtrListINS_9lduMatrix8smootherEEERNS_5FieldIdEERKS 8_S9_S9_S9_RNS1_IS8_EESD_h+0xb7a) [0x2a96c72eea]
[pbplnx11:18402] [10] /appl/openfoam/1.5/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM.so(_ZNK4Foam10GAMGSolver5solveERNS_5Fi eldIdEERKS2_h+0x426) [0x2a96c744b6]
[pbplnx11:18402] [11] /appl/openfoam/1.5/OpenFOAM-1.5/lib/linux64GccDPOpt/libfiniteVolume.so(_ZN4Foam8fvMatrixIdE5solveERNS_ 7IstreamE+0x167) [0x2a96084257]
[pbplnx11:18402] [12] rasCavitatingFoam [0x41f383]
[pbplnx11:18402] [13] /lib64/tls/libc.so.6(__libc_start_main+0xdb) [0x3fd2c1c40b]
[pbplnx11:18402] [14] rasCavitatingFoam(_ZN4Foam11regIOobject14readIfMod ifiedEv+0x1a9) [0x41aa09]
[pbplnx11:18402] *** End of error message ***
mpirun noticed that job rank 0 with PID 18401 on node pbplnx11 exited on signal 15 (Terminated).


thnX!!

[gschaider]

Originally Posted by vw.cfd

Hi,

Can anybody tell me what this error log means..

[1] #0 Foam::error:rintStack(Foam::Ostream&) in "/appl/openfoam/1.5/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM.so"
[1] #1 Foam::sigFpe::sigFpeHandler(int) in "/appl/openfoam/1.5/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM.so"
[1] #2 __restore_rt in "/lib64/tls/libc.so.6"
[1] #3 void Foam:rocessorLduInterface::compressedSend<double >(Foam::Pstream::commsTypes, Foam::UList<double> const&) const in "/appl/openfoam/1.5/OpenFOAM-1.5/lib/linux64GccDPOpt/libfiniteVolume.so"
[1] #4 Foam:rocessorGAMGInterfaceField::initInterfaceMa trixUpdate(Foam::Field<double> const&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::Field<double> const&, unsigned char, Foam::Pstream::commsTypes) const in "/appl/openfoam/1.5/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM.so"
[1] #5 Foam::lduMatrix::initMatrixInterfaces(Foam::FieldF ield<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, Foam::Field<double> const&, Foam::Field<double>&, unsigned char) const in "/appl/openfoam/1.5/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM.so"
[1] #6 Foam::lduMatrix::Amul(Foam::Field<double>&, Foam::tmp<Foam::Field<double> > const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, unsigned char) const in "/appl/openfoam/1.5/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM.so"
[1] #7 Foam::PCG::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in "/appl/openfoam/1.5/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM.so"
[1] #8 Foam::GAMGSolver::solveCoarsestLevel(Foam::Field<d ouble>&, Foam::Field<double> const&) const in "/appl/openfoam/1.5/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM.so"
[1] #9

Not much except that something seems to be going wrong with the AMG-Solver during a parallel run. I'm not really a specialist on that. Just some questions and pointers:
- Is this a vanilla-1.5 (the way it was released a year ago, without any bugfixes from 1.5.x)? If yes: the problem may or may not have been already solved in that version
- Does this happen during the first timestep and which field is being solved? Have you tried a different solver for that field?
- Look elsewhere on the board for discussion on "floatTransfer" and how to switch it of (it seems to cause problems)
- You could always compile a Debug-Version to get a more defailed view of what is happening, but I dobt that this helps in this case

Bernhard

[henrik]

Dear VW.CFD,

I agree with Bernhard's analysis and would like to add that the traceback can be quite misleading in the sense that the simulation might have gone wrong a long time ago. Please post the log of the last three iterations.

Henrik

[vw.cfd]

Hi,

Thank you Bernhard and Henrik for your feedbacks. As you suggested i looked at the field it is beeing solved, in the present case this is for the turbulence parameter 'omega'. But, earlier i had the same error when the solver was solving for the 'pressure'. It is not the first time-step it is after a certain random number of time-steps.

As requested i'm uploading the log file for the previous few time-steps. I am also attaching a file of the 'checkMesh' output for your considerations.


phiv Courant Number mean: 0.000380924308 max: 0.371157546 acoustic max: 5.82495376
deltaT = 2.35594662e-08
Time = 1.24556e-05
DILUPBiCG: Solving for rho, Initial residual = 0.0023438356, Final residual = 8.24817937e-10, No Iterations 3
max-min rho: 853.281786 -33.1438517
max-min gamma: 0.999998956 0
DILUPBiCG: Solving for Ux, Initial residual = 3.74124376e-05, Final residual = 5.3823599e-09, No Iterations 2
DILUPBiCG: Solving for Uy, Initial residual = 4.21000707e-05, Final residual = 5.18736189e-10, No Iterations 3
DILUPBiCG: Solving for Uz, Initial residual = 0.000266227609, Final residual = 1.17634205e-10, No Iterations 3
GAMG: Solving for p, Initial residual = 5.25665764e-05, Final residual = 3.82712155e-09, No Iterations 4
GAMG: Solving for p, Initial residual = 6.80353456e-06, Final residual = 3.21031189e-09, No Iterations 3
max-min p: 20561545.3 -8170302.45
max(U) 343.420096
max-min rho: 853.280749 0.001
max-min gamma: 0.999998956 0
GAMG: Solving for p, Initial residual = 2.41274035e-05, Final residual = 2.15307058e-09, No Iterations 4
GAMG: Solving for p, Initial residual = 3.36319336e-06, Final residual = 1.98127757e-09, No Iterations 3
max-min p: 20561546 -7177911.83
max(U) 344.737214
max-min rho: 853.280749 0.001
max-min gamma: 0.999998956 0
DILUPBiCG: Solving for omega, Initial residual = 3.23259431e-05, Final residual = 4.56596357e-10, No Iterations 2
DILUPBiCG: Solving for k, Initial residual = 0.000661195431, Final residual = 1.90450155e-10, No Iterations 2
ExecutionTime = 754.84 s ClockTime = 764 s


phiv Courant Number mean: 0.000430361638 max: 0.426024778 acoustic max: 6.57311296
deltaT = 2.62311583e-08
Time = 1.24818e-05
DILUPBiCG: Solving for rho, Initial residual = 0.00284606542, Final residual = 1.45435808e-09, No Iterations 3
max-min rho: 853.253409 -59.3930822
max-min gamma: 0.999998956 0
DILUPBiCG: Solving for Ux, Initial residual = 4.17829742e-05, Final residual = 7.9040332e-10, No Iterations 2
DILUPBiCG: Solving for Uy, Initial residual = 4.64698452e-05, Final residual = 5.89524353e-09, No Iterations 2
DILUPBiCG: Solving for Uz, Initial residual = 0.000296500127, Final residual = 7.69338574e-10, No Iterations 3
GAMG: Solving for p, Initial residual = 6.87716678e-05, Final residual = 8.14339585e-09, No Iterations 4
GAMG: Solving for p, Initial residual = 9.72950493e-06, Final residual = 6.44368233e-09, No Iterations 3
max-min p: 20504419.3 -8459374.5
max(U) 323.578958
max-min rho: 853.252186 0.001
max-min gamma: 0.999998956 0
GAMG: Solving for p, Initial residual = 3.18981244e-05, Final residual = 4.99989942e-09, No Iterations 4
GAMG: Solving for p, Initial residual = 4.59918881e-06, Final residual = 3.09157315e-09, No Iterations 3
max-min p: 20504420.7 -6459866.19
max(U) 319.938852
max-min rho: 853.252186 0.001
max-min gamma: 0.999998956 0
DILUPBiCG: Solving for omega, Initial residual = 3.61815272e-05, Final residual = 7.6978693e-10, No Iterations 2
DILUPBiCG: Solving for k, Initial residual = 0.000731965859, Final residual = 2.86899069e-09, No Iterations 2
ExecutionTime = 759.38 s ClockTime = 769 s


phiv Courant Number mean: 0.000479822444 max: 0.417382535 acoustic max: 7.31851752
deltaT = 2.92812929e-08
Time = 1.25111e-05
DILUPBiCG: Solving for rho, Initial residual = 0.00352277146, Final residual = 4.72396138e-09, No Iterations 3
max-min rho: 853.23056 -28.5267983
max-min gamma: 0.999998956 0
DILUPBiCG: Solving for Ux, Initial residual = 4.67640442e-05, Final residual = 5.60401146e-09, No Iterations 2
DILUPBiCG: Solving for Uy, Initial residual = 5.15667913e-05, Final residual = 1.78882879e-10, No Iterations 4
DILUPBiCG: Solving for Uz, Initial residual = 0.000331103855, Final residual = 5.43882975e-09, No Iterations 3
GAMG: Solving for p, Initial residual = 9.47929981e-05, Final residual = 3.90432687e-09, No Iterations 5
GAMG: Solving for p, Initial residual = 1.40666394e-05, Final residual = 2.92658361e-09, No Iterations 4
max-min p: 20458084.5 -10492652.3
max(U) 296.184552
max-min rho: 853.229018 0.001
max-min gamma: 0.999998956 0
GAMG: Solving for p, Initial residual = 4.36100877e-05, Final residual = 3.72547745e-09, No Iterations 5
GAMG: Solving for p, Initial residual = 6.42895012e-06, Final residual = 7.29378262e-09, No Iterations 3
max-min p: 20458086.2 -9626283.2
max(U) 334.641901
max-min rho: 853.229019 0.001
max-min gamma: 0.999998956 0
DILUPBiCG: Solving for omega, Initial residual = 4.05266552e-05, Final residual = 1.67471828e-10, No Iterations 2
bounding omega, min: -869.709373 max: 42128310.2 average: 189138.438
DILUPBiCG: Solving for k, Initial residual = 0.000811794653, Final residual = 4.00167814e-09, No Iterations 2
bounding k, min: -1.14326837 max: 95.0198859 average: 10.1747313
ExecutionTime = 764.21 s ClockTime = 773 s


phiv Courant Number mean: 0.00053635358 max: 0.47651801 acoustic max: 8.16950793
deltaT = 3.03525598e-08
Time = 1.25414e-05
DILUPBiCG: Solving for rho, Initial residual = 0.0041395039, Final residual = 1.33256684e-10, No Iterations 5
max-min rho: 864.551168 -65.0711859
max-min gamma: 0.999998956 0
DILUPBiCG: Solving for Ux, Initial residual = 4.87134868e-05, Final residual = 2.31588453e-09, No Iterations 3
DILUPBiCG: Solving for Uy, Initial residual = 5.37202582e-05, Final residual = 2.82646994e-10, No Iterations 4
DILUPBiCG: Solving for Uz, Initial residual = 0.000343392475, Final residual = 4.29377918e-13, No Iterations 6
GAMG: Solving for p, Initial residual = 0.000111569338, Final residual = 7.32231451e-09, No Iterations 6
GAMG: Solving for p, Initial residual = 1.70657195e-05, Final residual = 4.13042209e-09, No Iterations 4
max-min p: 20406783.2 -10849833.8
max(U) 449.864622
max-min rho: 853.203368 0.001
max-min gamma: 0.999998956 0
GAMG: Solving for p, Initial residual = 5.21365559e-05, Final residual = 5.13079653e-09, No Iterations 6
GAMG: Solving for p, Initial residual = 8.28114896e-06, Final residual = 9.41193553e-09, No Iterations 3
max-min p: 20406785.1 -9701693.22
max(U) 419.448286
max-min rho: 853.203369 0.001
max-min gamma: 0.999998956 0
DILUPBiCG: Solving for omega, Initial residual = 4.22738141e-05, Final residual = 1.18288164e-09, No Iterations 2
bounding omega, min: -1217516.55 max: 42250552.8 average: 189394.238
DILUPBiCG: Solving for k, Initial residual = 0.000835537576, Final residual = 2.05331868e-11, No Iterations 3
ExecutionTime = 769.42 s ClockTime = 779 s


phiv Courant Number mean: 0.000556802097 max: 0.664929262 acoustic max: 8.46839237
deltaT = 2.27644198e-08
Time = 1.25642e-05
DILUPBiCG: Solving for rho, Initial residual = 0.00354665096, Final residual = 1.25658862e-09, No Iterations 3
max-min rho: 869.402441 -24.5276701
max-min gamma: 0.999998956 0
DILUPBiCG: Solving for Ux, Initial residual = 3.66717904e-05, Final residual = 8.39855408e-09, No Iterations 2
DILUPBiCG: Solving for Uy, Initial residual = 4.03500128e-05, Final residual = 4.88353931e-10, No Iterations 3
DILUPBiCG: Solving for Uz, Initial residual = 0.000257824366, Final residual = 2.79214584e-09, No Iterations 3
GAMG: Solving for p, Initial residual = 7.22796351e-05, Final residual = 5.23238124e-09, No Iterations 4
GAMG: Solving for p, Initial residual = 9.20285287e-06, Final residual = 4.20702159e-09, No Iterations 3
max-min p: 20383897.8 -6772311.36
max(U) 346.543414
max-min rho: 853.191925 0.001
max-min gamma: 0.999998956 0
GAMG: Solving for p, Initial residual = 2.73424891e-05, Final residual = 2.51530559e-09, No Iterations 4
GAMG: Solving for p, Initial residual = 3.95192647e-06, Final residual = 2.1913065e-09, No Iterations 3
max-min p: 20383898.6 -5244772.98
max(U) 339.871797
max-min rho: 853.191925 0.001
max-min gamma: 0.999998956 0
DILUPBiCG: Solving for omega, Initial residual = 3.16516414e-05, Final residual = 4.14998638e-10, No Iterations 2
bounding omega, min: -57693.0699 max: 41123338.9 average: 189607.333
DILUPBiCG: Solving for k, Initial residual = 0.000622284076, Final residual = 7.82122104e-10, No Iterations 2
ExecutionTime = 774.03 s ClockTime = 783 s


phiv Courant Number mean: 0.000418043338 max: 0.359475769 acoustic max: 6.35129427
deltaT = 2.59126906e-08
Time = 1.25901e-05
DILUPBiCG: Solving for rho, Initial residual = 0.00435851381, Final residual = 2.97328137e-09, No Iterations 3
max-min rho: 853.18138 -18.9434031
max-min gamma: 0.999998956 0
DILUPBiCG: Solving for Ux, Initial residual = 4.17675039e-05, Final residual = 8.76460186e-10, No Iterations 4
DILUPBiCG: Solving for Uy, Initial residual = 4.56364255e-05, Final residual = 1.52260376e-09, No Iterations 5
DILUPBiCG: Solving for Uz, Initial residual = 0.000293445154, Final residual = 2.66404975e-10, No Iterations 4
GAMG: Solving for p, Initial residual = 9.57640968e-05, Final residual = 1.13728303e-07, No Iterations 1000
GAMG: Solving for p, Initial residual = 1.33085803e-05, Final residual = 4.90273077e-08, No Iterations 1000
max-min p: 20360590.4 -20260388.1
max(U) 925.025696
max-min rho: 853.180271 0.001
max-min gamma: 0.999998956 0
1 additional process aborted (not shown)


Please post your feedbacks and any criticisms..


Thank you

varun

[henrik]

Dear VW.CFD,

the density got unbounded (ie. negative) and then the pressure solver blew up (1000 iterations). There are many possible reasons for this (setup, boundary conditions, numerics of top-level-solver), but your chances to find the problem in the core library where you traceback emerged are very slim.

Henrik

[vw.cfd]

Dear Henrik,


Thank you again for your views. I have a few more questions since, i am a new user to CFD. I would like you to suggest a numerical scheme for my problem. Let me define my problem to you. I am simulating a two phase flow with one phase as 'oil' and the other as 'air'. There inlet boundary condition is mass inflow with a pressure outflow as exit BC. The geometry is a 20° sector of a cylinder( i am attaching a picture of the mesh file for visualisation). The two side walls are assumed as symmetry planes.

The numerics i have employed is the same as provided for the tutorial test case 'tutorials/rasCavitatingFoam/throttle' . I will repeat the schemes for you i'm using TVD schemes for the Divergence terms and an Unbounded , second order conservative scheme for the Laplacian terms. Turbelence is modelled by the Wilcox K-omega + SST scheme.

Would it be better to use an upwind scheme?? Please give me your suggestions.


Thank you,

Varun
Thank you,

Varun

[henrik]

Dear Varun,

did run the case single phase incompressible, yet?

If not, this is the first thing to do!

Henrik