Flux condition on an internal surface
16 days ago by
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:
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
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.