Difference between revisions of "DWSIM/C3/Heterogeneous-Catalytic-Reaction/English"

From Script | Spoken-Tutorial
Jump to: navigation, search
 
(One intermediate revision by one other user not shown)
Line 134: Line 134:
 
|-
 
|-
 
||  
 
||  
|| Before we start the simulation, we will simplify the rate expression.
+
|| Before we start the '''simulation''', we will simplify the '''rate expression'''.
 
|-
 
|-
 
||  
 
||  
|| We will substitute the coefficients in the rate for the numerator and denominator separately.
+
|| We will substitute the coefficients in the '''rate''' for the numerator and denominator separately.
  
  
This will be helpful when we have to enter the reaction rate for '''Heterogeneous Catalytic reaction'''.
+
This will be helpful when we have to enter the reaction '''rate''' for '''Heterogeneous Catalytic reaction'''.
 
|-
 
|-
 
||  
 
||  
Line 171: Line 171:
  
 
4.22E+15*exp(-240100/(8.314*T))/P1<sup>2.5</sup>*((R1*R2)-(P1<sup>3</sup>*P2)/(exp(30.42-27106/T)))
 
4.22E+15*exp(-240100/(8.314*T))/P1<sup>2.5</sup>*((R1*R2)-(P1<sup>3</sup>*P2)/(exp(30.42-27106/T)))
|| The numerator term in the rate expression is this.
+
|| The numerator term in the '''rate expression''' is this.
  
 
Substituting the values of '''k<sub>1'''</sub>, '''K<sub>1'''</sub> and '''partial pressure''' terms, the numerator becomes like this.
 
Substituting the values of '''k<sub>1'''</sub>, '''K<sub>1'''</sub> and '''partial pressure''' terms, the numerator becomes like this.
Line 184: Line 184:
  
 
(1+6.65E-4*exp(38280/(8.314*T))*R1+8.23E-5*exp(70650/(8.314*T))*P2+6.12E-9*exp(82900/(8.314*T))*P1+1.77E+5*exp(-88680/(8.314*T))*R2/P1)^2
 
(1+6.65E-4*exp(38280/(8.314*T))*R1+8.23E-5*exp(70650/(8.314*T))*P2+6.12E-9*exp(82900/(8.314*T))*P1+1.77E+5*exp(-88680/(8.314*T))*R2/P1)^2
|| The denominator term in the rate expression is this.
+
|| The denominator term in the '''rate expression''' is this.
  
 
Substituting the values of '''K<sub>CH4'''</sub>, '''K<sub>CO'''</sub>, '''K<sub>H2'''</sub>, '''K<sub>H2O'''</sub> and '''partial pressure''' terms, the denominator becomes like this.
 
Substituting the values of '''K<sub>CH4'''</sub>, '''K<sub>CO'''</sub>, '''K<sub>H2'''</sub>, '''K<sub>H2O'''</sub> and '''partial pressure''' terms, the denominator becomes like this.
 
|-
 
|-
 
||  
 
||  
|| We have simplified the rate expression to be entered in '''DWSIM'''.
+
|| We have simplified the '''rate expression''' to be entered in '''DWSIM'''.
  
 
Now we will start the simulation.
 
Now we will start the simulation.
Line 214: Line 214:
 
|| In the '''Compounds Search '''tab, type '''Methane.'''
 
|| In the '''Compounds Search '''tab, type '''Methane.'''
  
Select '''Methane '''from '''ChemSep '''database.
+
Select '''Methane '''from '''ChemSep database'''.
 
|-
 
|-
 
|| Type '''Water '''in the '''Search''' tab
 
|| Type '''Water '''in the '''Search''' tab
|| Next, add '''Water''' from '''ChemSep '''database.
+
|| Next, add '''Water''' from '''ChemSep database'''.
 
|-
 
|-
 
|| Type '''Hydrogen '''in the '''Search''' tab
 
|| Type '''Hydrogen '''in the '''Search''' tab
|| Then, add '''Hydrogen''' from '''ChemSep '''database.
+
|| Then, add '''Hydrogen''' from '''ChemSep database'''.
 
|-
 
|-
 
|| Type '''Carbon monoxide '''in the '''Search''' tab
 
|| Type '''Carbon monoxide '''in the '''Search''' tab
|| Next, add '''Carbon monoxide''' from '''ChemSep '''database.
+
|| Next, add '''Carbon monoxide''' from '''ChemSep database'''.
 
