[karkar]

I would like to use patchAverage to obtain the average pressure on an object that consists of multiple patches (due to different boundary conditions), however patchAverage only averages over single patches. As long as the constituent patches aren't shared with other objects, I could simply take the area-weighted average of each patch's average.

However, due to a dynamic geometry where the individual patches change area in time, this would become quite messy. So instead, my idea is to use topoSet to define a faceSet and turn that into a new object patch via createPatch (replacing the old ones, which are no longer needed after simulating).

Would createPatch correctly map the simulated fields to the new patch or would I end up with values being mapped to the wrong vertices? Or is there another way to achieve this? (I also considered using sample, which can obtain the values for multiple surfaces at once, but then I'd have to manually average, and I'm not sure this way I'd weigh the patch faces correctly.)

(Also asked at the Scientific Computing StackExchange: http://scicomp.stackexchange.com/q/2909/1669)