Divide a surface quantity by the patch area

crmccreary · August 26, 2020, 21:52

I would like to divide a surface quantity by the patch area.
Specifically, I am using MPPICFoam and calculating particle erosion, Q.
Q is in L^3 (m^3).

I would like to divide Q by the patch area to get a new variable in terms of L (m).
I could modify the particle erosion function to create a new variable, but I would then need to re-run the model.

How can I accomplish this post-processing task either using a post-processing function or paraview?

crmccreary · August 28, 2020, 16:08

To answer my own question, I accomplished my objective in two ways.
I modified ParticleErosion.C to divide by the cell face area so that I get wall loss due to erosion.

Using the python calculator in Paraview,
inputs[0].CellData['kinematic:Q']/area(inputs[0])

August 26, 2020, 21:52	Divide a surface quantity by the patch area	#1
crmccreary New Member Charles McCreary Join Date: Jun 2010 Posts: 12 Rep Power: 15	I would like to divide a surface quantity by the patch area. Specifically, I am using MPPICFoam and calculating particle erosion, Q. Q is in L^3 (m^3). I would like to divide Q by the patch area to get a new variable in terms of L (m). I could modify the particle erosion function to create a new variable, but I would then need to re-run the model. How can I accomplish this post-processing task either using a post-processing function or paraview?

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
[OLAFLOW] The OLAFLOW Thread	Phicau	OpenFOAM Community Contributions	457	March 27, 2024 00:59
y+ and u+ values with low-Re RANS turbulence models: utility + testcase	florian_krause	OpenFOAM	114	August 23, 2023 05:37
[Commercial meshers] Mesh conversion problem (fluent3DMeshToFoam)	Aadhavan	OpenFOAM Meshing & Mesh Conversion	2	March 8, 2018 01:47
[snappyHexMesh] Bad meshing result on corrugated metal sheet	UebertreibeR	OpenFOAM Meshing & Mesh Conversion	1	August 19, 2016 04:55
[GAMBIT] periodic faces not matching	Aadhavan	ANSYS Meshing & Geometry	6	August 31, 2013 11:25

August 28, 2020, 16:08		#2
crmccreary New Member Charles McCreary Join Date: Jun 2010 Posts: 12 Rep Power: 15	To answer my own question, I accomplished my objective in two ways. I modified ParticleErosion.C to divide by the cell face area so that I get wall loss due to erosion. Using the python calculator in Paraview, inputs[0].CellData['kinematic:Q']/area(inputs[0])