source: branches/newlanguage/sample/miscellaneous/sample_switcher.mso @ 141

Last change on this file since 141 was 141, checked in by Rafael de Pelegrini Soares, 15 years ago

Updated the switch sample, now working fine
  • Property svn:eol-style set to native
  • Property svn:keywords set to Id
File size: 1.7 KB
Line 
1
2using "types";
3
4Model BallWithSafetyValve
5        ATTRIBUTES
6        Info = "A fixed volume ball with a safety valve
7        for the pressure.";
8       
9    PARAMETERS
10    ValveMode as Switcher(Valid = ["closed", "opened"], Default = "closed");
11        R as Real(Unit="J/(K*mol)", Default=8.314472);
12        V as volume;
13        T as temperature;
14        Pout as pressure;
15        k as Real(Default=1e-3);
16
17        VARIABLES
18        flow_in as flow_mol;
19        flow_out as flow_mol;
20        P as pressure(Default=3);
21        n as mol;
22       
23    EQUATIONS
24
25        "Ideal gas relation"
26        P*V = n*R*T;
27       
28        "Mass balance"
29        diff(n) = flow_in - flow_out;
30       
31    switch ValveMode
32        case "closed":
33                flow_out = 0 * "mol/h";
34                when P > 2*"atm" switchto "opened";
35
36        case "opened":
37                flow_out = k*sqrt(P - Pout);
38                when P < Pout * 1.1 switchto "closed";
39        end
40end
41
42
43Model BallWithSafetyValve2
44        ATTRIBUTES
45        Info = "A fixed volume ball with a safety valve
46        for the pressure.";
47       
48    PARAMETERS
49    ValveMode as Switcher(Valid = ["closed", "opened"], Default = "closed");
50        R as Real(Unit="J/(K*mol)", Default=8.314472);
51        V as volume;
52        T as temperature;
53        Pout as pressure;
54        k as Real(Default=1e-4);
55
56        VARIABLES
57        flow_in as flow_mol;
58        flow_out as flow_mol;
59        P as pressure(Default=3);
60        n as mol;
61       
62    EQUATIONS
63
64        "Ideal gas relation"
65        P*V = n*R*T;
66       
67        "Mass balance"
68        diff(n) = flow_in - flow_out;
69       
70    #switch ValveMode
71        #case "closed":
72        #       flow_out = 0 * "mol/h";
73        #       when P > 2*"atm" switchto "opened";
74
75        #case "opened":
76                flow_out = k*sqrt(P - Pout);
77        #       when P < Pout * 1.1 switchto "closed";
78        #end
79end
80
81FlowSheet Test
82        DEVICES
83        ball as BallWithSafetyValve;
84       
85        SET
86        ball.V = 1 * "l";
87        ball.T = 300 * "K";
88        ball.Pout = 1 * "atm";
89        ball.ValveMode = "opened";
90       
91        SPECIFY
92        ball.flow_in = 0.1 * "kmol/h";
93       
94        INITIAL
95        ball.P = 1 * "atm";
96       
97        OPTIONS
98        TimeStep = 0.01;
99        TimeEnd = 10;
100end
Note: See TracBrowser for help on using the repository browser.