/[ascend]/branches/ksenija2/models/johnpye/solardynamics.a4c
ViewVC logotype

Diff of /branches/ksenija2/models/johnpye/solardynamics.a4c

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 2858 by jpye, Sat Mar 21 13:13:27 2015 UTC revision 2859 by jpye, Sun Mar 22 01:40:44 2015 UTC
# Line 35  MODEL solardynamics; Line 35  MODEL solardynamics;
35      *)      *)
36    
37      (* Solar collector *)      (* Solar collector *)
38      C IS_A factor;      (*C IS_A factor;
39      A IS_A area;      A IS_A area;*)
40      Qdot_abs IS_A energy_rate;      Qdot_abs IS_A energy_rate;
41      col_on: Qdot_abs = C*A*G;      col_on: Qdot_abs = 800.*100.*G;
42      col_off: Qdot_abs = 0;      col_off: Qdot_abs = 0;
43    
44      (* Storage *)      (* Storage *)
# Line 48  MODEL solardynamics; Line 48  MODEL solardynamics;
48      T ALIASES TA.T;      T ALIASES TA.T;
49      DERIVATIVE OF h;      DERIVATIVE OF h;
50      Qdot_elec IS_A energy_rate;      Qdot_elec IS_A energy_rate;
51      m * der(h) = Qdot_abs - Qdot_elec;      en_eq: m * der(h) = Qdot_abs - Qdot_elec;
52    
53      (* Power block *)      (* Power block *)
54      Wdot_elec, Qdot_elec_des IS_A energy_rate;      Wdot_elec, Qdot_elec_des IS_A energy_rate;
55      eta IS_A fraction;      eta IS_A fraction;
56      (* power cycle efficiency: chambadal novikov efficiency *)      (* power cycle efficiency: chambadal novikov efficiency *)
57      eta = 1 - sqrt(T_amb / T);      eta_eq: eta = 1 - sqrt(T_amb / T);
58      Qdot_elec*eta = Wdot_elec;      Wdot_eq: Qdot_elec*eta = Wdot_elec;
59      pow_on: Qdot_elec = Qdot_elec_des;      pow_on: Qdot_elec = Qdot_elec_des;
60      pow_off: Qdot_elec = 0;      pow_off: Qdot_elec = 0;
61    
# Line 91  METHODS Line 91  METHODS
91          Wdot_elec.obs_id := 4;          Wdot_elec.obs_id := 4;
92      END obs_init;      END obs_init;
93      METHOD specify;      METHOD specify;
94          FIX C := 2000;          (*FIX C := 2000;*)
95          FIX Qdot_elec_des := 1 {MW};          FIX Qdot_elec_des := 1 {MW};
96          FIX A := 2000 {m^2};          (*FIX A := 2000 {m^2};*)
97          FIX TA.p := 1 {bar};          FIX TA.p := 1 {bar};
98            FIX m := 20 {t};
99          t := 0 {s};          t := 0 {s};
100          T := 300 {K} + 273.15 {K};          T := 300 {K} + 273.15 {K};
101          is_G_high := FALSE;          is_G_high := FALSE;
# Line 102  METHODS Line 103  METHODS
103          is_T_start := FALSE;          is_T_start := FALSE;
104          is_T_stop := TRUE;          is_T_stop := TRUE;
105      END specify;      END specify;
106        METHOD nla_init;
107            FIX T;
108            FIX t := 20.5 {d};
109        END nla_init;
110      METHOD on_load;      METHOD on_load;
111          RUN specify;          RUN specify;
112          RUN obs_init;          RUN obs_init;

Legend:
Removed from v.2858  
changed lines
  Added in v.2859

john.pye@anu.edu.au
ViewVC Help
Powered by ViewVC 1.1.22