Modeling of Steam Methane Reforming (SMR) Over A Ni - AL O Catalyst
Modeling of Steam Methane Reforming (SMR) Over A Ni - AL O Catalyst
ISSN: 2458-9403
                                                                                                            Vol. 3 Issue 8, August - 2016
 Abstract—In present paper modeling of natural gas                 After desulfurization, natural gas and carbon dioxide
 reforming by steam over nickel-alumina catalyst is                mixed with water and flue gases and pre-heated to the
 investigated and key parameters such as pressure,                 reaction temperature. Inside the reactor, there are long
 temperature and mole ratio have been studied. The                 length and small diameter pipes with nickel-alumina
 results show that the reforming performance is                    catalyst contained. Reforming conditions with a high
 affected by temperature, pressure and ratio of                    temperature furnace are converted all hydrocarbons
 steam to methane. The equations for mass and                      into products. In this process, the required carbon
 energy balance equations and catalyst are                         dioxide usually achieved from returns all the carbon
 mentioned.                                                        dioxide and also a significant amount of outside unit [2].
 The combination of carbon monoxide and                            Process schematic is as follows:
 hydrogen are known as syngas and are used as a
 feedstock in the production of various products
 such as ammonia and liquid fuel applications.
 There are different techniques for the production of
 synthesis gas. One of these techniques is dry
 reforming that combination carbon dioxide and
 methane are converted into synthesis gas.
   The synthesis gas is odorless and colorless                                   Figure 1- Process schematic [3]
   poisonous gas that is a mixture of hydrogen and
   carbon monoxide and use in various industries                   Steam reforming technology was used from 1920
   including petrochemical and metallurgy industries.              onwards, but from the 1950s onwards was used as a
   There are several methods for the production of                 process for synthesis gas production of gas
   synthesis gas from steam reforming of natural gas is            compounds. The key point of these reactions is
   mostly done with hydrocarbons such as naphtha.                  endothermic that require high temperatures and
   Natural gas that primarily containing methane, reacts           decomposition reaction of H-C band is the first
   with water within tubular reactors at high temperature          controller. As shown in the following equations
   and absorb large amounts of heat, in the presence of            C2+reactions (reaction I) are irreversible, but the
   nickel-alumina catalyst and reaction final result is the        decomposition of methane (reaction II) equilibrium
   creation of reaction products and the ratio of                  relationship formed by the reaction of methane is
   hydrogen to carbon monoxide is as the most                      returned to the reaction [4].
   important criteria for reactor design process[1].
   Modeling gas production reactor unit in the                     I: CnH2n+2 + nH2O ↔ nCO + (2n+1) H2
   petrochemical industry is great importance.
                                                                   II: CH4 + H2O ↔ CO + 3H2
   II.    Process described
                                                                   Δ HII = +206 kJ/mol
                                                      www.jmest.org
 JMESTN42351725                                                                                                                    5394
                                                              Journal of Multidisciplinary Engineering Science and Technology (JMEST)
                                                                                                                        ISSN: 2458-9403
                                                                                                            Vol. 3 Issue 8, August - 2016
  Third reaction:                                                 𝑈
                                                                                             1          1
  CH4 + 2H2O ↔ CO2 + 4H2                                                𝐾𝑔         𝐷𝑝 𝐺 3 𝐶𝑝 𝜇 3
                                                                  = 0.4    [2.58 (     ) (    )
                                                                        𝐷𝑝          𝜇      𝐾𝑔
            𝑘3 𝑃𝐶𝐻4 𝑃𝐻22𝑂 𝑃𝐶𝑂2 𝑃𝐻0.5
                                  2
  𝑟3 = −        (        −           ) [7]                                                       0.4
           𝑑𝑒𝑛2   𝑃𝐻3.5      𝐾 3
                                                                              𝐷𝑝 𝐺 0.8 𝐶𝑝 𝜇
                     2                                            + 0.094 (       ) (       ) ]                    [24]
                                                                               𝜇       𝐾𝑔
                            −29336
  𝑘3 = 2.290 ∗ 1016 exp (         )        [8]
                              𝑇                                   Arrhenius equation:
  K3= K1 * K2                                                                  𝐸
                                                                  𝑘𝑖 = 𝐴𝑖 𝑒 −𝑅𝑇        , i = 1, 2, 3,…              [25]
   k3= kg.h/ Kmol. K Pa0.5
                                                                  Vant Hoff equation:
  K3= Pa2
                                                                               𝐸
                                                                  𝑘𝑗 = 𝐴𝑗 𝑒 −𝑅𝑇 , j = H2, CO , CO2 ,CH4               [26]
  𝑑𝑒𝑛 = 1 + 𝐾𝐶𝑂 𝑃𝐶𝑂 + 𝐾𝐻2 𝑃𝐻2 + 𝐾𝐶𝐻4 𝑃𝐶𝐻4
                          𝑃                                       VII.       The pressure drop in the reactor
                 + 𝐾𝐻2𝑂 𝐻2 𝑂⁄𝑃          [9]
                                𝐻2
                                                                  Because the reaction in a catalytic bed, the gas
                    −6
                               4604.28
  𝐾𝐶𝐻4 = 6.65 ∗ 10       exp (        )          [10]             pressure in the bed are facing. Pressure drop into bed
                                  𝑇                               it has an impact on the performance and efficiency [10].
                            −10666.356
  𝐾𝐻2𝑂 = 1.77 ∗ 103 exp (             )          [11]             𝛿𝑃𝑡           𝑒𝑔 𝑣𝑠2
                                𝑇                                     = −10−5 𝑓                                    [27]
                                                                  𝛿𝑧             𝑑𝑃
                             8497.71
  𝐾𝐶𝑂 = 8.23 ∗ 10−7 exp (           )            [12]             f from the following equation is obtained and also in
                                𝑇
                                                                  terms of atmospheric pressure:
                             9971.13
  𝐾𝐻2 = 6.12 ∗ 10−11 exp (          )            [13]                     1 − 𝜀𝐵     𝑏(1 − 𝜀𝐵 )
                                𝑇
                                                                  𝑓=             [𝑎+           ]                    [28]
                                                                            𝜀𝐵          𝑅𝑒
    V.      Mass balance
                                                                                               𝑅𝑒
  𝑑𝑋𝐶𝑂2    𝜋 𝑑𝑖2 𝜌𝑏 𝑟𝐶𝑂2                                          𝑎 = 1.75 , 𝑏 = 150 ,              < 500 [29]
        = (     )                     [14]                                                   1 − 𝜀𝐵
   𝑑𝑙       4       𝐹
                                                                                                         𝑅𝑒
   𝑋𝐶𝑂2 = 0 @ 𝑙 = 0                  [15]                         𝑎 = 1.24 , 𝑏 = 368 , 1000 <                 < 5000 [30]
                                                                                                       1 − 𝜀𝐵
  𝑟𝐶𝐻4 = 𝑟1 + 𝑟3                     [18]                         For industrial catalysts that are large in size due to the
                                                                  limited penetration velocity equation must be corrected
  𝑟𝐶𝑂2 = 𝑟2 + 𝑟3                    [19]                          that actual velocity of the particles obtained by
                                                                  multiplying the impact factor [11]. The impact factor is
  𝑑𝑋𝐶𝑂    𝜋 𝑑𝑖2 𝜌𝑏 𝑟𝐶𝑂                                            ratio of velocity in catalyst to velocity on the outer
       = (     )                      [20]
   𝑑𝑙      4      𝐹                                               surface of the catalyst.
                                                        www.jmest.org
  JMESTN42351725                                                                                                                  5396
                                                                   Journal of Multidisciplinary Engineering Science and Technology (JMEST)
                                                                                                                             ISSN: 2458-9403
                                                                                                                 Vol. 3 Issue 8, August - 2016
                                                                        λ𝐵 = (1 + 0.01P) ∑ 𝑦𝑖 𝜆𝑖                      [43]
  Ps,H2O–PH2O=(De,CO2/De,H2O)(PCO2-Ps,CO2)–(De,CH4/
                                                                                                  𝑖
  De,H2O) (Ps,CH4-PCH4)                  [34]                           𝜆𝑖 = 𝐴 + 𝐵𝑇 + 𝐶𝑇 2 + 𝐷𝑇 3                     [44]
                                                                     λB = Thermal conductivity of an ideal gas mixture
  Ps,H2 – PH2=(De,CO2/De,H2)(PCO2 – Ps,CO2) – (3De,CH4 /De,H2
                                                                         XI.      The viscosity of the ideal gas mixture
  ) (Ps,CH4 - PCH4)               [35]
      IX.Equations to calculate the physical properties                 The viscosity of the mixture can be obtained from the
         of the fluid:                                                  following equation:
  Assuming the gas will be ideal for mixed density:                     𝜇𝐵 = ∑𝑁 1 𝑦𝑖 𝜇𝑖                    [45]
  𝑃𝑉 = 𝑛𝑅𝑇                            [36]                              Viscosity is a function of temperature and pressure, but
      𝑚                                                                 its dependence to the temperature is much more than
  𝑛=                                  [37]                              pressure and only in very high pressures change
      𝑀
                                                                        noticeable [13].
        𝑅
  𝑃 = 𝑐𝑅 𝑇                           [38]                               From the following equation can be found on the
        𝑀
                                                                        temperature dependence of viscosity:
             𝑁
  𝑀 = ∑ 𝑦𝑖 𝑀𝑖                                   [39]
                                                                        𝜇𝑖 × 𝐴𝐽 = [0.807 𝑇𝑟𝑖0.618 − 0.357 exp(−0.449𝑇𝑟𝑖 ) +
            𝑖=1
                                                                        0.34 exp(−4.035𝑇𝑟𝑖 ) + 0.18]𝐹𝑄∗ 𝐹𝑃∗           [46]
  For an ideal gas specific heat the mixture is calculated                                            0.166
                                                                                         𝑇𝑐𝑖
  as follows:                                                           𝐴𝐽 = 0.176 {             }                  [47]
                  𝑛
                                                                                        𝑀𝑖3 𝑃𝑐𝑖4
      𝑖𝑔               𝑖𝑔
  𝐶𝑝𝑅 = ∑ 𝑦𝑖 𝐶𝑝𝑖                           [40]                         𝐹𝑃∗ = 1        0 < 𝑑𝑟 < 0.022                [48]
              𝑖=1
   𝑖𝑔
                                                                                        𝑑𝑖2 𝑃𝑐𝑖
  𝐶                                                                     𝑑𝑟𝑖 = 52.46                                [49]
   𝑝ه
           = 𝐴 + 𝐵𝑇 + 𝐶𝑇 2 + 𝐷𝑇 −2         [41]                                          𝑇𝑐𝑖2
  𝑅𝑇
                                                            www.jmest.org
  JMESTN42351725                                                                                                                        5397
                                                                                       Journal of Multidisciplinary Engineering Science and Technology (JMEST)
                                                                                                                                                 ISSN: 2458-9403
                                                                                                                                     Vol. 3 Issue 8, August - 2016
 entrance of the tube and produce methane.                                                  Because water gas shift reaction is exothermic
                                                                                            thereupon in a high-temperature reaction suppresses
                                                                                            so increase in temperature causes more methane into
                                                                                            H2 and CO, CO2 and convert CO into CO2 is less.
                                                                                  www.jmest.org
 JMESTN42351725                                                                                                                                                 5398
                                                                  Journal of Multidisciplinary Engineering Science and Technology (JMEST)
                                                                                                                            ISSN: 2458-9403
                                                                                                                Vol. 3 Issue 8, August - 2016
                                                                           PH2( bar)
                                                                                       7,9
