[mikeP]

Hi,
I use following two methods in order to get the multiplication of Reynolds number and drag coefficient. (Drag coefficient, Cd, here is assumed constant scalar here for convenience)

Re = (magUr*dp*rhoa)/mua;
CdRe = 0.5*Re;

Re = (magUr*dp*rhoa)/mua;
forAll (Re, celli) 
{
     CdRe[celli] = 0.5*Re[celli];
}

These should give the same volScalarField, CdRe.
However, when I use these two code parts, I get different velocity field results at the end.
Are these two code parts exactly the same? or can they give different results when I, let's say, interpolate them?