|-
 
|-
 
|| Type '''Carbon dioxide '''in the '''Search''' tab
 
|| Type '''Carbon dioxide '''in the '''Search''' tab
|| Next, add '''Carbon dioxide''' from '''ChemSep '''database.
+
|| Next, add '''Carbon dioxide''' from '''ChemSep database'''.
  
  
Line 240: Line 240:
  
 
Double click on '''Peng-Robinson (PR)'''
 
Double click on '''Peng-Robinson (PR)'''
|| From '''Available Property Packages''' list, double click on '''Peng-Robinson (PR) '''option.
+
|| From '''Available Property Packages''' list, double-click on '''Peng-Robinson (PR) '''option.
 
|-
 
|-
 
|| Click on '''Next''' button at the bottom.
 
|| Click on '''Next''' button at the bottom.
Line 255: Line 255:
 
|-
 
|-
 
|| Click on '''Maximize''' button
 
|| Click on '''Maximize''' button
|| Let us maximize the simulation window.
+
|| Let us maximize the '''simulation''' window.
 
|-
 
|-
 
|| Cursor on the simulation window.
 
|| Cursor on the simulation window.
Line 261: Line 261:
 
|-
 
|-
 
|| Point to '''Streams.'''
 
|| Point to '''Streams.'''
|| At the bottom of the main simulation window, go to '''Streams.'''
+
|| At the bottom of the main '''simulation''' window, go to '''Streams.'''
 
|-
 
|-
 
|| Click and drag '''Material Stream '''to the flowsheet
 
|| Click and drag '''Material Stream '''to the flowsheet
Line 267: Line 267:
 
|-
 
|-
 
|| Click '''MSTR-000'''
 
|| Click '''MSTR-000'''
|| Click on the material stream named '''MSTR-000.'''
+
|| Click on the '''material stream''' named '''MSTR-000.'''
 
|-
 
|-
 
|| Type '''Feed'''
 
|| Type '''Feed'''
Line 273: Line 273:
 
|-
 
|-
 
||  
 
||  
|| Now we will specify the Feed''' stream''' properties.
+
|| Now we will specify the '''Feed stream''' properties.
 
|-
 
|-
 
|| '''Stream Conditions >> Flash Spec >> Temperature and Pressure (TP)'''
 
|| '''Stream Conditions >> Flash Spec >> Temperature and Pressure (TP)'''
Line 299: Line 299:
 
|-
 
|-
 
|| '''Input Data '''>>''' Compound Amounts'''
 
|| '''Input Data '''>>''' Compound Amounts'''
|| Now let us specify the '''Feed stream compositions.'''
+
|| Now let us specify the '''Feed stream''' compositions.
  
 
For this under '''Input Data''' click on '''Compound Amounts''' tab.
 
For this under '''Input Data''' click on '''Compound Amounts''' tab.
Line 413: Line 413:
  
  
Then press Enter.
+
Then press '''Enter'''.
 
|-
 
|-
 
|| Point to '''Stoichiometry '''field
 
|| Point to '''Stoichiometry '''field
Line 419: Line 419:
 
|-
 
|-
 
|| Point to '''Equation '''field
 
|| Point to '''Equation '''field
|| Here''' '''the '''Equation '''field shows the reaction equation'''.'''
+
|| Here''' '''the '''Equation '''field shows the reaction equation.
 
|-
 
|-
 
|| Point to '''Heterogeneous Kinetic Reaction Parameters'''
 
|| Point to '''Heterogeneous Kinetic Reaction Parameters'''
Line 434: Line 434:
 
|-
 
|-
 
||  
 
||  
|| Now we will specify the reaction rate.
+
|| Now we will specify the '''reaction rate'''.
 
|-
 
|-
 
|| Point to '''Reaction Rate(Base Component)'''
 
|| Point to '''Reaction Rate(Base Component)'''
Line 440: Line 440:
 
|-
 
|-
 
|| Point to '''Numerator''' and '''Denominator''' field.
 
|| Point to '''Numerator''' and '''Denominator''' field.
|| The reaction rate is to be entered separately as '''Numerator''' and '''Denominator.'''
+
|| The '''reaction rate''' is to be entered separately as '''Numerator''' and '''Denominator.'''
  
  
Line 469: Line 469:
 
|-
 
|-
 
|| '''Amount Unit >> atm'''
 
