Implementation of Reaction Kinetics into Reactor Model Followed by Validation of the Reactor Model

Sammanfattning: This master thesis project was performed in collaboration with Johnson Matthey. Johnson Matthey is a worldwide company working with sustainable technologies. The site in Perstorp, where this work was performed, focus on formaldehyde production. In order to continuously improve chemical processes, simulation and modeling are advantageous tools. They can give a better understanding of the processes and the optimum operation conditions. By using reactor modeling, design improvements can be investigated, and the efficiency of the plant can be increased. The aim with this project was to implement a kinetic model into an existing reactor model and to further on validate the reactor model against data from the pilot reactor. The model was developed in Aspen Custom Modeler. The kinetic model was described by the power law, but at an early stage, numerical errors were obtained with this kinetic model. The power law equations were replaced with another kinetics model from literature. The kinetic model from literature was more complex and more equations were added, including the surface fraction of the components. The development of the kinetic model was an ongoing project and two kinetic models were developed from the literature kinetic model. The rebased model including a reference temperature and the second one with refitted parameters to the experimental data provided from Johnson Matthey. The rebased model was validated against the pilot reactor and the outlet composition and temperature were compared. The results showed deviation when compared with the pilot plant data, with relative deviation for the temperature profile of +15% and relative deviation for the outlet composition of +2.6% for formaldehyde. A possible explanation for the deviation could be higher reaction rates for reaction 1 and 2. However, the reactor model showed similar behaviour of the pilot plant data regarding the hot spots and similar range of concentration for the main components in the outlet, which shows great potential describing the pilot plant in the future after final validation. The different layers of catalyst along the tube were implemented in the reactor model and the inert dilution was taking into consideration. When analysing the result in the radial direction, the reactor model showed expected results for the composition and temperature, i.e. constant values for the composition in the radial direction and warmer temperature in the center of the tube in comparison to the wall. This result indicates that the kinetic model was successfully implemented in 2D. To improve the reactor model, the refitted kinetic model can be evaluated further. This kinetic model with refitted parameters showed problems with convergence and due to the time limit of this project, the convergence problem could not be solved. The mass transfer limitation within the catalyst particle should preferably be implemented as well in order to describe the correct behaviour in the reactor tube.