feed advances gas-water shift reaction and reduces in
                                                                                       7,8
CO/H2 production. By increasing the ratio of H2O/CH4 in
input, the ratio of CO2/CH4 constant in input and                                      7,7
temperature is constant and increases the ratio of CO/H2                               7,6
production. By increasing of CO2/CH4 input in H2O/CH4                                  7,5
constant input and constant-temperature the CO/H2                                      7,4
ratio will reduce.                                                                                                   0             10           20              30
                                                                                                                               (E=6m) Length of catalyst
      Table 6 -values of E and H vant Hoff and Arrhenius               Figure 3 – H2 pressure on the length of the catalyst for
                      equation( kJ / mol)                                                    the (E= 6m)
    E1     E2    E3     ∆HCO     ∆HH2      ∆HCH4    ∆HH2O
   241    66     23      -70      -81       -38       89               Becomes smaller catalyst particles in the pipes are
 XIV.     Conclusions                                                  reduced catalyst loading and increases the pressure
 In the process, concentration gradient of gas                         drop across the bed.
 component is limited to a very thin layer of the active               In industry has always faced pressure restriction so can
 catalyst. Catalyst component concentration gradient is                built a catalyst with complex shapes in order to create
 very steep in grain and mass transfer resistance                      more internal level, low pressure drop, or active
 components in the catalyst grain are process controller.              substances in the outer layer of the catalyst, on the one
                                                                       hand improve the impact and also saving the use of
                                                                       active metal that is usually expensive.
                                                                                                           2
                                                                          PCO2( bar)
