CFD Online Discussion Forums - Viscosity UDF Gives different Results to fluent Model?

Hi All.

I am examining many different viscosity models through UDF's in Fluent. I have written a UDF for the Carreau model even though it is in Fluent. When I run 2 identical simulations, one using my UDF for Carreau Viscosity, the other using Fluents inbuilt model, they give different results. This is implying my UDF's are not working as intended, and are likely wrong for the other models. Is there a reaosn why this UDF doesnt implement the Carreau model? (both models use the same values for max/min shear visco, lambda and n.)

Could this be because of the rounding errors in my code? Would it be more accurate to use "double" instead of "float"?

/* Carreau (p=2) Viscosity Model */

#include "udf.h"

float muinf=0.00345;
float mu0=0.056;
float lambda=3.313;
float p=2;
float n=0.3568;

DEFINE_PROPERTY(cell_viscosity,c,t)
{
real diff;
real shear;
real prod;
real power;
real visco;
diff=(mu0-muinf);
prod=((lambda)*(C_STRAIN_RATE_MAG(c,t)));
shear=(1+(pow(prod,p)));
power=((n-1)/(p));
visco=((muinf)+(diff)*(pow(shear,power)));
return visco;
}