We have developed a kinetics model that incorporates the essential mechanisms of the enzymatic depolymerization of cellulose but retains the computational speed necessary for process design and optimization. The model incorporates the distinct modes of action an endoglucanase (EGI) and an exoglucanase (cellobiohydrolase, CBHI) for depolymerization of insoluble cellulose polymer chains, and beta-glucosidase acting on soluble oligosaccharides. The processive action of CBHI is modeled using a population-balance approach to describe the depolymerization of cellulose chains, which incorporates distinct steps for adsorption, complexation, processive hydrolysis, and desorption. Random-chain scission of insoluble cellulose by EGI is also handled using population-balance equations, which allows for a kinetic description of the transformation of a polydisperse distribution of cellulose polymer lengths. Utilizing population balance models, which track the evolution of the distribution of polymer lengths, does not require solving equations for all chemical species present in the reacting mixture and affords computationally efficient simulations.