This technique was originally introduced in a series of seminal papers by Gajewski and Gröger [1]}, [2]}, [3]}. The main differences with respect to those works consist in a different expression of the total free energy and in the presence of nonlinear Robin boundary conditions we have to deal with.
