[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: [TRNSYS-users] PSYCHROMETRICS routine issue



Dear all

 

Eventually the same kind of problem appears with standard Type 33.a.

 

*TRACE*  UNIT   3  TYPE  33  AT TIME   1.3750000000000000E+00  ITERATION    1  CALL      25  POST-CONVERGENCE CALL

P  6.0000000000000000E+00  1.0000000000000000E+00  2.0000000000000000E+00

I  8.5248000000000025E-03  5.1599257000000001E+01  1.0000000000000000E+00

O  8.5248000000000025E-03  8.5248000000000025E-03  5.1599257000000001E+01  1.2856158812905980E+00  1.2747489018520892E+00

O  1.0000000000000000E+02  8.5248000000000025E-03  8.5248000000000025E-03  0.0000000000000000E+00  1.0000000000000000E+00

 

T = O[7] = 8.5E-3 [°C]

x = O[1] = 8.5E-3 []

h = O[3] = 51.6 [kJ/kg]

(No warning printed in lst file.)

 

I am not sure if it is related with the issue of PSYCHROMETICS routine below.

Thank you for your opinion.

 

Regards

--

Antoine Gautier
Solam
en
+33 6 98 78 56 46

antoine.gautier@solamen.fr

 

De : Antoine Gautier [mailto:antoine.gautier@solamen.fr]
Envoyé : mardi 31 mai 2011 09:52
À : 'trnsys-users@cae.wisc.edu'
Objet : RE: PSYCHROMETRICS routine issue

 

Dear all

 

Same problem illustrated in a more simplified way.

 

source code

 

MODE=6

      PSYDAT(6)=0.8524754597715760D-02

      PSYDAT(7)=0.5159925701758463D+02

      CALL PSYCHROMETRICS(TIME,INFO,IUNITS,MODE,WBMODE,

     &                                     PSYDAT,EMODE,STATUS)

      write(myMessage,'(" PSYDAT(2)=",D," PSYDAT(6)=",D,

     &            " PSYDAT(7)=",D)') PSYDAT(2),PSYDAT(6),PSYDAT(7)

      CALL MESSAGES(-1,trim(myMessage),

     &                  'Notice',info(1),info(2))

 

simulation response

 

*** Notice at time    :   9.125000

      Generated by Unit     :     2

      Generated by Type     :   206

Message               :  PSYDAT(2)=   0.1732587158021238D+02 PSYDAT(6)=   0.8524754597715760D-02 PSYDAT(7)=   0.5159925701758463D+02

 

DB Temperature = 17.3°C instead of ~ 30°C for that state point.

Is there any explanation?

 

Thanks a lot.

Regards

--

Antoine Gautier
Solam
en
+33 6 98 78 56 46

antoine.gautier@solamen.fr

 

De : Antoine Gautier [mailto:antoine.gautier@solamen.fr]
Envoyé : lundi 30 mai 2011 18:50
À : 'trnsys-users@cae.wisc.edu'
Objet : PSYCHROMETRICS routine issue

 

Dear all

 

Please have a look below at the short excerpt of a source code.

 

      MODE=6

      PSYDAT(6)=wAirEnt

      PSYDAT(7)=hAirEnt

      CALL PSYCHROMETRICS(TIME,INFO,IUNITS,MODE,WBMODE,

     &                                     PSYDAT,EMODE,STATUS)

      write(myMessage,'(" PSYDAT(2)=",D," PSYDAT(6)=",D,

     &            " PSYDAT(7)=",D)') PSYDAT(2),PSYDAT(6),PSYDAT(7)

      CALL MESSAGES(-1,trim(myMessage),

     &                  'Notice',info(1),info(2))

 

      PSYDAT(2)=29.82D0

 

      CALL PSYCHROMETRICS(TIME,INFO,IUNITS,MODE,WBMODE,

     &                          PSYDAT,EMODE,STATUS)

      write(myMessage,'(" PSYDAT(2)=",D," PSYDAT(6)=",D,

     &            " PSYDAT(7)=",D)') PSYDAT(2),PSYDAT(6),PSYDAT(7)

      CALL MESSAGES(-1,trim(myMessage),

     &                  'Notice',info(1),info(2))               

 

I know PSYDAT(2) is not used in MODE 6 but please have a look now at the response I get…

 

*** Notice at time        :        16.375000

    Generated by Unit     :     2

    Generated by Type     :   206

    Message               :  PSYDAT(2)=   0.1732587158021238D+02 PSYDAT(6)=   0.8524754597715760D-02 PSYDAT(7)=   0.5159925701758463D+02

 

 

*** Notice at time        :        16.375000

    Generated by Unit     :     2

    Generated by Type     :   206

    Message               :  PSYDAT(2)=   0.2981999969482422D+02 PSYDAT(6)=   0.8524754597715760D-02 PSYDAT(7)=   0.5159925701758463D+02

 

Same humidity ratio, same enthalpy and quite different dry bulb temperature. Why the second PSYCHROMETRICS statement does not rewrite on PSYDAT(2)?

Any help would be quite appreciated as I am slowly but surely getting crazy!

 

Thanks a lot.

Regards

--

Antoine Gautier
Solam
en
+33 6 98 78 56 46

antoine.gautier@solamen.fr