|| '''Amount Unit >> atm'''
|| Select '''Amount Unit''' as '''atm''' from the dropdown.
+
|| Select '''Amount Unit''' as '''atmosphere''' from the drop-down.
  
The '''Reaction Rate '''is expressed in terms of '''Partial Pressure, '''whose unit is '''atm.'''  
+
The '''Reaction Rate '''is expressed in terms of '''Partial Pressure, '''whose unit is '''atmosphere.'''  
 
|-
 
|-
 
|| '''Velocity Unit >> kmol/[kg.h]'''
 
|| '''Velocity Unit >> kmol/[kg.h]'''
Line 488: Line 488:
 
|-
 
|-
 
|| Cursor to''' Reactors''' tab
 
|| Cursor to''' Reactors''' tab
|| At the bottom of the main simulation window, go to '''Reactors '''tab.
+
|| At the bottom of the main '''simulation''' window, go to '''Reactors '''tab.
 
|-
 
|-
 
|| Click and drag a '''Plug-Flow Reactor(PFR)''' to the flowsheet
 
|| Click and drag a '''Plug-Flow Reactor(PFR)''' to the flowsheet
Line 503: Line 503:
 
|-
 
|-
 
||  
 
||  
|| Let us once again arrange it in the flowsheet.
+
|| Let us once again arrange it in the '''flowsheet'''.
 
|-
 
|-
 
||  
 
||  
|| Leave that stream as unspecified.
+
|| Leave that '''stream''' as unspecified.
 
|-
 
|-
 
|| Type '''Product'''
 
|| Type '''Product'''
|| We will change the name of this stream to '''Product.'''
+
|| We will change the name of this '''stream''' to '''Product.'''
 
|-
 
|-
|| Click and drag '''Energy Stream''' to the flowsheet
+
|| Click and drag '''Energy Stream''' to the '''flowsheet'''
 
|| Let us now insert one '''Energy Stream.'''
 
|| Let us now insert one '''Energy Stream.'''
 
|-
 
|-
 
|| Type '''Energy.'''
 
|| Type '''Energy.'''
|| Name this stream as '''Energy.'''
+
|| Name this '''stream''' as '''Energy.'''
 
|-
 
|-
 
||Click '''Plug-Flow Reactor (PFR)'''
 
||Click '''Plug-Flow Reactor (PFR)'''
Line 547: Line 547:
 
|-
 
|-
 
|| Hover mouse at '''Calculation Parameters'''
 
|| Hover mouse at '''Calculation Parameters'''
|| Now, Go to the '''Calculation Parameters.'''
+
|| Now, go to the '''Calculation Parameters.'''
 
|-
 
|-
 
|| '''Reaction Set''' >> '''Default Set'''
 
|| '''Reaction Set''' >> '''Default Set'''
|| In this section, first option is '''Reaction Set.'''
+
|| In this section, the first option is '''Reaction Set.'''
  
 
By default, it is '''Default Set'''.
 
By default, it is '''Default Set'''.
Line 559: Line 559:
  
 
Select''' Isothermic'''
 
Select''' Isothermic'''
|| Click on the drop down against '''Calculation Mode'''
+
|| Click on the drop down against '''Calculation Mode''' and select '''Isothermic.'''
 
+
and select '''Isothermic.'''
+
 
|-
 
|-
 
|| '''Reactor Volume''' >> 1''' m3'''
 
|| '''Reactor Volume''' >> 1''' m3'''
Line 587: Line 585:
  
 
Press '''Enter'''
 
Press '''Enter'''
|| Click on the field against '''Catalyst Particle Diameter''' and enter it as 2''' mm.'''
+
|| Click on the field against '''Catalyst Particle Diameter''' and enter it as 2''' mm'''
  
 
and press '''Enter.'''
 
and press '''Enter.'''
Line 601: Line 599:
 
|| Wait for a few seconds.
 
|| Wait for a few seconds.
  
The flowsheet is getting simulated.
+
The flowsheet is getting '''simulated'''.
 
|-
 
|-
 
||  
 
||  
|| Once again, we will run the simulation.
+
|| Once again, we will run the '''simulation'''.
 
|-
 
|-
 
|| Click '''Solve Flowsheet'''
 
|| Click '''Solve Flowsheet'''
Line 610: Line 608:
 
|-
 
|-
 
|| Click '''Plug-Flow Reactor'''
 
