The present paper reports the results of an extensive numerical analysis of electromagnetic fields in random dielectric materials. The effective permittivity of a two-dimensional (2-D) dielectric mixture is calculated by FDTD simulations of such a sample in a TEM waveguide, Various theoretical bounds are tested in light of the numerical simulations. The results show how the effective permittivity of a mixture with random inclusion positionings is distributed, All possible permittivity values lie between Wiener limits, and according to FDTD simulations the values are almost always between Hashin-Shtrikman limits, Calculated permittivity distribution is also compared with theoretical mixture models. No model seems to be able to predict the simulated behavior over the whole range of volume fraction.