modelSteamMethaneReformer_L4

Discretized pseudohomogeneous PFR model of a fixed-bed steam methane reformer

Extends from Base.PartialFixedbedReactorIdealGas_L4 (Discretized model of a fixed-bed reactor using ideal gas models).

Information

1. Purpose of model

This model represents a discretized fixed bed steam methane reformer with reaction kinetics with constant effectiveness factors.

2. Level of detail, physical effects considered, and physical insight

The reactor is discretized and in each volume mass, impulse and energy balances as well as heat transfer between furnace gas and synthesis gas in the tubes and reaction rate equations are solved. These equations are taken from Nandasana et al. [1] except that the pressure loss is assumed to be constant and the effective

reaction rates are

calculated using constant effectiveness factors. Also the mass balances are stationary so changes in density are neglected.

3. Limits of validity

The model is valid if the changes of the effectiveness factors and the pressure loss are negligible.

4. Interfaces

gasPortIn: ideal gas inlet

gasPortOut: ideal gas outlet

heat: heat port

5. Nomenclature

(no elements)

6. Governing Equations

The used equations are described in [1] except for the changes described in 2. The pressure calculation for each volume can be done either using the pressure in the middle or at the end of the volume.

7. Remarks for Usage

(no remarks)

8. Validation

The results of the model were compared to results from Rajesh et al. [2] and the outlet variables fit the results satisfyingly.

9. References

[1] Nandasana, Anjana D.; Ray, Ajay K.; Gupta, Santosh K. (2003): Dynamic Model of an Industrial Steam Reformer and Its Use for Multiobjective Optimization. In: Ind. Eng. Chem. Res. 42 (17), S. 4028–4042. DOI: 10.1021/ie0209576.

[2] Rajesh, J. K.; Gupta, Santosh K.; Rangaiah, G. P.; Ray, Ajay K. (2000): Multiobjective Optimization of Steam Reformer Performance Using Genetic Algorithm. In: Ind. Eng. Chem. Res. 39 (3), S. 706–717. DOI: 10.1021/ie9905409.

10. Version History


Model created by Carsten Bode (c.bode@tuhh.de) on Tue Apr 05 2016

Parameters

