Non-isothermal flow is simulated in the screw and kneading disc elements in a modular co-rotating twin screw extruder. Methods of calculating mean temperature rises for individual elements are discussed and results are presented. The implications of non-isothermal operation for scale-up is discussed. A method is then described for making calculations in a modular co-rotating machine, that contains many different elements. Example calculations are given showing the tendency of larger machines to buildup greater temperatures when viscous heating is included.