Simul 5, the simulation program based on the 1-D model of electrophoresis has been extended to simulate isoelectric trapping (IET) separations that take place in recirculating multicompartmental electrolyzers (MCEs). In the extended Simul 5, the simulated separation space between the anode and cathode can be divided into multiple segments to represent the anode compartment, separation compartment(s) and the cathode compartment. The compartments may have identical or different cross sections. A new algorithm simulates convective mixing that occurs in the recirculating MCEs where the distances between the buffering membranes are short and the velocities of tangential flows through the compartments, orthogonal to the electric field, are high. The intensity of simulated mixing can be independently controlled in each compartment. pH transients that were reported to occur during the desalting step in IET separations were simulated with the extended Simul 5 program: the main features of the experimental results were reproduced by the simulations. The simulations shed light on the possible causes of uneven anion and cation removal rates, pH transients and the transient invasion of the separation compartment by components of the electrode solutions that might occur during the desalting step.