This work investigates the mixed site occupancy of aluminium and silicon atoms in the β-AlFeSi phase. For this purpose, the six mixed Al/Si sites of the β-AlFeSi structure were considered independent and alternatively substituted by Al or Si, thus generating 64 ordered structures or end-members. The enthalpy of formation of each end-member was calculated by DFT. These calculations allowed us to derive the enthalpy of mixing of the solid solution at 0 K, over a wide range of chemical compositions, from the Al-Fe binary system to the Si-Fe binary system. In addition, the heat capacities of the solid solution were determined using a Debye model based on the calculation of the elastic constants and the equation of state of each end-member. These heat capacity values were used along with the enthalpy of formation we calculated to determine the Gibbs free energies of all the end-members of the β-AlFeSi structure. Finally, the configurational entropy of mixing from the Compound Energy Formalism (CEF) for the configurational entropy of mixing was subsequently used to calculate the occupation fractions of the Si sites on the Al sites of the β-AlFeSi structure, at 300 K and 938 K, the latter being the thermal decomposition temperature of this compound. These original site occupancy data were used to quantify the chemical ordering of the solid solution and to compare different sublattice (SL) models. We thus highlight that the SL model of the β-AlFeSi solution most commonly accepted in the literature generates considerable errors in its thermodynamic description, contrary to the model proposed in this paper, which is both simple and particularly accurate, consisting in merging the sites Al(1)-Al(6), the sites Al(2)-Al(3) as well as the sites Al(4)-Al(5).