1,75
                                                                                        1,5
                                                                                                                      0              10            20           30
                                                                                                                               (E=6m) Length of catalyst
  Figure 1 - CO pressure on the length of the catalyst for the
                                                                                        Figure 4 – CO2 pressure on the length of the catalyst
                          (E=6m)
                                                                                                          for the (E=6m)
 Figure 2.Experimental data and model prediction
 against the input oxygen flow for changes of methane                                                                  3
                                                                                        Ration of water to methane
                                                           www.jmest.org
 JMESTN42351725                                                                                                                                                       5399
                                                                                                                                    Journal of Multidisciplinary Engineering Science and Technology (JMEST)
                                                                                                                                                                                              ISSN: 2458-9403
                                                                                                                                                                                  Vol. 3 Issue 8, August - 2016
80 0,015
                                                                                                                                           methan concentration
              The conversion of methane
                                                     70
                                                                                                                                                                              0,01
                                                     60
                                                     50
                                                                                                                                                                     0,005
                                                     40
                                                     30
                                                                                                                                                                                      0
                                                                        0        5         10         15         20       25
                                                                                                                                                                       2           4      06
                                                                                       input flow oxygen                                                             length of reactor
                                                                                                                                                    Figure10. Methane concentration changes in the
                                                                    Figure 6 - Oxygen flow at input to the conversion of                                            reactor length
                                                                                         methane
                                                                                                                                                                                       0,001
                                                                                                                                                                  CO2 concentration
                                                                    3                                                                                                                 0,0008
   Ratio of carbon dioxide to methane
                                          2,5                                                                                                                                         0,0006
                                                                    2                                                                                                                 0,0004
                                          1,5                                                                                                                                         0,0002
                                                                    1                                                                                                                         0
                                                                                                                                                                                                  0            2          4
                                          0,5                                                                                                                                                         length of reactor
                                                                        0            0,5         1           1,5              2
                                                                         The ratio of hydrogen to carbon monoxide in output
                                                                                                                                          Figure9. CO2 concentration changes in the reactor
                                                                        Figure 7 - The ratio of in output of carbon dioxide                                    length
                                                                                            to methane                                   XV.    Reference
                                                                                                                                         [1].Harrison. M.R,shires,T.M. wessels, jk,"methane
                                  Table 7. Specifications Feed and conditions                                                            emissions from the natural Gas industry" (1997)
                                                                                                                                         [2]. padro. CEG and v.putche," survey of the economics
   Feed flow rate (kmol/hr)                                                                                      3392                    of hydrogen technologies".national renewable energy
                                                                                                                                         laboratory, (September 1999)
   Temperature ( K)                                                                                               586                    [3].Hart stein, Arthur “Hydrogen production from
   Pressure ( KPa)                                                                                               2074                    Natural Gas".( 2003(
                                                                                                                                         [4].m.karimi–Golpayegani,                   A.Akhavan–
                                                                                                                                         Abdolahian,N.Mostoufi" simulation of a fluidized–bed
