[jipaz]

Hi Everybody:

I'm simulating a 2D Fluidized Bed Riser in Ansys Fluent 14.5 using an Eulerian approach, k-e turbulent model, Johnson-Jackson wall boundary condition. When I use some drag model included into fluent, I havo no problem about the convergence, but I can not obtain the proper mass flux profiles into the riser, so I tried to use the EMMS drag model as an UDF, I compiled it well, without any error message, but when I started the calculation procedure, before the first iteration I got the following error messages:

Updating solution at time levels N and N-1.
done.
iter continuity u-aire u-solidos v-aire v-solidos gt-solidos vf-solidos time/iter

mp-x-momentum equations:

Divergence detected in AMG solver: pressure correction
Divergence detected in AMG solver: vof-1
Divergence detected in AMG solver: pressure correction
Divergence detected in AMG solver: vof-1
Divergence detected in AMG solver: pressure correction
Divergence detected in AMG solver: vof-1
Divergence detected in AMG solver: pressure correction
Divergence detected in AMG solver: vof-1
Primitive Error at Node 0: floating point exception

Primitive Error at Node 1: floating point exception

Primitive Error at Node 2: floating point exception

Primitive Error at Node 3: floating point exception

I test my problem using the modified Syamlal-O'Brien drag model example, and the run progress properly, so I understand that the problem is the EMMS UDF, but I don't know what I should do. Any help is very welcome.

Regards