In the present study, a comprehensive mathematical model is developed to simulate the dynamic behaviour of an industrial slurry-phase olefin catalytic polymerization loop reactor series. More specifically, the effects of various operating conditions on the dynamic reactor behaviour (i.e., reaction temperature and pressure, inflow rates of catalyst, monomers and diluent, etc.) as well as the on the molecular and rheological polyolefin properties (i.e., M-n, M-w, MWD, complex viscosity, etc.) are fully assessed. According to the proposed modeling approach, each loop reactor (i.e., consisting of the loop reactor and the settling legs) is modeled as an ideal CSTR in series with a semi-continuous product removal unit. Dynamic macroscopic mass species and energy balances are derived to calculate the dynamic evolution of the concentrations of the various molecular species as well as of temperature profiles and heat removal in the two loop reactors. The polymer molecular properties (i.e., number- and weight-average molecular weights and molecular weight distribution) are determined by employing a generalized multi-site, Ziegler-Natta (Z-N) kinetic scheme in conjunction with the well-known method of moments. All the thermodynamic calculations, regarding the equilibrium species concentration in the various phases (i.e., solids and liquid), are carried out using the Sanchez-Lacombe Equation of State (S-L EOS). It is shown that the proposed comprehensive model is capable of simulating the dynamic operation of an industrial slurry-phase cascade-loop reactor series under different plant operating policies (i.e., start-up, grade transition, etc.). (C) 2010 Elsevier Ltd. All rights reserved.