Almost all catalytic reactions in the petrochemical industry are reversible and therefore their conversion is limited by the thermodynamic equilibrium. This conservative limitation can be broken by using selective membranes to remove one of the products. In this paper this revolutionary concept is used for the dehydrogenation reaction where the selective membranes are used for the perm-selective removal of hydrogen.1-3 These membranes have 100% selectivity for the removal of hydrogen. Most efficient configuration is when in the other side of the membrane is a hydrogenation reaction and the flows in the two sides of the membrane are counter-current.
International Journal of Petrochemical Science & Engineering
Research ArticleOpen Access
Simulation investigation of novel membrane reactors for catalytic reactions
Governing equations for membrane reactorTo obtain the mole balance equation and the energy balance equation, a differential element inside the membrane reactor was considered. After writing the two balances under the steady state assumption, both sides of the resulting equations were divided by the thickness of the differential element, which was then forced to approach zero. The resulting balances equations of the shell side can be expressed as:Mole balance:Energy balance:( ) ( )Pressure drop:The corresponding balance equations for the tube side can be expressed as:Mole balance:Energy balance:The hydrogen flux across the membrane surface obeys Sievert's law,The pre-exponential constant, permeation activation energy, and the thickness of the hydrogen permeation membrane are taken as 6.33×10 , 15700 J/mole, and 1~2×10 -5 m, respectively. 8 Heat transfer across the membrane involves both convection from the gas mixture to the membrane, conduction across the membrane layer, and finally convection from the membrane to the second gas mixture. Radiation of heat is neglected. The membrane tube is considered to be The convective heat transfer coefficients in equation (24) are calculated using Leva's correlation (1949).11 For the shell side in which the reacting mixture is heated up, the convective heat transfer coefficient is calculated by: 11 0.9 6 0.813 exp…………………… (25) In contrast, the reacting mixture in the tube side is cooled and consequently the convective heat transfer coefficient is calculated by: 11 0.7 4.6 3.50 expPhysical properties such as, thermal conductivity, gas density and viscosity, and heat capacities are taken as functions of temperature from Yaws.
Boundary conditionsIn the case of the co-current operation ( ) 2 b= , the above system of differential equations gives an initial value problem which can be solved by a Runge-Kutta Verner fifth and sixth order method with an automatic step size, double precision calculation, and a relative error of 1×10 -12 to ensure high accuracy. The initial conditions are:For the counter-current operation case ( ) 1 b= , the above system of differential equations results in a split two-point boundary value problem which can be solved by an orthogonal c...