
[Sponsors] 
December 6, 2010, 07:37 
calculate surface area

#1 
Member

Deal OpenFOAM User,
I am implementing PLIC in interFOAM. Then I need to find the way to calculate the volume of each liquid occupied in each cell. I found in the primitiveMesh class two functions that can help me: void Foam: primitiveMesh::makeFaceCentresAndAreas ( const pointField& p, vectorField& fCtrs, vectorField& fAreas ) and void Foam: primitiveMesh::makeCellCentresAndVols ( const vectorField& fCtrs, const vectorField& fAreas, vectorField& cellCtrs, scalarField& cellVols ) However, I have a problem with construct pointField p in function makeFaceCentresAndAreas. Let's say we have four points: (1 0 0) (2 0 0) (2 1 0) (1 2 0) Does anyone have ideas about how to construct variable p from these four points for the input of this function? Regards, Duong 

December 6, 2010, 13:26 

#2 
Senior Member
Santiago Marquez Damian
Join Date: Aug 2009
Location: Santa Fe, Santa Fe, Argentina
Posts: 420
Rep Power: 15 
Hi Duong, what about doing,
volScalarField alpha1Vol=mesh.V()*alpha1; volScalarField alpha2Vol=mesh.V()*(scalar(1)alpha1); or something like this. Best.
__________________
Santiago MÁRQUEZ DAMIÁN, Ph.D. Postdoctoral Fellow Research Center for Computational Mechanics (CIMEC)  CONICET/FICHUNL T.E.: 543424511594 Ext. 1005 Güemes 3450  (3000) Santa Fe Santa Fe  Argentina http://www.cimec.org.ar 

December 7, 2010, 11:03 

#3 
Member

Hi Santiago,
Thanks for the suggestion. However, since I want to construct the interface which has a form n.x+const = 0 (n is the normal vector and x is the point), I should find value of the constant to ensure the conservation of volume fraction. Then, I must find a way to use those function. Btw, I found how to construct this pointField p. What you can do is: pointField& p(4,point(0,0,0)); p[0] = point(1,0,0); p[1] = point(2,1,0); etc.. That's how we define pointField in OpenFOAM. Regards, Duong 

December 7, 2010, 14:00 

#4 
Senior Member
Santiago Marquez Damian
Join Date: Aug 2009
Location: Santa Fe, Santa Fe, Argentina
Posts: 420
Rep Power: 15 
Aha, I see. My answer was related to "then I need to find the way to calculate the volume of each liquid occupied in each cell". I understand that you need to calculate a plane within a cell such as the subvolumes would be equal to those are calculated by the formulas I exposed before...
Bye.
__________________
Santiago MÁRQUEZ DAMIÁN, Ph.D. Postdoctoral Fellow Research Center for Computational Mechanics (CIMEC)  CONICET/FICHUNL T.E.: 543424511594 Ext. 1005 Güemes 3450  (3000) Santa Fe Santa Fe  Argentina http://www.cimec.org.ar 

Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Problem with Gmsh  nishant_hull  Open Source Meshers: Gmsh, Netgen, CGNS, ...  23  August 5, 2015 02:09 
Export Cell Surface Area  Thierros  CFX  0  April 28, 2010 06:01 
cell surface area in boundary..  Chiar  FLUENT  0  March 7, 2007 04:53 
How to calculate the cell area  Le  FLUENT  0  February 18, 2007 23:15 
Compute total surface area of Dpm  figo  FLUENT  0  September 8, 2004 21:44 