source: trunk/sample/auto/cstr_dae.mso @ 388

Last change on this file since 388 was 388, checked in by Argimiro Resende Secchi, 15 years ago

Adding AUTO DAE examples.

File size: 1.2 KB
Line 
1using "types";
2
3FlowSheet cstr2_dae
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;
21
22        VARIABLES
23        u1 as Real;
24        u2 as Real;
25        u3 as Real;
26        u4 as Real;
27        u5 as Real;
28        u6 as Real;
29
30        SET
31        fi              = 0.04;
32
33        beta    = 8;
34
35        delta   = 0.78;                                         
36
37        q               = 1;
38
39        alfa    = 0.19;
40
41        S               = 1.015;
42
43        csi             = 0.32;
44
45        delta1  = 10;
46
47        delta2  = 0.952;
48
49        gama    = 27.85;                       
50
51        x1f             = 1;                   
52
53        x2f             = 0;
54
55        x3f             = 0;
56
57        x4f             = -1;
58
59        qc              = 0;
60
61        EQUATIONS
62 
63
64        q*(x1f-u1)-u1*u5*fi - diff(u1)*'s' = 0;
65
66        q*(x2f-u2)-u2*fi*S*u6+u1*fi*u5 - diff(u2)*'s' = 0;
67
68        q*(x3f-u3)+delta*(u4-u3)+beta*fi*(u1*u5+alfa*u2*u6*S) - diff(u3)*'s' = 0;
69
70        delta1*(qc*(x4f-u4)+delta*delta2*(u3-u4)) - diff(u4)*'s' = 0;
71
72        -u5 + exp(u3/(1+(u3/gama))) = 0;
73
74        -u6 + exp(csi*u3/(1+(u3/gama))) = 0;
75
76        GUESS
77        u1 = 0.6419;
78
79        u2 = 0.3270;
80
81        u3 = 2.9101;
82
83        u4 = 2.9101;
84
85        u5 = 13.9403;
86
87        u6 = 2.3236;                                     
88
89
90    INITIAL
91        u1 = 0.6;
92
93        u2 = 0.3;
94
95        u3 = 3;
96
97        u4 = 3;
98
99
100        OPTIONS
101        TimeEnd = 10;
102        TimeStep = 0.2;
103        Dynamic = false;
104end
Note: See TracBrowser for help on using the repository browser.