[TRNSYS-users] PSYCHROMETRICS routine issue

[Antoine Gautier]

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
Solamen
+33 6 98 78 56 46

antoine.gautier at solamen.fr

 

De : Antoine Gautier [mailto:antoine.gautier at solamen.fr] 
Envoyé : mardi 31 mai 2011 09:52
À : 'trnsys-users at 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
Solamen
+33 6 98 78 56 46

antoine.gautier at solamen.fr

 

De : Antoine Gautier [mailto:antoine.gautier at solamen.fr] 
Envoyé : lundi 30 mai 2011 18:50
À : 'trnsys-users at 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
Solamen
+33 6 98 78 56 46

antoine.gautier at solamen.fr

 

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.onebuilding.org/pipermail/trnsys-users-onebuilding.org/attachments/20110531/4886afda/attachment-0005.htm>