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 The partial derivative of a variable x with respect to time t can be thus written with respect to the angular position θ following: This yields the following system: ∂x ∂t π ∂x = τ ∂θ ro Mass conservation equation M d π ∂W + m τ ∂θ ro a ⎛ ⎜ ⎝ π τ ro Mass transfer equation: M d π ∂W = h τ ∂θ ro m S ∂wa ∂w + u ∂θ ∂z ( w − w ) Heat conservation equation a eq a ⎞ ⎟ = 0 ⎠ (31) (32) M d π ∂H + m τ ∂θ ro a ⎛ ⎜ ⎝ π τ ro Heat transfer equation M d π ∂H = h τ ∂θ ro m S ∂ha ∂ha + u ∂θ ∂z ⎞ ⎟ = 0 ⎠ ( w − w )( h + c T ) + h S( T − T ) a eq fg pv a t The obtained equation system is coupled non linear hyperbolic, since the enthalpy of the moist air and of the desiccant depend on the water content and on the temperature. In the bibliography the equations were solved by Macalaine-cross, [17] and Stabat, [11] by neglecting the terms ∂w a ∂h and a ∂θ ∂θ a d (33) (34) in front of transport term which yields a simplification of the equations and of the matrix to be solved. In this work we conserved these terms. Transformation to a discrete problem and restriction to 2 dimensions The wheel is divided into elementary domains following the figure below R ∆z ∆θ Figure 20 : discretisation of the desiccant wheel page 48

IEA SHC Task 38 Solar Air Conditioning and Refrigeration Subtask C2-A, November 9, 2009 Since in equations the radius R of the wheel does not appear since the variables are assumed to be constant in the radial direction the problem can thus be treated as two dimensional problem following z and θ. The considered steps are then: 2π ∆θ = and Na ∆ z = L Np Where Na is the number of steps in the angular position θ direction and Np is the number of steps in the width z direction. The domain can be represented on Cartesian scale: z=0 z=L θ=0 i=1 j=1 j=Np j-1 j j+1 rotation de la roue θ=π i=Na/2 air process i-1 i i+1 ∆θ 2∆θ θ=2π Roue complète air régénération i=Na j=Np j=1 z=L z=0 ∆z 2∆z Process ou régénération Figure 21: Cartesian representation of the grid For each node in the center of the elementary domain is identified by the coordinate (i,j) where i corresponds to the angular position and j to the width position: For the process: Na 0 < θ < π ⇔ i = 1 to 2 For the regeneration: ⎧ j = 1 ⎨ ⎩ j = Np Na π < θ < 2π ⇔ i = + 1 to Na 2 process process ⎧ j = 1 ⎨ ⎩ j = Np inlet outlet regeneration regeneration ( z = 0) ( z = L) inlet outlet ( z = 0) ( z = L) Estimation of the partial derivatives page 49

<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 />

The partial deriv<strong>at</strong>ive of a variable x with respect to time t can be thus written with respect to<br />

the angular position θ following:<br />

This yields the following system:<br />

∂x<br />

∂t<br />

π ∂x<br />

=<br />

τ ∂θ<br />

ro<br />

Mass conserv<strong>at</strong>ion equ<strong>at</strong>ion<br />

M<br />

d<br />

π ∂W<br />

+ m<br />

τ ∂θ<br />

ro<br />

a<br />

⎛<br />

⎜<br />

⎝<br />

π<br />

τ<br />

ro<br />

Mass transfer equ<strong>at</strong>ion:<br />

M<br />

d<br />

π ∂W<br />

= h<br />

τ ∂θ<br />

ro<br />

m<br />

S<br />

∂wa<br />

∂w<br />

+ u<br />

∂θ<br />

∂z<br />

( w − w )<br />

He<strong>at</strong> conserv<strong>at</strong>ion equ<strong>at</strong>ion<br />

a<br />

eq<br />

a<br />

⎞<br />

⎟ = 0<br />

⎠<br />

(31)<br />

(32)<br />

M<br />

d<br />

π ∂H<br />

+ m<br />

τ ∂θ<br />

ro<br />

a<br />

⎛<br />

⎜<br />

⎝<br />

π<br />

τ<br />

ro<br />

He<strong>at</strong> transfer equ<strong>at</strong>ion<br />

M<br />

d<br />

π ∂H<br />

= h<br />

τ ∂θ<br />

ro<br />

m<br />

S<br />

∂ha<br />

∂ha<br />

+ u<br />

∂θ<br />

∂z<br />

⎞<br />

⎟ = 0<br />

⎠<br />

( w − w )( h + c T ) + h S( T − T )<br />

a<br />

eq<br />

fg<br />

pv<br />

a<br />

t<br />

The obtained equ<strong>at</strong>ion system is coupled non linear hyperbolic, since the enthalpy of the<br />

moist air <strong>and</strong> of the desiccant depend on the w<strong>at</strong>er content <strong>and</strong> on the temper<strong>at</strong>ure.<br />

In the bibliography the equ<strong>at</strong>ions were solved by Macalaine-cross, [17] <strong>and</strong> Stab<strong>at</strong>, [11] by<br />

neglecting the terms<br />

∂w a<br />

∂h<br />

<strong>and</strong> a<br />

∂θ ∂θ<br />

a<br />

d<br />

(33)<br />

(34)<br />

in front of transport term which yields a simplific<strong>at</strong>ion of<br />

the equ<strong>at</strong>ions <strong>and</strong> of the m<strong>at</strong>rix to be solved. In this work we conserved these terms.<br />

Transform<strong>at</strong>ion to a discrete problem <strong>and</strong> restriction to 2 dimensions<br />

The wheel is divided into elementary domains following the figure below<br />

R<br />

∆z<br />

∆θ<br />

Figure 20 : discretis<strong>at</strong>ion of the desiccant wheel<br />

page 48

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

Saved successfully!

Ooh no, something went wrong!