1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
|
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
|
+
+
+
|
# -*-octave-*- Put Emacs into octave-mode
# State specification (Reactor_state.txt)
# Generated by MTT at Fri Mar 3 11:52:23 GMT 2000
###############################################################
## Version control history
###############################################################
## $Id$
## $Log$
## Revision 1.1 2000/12/28 17:09:55 peterg
## To RCS
##
###############################################################
## Reduce steady-state parameter file (Reactor_sspar.r)
## as siso_sspar ecxept that states/states have different meaning
## Steady state for constant c_a, c_b and t=t_s and f=f_s
## Unit volume Reactor:
|
45
46
47
48
49
50
51
52
53
54
55
56
57
58
|
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
|
-
-
-
+
+
+
+
+
+
|
q_s = -((t_0-t_s)*c_p*f_s + e_1*h_1*k_1*x1 + e_2*h_2*k_2*x2 + e_3*h_3*k_3*x1^2);
## The two inputs at steady-state
u1 = f_s;
u2 = q_s;
## Load up the states
mttx(1) = x1;
mttx(2) = x2;
mttx(3) = x3;
## Removed by MTT on Thu Dec 28 18:46:20 GMT 2000: mttx(1) = x1;
## Removed by MTT on Thu Dec 28 18:46:20 GMT 2000: mttx(2) = x2;
## Removed by MTT on Thu Dec 28 18:46:20 GMT 2000: mttx(3) = x3;
reactor_h_r = x3; # Added by MTT on Thu Dec 28 18:46:25 GMT 2000
reactor_m_a = x1; # Added by MTT on Thu Dec 28 18:46:25 GMT 2000
reactor_m_b = x2; # Added by MTT on Thu Dec 28 18:46:25 GMT 2000
|