A non-isothermal flow simulation involving non-conformal interface and (discontinuous) temperature imposed along outside walls may produce unrealistic (unphysical) temperature distribution while the flow seems properly calculated. In fact, the temperature at the non-conformal interface exhibits overshoots and undershoots. Changing the dilation coefficient or the volume coefficient for non-conformal boundary conditions does not help. Solution: We suspect that there is a conflict of boundary conditions. Indeed, the set of constraints created at the non-conformal interface may be affected by the additional temperature condition imposed. This may become worse when there is a discontinuity in temperature boundary conditions on both sides of the non-conformal interface. The case is presumably over constrained. Actually, a discontinuity in temperature boundary conditions may cause troubles for the set of constraints between variables and this adds to the fact that discontinuity in temperature boundary conditions usually produce oscillations which are also propagated with the flow. In order to circumvent this, temperature boundary conditions along walls can be replaced by forced convection condition alpha*(T-Talpha) (by selecting “Flux density imposed” under Thermal boundary conditions) where a reasonable value for alpha could be 250 W/m/°C while Talpha corresponds to the initially imposed temperature. For a heat conduction calculation, this has produced a significantly better temperature distribution.