### Flux condition on an internal surface

46

views

0

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$

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!

Rem

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$

`T`^{−}−`T`^{+}=`C``T``R`·`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!

Rem

Community: FEniCS Project

What does your code look like, and where are you stuck?

written
7 days ago by
lhdamiani

### 1 Answer

0

Hi,

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

Remi

Attached is my code.

File attached: 3-MOT.py (13.31 KB)

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

Remi

Attached is my code.

File attached: 3-MOT.py (13.31 KB)

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