# Diff of /trunk/models/johnpye/absorber.a4c

revision 319 by johnpye, Fri Feb 10 04:19:11 2006 UTC revision 320 by johnpye, Thu Feb 23 13:38:40 2006 UTC
# Line 1  Line 1
1  REQUIRE "johnpye/iapws95.a4c";  REQUIRE "johnpye/iapws95.a4c";
2  REQUIRE "johnpye/iapws_sat_curves.a4c";  REQUIRE "johnpye/iapws_sat_curves.a4c";
3
4
5  MODEL absorber;  MODEL absorber;
6      (* assumptions:      (* assumptions:
7          outlet is saturated steam          outlet is saturated steam
# Line 78  METHOD values; Line 79  METHOD values;
79  END values;  END values;
80
81  END absorber;  END absorber;
82
83
84    (*
85        This model seems completely correct but it won't converge.
86        It's a problem with the S_out converging from defined (p,h).
87
88        Need to investivate
89    *)
90    MODEL absorber2;
91        S_in IS_A iapws95_2phase;
92        S_out IS_A iapws95_2phase;
93        mdot_in IS_A mass_rate;
94        mdot_out IS_A mass_rate;
95        Q IS_A energy_rate;
96        m_water IS_A mass;
97        V_total IS_A volume;
98
99        H_in IS_A energy_rate;
100        H_in = mdot_in*S_in.h;
101        H_out IS_A energy_rate;
102        H_out = mdot_out*S_out.h;
103
104        Q = H_out - H_in;
105
106        mdot_out = mdot_in;
107        S_out.p = S_in.p;
108
109        x_avg IS_A fraction;
110        alpha IS_A fraction;
111        x_avg = (S_in.x + S_out.x) / 2;
112        alpha * S_out.rho * (1-x_avg) = 1000{kg/m^3} * x_avg * (1-alpha);
113
114        m_water = S_in.rhol * (1-alpha) * V_total;
115
116    METHODS
117    METHOD default_self;
118        RUN reset; RUN values;
119        RUN scale_self;
120    END default_self;
121
122    METHOD scale_self;
123        S_out.Sl.rho.nominal := 800 {kg/m^3};
124    END scale_self;
125
126    METHOD specify;
127        FIX V_total, Q;
128        FIX mdot_in, S_in.T, S_in.rho;
129    END specify;
130
131    METHOD values;
132        V_total := 300{m} * 16 * 1{PI}*( 40{mm} )^2;
133        mdot_in := 0.4 {kg/s};
134        Q := 800 {W/m^2} * 27(*concentration*) * 500{mm} * 60{m};
135        S_in.T := 175 {K} + 273.15 {K};
136        S_in.rho := 892 {kg/m^3};
137
138        (* free *)
139        S_out.T := S_in.T;
140        S_out.rho := S_in.rho;
141    END values;
142
143    END absorber2;

Legend:
 Removed from v.319 changed lines Added in v.320