So it is trying this when creating the thermophysical model:

h.boundaryField()[fluid_internal]

but then it crashes because it is a mixedPatch.