|| Click '''Plug-Flow Reactor'''
|| When calculations are completed, click on the '''PFR''' in the Flowsheet.  
+
|| When calculations are completed, click on the '''PFR''' in the flowsheet.  
 
|-
 
|-
 
|| Hover mouse at '''Results'''
 
|| Hover mouse at '''Results'''
Line 616: Line 614:
 
|-
 
|-
 
|| '''Results >> General'''
 
|| '''Results >> General'''
|| Under '''General''' tab, Check '''Residence time.'''
+
|| Under '''General''' tab, check '''Residence time.'''
  
 
It is '''0.00068 h.'''
 
It is '''0.00068 h.'''
Line 624: Line 622:
  
 
For '''Methane,''' the conversion is '''80.67%''' and for '''Water,''' it is '''80.67%.'''
 
For '''Methane,''' the conversion is '''80.67%''' and for '''Water,''' it is '''80.67%.'''
|-
 
|| '''Insert >> Master Property Table'''
 
|| To check the material balances, go to '''Insert''' option from the toolbar.
 
 
Select '''Master Property Table.'''
 
|-
 
|| Double click on '''Master Property Table'''
 
 
Point to '''Configure Master Property Table'''
 
|| Double click on the '''Master Property Table '''to edit it.
 
 
'''Configure Master Property Table''' window opens.
 
|-
 
|| Type '''Stream Wise Results – Heterogeneous Catalytic Reaction'''
 
|| Enter the '''Name''' as '''Stream Wise Results – Heterogeneous Catalytic Reaction.'''
 
|-
 
|| Type '''Material Stream'''
 
|| Type '''Object Type''' as '''Material Stream'''.
 
 
By default, '''Material Stream''' is already selected.
 
 
So we will not change it.
 
|-
 
|| '''Object >> Feed''' and '''Product'''
 
|| Under '''Properties to display''', select '''Object '''as '''Feed''' and '''Product.'''
 
|-
 
|| '''Configure Master Property Table>> Property'''
 
|| Under '''Property''', select the properties as
 
 
'''Temperature'''
 
 
'''Pressure'''
 
 
'''Mass Flow'''
 
 
'''Molar Flow'''
 
 
'''Molar Flow (Mixture) / Methane'''
 
 
'''Mass Flow (Mixture) / Methane'''
 
 
'''Molar Flow (Mixture) / Water'''
 
 
'''Mass Flow (Mixture) / Water'''
 
 
'''Molar Flow (Mixture) / Hydrogen'''
 
 
'''Mass Flow (Mixture) / Hydrogen'''
 
 
'''Molar Flow (Mixture) /Carbon monoxide'''
 
 
'''Mass Flow (Mixture) / Carbon monoxide'''
 
 
'''Molar Flow (Mixture) / Carbon dioxide''' and
 
 
'''Mass Flow (Mixture) / Carbon dioxide'''
 
|-
 
|| Close this window.
 
|| Close this window.
 
|-
 
|| Point to the results.
 
|| Move the '''Master Property Table '''for better visibility.
 
  
Here we can see the corresponding results for '''Product''' and''' Feed.'''
 
 
|-
 
|-
 
||  
 
||  
Line 718: Line 653:
 
||As an assignment,
 
||As an assignment,
  
Add two more Heterogeneous Catalytic Reactions for the existing PFR system.
+
Add two more '''Heterogeneous Catalytic reactions''' for the existing '''PFR system'''.
  
 
'''Feed Stream''' conditions and '''Plug-Flow Reactor''' parameters remain unchanged.
 
'''Feed Stream''' conditions and '''Plug-Flow Reactor''' parameters remain unchanged.
Line 736: Line 671:
 
K<sub>3</sub> = exp (34.218 - 31266/T)
 
K<sub>3</sub> = exp (34.218 - 31266/T)
  
|| Here we give the reaction rates and its coefficients.  
+
|| Here we give the '''reaction rates''' and its coefficients.  
 
|-
 
|-
 
|| '''Slide Number 16'''
 
|| '''Slide Number 16'''
Line 742: Line 677:
 
'''About the Spoken Tutorial Project'''
 
'''About the Spoken Tutorial Project'''
 
|| Watch the video available at following link.
 
|| Watch the video available at following link.
 
'''http://spoken-tutorial.org/ '''
 
  
 
It summarizes the Spoken Tutorial project.
 
