Modeling the Drying Process of Porous Catalysts: Impact of the Pore Size Distribution.

The distribution of catalytically active species in heterogeneous porous catalysts strongly influences their performance and durability in industrial reactors. A drying model for investigating this redistribution was developed and implemented using the finite volume method. This model embeds an analytical approach regarding the permeability and capillary pressure from arbitrary pore size distributions. Subsequently, a set of varying pore size distributions are investigated, and their impact on the species redistribution during drying is quantified. It was found that small amounts of large pores speed up the drying process and reduce internal pressure build up significantly while having a negligible impact on the final distribution of the catalytically active species. By further increasing the amount of large pores, the accumulation of species at the drying surface is facilitated.

CFD-DEM Evaluation of the Clustering Behavior in a Riser-the Effect of the Drag Force Model.

Riser reactors are frequently applied in catalytic processes involving rapid catalyst deactivation. Typically heterogeneous flow structures prevail because of the clustering of particles, which impacts the quality of the gas-solid contact. This phenomenon results as a competition between fluid-particle interaction (i.e., drag) and particle-particle interaction (i.e., collisions). In this study, five drag force correlations were used in a combined computational fluid dynamics-discrete element method Immersed Boundary Model to predict the clustering. The simulation results were compared with experimental data obtained from a pseudo-2D riser in the fast fluidization regime. The clusters were detected on the basis of a core-wake approach using constant thresholds. Although good predictions for the global (solids volume fraction and mass flux) variables and cluster (spatial distribution, size, and number of clusters) variables were obtained with two of the approaches in most of the simulations, all the correlations show significant deviations in the onset of a pneumatic transport regime. However, the correlations of Felice (Int. J. Multiphase Flow1994, 20, 153-159) and Tang et al. [AIChE J.2015, 61 ( (2), ), 688-698] show the closest correspondence for the time-averaged quantities and the clustering behavior in the fast fluidization regime.

Direct Numerical Simulation of Reactive Fluid-Particle Systems Using an Immersed Boundary Method.

In this paper, direct numerical simulation (DNS) is performed to study coupled heat and mass-transfer problems in fluid-particle systems. On the particles, an exothermic surface reaction takes place. The heat and mass transport is coupled through the particle temperature, which offers a dynamic boundary condition for the thermal energy equation of the fluid phase. Following the case of the unsteady mass and heat diffusion in a large pool of static fluid, we consider a stationary spherical particle under forced convection. In both cases, the particle temperatures obtained from DNS show excellent agreement with established solutions. After that, we investigate the three-bead reactor, and finally a dense particle array composed of hundreds of particles distributed in a random fashion is studied. The concentration and temperature profiles are compared with a one-dimensional heterogeneous reactor model, and the heterogeneity inside the array is discussed.

Deactivation Kinetics of Solid Acid Catalyst with Laterally Interacting Protons.

Differences in catalyst deactivation kinetics in solid acid catalysis are studied with catalyst models that allow for lateral interaction between protons. Deactivation of a solid acid catalyst with laterally interacting protons induces inhomogeneity of proton reactivity that develops with time. As a consequence, product selectivity changes and deactivation will accelerate. This is demonstrated by simulations of the deactivation kinetics of the alkylation reaction of propylene with isobutane. The effect of lateral interactions between protons arises because initial catalyst deactivation is not caused by pore blocking or coke deposition but by a molecular mechanism where protons are consumed due to the formation of stable nonreactive carbenium ions. High selectivity to alkylate requires a catalyst with protons of high reactivity. When protons become consumed by formation of stable deactivating carbenium ions, initially reactive protons are converted into protons of lower reactivity. The latter only catalyze deactivating oligomerization reactions. Simulations that compare the deactivation kinetics of a catalyst model with laterally interacting protons and a catalyst model that contains protons of similar but different reactivity, but that do not laterally interact, are compared. These simulations demonstrate that the lateral interaction catalyst model is initially more selective but also has a lower stability. Catalyst deactivation of the alkylation reaction occurs through two reaction channels. One reaction channel is due to oligomerization of reactant propylene. The other deactivation reaction channel is initiated by deprotonation of intermediate carbenium ions that increase alkene concentration. By consecutive reactions, this also leads to deactivation. The hydride transfer reaction competes with oligomerization reactions. It is favored by strongly acid sites that also suppress the deprotonation reaction. Within the laterally interacting proton catalyst model, when reactive protons become deactivated, weakly reactive protons are generated that only catalyze the deactivating alkene oligomerization and consecutive reactions. This rapid formation of the weakly reactive protons is the cause of decreasing selectivity with reaction time and increased rate of deactivation. Solutions of the mean field kinetic equations as well as stochastic simulations are presented. Comparative simulations with a reduced number of neighbors of the protons illustrate decreased deactivation rates when the proton density decreases. Island formation of adsorbed reaction intermediates on the catalyst surface is observed in stochastic kinetics simulations. When alkylation selectivity is high, this island formation increases the rate of catalyst deactivation in comparison to the rate of deactivation according to the mean field studies. A nonlinear dynamics model of proton dynamics is provided, which shows that the differences between stochastic and mean field simulations are due to frustrated proton state percolation.

Direct Numerical Simulation of Fluid Flow and Mass Transfer in Particle Clusters.

In this paper, an efficient ghost-cell based immersed boundary method is applied to perform direct numerical simulation (DNS) of mass transfer problems in particle clusters. To be specific, a nine-sphere cuboid cluster and a random-generated spherical cluster consisting of 100 spheres are studied. In both cases, the cluster is composed of active catalysts and inert particles, and the mutual influence of particles on their mass transfer performance is studied. To simulate active catalysts the Dirichlet boundary condition is imposed at the external surface of spheres, while the zero-flux Neumann boundary condition is applied for inert particles. Through our studies, clustering is found to have negative influence on the mass transfer performance, which can be then improved by dilution with inert particles and higher Reynolds numbers. The distribution of active/inert particles may lead to large variations of the cluster mass transfer performance, and individual particle deep inside the cluster may possess a high Sherwood number.

Effect of Superficial Gas Velocity on the Solid Temperature Distribution in Gas Fluidized Beds with Heat Production.

The hydrodynamics and heat transfer of cylindrical gas-solid fluidized beds for polyolefin production was investigated with the two-fluid model (TFM) based on the kinetic theory of granular flow (KTGF). It was found that the fluidized bed becomes more isothermal with increasing superficial gas velocity. This is mainly due to the increase of solids circulation and improvement in gas solid contact. It was also found that the average Nusselt number weakly depends on the gas velocity. The TFM results were qualitatively compared with simulation results of computational fluid dynamics combined with the discrete element model (CFD-DEM). The TFM results were in very good agreement with the CFD-DEM outcomes, so the TFM can be a reliable source for further investigations of fluidized beds especially large lab-scale reactors.