Convergence and discretization
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 (k-w SST) and standard for pressure. The solution converged with residual set to 1e-4 for all except 1e-6 for energy. I also have mass avrged pressure, velocity and temperature monitors at outlet.
After the solution converged (just with residual criteria, surface monitors are still varying), to get more accurate results, I switched to Second order upwind for all except pressure which was left as standard and reducing the relaxation factors. this time monitoring both residuals and surface monitors. I ran the simulation for a while with 'none' as convergence criteria.
I have a couple of questions:
Firstly, I plotted the residuals and surface monitors. In the velocity monitor, the velocity has almost stabilized, but there exists a few wiggles (extremely minor oscillations). Can I consider the solution to be converged (all other monitors are flat, and converged) ? Is this expected for the discretization method chosen? Also, the residual for energy did not satisfy the 1e-6 mark, can i still say the solution is converged? (I have added a few pics for the grid and convergence history)
Secondly, Could doing this help? ->
Instead of changing the discretization to 2nd order upwind, I first change it to Power law from 1st order upwind. As far as i remember , power law is 1st order accurate but it is way less diffusive and it seems quite similar to hybrid differencing scheme. Is the solution obtained here good enough?
After I switch to power law, I come back to 2nd order, which is a higher order scheme. Could this reduce the wiggles, or help me approach a converged solution? I thought of going for QUICK, but i have a lot of tet meshes.
The solver is Steady state -SIMPLE scheme for the staggered grid. 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. Entire meshed domain is fluid. Grid is conformal.
Images for grid: (coarse mesh)
Images for convergence (a bit refined mesh and better transition between mesh regions-i didn't add images for them, but you get the basic idea from above)
Could anyone please help?
More often than not physical flow is inherently transient, even though steady solver is used.
I would say that you can call your solution "converged" after the 6000th iteration. It's probably still going to do the wiggles for the next 4000 iterations if you let it run on. You may have an issue with your mesh. There are sharp transitions in your mesh density that I don't like. See if you can improve those transitions and run again to 6000 iterations and compare. What's your current mesh size? Basically you might want to do a grid convergence study.
|All times are GMT -4. The time now is 19:46.|