[Pagoda]

Dear Foamers,
I am using the cuttingPlane class to evaluate integrals over a slice of the three-dimensional domain.

I need to use cuttingPlane itself, not function objects, because I have to take the slice only on a subset of the whole domain.

My problem is how to measure the surface of the portion of plane inside a prescribed cell.

To be more clear, given the cell indexes of the domain subset "subSetCells", and a field to be integrated "U", this is the code I use:

\\ plane base point and vector
point pt(0,0,1);                  
dir (0,0,1);
plane pln(pt, dir);

\\ cutting plane triangulated (true), over a mesh subset
cuttingPlane cutPln(pln, mesh, true, subSetCells); 
labelList cutLabels = cutPln.cutCells();


\\ loop over all the cut faces
scalar US(0);
scalar S(0);
forAll(cutLabels, labelI)
{
    scalar ui = U[cutLabels[labelI]];
    scalar Si = ???????     // this should be the surface of the triangulated plane
    US += ui*Si;
    S += Si;
}
US /= S;

Any idea?

Kind regards,
Matteo

[ashish.svm]

I want to calculate the area of the plane of intersection in a grid block.

For example, if blockMesh is used then I want to find out the area of the plane of intersection with a cube.

[Gary]

Hi Have you found the solution to your question?


Regards, Gary

[ashish.svm]

Originally Posted by Gary

Hi Have you found the solution to your question?


Regards, Gary

Not yet

For the time being, I have simplified my problem to focus on simple geometry. But I would be happy to get some idea to solve this