Difference between revisions of "ASCEND/C2/Solving-Cubic-EOS/English"

Title slide

We will learn how to

• Create a reusable model
• Reuse the model for different components
• Import data from the Thermodynamics library

• Solve the model

System Requirement

Ubuntu Linux OS v. 12.04

ASCEND v. 0.9.8

Pre-requisites

Linux

ASCEND and

Peng Robinson cubic EOS

Pre-requisites

For relevant tutorials please visit our website

http://spoken-tutorial.org

Peng Robinson Cubic EOS

Peng Robinson equation of state is a cubic equation of state which contains volume terms to the third power.

It is usually expressed to give pressure in terms of temperature and molar volume. These are the equation for Peng Robinson equation of state

Peng Robinson Cubic EOS

These are the equations which we will be using in our model to compute the value of Z

I already have the code for Peng Robinson model in a text file.

Open the file Peng Robinson.a4c in a text editor

Let me explain the code.

So, add REQUIRE “thermodynamics.a4l”.

After this define variables like P, T, V etc. using WILL_BE operator.

WILL_BE operator is used to make our model reusable

The variables defined using WILL_BE are used as parameters in our test model

I will come back to the explanation of WILL_BE as we move ahead in our tutorial

R IS_A molar_gas_constant;

Pc IS_A pressure;

Tc IS_A temperature;

T_degC IS_A factor;

alpha, beta, q,Tr, Pr, omega IS_A factor;

PSI, OMEGA, sigma, eps IS_A real_constant;

Note that factor denotes a dimensionless variable, Thus alpha, beta, q etc. are defined as factor

PSI :== 0.45724;

OMEGA :== 0.07780;

sigma :== 2.414213562;

eps :== 0.414213562;

Please note that the values of these constants will be reused in our test model.

So we use a colon and is equal to sign twice to assign values to constants

Pc = data.Pc;

Tc = data.Tc;

omega = data.omega;

Tr = T/Tc;

Pr = P/Pc;

T_degC = T/1{K} - 273.15;

And so we write

Pc=data.Pc and so on.

Now we define the reduced temperature and pressure in terms of Tc and Pc.

Now we mention the Peng Robinson equation as shown in the previous slide

Note that labels like eq 4, eq 5 etc. are not compulsory to write but they help in debugging the problem

We use the method default to solve this model and run the METHODS specify and values under it.

Under specify we fix the variables by writing , P.fixed := TRUE;and T.fixed := TRUE;

Note that True implies fixed variable whereas all other variables that are not mentioned here are by default set to False, i.e. Free variables

Now under values, we define values to the fixed variables

For Z, we set the upper and lower bound and also provide a nominal value to initiate the solution.

The syntax used is

Z.lower_bound := 0.0;

Z.nominal := 0.1;

Now we end Peng Robinson model

Parameters defined in this model, can now be reused by any model to calculate Z for specific components

Let's now see how to use these parameters and calculate Z for ethylene in a test model.

Let's now look at the code for test model

Note that the test model is written below the Peng Robinson model in the same document

We name this model as test

We define a symbol constant – c1.

And assign c1 is equal to ethylene.

We use a colon and equal to sign twice as c1 is passed as a parameter to another variable components_data which is present in thermodynamics library

Here we define cd as components_data of ethylene

Note that in the parenthesis,c1 in square bracket implies such data for ethylene from thermodynamics library whereas c1 outside square bracket is set as reference component

More on this will be explained in future tutorials

Note that we use WILL_BE to define P,V,Z etc

To reuse these parameters from previous model we need to define them here using IS_A operator

Thus we define

P IS_A pressure;

T IS_A temperature and so on;

We define tpeng IS_A pengrobinson

Now tpeng includes all the equations and methods

defined in the peng robinson model

We have passed the required parameters to pengrobinson model in parenthesis,so we don't have to write the variables and the equations again

the values of constants and other variables will be called in the test model from peng robinson model

This is how we make our model reusable

We define equations for Tr and Pr as

Tr = T/cd.data[c1].Tc; and similarly for Pr;

We reuse this model in the test model to find out Z for ethylene.

I already have the code for testpengrobinson model in a text file.

Let me explain the code.

We name this model as testpengrobinson.

We define a symbol constant – c1.

And assign c1 is equal to ethylene.

This c1 is passed as a parameters to another model components_data which is present in one of our ‘REQUIRED’ libraries.

Here we define cd as the component data of ethylene.

To reuse the parameters from previous model like P, V, Z etc we need to define them here using the IS_A operator

We define tpeng as a pengrobinson.

Now tpeng includes all the equations and methods defined in pengrobinson model.

We have passed the required parameters to pengrobinson

So we do not have to write the equations and variables again.

This is how we make our model reusable.

We define equations for Tr and Pr as :

Tr = T/cd.data[c1].Tc;

Pr = P/cd.data[c1].Pc;

cd.dada[c1].Tc imports the value of Tc for ethylene from thermodynamics library.

In METHODS we define the standard methods again

But now under METHOD specify we RUN tpeng.specify

Under METHOD values we RUN tpeng.values. This is done to run the METHODS from pengrobinson model in our test model.

In the end we define the METHOD on_load in which we run specify and values;

Now end the test model

Under Modules section, you will see two names

peng robinsing and test model

In the panel below you will see the message

'won’t auto-initiate the model peng robinson

model requires parameters'

As we have used WILL_BE in peng robinson, we can’t execute this model but the parameters defined here will be executed in our test model

Now double click on test, the model opens under simulation

Now run the on_load method and click on solve

You will see pressure, temperature are the fixed variables and ASCEND has successfully calculated values for z, volume, reduced Temperature and reduced Pressure

Summary

• Create a reusable model
• Reuse the model for different components
• Import data from the Thermodynamics library
• Solve the model

Assignment

Also re-use this model in a test model to compute value of V for ethylene.

About the Spoken Tutorial Project

spoken-tutorial.org

Spoken Tutorial Workshops

Acknowledgements

Thank You

Thank You for joining