IEA Solar Heating and Cooling Programm - NachhaltigWirtschaften.at

IEA Solar Heating and Cooling Programm - NachhaltigWirtschaften.at IEA Solar Heating and Cooling Programm - NachhaltigWirtschaften.at

nachhaltigwirtschaften.at
from nachhaltigwirtschaften.at More from this publisher
04.11.2013 Views

IEA SHC Task 38 Solar Air Conditioning and Refrigeration Subtask C2-A, November 9, 2009 As stated in equation (14) above, the mass flow of weak solution is set constant. In equation (12) it is assumed that the mass flow entering the absorber tube bundle at time interval i equals the mass flow which has left the generator c 1 time intervals ago. Doing so, the physical incompressibility of the solution in the piping between generator and absorber is not modelled correctly. In reality, the mass flow entering the tube at time interval i should equal the mass flow leaving the tube at time interval i. The assumption, however, had to be made in order to achieve a correct salt balance in equation (22). The result is a hidden mass storage term in the solution tube. However, the amount of stored solution in the tube is only 0.6 % of the total amount of stored solution in the sumps [3]. The error introduced by this assumption is therefore negligible. Mass balances The vessel sumps (see Figure 2) are assumed to be fully mixed at each time interval. Thus vessel sump and solution leaving the vessel are assumed to have the same salt concentration, x sA,i and x wG,i , respectively. The outlet of each tube bundle is assumed to exhibit the equilibrium concentration x sG,i and x wA,i . The total solution mass stored in generator and absorber sump at time i is expressed by M sol st, G, i M , , and sol, st, A i , respectively. We know that the total mass (or salt mass, respectively) in the sump at time interval i equals the total mass (or salt mass) at time interval (i-1) plus the difference of ingoing and outgoing solution (or salt) flows at time interval i. The balances for the total solution contents in the sumps are given in equations (15) and (16). ∆ t is the time between two simulation intervals. M sol, , G, i − M sol, , G, i−1 Generator: m& , , , − m& , , − st st sol tb G i sol sG i = 0 (15) ∆t M sol, st, A, i − M sol, st, A, i−1 Absorber: m& , , , − m& sol tb A i sol, wA, i − = 0 (16) ∆t The fraction term in equations (15 and 16) is the amount of solution which is added to the solution stored in either generator or absorber sump during the time simulation intervals i and (i-1). ∆ t between consecutive The salt flow balance, analogously, in both sumps can be written as: Generator M sol, , G, i ⋅ xsA, i − M sol, , G, i−1 ⋅ xsA, i−1 m& , , , ⋅ x , − m& , , ⋅ x , − st st = 0 (17 sol tb G i sG i sol sG i sA i ∆t : ) Absorber : M sol, st, A, i ⋅ xwG , i − M sol, st, A, i−1 ⋅ xwG , i−1 m& , , , ⋅ x , − m& sol tb A i wA i sol, wA, i ⋅ xwG , i − = 0 (18 ∆t ) page 68

IEA SHC Task 38 Solar Air Conditioning and Refrigeration Subtask C2-A, November 9, 2009 Proceeding from the sumps to the tube bundles, the mass flow balance yields Generator: m & − m& − m& 0 (19) sol, wA, i−c2 vG, i sol, tb, G, i = Absorber: m & + m& − m& 0 (20) sol, sG, i−c1 vA, i sol, tb, A, i = The salt flow balances of the tube bundles read Generator: m & ⋅ x − m& ⋅ x 0 (21) sol, wA, i−c2 wG , i−c2 sol, tb, G, i sG, i = Absorber: m & ⋅ x − m& ⋅ x 0 (22) sol, sG, i−c1 sA, i−c1 sol, tb, A, i wA, i = Pressure drop The strong solution flow from generator to absorber is driven by gravity and a pressure difference. The pressure loss of the solution heat exchanger and its adjacent piping is assumed constant which is a good approximation as long as the flow variations are not too large; hence the mass flow of the strong solution depends only on the pressure difference between generator/absorber and the liquid solution column at heat exchanger inlet. This assumption, of course, can be refined in future works. During operation, the static pressure of the liquid column has to equal the dynamic pressure loss caused by flow restrictions through piping and solution heat exchanger. The mass flow of the strong solution can therefore be calculated as ( p − p + ρ ⋅ g ⋅ ( h + z )) 2 ⋅ ρsol, s ⋅ G, i A, i sol, s i m& sol, s, i = At ⋅ , (23) ς where generator/condenser and evaporator/absorber pressure are calculated using the equation of Clausius-Clapeyron as shown in equation (24). There, the pressure and temperature values of solution interval (i-1) are being used to calculate the values of simulation interval i. This is a valid approximation as long as the time between two simulation intervals, ∆ t , is not too long. ⎛ p X , i ⎞ r ⎛ ⎞ ⎜ ⎟ 0 = ⋅⎜ 1 1 ln ⎟ − (24) ⎝ p X , i−1 ⎠ R ⎝ TX , i−1 TX , i ⎠ Thermal storage page 69

