Flux condition on an internal surface

16 days ago by
Hi every one

I'm currently solving the poisson equation and I would like to add an internal condition.

Let's say I'm solving the heat equation and I want to add a contact thermal resistance (CTR). This CTR is described as a temperature discontinuity on a surface as follow:
$T^--T^+=CTR\cdot q$TT+=CTR·q
where q is the thermal flux.

Do you have any idea how to implement this kind of condition in FEniCS (let's take the poisson.py demo)?

Thanks a lot for your precious help!

Community: FEniCS Project
What does your code look like, and where are you stuck?
written 7 days ago by lhdamiani  

1 Answer

7 days ago by

My code is able to simulate heat transfer, identify subdomains (volumes and surfaces) based on XDMF attributes. I can easily apply different materials properties and boundary conditions of any kind on the external surfaces. But I would like now to apply a Contact Thermal Resistance (CTR) on an internal surface between two subdomains to model imperfect contact between them.

I don't know how to implement this in the weak formulation.

Thank you very much for your help


Attached is my code.
File attached: 3-MOT.py (13.31 KB)

Please login to add an answer/comment or follow this question.

Similar posts:
Search »