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/)
-   -   RNG diverged during the analysis the flow over a multi element airfoil, why? (https://www.cfd-online.com/Forums/openfoam-solving/121777-rng-diverged-during-analysis-flow-over-multi-element-airfoil-why.html)

s.m August 5, 2013 08:39
RNG diverged during the analysis the flow over a multi element airfoil, why?
 
4 Attachment(s)
i am analysing the flow over multi element airfoil, i use free stream boundary condition, RNGkEpsilon turbulence model, for this analysi, but it become diverged, after 10 iteration,
what is wrong with this setting that the "RNGkEpsilon" model decome diverged?

Create mesh for time = 300

Reading field p

Reading field U

Reading/calculating face flux field phi

Selecting incompressible transport model Newtonian
Selecting RAS turbulence model RNGkEpsilon
--> FOAM Warning :
From function Field<Type>::Field(const word& keyword, const dictionary&, const label)
in file /home/alberto/OpenFOAM/OpenFOAM-2.2.x/src/OpenFOAM/lnInclude/Field.C at line 262
Reading "/home/saeide/of/f15/3element-good/noghatekhodam-simpleFoam/stl-riztar/free-stream-very-good-results/RNG-freestream-diverge/300/epsilon.boundaryField.airfoil_solid" from line 46 to line 47
expected keyword 'uniform' or 'nonuniform', assuming deprecated Field format from Foam version 2.0.
--> FOAM Warning :
From function Field<Type>::Field(const word& keyword, const dictionary&, const label)
in file /home/alberto/OpenFOAM/OpenFOAM-2.2.x/src/OpenFOAM/lnInclude/Field.C at line 262
Reading "/home/saeide/of/f15/3element-good/noghatekhodam-simpleFoam/stl-riztar/free-stream-very-good-results/RNG-freestream-diverge/300/epsilon.boundaryField.airfoill_solid" from line 52 to line 53
expected keyword 'uniform' or 'nonuniform', assuming deprecated Field format from Foam version 2.0.
--> FOAM Warning :
From function Field<Type>::Field(const word& keyword, const dictionary&, const label)
in file /home/alberto/OpenFOAM/OpenFOAM-2.2.x/src/OpenFOAM/lnInclude/Field.C at line 262
Reading "/home/saeide/of/f15/3element-good/noghatekhodam-simpleFoam/stl-riztar/free-stream-very-good-results/RNG-freestream-diverge/300/epsilon.boundaryField.airfoils_solid" from line 58 to line 59
expected keyword 'uniform' or 'nonuniform', assuming deprecated Field format from Foam version 2.0.
RNGkEpsilonCoeffs
{
Cmu 0.0845;
C1 1.42;
C2 1.68;
sigmak 0.71942;
sigmaEps 0.71942;
eta0 4.38;
beta 0.012;
}

No finite volume options present


SIMPLE: convergence criteria
field p tolerance 1e-06
field U tolerance 1e-05
field nuTilda tolerance 1e-05
field omega tolerance 1e-05
field k tolerance 1e-05


Starting time loop

Time = 301

smoothSolver: Solving for Ux, Initial residual = 0.00800557, Final residual = 8.03219e-07, No Iterations 11
smoothSolver: Solving for Uy, Initial residual = 0.0123494, Final residual = 5.67371e-07, No Iterations 12
GAMG: Solving for p, Initial residual = 0.0914344, Final residual = 7.85428e-08, No Iterations 21
GAMG: Solving for p, Initial residual = 0.011361, Final residual = 7.82593e-08, No Iterations 14
GAMG: Solving for p, Initial residual = 0.000954229, Final residual = 9.41026e-08, No Iterations 10
time step continuity errors : sum local = 5.69478e-08, global = 6.11691e-15, cumulative = 6.11691e-15
smoothSolver: Solving for epsilon, Initial residual = 0.353798, Final residual = 0.00918427, No Iterations 100
smoothSolver: Solving for k, Initial residual = 1, Final residual = 7.01656e-07, No Iterations 11
ExecutionTime = 2.62 s ClockTime = 2 s

forceCoeffs output:
Cm = 1.36917
Cd = 0.217473
Cl = 2.83816
Cl(f) = 2.78825
Cl(r) = 0.0499078

forceCoeffs output:
Cm = 0.42901
Cd = 0.301073
Cl = 0.393971
Cl(f) = 0.625995
Cl(r) = -0.232024

forceCoeffs output:
Cm = -0.00911897
Cd = -0.0815906
Cl = 0.212466
Cl(f) = 0.0971139
Cl(r) = 0.115352

forces output:
forces(pressure,viscous)((-26.7641 212.265 -5.47395e-19),(0.303963 0.121734 1.8796e-21))
moment(pressure,viscous)((-42.453 -5.35282 54.6227),(-0.0243469 0.0607927 0.000624143))



Time = 305

smoothSolver: Solving for Ux, Initial residual = 0.000527167, Final residual = 5.34816e-07, No Iterations 10
smoothSolver: Solving for Uy, Initial residual = 0.00129551, Final residual = 8.2923e-07, No Iterations 11
GAMG: Solving for p, Initial residual = 0.265218, Final residual = 8.90272e-08, No Iterations 59
GAMG: Solving for p, Initial residual = 0.0116877, Final residual = 9.93816e-08, No Iterations 37
GAMG: Solving for p, Initial residual = 0.000903402, Final residual = 9.85498e-08, No Iterations 24
time step continuity errors : sum local = 9.27064e-08, global = 2.28966e-16, cumulative = 2.23089e-14
smoothSolver: Solving for epsilon, Initial residual = 0.231188, Final residual = 3.96873e-06, No Iterations 100
bounding epsilon, min: -0.250843 max: 1.33142e+07 average: 81719.6
smoothSolver: Solving for k, Initial residual = 0.00195982, Final residual = 5.89048e-07, No Iterations 6
ExecutionTime = 15.05 s ClockTime = 15 s

forceCoeffs output:
Cm = 1.35544
Cd = 0.198891
Cl = 2.8304
Cl(f) = 2.77063
Cl(r) = 0.0597623

forceCoeffs output:
Cm = 0.42617
Cd = 0.297069
Cl = 0.392494
Cl(f) = 0.622417
Cl(r) = -0.229923

forceCoeffs output:
Cm = -0.0101199
Cd = -0.0873238
Cl = 0.210361
Cl(f) = 0.0950609
Cl(r) = 0.115301

forces output:
forces(pressure,viscous)((-27.8153 211.754 -5.68802e-19),(0.562578 0.130968 6.85759e-20))
moment(pressure,viscous)((-42.3507 -5.56306 54.0607),(-0.0261936 0.112516 -0.00682002))

T
Time = 309

smoothSolver: Solving for Ux, Initial residual = 0.49735, Final residual = 6.23036e-07, No Iterations 17
smoothSolver: Solving for Uy, Initial residual = 0.459138, Final residual = 5.22546e-07, No Iterations 16
GAMG: Solving for p, Initial residual = 0.00256564, Final residual = 8.10012e-08, No Iterations 43
GAMG: Solving for p, Initial residual = 2.26878e-08, Final residual = 2.26878e-08, No Iterations 0
GAMG: Solving for p, Initial residual = 2.26878e-08, Final residual = 2.26878e-08, No Iterations 0
time step continuity errors : sum local = 1.82335e+09, global = -2.59111e-07, cumulative = -2.62237e-07
smoothSolver: Solving for epsilon, Initial residual = 0.973025, Final residual = 8.42558e-08, No Iterations 1
bounding epsilon, min: -1.45831e+15 max: 3.20811e+32 average: 7.54731e+27
smoothSolver: Solving for k, Initial residual = 0.309414, Final residual = 4.97484e-07, No Iterations 10
ExecutionTime = 29.25 s ClockTime = 29 s

forceCoeffs output:
Cm = 2.46333e+17
Cd = 3.20069e+16
Cl = 4.98601e+17
Cl(f) = 4.95633e+17
Cl(r) = 2.96779e+15

forceCoeffs output:
Cm = 2.52658e+12
Cd = 8.6173e+12
Cl = 1.88077e+12
Cl(f) = 3.46697e+12
Cl(r) = -1.58619e+12

forceCoeffs output:
Cm = -5.43653e+11
Cd = -7.79626e+12
Cl = 5.91847e+12
Cl(f) = 2.41558e+12
Cl(r) = 3.50289e+12

forces output:
forces(pressure,viscous)((-2.82209e+18 4.78317e+19 -0.540524),(-4.26067e+14 3.43531e+13 -0.000406198))
moment(pressure,viscous)((-9.56634e+18 -5.64419e+17 1.41743e+19),(-6.87062e+12 -8.52133e+13 -1.75202e+12))

Time = 310
Duration of pickling 0.0778450965881 too long. Extending frequency from 1.0 to 3.89225482941

smoothSolver: Solving for Ux, Initial residual = 0.139675, Final residual = 5.12311e-07, No Iterations 15
smoothSolver: Solving for Uy, Initial residual = 0.364373, Final residual = 7.20488e-07, No Iterations 16
GAMG: Solving for p, Initial residual = 0.000169534, Final residual = 0.000536, No Iterations 100
GAMG: Solving for p, Initial residual = 4.52632e-11, Final residual = 4.52632e-11, No Iterations 0
GAMG: Solving for p, Initial residual = 4.52632e-11, Final residual = 4.52632e-11, No Iterations 0
time step continuity errors : sum local = 4.56399e+14, global = -0.678477, cumulative = -0.678477
smoothSolver: Solving for epsilon, Initial residual = 0.32509, Final residual = 2.91165e-07, No Iterations 1
bounding epsilon, min: -3.22545e+29 max: 2.97866e+33 average: 1.54415e+29
smoothSolver: Solving for k, Initial residual = 0.696664, Final residual = 3.17524e-07, No Iterations 11
bounding k, min: -3.26008e+18 max: 7.62899e+23 average: 1.04585e+20
ExecutionTime = 31.73 s ClockTime = 32 s

forceCoeffs output:
Cm = 4.92094e+32
Cd = 3.50449e+31
Cl = 9.11494e+32
Cl(f) = 9.47841e+32
Cl(r) = -3.6347e+31

forceCoeffs output:
Cm = 7.12248e+14
Cd = 6.87955e+14
Cl = 6.19634e+14
Cl(f) = 1.02206e+15
Cl(r) = -4.02431e+14

forceCoeffs output:
Cm = -9.69654e+16
Cd = -5.61894e+18
Cl = 6.27121e+18
Cl(f) = 3.03864e+18
Cl(r) = 3.23257e+18

forces output:
forces(pressure,viscous)((-7.39342e+33 8.71665e+34 5.56692e+14),(3.01979e+18 1.10144e+17 0.726169))
moment(pressure,viscous)((-1.74333e+34 -1.47868e+33 2.83159e+34),(-2.20288e+16 6.03957e+17 1.2933e+17))

Time = 311

smoothSolver: Solving for Ux, Initial residual = 0.419499, Final residual = 5.49432e-07, No Iterations 17
smoothSolver: Solving for Uy, Initial residual = 0.374946, Final residual = 8.83925e-07, No Iterations 16
GAMG: Solving for p, Initial residual = 0.997043, Final residual = 0.000130123, No Iterations 100
GAMG: Solving for p, Initial residual = 0.981997, Final residual = 1.0727e-06, No Iterations 100
GAMG: Solving for p, Initial residual = 0.0648716, Final residual = 2.2693e-07, No Iterations 100
time step continuity errors : sum local = 3.6604e+19, global = 4.86138e+06, cumulative = 4.86138e+06
smoothSolver: Solving for epsilon, Initial residual = 1, Final residual = 0.732293, No Iterations 100
smoothSolver: Solving for k, Initial residual = 0.430025, Final residual = 1.12638e-09, No Iterations 1
ExecutionTime = 38.25 s ClockTime = 38 s

forceCoeffs output:
Cm = 1.3753e+33
Cd = 1.61358e+32
Cl = 2.71282e+33
Cl(f) = 2.73171e+33
Cl(r) = -1.88959e+31

forceCoeffs output:
Cm = -8.41831e+24
Cd = -2.16018e+25
Cl = -1.13321e+25
Cl(f) = -1.40843e+25
Cl(r) = 2.75228e+24

forceCoeffs output:
Cm = -4.71304e+22
Cd = -1.00122e+24
Cl = 1.70514e+24
Cl(f) = 8.05441e+23
Cl(r) = 8.99702e+23

forces output:
forces(pressure,viscous)((-1.69807e+34 2.63097e+35 2.56593e+14),(4.06702e+32 -2.99769e+33 -6.32977e+12))
moment(pressure,viscous)((-5.26194e+34 -3.39613e+33 8.01073e+34),(5.99538e+32 8.13403e+31 -9.69927e+32))

Time = 312
Duration of pickling 0.0470168590546 too long. Extending frequency from 1.0 to 2.35084295273

smoothSolver: Solving for Ux, Initial residual = 3.84897e-05, Final residual = 4.89998e-07, No Iterations 4
smoothSolver: Solving for Uy, Initial residual = 0.000175971, Final residual = 5.12298e-07, No Iterations 6
GAMG: Solving for p, Initial residual = 0.923463, Final residual = 2.76536e+17, No Iterations 100
GAMG: Solving for p, Initial residual = 1, Final residual = 8.80639e-08, No Iterations 4
GAMG: Solving for p, Initial residual = 0.573316, Final residual = 4.85976e-08, No Iterations 3
time step continuity errors : sum local = 9.33282e+49, global = 2.35406e+34, cumulative = 2.35406e+34
smoothSolver: Solving for epsilon, Initial residual = 1, Final residual = 1.5617e+11, No Iterations 100
bounding epsilon, min: -1.0805e+81 max: 2.15564e+103 average: 1.90597e+101
smoothSolver: Solving for k, Initial residual = 2.11034e-17, Final residual = 2.11034e-17, No Iterations 0
ExecutionTime = 43.63 s ClockTime = 44 s

forceCoeffs output:
Cm = -2.28422e+59
Cd = -1.5873e+59
Cl = -1.48431e+59
Cl(f) = -3.02638e+59
Cl(r) = 1.54207e+59

forceCoeffs output:
Cm = -3.21871e+59
Cd = -2.03979e+59
Cl = -2.63668e+59
Cl(f) = -4.53705e+59
Cl(r) = 1.90037e+59

forceCoeffs output:
Cm = 5.88367e+57
Cd = 3.68239e+58
Cl = -7.32437e+57
Cl(f) = 2.22149e+57
Cl(r) = -9.54586e+57

forces output:
forces(pressure,viscous)((-6.40749e+59 1.17643e+61 8.63515e+40),(1.06485e+50 1.57963e+49 1.30201e+34))
moment(pressure,viscous)((-2.35286e+60 -1.2815e+59 5.03863e+60),(-3.15925e+48 2.12969e+49 8.52946e+48))

Time = 313
Duration of pickling 0.117083072662 too long. Extending frequency from 1.0 to 5.85415363312

smoothSolver: Solving for Ux, Initial residual = 0.658273, Final residual = 5.32991e-07, No Iterations 20
smoothSolver: Solving for Uy, Initial residual = 0.122935, Final residual = 7.89541e-07, No Iterations 17
#0 Foam::error::printStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 in "/lib64/libc.so.6"
#3 Foam::GaussSeidelSmoother::smooth(Foam::word const&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::Field<double> const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, unsigned char, int) at ??:?
#4 Foam::GaussSeidelSmoother::smooth(Foam::Field<doub le>&, Foam::Field<double> const&, unsigned char, int) const at ??:?
#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 at ??:?
#6 Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
#7 Foam::fvMatrix<double>::solveSegregated(Foam::dict ionary const&) at ??:?
#8 Foam::fvMatrix<double>::solve(Foam::dictionary const&) at ??:?
#9
at ??:?
#10 __libc_start_main in "/lib64/libc.so.6"
#11
at /home/abuild/rpmbuild/BUILD/glibc-2.17/csu/../sysdeps/x86_64/start.S:126


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