[zmester]

Hi

Is there anyone who can explain me how do the cell centroid(c_centroid) and cell volume(C_Volume) functions of Fluent work? What method they use to calculate the centroid and the volume of a cell. Im writing my code, and would like to create functions similar to these.Thanks.

[esozer]

I don't how fluent does this but the following will work for arbitrary cell geometries:

Loop faces
   Loop face nodes
     Form a triangle using two subsequent face nodes and a point within the face (e.g  average of face nodes)
     Calculate the resulting triangle area
     Add to face area
     Form a tetra using the triangle and a point within the cell volume
     Calculate tetra volume
     Add to cell volume
     Tetra centroid is known (google it)
     cell_centroid=cell_centroid+tetra_volume*tetra_centroid
   End face node loop
End face loop

cell_centroid=cell_centroid/cell_volume

Hope this helps.

[zmester]

Thank you for your answer. Since i asked i have already find an other method using the divergence theorem (Gauss) to transform volume integrals to surface integrals.

[esozer]

Right, that will work too. I preferred the method I described since it will give the face areas, face centroids, cell volumes and cell centroids in one swipe.

Good luck with your code.