IEA Solar Heating and Cooling Programm - NachhaltigWirtschaften.at
IEA Solar Heating and Cooling Programm - NachhaltigWirtschaften.at IEA Solar Heating and Cooling Programm - NachhaltigWirtschaften.at
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
- Page 474 and 475: IEA SHC Task 38 Solar Air Condition
- Page 476 and 477: IEA SHC Task 38 Solar Air Condition
- Page 478 and 479: IEA SHC Task 38 Solar Air Condition
- Page 480 and 481: IEA SHC Task 38 Solar Air Condition
- Page 482 and 483: IEA SHC Task 38 Solar Air Condition
- Page 484 and 485: IEA SHC Task 38 Solar Air Condition
- Page 486 and 487: IEA SHC Task 38 Solar Air Condition
- Page 488 and 489: IEA SHC Task 38 Solar Air Condition
- Page 490 and 491: IEA SHC Task 38 Solar Air Condition
- Page 492 and 493: IEA SHC Task 38 Solar Air Condition
- Page 494 and 495: IEA SHC Task 38 Solar Air Condition
- Page 496 and 497: IEA SHC Task 38 Solar Air Condition
- Page 498 and 499: IEA SHC Task 38 Solar Air Condition
- Page 500 and 501: IEA SHC Task 38 Solar Air Condition
- Page 502 and 503: IEA SHC Task 38 Solar Air Condition
- Page 504 and 505: IEA SHC Task 38 Solar Air Condition
- Page 506 and 507: IEA SHC Task 38 Solar Air Condition
- Page 508 and 509: IEA SHC Task 38 Solar Air Condition
- Page 510 and 511: IEA SHC Task 38 Solar Air Condition
- Page 512 and 513: IEA SHC Task 38 Solar Air Condition
- Page 514 and 515: IEA SHC Task 38 Solar Air Condition
- Page 516 and 517: IEA SHC Task 38 Solar Air Condition
- Page 518 and 519: IEA SHC Task 38 Solar Air Condition
- Page 520 and 521: IEA SHC Task 38 Solar Air Condition
- Page 522 and 523: IEA SHC Task 38 Solar Air Condition
- Page 526 and 527: IEA SHC Task 38 Solar Air Condition
- Page 528 and 529: IEA SHC Task 38 Solar Air Condition
- Page 530 and 531: IEA SHC Task 38 Solar Air Condition
- Page 532 and 533: IEA SHC Task 38 Solar Air Condition
- Page 534 and 535: IEA SHC Task 38 Solar Air Condition
- Page 536 and 537: IEA SHC Task 38 Solar Air Condition
- Page 538 and 539: IEA SHC Task 38 Solar Air Condition
- Page 540 and 541: IEA SHC Task 38 Solar Air Condition
- Page 542 and 543: IEA SHC Task 38 Solar Air Condition
- Page 544 and 545: IEA SHC Task 38 Solar Air Condition
- Page 546 and 547: IEA SHC Task 38 Solar Air Condition
- Page 548 and 549: IEA SHC Task 38 Solar Air Condition
- Page 550 and 551: IEA SHC Task 38 Solar Air Condition
- Page 552 and 553: IEA SHC Task 38 Solar Air Condition
- Page 554 and 555: IEA SHC Task 38 Solar Air Condition
- Page 556 and 557: IEA SHC Task 38 Solar Air Condition
- Page 558 and 559: IEA SHC Task 38 Solar Air Condition
- Page 560 and 561: IEA SHC Task 38 Solar Air Condition
- Page 562 and 563: IEA SHC Task 38 Solar Air Condition
- Page 564 and 565: IEA SHC Task 38 Solar Air Condition
- Page 566 and 567: IEA SHC Task 38 Solar Air Condition
- Page 568 and 569: IEA SHC Task 38 Solar Air Condition
- Page 570 and 571: IEA SHC Task 38 Solar Air Condition
- Page 572 and 573: Task 38 Solar Air-Conditioning and
<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