modelMethanatorBlock_equilibrium_L2

model of a methanation block
Diagram of MethanatorBlock_equilibrium_L2

Extends from TransiEnt.Basics.Icons.Reactor.

Information

1. Purpose of model

Equilibrium model of reactor block for methanation of hydrogen.

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

The model is based on the equilibrium constants for the reaction equations of CO2-methanation, CO-methanation und CO-Shift. The gas components in the output are in an equilibrium. The model considers the heat of the exothermal reaction and therewith the temperature increase of the product gas and of the temperature increase of the reactor. A detailed description is found in [1].

3. Limits of validity

The calculation of the equilibrium is validated in [1].

4. Interfaces

gasPortIn: Input port for methanation

gasPortOut: Output port of product gas

5. Nomenclature


6. Governing Equations

The chemical balance considers the equilibrium constants. These are calculated via the following formulas:

1. equilibrium constant of CO-methanation: p_CH4*p_H2O/(p_CO*p_H2^3)=9.74*10^(-11)*exp(26830/T_1-30.11)

2. equilibirum constant for CO-shift: p_H2*p_CO2/(p_CO*p_H2O)=exp(4400/T_1-4.063)

The equilibrium constant for CO2-methanation is lineraly dependent on those two equilibrium constants such that the equation can be omited.

Moreover a substance balance and an thermal balance govern the calculation. The substance balance ensures that the amount of subtance of each element entering the reactor block equals the amount of substance leaving the reactor block.

For the thermal balance a simplified model of a fixed-bed reactor is considered with the lenght L, the inner diameter di, a thickness of the reactor wall of delta_wall and an isolation layer with the thickness delta_iso. The cross section of the reactor block is shown in graphic 1.


The following differential equation is needed for the calculation of the temperature of the reactor mass T_reactor:

C* (d/dt *T_reactor)=Q_flow_loss+Q_flow_inner


Here Q_flow_loss describes the heat losses towards the environment and Q_flow_inner is the heat transfer from the process gas onto the reactor mass. Q_flow_loss consists of the heat loss through radiation and convection whereby the necessary values for thermal conductivities and the heat-transfer coefficient can be defined as parameters. The formula for heat-transfer coefficient alpha_i for the heat transfer from the process gas onto the reactor mass (Q_flow_inner) is a simplified correlation which results from a more complex correlation:

alpha_i=[54.2853*(m_flow_total/m_flow_nom)+0.209407]*[-0.000747*T_0+1.465437]


The temperature increase of the process gas is calculated via the following formula:

F_average*c_pgas*(T_in-T_out)=O_flow_inner+Q_flow_reaction


Here F_average is the averaged molar flow between output and input, c_p,gas is the specific heat capacity of the gas, T_out and T_in are the input, respectively the output temperature of the process gas and Q_flow_reaction is the reaction heat that is released during the reaction.

Q_flow_reaction is calculated via the specific molar reaction heat:

Q_flow_reaction=(F_out_CO2-F_in_CO2)*h_3+(F_out_CH4-F_0_CH4)*h_1

Here F_i stands for the molar flow of the respective substance and h1 and h3 stand for the specific molar reaction heat. These specifc molar reaction heat depend on the temperature. As a simplification the reaction heat are calculated with the input temperature T_in:

specific molar reaction heat for CO-methanation: h_1=0.0266*T_0^2-47.7331*T_0-205094.5788

specific molar reaction heat for CO-shift: h_3=0.0026*T_0^2-7.4437*T_0-41557.3842


7. Remarks for Usage

8. Validation

Validated in [1].

9. References

[1] Schülting, Oliver - Vergleich von Power-to-Gas-Speichern mit Ziel der Rückverstromung unter derzeit gültigen technischen Restriktionen (Masterarbeit), Technische Universität Hamburg - Institut für Energietechnik, 2016

10. Version History

Model created by Oliver Schülting (oliver.schuelting@tuhh.de) in Nov 2019

Parameters

TypeNameDefaultDescription
BooleanSteadyStatefalseif '=true' no heating up/cooling down of reactor considered
Modelica.Units.SI.MassFlowRatem_flow_nominal12nominal mass flow
Modelica.Units.SI.AbsolutePressureDelta_p_nominal0.5e5pressure loss for nominal mass flow
BooleanuseHomotopytrueTrue, if homotopy method is used during initialisation
Parameteres for heat transfer
Modelica.Units.SI.CoefficientOfHeatTransferalpha_o25.73outer coefficent of heat transfer
Modelica.Units.SI.ThermalConductivitylambda50thermal conductivity of reactor wall
Modelica.Units.SI.ThermalConductivitylambda_insulation0.04thermal condructivity of insulation layer
Modelica.Units.SI.HeatCapacitymCp_Nenn15E6nominal heat capacity for nominal mass flow 1kg/s
Modelica.Units.SI.Lengthdelta_insulation0.1thickness of insulation layer
Modelica.Units.SI.TemperatureT_ambient283.15constant ambient temperature
Modelica.Units.SI.TemperatureT_reactor_startT_ambientaverage start temperature for reactor mass

