C C TYPE 205 : PIPE C C CALL THE SUBROUTINE PIPE WHICH CALCULATE THE HEAT TRANSFER C COEFFICENT USING SIPLIFIED LAWS, THE HEAT TRANSFER OF THE C PIPE AND THE EXHAUST TEMPERATURE OF THE FLUID GOING THROUGH C THE PIPE. C C C C SI UNITS ! C C PARAMETERS : PAR(1)=XL : LENGTH OF THE PIPE C PAR(2)=D : DIAMETER OF THE PIPE C PAR(3)=HV : HV >= 0 THEN VERTICAL PIPE ELSE HORIZONTAL C PAR(4)=EPS: EMISION COEFFICIENT FOR OUTER SURFACE C PAR(5)=CPW: MASSIC HEAT OF FLUID C C INPUTS : XIN(1)=TSU : SUPPLY FLUID TEMPERATURE C XIN(2)=XMW : FLUID MASSFLOW RATE C XIN(3)=TAMB: TEMPERATURE OF SURROUNDING WALLS AND AIR C C OUTPUTS : OUT(1)=TEX : EXHAUST TEMPERATURE OF THE FLUID C OUT(2)=XMW : FLUID MASSFLOW RATE C OUT(3)=Q : HEAT TRANSFER C CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC C C SAMPLE INPUT AND OUTPUT VALUES C C PARAMETERS SAMPLE VALUES C PAR(1) XL 3.5 C PAR(2) D 0.07 C PAR(3) HV -1.0 C PAR(4) EPS 0.5 C PAR(5) CPW 4187.0 C C INPUTS C XIN(1) TSU 68.65 C XIN(2) XMW 2.4 C XIN(3) TAMB 15.0 C C OUTPUTS C OUT(1) TEX 68.61 C OUT(2) XMW 2.4 C OUT(3) Q 430.3 C CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC C ! This component has been assigned Type Number 205. If that number conflicts with ! another user Type number, you will need to change it and recompile the appropriate ! dll. SUBROUTINE TYPE205(TIME,XIN,OUT,T,DTDT,PAR,INFO,ICNTRL,*) !DEC$ATTRIBUTES DLLEXPORT :: TYPE205 DOUBLE PRECISION XIN,OUT DIMENSION XIN(3),OUT(3),PAR(5),INFO(15) ! Set the version information for TRNSYS IF (INFO(7).EQ.-2) THEN INFO(12) = 15 RETURN 1 ENDIF C C FIRST CALL C IF(INFO(7).EQ.-1) CALL TYPECK(1,INFO,3,5,0) C C INPUTS C TSU=XIN(1) XMW=XIN(2) TAMB=XIN(3) C C PARAMETERS C XL=PAR(1) D=PAR(2) HV=PAR(3) EPS=PAR(4) CPW=PAR(5) C C CALL SUBROUTINE PIPE C CALL PIPE(TSU,XMW,TAMB,XL,D,HV,EPS,CPW,TEX,Q) C C OUTPUTS C OUT(1)=TEX OUT(2)=XMW OUT(3)=Q C RETURN 1 END C C C SUBROUTINE PIPE C C CALCULATE THE HEAT TRANSFER COFFEFICIENT (SIMPLIFIED C FORMULA), THE HEAT TRANSFER FRO THE FLUID AND THE EXHAUST C TEMPERATURE OF THE FLUID. C C SEE DOCUMENT AN10-880603-02 OF THE ANNEX10 DATA BANK. C SUBROUTINE PIPE(TSU,XMW,TAMB,XL,D,HV,EPS,CPW,TEX,Q) C C C "PARAMETERS": XL : LENGTH OF THE PIPE C D : DIAMETER OF THE PIPE C HV : HV >= 0 THEN VERTICAL PIPE ELSE HORIZONTAL C EPS: EMISION COEFFICIENT FOR OUTER SURFACE C CPW: MASSIC HEAT OF FLUID C C "INPUTS" : TSU : SUPPLY FLUID TEMPERATURE C XMW : FLUID MASSFLOW RATE C TAMB: TEMPERATURE OF SURROUNDING WALLS AND AIR C C "OUTPUTS" : TEX : EXHAUST TEMPERATURE OF THE FLUID C Q : HEAT TRANSFER C C OTHER VARIABLES C C ARE THE VALUES ACCEPTABLE ? C C IS THE MASSFLOW RATE POSITVE ? COMMON /LUNITS/LUR,LUW,LUK,IFORM IF(XMW.LT.0.0) THEN WRITE(LUW,*) 'NEGATIVE MASSFLOW RATE IN SUBROUTINE PIPE' CALL MYSTOP(1001) RETURN ENDIF C ARE THE DIMENSIONS POSITIVE ? IF((D.LT.0.0).OR.(XL.LT.0.0)) THEN WRITE(LUW,*) 'NEGATIVE DIMENSION FOR PIPE IN SUBROUTINE PIPE' CALL MYSTOP(1001) RETURN ENDIF C ARE THE TEMPERATURES HIGHER THAN 0 ? IF((TAMB.LT.0.0).OR.(TSU.LT.0.0)) THEN WRITE(LUW,*) 'TEMPERATURE LOWER THAN 0' C STOP ENDIF C VERY LOW MASSFLOW RATE IF(XMW.LT.1.E-5) THEN TEX=TAMB GO TO 111 ENDIF C C AREA CALCULATION C A=3.1416*D*XL C C HEAT TRANSFER COEFFICIENT DUE TO RADIATION C HR=5.7E-8*EPS*((TSU+273)**2+(TAMB+273)**2)*(TSU+273+TAMB+273) C C HEAT TRANSFER COEFFICIENT DUE TO CONVECTION C DELTAT=TSU-TAMB IF(DELTAT.LT.0.0) DELTAT=-DELTAT IF(HV.GE.0.0) THEN C VERTICAL PIPE HC=1.3*DELTAT**(1./3.) ELSE C HORIZONTAL PIPE HC=1.3*(DELTAT/D)**0.25 ENDIF C C TOTAL HEAT TRANSFER COEFFICIENT C U=HR+HC C C NO AREA OR HEAT TRANSFER COEFFICIENT = 0 C IF((U.EQ.0.0).OR.(A.EQ.0.0)) THEN TEX=TSU GO TO 111 ENDIF C C EXHAUST TEMPERATURE C PROV1=(U*A/XMW/CPW) IF(PROV1.GT.20) THEN TEX=TAMB ELSE PROV=EXP(PROV1) TEX=(TSU-TAMB*(1-PROV))/PROV ENDIF C C HEAT TRANSFER C 111 Q=XMW*CPW*(TSU-TEX) C RETURN END