Thermal equilibrium is given when the temperature is constant in the entire system. In
the present work, local thermal equilibrium is assumed, i.e. the temperature of all phases
including the matrix are the same at a finite volume in space. Therefore, local temperature
differences, e.g. as a result of slow heat transfer from the fluid to the matrix, are neglected.
If flow velocities are slow, this is a reasonable assumption