<strong>IEA</strong> SHC Task 38 <strong>Solar</strong> Air Conditioning <strong>and</strong> Refriger<strong>at</strong>ion Subtask C2-A, November 9, 2009<br />

Proceeding from the sumps to the tube bundles, the mass flow balance yields<br />

Gener<strong>at</strong>or: m & − m&<br />

− m&<br />

0<br />

(19)<br />

sol, wA,<br />

i−c2<br />

vG,<br />

i sol,<br />

tb,<br />

G,<br />

i<br />

=<br />

Absorber: m & + m&<br />

− m&<br />

0<br />

(20)<br />

sol, sG,<br />

i−c1<br />

vA,<br />

i sol,<br />

tb,<br />

A,<br />

i<br />

=<br />

The salt flow balances of the tube bundles read<br />

Gener<strong>at</strong>or: m & ⋅ x − m&<br />

⋅ x 0<br />

(21)<br />

sol, wA,<br />

i−c2<br />

wG , i−c2<br />

sol,<br />

tb,<br />

G,<br />

i sG,<br />

i<br />

=<br />

Absorber: m & ⋅ x − m&<br />

⋅ x 0<br />

(22)<br />

sol, sG,<br />

i−c1<br />

sA,<br />

i−c1<br />

sol,<br />

tb,<br />

A,<br />

i wA,<br />

i<br />

=<br />

Pressure drop<br />

The strong solution flow from gener<strong>at</strong>or to absorber is driven by gravity <strong>and</strong> a pressure<br />

difference. The pressure loss of the solution he<strong>at</strong> exchanger <strong>and</strong> its adjacent piping is<br />

assumed constant which is a good approxim<strong>at</strong>ion as long as the flow vari<strong>at</strong>ions are not too<br />

large; hence the mass flow of the strong solution depends only on the pressure difference<br />

between gener<strong>at</strong>or/absorber <strong>and</strong> the liquid solution column <strong>at</strong> he<strong>at</strong> exchanger inlet. This<br />

assumption, of course, can be refined in future works. During oper<strong>at</strong>ion, the st<strong>at</strong>ic pressure<br />

of the liquid column has to equal the dynamic pressure loss caused by flow restrictions<br />

through piping <strong>and</strong> solution he<strong>at</strong> exchanger. The mass flow of the strong solution can<br />

therefore be calcul<strong>at</strong>ed as<br />

( p − p + ρ ⋅ g ⋅ ( h + z ))<br />

2 ⋅ ρsol,<br />

s<br />

⋅<br />

G,<br />

i A,<br />

i sol,<br />

s<br />

i<br />

m& sol,<br />

s,<br />

i<br />

= At<br />

⋅<br />

, (23)<br />

ς<br />

where gener<strong>at</strong>or/condenser <strong>and</strong> evapor<strong>at</strong>or/absorber pressure are calcul<strong>at</strong>ed using the<br />

equ<strong>at</strong>ion of Clausius-Clapeyron as shown in equ<strong>at</strong>ion (24). There, the pressure <strong>and</strong><br />

temper<strong>at</strong>ure values of solution interval (i-1) are being used to calcul<strong>at</strong>e the values of<br />

simul<strong>at</strong>ion interval i. This is a valid approxim<strong>at</strong>ion as long as the time between two simul<strong>at</strong>ion<br />

intervals,<br />

∆ t , is not too long.<br />

⎛ p<br />

X , i<br />

⎞ r ⎛<br />

⎞<br />

⎜ ⎟<br />

0<br />

= ⋅⎜<br />

1 1<br />

ln ⎟<br />

−<br />

(24)<br />

⎝ p<br />

X , i−1<br />

⎠ R ⎝ TX<br />

, i−1<br />

TX<br />

, i ⎠<br />

Thermal storage<br />

page 69

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!