DWSIM/C3/Custom-Unit-Operation-using-Scilab/English

From Script | Spoken-Tutorial
Revision as of 17:09, 30 April 2019 by Kaushik Datta (Talk | contribs)

(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)
Jump to: navigation, search


Visual Cue Narration
Slide Number 1

Title Slide

Welcome to this tutorial on simulating a Custom Unit Operation using Scilab.
Slide Number 2

Learning Objective

In this tutorial, we will learn to:


  • Create a custom unit operation using Scilab
  • Calculate the Outlet pressure of product stream
  • Calculate Molar flow rate and Enthalpy
Slide Number 3

System Requirements

To record this tutorial, I am using


  • DWSIM 5. 6 update 8 (Classic UI) and
  • Windows 10

But, this process is identical in 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 streams and specify their properties
Slide Number 5


Prerequisite Tutorials and Files


The prerequisite tutorials are available on our website.

Spoken-tutorial.org.


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


You should also have Scilab 5.02 (or higher) installed on your system.

Slide Number 6


Scilab CAPE-OPEN Unit Operation

First let us learn what is Scilab CAPE-OPEN Unit Operation.


It is an unit operation implementation for which the calculations can be entered in Scilab.

More information on this is given in the Additional reading material link of this tutorial.

Slide Number 7


Components: Water, methanol and Ethanol


Property Package: Raoult’s Law


Outlet Pressure: Average of Inlet Stream

We will develop a flowsheet to determine the product stream Temperature, Pressure,


Molar flow rate and Molar enthalpy from a Mixer model created using Scilab.

Slide Number 8


Inlet Stream 1:


Component: Pure Water

Temperature: 300 K

Pressure: 100000 Pa

Mass flow: 10 kg/s


Inlet Stream 2:


Component: Pure Methanol

Temperature: 305 K

Pressure: 150000 Pa

Mass flow: 15 kg/s


Inlet Stream 3:


Component: Pure Ethanol

Temperature: 310 K

Pressure: 200000 Pa

Mass flow: 20 kg/s

Here we give Inlet Stream Conditions of the product stream.


Now let me open DWSIM (Classic UI).
Cursor to Desktop


DWSIM (Classic UI) >> Run as Administrator.

To open DWSIM, go to the Desktop.


Right-click on DWSIM (Classic UI) and click Run as Administrator.

User Control Window >> Yes Click Yes on User Access Control window.
Point to Welcome to DWSIM window Welcome to DWSIM window opens.


Close this window.

File >> New Steady-state Simulation. Now go to the 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 the Next button.
Type Water in the Search tab


ChemSep database >> Water

In the Compounds search tab, type Water.


Select Water from the ChemSep database.

Type Methanol in the Search tab Similarly, add Methanol.
Type Ethanol in the Search tab Followed by Ethanol.
Click on Next. And then click on the Next button at the bottom.
Point to Property Packages The Property Packages window opens.
Property Packages >> Available Property Package


Double click on Raoult’s Law

From Available Property Package list, double-click on Raoult’s Law.
Click on Next. Then click on the Next button.
Point to System of Units The next option after this, is System of Units.
System of Units >> SI Under System of Units, select SI.
Click on Finish. Click on the Finish button at the bottom.
Click on Maximize button. Let us now maximize the simulation window.
Cursor on the simulation window. Let’s insert the feed streams that enter the Custom Mixer.
Point to Streams. Go to Streams tab, at the bottom of the main simulation.
Click and drag Material Streams to the flowsheet From the displayed list, drag and drop three Material Streams to the Flowsheet.
Click MSTR-000 Click on the first material stream named as MSTR-000
Type Water Let’s change the name of this stream to Water.
Now we will specify the Water stream properties.
Point to input data. Go to Input Data.
Stream Conditions >> Flash Spec >> Temperature and Pressure (TP) Under Stream Conditions tab, select Flash Spec as Temperature and Pressure (TP).


By default, Temperature and Pressure should already be selected as Flash Spec.

Stream Conditions >>Temperature >> 300 K

Press Enter

Change Temperature to 300 K and press Enter.
Stream Conditions >> Pressure >> 100000 Pa

Press Enter

Change Pressure to 100000 Pa and press Enter.
Stream Conditions >> Mass Flow >> 10 kg/s

Press Enter

Change Mass Flow to 10 kg/s and press Enter.
Input Data >> Compound Amounts Now let us specify the Water stream compositions.


For this, click on Compound Amounts tab under Input Data.

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


By default, Mole Fractions is already selected as Basis.

Water: 1 Now for Water, enter the Amount as 1 and press Enter.
Methanol: 0 For Methanol, enter 0 and press Enter.
Ethanol: 0 For Ethanol, enter 0 and press Enter.
Click Accept Changes button On the right, click on Accept Changes button.
Similarly, we will define the other two material streams.
Click MSTR-001 Click on the second material stream named MSTR-001.
Type Methanol Let’s change the name of this stream to Methanol.
Now we will specify the Methanol 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).

Stream Conditions >>Temperature >> 305 K

Press Enter

