|
[Sponsors] |
November 20, 2009, 10:03 |
pisoFoam Error
|
#1 |
Member
Join Date: Nov 2009
Location: Munich
Posts: 43
Rep Power: 16 |
Hej all,
I recieve this error when I execute pisoFoam: __________________________________________________ ____________________ Calculating averages Time = 0.00142 Courant Number mean: 3.72507e+52 max: 1.36418e+54 DILUPBiCG: Solving for Ux, Initial residual = 0.555198, Final residual = 1.81207e-06, No Iterations 35 DILUPBiCG: Solving for Uy, Initial residual = 0.882337, Final residual = 2.82071e-06, No Iterations 36 DILUPBiCG: Solving for Uz, Initial residual = 0.959453, Final residual = 2.3849e-06, No Iterations 29 DICPCG: Solving for p, Initial residual = 0.685706, Final residual = 0.0341773, No Iterations 70 time step continuity errors : sum local = 2.65701e+50, global = 1.48284e+49, cumulative = 1.38529e+49 DICPCG: Solving for p, Initial residual = 0.684507, Final residual = 9.50829e-07, No Iterations 171 time step continuity errors : sum local = 2.06627e+46, global = -1.52172e+44, cumulative = 1.38528e+49 DILUPBiCG: Solving for k, Initial residual = 0.88271, Final residual = 7.51988e-06, No Iterations 6 bounding k, min: -3.72039e+117 max: 7.6031e+117 average: 2.9493e+116 ExecutionTime = 164.59 s ClockTime = 165 s Calculating averages Time = 0.00143 Courant Number mean: 8.88307e+51 max: 2.84151e+53 DILUPBiCG: Solving for Ux, Initial residual = 0.856534, Final residual = 7.94225e-06, No Iterations 39 DILUPBiCG: Solving for Uy, Initial residual = 0.826641, Final residual = 9.82446e-06, No Iterations 37 DILUPBiCG: Solving for Uz, Initial residual = 0.6356, Final residual = 6.89001e-06, No Iterations 35 DICPCG: Solving for p, Initial residual = 0.702043, Final residual = 0.0343882, No Iterations 78 time step continuity errors : sum local = 6.28092e+51, global = -2.13804e+49, cumulative = -7.5276e+48 DICPCG: Solving for p, Initial residual = 0.778709, Final residual = 8.64955e-07, No Iterations 171 time step continuity errors : sum local = 5.07704e+47, global = -7.17519e+44, cumulative = -7.52832e+48 #0 Foam::error:rintStack(Foam::Ostream&) in "/opt/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libOpenFOAM.so" #1 Foam::sigFpe::sigFpeHandler(int) in "/opt/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libOpenFOAM.so" #2 ?? in "/lib/libc.so.6" #3 Foam::PBiCG::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in "/opt/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libOpenFOAM.so" #4 Foam::fvMatrix<double>::solve(Foam::dictionary const&) in "/opt/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libfiniteVolume.so" #5 Foam::incompressible::LESModels:neEqEddy::correc t(Foam::tmp<Foam::GeometricField<Foam::Tensor<doub le>, Foam::fvPatchField, Foam::volMesh> > const&) in "/opt/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libincompressibleLESModels.so" #6 Foam::incompressible::LESModel::correct() in "/opt/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libincompressibleLESModels.so" #7 main in "/opt/OpenFOAM/OpenFOAM-1.6/applications/bin/linux64GccDPOpt/pisoFoam" #8 __libc_start_main in "/lib/libc.so.6" #9 _start at /usr/src/packages/BUILD/glibc-2.9/csu/../sysdeps/x86_64/elf/start.S:116 Gleitkomma-Ausnahme __________________________________________________ ____________________ It's running for a while and then it crashes. Could it be, that my Courant Number is much to hight? Or isn't there enough memory? Have I nice weekend, Erik PS: I'll only be back on monday :-) |
|
November 20, 2009, 14:23 |
|
#2 |
Senior Member
Anton Kidess
Join Date: May 2009
Location: Germany
Posts: 1,377
Rep Power: 29 |
Yes, your simulations are diverging and the Courant number gets way out of hand. As usual, check your boundary and initial conditions, and then your numerical parameters. Time step looks like it's pretty small already, so I'm assuming that's not your problem.
|
|
November 23, 2009, 07:49 |
|
#3 |
Member
Join Date: Nov 2009
Location: Munich
Posts: 43
Rep Power: 16 |
Hi akidess,
thx for your help. The courant number depends on the speed, the timestep and the cell size, right? But how can I improve the courant number? I have already tried even smaler time steps and also a smaler inital speed U, but with the same result. What else could I do to make the calculation converging? |
|
November 23, 2009, 18:13 |
|
#4 | |
Senior Member
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,912
Rep Power: 36 |
Quote:
About your case, it would be useful to know some more details on the case and the numerical setup (fvSchemes, essentially). Best,
__________________
Alberto Passalacqua GeekoCFD - A free distribution based on openSUSE 64 bit with CFD tools, including OpenFOAM. Available as in both physical and virtual formats (current status: http://albertopassalacqua.com/?p=1541) OpenQBMM - An open-source implementation of quadrature-based moment methods. To obtain more accurate answers, please specify the version of OpenFOAM you are using. Last edited by alberto; November 23, 2009 at 18:14. Reason: Corrected link |
||
November 23, 2009, 19:53 |
|
#5 | |
Senior Member
Anton Kidess
Join Date: May 2009
Location: Germany
Posts: 1,377
Rep Power: 29 |
Quote:
|
||
November 24, 2009, 03:33 |
|
#6 | ||||
Member
Join Date: Nov 2009
Location: Munich
Posts: 43
Rep Power: 16 |
Thank you very much, I really appreciate your help!
My case is about a pipe flow with a tube in the middle. This is my fvSchemes file: Quote:
There are the initial conditions for U: Quote:
Quote:
Quote:
Sorry for the long post. Last edited by sErik; November 24, 2009 at 05:07. |
|||||
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Compile problem | ivanyao | OpenFOAM Running, Solving & CFD | 1 | October 12, 2012 09:31 |
Problem with compile the setParabolicInlet | ivanyao | OpenFOAM Running, Solving & CFD | 6 | September 5, 2008 20:50 |
Compiling problems with hello worldC | fw407 | OpenFOAM Installation | 21 | January 6, 2008 17:38 |
DecomposePar links against liblamso0 with OpenMPI | jens_klostermann | OpenFOAM Bugs | 11 | June 28, 2007 17:51 |
user subroutine error | CFDUSER | CFX | 2 | December 9, 2006 06:31 |