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

RNG diverged during the analysis the flow over a multi element airfoil, why?

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

Reply
 
LinkBack Thread Tools Display Modes
Old   August 5, 2013, 08:39
Default RNG diverged during the analysis the flow over a multi element airfoil, why?
  #1
s.m
Senior Member
 
saeideh mohamadi
Join Date: Aug 2012
Posts: 229
Rep Power: 6
s.m is on a distinguished road
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:rintStack(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
Attached Images
File Type: jpg residual194.jpg (24.0 KB, 2 views)
Attached Files
File Type: txt fvSchemes.txt (3.5 KB, 2 views)
File Type: txt fvSolution.txt (4.6 KB, 2 views)
File Type: gz 0.tar.gz (1.3 KB, 2 views)
s.m 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
Supersonic flow over an airfoil. FLUENT Lift don't match up with Shock Expansion Tech enriccasas FLUENT 1 August 21, 2013 06:43
Re: Meshing in 2d analysis of incompressible flow over airfoil yashganatra FLUENT 4 June 19, 2013 07:34
Problem with airfoil shape optimization robyTKD SU2 Shape Design 6 June 13, 2013 19:13
Problem with restart solution in shape_optimization.py robyTKD SU2 Shape Design 21 May 29, 2013 09:26
problem in modeling flow over a 3d airfoil guess Main CFD Forum 8 March 20, 2010 14:54


All times are GMT -4. The time now is 04:24.