
[Sponsors] 
March 18, 2013, 18:32 
Zero iteration step, why??

#1 
Senior Member
Join Date: Nov 2012
Posts: 168
Rep Power: 5 
Hi All,
In one of my simulations, I used kepsilon model, and got the following output: ............... pEqn.H_1: rho max/min : 1.1966 1.19659 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 time step continuity errors (Low Ma): sum local = 8.61908e16, global = 2.16302e17, cumulative = 2.87666e15 pEqn.H_2: rho max/min : 1.19661 1.1966 GAMG: Solving for phd, Initial residual = 0.00121095, Final residual = 8.07744e06, No Iterations 8 GAMG: Solving for phd, Initial residual = 0.00275727, Final residual = 2.14566e05, No Iterations 5 GAMG: Solving for phd, Initial residual = 0.00166235, Final residual = 9.90116e06, No Iterations 4 pEqn.H_1: rho max/min : 1.19661 1.1966 diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 time step continuity errors (Low Ma): sum local = 7.21639e16, global = 2.3332e17, cumulative = 2.89999e15 pEqn.H_2: rho max/min : 1.1966 1.19659 DILUPBiCG: Solving for epsilon, Initial residual = 5.56717e07, Final residual = 5.56717e07, No Iterations 0 DILUPBiCG: Solving for k, Initial residual = 4.73169e07, Final residual = 4.73169e07, No Iterations 0 Pimple Loop End: rho max/min : 1.1966 1.1966 Pimple Loop End: nu max/min : 1.528e05 1.528e05 ExecutionTime = 1907.55 s ClockTime = 1915 s .................... From above, we can see the number of the iterations is zero. What does it mean? Does it mean k and epsilon are solved using explicit method? (I have this think because in the compressible solvers, the density is solved using explicit method with zero iteration step). Why are they (k and epsilon) are solved with zero iteration steps? In the beginning of the simulations, they are not zero, but 35 steps. Does anybody give me some hints about this? best H 

March 18, 2013, 19:05 

#2 
Senior Member

Hi H,
The number of iteration is 0 since the initial residual is less than the convergence criteria you set. So for all intent and purpose, the solver is assuming that the epsilon and k fields are "converged" and do not need iterations. If you do want to have more iteration, decrease the convergence criteria for k and epsilon. Julien
__________________
 Julien de Charentenay 

March 18, 2013, 19:57 

#3  
Senior Member
Join Date: Nov 2012
Posts: 168
Rep Power: 5 
Hi Julien,
When I change the tolerance in system/fvSolution as follows: k { solver smoothSolver; smoother GaussSeidel; tolerance 1e9; relTol 0; nSweeps 1; /* solver PBiCG; preconditioner DILU; tolerance 1e010; relTol 0.1; */ }; epsilon { solver smoothSolver; smoother GaussSeidel; tolerance 1e9; relTol 0; nSweeps 1; /* solver PBiCG; preconditioner DILU; tolerance 1e010; relTol 0.1; */ }; The output from the solver is like the followings: diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 time step continuity errors (Low Ma): sum local = 7.207e16, global = 8.28698e18, cumulative = 7.88457e15 pEqn.H_2: rho max/min : 1.1966 1.19659 DILUPBiCG: Solving for epsilon, Initial residual = 4.13716e07, Final residual = 4.13716e07, No Iterations 0 DILUPBiCG: Solving for k, Initial residual = 5.42019e07, Final residual = 5.42019e07, No Iterations 0 It is a little strang that changing the tolerance does not change the iteration number. Besides, I found that linear algebraic solver DILUPBiCG is not the one I specified in the fvSolution file (smoothSolver). When I change the linear algebraic solver for the velocity U, it indeed changes. So I am a little confused about the controling parameter for k and epsilon. Why can they not be changed in system/fvSolution although they should be? Thank you so much! hz283 Quote:


March 18, 2013, 20:08 

#4  
Senior Member
Join Date: Nov 2012
Posts: 168
Rep Power: 5 
Hi Julien,
If during the simulation the iteration number of k and epsilon is always zero, does it mean that actually the values of k and epsilon are not updated by the solver? Thanks. Quote:


Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Foam::error::PrintStack  almir  OpenFOAM Running, Solving & CFD  62  June 16, 2016 03:01 
InterFoam negative alpha  karasa03  OpenFOAM  7  December 12, 2013 04:41 
why PIMPLE doesn't converge at each time step but no diverging?  immortality  OpenFOAM Running, Solving & CFD  7  May 19, 2013 14:16 
time step and iteration  tony25800  FLUENT  2  February 3, 2012 05:33 
Heat exchanger problem  chiseung  FLUENT  16  October 20, 2001 04:36 