457943 Multiscale Simulation of Heterogeneous Catalysis: Butadiene Hydrogenation from Molecular Reaction Steps to the Fluid Flow

Sunday, November 13, 2016: 4:45 PM
Franciscan D (Hilton San Francisco Union Square)
Phillippe Sautet1, Ludovic Briquet2, Céline Chizallet3, Pascal Raybaud4, Jean-Marc Schweitzer5, Jan J. Verstraete6, Anthony Wachs7, Matthieu Rolland4 and Abdelkader Hammouti4, (1)Chemical and Biomolecular Engineering, University of California Los Angeles, Los Angeles, CA, (2)ENS Lyon, Université de Lyon, Lyon, France, (3)Catalysis, IFP Energies Nouvelles, Solaize, France, (4)IFP Energies Nouvelles, Solaize, France, (5)Catalysis and Separation Division, IFP Energies nouvelles, 69360 Solaize, France, (6)Process Design and Modeling Division, IFP Energies nouvelles, 69360 Solaize, France, (7)University of British Columbia, Vancouver, Canada

The selective hydrogenation of poly-unsaturated organic molecules is of high industrial interest in the area of fine and polymer chemistry. In order to have a good control of the reaction, it is necessary to gain a detailed understanding of all its stages and scales. On the smallest atomic scale, the reaction is governed by the interactions of the olefin with the metal catalyst surface. At a larger scale, the reaction is governed by temperature and pressure conditions and reactant and product flow in the reactor. Nevertheless the two scales are closely interlinked since local conditions at the surface of the catalyst's particle have a direct impact on the atomistic environment of the metal surface. In particular, hydrogen is known to easily cover the palladium and platinum surfaces that are used in such reactions.1The coverage may depend on the local temperature and pressure condition but also on the kinetics of the reaction in those conditions. In this work we have studied the reaction with a multi-scale approach, the atomic scale with electronic quantum chemical methods, and the macroscopic scale with Computational Fluid Dynamics (CFD) simulations, both scales being coupled by a mesoscopic stage using mean-field kinetic simulations.

Using periodic plane wave density functional theory calculations,2 the hydrogenation reactivity of the simplest poly-unsaturated hydrocarbon, butadiene, on the Pd(111) surface has been investigated. The effect of the hydrogen coverage on the Pd(111) surface has been accounted for in order to model different macroscopic reaction conditions. The most important energy barriers have been exactly determined using the Nudge Elastic Band method3 while transition states suspected to have a higher potential energy have been estimated using a Bronsted-Evan-Polanyi relationship.4Using transition state theory, the hydrogenation energy profile at low hydrogen coverage, developed in an earlier stage of the project, has been used to build a first implementation of a micro-kinetic model. The PeliGRIFF CFD code developed at IFP Energies Nouvelles has been upgraded in order to include for the first time a reaction term in the catalyst grain. The obtained kinetic law has been included in PeliGRIFF and it was shown that with intra-grain diffusion coefficient one order of magnitude lower than the gas phase diffusion, the reaction mainly happens in a shell surrounding the grain.


  1. C. Chizallet, G. Bonnard, E. Krebs, C. Thomazeau, P. Raybaud, J. Phys. Chem. B 115 (2011) 12135-12149.
  2. G. Kresse, J. Furthmüller, Phys. Rev. B 54 (1996) 11169-11186.
  3. D. Sheppard, P. Xiao, W. Chemelewski, D. D. Johnson, G. Henkelman, J. Chem. Phys. 136 (2012) 074103.
  4. D. Loffreda, F. Delbecq, F. Vigne, P. Sautet, Angew. Chem. Int. Ed. 48 (2009) 8978-8980.

Extended Abstract: File Not Uploaded
See more of this Session: Multiscale Modeling
See more of this Group/Topical: Catalysis and Reaction Engineering Division