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/)
-   -   Foam:: Error:: PrintStack (https://www.cfd-online.com/Forums/openfoam-solving/108300-foam-error-printstack.html)

jr33 October 18, 2012 18:03

Foam:: Error:: PrintStack
 
I am running an example of flow in pipe using openfoam (Ubuntu), I get following error, although it works fine till time point 16, How to fix the error.
Files are here:
http://dl.dropbox.com/u/68983831/files.zip


// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0

Reading transportProperties

Reading field p

Reading field U

Reading/calculating face flux field phi


Starting time loop

Time = 15

Courant Number mean: 2.30651e+31 max: 4.76733e+34
DILUPBiCG: Solving for Ux, Initial residual = 0.998836, Final residual = 0.0837974, No Iterations 14
DILUPBiCG: Solving for Uy, Initial residual = 0.999846, Final residual = 0.079946, No Iterations 5
DILUPBiCG: Solving for Uz, Initial residual = 0.999593, Final residual = 0.0819008, No Iterations 13
DICPCG: Solving for p, Initial residual = 1, Final residual = 0.000810451, No Iterations 2
time step continuity errors : sum local = 8.6096e+60, global = -1.31187e+44, cumulative = -1.31187e+44
DICPCG: Solving for p, Initial residual = 0.997834, Final residual = 0.00011985, No Iterations 1
time step continuity errors : sum local = 3.07997e+62, global = 9.1217e+45, cumulative = 8.99052e+45
ExecutionTime = 1.84 s ClockTime = 3 s

Time = 16

Courant Number mean: 3.06096e+62 max: 5.97395e+65
DILUPBiCG: Solving for Ux, Initial residual = 0.585237, Final residual = 0.0331803, No Iterations 23
DILUPBiCG: Solving for Uy, Initial residual = 0.58276, Final residual = 0.0407138, No Iterations 12
DILUPBiCG: Solving for Uz, Initial residual = 0.908282, Final residual = 0.0424477, No Iterations 1
#0 Foam::error:rintStack(Foam::Ostream&) in "/opt/openfoam211/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
#1 Foam::sigFpe::sigHandler(int) in "/opt/openfoam211/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
#2 Uninterpreted:
#3 Foam:ICPreconditioner::calcReciprocalD(Foam::Field <double>&, Foam::lduMatrix const&) in "/opt/openfoam211/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
#4 Foam:ICPreconditioner:ICPreconditioner(Foam::lduMa trix::solver const&, Foam::dictionary const&) in "/opt/openfoam211/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
#5 Foam::lduMatrix:reconditioner::addsymMatrixConstru ctorToTable<Foam:ICPreconditioner>::New(Foam::lduM atrix::solver const&, Foam::dictionary const&) in "/opt/openfoam211/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
#6 Foam::lduMatrix:reconditioner::New(Foam::lduMatrix ::solver const&, Foam::dictionary const&) in "/opt/openfoam211/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
#7 Foam::PCG::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in "/opt/openfoam211/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
#8 Foam::fvMatrix<double>::solve(Foam::dictionary const&) in "/opt/openfoam211/platforms/linuxGccDPOpt/lib/libfiniteVolume.so"
#9
in "/opt/openfoam211/platforms/linuxGccDPOpt/bin/icoFoam"
#10 __libc_start_main in "/lib/i386-linux-gnu/libc.so.6"
#11
in "/opt/openfoam211/platforms/linuxGccDPOpt/bin/icoFoam"
Floating point exception (core dumped)

//************************************************** ************************************************** *****//

gschaider October 18, 2012 18:19

Quote:

Originally Posted by jr33 (Post 387407)
I am running an example of flow in pipe using openfoam (Ubuntu), I get following error, although it works fine till time point 16,

No it doesn't. Have a look at time 15: Courant numbers of the order 1e34 are not fine in this universe

Problems start probably even further up (there is probably a problem with your case)

Quote:

Originally Posted by jr33 (Post 387407)
How to fix the error.
Files are here

Best way to find the error is write out every time-step (with 16 timesteps to failure this should be possible), have a look at them and find out what is going wrong

jr33 October 18, 2012 18:23

Quote:

Originally Posted by gschaider (Post 387408)
No it doesn't. Have a look at time 15: Courant numbers of the order 1e34 are not fine in this universe

Problems start probably even further up (there is probably a problem with your case)



Best way to find the error is write out every time-step (with 16 timesteps to failure this should be possible), have a look at them and find out what is going wrong


I saw at time point 1:
Courant Number mean: 5.80776 max: 1368.61
Thanks

jr33 October 18, 2012 18:29

Above error came when I fixed keyword laplacian
 
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 10

Reading transportProperties

Reading field p

Reading field U

Reading/calculating face flux field phi


Starting time loop

Time = 11

Courant Number mean: 1.54422e+12 max: 7.57587e+14


--> FOAM FATAL IO ERROR:
keyword laplacian(nu,U) is undefined in dictionary "/home/jay/OpenFOAM/jay-2.1.1/run/tutorials/incompressible/simpleFoam/elbow/system/fvSchemes::laplacianSchemes"

file: /home/jay/OpenFOAM/jay-2.1.1/run/tutorials/incompressible/simpleFoam/elbow/system/fvSchemes::laplacianSchemes from line 44 to line 50.

From function dictionary::lookupEntry(const word&, bool, bool) const
in file db/dictionary/dictionary.C at line 400.

FOAM exiting


All times are GMT -4. The time now is 06:40.