CFD Online Discussion Forums - Getting volume of cells with mesh.V does not work

Hello everyone!

For the last few days, I was toying around with some Algorithms I found in a Paper and tried to put them in an OpenFOAM Solver. Some Calculations require the Volume of Cells. Since all of the other Field Variables are of Type volScalarField or volVectorField, I thought of creating a volScalarField and initializing it with mesh.V().

Google turned up this Thread: http://http://www.cfd-online.com/Forums/openfoam/82866-write-cell-volumes.html
The second Answer provided a Code-Snippet, which seemed to be exactly what I was looking for:

Code:

volScalarField cv     

    (             

    IOobject 

               (                     

        "cv",

        runTime.timeName(),

        mesh,

        IOobject::NO_READ,                     

        IOobject::AUTO_WRITE 

    ), 

    mesh,

    dimensionedScalar("zero",dimVolume,0.0) 

    ); 



cv.internalField() = mesh.V();

However, wmake seems not to like this:

Code:

optimVariables.H:77:24: error: passing ‘const Internal {aka const Foam::DimensionedField<double, Foam::volMesh>}’ as ‘this’ argument of ‘void Foam::DimensionedField<Type, GeoMesh>::operator=(const Foam::DimensionedField<Type, GeoMesh>&) [with Type = double; GeoMesh = Foam::volMesh]’ discards qualifiers [-fpermissive]

     cv.internalField() = mesh.V();

                        ^

I am running OpenFOAM 4.0 under Ubuntu Mate 15.04.

It feels like I am missing something pretty obvious about the usage of someField.internalField() and mesh.V(). Some additional googling told me that mesh.V() returns a scalarField, which resembles a volScalarField without boundaryField and should be able to be written to a volScalarField using the internalField(). Is this correct? Also, is it possible to write to a internalField()? I tried grepping for "internalField() =" in $FOAM_SRC, but found nothing (possibly my grep skills are not the best:)). What would be the correct way to do this?

Also I tried to initialize the internalField to a constant:

Code:

cv.internalField() = 1.0;

which yields using wmake

Code:

optimVariables.H:77:24: error: passing ‘const Internal {aka const Foam::DimensionedField<double, Foam::volMesh>}’ as ‘this’ argument of ‘void Foam::DimensionedField<Type, GeoMesh>::operator=(const Foam::dimensioned<Type>&) [with Type = double; GeoMesh = Foam::volMesh]’ discards qualifiers [-fpermissive]

     cv.internalField() = 1.0;

Googling the error message told me something about the Implementation of the underlaying classes and operators might be off, but I honestly can't believe that. I'd also like to avoid digging deeper in the code than necessary...

Is there any other way to achieve this? Most of my input for those Equations are volScalar/VectorFields, but as long as I get a volScalarField in the end that I can write to file, I would not mind doing things an other way.

confused Greetings,

Juggler