[anger]

Hello Foamers,

I want to model a flow with a porous zone which works well using rhoPorousSimpleFoam in OF 1.6.x.
In a next step, I want to add a heat sink to the porous medium. Is it a good idea to use (and probably extend) the porousZone machinery to this end or is it a better idea to do something on the solver level?
A possible approach is to modify the hEqn with an appropriate source term and multiply it with a field value which corresponds to 1 in the porous region and 0 elsewhere, but this seems not to be very elegant.
What do you think?

Best regards,
-Thomas