source: trunk/sample/auto/cstr_ode.mso @ 984

Last change on this file since 984 was 984, checked in by Argimiro Resende Secchi, 6 years ago
File size: 1.2 KB
Line 
1using "types";
2
3FlowSheet cstr_ode
4
5        PARAMETERS
6        fi as Real;
7        beta as Real;
8        delta as Real;
9        q as Real;
10        alfa as Real;
11        S as Real;
12        csi as Real;
13        delta1 as Real;
14        delta2 as Real;
15        gama as Real;
16        x1f as Real;
17        x2f as Real;
18        x3f as Real;
19        x4f as Real;
20        qc as Real (Lower=-1,Upper=20);
21
22        VARIABLES
23        u1 as Real (Lower=0,Upper=1);
24        u2 as Real (Lower=0,Upper=1);
25        u3 as Real (Lower=0,Upper=20);
26        u4 as Real (Lower=0,Upper=20);
27
28        SET
29
30        fi              = 0.133;
31        beta    = 8;
32        delta   = 1;                                           
33        q               = 1;
34        alfa    = 1;
35        S               = 0.01;
36        csi             = 1;
37        delta1  = 10;
38        delta2  = 1;
39        gama    = 1000;                 
40        x1f             = 1;   
41        x2f             = 0;
42        x3f             = 0;
43        x4f             = -1;
44        qc              = 0;
45       
46        EQUATIONS
47 
48
49        q*(x1f-u1)-u1*exp(u3/(1+(u3/gama)))*fi - diff(u1)*'s' = 0;
50
51        q*(x2f-u2)-u2*fi*S*exp(csi*u3/(1+(u3/gama)))+u1*fi*exp(u3/(1+(u3/gama))) - diff(u2)*'s' = 0;
52
53        q*(x3f-u3)+delta*(u4-u3)+beta*fi*(u1*exp(u3/(1+(u3/gama)))+alfa*u2*exp(csi*u3/(1+(u3/gama)))*S) - diff(u3)*'s' = 0;
54
55        delta1*(qc*(x4f-u4)+delta*delta2*(u3-u4)) - diff(u4)*'s' = 0;
56
57        GUESS
58        u1 = 0;
59        u2 = 0;
60        u3 = 16;
61        u4 = 16;
62
63
64    INITIAL
65        u1 = 0.6;
66        u2 = 0.3;
67        u3 = 3;
68        u4 = 3;
69
70
71        OPTIONS
72        TimeEnd = 10;
73        TimeStep = 0.2;
74        Dynamic = false;
75end
Note: See TracBrowser for help on using the repository browser.