In this paper, we present a fundamental study on the stationarity and ergodicity of eight classes of sum-of-cisoids (SOC) processes for the modeling and simulation of frequency-nonselective mobile Rayleigh fading channels. The purpose of this study is to determine which classes of SOC models enable the design of channel simulators that accurately reproduce the channel's statistical properties without demanding information on the time origin or the time-consuming computation of an ensemble average. We investigate the wide-sense stationarity, first-order stationarity of the envelope, mean ergodicity, and autocorrelation ergodicity of the underlying random processes characterizing the different classes of stochastic SOC simulators. The obtained results demonstrate that only the class of SOC models comprising cisoids with constant gains, constant frequencies, and random phases is defined by a set of stationary and ergodic random processes. The analysis presented here can easily be extended with respect to the modeling and simulation of frequency-selective single-input single-output (SISO) and multiple-input multiple-output channels. For the case of frequency-selective SISO channels, we investigate the stationarity and ergodicity in both time and frequency of 16 different classes of SOC simulation models. The findings presented in this paper can be used in the laboratory as guidelines to design efficient simulation platforms for the performance evaluation of modern mobile communication systems.