We consider the sample to sample fluctuations that occur in the value of a thermodynamic quantity
P in an ensemble of finite systems with quenched disorder, at equilibrium. The variance of
P,
VP, which characterizes these fluctuations is calculated as a function of the systems' linear size
l, focusing on the behavior at the critical point. The specific model considered is the bond-disordered Ashkin-Teller model on a square lattice [Phys. Rev. 64, 178 (1943)]. Using extensive Monte Carlo simulations, several bond-disordered Ashkin-Teller models were examined, including the bond-disordered Ising model and the bond-disordered four-state Potts model. It was found that far from criticality all thermodynamic quantities which were examined (energy, magnetization, specific heat, susceptibility) are strongly self-averaging, that is ${V}_{P}\ensuremath{\sim}{l}^{\ensuremath{-}d}$ (where
d=2 is the dimension). At criticality though, the results indicate that the magnetization
M and the susceptibility $\ensuremath{\chi}$ are non-self-averaging, i.e., $\frac{{V}_{\ensuremath{\chi}}}{{\ensuremath{\chi}}^{2}}, \frac{{V}_{M}}{{M}^{2}}\ensuremath{\nrightarrow}0$. The energy
E at criticality is clearly weakly self-averaging, that is ${V}_{E}\ensuremath{\sim}{l}^{\ensuremath{-}{y}_{v}}$ with
0<yv<d. Less conclusively, and possibly only as a transient behavior, the specific heat too is found to be weakly self-averaging. A phenomenological theory of finite size scaling for disordered systems is developed, based on physical considerations similar to those leading to the Harris criterion. Its main prediction is that when the specific heat exponent $\ensuremath{\alpha}<0$ ($\ensuremath{\alpha}$ of the disordered model) then, for a quantity
P which scales as ${l}^{\ensuremath{\rho}}$ at criticality, its variance
VP will scale asymptotically as ${l}^{2\ensuremath{\rho}+\frac{\ensuremath{\alpha}}{\ensuremath{\nu}}}$. The theory is not applicable in the asymptotic limit ($l\ensuremath{\rightarrow}\ensuremath{\infty}$) to the bond-disordered Ashkin-Teller model where $\frac{\ensuremath{\alpha}}{\ensuremath{\nu}}={0}_{+}$. Nonetheless in the accessible range of lattice sizes we found very good agreement between the theory and the data for ${V}_{\ensuremath{\chi}}$ and
VE. The theory may also be compatible with the data for the variance of the magnetization
VM and the variance of the specific heat
VC, but evidence for this is less convincing.
The relative variances of the magnetization Rm, heat capacity Rc, and susceptibility Rχ are calculated by the Monte Carlo method in the spin lattice four-component weakly diluted Potts model on a square lattice at a spin density p = 0.80. It is shown that introducing a disorder in the form of nonmagnetic impurities into the 2D Potts model leads to nonzero values of Rm, Rc, and Rχ at the critical point. It has been found that these values decrease noticeably for systems with linear sizes L ≥ 120 interatomic distances.