C C C LIST OF SUBROUTINES IN THIS FILE C C HC_THLEQN : LIQUID THERMAL CONDUDCITIVITY EQUATIONS C HC_THLSATINI : LATINI ET AL. METHOD C C******************************************************************* C C KDB THEMOPHYSICAL PROPERTIES CALCULATION FORTRAN ROUTINE LIBRARY C C [NAME ] HC_THLEQN C C [TYPE ] FORTRAN SUBROUTINE C C [PURPOSE] LIQUID THERMAL CONDUCTIVITY CALCULATION USING INTERNAL CORRELATION EQUATION C C [USAGE ] CALL HC_THLEQN(ICN,T,THL,IST) C C [ARGUMENTS] C ICN : COMPONENT NUMBER (1-50) TO CALCULATE LIQUID THERMAL CONDUCTIVITY C (INTEGER, INPUT) C T : TEMPERATURE IN KELVIN (REAL*8, INPUT) C C THL : LIQUID THERMAL CONDUCTIVITY IN W/m.K (REAL*8, OUTPUT) C IST : STATUS OF CALCULATION (INTEGER, OUTPUT) C = 0 : NORMAL TERMINATION C = 801 : LIQUID THERMAL CONDUCTIVITY COEFFICIENT NOT AVAILABLE C = 802 : OUT OF RANGE FOR THE APPLICATION C C [COMMENTS] C C [REQUIRED COMMON BLOCKS] C COMMON /HC_KTHL/ LIQUID THERMAL CONDUCTIVITY COEFFICIENTS C C [REQUIRED SUBROUTINES OR FUNCTIONS] C NONE C C [REFERENCE] C NONE C C [REVISION INFORMATION] C 1.PROGRMAMMED BY J.W.KANG, KOREA UNIVERSITY, 1996 C 2.REVISED BY J.W.KANG, KOREA UNIVERSITY, 2000 C 3.REVISED BY J.W.KANG, KOREA UNIVERSITY, 2001 C 4.REVISED BY Y.S.KIM, KOREA UNIVERSITY, 2002 C******************************************************************* SUBROUTINE HC_THLEQN(ICN,T,THL,IST) IMPLICIT DOUBLE PRECISION (A-H,O-Z) INTEGER IST,ICN REAL*8 T,THL COMMON/HC_KTHL/IEQ_THL(50),THL_A(50),THL_B(50),THL_C(50), 1 THL_D(50),THL_E(50),THL_F(50),THL_G(50),THL_T1(50),THL_T2(50) C --- INITIALIZATION IST = 0 THL = 0.0 C --- CHECK IF EQUATION COEFF.S AVAILABLE IF(IEQ_THL(ICN).LE.0) IST = 801 C --- CHECK IF THE RANGE OF APPLICATION IF(IST.NE.0) RETURN C --- LIQUID HEAT CAPAICITY CALCULATION THL =THL_A(ICN)+THL_B(ICN)*T+THL_C(ICN)*T*T+THL_D(ICN)*T*T*T RETURN END C******************************************************************* C C KDB THEMOPHYSICAL PROPERTIES CALCULATION FORTRAN ROUTINE LIBRARY C C [NAME ] HC_THLLATINI C C [TYPE ] FORTRAN SUBROUTINE C C [PURPOSE] LIQUID THERMAL CONDUCTIVITY CALCULATION USING LATINI ET AL. METHOD C C [USAGE ] CALL HC_THLLATINI(ICN,T,THL,IST) C C [ARGUMENTS] C ICN : COMPONENT NUMBER (1-50) TO CALCULATE LIQUID THERMAL CONDUCTIVITY C (INTEGER, INPUT) C T : TEMPERATURE IN KELVIN (REAL*8, INPUT) C C THL : LIQUID THERMAL CONDUCTIVITY IN W/m.K (REAL*8, OUTPUT) C IST : STATUS OF CALCULATION (INTEGER, OUTPUT) C = 0 : NORMAL TERMINATION C = 811 : CRITICAL TEMPERATURE DATA NOT AVAILABLE C = 812 : NORMAL BOILING POINT DATA NOT AVAILABLE C = 813 : MOLECULAR WEIGHT DATA NOT AVAILABLE C = 814 : GIVEN T EXCEEDS CRITICAL TEMPERATURE C = 815 : OUT OF APPLICATION C C [COMMENTS] C C [REQUIRED COMMON BLOCKS] C COMMON /HC_KTHL/ LIQUID THERMAL CONDUCTIVITY COEFFICIENTS C C [REQUIRED SUBROUTINES OR FUNCTIONS] C NONE C C [REFERENCE] C NONE C C [REVISION INFORMATION] C 1.PROGRMAMMED BY J.W.KANG, KOREA UNIVERSITY, 1996 C 2.REVISED BY J.W.KANG, KOREA UNIVERSITY, 2000 C 3.REVISED BY J.W.KANG, KOREA UNIVERSITY, 2001 C 4.REVISED BY Y.S.KIM, KOREA UNIVERSITY, 2002 C******************************************************************* SUBROUTINE HC_THLLATINI(ICN,T,THL,IST) IMPLICIT DOUBLE PRECISION (A-H,O-Z) REAL*8 THL INTEGER ICN,K CHARACTER*50 CNAME CHARACTER*20 FORMULA COMMON/HC_NAME/ICMP(50),CNAME(50),FORMULA(50),ICLASS(50) COMMON/HC_PROP/WT(50),TB(50),TF(50),TC(50),PC(50),VC(50),ZC(50) 1 ,ACCF(50),WSRK(50),VEST(50),ZRA(50),SOLP(50),VOLP(50) 2 ,QI(50),RI(50),DM(50) DIMENSION AA(11),AL(11),BE(11),GAM(11) DATA (AA(J),J=1,11)/3.5D-3,3.61D-2,3.1D-2,3.46D-2,3.39D-3, 1 3.19D-3,3.83D-3,4.15D-2,3.85D-2,5.62D-1, 1 4.94D-1/ DATA (AL(J),J=1,11)/9*1.2D0,2*0.D0/ DATA (BE(J),J=1,11)/0.5D0,3*1.D0,3*0.5D0,2*1.D0,2*0.5D0/ DATA (GAM(J),J=1,11)/9*0.167D0,-0.167D0,-0.167D0/ C-------------------------- INITIALIZATION -------------------------C IST = 0 I = ICN C---------------- SCAN IF INSUFFICIENT DATA AVAILABLE --------------C IF(TC(I).LE.0.D0) THEN IST = 811 RETURN ENDIF IF(TB(I).LE.0.D0) THEN IST = 812 RETURN ENDIF IF(WT(I).LE.0.D0) THEN IST = 813 RETURN ENDIF C----------------------- CHECKING TEMP. RENAGE ---------------------C IF(T.GE.TC(I)) THEN IST = 814 RETURN ENDIF C--- IDENTIFICATION OF PARAMETERS FOR LATINI ET AL. METHOD IF(ICLASS(I).GE.1 .AND. ICLASS(I).LE.2) K = 1 ! SATURATED HYDROCARBONS IF(ICLASS(I).GE.3 .AND. ICLASS(I).LE.5) K = 2 ! OLEFINS IF(ICLASS(I).GE.6 .AND. ICLASS(I).LE.7) K = 3 ! CYCLOPARAFFINS IF(ICLASS(I).GE.8 .AND. ICLASS(I).LE.10) K = 4 ! AROMATICS IF(ICLASS(I).GE.11 .AND. ICLASS(I).LE.14) K = 5 ! ALCOHOLS IF(ICLASS(I).GE.15 .AND. ICLASS(I).LE.19) K = 6 ! ACIDS(ORGANICS) IF(ICLASS(I).GE.28 .AND. ICLASS(I).LE.30) K = 7 ! KETONS IF(ICLASS(I).GE.25 .AND. ICLASS(I).LE.27) K = 8 ! ESTERS IF(ICLASS(I).GE.20 .AND. ICLASS(I).LE.24) K = 9 ! ETHERS IF(ICLASS(I).EQ.48) THEN ! HALOGENATED PARAFFINS K = 11 IF( ICMP(I).EQ.1514 .OR. ICMP(I).EQ.1516) K = 10 ! R-20 IF( ICMP(I).EQ.1517 ) K =10 ! R-21 IF( ICMP(I).EQ.1515 ) K =10 ! R-22 IF( ICMP(I).EQ.1521 ) K =10 ! R-23 ENDIF C----------------------- LATINI et al. METHOD ----------------------C IF(K.GE.1.AND.K.LE.11) THEN TR=T/TC(I) VALA = AA(K)*TB(I)**AL(K)/WT(I)**BE(K)/TC(I)**GAM(K) THL = (VALA*(1.D0-TR)**0.38D0)/(TR**(1.D0/6.D0)) ELSE IST = 815 RETURN ENDIF RETURN END C******************************************************************* C END OF FILE