/[ascend]/trunk/models/johnpye/fprops/brayton_fprops.a4c
ViewVC logotype

Diff of /trunk/models/johnpye/fprops/brayton_fprops.a4c

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

revision 2673 by jpye, Wed Dec 12 12:39:25 2012 UTC revision 2674 by jpye, Thu Jan 24 02:27:41 2013 UTC
# Line 26  REQUIRE "johnpye/airprops.a4c"; Line 26  REQUIRE "johnpye/airprops.a4c";
26  REQUIRE "johnpye/fprops/rankine_fprops.a4c";  REQUIRE "johnpye/fprops/rankine_fprops.a4c";
27  REQUIRE "johnpye/fprops/heatex_pinch.a4c";  REQUIRE "johnpye/fprops/heatex_pinch.a4c";
28    
29    IMPORT "sensitivity/solve";
30  IMPORT "johnpye/extpy/extpy";  IMPORT "johnpye/extpy/extpy";
31  IMPORT "johnpye/fprops/cycle_plot_rachel";  IMPORT "johnpye/fprops/cycle_plot_rachel";
32    
# Line 146  END brayton_common; Line 147  END brayton_common;
147    
148  MODEL brayton_co2 REFINES brayton_common;  MODEL brayton_co2 REFINES brayton_common;
149      cd.component :== 'carbondioxide';      cd.component :== 'carbondioxide';
150        cd.type :== 'helmholtz';
151  METHODS  METHODS
152  METHOD on_load;  METHOD on_load;
153      RUN ClearAll;      RUN ClearAll;
# Line 209  END regenerator; Line 211  END regenerator;
211    
212  MODEL regenerator_test REFINES regenerator;  MODEL regenerator_test REFINES regenerator;
213      inlet.cd.component :== 'carbondioxide';      inlet.cd.component :== 'carbondioxide';
214      n :== 8;      inlet.cd.type :== 'helmholtz';
215      inlet_hot.cd, inlet.cd ARE_THE_SAME;      inlet_hot.cd, inlet.cd ARE_THE_SAME;
216        n :== 8;
217  METHODS  METHODS
218  METHOD on_load;  METHOD on_load;
219        RUN heatex_pinch::default_self;
220      RUN bound_self;      RUN bound_self;
221      FIX inlet.mdot, inlet.p, inlet.T;      RUN heatex_pinch::default_self;
222      FIX inlet_hot.mdot, inlet_hot.p, inlet_hot.T;      (* heat transfer for a particular value of temperature drop at the pinch *)
223      inlet.mdot := 1 {kg/s};      FIX inlet.p := 200 {bar};
224      inlet.p := 10 {bar};      FIX inlet_hot.p := 70 {bar}; (* doesn't currently solve for CO2 if the pressure is too low... needed to figure out why *)
225      inlet.T := 300 {K};      FIX inlet_hot.T := 100 {K} + 273.15 {K};
226      inlet_hot.mdot := 1.05 {kg/s};      FIX inlet.T := 31 {K} + 273.15 {K};
227      inlet_hot.p := 15 {bar};      FIX mdot_hot := 1.1 {kg/s};
228      inlet_hot.T := 500 {K};      FIX mdot := 0.81 {kg/s};
229      FIX epsilon;      FIX DT_pinch := 10 {K};
     epsilon := 1;  
230  END on_load;  END on_load;
231  END regenerator_test;  END regenerator_test;
232    
# Line 248  MODEL brayton_regen; Line 251  MODEL brayton_regen;
251      RE IS_A heatex_pinch;      RE IS_A heatex_pinch;
252      cd ALIASES BO.cd;      cd ALIASES BO.cd;
253      cd.component :== 'carbondioxide';      cd.component :== 'carbondioxide';
254        cd.type :== 'helmholtz';
255      RE.n :== 8;      RE.n :== 8;
256            
257      PU.outlet, RE.inlet ARE_THE_SAME;      PU.outlet, RE.inlet ARE_THE_SAME;
# Line 286  MODEL brayton_regen; Line 290  MODEL brayton_regen;
290      eta_PU ALIASES PU.eta;      eta_PU ALIASES PU.eta;
291      DT_RE ALIASES RE.DT_pinch;      DT_RE ALIASES RE.DT_pinch;
292  METHODS  METHODS
293    METHOD default_self;
294        RUN BO.default_self;
295        RUN TU.default_self;
296        RUN CO.default_self;
297        RUN PU.default_self;
298        RUN RE.default_self;
299    END default_self;
300  METHOD on_load;  METHOD on_load;
301      RUN RE.bound_self;      RUN default_self;
302      FIX PU.eta, TU.eta;      FIX PU.eta, TU.eta;
303      PU.eta := 0.8;      PU.eta := 0.8;
304      TU.eta := 0.85;      TU.eta := 0.85;
# Line 349  MODEL brayton_regen_reheat_intercool; Line 360  MODEL brayton_regen_reheat_intercool;
360      eta IS_A fraction;      eta IS_A fraction;
361      eta = Wdot / Qdot_H;      eta = Wdot / Qdot_H;
362  METHODS  METHODS
363    METHOD default_self;
364        RUN BU1.default_self;
365        RUN BU2.default_self;
366        RUN TU1.default_self;
367        RUN TU2.default_self;
368        RUN CO1.default_self;
369        RUN CO2.default_self;
370        RUN RE.default_self;
371        RUN DI1.default_self;
372        RUN DI2.default_self;
373    END default_self;
374  METHOD on_load;  METHOD on_load;
375      RUN RE.bound_self;      RUN RE.bound_self;
376        RUN default_self;
377      FIX BU1.eta := 1;      FIX BU1.eta := 1;
378      FIX BU2.eta := 1;      FIX BU2.eta := 1;
379      FIX TU1.eta := 0.85;      FIX TU1.eta := 0.85;
# Line 368  END brayton_regen_reheat_intercool; Line 391  END brayton_regen_reheat_intercool;
391    
392  MODEL brayton_regen_reheat_intercool_co2 REFINES brayton_regen_reheat_intercool;  MODEL brayton_regen_reheat_intercool_co2 REFINES brayton_regen_reheat_intercool;
393      cd.component :== 'carbondioxide';      cd.component :== 'carbondioxide';
394        cd.type :== 'helmholtz';
395      RE.n :== 8;      RE.n :== 8;
396  METHODS  METHODS
397  METHOD on_load;  METHOD on_load;
# Line 378  METHOD on_load; Line 402  METHOD on_load;
402      FIX p_2a := 150 {bar};      FIX p_2a := 150 {bar};
403      FIX p_2b := p_2a;      FIX p_2b := p_2a;
404      FIX p_3 := 300 {bar};      FIX p_3 := 300 {bar};
405      FIX RE.epsilon := 0.77;      (* FIX RE.epsilon := 0.77; *)
406        FIX RE.DT_pinch := 15 {K};
407  END on_load;  END on_load;
408  END brayton_regen_reheat_intercool_co2;  END brayton_regen_reheat_intercool_co2;
409    

Legend:
Removed from v.2673  
changed lines
  Added in v.2674

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