Surfaces output for two regions.
 

Hi,

I am doing conjugate heat transfer simulations, where I have two different regions. At those separate regions, I want to write temperature fields for specified surfaces. The problem is, that no matter how I write my code, OpenFOAM outputs only surfaces of one regions.

Please find the surfaces files in the attachment.

In controldict I simply put:
functions
{
include "surfaces"
include "surfaces2"
}

Thanks for the suggestions.

EDIT:
The above code works. If anyone needs it, take it.

What you've noted is quite true and is a missing feature in OpenFOAM. Supporting it would take some additional work, here's my thoughts or explanation.
The surface sampling are built on top of region/mesh function object, which means their "natural" context is a single mesh region. For sampled patches (for an example) sampling on multiple regions doesn't make much sense. For samplers such as an isoSurface, cutting plane or a triSurface, covering multi-regions makes sense, but still has other issues. Probably cannot sample quantities that only exist in one region (eg, U) - need a zero of undefined for the other region? The various mesh searchers are for a single domain only.

So with these various issues, which sampler types do you think are absolutely needed to handle multi-domain?