Thermal resistance layer(s) not accounted for in turbulentTemperatureRadCoupledMixedFvPatchScalarField

Summary

The thermal resistance layers are not properly taken into account in the turbulentTemperatureRadCoupledMixed boundary condition (possibly the same problem in other boundary conditions as well).

Steps to reproduce

Add a thermal resistance, either as contactLayer or as contactLayers, to a multi-region heat heat transfer case.

Example case

Any multi-region tutorial will do.

What is the current bug behaviour?

The temperatures do not change when adding the resistance layer(s).

What is the expected correct behavior?

A temperature difference between the two regions at the interface. A change in temperatures of the interface and regions

Environment information

• OpenFOAM version : v2112
• Operating system : CentOS
• Compiler : gcc

Possible fixes

Possible fix attached as patch file (only tested in non-implicit mode). It was set up based on the v2012 version of the code.turbulentTemperatureRadCoupledMixedFvPatchScalarField.patch