Particle size distributions of obtained samples from several sampling campaigns were determined and raw data were mass balanced before being used in simulation studies. After determination of breakage function, selection function, Bond work index, residence time distribution parameters, and Whiten’s model parameters for air separators and diaphragms between the two compartments of tube ball mills, performance of the circuits was simulated for given throughputs and feed particle size distribution. Whiten’s model parameters were determined by GA (genetic algorithm) toolbox of MATLAB software. Based on implemented models for modeling and simulation, optimization of circuits was carried out. It increased nearly 10.5% and 15.8% in fresh feed capacity input to each tube ball mill. In addition, circulating load ratios of circuits are modified to 118% and 127% from low level of 57% and 22%, respectively, and also cut points of air separators are adjusted at 30 and 40 μm from high range of 53 and 97 μm, respectively. All applications helped in well-operation and energy consumption reduction of equipments.