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

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

Adding AUTO DAE examples.

File size: 1.1 KB
Line 
1using "types";
2
3FlowSheet cstr_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        beta    = 8;
33        delta   = 0.78;                                         
34        q               = 1;
35        alfa    = 0.19;
36        S               = 1.015;
37        csi             = 0.32;
38        delta1  = 10;
39        delta2  = 0.952;
40        gama    = 27.85;                       
41        x1f             = 1;                   
42        x2f             = 0;
43        x3f             = 0;
44        x4f             = -1;
45        qc              = 0;
46
47        EQUATIONS
48 
49        q*(x1f-u1)-u1*u5*fi - diff(u1)*'s' = 0;
50        q*(x2f-u2)-u2*fi*S*u6+u1*fi*u5 - diff(u2)*'s' = 0;
51        q*(x3f-u3)+delta*(u4-u3)+beta*fi*(u1*u5+alfa*u2*u6*S) - diff(u3)*'s' = 0;
52        delta1*(qc*(x4f-u4)+delta*delta2*(u3-u4)) - diff(u4)*'s' = 0;
53        -u5 + exp(u3/(1+(u3/gama))) = 0;
54        -u6 + exp(csi*u3/(1+(u3/gama))) = 0;
55
56        GUESS
57        u1 = 0.6419;
58        u2 = 0.3270;
59        u3 = 2.9101;
60        u4 = 2.9101;
61        u5 = 13.9403;
62        u6 = 2.3236;                                     
63
64    INITIAL
65        u1 = 0.6;
66        u2 = 0.3;
67        u3 = 3;
68        u4 = 3;
69
70        OPTIONS
71        TimeEnd = 10;
72        TimeStep = 0.2;
73        Dynamic = false;
74end
Note: See TracBrowser for help on using the repository browser.