
[Sponsors] 
July 15, 2010, 23:47 
Convergence question with regards to discretization

#1 
New Member
Karan Anand
Join Date: Feb 2010
Posts: 23
Rep Power: 9 
I have a model where the turbine blades are cooled internally. The discretization to NS equation was first performed with 1st order upwind for momentum, energy and viscous model chosen and standard for pressure. The gradient was calculated with Least squares Cell based. the solution converged with residual set to 1e4 for all except 1e6 for energy.
After the solution converged, to get more accurate results, I switched to Second order upwind for all except pressure which was left as standard with gradient calculated with Least squares Cell based. This time however the solution seems stagnant with continuity at 1e2, turbulence at around 5e4 and momentum around 2e4 energy is at 1e5. I started with the converged solution to begin with. Why is it that the solution is not getting converged for second order upwind? Is there a way I can force the solution to converge for this discretization method? 

July 16, 2010, 08:45 

#2 
Senior Member
karine
Join Date: Nov 2009
Posts: 158
Rep Power: 9 
Hi
1u are doing a DNS (laminar model in Fluent)?? if yes how much cells do u have, whar is ure Re 2u are using steady solver?? if yes is the flow really steady??? I need to know ure answers so i can know what is the problem. FOr convergence, it is not like this that u must do: try to monitor variables in the flow by exemple Vz in point P. For a steady solver, u will get convergence when ure variables dont change anymore. 

July 16, 2010, 12:10 

#3  
New Member
Karan Anand
Join Date: Feb 2010
Posts: 23
Rep Power: 9 
Quote:
Turbulence model: kw SST. At inlet (which is the hub of the blade) the velocity is 1.2266 m/s. Other operating conditions/Boundary conditions: Operating pressure: 1.435 MPa Air temp: 500 K Airfoil temp: 300 K (instead of cooling, heating is applied as it is easy to converge due to change in viscosity of air, this is done to study the effect of heat transfer coefficient) Outlet: OUTFLOW Air properties: Viscosity (µ): 2.6375 x 105 kg/ms Thermal Conductivity: 0.040284 W/mK Specific heat: 1030.305 J/kgK Density of air: incompressibleideal gas (using ideal gas law, around the inlet, density is 10 kg/m3) Quote:
I need the solution for a steady state. The model has no significant curvature , nor high values of natural convection. Grid is Hybrid with hex core and tetra around. Have a prism boundary layer at the airfoil surface. Quote:
Is absolute pressure or velocity at outlet a good option? should I also monitor local surface temperature or Nu number or heat transfer coefficient at point? are these good enough variables to check for convergence? Thanks for your help. 

July 16, 2010, 13:02 

#4 
Senior Member
karine
Join Date: Nov 2009
Posts: 158
Rep Power: 9 
can u send a photo of ure geometry? tell me also if u have highly deformed mesh.
for monitoring, i suggest that u monitor 23 variables, and make at least one of these in the interior domain. In fact, it is possible that ure 'steady  state' include some coherent motion that persists. While using hgh order descretisatin scheme (and of course a better turbulence model, lessly diffusive) ure simulation can better reproduce those strctures. So ure residuals and monitors will include a cyclic behavior and it will not conerge. Tell me what do u get when u monitor variables, and u can do the same thing for ure simulation with 1st upwind, that maybe did not converge too 

July 16, 2010, 13:16 

#5 
New Member
Karan Anand
Join Date: Feb 2010
Posts: 23
Rep Power: 9 
As of now I am using quite a coarse mesh:
http://img85.imageshack.us/gal.php?g...meshstrctu.jpg When i tried solving without monitoring variables except the default ones in fluent and put it to 104, the first order converged. (I didnt see the variables as such so ill plot them to check.) However, after I changed to 2nd order, I did not see any cyclic behavior (again for momentum, energy and cont). It was just that it was steady and converging towards the residuals. Ill check using individual variables anyway. Also could the grid`s coarseness be one issue? 

July 16, 2010, 13:33 

#6 
Senior Member
karine
Join Date: Nov 2009
Posts: 158
Rep Power: 9 
Even if residuals drop under 10^4, it is not sure that u converged. If residuals of continuty are >10^2, it does not mean that u didnt converge.
Check the variables. For the mesh, where there is coarse elements is it solid or fluid elements??? Anyway i dont think the mesh is a problem. Generally u have more chance to converge if the mesh is coarse (but with a bad solution...U must refine to verify that ure solution is independant of grid size). And ure mesh seems to be reasonnably good. I advice u to use hexa in place of tetra elements (use non conformal grid) 

July 16, 2010, 13:46 

#7 
Senior Member
Chris
Join Date: Jul 2009
Location: Ohio, USA
Posts: 169
Rep Power: 9 
I actually think that the mesh could be problem, particularly where the cell size seems to change very rapidly. I've seen bad grids that run with first order and blow up with second order in fluent. Maybe it's because first order is very dissipative and knocks out everything that could blow out the code. I'm not really sure, but that might have something to do with it.


July 16, 2010, 14:44 

#8  
New Member
Karan Anand
Join Date: Feb 2010
Posts: 23
Rep Power: 9 
Quote:
Quote:
Also in FLUENT, to have non conformal grid: "If you create a single grid with multiple cell zones separated by a nonconformal boundary, you must be sure that each cell zone has a distinct face zone on the nonconformal boundary" Fluent Manual: http://cdlab2.fluid.tuwien.ac.at/LEH...ug/node175.htm Accordingly, I will have to divide the entire zone where I have hexcore to separate zones for each little cylinder and then have non conformal grid. Just to let you know. I am validating the results to published numerical work. The paper seems to have similar grid (conformal) but just had tet elements. I added hexcore to reduce the number of elements and seed things up. 

July 16, 2010, 14:49 

#9  
New Member
Karan Anand
Join Date: Feb 2010
Posts: 23
Rep Power: 9 
Quote:
Also had a question about monitoring variables inside the domain in fluent. I have added surface monitors (monitoring abs pressure using mass avg, velocity using mass avg at outlet and heat transfer coeff at the required walls) I did not get the volume monitoring option. Is the suggestion to pick a point in the domain and monitor that? or is it to monitor averaged values for the domain? 

Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Problems with convergence with an easy system  franzdrs  Main CFD Forum  0  June 15, 2009 18:17 
Question about convergence criteria?  Thomas  FLUENT  7  April 8, 2005 19:38 
help! cyclone convergence question  Lcw  FLUENT  1  January 26, 2005 16:29 
Discretisation / Convergence Question  Johnny B  FLUENT  1  November 15, 2003 15:27 
Convergence Question  Colin  FLUENT  13  May 16, 2003 11:41 