/[ascend]/trunk/models/johnpye/airprops2.a4c
ViewVC logotype

Diff of /trunk/models/johnpye/airprops2.a4c

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

revision 2801 by jpye, Fri Jul 18 03:48:18 2014 UTC revision 2802 by jpye, Fri Jul 18 08:22:02 2014 UTC
# Line 20  REQUIRE "atoms.a4l"; Line 20  REQUIRE "atoms.a4l";
20  REQUIRE "johnpye/moltensalt.a4c";  REQUIRE "johnpye/moltensalt.a4c";
21    
22  (*  (*
23      Ideal-gas air, together with some power-law correlations for      Ideal-gas air, together with some simplified correlations for
24      viscosity and conductivity of uncertain accuracy, obtained from      viscosity and thermal conductivity.
     http://dx.doi.org/10.1615/AtoZ.a.air_properties_of  
   
     THIS MODEL IS STILL IN DEVELOPMENT AND DOESN'T WORK YET!!!  
25  *)  *)
26  MODEL airprops2 REFINES thermophysical_fluid_base;  MODEL air_fluid REFINES thermophysical_fluid_base;
27      M IS_A real_constant;      M IS_A real_constant;
28      M :== 28.97{kg/kmol}; (* Cengel *)      M :== 28.97{kg/kmol}; (* Cengel *)
29      R IS_A real_constant;      R IS_A real_constant;
# Line 61  MODEL airprops2 REFINES thermophysical_f Line 58  MODEL airprops2 REFINES thermophysical_f
58      z2_eq: mu - mustar = 5.76e-3{cP} * (dpdT_v/1{atm})^1.126;      z2_eq: mu - mustar = 5.76e-3{cP} * (dpdT_v/1{atm})^1.126;
59      dpdT_v = rho * R;      dpdT_v = rho * R;
60    
61      (* conductivity *)      (* conductivity, from Kadoya et al, doi:10.1063/1.555744, using ideal EOS because of use of rho in calc *)
62      cc[-4..0,1,5] IS_A real_constant;      cc[-4..0,1,5] IS_A real_constant;
63      cc[1] :== 0.239503;      cc[1] :== 0.239503;
64      cc[5] :== 0.00649768;      cc[5] :== 0.00649768;
# Line 95  METHODS Line 92  METHODS
92          h.lower_bound := -1e5 {kJ/kg};          h.lower_bound := -1e5 {kJ/kg};
93          s.lower_bound := -1e5 {kJ/kg/K};          s.lower_bound := -1e5 {kJ/kg/K};
94          u.lower_bound := -1e5 {kJ/kg};          u.lower_bound := -1e5 {kJ/kg};
       
95      END bound_self;      END bound_self;
96      METHOD on_load;      METHOD default;
         RUN ClearAll;  
         RUN bound_self;  
97          FIX h_ref := -295.504033483 {kJ/kg};          FIX h_ref := -295.504033483 {kJ/kg};
98          FIX s_ref := 108.501816216 {kJ/kg/K};          FIX s_ref := 108.501816216 {kJ/kg/K};
99            h := 10 {kJ/kg};    
100        END default;
101        METHOD on_load;
102            RUN ClearAll;
103            RUN default;
104          FIX T := 650{K};          FIX T := 650{K};
105          FIX p := 10 {MPa};          FIX p := 10 {MPa};
106      END on_load;      END on_load;
107  END airprops2;  END air_fluid;
108    

Legend:
Removed from v.2801  
changed lines
  Added in v.2802

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