Change Temperature to 305 K and press Enter.
Stream Conditions >> Pressure >> 150000 Pa

Press Enter

Change Pressure to 150000 Pa and press Enter.
Stream Conditions >> Mass Flow >> 15 kg/s

Press Enter

Change Mass Flow to 15 kg/s and press Enter.
Input Data >> Compound Amounts Now click on the Compound Amounts tab under Input Data.
Basis >> Mole Fractions Choose Basis as Mole Fractions.
Water: 0 Now for Water, enter the Amount as 0 and press Enter.
Methanol: 1 For Methanol, enter 1 and press Enter.
Ethanol: 0 For Ethanol, enter 0 and press Enter.
Click Accept Changes button. On the right, click on Accept Changes button.
Click MSTR-002 Now, click on the third material stream named as MSTR-002.
Type Ethanol Let’s change the name of this stream to Ethanol.
Now we will specify the Ethanol 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)

Stream Conditions >>Temperature >> 310 K

Press Enter

Change Temperature to 310 K and press Enter.
Stream Conditions >> Pressure >> 200000 Pa

Press Enter

Change Pressure to 200000 Pa and press Enter.
Stream Conditions >> Mass Flow >> 20 kg/s

Press Enter

Change Mass Flow to 20 kg/s and press Enter.
Input Data >> Compound Amounts Then under Input Data click on Compound Amounts tab.
Basis >> Mole Fractions Choose Basis as Mole Fractions.
Water: 0 For Water, enter the Amount as 0 and press Enter.
Methanol: 0 For Methanol, enter 0 and press Enter.
Ethanol: 1 For Ethanol, enter 1 and press Enter.
Click Accept Changes button. On the right side, click on Accept Changes button.
Now, let us insert another material stream that exits the custom mixer.
Click and Drag Material Stream to the flowsheet To do that, let us drag one material stream to the flowsheet.
Let us arrange the streams.
We will leave that stream as unspecified.
Type Mixed Product Then we will change the name of this stream to Mixed Product.
Cursor to CAPE-OPEN tab At the bottom of the main simulation window, go to CAPE-OPEN tab.
Click and drag CAPE-OPEN Unit Operation to the flowsheet Drag and drop CAPE-OPEN Unit Operation to the flowsheet area.
Cursor to CAPE-OPEN Unit Operation window Add CAPE-OPEN Unit Operation window opens.
Point to SciLab Unit Operation.


Click on OK

From the displayed list, select SciLab Unit Operation.


Click on OK at the bottom.

Type Custom Mixer. Let us change the name of this Unit Operation to Custom Mixer.
Now, we will create the inlet and outlet connection ports for the Custom Mixer.
Click Custom Mixer For this, click on Custom Mixer.
Click Open CAPE-OPEN Object Editor From Property Editor window, click on Open CAPE-OPEN Object Editor.
Point to Scilab CAPE-OPEN Unit Operation window Scilab CAPE-OPEN Unit Operation window opens.
Click on Ports tab. Click on Ports tab.
Feed ports >> Add Under Feed ports, there is an Add button at the bottom.
Click on the Add button Click on Add button thrice to create three inlet ports for the three input material streams.
Product ports >> Add Under Product ports, click on Add button at the bottom.


This will create an outlet port for the output material stream.

Click on Close button. Click on Close button at the bottom.
Point to Property Editor Window Go to the Property Editor Window.
Go to Connections. Under Connections, go to inlet tab.
Click on drop down arrow against Feed 1 Select Water. Click on the drop-down against Feed 1 and select Water.
Click on drop down arrow against Feed 2


Select Methanol.

Next, click on the drop-down against Feed 2 and select Methanol.
Click on drop down arrow against Feed 3


Select Ethanol.

Next, click on the drop-down against Feed 3 and select Ethanol.
Click on the Outlet tab. Now click on the Outlet tab.
Click on drop down arrow against Product 1


Select Mixed Product.

Click on the drop-down against Product 1 and

select Mixed Product.

Point to all the connections. All the connections now are complete.


Now we will write the code to perform the computation.

Click on Open CAPE-OPEN Object Editor Click on Open CAPE-OPEN Object Editor.
Click on SciLab tab. Click on SciLab tab.
Cursor on the interface. First we will get the feed properties of the material streams.


Then we will write the calculation routines.

Type:


f1=getFeedProp(1,"flow")

f2=getFeedProp(2,"flow")

f3=getFeedProp(3,"flow")

Type this code in the Script tab.


This is to read the molar flow rate of the compounds from the inlet material streams.

Cursor to f1, f2 and f3 variables. Here, f1, f2 and f3 are the variables.


They store molar flow rate of the compounds from the input material streams.

Cursor to getFeedProp command getFeedProp is the command used to read the specified feed property.
Cursor to 1, 2 and 3 feed port numbers 1, 2 and 3 are the feed port numbers in SciLab to which the material streams in DWSIM is connected.
Cursor to flow keyword flow is the keyword.


flow indicates that components in the streams, has to be stored in f1,f2,f3 in form of an array.

Type:


h1=getFeedProp(1,"enthalpy")

h2=getFeedProp(2,"enthalpy")

