We developed a mathematical model to describe the solution polymerization of olefins with two single‐site catalysts in a series of two CSTRs. The model was used to simulate processes where semi‐crystalline macromonomers produced in the first reactor are incorporated as long chain branches onto amorphous (or lower crystallinity) chains in the second reactor (cross‐products). The simulation results show that CSTRs are more efficient to make chains with high LCB density and high weight percent of cross‐products. The model can also predict the polydispersity index, average chain lengths, and fractions of the different polymer populations, and help the polymer reactor engineer formulate new products with complex microstructures.