Pac V2-1
Pac V2-1
Pac V2-1
o :
Nom
se
ei
se
ei
t
t
t
t
t
=
( ) ( ) ( )
seNom eiNom
t t + + =
se
F2
ei
F1 fNom ffl t D t D 1 P P
P
a
=P
afl
. P
f
/ P
ffl
+P
ventNom
La qualit de la modlisation dpend du nombre et du choix des donnes constructeur
disponibles. En plus, la forme adimensionnelle du modle rend possible son utilisation avec les
mmes valeurs des coefficients de pondration pour d'autres types d'installations appartenant
la mme gamme.
Les coefficients de rgressions C
i
et D
i
sont dtermins par le module PREPGF partir de 1 ou
au moins 3 points de fonctionnements
5.2 EN MODE POMPE A CHALEUR
Ce modle caractrise le groupe frigorifique dans son ensemble par ses tempratures de
sources chaudes et froides et permet dobtenir la puissance appele au compresseur ainsi que
la puissance fournie. Le modle polynomial est bas sur une quation de modlisation
quadratique adimensionnelle tablie partir de l'expression du coefficient de performance
thorique et d'une rgression partir des donnes constructeur.
La puissance appele au compresseur, P
a
, est calcule par un polynme de deuxime degr
en fonction de la temprature de l'eau la sortie de la batterie intrieure, de la temprature
d'entre de l'air sur la batterie extrieure et du rapport des puissances appele et fournie au
nominal. La puissance fournie en conditions non nominales a t modlise d'une manire
similaire. Les quations de modlisation sont les suivantes :
( )
2
C2 C1 CO
Nom
c
a
cfl
afl
t C t C C
P
P
P
P
+ +
o:
Nom
se
ei
se
ei
t
t
t
t
t
=
( )
2
C2 C1 C0
Nom fl t D t D D Pc Pc + + =
( ) ( ) ( )
seNom eiNom
t t + + =
se
C2
ei
C1 cNom fl c t D t D 1 P P
P
a
=P
afl
. P
c
/ P
cfl
+P
ventNom
rf.00.135R PAC V.2.1 Mai 2004 page 5
rf.00.135R PAC V.2.1 Mai 2004 page 6
La qualit de la modlisation dpend du nombre et du choix des donnes constructeur
disponibles. En plus, la forme adimensionnelle du modle rend possible son utilisation avec les
mmes valeurs des coefficients de pondration pour d'autres types d'installations appartenant
la mme gamme.
6 REFERENCES
Catalin BRANESCU, Comparaison des modles de groupes frigorifiques en vue dun calcul des
consommations des installations de climatisation, rapport de stage de DEA, bibliothque de
lEcole des Mines de Paris, 1996
Alexandra NITULESCU, Modlisation dun groupe frigorifique volumtrique en vue du calcul
des consommations dnergie, rapport de stage TEMPUS, Ecole des Mines de Paris
H.C. PEITSMANN, H.J . NICOLAAS, IEA Annex 10 "Liquid Chiling System", 1988, TNO Institute
of Applied Physics. The Netherlands
7 ALGORITHMES
En mode froid :
conversion des tempratures en Kelvin
calcul de la diffrence de temprature adimensionnelle
calcul de la puissance frigorifique demande
dtermination de la puissance fournie pleine charge en conditions non nominales
calcul de la puissance absorbe pleine charge en conditions non nominales
calcul de la puissance absorbe
En mode pompe chaleur :
conversion des tempratures en Kelvin
calcul de la diffrence de temprature adimensionnelle
calcul de la puissance frigorifique demande
dtermination de la puissance fournie pleine charge en conditions non nominales
calcul de la puissance absorbe pleine charge en conditions non nominales
calcul de la puissance absorbe
rf.00.135R PAC V.2.1 Mai 2004 page 7
8 PROGRAMME (CODES SOURCE FORTRAN)
!**************************************************************************
!14 - SUBROUTINE PAC () [POMPE A CHALEUR REVERSIBLE]
!***********************************************************************
! REFERENCE DE LA METHODE : CSTB/ENEA/CVA - 00.000R
! REFERENCE DE L'ALGORITHME : FICHE PAC - version 1.0
! PAC : POMPE A CHALEUR REVERSIBLE
! DEVELOPER : J R. MILLET E. FLEURY
! INFORMATISATION : A. BOLHER
! VERSION : 1.0
! DATE : 11 sept 2000
!**************************************************************************
! REFERENCE : BRANESCU C. Comparaison des modles de groupes
! frigorifiques en vue du calcul des consommations des installations de
! climatisation
! ROOF-TOP : Calculate roof-top performance at part-load and at
! no-rating temperature. Determination of sensible and
! latente energy rate.
! DEVELOPER : Olivier MORISOT PhD Studient Ecole des Mines de Paris,
! VERSION : 2.1
! DATE : 2 septembre 1999
!
! Adapted by Millet, Fleury, Bolher, CSTB
! Calculate heat pump performances at part-load and at
! no-rating temperature.
!**************************************************************************
!***********************************************************************
! INPUT VARIABLES :
!***********************************************************************
! ModePAC heating or cooling
! 0 arret
! 1 Chauffage
! 2 Refroidissemnt
! Texte outside air temperature (C)
! Tse chilled ot heated water temperature (C)
! Tee inlet water temperature (C)
!CCi,DCi i=0,1,2 regression coefficients for heating process
!CFi,DFi i=0,1,2 regression coefficients for cooling process
!
!QaRatF rating consumed energy rates for cooling process (W)
!QaRatC rating consumed energy rates for heating process (W)
!QfRat rating required energy rates for cooling process (W)
!QcRat rating required energy rates for heating process (W)
!QfanextRat rating outside fan energy rate (W)
! if QfanextRat =0, fan energy rate included in consumed energy rate
!
!TexteRatF rating outside air temperature, inlet condenser in cooling process (C)
!TexteRatC rating outside air temperature, inlet evaporator in heating process (C)
!TseRatF rating outlet evaporator water temperature in cooling process (C)
!TseRatC rating outlet condensor water temperature, in heating process (C)
!MaRat rating outside air flow rate (kg/s)
!MegRat rating inside water flow rate (kg/s)
!***********************************************************************
! OUTPUT VARIABLES :
!***********************************************************************
!Qa consumed energy rates (W)
!Qc required energy rates given by the heat pump (cooling) (W)
!Qf required energy rates given by the heat pump (heating) (W)
!***********************************************************************
! INTERNAL VARIABLES:
!***********************************************************************
! Qaffl consumed energy rates at full load for no-rating conditions (W)
! Qffl required energy rates at full load for no-rating conditions (W)
! Qacfl consumed energy rates at full load for no-rating conditions (W)
! Qcfl required energy rates at full load for no-rating conditions (W)
! Dt non-dimensionnal temperature difference (K/K)
!***********************************************************************
!
!**************************************************************************
! NOTE Inputs must be column vectors of same size.
! expect rating values which are scalair values.
rf.00.135R PAC V.2.1 Mai 2004 page 8
!**************************************************************************
SUBROUTINE PAC &
& (ModePAC,&
& Texte,&
& Tse,&
& Tee,&
& CC0,&
& CC1,&
& CC2,&
& DC0,&
& DC1,&
& DC2, &
& CF0,&
& CF1,&
& CF2,&
& DF0,&
& DF1,&
& DF2,&
& QaRatF,&
& QaRatC,&
& QfRat,&
& QcRat,&
& QfanextRat, &
& TexteRatF,&
& TexteRatC,&
& TseRatF,&
& TseRatC,&
& MegRat,&
& Qf,&
& Qc,&
& Cpeg
& Qa)
USE CONSTANTES
IMPLICIT NONE
Real Small,large,CDefaut
Integer ModePAC
Real Texte,Tse,Tee
Real QfRat,QaRatF
Real CF0,CF1,CF2,DF0,DF1,DF2,CC0,CC1,CC2,DC0,DC1,DC2
Real QfanextRat,Tcomp,CdegradF,CdegradC,TexteRatC,TseRatC
Real QcRat,QaRatC
!internes
Real texteRatF,tseRatF,Dt,Qffl,Qaffl,MegRat
Real Qcfl,Qacfl
!sorties
Real Qa,Qf,Qc
!---- Initialize constants--------------------------------------------
small =1.E-15
large =1.E15
!--------------------------------------------------------------------------
Qc=0
Qf=0
Qa=0
if (ModePAC.EQ.2) then !cooling process
!*** Convert temperatures in Kelvin
texte =Texte +273.15
tse =Tse +273.15
texteRatF =TexteRatF +273.15
tseRatF =TseRatF +273.15
tee =Tee +273.15
!*** Calculate non-dimensionnal temperature difference
Dt =texte/tse - texteRatF/tseRatF
!*** Calculate required energy rates at full load for no-rating conditions
Qffl =QfRat * (DF0 +DF1*Dt +DF2*Dt**2)
Qffl =QfRat * (DF0 +DF1*(texte-texteRatF) +DF2*(tse-tseRatF))
!*** Calculate consumed energy rates at full load for no-rating conditions
Qaffl =Qffl* QaRatF/QfRat* (CF0 +CF1*Dt +CF2*Dt**2)
!*** Calculate required energy rates
Qf =MegRat*Cpeg*(tee - tse)
if ((Qf /Qffl).GT.1) then
rf.00.135R PAC V.2.1 Mai 2004 page 9
Qf =Qffl
Qa =(Qaffl*Qf/Qffl) +QfanextRat*Qf/Qffl
else
Qa =(Qaffl*Qf/Qffl) +QfanextRat*Qf/Qffl
endif
else
! heating process
!*** Convert temperatures in Kelvin
texte =Texte +273.15
tse =Tse +273.15
texteRatC =TexteRatC +273.15
tseRatC =TseRatC +273.15
tee =Tee +273.15
!*** Calculate non-dimensionnal temperature difference
Dt =texte/tse - texteRatC/tseRatC
!*** Calculate required energy rates at full load for no-rating conditions
Qcfl =QcRat * (DC0 +DC1*Dt +DC2*Dt**2)
Qcfl=QcRat*(DC0+DC1*(texte-texteRatC)+DC2*(tse-tseRatC))
!*** Calculate consumed energy rates at full load for no-rating conditions
Qacfl =Qcfl* QaRatC/QcRat* (CC0 +CC1*Dt +CC2*Dt**2)
!*** Calculate required energy rates
Qc =MegRat*Cpeg*(tse - tee)
if ((Qc /Qcfl).GT.1) then
Qc =Qcfl
Qa =(Qacfl*Qc/Qcfl) +QfanextRat*Qc/Qcfl
else
Qa =(Qacfl*Qc/Qcfl) +QfanextRat*Qc/Qcfl
endif
endif
return
end subroutine PAC
9 EXEMPLES DE FONCTIONNEMENT
Determination des caractristiques en fonctionnement froid et chaud :
CIAT COOLER ILK ILKH 600A(3 compresseurs, 2 circuits frigorifiques)
Mode froid : on connat un point de fonctionnement nominal et deux autres points.
Temperature d'entre d'air au condenseur Texte =32C
Temprature sortie d'eau l'vaporateur Tsen =7C
Puissance absorbe compresseurs Pan =56200 W
Puissance frigorifique Pfn =146000 W
Texte (1)=28
Tse(1)=12
Pa(1)=58400
Pf(1)=180000
Texte (2)=40
Tse(2)=5
Pa(2)=58500
Pf(2)=122000
Sorties PREPGF :
C0 1.000000
C1 5.679862
C2 28.26681
C3 0.0000000
D0 1.000000
D1 -5.824059
rf.00.135R PAC V.2.1 Mai 2004 page 10
D2 36.39729
D3 0.0000000
Mode chaud : on connat un point de fonctionnement nominal et deux autres points de
fonctionnement. La temprature de dpart deau est fonction de la temprature extrieure.
Texte =0
Tsen =40
Pan=43700
Pfn=143500
Texte (1)=10
Tse(1)=30
Pa(1)=42700
Pf(1)=209000
Texte (2)=-10
Tse(2)=50
Pa(2)=39000
Pf(2)=96900
Sorties PREPGF :
C0 1.000000
C1 -5.444
C2 1.859
C3 0.0000000E+00
D0 1.000000
D1 6.469
D2 14.918
D3 0.0000000E+00
Exemple de fonctionnement
Dfinition des paramtres PAC :
Cpeg=3800
Fonctionnement en mode chaud :
ModePAC =1
Texte=7 ! Temprature d' air extrieur (C)
Tse=34 ! Temprature dpart d'eau (C)
Tee=29 ! temperature de retour d'eau(C)
!CCi,DCi i=0,1,2 coefficients de regression :
C0 1.000000
C1 -5.444
C2 1.859
C3 0.0000000E+00
D0 1.000000
D1 6.469
D2 14.918
D3 0.0000000E+00
TexteRatC=0 C
TseRatC=40 C
QaRatC=43700 W
QcRat=143500 W
rf.00.135R PAC V.2.1 Mai 2004 page 11
Fonctionnement en mode froid :
ModePAC =2
Texte=35 ! Temprature d' air extrieur (C)
Tse=7 ! Temprature dpart d'eau (C)
Tee= 11 ! temperature de retour d'eau(C)
!CFi,DFi i=0,1,2 regression coefficients for cooling process :
C0 1.000000
C1 5.679862
C2 28.26681
C3 0.0000000
D0 1.000000
D1 -5.824059
D2 36.39729
D3 0.0000000
TexteRatF=32 C
TseRatF=7 C
QaRatF=56200 W
QfRat=146000 W
QfanextRat= 3346W
MegRat=7.68 kg/s
Sorties en mode Chaud:
Qa =38279 W
Qf=0
Qc=145920 W
Sorties en mode froid:
Qa=51159 W
Qf=- 116736 W
Qc=0
10 SUIVI DES MODIFICATIONS
Auteur Date Version Modification
J .M.
ALESSANDRINI
28/05/04 V2.0->2.1 Description gnrale : pas de prise en compte du
dgivrage
Correction de Qffl et de Qcfl mise en cohrence avec
les prprocesseurs
Les exemples de fonctionnement ne sont pas corrigs