Heat flux boundary condition

[Martin_]

Dear Elmer users,

I read some posts on the forum regarding the flux boundary condition for the heat equation and this is what I've understood so far:

- when I solve the heat equation in 2 bodies and I don't specify and boundary condition on the interface between them, then continuity of heat flux and temperature is automatically imposed
- when I want to have discontinuity in temperature at the interface, I need two sets of nodes. Say the interface has label 10, I need to run ElmerGrid with -discont 10. The heat flux on the interface will remain continuous.

Can you please confirm this? I would like to know, for example, if it's necessary to specify

Code: Select all

Heat Flux BC = true

(with no further parameters) to have continuous heat flux? Can/should this be completely omitted?

Thank you.

Best regards,

Martin Vymazal

[raback]

Yes, confirmed. The "Heat Flux BC" as the only definition does not have any effect. The continuity of fields and fluxes (in weak sense) is always striven by the standard FE formulation. -Peter