How can I add a thermal contact resistance (or surface thermal conductance) in heat PDEs?

[ChenghanKevinWang]

I am a beginner of FEniCSx, and I am trying to use FEniCSx to solve heat PDEs. But I have a question:
Imagine that I have two stacked cuboids and they have an horizontal interface. How can I add a thermal contact resistance (or surface thermal conductance) between the two cuboids? Can this be done by modifying the weak form? For example, Can I create two function spaces and using Robin-like weak form to link the two surfaces. Thanks!!!

[dokken]

You can do this with something along the lines of

Discontinuity at interface using mixed domains

Hi all, Exactly like in this topic I am trying to impose a discontinuity at an interface between two subdomains. This time having two subdomains with only CG elements and coupling them with a dS term. @dokken published this gist that I am trying to adapt in order to have: \frac{u(-)}{K(-)} = \frac{u(+)}{K(+)} \nabla u(-) - \nabla u(+) = 0 However I am struggling to adapt the interface terms to achieve the desired result. This is the DG formulation we had earlier: # Interface F += - dot…

For instance a modified version of this for is available at Thermal contact resistance problem · GitHub

[ChenghanKevinWang]

I tried the rtc_dg and it works very very well. Thank you man!!!