It summarizes the Spoken Tutorial project.
Line 750: Line 683:
  
 
'''Spoken Tutorial Workshops'''
 
'''Spoken Tutorial Workshops'''
|| The Spoken Tutorial Project Team
+
|| The Spoken Tutorial Project Team conducts workshops and gives certificates.
  
* Conducts workshops and
+
For more details, please write to us.
* Gives certificates.
+
* For more details, please write to us.
+
 
|-
 
|-
 
|| '''Slide Number 18'''
 
|| '''Slide Number 18'''

Latest revision as of 11:32, 6 August 2019

Visual Cue Narration
Slide Number 1

Title Slide

Welcome to this tutorial on Simulating a PFR using Heterogeneous Catalytic Reaction in DWSIM.
Slide Number 2

Learning Objective

In this tutorial, we will learn to:
  • Define a Heterogeneous Catalytic reaction
  • Define parameters for Plug Flow Reactor required to simulate a HCR
  • Calculate Conversion and Residence time for HCR in a PFR
Slide Number 3

System Requirements

To record this tutorial, I am using
  • DWSIM 5.6 (Classic UI) Update 8 and
  • Windows 10

The process demonstrated in this tutorial is identical in other OS also such as-

  • Ubuntu Linux,
  • Mac OS X or
  • FOSSEE OS on ARM.
Slide Number 4


Pre-requisites

To practice this tutorial, you should know to-
  • add components to a flowsheet
  • select thermodynamic packages
  • add material and energy streams and specify their properties.
Slide Number 5

Prerequisite Tutorials and Files

https:\\spoken-tutorial.org

The prerequisite tutorials are available on this website.

You can access these tutorials and all the associated files from this site.

Slide Number 6

Problem Statement

Reaction: CH4 + H2O ⇔ 3H2 + CO

Inlet Stream:

Mass Flow: 300 kg/h

Mole Fraction(CH4): 0.4975

Mole Fraction(H2O): 0.4975

Mole Fraction(H2): 0.005

Temperature: 730 degree C

Pressure: 1.01325 bar


Property Package: Peng Robinson(PR)

PFR Dimensions:

Volume: 1 m3, length: 1 m

We will develop a flowsheet to determine the exit composition from an Isothermal PFR with Heterogeneous Catalytic Reaction.



Here we give Inlet stream conditions, Property package and Reactor Parameters.

Slide Number 7


Catalyst Properties:

Loading: 0.386 kg/m3

Particle Diameter: 0.002 m

Void Fraction: 0.4


Gas Constant: R = 8.314 J/mol-K

Next, we give the Catalyst properties and Gas constant.
Slide Number 8

Reaction Rate and Coefficients:

R1 = [k1/PH22.5 {PCH4 PH2O - PH23 PCO / K1}] / DEN2


Reaction Rate

k1 = 4.22 * 10 15 exp (- 240100/RT)

K1 = exp (30.42 - 27106/T)

DEN = 1 + KCH4 PCH4 + KCO PCO + KH2 PH2 + KH2O PH2O/PH2

Now we give the reaction rate and its coefficients.
Slide Number 9

Reaction Rate Coefficients

k1 = 6.65 * 10-4 exp (38280/RT)


k1 = 1.77 * 105 exp (- 88680/RT)


k1 = 6.12 * 10-9 exp (82920/RT)


k1 = 8.23 * 10-5 exp (70650/RT)

Here we give the reaction rate coefficients.
Before we start the simulation, we will simplify the rate expression.
We will substitute the coefficients in the rate for the numerator and denominator separately.


This will be helpful when we have to enter the reaction rate for Heterogeneous Catalytic reaction.

In the rate expression, we will assign the variables for partial pressure of reactants and products.
Slide Number 10

Reactants: R

PCH4: R1, PH2O: R2

Products: P

PH2: P1, PCO: P2

R corresponds to reactants.

Partial pressure of Methane(CH4) is R1 and Water(H2O) is R2.

P corresponds to products.

So the partial pressure of Hydrogen(H2) is P1 and Carbon monoxide(CO) is P2.

Slide Number 11

Numerator of R1 = k1/PH22.5 {PCH4 PH2O - PH23 PCO / K1}


Substituting for k1, PH2, PCH4, PH3, PCO and K1, we get


4.22E+15*exp(-240100/(8.314*T))/P12.5*((R1*R2)-(P13*P2)/(exp(30.42-27106/T)))