TypeNameDefaultDescription
IntegerN_comp (from PartialFixedbedReactorIdealGas_L4)medium.ncnumber of components
SI.CrossSectionA_c (from PartialFixedbedReactorIdealGas_L4)N_tube*pi/4*dia_tube_i^2total crosssectional area of all reactor tubes
SI.LengthDelta_x (from PartialFixedbedReactorIdealGas_L4)l/N_cvLength of one control volume
TransiEnt.Basics.Media.Gases.Gas_VDIWA_SG6_vargas_sg6Medium used in the steam methane reformer
SI.AreaA_tube_oN_tube*pi*dia_tube_o*lOuter surface of all tubes
Real[4]AK_j{6.65e-4, 8.23e-5, 6.12e-9, 1.77e5}Pre-exponential factors for adsorption energies {CH4, CO, H2, H2O}
SI.MolarEnthalpy[4]dH_j{-38280, -70650, -82900, 88680}Adsorption enthalpies for components j {CH4, CO, H2, H2O}
Real[N_reac]Ak_i1/3.6*{4.225e15, 1.955e6, 1.02e15}Pre-exponential factors for activation energies
Real[N_reac]AK_i{4.707e12, 1.142e-2, 5.375e10}Pre-exponential factors for enthalpy changes of reactions
Fundamental Definitions
TILMedia.GasTypes.BaseGasmedium (from PartialFixedbedReactorIdealGas_L4)simCenter.gasModel2Medium model
BooleanuseHomotopy (from PartialFixedbedReactorIdealGas_L4)simCenter.useHomotopytrue if homotopy should be used
IntegerN_cv (from PartialFixedbedReactorIdealGas_L4)1Number of control volumes
IntegerN_reac (from PartialFixedbedReactorIdealGas_L4)Number of reactions
Integer[N_reac,N_comp - 1]nu (from PartialFixedbedReactorIdealGas_L4)Matrix with stochiometric coefficients of all components-1 in all reactions
SI.MolarEnergy[N_reac]E_i (from PartialFixedbedReactorIdealGas_L4)Activation energies
SI.MolarEnthalpy[N_reac]dH_R_i (from PartialFixedbedReactorIdealGas_L4)Reaction enthalpies
SI.Efficiency[N_cv,N_reac]eff (from PartialFixedbedReactorIdealGas_L4)Effectiveness factors for the reactions
SI.PressureDifferenceDelta_p (from PartialFixedbedReactorIdealGas_L4)Total pressure loss over the reactor
IntegerpressureCalculation (from PartialFixedbedReactorIdealGas_L4)1Method of pressure calculation
Geometry
RealN_tube (from PartialFixedbedReactorIdealGas_L4)Number of tubes
SI.Diameterdia_tube_i (from PartialFixedbedReactorIdealGas_L4)Inner tube diameter
SI.Lengthl (from PartialFixedbedReactorIdealGas_L4)Length of reactor
SI.Diameterdia_tube_o0.1020Outer tube diameter
Catalyst
SI.Diameterdia_part (from PartialFixedbedReactorIdealGas_L4)Equivalent pellet diameter
SI.VolumeFractioneps_bed (from PartialFixedbedReactorIdealGas_L4)Bed porosity
SI.VolumeFractioneps_cat (from PartialFixedbedReactorIdealGas_L4)Catalyst porosity
SI.Densityd_cat (from PartialFixedbedReactorIdealGas_L4)Density of catalyst particle
SI.Densityd_bed (from PartialFixedbedReactorIdealGas_L4)d_cat*(1 - eps_bed)Density of fixed-bed
SI.SpecificHeatCapacitycp_cat (from PartialFixedbedReactorIdealGas_L4)Specific heat capacity of catalyst
Nominal Values
SI.Temperature[N_cv]T_nom (from PartialFixedbedReactorIdealGas_L4)Nominal gas and catalyst temperature in the control volumes
SI.Pressure[N_cv]p_nom (from PartialFixedbedReactorIdealGas_L4)Nominal pressure in the control volumes
SI.MassFraction[N_cv,N_comp - 1]xi_nom (from PartialFixedbedReactorIdealGas_L4)Nominal values for mass fractions
Heat Transfer
SI.AreaA_refractor1164*N_tube/176Surface area of refractor
SI.AreaA_flame0.01Surface area of a flame of one burner
SI.Emissivityem_gas0.1Emissivity of furnace gas
SI.Emissivityem_flame0.1Emissivity of flames
SI.Emissivityem_tube0.95Emissivity of reformer tubes
SI.TemperatureT_adi2200Adiabatic flame temperature
RealN_burner112*N_tube/176Number of burners

Connectors

TypeNameDefaultDescription
TransiEnt.Basics.Interfaces.Gas.IdealGasEnthPortIngasPortIn (from PartialFixedbedReactorIdealGas_L4)
TransiEnt.Basics.Interfaces.Gas.IdealGasEnthPortOutgasPortOut (from PartialFixedbedReactorIdealGas_L4)
Modelica.Thermal.HeatTransfer.Interfaces.HeatPort_a[N_cv]heat (from PartialFixedbedReactorIdealGas_L4)

Components

TypeNameDefaultDescription
TransiEnt.SimCentersimCenter (from PartialFixedbedReactorIdealGas_L4)
SI.Temperature[N_cv]T (from PartialFixedbedReactorIdealGas_L4)Gas and catalyst temperature in the control volumes
SI.Pressure[N_cv]p (from PartialFixedbedReactorIdealGas_L4)Pressure in the control volumes
SI.MassFraction[N_cv,N_comp - 1]xi (from PartialFixedbedReactorIdealGas_L4)Mass fraction in the control volumes
Summarysummary

Contents

NameDescription
Heatprotected
Summaryprotected