[kamakura117]

Hello,

I am using the turbulentTemperatureCoupledBaffleMixed boundary condition with chtMultiRegionFoam and OpenFOAM v 2.3.1. I am running the solution until a steady-state condition is reached to verify with steady-state solutions.

I have two blocks in contact; one block contains a constant volumetric heat generation (first image, block aligned with x=0 has a source of 1e6 w/m^3). In the first case there is no contact resistance and the results match with analytical calculations (second image). At the x = 2*L (with L being the length of one block), I have a constant wall temperature. The other boundaries effectively make this a 1D problem.

However, now I apply a contact resistance:

box_1_to_box_2
    {
        type            compressible::turbulentTemperatureCoupledBaffleMixed;
        value           uniform 310.0;
        Tnbr            T;
        kappa           solidThermo;
        kappaName       none;
        thicknessLayers (1e-5);
        kappaLayers     (1e-2);
    }

Such that the resulting temperature drop, with my source, should be 20 K. When I run the case, I get this temperature drop, roughly. However, the peak temperature of the system is significantly higher than it should be. MOREOVER, when I refine the mesh, the temperature distribution just climbs(third image). In this image, two of the many meshes are shown, a tet mesh and hex mesh. The high temperature distribution is for the tet mesh. The "peak" temperature should be 686 K.

I am just starting to diagnose why this is happening, but I won't be able to get to it in the near term. I just wanted to put this out there.

[kamakura117]

I just submitted a bug report:

http://www.openfoam.org/mantisbt/view.php?id=1682