/[ascend]/branches/fprops-incomp/models/johnpye/fprops/towerarray.py
ViewVC logotype

Diff of /branches/fprops-incomp/models/johnpye/fprops/towerarray.py

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

revision 3468 by jpye, Wed Feb 12 04:03:13 2020 UTC revision 3469 by jpye, Thu Feb 13 01:45:55 2020 UTC
# Line 305  MODEL towerarray{suffix} REFINES layout{ Line 305  MODEL towerarray{suffix} REFINES layout{
305    
306    
307      (* total materials *)      (* total materials *)
308      m_pipe, m_insul IS_A mass;      m_pipe, m_insul, m_sodium IS_A mass;
309        V_insul IS_A volume;
310      rho_pipe, rho_insul IS_A mass_density;      rho_pipe, rho_insul IS_A mass_density;
311      L_pipe_tot IS_A distance;      L_pipe_tot IS_A distance;
312      m_pipe = rho_pipe * SUM[PC[i].solid_pipe.V + PH[i].solid_pipe.V | i IN [0..n_PH-1]];      m_pipe = rho_pipe * SUM[PC[i].solid_pipe.V + PH[i].solid_pipe.V | i IN [0..n_PH-1]];
313      m_insul = rho_insul * SUM[PC[i].solid_insul.V + PH[i].solid_insul.V | i IN [0..n_PH-1]];      V_insul = SUM[PC[i].solid_insul.V + PH[i].solid_insul.V | i IN [0..n_PH-1]];
314      L_pipe_tot = SUM[PC[i].L + PH[i].L | i IN [0..n_PH-1]];      L_pipe_tot = SUM[PC[i].L + PH[i].L | i IN [0..n_PH-1]];
315        m_sodium = SUM[PC[i].solid.V / PC[i].inlet.v + PH[i].solid.V / PH[i].inlet.v | i IN [0..n_PH-1]];
316        m_insul = rho_insul * V_insul;
317    
318      (* costs *)      (* costs *)
319      C_pipe, C_insul, C_inst_pipe, C_inst_insul, C_supp, C_tot IS_A monetary_unit;      C_pipe, C_insul, C_inst_pipe, C_inst_insul, C_supp, C_sodium, C_tot IS_A monetary_unit;
320      c_pipe, c_insul IS_A cost_per_mass;      c_pipe IS_A cost_per_mass;
321        c_insul IS_A cost_per_volume;
322      C_pipe = c_pipe * m_pipe;      C_pipe = c_pipe * m_pipe;
323      C_insul = c_insul * m_insul;      C_insul = c_insul * V_insul;
324        C_sodium = c_sodium * m_sodium;
325    
326      c_inst_pipe_0, c_inst_insul_0, c_supp_0 IS_A cost_per_length;      c_inst_pipe_0, c_inst_insul_0, c_supp_0 IS_A cost_per_length;
327      c_inst_pipe_1, c_inst_insul_1 IS_A cost_per_area;      c_inst_pipe_1, c_inst_insul_1 IS_A cost_per_area;
328        c_sodium IS_A cost_per_mass;
329    
330      C_inst_pipe = SUM[PC[i].L * (c_inst_pipe_1*PC[i].D_o + c_inst_pipe_0)      C_inst_pipe = SUM[PC[i].L * (c_inst_pipe_1*PC[i].D_o + c_inst_pipe_0)
331                      + PH[i].L * (c_inst_pipe_1*PH[i].D_o + c_inst_pipe_0) | i IN [0..n_PH-1]];                      + PH[i].L * (c_inst_pipe_1*PH[i].D_o + c_inst_pipe_0) | i IN [0..n_PH-1]];
# Line 381  METHOD on_load; Line 387  METHOD on_load;
387      FIX k_insul := 0.05 {W/m/K}; (* 'microporous insulation board' @ 700°C, https://is.gd/5j9Gkw *)      FIX k_insul := 0.05 {W/m/K}; (* 'microporous insulation board' @ 700°C, https://is.gd/5j9Gkw *)
388      FIX k_pipe := 22.4 {W/m/K}; (* Haynes 230 @ 700°C: https://is.gd/F0RICh *)      FIX k_pipe := 22.4 {W/m/K}; (* Haynes 230 @ 700°C: https://is.gd/F0RICh *)
389      FIX rho_pipe := 9 {g/cm^3}; (* Haynes 230 @ 700°C: https://is.gd/F0RICh *)      FIX rho_pipe := 9 {g/cm^3}; (* Haynes 230 @ 700°C: https://is.gd/F0RICh *)
390      FIX rho_insul := 380 {g/m^3}; (* 'microporous insulation board', approx value from Felix *)      FIX rho_insul := 380 {kg/m^3}; (* 'microporous insulation board', approx value from Felix *)
391    
392      FIX c_pipe := 8 {USD/kg} (* 316 stainless steel assumed *);      FIX c_pipe := 8 {USD/kg} (* 316 stainless steel assumed *);
393      FIX c_insul := 10000 {USD/kg}; (* old quote from brandname product...or maybe 250 USD/kg if alibaba.com *)  
394        FIX c_insul := 10000 {USD/m^3}; (* old quote from brandname product...or maybe 250 USD/kg if alibaba.com *)
395      FIX c_inst_pipe_0 := 22 {AUD/m};      FIX c_inst_pipe_0 := 22 {AUD/m};
396      FIX c_inst_pipe_1 := 0.3734 {AUD/m/mm};      FIX c_inst_pipe_1 := 0.3734 {AUD/m/mm};
397      FIX c_inst_insul_0 := 29 {AUD/m};      FIX c_inst_insul_0 := 29 {AUD/m};
398      FIX c_inst_insul_1 := 0.6257 {AUD/m/mm};      FIX c_inst_insul_1 := 0.6257 {AUD/m/mm};
399      FIX c_supp_0 := 100 {USD/m}; (* this one's in yankee dollars *)      FIX c_supp_0 := 100 {USD/m}; (* this one's in yankee dollars *)
400        FIX c_sodium := 3000 {USD/t}; (* based on a quick look at alibaba https://is.gd/pPKCYt *)
401    
402      (*      (*
403      installation cost for pipe, per length installed      installation cost for pipe, per length installed
# Line 425  METHOD correct_dp_and_Tout; Line 433  METHOD correct_dp_and_Tout;
433      FREE mdot_onecoll;      FREE mdot_onecoll;
434      FIX T_out := 740 {K} + 273.15{K};      FIX T_out := 740 {K} + 273.15{K};
435  END correct_dp_and_Tout;  END correct_dp_and_Tout;
436    METHOD show_temperatures;
437        EXTERNAL disharray_temperature_list(SELF);
438    END show_temperatures;
439      """ % (Qdot_onecoll_MW,mdot_onecoll));      """ % (Qdot_onecoll_MW,mdot_onecoll));
440    
441      f.write("END towerarray{suffix};\n".format(suffix=suffix));      f.write("END towerarray{suffix};\n".format(suffix=suffix));

Legend:
Removed from v.3468  
changed lines
  Added in v.3469

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