The numerator term in the rate expression is this.

Substituting the values of k1, K1 and partial pressure terms, the numerator becomes like this.

Slide Number 12

Denominator of R1 = [1 + KCH4 PCH4 + KCO PCO + KH2 PH2 + KH2O PH2O/PH2] 2


Substituting for KCH4, PCH4, KCO, PCO, KH2, PH2, KH2O and PH2O, we get


(1+6.65E-4*exp(38280/(8.314*T))*R1+8.23E-5*exp(70650/(8.314*T))*P2+6.12E-9*exp(82900/(8.314*T))*P1+1.77E+5*exp(-88680/(8.314*T))*R2/P1)^2

The denominator term in the rate expression is this.

Substituting the values of KCH4, KCO, KH2, KH2O and partial pressure terms, the denominator becomes like this.

We have simplified the rate expression to be entered in DWSIM.

Now we will start the simulation.

File >> New Steady-state Simulation I have already opened DWSIM.

Go to File menu and select New Steady-state Simulation.

Point to Simulation Configuration Wizard window Simulation Configuration Wizard window appears.
Click on Next button at the bottom. Click on Next button at the bottom.
First, we will add the compounds.
Type Methane in the Search tab

ChemSep database >> Methane

Check the checkbox under column Added

In the Compounds Search tab, type Methane.

Select Methane from ChemSep database.

Type Water in the Search tab Next, add Water from ChemSep database.
Type Hydrogen in the Search tab Then, add Hydrogen from ChemSep database.
Type Carbon monoxide in the Search tab Next, add Carbon monoxide from ChemSep database.
Type Carbon dioxide in the Search tab Next, add Carbon dioxide from ChemSep database.


Now, all the compounds are added.

Click on Next button at the bottom. Click on Next button at the bottom.
Point to Property Packages Now comes Property Packages.
Property Packages >> Available Property Package

Double click on Peng-Robinson (PR)

From Available Property Packages list, double-click on Peng-Robinson (PR) option.
Click on Next button at the bottom. Then click on the Next button at the bottom.
Point to System of Units Next option is System of Units.
System of Units >> C5 Under System of Units, we will select C5.
Click on Finish Then at the bottom, click on the Finish button.
Click on Maximize button Let us maximize the simulation window.
Cursor on the simulation window. Now let's insert a feed stream that enters the Reactor.
Point to Streams. At the bottom of the main simulation window, go to Streams.
Click and drag Material Stream to the flowsheet From the displayed list, drag and drop a Material Stream to the Flowsheet.
Click MSTR-000 Click on the material stream named MSTR-000.
Type Feed Let’s change the name of this stream to Feed.
Now we will specify the Feed stream properties.
Stream Conditions >> Flash Spec >> Temperature and Pressure (TP) Go to Input Data.

Under Stream Conditions tab, select Flash Spec as Temperature and Pressure (TP), if not already selected.

Point to the default temperature and pressure. By default, Temperature and Pressure are already selected as Flash Spec.
Stream Conditions >>Temperature >> 730 C

Press Enter

Change Temperature to 730 degree C and press Enter.
Stream Conditions >> Pressure >> 1.01325 bar

Press Enter

Change Pressure to 1.01325 bar and press Enter.
Stream Conditions >> Mass Flow >> 300 kg/h

Press Enter

Change Mass Flow to 300 kg/h and press Enter.
Input Data >> Compound Amounts Now let us specify the Feed stream compositions.

For this under Input Data click on Compound Amounts tab.

Basis >> Mole Fractions Choose the Basis as Mole Fractions, if not already selected.

By default, Mole Fractions is already selected as Basis.

Methane: 0.4975 Now for Methane, enter the Amount as 0.4975 and press Enter.
Water: 0.4975 For Water, enter 0.4975 and press Enter.
Hydrogen: 0.005 For Hydrogen, enter 0.005 and press Enter.
Carbon monoxide: 0 For Carbon monoxide, enter 0 and press Enter.
Carbon dioxide: 0 For Carbon dioxide, enter 0 and press Enter.
Click Accept Changes button On the right, click on the Accept Changes button.
Now we will define the Heterogeneous Catalytic Reaction.
Tools >> Reaction Manager

Point to Chemical Reactions Manager

Under Tools, click on Reactions Manager.

Chemical Reactions Manager window opens.

