CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM (http://www.cfd-online.com/Forums/openfoam/)
-   -   print stack error (http://www.cfd-online.com/Forums/openfoam/75934-print-stack-error.html)

CFT May 10, 2010 05:25

print stack error
 
Hi all

I am running a rans simulation, however in the first step i get the following print stack error

Starting time loop

Time = 0.0001

DILUPBiCG: Solving for Ux, Initial residual = 1, Final residual = 5.57e-06, No Iterations 1
DILUPBiCG: Solving for Uy, Initial residual = 1, Final residual = 5.7e-06, No Iterations 1
DILUPBiCG: Solving for Uz, Initial residual = 0, Final residual = 0, No Iterations 0
DICPCG: Solving for p, Initial residual = 1, Final residual = 0.00966, No Iterations 41
time step continuity errors : sum local = 0.0718, global = 0.0712, cumulative = 0.0712
#0 Foam::error::printStack(Foam::Ostream&) in "/home/steven/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libOpenFOAM.so"
#1 Foam::sigFpe::sigFpeHandler(int) in "/home/steven/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libOpenFOAM.so"
#2 Uninterpreted: [0xffffe400]
#3 Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) in "/home/steven/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libOpenFOAM.so"
#4 void Foam::divide<Foam::fvPatchField, Foam::volMesh>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in "/home/steven/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libincompressibleRASModels.so"
#5 Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::operator/<Foam::fvPatchField, Foam::volMesh>(Foam::tmp<Foam::GeometricField<doub le, Foam::fvPatchField, Foam::volMesh> > const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in "/home/steven/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libincompressibleRASModels.so"
#6 Foam::incompressible::RASModels::RNGkEpsilon::corr ect() in "/home/steven/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libincompressibleRASModels.so"
#7 main in "/home/steven/OpenFOAM/OpenFOAM-1.5/applications/bin/linuxGccDPOpt/simpleFoam"
#8 __libc_start_main in "/lib/libc.so.6"
#9 Foam::regIOobject::readIfModified() in "/home/steven/OpenFOAM/OpenFOAM-1.5/applications/bin/linuxGccDPOpt/simpleFoam"
Floating point exception

any suggestions will be appreciated


CFT

gschaider May 11, 2010 07:13

Quote:

Originally Posted by CFT (Post 258147)
#3 Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) in "/home/steven/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libOpenFOAM.so"
#4 void Foam::divide<Foam::fvPatchField, Foam::volMesh>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in "/home/steven/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libincompressibleRASModels.so"

It's the old "k or epsilon zero somewhere (internal or boundary)"-problem

It gets asked for every second week (on the other weeks it's "pressure 0 in compressible solver")

CFT May 11, 2010 22:41

problem solved thanks for help

immortality January 9, 2013 14:29

please help Bernhard.what does make this error?
Code:

Time = 2.1641e-06

diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG:  Solving for Ux, Initial residual = 0.003904074052, Final residual = 2.005454408e-07, No Iterations 1
DILUPBiCG:  Solving for Uy, Initial residual = 0.003938431751, Final residual = 1.762807354e-07, No Iterations 1
DILUPBiCG:  Solving for h, Initial residual = 0.002436468979, Final residual = 1.507492028e-07, No Iterations 1
DICPCG:  Solving for p, Initial residual = 0.004935048413, Final residual = 6.446599636e-09, No Iterations 1
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 0.2036037114, global = -0.2036037114, cumulative = -57.13916748
rho max/min : 2 0.6061772378
DICPCG:  Solving for p, Initial residual = 4.573658842e-06, Final residual = 1.200479079e-11, No Iterations 1
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 0.203603714, global = -0.203603714, cumulative = -57.34277119
rho max/min : 2 0.6061337361
DILUPBiCG:  Solving for omega, Initial residual = 0.0006665378688, Final residual = 4.023732082e-08, No Iterations 1
DILUPBiCG:  Solving for k, Initial residual = 0.002765522059, Final residual = 3.385841907e-07, No Iterations 1
ExecutionTime = 18.45 s  ClockTime = 18 s

Courant Number mean: 6.378037208e-05 max: 0.09956912126
deltaT = 7.530598113e-09
Time = 2.17163e-06

diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG:  Solving for Ux, Initial residual = 0.003889209539, Final residual = 2.070157988e-07, No Iterations 1
DILUPBiCG:  Solving for Uy, Initial residual = 0.003914010891, Final residual = 1.713622134e-07, No Iterations 1
DILUPBiCG:  Solving for h, Initial residual = 0.002426461002, Final residual = 1.448322811e-07, No Iterations 1
DICPCG:  Solving for p, Initial residual = 0.004904283451, Final residual = 6.405767558e-09, No Iterations 1
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 0.2036078435, global = -0.2036078435, cumulative = -57.54637903
rho max/min : 2 0.6059568975
DICPCG:  Solving for p, Initial residual = 4.317143549e-06, Final residual = 1.144976892e-11, No Iterations 1
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 0.2036078461, global = -0.2036078461, cumulative = -57.74998688
rho max/min : 2 0.6059136717
DILUPBiCG:  Solving for omega, Initial residual = 0.0006663067204, Final residual = 4.144346607e-08, No Iterations 1
DILUPBiCG:  Solving for k, Initial residual = 0.002749633138, Final residual = 3.045808594e-07, No Iterations 1
ExecutionTime = 18.57 s  ClockTime = 18 s

Courant Number mean: 6.422498588e-05 max: 0.09991615161
deltaT = 7.530598113e-09
Time = 2.17916e-06

diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG:  Solving for Ux, Initial residual = 0.003887642308, Final residual = 2.368225887e-07, No Iterations 1
DILUPBiCG:  Solving for Uy, Initial residual = 0.003887776118, Final residual = 1.508081681e-07, No Iterations 1
DILUPBiCG:  Solving for h, Initial residual = 0.002418835722, Final residual = 1.453036584e-07, No Iterations 1
DICPCG:  Solving for p, Initial residual = 0.004871840268, Final residual = 6.358455224e-09, No Iterations 1
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 0.2036118998, global = -0.2036118998, cumulative = -57.95359878
rho max/min : 2 0.6057834604
DICPCG:  Solving for p, Initial residual = 4.046935923e-06, Final residual = 1.087937439e-11, No Iterations 1
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 0.2036119023, global = -0.2036119023, cumulative = -58.15721068
rho max/min : 2 0.605740617
DILUPBiCG:  Solving for omega, Initial residual = 0.0006657876225, Final residual = 4.224787986e-08, No Iterations 1
DILUPBiCG:  Solving for k, Initial residual = 0.002735420898, Final residual = 2.776874161e-07, No Iterations 1
ExecutionTime = 18.69 s  ClockTime = 19 s

Courant Number mean: 6.466661485e-05 max: 0.1002717423
deltaT = 7.46213813e-09
Time = 2.18663e-06

diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG:  Solving for Ux, Initial residual = 0.003843906995, Final residual = 2.437691487e-07, No Iterations 1
DILUPBiCG:  Solving for Uy, Initial residual = 0.003845199722, Final residual = 1.536313099e-07, No Iterations 1
DILUPBiCG:  Solving for h, Initial residual = 0.002384056736, Final residual = 1.424960386e-07, No Iterations 1
DICPCG:  Solving for p, Initial residual = 0.004789183743, Final residual = 6.011477139e-09, No Iterations 1
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 0.2036158846, global = -0.2036158846, cumulative = -58.36082657
rho max/min : 2 0.6056450947
DICPCG:  Solving for p, Initial residual = 3.693091856e-06, Final residual = 9.695100434e-12, No Iterations 1
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 0.2036158871, global = -0.2036158871, cumulative = -58.56444245
rho max/min : 2 0.6056032962
DILUPBiCG:  Solving for omega, Initial residual = 0.0006652707718, Final residual = 3.054000263e-08, No Iterations 1
DILUPBiCG:  Solving for k, Initial residual = 0.002696792936, Final residual = 2.498247668e-07, No Iterations 1
ExecutionTime = 18.82 s  ClockTime = 19 s

Courant Number mean: 6.451143786e-05 max: 0.09971726235
deltaT = 7.46213813e-09
Time = 2.19409e-06

diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG:  Solving for Ux, Initial residual = 0.003832159235, Final residual = 2.667268357e-07, No Iterations 1
DILUPBiCG:  Solving for Uy, Initial residual = 0.003840701175, Final residual = 1.537092273e-07, No Iterations 1
DILUPBiCG:  Solving for h, Initial residual = 0.002366988328, Final residual = 1.442548449e-07, No Iterations 1
#0  Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1  Foam::sigFpe::sigHandler(int) in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2  in "/lib/x86_64-linux-gnu/libc.so.6"
#3  in "/lib/x86_64-linux-gnu/libm.so.6"
#4  Foam::pow(Foam::Field<double>&, Foam::UList<double> const&, double const&) in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#5  Foam::pow(Foam::UList<double> const&, double const&) in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#6  parserPatch::PatchValueExpressionParser::parse() in "/home/thesis/OpenFOAM/thesis-2.1.0/platforms/linux64GccDPOpt/lib/libswak4FoamParsers.so"
#7  Foam::PatchValueExpressionDriver::parse(std::string const&) in "/home/thesis/OpenFOAM/thesis-2.1.0/platforms/linux64GccDPOpt/lib/libswak4FoamParsers.so"
#8  Foam::groovyBCFvPatchField<double>::updateCoeffs() in "/home/thesis/OpenFOAM/thesis-2.1.0/platforms/linux64GccDPOpt/lib/libgroovyBC.so"
#9  at gaussLaplacianSchemes.C:0
#10  Foam::fv::gaussLaplacianScheme<double, double>::fvmLaplacianUncorrected(Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#11  Foam::fv::gaussLaplacianScheme<double, double>::fvmLaplacian(Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#12  Foam::fv::laplacianScheme<double, double>::fvmLaplacian(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#13  Foam::tmp<Foam::fvMatrix<double> > Foam::fvm::laplacian<double, double>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::word const&) in "/opt/openfoam210/platforms/linux64GccDPOpt/bin/rhoPimpleFoam"
#14  Foam::tmp<Foam::fvMatrix<double> > Foam::fvm::laplacian<double, double>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in "/opt/openfoam210/platforms/linux64GccDPOpt/bin/rhoPimpleFoam"
#15 
 in "/opt/openfoam210/platforms/linux64GccDPOpt/bin/rhoPimpleFoam"
#16  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#17 
 in "/opt/openfoam210/platforms/linux64GccDPOpt/bin/rhoPimpleFoam"
thesis@thesis-X58A-UD7:~/Desktop/WR_Pimple$


gschaider January 9, 2013 18:00

Quote:

Originally Posted by immortality (Post 401017)
please help Bernhard.what does make this error?
Code:

Courant Number mean: 6.451143786e-05 max: 0.09971726235
deltaT = 7.46213813e-09
Time = 2.19409e-06

diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG:  Solving for Ux, Initial residual = 0.003832159235, Final residual = 2.667268357e-07, No Iterations 1
DILUPBiCG:  Solving for Uy, Initial residual = 0.003840701175, Final residual = 1.537092273e-07, No Iterations 1
DILUPBiCG:  Solving for h, Initial residual = 0.002366988328, Final residual = 1.442548449e-07, No Iterations 1
#0  Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1  Foam::sigFpe::sigHandler(int) in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2  in "/lib/x86_64-linux-gnu/libc.so.6"
#3  in "/lib/x86_64-linux-gnu/libm.so.6"
#4  Foam::pow(Foam::Field<double>&, Foam::UList<double> const&, double const&) in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#5  Foam::pow(Foam::UList<double> const&, double const&) in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#6  parserPatch::PatchValueExpressionParser::parse() in "/home/thesis/OpenFOAM/thesis-2.1.0/platforms/linux64GccDPOpt/lib/libswak4FoamParsers.so"
#7  Foam::PatchValueExpressionDriver::parse(std::string const&) in "/home/thesis/OpenFOAM/thesis-2.1.0/platforms/linux64GccDPOpt/lib/libswak4FoamParsers.so"
#8  Foam::groovyBCFvPatchField<double>::updateCoeffs() in "/home/thesis/OpenFOAM/thesis-2.1.0/platforms/linux64GccDPOpt/lib/libgroovyBC.so"
#9  at gaussLaplacianSchemes.C:0
#10  Foam::fv::gaussLaplacianScheme<double, double>::fvmLaplacianUncorrected(Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#11  Foam::fv::gaussLaplacianScheme<double, double>::fvmLaplacian(Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#12  Foam::fv::laplacianScheme<double, double>::fvmLaplacian(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#13  Foam::tmp<Foam::fvMatrix<double> > Foam::fvm::laplacian<double, double>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::word const&) in "/opt/openfoam210/platforms/linux64GccDPOpt/bin/rhoPimpleFoam"
#14  Foam::tmp<Foam::fvMatrix<double> > Foam::fvm::laplacian<double, double>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in "/opt/openfoam210/platforms/linux64GccDPOpt/bin/rhoPimpleFoam"
#15 
 in "/opt/openfoam210/platforms/linux64GccDPOpt/bin/rhoPimpleFoam"
#16  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#17 
 in "/opt/openfoam210/platforms/linux64GccDPOpt/bin/rhoPimpleFoam"
thesis@thesis-X58A-UD7:~/Desktop/WR_Pimple$


You explained your situation in the other thread and I see that you are under pressure, but PLEASE don't post the same question 3 times in the range of an hour.

Well according to your stacktrace the function pow fails. Why that is I don't know. Either it became undefined (for instance with something like pow(-1,0.5) - square root of a negative) or you got an overflow because the number became to big. Have a look at all your expressions that use pow and see what could have gone wrong (mathematically speaking)

rosswin February 28, 2013 05:34

Hi,

I think I am getting a similar pow error.
I have not had much experience with OpenFOAM.
Could you please explain to me the meaning of this error and if possible how I should try to resolve it.

Code:

Time = 114

smoothSolver:  Solving for Ux, Initial residual = 0.503878, Final residual = 0.029159, No Iterations 4
smoothSolver:  Solving for Uz, Initial residual = 0.99092, Final residual = 0.0338131, No Iterations 4
GAMG:  Solving for p, Initial residual = 0.982744, Final residual = 0.0612016, No Iterations 8
time step continuity errors : sum local = 1.37064e+99, global = -1.63939e+85, cumulative = -1.95407e+85
#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::sigFpe::sigHandler(int) at ??:?
#2  Uninterpreted:
#3  Foam::pow3(Foam::Field<double>&, Foam::UList<double> const&) at ??:?
#4  Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::pow3<Foam::fvPatchField, Foam::volMesh>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) at ??:?
#5  Foam::incompressible::RASModels::SpalartAllmaras::fv1(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) const at ??:?
#6  Foam::incompressible::RASModels::SpalartAllmaras::correct() at ??:?
#7 
 at ??:?
#8  __libc_start_main in "/lib/i386-linux-gnu/libc.so.6"
#9 
 at ??:?
Floating point exception (core dumped)

Thank you very much for your time and effort.

Ross

gschaider February 28, 2013 17:52

Quote:

Originally Posted by rosswin (Post 410597)
Hi,

I think I am getting a similar pow error.
I have not had much experience with OpenFOAM.
Could you please explain to me the meaning of this error and if possible how I should try to resolve it.

Code:

Time = 114

smoothSolver:  Solving for Ux, Initial residual = 0.503878, Final residual = 0.029159, No Iterations 4
smoothSolver:  Solving for Uz, Initial residual = 0.99092, Final residual = 0.0338131, No Iterations 4
GAMG:  Solving for p, Initial residual = 0.982744, Final residual = 0.0612016, No Iterations 8
time step continuity errors : sum local = 1.37064e+99, global = -1.63939e+85, cumulative = -1.95407e+85
#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::sigFpe::sigHandler(int) at ??:?
#2  Uninterpreted:
#3  Foam::pow3(Foam::Field<double>&, Foam::UList<double> const&) at ??:?
#4  Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::pow3<Foam::fvPatchField, Foam::volMesh>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) at ??:?
#5  Foam::incompressible::RASModels::SpalartAllmaras::fv1(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) const at ??:?
#6  Foam::incompressible::RASModels::SpalartAllmaras::correct() at ??:?
#7 
 at ??:?
#8  __libc_start_main in "/lib/i386-linux-gnu/libc.so.6"
#9 
 at ??:?
Floating point exception (core dumped)

Thank you very much for your time and effort.

Ross

Have a look at the time-step continuity error: things were going terribly wrong before that (the error in pow3 probably is just a symptom)

There also is an explanation on the meaning of stack-traces on openfoamwiki.net


All times are GMT -4. The time now is 15:35.