FLUENT 6  Slow convergence along centerline for polar grids (Ogrids)
It is known that FLUENT converges very slowly along the centerline of 3D polar grids. These types of grids are generated in some preprocessors such as GridGen and ICEM. This problem represents some limitations of the FLUENT segregated solver for Otype grids with a high degree of anisotropy. With the standard settings, FLUENT is not able to obtain convergence near the axis. The explanation is the following: FLUENT's segregated solver in STEADY STATE mode requires underrelaxation for the momentum equations to facilitate the work of pressurevelocity coupling schemes like SIMPLE, SIMPLEC, or PISO. This causes very high diagonal dominance in the matrix on coarse levels of the algebraic multigrid (AMG) solver for the momentum equations, and the nearaxis cells are almost completely detached from the other cells inside the AMG. The workaround is to modify the conditions in the solver such that the segregated pressurevelocity coupling algorithm will allow application of momentum underrelaxation equal to 1.0. The simplest method of overcoming this problem is to use the transient solver with the SIMPLEC or PISO pressurevelocity coupling algorithms and an underrelaxation factor of 1.0 for momentum. A more computationally efficient method for steady flow simulations is to include the pseudotransient terms in the transport equations for momentum and turbulence quantities in the form of a userdefined source term. This allows the pseudotransient term to act as some sort of underrelaxation to prevent divergence, instead of the underrelaxation factors in the individual transport equations. To set up such a problem, 1. Establish the flow field everywhere except near the centerline 2. compile the UDF below 3. hook the UDF to the transport equations for U, V, W, k, and epsilon 4. change the PV coupling to SIMPLEC 5. Set all underrelaxation factors to 1 6. Iterate The value of delta_t (the pseudotimestep) may need to be varied to suit the specific conditions of the problem. Also, to speed up convergence, the AMG can be forced to visit more coarse levels by applying the V or F type of AMG cycles for the velocity components (optional).  #include "udf.h" static real delta_t = 1e5; DEFINE_SOURCE(unsteady_source,c,t,dS,eqn) { real rho = C_R(c,t); real source; source = 0.0; dS[eqn] = rho/delta_t; return source; } 