Chemical Reactions >> Add Reaction Under Chemical Reactions tab, click on the green coloured Add Reaction button.
Click on Heterogeneous Catalytic Click on Heterogeneous Catalytic.
Point to Heterogeneous Catalytic Reaction Heterogeneous Catalytic Reaction window opens.
Identification >> Name >> Steam Reforming Under Identification, enter the Name as Steam Reforming.
Description >> Catalytic Reaction for Steam Reforming Next, type the Description as,

Catalytic Reaction for Steam Reforming

Point to Components and Stoichiometry Next part is a table of Components and Stoichiometry.
Point to Name field First column Name, shows the available components here.
Point to Molar Weight Second column corresponds to its Molar Weight.
Point to Heat of Formation Third column corresponds to its Heat of Formation.
Point to Include Next column is Include.


Under Include, check the check boxes against Methane, Water, Hydrogen and Carbon monoxide.

Carbon dioxide doesn’t take part in this reaction.


We will keep it unchecked.

Point to BC


Check Methane check box

Fourth column is BC.

Under BC, check the Methane check box as Methane is the base component.

Point to SC Next column is SC (stoichiometric coefficient)

Stoich. Coeff >> Methane: -1,

Water: -1

Hydrogen: 3,

Carbon monoxide: 1,

Carbon dioxide: 0'

Under SC column, enter


-1 for Methane

-1 for Water

3 for Hydrogen

1 for Carbon monoxide

0 for Carbon dioxide


Then press Enter.

Point to Stoichiometry field In the Stoichiometry field, we can see it is showing OK.
Point to Equation field Here the Equation field shows the reaction equation.
Point to Heterogeneous Kinetic Reaction Parameters Next part is Heterogeneous Kinetic Reaction Parameters.
Basis >> Partial Pressures Select Basis as Partial Pressures.
Phase >> Vapor Select Phase as Vapor.
Leave Tmin and Tmax unchanged. We will leave Tmin and Tmax field as unchanged.
Now we will specify the reaction rate.
Point to Reaction Rate(Base Component) Go to Reaction Rate (Base Component) = Numerator/Denominator.
Point to Numerator and Denominator field. The reaction rate is to be entered separately as Numerator and Denominator.


We have already simplified the Numerator and Denominator part separately.

Point to Numerator


Type

4.22E+15*exp(-240100/(8.314*T))/P1^2.5*((R1*R2)-(P1^3*P2)/(exp(30.42-27106/T)))

We will enter the Numerator part first.


Type the equation as shown here.

Point to Denominator

Type

(1+6.65E-4*exp(38280/(8.314*T))*R1+8.23E-5*exp(70650/(8.314*T))*P2+6.12E-9*exp(82900/(8.314*T))*P1+1.77E+5*exp(-88680/(8.314*T))*R2/P1)^2

Now we will enter the Denominator part.

Type the equation as shown here.

Now we will specify the units for Amount and Velocity.
Amount Unit >> atm Select Amount Unit as atmosphere from the drop-down.

The Reaction Rate is expressed in terms of Partial Pressure, whose unit is atmosphere.

Velocity Unit >> kmol/[kg.h] From the drop-down select Velocity Unit as kmol/[kg.h].

This is the unit for the Reaction Rate.

Click OK Click on the OK button at the bottom.
Click on Flowsheet Now go to the Flowsheet area.
We will use a Plug Flow Reactor for the simulation.
Cursor to Reactors tab At the bottom of the main simulation window, go to Reactors tab.
Click and drag a Plug-Flow Reactor(PFR) to the flowsheet Drag and drop a Plug-Flow Reactor(PFR) to the flowsheet area.
Let us arrange it as required.
And then let us insert one Output Stream.
Click and drag Material Stream to the flowsheet To do that let us drag one Material Stream.
Let us once again arrange it in the flowsheet.
Leave that stream as unspecified.
Type Product We will change the name of this stream to Product.
Click and drag Energy Stream to the flowsheet Let us now insert one Energy Stream.
Type Energy. Name this stream as Energy.
Click Plug-Flow Reactor (PFR) We are now ready to specify the Plug-Flow Reactor.

So let’s click on it.

Point to the tab on the left. On the left, we can see a tab displaying properties related to the PFR.
Go to aConnections


Click on drop down arrow against Inlet Stream


Select Feed.

Under Connections, click on the drop down against Inlet Stream.

Select Feed.

