The aim of the present paper is to develop a reliable and accurate model of the wastewater biochemical treatment process and to explore the behaviour through a general dynamic simulation environment, namely the INtegrated Simulation Environment Language (INSEL), for the analysis of the energy demand of the whole wastewater treatment plant. In particular, the presented model pays special attention to the chemical kinetics involved in the activated sludge process for the reduction of nitrogen and carbon compounds. According to the best practices, the plant configuration considered in this work includes the denitrification-nitrification process, performed by completely mixed reactors. In particular, the process analysed in this paper is based on the Ludzak-Ettinger process. The biological process is simulated according to the well-known method widely used in the literature, namely the Activated Sludge Model No 1 (ASM1). The model includes a set of equations for the calculation of aerobic growth of heterotrophs, anoxic growth of heterotrophs, aerobic growth of autotrophs, decay of autotrophs, ammonification of soluble nitrogen, hydrolysis of entrapped organics, and hydrolysis of entrapped organic nitrogen. All these equations, along with energy and mass balances, are solved by the explicit Euler method. The developed model is validated using literature data, showing a great accuracy (deviation below 1%). As for the temperature, results show that, between 15 and 25 °C, in the initial part of the process, transport effects dominate the consumption ones. When the temperature is higher than 30 °C, nitrate consumption is so fast that biomass growth is limited by this effect. Conversely, in case of low temperatures (5–10 °C), biomass growth is not limited by nitrate availability. Finally, results also showed that temperature significantly affects the denitrification process, whereas the effect on the oxygen is lower.