In this collection of methane in the reactor pressure                                                                                    steam reformer" university of tehran,iran( 2005)
decreased with increasing pressure. This process was                                                                                     [5].Trimm, D.L. and Wainwright, M.S. Catal. Today, 6,
assessed as a thermodynamically unfavorable reaction                                                                                     261, 1990
system.                                                                                                                                  [6]. John, N. A., Applied Catalysis A: General 176, 159,
                                                                                                                                         1999.
                                                                                                                                         [7]. Solymosi, F., Tolmacsov, P., K. Kedves, J. Catal.
                                                                    70                                                                   216, 377, 2003.
                                        The conversion of methane
                                                                                                                               www.jmest.org
JMESTN42351725                                                                                                                                                                                                            5400
                                                            Journal of Multidisciplinary Engineering Science and Technology (JMEST)
                                                                                                                      ISSN: 2458-9403
                                                                                                          Vol. 3 Issue 8, August - 2016
hydrotalcite-like anionic clay”, J. Catal., Vol. 221, pp.        [15].Jack  H.Lunferd,"Catalytic    Conversion    of
43-54, 2004.                                                     Methane to more useful chemical s and fuel s: a
[14]. Hoang D.L., Chan S.H., Ding O.L., Kinetic and              challenge   for  the    21st    century",  Catalysi
modeling Study of Methane Steam Reforming over                   sToday,63(2000)165-174.
Sulfide Nickel Catalyst on a Gamma Alumina Support,
Chemical Engineering Journal.
                                                     www.jmest.org
JMESTN42351725                                                                                                                   5401