[cfdmarkus]

Hi Tammo,

I am having exactly the same problem with prescribing a mass flux instead of Ubar.

The mas flux through cyclic patches is obviously zero.

Have you found a way to compute the mass flux through your domain?



Regards,
Markus

[ngj]

Hi

I have done something similar once upon a time. As you are having the flowDirection, I believe something like this would give you the flux:

scalar flux = sum(pos(flowDirection & mesh.boundary()[patchID].Sf()) * phi.boundaryField()[patchID])

where patchID is referring to the cyclic patch in question. As the face-normal is pointing in opposite directions on the two cyclic patches the inner product will yield opposite sign, thus the "pos" part of the entry filters all opposite phis out of the summation, and voila

Have fun,

Niels

[cfdmarkus]

Brilliant, it does exactly what I need.

Thank you very much!