Matthew Knepley <knep...@gmail.com> writes: > Hmm, let's start here because that does not make sense to me. Can you write > down the mathematical model? Usually I think of temperature as being > a continuous field. I have only seen discontinuities for rareified gases.
Well, you can still model it as a discontinuous field with a flux derived from a jump condition similar to DG. I think that's a reasonable choice when the fluid uses conservative variables and/or a function space that does not have boundary nodes (like FV or DG). But you should check the literature to identify a formulation you're comfortable with and then we can advise on how to implement it with DMPlex.