[froberto]

Hi all,
I am trying to model a problem consisting of a solid (sub-divided into two regions) and a fluid.
The two solid regions produce heat at different rates, and they have the same material properties.

I dug into this forum for quite a while and found two posts:
https://www.cfd-online.com/Forums/openfoam-pre-processing/211078-solid-solid-thermal-resistance-layer-thickness.html
https://www.cfd-online.com/Forums/openfoam-solving/151469-boundary-condition-conjugate-thermal-heat-conduction.html

They seem to agree on the fact that solid-to-solid interfaces should be modeled using compressible::turbulentTemperatureCoupledBaffleMix ed.
However, my result doesn't look good at the interface (see attached figure).

I have attached all the files.
Any ideas of what is going on/suggestions are welcome!
Thank you,


Here is a snippet of the solid "heater0" changeDictionaryDict:

        heater0_to_heater1
        {
            type            compressible::turbulentTemperatureCoupledBaffleMixed;
            Tnbr            T;
            kappaMethod     solidThermo;
            value           uniform 298;
        }

Here is a snippet of the solid "heater1" changeDictionaryDict:

        heater1_to_heater0
        {
            type            compressible::turbulentTemperatureCoupledBaffleMixed;
            Tnbr            T;
            kappaMethod     solidThermo;
            value           uniform 298;
        }

And here is a snippet of the topoSetDict:

    // Heater 0
    {
        name    heater0CellSet;
        type    cellSet;
        action  new;
        source  boxToCell;
        box     (-0.004 -0.05 -0.1) (0.004 -0.025 0.1);
    }
    {
        name    heater0;
        type    cellZoneSet;
        action  new;
        source  setToCellZone;
        set     heater0CellSet;
    }
    // Heater 1
    {
        name    heater1CellSet;
        type    cellSet;
        action  new;
        source  boxToCell;
        box     (-0.004 -0.025 -0.1) (0.004 0.0 0.1);
    }
    {
        name    heater1;
        type    cellZoneSet;
        action  new;
        source  setToCellZone;
        set     heater1CellSet;
    }

Attachments:
of-T-final-line.png
post.zip

[froberto]

When using the sample utility I have the possibility to select "cellPointFace" interpolationScheme.
The plots now look right.
I am assuming, paraview plots the cell values by default or something like that.

Anyway, I thought it was mentioning that the following post is related:
chtMultiregionFoam--unequal temperature at coupled patches