Connectors

TypeNameDefaultDescription
TransiEnt.Basics.Interfaces.Gas.RealGasPortIngasPortIn
TransiEnt.Basics.Interfaces.Gas.RealGasPortOutgasPortOut

Components

TypeNameDefaultDescription
TransiEnt.SimCentersimCenter
TILMedia.Internals.VLEFluidConfigurations.FullyMixtureCompatible.VLEFluid_phgasIn_properties
TILMedia.Internals.VLEFluidConfigurations.FullyMixtureCompatible.VLEFluid_pTgasOut_properties
TransiEnt.Basics.Media.Gases.VLE_VDIWA_SG6_varmedium
Realkp1equilibrium constant of methanation reaction
Realkp2equilibrium constant of CO-shift reaction
Modelica.Units.SI.Lengthdelta
Modelica.Units.SI.AbsolutePressureptotal_0total pressure
Modelica.Units.SI.AbsolutePressureptotal_1total pressure
Modelica.Units.SI.AbsolutePressureDelta_p
Modelica.Units.SI.MolarFlowRateFtotal_1output total molar flow
Modelica.Units.SI.MolarFlowRateFtotal_0input total molar flow
Modelica.Units.SI.MolarFlowRate[6]F_0input molar flow
Modelica.Units.SI.MolarFlowRate[6]F_1output molar flow
Modelica.Units.SI.MolarMass[6]Mvector of molar masses
Modelica.Units.SI.MassFraction[6]wt_1weight fraction of output
Modelica.Units.SI.MassFraction[6]wt_0weight fraction of input
Modelica.Units.SI.MassFraction[5]wt_0_5
Modelica.Units.SI.MassFraction[5]wt_1_5
Modelica.Units.SI.MassFlowRatem_flow_totaltotal input mass flow
Modelica.Units.SI.MolarEnergyEnthalpy1enthalpy of methanation reaction 1
Modelica.Units.SI.MolarEnergyEnthalpy2
Modelica.Units.SI.MolarEnergyEnthalpy3
Modelica.Units.SI.EnergyQ_flow_reactionenergy of complete reaction
Modelica.Units.SI.TemperatureT_outoutput Temperature
Modelica.Units.SI.TemperatureT_ininput Temperature
Modelica.Units.SI.TemperatureT_averageaverage Temperature
Modelica.Units.SI.TemperatureT_reactor
Modelica.Units.SI.TemperatureT_outer
Modelica.Units.SI.TemperatureT_out_min200minimum outlet temperature
Modelica.Units.SI.HeatCapacitymCp
Real[4,6]A2[{1.925E1, 5.213E-2, 1.197E-5, -1.132E-8}, {1.980E1, 7.344E-2, -5.602E-5, 1.715E-8}, {3.224E1, 1.924E-3, 1.055E-5, -3.596E-9}, {2.714E1, 9.274E-3, -1.381E-5, 7.645E-9}, {3.087E1, -1.285E-2, 2.789E-5, -1.272E-8}, {28.3, 2.537/1000, 0.5443/1000^2, 0}]coefficient for calculation of heat capacities
Modelica.Units.SI.MolarHeatCapacityCp_averageaverage overall heat capacity
Modelica.Units.SI.MolarHeatCapacityCp_average_2
Modelica.Units.SI.MolarHeatCapacity[6]Cp_0
Modelica.Units.SI.HeatFlowRateQ_flow_transfer
Modelica.Units.SI.HeatFlowRateQ_flow_loss
Modelica.Units.SI.HeatFlowRateQ_flow_loss_convection
Modelica.Units.SI.HeatFlowRateQ_flow_loss_radiation
Modelica.Units.SI.HeatFlowRateQ_flow_inner
Modelica.Units.SI.HeatQ_stored
Modelica.Units.SI.CoefficientOfHeatTransferalpha_i
Modelica.Units.SI.Lengthd_i_ref2.686
Modelica.Units.SI.Lengthd_i
Modelica.Units.SI.Lengthd_o
Modelica.Units.SI.Lengthd_average
Modelica.Units.SI.LengthL_ref4.8
Modelica.Units.SI.LengthL
Modelica.Units.SI.Lengths_min
Modelica.Units.SI.AreaA_o
Modelica.Units.SI.AreaA_i
Modelica.Units.SI.MoleFraction[6]Fshare_1
Modelica.Units.SI.MoleFraction[6]Fshare_0
BooleanOn
Reals1
Reals2
Modelica.Units.SI.Stresssigma_zul