[sharonyue]

Hi guys,

I have a volScalarField,

volScalarField epsilonD(pow(epsilon/pow(d,2),1/3));

you can see that epsilon's dimension is [0 2 -3 0 0 0 0] and d's dimension is [0 1 0 0 0 0 0], So epsilonD's dimension should be [0 0 -1 0 0 0 0]

But openfoam tell me its dimention is [0 0 -0 0 0 0 0]?

Any ideas?

[sharonyue]

volScalarField epsilonD(pow(epsilonF,1/3));

This is very simple, epsilonF's dimension is [0 0 -3 0 0 0 0] , but epsilonD's dimention is [0 0 -0 0 0 0 0]

[sharonyue]

Solved.

volScalarField epsilonD(pow(epsilonF,1.0/3.0));

I checked it out in dimensionSet.C and I got this piece of code:

Foam::dimensionSet Foam::cbrt(const dimensionSet& ds)
{
    return pow(ds, 1.0/3.0);
}
..........................

Foam::dimensionSet Foam::pow(const dimensionSet& ds, const scalar p)
{
    dimensionSet dimPow
    (
        ds[dimensionSet::MASS]*p,
        ds[dimensionSet::LENGTH]*p,
        ds[dimensionSet::TIME]*p,
        ds[dimensionSet::TEMPERATURE]*p,
        ds[dimensionSet::MOLES]*p,
        ds[dimensionSet::CURRENT]*p,
        ds[dimensionSet::LUMINOUS_INTENSITY]*p
    );

    return dimPow;
}

I dont know if its a bug. But its solved right now. Any comments are welcomed.

[bieshuxuhe]

int a=1;
int b=3;
float c=a/b;

when a and b both are int, c=0

best wishes!