An electrocoagulation (EC) model is developed for hexavalent chromium reduction and precipitation, using iron electrodes. Parallel removal mechanisms such as adsorption of chromium on ferrihydrite and direct reduction at the cathode is assumed negligible due to low concentration of Cr(VI). The reaction model presented for batch system represents species complexation, precipitation/dissolution, acid/base, and oxidation-reduction reactions. Batch reactor simulation is verified using experimental data obtained by Sarahney et al. (2012), where the effect of initial chromium concentration, pH, volumetric current density, and ionic strength is considered (Sarahney et al., 2012). The model couples multicomponent ionic transport in MATLAB with chemical reaction model in PHREEQC, as a widely used computational programming tool and a geochemical reaction simulator with comprehensive geochemistry databases. The suggested current density is 0.05−0.3mA/cm2 and the surface to volume ratio in batch reactor is considered 0.017 1/cm. Design parameters are presented for operation of a flow-through hexavalent chromium removal using electrocoagulation by iron electrode to treat Cr(VI) in range of 10–50 mg/L. The operational parameters for a flow-through EC reactor for Cr(VI) removal is suggested to follow 0.05mA/cm2≤3nFeQcCrVIinlet≤0.3mA/cm2.