We develop a phenomenological, computation model based on a system of interacting spin-pairs to examine the thermodynamics of magnetization kinetics in permanent magnets. In particular, we examine the evolution of the system across the magnetic Gibbs free energy surface, at constant field and temperature (viscosity conditions), from different non-equilibrium starting configurations. Under viscosity conditions, the system will converge to either the (single) local free energy minimum, if it exists, or the global minimum. The Gibbs free energy will decrease monotonically along a convergence path, but the model shows clearly that there is no such constraint on the internal energy, the entropy or the magnetization.