CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

MRFSimpleFOAM goes divergenced!

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   November 19, 2009, 03:11
Smile MRFSimpleFOAM goes divergenced!
  #1
Member
 
任芸
Join Date: Jun 2009
Posts: 75
Rep Power: 8
renyun0511 is on a distinguished road
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:rintStack(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
renyun0511 is offline   Reply With Quote

Reply

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
MRFSimpleFoam Tutorial bastil OpenFOAM Running, Solving & CFD 48 August 1, 2012 10:00
GGI implementation in MRFSimpleFoam amgode OpenFOAM Running, Solving & CFD 8 August 5, 2011 06:03
MRFSimpleFoam with cyclicGgi in parallel hannes OpenFOAM Bugs 7 April 29, 2010 12:50
Difference between OF1.5x and OF1.5-dev for the MRFSimpleFoam tutorial mai OpenFOAM 0 November 15, 2009 18:09
Building error with simpleSRFFoam and MRFSimpleFOAM in 15 waynezw0618 OpenFOAM Running, Solving & CFD 4 August 27, 2008 10:31


All times are GMT -4. The time now is 21:12.