|
[Sponsors] |
November 19, 2009, 02:11 |
MRFSimpleFOAM goes divergenced!
|
#1 |
Member
任芸
Join Date: Jun 2009
Posts: 75
Rep Power: 17 |
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 |
|
Thread Tools | Search this Thread |
Display Modes | |
|
|
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 17:09 |
Building error with simpleSRFFoam and MRFSimpleFOAM in 15 | waynezw0618 | OpenFOAM Running, Solving & CFD | 4 | August 27, 2008 10:31 |