CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (http://www.cfd-online.com/Forums/openfoam-solving/)
-   -   MRFSimpleFOAM goes divergenced! (http://www.cfd-online.com/Forums/openfoam-solving/70276-mrfsimplefoam-goes-divergenced.html)

renyun0511 November 19, 2009 03:11

MRFSimpleFOAM goes divergenced!
 
hi all
i am calculating an impeller of pump,which was meshed by Gambit,using MRFSimpleFoam (OF1.5),and after several steps,it divergenced,the log and error are as follow:
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 1.5 |
| \\ / A nd | Web: http://www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Exec : MRFSimpleFoam
Date : Nov 19 2009
Time : 10:06:21
Host : linux-pw3p
PID : 4721
Case : /home/ry/OpenFOAM/ry-1.5/tutorials/MRFSimpleFoam/091119
nProcs : 1

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

Create mesh for time = 0

Reading field p

Reading field U

Reading/calculating face flux field phi

Selecting incompressible transport model Newtonian
Selecting RAS turbulence model kEpsilon
kEpsilonCoeffs
{
Cmu 0.09;
C1 1.44;
C2 1.92;
alphaEps 0.76923;
}


Starting time loop

Time = 1

smoothSolver: Solving for Ux, Initial residual = 1, Final residual = 0.0200145, No Iterations 2
smoothSolver: Solving for Uy, Initial residual = 1, Final residual = 0.0201876, No Iterations 2
smoothSolver: Solving for Uz, Initial residual = 1, Final residual = 0.0220911, No Iterations 2
GAMG: Solving for p, Initial residual = 1, Final residual = 0.0458254, No Iterations 5
time step continuity errors : sum local = 32.1584, global = -15.972, cumulative = -15.972
smoothSolver: Solving for epsilon, Initial residual = 0.377664, Final residual = 0.0176228, No Iterations 2
bounding epsilon, min: -20.6763 max: 328.148 average: 8.48132
smoothSolver: Solving for k, Initial residual = 1, Final residual = 0.0657878, No Iterations 2
bounding k, min: -3.18187 max: 28.4237 average: 0.716773
ExecutionTime = 57.54 s ClockTime = 64 s

Time = 2

smoothSolver: Solving for Ux, Initial residual = 0.184681, Final residual = 0.00751509, No Iterations 2
smoothSolver: Solving for Uy, Initial residual = 0.184377, Final residual = 0.00751749, No Iterations 2
smoothSolver: Solving for Uz, Initial residual = 0.319056, Final residual = 0.00961937, No Iterations 2
GAMG: Solving for p, Initial residual = 0.51407, Final residual = 0.0136539, No Iterations 5
time step continuity errors : sum local = 30.467, global = -10.0327, cumulative = -26.0047
smoothSolver: Solving for epsilon, Initial residual = 0.46798, Final residual = 0.0345097, No Iterations 2
bounding epsilon, min: -107.641 max: 34135.9 average: 73.8284
smoothSolver: Solving for k, Initial residual = 0.693379, Final residual = 0.0663071, No Iterations 2
bounding k, min: -6.09496 max: 201.319 average: 3.88779
ExecutionTime = 94.54 s ClockTime = 102 s

Time = 3

smoothSolver: Solving for Ux, Initial residual = 0.179685, Final residual = 0.0108861, No Iterations 2
smoothSolver: Solving for Uy, Initial residual = 0.180145, Final residual = 0.0111631, No Iterations 2
smoothSolver: Solving for Uz, Initial residual = 0.333866, Final residual = 0.0100125, No Iterations 2
GAMG: Solving for p, Initial residual = 0.18048, Final residual = 0.0084598, No Iterations 4
time step continuity errors : sum local = 68.6218, global = 16.4865, cumulative = -9.51821
smoothSolver: Solving for epsilon, Initial residual = 0.153731, Final residual = 0.00253116, No Iterations 2
bounding epsilon, min: -12055.4 max: 665628 average: 775.906
smoothSolver: Solving for k, Initial residual = 0.404658, Final residual = 0.0358232, No Iterations 2
bounding k, min: -7.47208 max: 385.105 average: 8.65698
ExecutionTime = 129.72 s ClockTime = 138 s

Time = 4

smoothSolver: Solving for Ux, Initial residual = 0.158165, Final residual = 0.0139721, No Iterations 2
smoothSolver: Solving for Uy, Initial residual = 0.164707, Final residual = 0.0138668, No Iterations 2
smoothSolver: Solving for Uz, Initial residual = 0.219469, Final residual = 0.00927912, No Iterations 2
GAMG: Solving for p, Initial residual = 0.205141, Final residual = 0.00886494, No Iterations 2
time step continuity errors : sum local = 129.77, global = 2.40719, cumulative = -7.11102
smoothSolver: Solving for epsilon, Initial residual = 0.9131, Final residual = 0.055472, No Iterations 2
bounding epsilon, min: -4.83352e+07 max: 1.96105e+09 average: 19641.1
smoothSolver: Solving for k, Initial residual = 0.874109, Final residual = 0.0143231, No Iterations 2
bounding k, min: -68307.1 max: 358792 average: 27.0937
ExecutionTime = 161.08 s ClockTime = 170 s

Time = 5

smoothSolver: Solving for Ux, Initial residual = 0.198946, Final residual = 0.00370692, No Iterations 2
smoothSolver: Solving for Uy, Initial residual = 0.495621, Final residual = 0.000454703, No Iterations 2
smoothSolver: Solving for Uz, Initial residual = 0.0174051, Final residual = 0.000303364, No Iterations 2
GAMG: Solving for p, Initial residual = 0.192162, Final residual = 0.00442174, No Iterations 5
time step continuity errors : sum local = 50.1274, global = 1.48055, cumulative = -5.63047
smoothSolver: Solving for epsilon, Initial residual = 4.12832e-05, Final residual = 3.26667e-07, No Iterations 2
bounding epsilon, min: -2.7054e+09 max: 4.57078e+10 average: 813456
smoothSolver: Solving for k, Initial residual = 0.00682609, Final residual = 0.000230793, No Iterations 2
bounding k, min: -49889.8 max: 6.10063e+06 average: 136.631
ExecutionTime = 198.18 s ClockTime = 208 s

Time = 6

smoothSolver: Solving for Ux, Initial residual = 0.0855071, Final residual = 0.00576353, No Iterations 2
smoothSolver: Solving for Uy, Initial residual = 0.0395296, Final residual = 0.000789192, No Iterations 2
smoothSolver: Solving for Uz, Initial residual = 0.224219, Final residual = 0.00988882, No Iterations 2
GAMG: Solving for p, Initial residual = 0.985461, Final residual = 0.0266543, No Iterations 2
time step continuity errors : sum local = 6.5826e+09, global = 5.5398e+08, cumulative = 5.5398e+08
smoothSolver: Solving for epsilon, Initial residual = 0.995357, Final residual = 0.00405112, No Iterations 4
bounding epsilon, min: -2.16538e+14 max: 2.24754e+16 average: 1.02789e+11
smoothSolver: Solving for k, Initial residual = 0.999998, Final residual = 0.0255665, No Iterations 2
bounding k, min: -2.50162e+13 max: 2.77207e+15 average: 9.55502e+09
ExecutionTime = 230.17 s ClockTime = 241 s

Time = 7

smoothSolver: Solving for Ux, Initial residual = 0.495146, Final residual = 0.0162337, No Iterations 2
smoothSolver: Solving for Uy, Initial residual = 0.560682, Final residual = 0.0162687, No Iterations 2
smoothSolver: Solving for Uz, Initial residual = 0.30485, Final residual = 0.00376747, No Iterations 2
GAMG: Solving for p, Initial residual = 0.996045, Final residual = 0.0263017, No Iterations 2
time step continuity errors : sum local = 3.31339e+09, global = 6.40301e+08, cumulative = 1.19428e+09
smoothSolver: Solving for epsilon, Initial residual = 0.160724, Final residual = 3.38554e-12, No Iterations 2
bounding epsilon, min: -1.57739e+21 max: 4.989e+25 average: 7.92797e+20
smoothSolver: Solving for k, Initial residual = 0.667466, Final residual = 0.0239193, No Iterations 2
bounding k, min: -3.40807e+17 max: 4.20692e+24 average: 3.59773e+19
ExecutionTime = 261.62 s ClockTime = 273 s

Time = 8

smoothSolver: Solving for Ux, Initial residual = 0.830768, Final residual = 0.04517, No Iterations 2
smoothSolver: Solving for Uy, Initial residual = 0.69165, Final residual = 0.00968096, No Iterations 2
smoothSolver: Solving for Uz, Initial residual = 0.513691, Final residual = 0.00846947, No Iterations 2
GAMG: Solving for p, Initial residual = 0.000344709, Final residual = 9.67256e-06, No Iterations 2
time step continuity errors : sum local = 9.57166e+22, global = -8.20482e+21, cumulative = -8.20482e+21
smoothSolver: Solving for epsilon, Initial residual = 0.999659, Final residual = 0.0158262, No Iterations 2
bounding epsilon, min: -1.50016e+40 max: 6.40414e+43 average: 2.56412e+38
smoothSolver: Solving for k, Initial residual = 0.00843413, Final residual = 0.000557551, No Iterations 2
bounding k, min: -6.14028e+34 max: 4.24e+35 average: 1.87793e+30
ExecutionTime = 292.88 s ClockTime = 305 s

Time = 9

smoothSolver: Solving for Ux, Initial residual = 0.000252193, Final residual = 5.32787e-07, No Iterations 2
smoothSolver: Solving for Uy, Initial residual = 4.39356e-05, Final residual = 3.30549e-07, No Iterations 2
smoothSolver: Solving for Uz, Initial residual = 0.000349894, Final residual = 6.32743e-07, No Iterations 2

ry@linux-pw3p:~/OpenFOAM/ry-1.5/tutorials/MRFSimpleFoam/091119> #0 Foam::error::printStack(Foam::Ostream&) in "/home/ry/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libOpenFOAM.so"

#1 Foam::sigFpe::sigFpeHandler(int) in "/home/ry/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libOpenFOAM.so"

#2 Uninterpreted: [0xffffe420]

#3 Foam::GAMGSolver::scalingFactor(Foam::Field<double >&, Foam::Field<double> const&, Foam::Field<double> const&, Foam::Field<double> const&) const in "/home/ry/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libOpenFOAM.so"

#4 Foam::GAMGSolver::scalingFactor(Foam::Field<double >&, Foam::lduMatrix const&, Foam::Field<double>&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, Foam::Field<double> const&, unsigned char) const in "/home/ry/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libOpenFOAM.so"

#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::PtrList<Foam::Field<double> >&, Foam::PtrList<Foam::Field<double> >&, unsigned char) const in "/home/ry/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libOpenFOAM.so"

#6 Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in "/home/ry/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libOpenFOAM.so"

#7 Foam::fvMatrix<double>::solve(Foam::Istream&) in "/home/ry/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libfiniteVolume.so"

#8 main in "/home/ry/OpenFOAM/ry-1.5/applications/bin/linuxGccDPOpt/MRFSimpleFoam"

#9 __libc_start_main in "/lib/libc.so.6"

#10 Foam::regIOobject::readIfModified() in "/home/ry/OpenFOAM/ry-1.5/applications/bin/linuxGccDPOpt/MRFSimpleFoam"


i don't know why the k and epsilon are changed so large! can someone help me?thanks!
yours jennifer


All times are GMT -4. The time now is 07:23.