CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Meshing & Mesh Conversion

How to get mesh total volume quickly?

Register Blogs Members List Search Today's Posts Mark Forums Read

Like Tree1Likes
  • 1 Post By alexeym

Reply
 
LinkBack Thread Tools Display Modes
Old   August 1, 2012, 15:24
Default How to get mesh total volume quickly?
  #1
Member
 
Join Date: Sep 2011
Posts: 45
Rep Power: 7
ic3wall is on a distinguished road
Hi,

What's the fastest way to get the total volume of a mesh precisely (sum of volume of each cell)?

I know checkMesh is giving Mesh total volume but i'm not sure if it's really mesh total volume or volume of bounding box ..

Thank you!
ic3wall is offline   Reply With Quote

Old   August 1, 2012, 17:24
Default
  #2
Senior Member
 
Travis Carrigan
Join Date: Jul 2010
Location: Arlington, TX
Posts: 144
Rep Power: 8
tcarrigan is on a distinguished road
What software are you using to generate the mesh? It may provide the appropriate utilities for doing a volume calculation.
tcarrigan is offline   Reply With Quote

Old   August 1, 2012, 19:20
Default
  #3
Member
 
Join Date: Feb 2012
Posts: 57
Rep Power: 6
matt.mech.eng is on a distinguished road
checkMesh gives the volume of the mesh not the bounding box
you can easily calculate the volume of the bounding box yourself.. compare with checkMesh volume given
matt.mech.eng is offline   Reply With Quote

Old   January 21, 2016, 13:34
Default
  #4
Member
 
Peter
Join Date: Nov 2015
Location: Hamburg, Germany
Posts: 37
Rep Power: 3
potentialFoam is on a distinguished road
Dear Foamers,

I would like to get the total volume of a solution domain.
I tried
Code:
scalar totalVolume(0.0);
forAll(mesh_.cells(),i) //cellI)
{
totalVolume += mesh_.V()[i]; //cellI];
}
without success. (Wrong result, but I don't know why.)

Then I found this last answer and hence I tried to get it with 'checkMesh'.
Unfortunately, I can't find the code snippet in the folder
/OpenFOAM-3.0.1/applications/utilities/mesh/manipulation/checkMesh/
using e.g.:
grep -nr "volume"

Can you please show me, where they hid it?

Regards,
Peter
potentialFoam is offline   Reply With Quote

Old   January 21, 2016, 16:14
Default
  #5
Senior Member
 
Alexey Matveichev
Join Date: Aug 2011
Location: Nancy, France
Posts: 1,438
Rep Power: 25
alexeym will become famous soon enoughalexeym will become famous soon enough
Hi,

Why have you decided to iterate over strange cells() list? This does the job:

Code:
scalar vol = 0.0;
forAll(mesh.V(), idx)
{
    vol += mesh.V()[idx];
}
reduce(vol, sumOp<scalar>());
Info<< vol << endl;
The same thing:

Code:
scalar vol = gSum(mesh.V());
Info<< vol << endl;
(in fact gSum is just short way to say
Code:
scalar vol = sum(mesh.V());
reduce(vol, sumOp<scalar>());
)

Also it is up to you to check if you run 2D axisymmetric case and multiply result to get real mesh volume.

Answering you last question: it is in primitiveMeshCheck.C (http://foam.sourceforge.net/docs/cpp...ce.html#l00318).
potentialFoam likes this.

Last edited by alexeym; January 22, 2016 at 18:27.
alexeym is offline   Reply With Quote

Old   January 22, 2016, 04:57
Default
  #6
Member
 
Peter
Join Date: Nov 2015
Location: Hamburg, Germany
Posts: 37
Rep Power: 3
potentialFoam is on a distinguished road
Thanks Alexey,

it works
(I use
Code:
scalar vol = gSum( mesh_.V() );
)
potentialFoam is offline   Reply With Quote

Reply

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
[ICEM] Negative volume error in hybrid mesh siw ANSYS Meshing & Geometry 4 September 3, 2014 05:25
Moving mesh Niklas Wikstrom (Wikstrom) OpenFOAM Running, Solving & CFD 122 June 15, 2014 06:20
3D Hybrid Mesh Errors DarrenC ANSYS Meshing & Geometry 11 August 5, 2013 06:42
channelFoam for a 3D pipe AlmostSurelyRob OpenFOAM 3 June 24, 2011 13:06
ICEM 12 CFD help creating volume mesh from stl EmpError ANSYS 0 November 13, 2010 07:38


All times are GMT -4. The time now is 22:12.