A Grand-canonical Monte-Carlo simulation method is investigated. Due to charge neutrality requirement of electrolyte solutions, ions must be added to or removed from the system in groups. It is then implemented to simulate solution of 1:1, 2:1 and 2:2 salts at different concentrations using the primitive ion model. | Size effect in Grand-canonical Monte-Carlo simulation of solutions of electrolyte