Ion distribution models for defect fluorite ZrO2 - AO1.5 (A = Ln, Y) solid solutions: II. Thermodynamics of mixing and ordering