[Elol]

Hi everybody,

I have a problem understanding this issue, I declared volScalarField as follow;

volScalarField klim
(
    IOobject
    (
        "klim",
        runTime.timeName(),
        mesh_
       // IOobject::MUST_READ,
       // IOobject::AUTO_WRITE
    ),
    mesh_
   //deltaX[i] *deltaX[i]/runTime.deltaT() 
//totalDeltaX1/runTime.deltaT()

);

and calculated it as follows;

scalarField deltaX = cbrt(mesh_.V());
forAll (deltaX, i)
{
//deltaX[i] = mag(deltaX[i]);
klim = deltaX[i] * deltaX[i]/runTime.deltaT();

}

But then I have this error;

--> FOAM FATAL ERROR: 
cannot find file "/home/Elol/OpenFOAM/Elol-v1712/run/tests_1/0.011976/klim"

I am using klim in a transport equation.

Anybody knows why this is happening ?

Thanks in advance.

[geth03]

so klim is a volScalarField and is accessed incorrectly in the loop?

klim = deltaX[i] * deltaX[i]/runTime.deltaT(); //klim[i] is correct access

[Elol]

Thanks for you reply,

Yes klim is volScalarField,

when I tried klim[i] it gives me a compilation error as follow

error: cannot convert ‘Foam::dimensionedScalar {aka Foam::dimensioned<double>}’ to ‘double’ in assignment
 klim[i] = deltaX[i] * deltaX[i]/runTime.deltaT();

[geth03]

Originally Posted by Elol

error: cannot convert ‘Foam::dimensionedScalar {aka Foam::dimensioned<double>}’ to ‘double’ in assignment
 klim[i] = deltaX[i] * deltaX[i]/runTime.deltaT();

obviously you didn't define the dimension of klim, so openfoam checks if dimensions match before doing calculations. so when declaring klim, also define the dimension!

[Elol]

Thanks.

It worked, it needed also to change the runTime.deltaT() to runTime.deltaT().value()

Have a nice day.

Best Regards
Elol