|
[Sponsors] |
RNG diverged during the analysis the flow over a multi element airfoil, why? |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
August 5, 2013, 08:39 |
RNG diverged during the analysis the flow over a multi element airfoil, why?
|
#1 |
Senior Member
saeideh mohamadi
Join Date: Aug 2012
Posts: 229
Rep Power: 15 |
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 |
|
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Problem with airfoil shape optimization | robyTKD | SU2 Shape Design | 7 | March 7, 2022 16:18 |
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 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 13:54 |