[mikemech]

Hello,

I'm facing this problem:

I want to calculate the volume of a certain cell that is contained in my grid.
Is there any straightforward method for doing that?

Thank you.

[romant]

you could create a field for all the cell volumes and then access your specific cell the field is created with

    volScalarField cellVolu
    (
    IOobject
       (
        "cellVolu",
        runTime.timeName(),
        mesh,
        IOobject::NO_READ,
        IOobject::AUTO_WRITE
       ),
        mesh,
    dimensionedScalar("zero", dimVolume, 0.0)
    );
    cellVolu.internalField() = mesh.V();

and access to the volume of a specific cell is then given by

cellVolu()[cellIndex];

or right away by

mesh.V()[cellIndex]

have a look at similar things
http://www.cfd-online.com/Forums/ope...ne-volume.html

[mikemech]

Hi romant,

You mean that if I type in the terminal window "mesh.V()[34332]" I ll get the volume of cell number 34332 ?

[romant]

no this should have been in C++ in the source code.

[tomislav_maric]

Originally Posted by mikemech

Hello,

I'm facing this problem:

I want to calculate the volume of a certain cell that is contained in my grid.
Is there any straightforward method for doing that?

Thank you.

There is a Foam::cell::mag() method for this. It is local and approximative:

cell.C:

00238 Foam::scalar Foam::cell::mag
00239 (
00240     const pointField& p,
00241     const faceUList& f
00242 ) const
00243 {
00244     // When one wants to access the cell centre and magnitude, the
00245     // functionality on the mesh level should be used in preference to the
00246     // functions provided here. They do not rely to the functionality
00247     // implemented here, provide additional checking and are more efficient.
00248     // The cell::centre and cell::mag functions may be removed in the future.
00249 
00250     // WARNING! See cell::centre

There are more accurate ways of computing cell magnitude (polyhedron volume)
than this one, and you should check out the literature if you absolutely need local computations. Do you need the volume of a single cell (set of cells), or do you need mesh level operations?

If you just want a volume of a single cell, or a zone of cells defined by some
criterion, you can go about it like this:

const cellList& cells = mesh.cells();
const faceList& faces = mesh.faces();
const pointField& points = mesh.points();

// Let's say that you want a volume of the 100th cell..

scalar cellV = cells[100].mag(points,faces);

Of course, finding out the proper label for the cell(s) will depend on what you
actually want to do.

T.

[mikemech]

Thank you tomislav_maric,
Well I need to find the volume of one cell, and I know its label.
So I ll try the second idea you mentioned.

[tomislav_maric]

Originally Posted by mikemech

Thank you tomislav_maric,
Well I need to find the volume of one cell, and I know its label.
So I ll try the second idea you mentioned.

Cool, I'm glad to be of help!

T.

[oort]

Hi!
If i put

cellVolu.write();

bellow the "cellVolu.internalField() = mesh.V();" is written the cellVolu field.
How can I sum all the cell volumes? I want calculate the global volume of my geometry.

Thanks.

Originally Posted by romant

you could create a field for all the cell volumes and then access your specific cell the field is created with

    volScalarField cellVolu
    (
    IOobject
       (
        "cellVolu",
        runTime.timeName(),
        mesh,
        IOobject::NO_READ,
        IOobject::AUTO_WRITE
       ),
        mesh,
    dimensionedScalar("zero", dimVolume, 0.0)
    );
    cellVolu.internalField() = mesh.V();

and access to the volume of a specific cell is then given by

cellVolu()[cellIndex];

or right away by

mesh.V()[cellIndex]

have a look at similar things
http://www.cfd-online.com/Forums/ope...ne-volume.html

[romant]

Hej,

it is easiest calculated

scalar volumeSum = sum(cellVolu).value();

or if you don't need it in a program, just run

checkMesh

which will give you the total volume as one of the outputs.

[oort]

Thanks romant.
The total cell's volume given by your method have deviation from "geometric" calculations of 0.002% in my mesh.

[pvpnr]

Originally Posted by oort

Hi!
If i put

cellVolu.write();

bellow the "cellVolu.internalField() = mesh.V();" is written the cellVolu field.
How can I sum all the cell volumes? I want calculate the global volume of my geometry.

Thanks.

in which source code file do we have to make these changes