Click on drop down arrow against Outlet Stream

Select Product.

Click on the drop down against Outlet Stream and select Product.
Click on the drop down against Energy Stream

Select Energy.

Click on the drop down against Energy Stream and select Energy.
Hover mouse at Calculation Parameters Now, go to the Calculation Parameters.
Reaction Set >> Default Set In this section, the first option is Reaction Set.

By default, it is Default Set.

Since we have only one reaction, we leave it as it is.

Click drop down against Calculation Mode

Select Isothermic

Click on the drop down against Calculation Mode and select Isothermic.
Reactor Volume >> 1 m3

Press Enter

Click on the field against Reactor Volume and enter 1 m3.

and press Enter.

Reactor Length >> 1 m

Press Enter

Click on the field against Reactor length and enter it as 1 m.

Then press Enter.

Catalyst Loading >> 0.386 kg/m3

Press Enter

Click on the field against Catalyst Loading and enter it as 0.386 kg/m3.

and press Enter.

Catalyst Particle Diameter >> 2 mm

Press Enter

Click on the field against Catalyst Particle Diameter and enter it as 2 mm

and press Enter.

Catalyst Void Fraction >> 0.4

Press Enter

Click on the field against Catalyst Void Fraction and enter it as 0.4.

Then press Enter.

Wait for a few seconds.

The flowsheet is getting simulated.

Once again, we will run the simulation.
Click Solve Flowsheet To do this, click on Solve Flowsheet button.
Click Plug-Flow Reactor When calculations are completed, click on the PFR in the flowsheet.
Hover mouse at Results Locate Results section.
Results >> General Under General tab, check Residence time.

It is 0.00068 h.

Results >> Conversions Under Conversions tab, check conversion for both the reacting compounds.

For Methane, the conversion is 80.67% and for Water, it is 80.67%.

Let's summarize.
Slide Number 13

Summary

In this tutorial, we have learnt to
  • Define a Heterogeneous Catalytic reaction
  • Define parameters for Plug Flow Reactor required to simulate a HCR
  • Calculate Conversion and Residence time for HCR in a PFR
Slide Number 14

Assignment

Water Gas Shift Reaction:

CO + H2O ⇔ H2 + CO2

R2 = [k2/PH2 {PCO PH2O - PH2 PCO2 / K2}] / DEN2

k2 = 1.96 * 10 6 exp (- 67130/RT)

K2 = exp (-3.798 + 4160/T)

As an assignment,

Add two more Heterogeneous Catalytic reactions for the existing PFR system.

Feed Stream conditions and Plug-Flow Reactor parameters remain unchanged.

Slide Number 15

Assignment

Overall Reaction:

CH4 + 2H2O ⇔ 4H2 + CO2

R3 = [k3/PH23.5 {PCH4 PH2O2 - PH24 PCO2 / K3}] / DEN2

k3 = 1.02 * 10 15 exp (- 243900/RT)

K3 = exp (34.218 - 31266/T)

Here we give the reaction rates and its coefficients.
Slide Number 16

About the Spoken Tutorial Project

Watch the video available at following link.

It summarizes the Spoken Tutorial project.

Slide Number 17

Spoken Tutorial Workshops

The Spoken Tutorial Project Team conducts workshops and gives certificates.

For more details, please write to us.

Slide Number 18

Forum slide


Do you have questions in this Spoken Tutorial?

Please visit this site

Choose the minute and second where you have the question.

Explain your question briefly.

Someone from the FOSSEE team will answer them.

Please post your times queries in this forum.
Slide Number 19

DWSIM Flowsheeting Project

The FOSSEE team coordinates conversion of existing flow sheets into DWSIM.

We give honorarium and certificates.

For more details, please visit this site.

Slide Number 20

TextBook Companion Project

The FOSSEE team coordinates coding of solved examples of popular books.

We give honorarium and certificates.

For more details, please visit this site.

Slide Number 21

Lab Migration Project

The FOSSEE team helps migrate commercial simulator labs to DWSIM.

We give honorarium and certificates.

For more details, please visit this site.

Slide Number 22

Acknowledgements

Spoken Tutorial and FOSSEE projects are funded by NMEICT, MHRD, Government of India.
Slide Number 23

Thanks

This tutorial is contributed by Kaushik Datta and Priyam Nayak.

Thanks for joining.

Contributors and Content Editors

Kaushik Datta, Nancyvarkey