Re: Problem with a Fuel Cell

[Fernando Rivas <ferriv@xxxxxxxxxxxxxxx>]

Dear TRNSYS/HYDROGEMS users,

Although the solution given by Michaël Kummert and the HYDROGEMS 
developers  allows me to introduce the H_PEM and W_PEM parameters 
without errors I still had the same problem:  H_PEM and W_PEM were not 
used to calcuate the main electrical equations; these are calculated 
with A_PEM (A_PEM was 232 cm2 always).
I have found a provisional solution adding the next lines to the source 
code (in the MAIN_EQN subroutine):
   ...
     INTEGER RTCTMODE                        !NEWCODE
     REAL H_PEM,W_PEM                        !NEWCODE
   ...
    RTCTMODE=INT(PARAM(10))           !NEWCODE
    IF (RTCTMODE.EQ.2) THEN               !NEWCODE
        H_PEM=PARAM(12)                        !NEWCODE
        W_PEM=PARAM(13)                       !NEWCODE
        A_PEM=W_PEM*H_PEM                !NEWCODE
    ENDIF                                                    !NEWCODE
...

I have only modified the MAIN_EQN subroutine (I have added a remark to 
find it easily: !NEWCODE ). I don´t know if  the rest of  the code needs 
more modifications. I have recompiled this code and it runs correctly. 
These changes should affect to the IISiBat proformas TYPE170b, TYPE170e, 
TYPE170h and TYPE170k (RTCTMODE=2). I have attached a copy of  the 
provisional source code.
 
Regards
Fernando Rivas

Michaël Kummert escribió

>Dear TRNSYS/HYDROGEMS users,
>
>Fernando received a "personal" answer from the HYDROGEMS developers but may
>be some of you are facing the same problem, so here is a copy:
>
>This is indeed a bug in Type170. The code should not be checking the value
>of A_PEM in that configuration. 
>You can simply comment out the following lines in the code:
>
>C	  AREA=H_PEM*W_PEM
>C	  IF((A_PEM.LE.(0.99*AREA)).OR.(A_PEM.GE.(1.01*AREA))) THEN
>C	    WRITE(LUW,150) ITYPE,TIME,A_PEM,H_PEM,W_PEM
>C	    STOP
>C        ENDIF
>
>An update of the HYDROGEMS library will be available in March from the
>TRNSYS website. It will include a few bug fixes, including this one.
>
>Regards,
>
>Michaël Kummert
>
>_________________________________________________________
>
>Michaël Kummert
>
>Solar Energy Laboratory - University of Wisconsin-Madison
>1303 Engr Res Bldg, 1500 Engineering Drive
>Madison, WI 53706
>
>Tel: +1 (608) 263-1589
>Fax: +1 (608) 262-8464
>E-mail: kummert@engr.wisc.edu
>
>SEL Web Site: http://sel.me.wisc.edu 
>TRNSYS Web Site: http://sel.me.wisc.edu/trnsys
>
>
>
>
>
>-----Original Message-----
>From: owner-trnsys@relay.doit.wisc.edu
>[mailto:owner-trnsys@relay.doit.wisc.edu]On Behalf Of Fernando Rivas
>Sent: Friday, 21 February 2003 6:38 PM
>To: trnsys@relay.doit.wisc.edu
>Subject: Problem with a Fuel Cell
>
>
>Dear TRNSYS users,
>
>I am trying to simulate a PEM Fuel Cell with IISiBat 3 and I have a
>problem with TYPE170b of Hydrogems
>(OXMODE=1, TMODE=1, RTCTMODE=2).  I have the following problem:  the
>A_PEM parameter is "dummy"
>and it can not be modified and this value is provided by H_PEM and W_PEM
>parameters (H_PEM * W_PEM = A_PEM).
>However if I modify those parameters (H_PEM and W_PEM) and run the
>simulation I get this error message in the list file:
>               ********ERROR**********
>               TYPE170  TIME =    1.000
>       INVALID PEMFC GEOMETRY:
>       LENGTH AND/OR HEIGHT NOT IN AGREEMENT WITH AREA!
>       A_PEM =  232.00[cm2]
>       H_PEM =   16.00 cm
>       W_PEM =   16.00 cm
>               SIMULATION 'STOPPED'
>               ***********************
>
>If A_PEM can not be modified (it has a value by default of 232.00 [cm2]
>always), How can I change the A_PEM value
>(in example 16.00*16.00 = 256.00 [cm2])?
>
>Thank you in advance
>Fernando Rivas
>