Divergence detected in AMG solver
1 Attachment(s)
Hi everyone,
I'm trying to simulate a rotating machine and i'm having some dificulties to get the convergence in the solution. I simulate the rotation using Multiple Reference Frame (MRF). The rotor has an angular velocity of 4000 rpm, and the stator is stationary. The mesh hasnīt shown any problems when i do the mesh check. The boundaries are defined as pressureinlet and pressureoutlet, and the walls are the following: Wallstator: The walls of the stator are defined as stationary Wallrotorstatic: Walls of the rotor that don`t rotate, defined as moving walls, absolute, with 0 velocity. Wallrotordinamic: Rotating walls of the rotor, defined as moving walls, relative to adjacent cell zone, with 0 velocity. When I do the simulation i got the following message: Error: Divergence detected in AMG solver: epsilon Error Object: #f I've uploaded an screenshot of the final result (as you can see, the divergence of the continuity starts growing in the 30th iteration more or less). Could anybody help me? Thanks. 
Hi:
First of all: which solver are you using? Segregated or Coupled? Which numerical schemes? First Order or higher order? What kind of initialization did you do? If you are usin second order or higher order schemes of discretization try to use first order schemes and then switch to higher order schemes. Second, if you are using standard or hybrid initialization try to use FMG initialization. Hope this helps you. Regards. Gonzalo 
Quote:
I'm using SIMPLE method and LEAST SQUARE CELL BASED, with all the parameters in First order. If tried with hybrid initialization and FMG initialization, but both of them got me to a similar result (the one that it's shown in the figure). I'm a bit desperate because i started this simulation some weeks ago and it`s been imposible to get a result. Thanks for the support. 
OK, try to use a Coupled solver with the pseudotransient algorithm and play with the timescale (rise it) if you have problems with convergence. Be carefull with it because it requieres more memory than the segregated solver. Reorder the domain to improve the memory usage. Regards
Gonzalo 
Your continuity residual is going up when it is diverging which means you will have high velocity spikes in your calculation. Stop the calculation before it dies and figure out where that is occuring. You might be able to see whether it is a boundary condition problem or mesh problem or, if it occurs middle of the domain and mesh is fine, selection of solver problem.
Also you have high number of cells hitting the viscosity limiter. Check the material property you are using. Try using (incompressible) ideal gas if you are using a constant density. If the density is very much off the solution might diverge. 
1 Attachment(s)
I've done some changes as you recomended me. And I`ve obtained better solution, but i don't reach the convergence criteriums i've set (103 in all the parameters).
My final solution is shown in the picture i'm attaching. (I am now trying to solve it with Coupled method, i hope this could get me to the solution). Thanks for the previous answers, they helped me a lot. 
Quote:

I see that you are using fluent 3d solver rather than 3ddp (3d double precision). Try that. Which means when you launch the Fluent you have to select 3ddp rather than 3d.
Also once the solution has converged with the firstorder discretization, rerun with the secondorder discretization. By the way, you are using standard Kepsilon turbulence model but unless you have a particular reason the Realizable kepsilon should be a better option for you. 
Quote:
Gonzalo 
Hi,
I've tried with double precision, but it still do not converge. I also tried with Standard kepsilon and realizable, but both of them get me to a similar solution. The thing is that my model is a machine rotating at 4000 rpm. To see the convergence, i change the rotating velocity to 100 rpm and it converged. But when i raise the velocity (even to 200 rpm) it doesn't converge. It stays near the convergence with 200 rpm. With 4000 rpm i got 100 m/s in some cells, should it be a problem of air compresibility? Does anyone have an idea? Thanks for the support. PD: Sorry for my english. 
Quote:
mmmmm, 4000rpm is a lot for me. What's the external radius of your compressor? With it you can calculate the velocity at the tip of the blades and vectorially adding the inlet velocity you can calculate the total velocity at these points and whit the local properties of the air calculate the mach number and decide either if the flow is compressible or not. with regard with your convergence problem, what is the minimal Orthogonal Quality of your mesh? There exis big diferences betwen adyacent cell sizes? If you're using BL meshes, firs try to get a converged solution with a mesh without it or a coarse mesh, then interpolate the data to the finner mesh. Regards Gonzalo 
Quote:
The max skeweness is 0.94 (in a few cells), with a skeweness average of 0.24 more or less. I think the mesh is not the problem, but i will check it again to be sure. Should i put the air in a compresible state? Thanks for your answer 
I don't think so, one more thing: what is the turbulence level and the lenght scale you're using at the inlet?

Quote:
And length scale: 0.01 m 
One thing you can do is run the case until it diverges, then stops it and make contours of velocity and p and look for zones where exits singularities. Then you can try to refine the mesh on that zones or improve the quality. I don't know, it's a little difficult to see what is happening without more specifications. Sorry I can't help you.
Gonzalo 
I hope I haven't been too late for replying you :p
You have to use a mixing plane model for the contact region of rotor ans stator. Needless to say that you need two separate fluid zones. You can also go to ANSYS customer portal and get the official tutorial for precisely what you want. There's a tutorial. 
All times are GMT 4. The time now is 15:53. 