h3=getFeedProp(3,"enthalpy")

Type this code in the Script tab.


This is to read the molar enthalpy of the inlet streams from the added inlet material stream.


Here, variables h1, h2 and h3 store the molar enthalpy of the added input material streams.

Cursor to enthalpy keyword enthalpy is the keyword.


It indicate that molar enthalpy of material streams has to be stored in h1, h2 and h3.

Type:


p1=getFeedProp(1,"pressure")

p2=getFeedProp(2,"pressure")

p3=getFeedProp(3,"pressure")

Next, type this code.


This is to read the pressure of the inlet streams from the added inlet material stream.


Here, variables p1, p2 and p3 store the pressures of Water, Methanol and Ethanol streams.

Cursor to pressure keyword pressure is the keyword.


It indicates that the pressure of material streams has to be stored in the variables p1, p2 and p3.

Now we will code the calculation routine
First we will calculate the outlet pressure.


For this let us consider the outlet pressure to be average pressure of all the 3 inlet streams.

p=(p1+p2+p3)/3 Now, type this code.

Outlet pressure is calculated and stored in the variable p.

Now we will calculate the molar flow rate of the compounds in the output stream.
f=f1+f2+f3 Type this code in the Script tab.


Molar flow rate of compounds in the product stream is calculated.


It is stored in the variable f in form of an array.

Now we will calculate the total molar flow rate of the product stream.
totF=sum(f) Type this code in the Script tab.


Molar flow rate of product stream is calculated and stored in the variable totF.

Now we will calculate the enthalpy of the product stream.
h=((h1*f1)+(h2*f2)+(h3*f3))/f Type this code.


Enthalpy of the product stream is calculated and stored in the variable h.

Now we will pass the calculated variables to the Mixed Product material stream.
setProduct(1,totF,f/totF,"pressure",p,"enthalpy",h) Type this code.


setProduct indicates that, variables in the braces are to be passed to the stream in the product.

Click Test At the bottom, click Test.
Point to the Output window Wait for a few seconds for Scilab script to be executed and display the results in the Output tab.
Point to the tabs. Values stored and calculated in all the variables are displayed in the Output tab.
Click Close button. Click Close button at the bottom.
Now we will run the simulation.
Click Solve Flowsheet So, from the toolbar, click on Solve Flowsheet button.
Click Mixed Product When the calculations are completed, click on the Mixed Product in the flowsheet.
Point to Property Editor Window


Hover mouse at Input Data

Under Stream Conditions, check Temperature, Pressure, Molar Flow.
Insert >> Master Property Table Now, go to Insert menu and select Master Property Table.
Double click on Master Property Table Double-click on the Master Property Table.
Point to Configure Master Property Table Configure Master Property Table window opens.
Type Stream Wise Results – Custom Mixer. Enter Name as Stream Wise Results – Custom Mixer.
Type Material Stream Enter Object Type as Material Stream.


By default, Material Stream is already selected.


So we will not change it.

Object >> Liquid Product, Vapour Product and Feed Under Properties to display, select Object as Water, Methanol, Ethanol and Mixed Product.
Configure Master Property Table>> Property Under Property, scroll down to see all the parameters.


Now select the properties as:

Temperature

Pressure

Mass Flow

Molar Flow

Vapor Phase Molar Fraction

Molar Fraction (Mixture) / Water

Molar Flow (Mixture) / Water

Molar Fraction (Mixture) / Methanol

Molar Flow (Mixture) / Methanol

Molar Fraction (Mixture) / Ethanol

Molar Flow (Mixture) / Ethanol

Close Configure Master Property Table window. Let’s close this window.
Point to the Master Property Table


Point to Liquid Product stream

Move the Master Property Table for better visibility.


Here we can see the corresponding results for Water, Methanol, Ethanol and Mixed Product.

Let's summarize.
Slide Number 9


Summary

In this tutorial, we have learnt to


  • Create a custom unit operation using Scilab
  • Calculate the Outlet pressure of product stream
  • Calculate Molar flow rate and Enthalpy
Slide Number 10


Assignment


Compounds:

Water – 100 kmol/h

Methanol – 100 kmol/h


Pressure – 1.25 atm


Property Package: Raoult’s Law


Inlet Stream Temperature: 300 K

Temperature Difference: 350 K

Pressure Drop: 1.1 atm

As an assignment,


Create a custom model for a heater to heat the given mixture of compounds.


Calculate the amount of heat required for heating.


Compare the results with in-built heater unit operation available.

Slide Number 11

About the Spoken Tutorial Project

Watch the video available at following link.

http://spoken-tutorial.org/

It summarizes the Spoken Tutorial project.

Slide Number 12

Spoken Tutorial Workshops

The Spoken Tutorial Project Team
  • Conducts workshops and
  • Gives certificates.
  • For more details, please write to us.
Slide Number 13


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 14

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 15

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 16

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 17

Acknowledgements

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

Thanks

This tutorial is contributed by Kaushik Datta and Priyam Nayak. Thanks for joining.

Contributors and Content Editors

Kaushik Datta, Nancyvarkey