The following lines appear in the elways-called 014.htm CHAR(ABT,ABP) ! absolute temperature and pressure for 089 ABT=(TEM1+TEMP0); ABP=(P1+PRESS0) ! default setting mesg(print-out from head of library case 089 abt temp0 abp press0 Fluid properties of several fluids according to: KD Hagen,"Heat Transfer with applications", Prentice-Hall, 1999 ** Ideal Gases Air Ammonia [NH3] Argon [Ar] Carbon Dioxide [CO2] Carbon Monoxide [CO] Helium [He] Hydrogen [H2] Methane [CH4] Nitrogen [N2] Oxygen [O2] Superheated Water vapour [H2O] ** saturated vapours saturated_water_vapor [H2O] saturated vapor Refrigerant-12 [CCl2F2] saturated vapor Refrigerant-134a [CH2FCF3] --------------- *** LIQUIDS *** --------------- ** saturated liquids saturated_water SAE 5W-30 engine oil saturated liquid SAE 10W-30 engine oil saturated liquid SAE 20W-20 engine oil saturated liquid Ethylene Glycol [C2H4(OH)2] saturated liquid Ethylene Glycol 50% by volume aqueous solution saturated liquid Gasoline saturated liquid Glycerin [C3H3(OH)3] saturated liquid Refrigerant-12 [CCl2F2] saturated liquid Refrigerant-134a [CH2FCF3] saturated liquid Therminol 59 saturated liquid Therminol 66 saturated liquid Dowtherm A saturated liquid Syltherm 800 saturated liquid FC-72 saturated liquid HFE-7100 saturated liquid Important note: Temperature MUST be expressed in degrees Kelvin char(rho_expression);rho_expression=XXX char(emu_expression);emu_expression=XXX char(enu_expression);enu_expression=XXX char(cp_expression);cp_expression=XXX char(cond_expression);cond_expression=XXX ABT is now declared as CHAR in 014.htm and set = TEM1+TEMP0 ABP is now declared as CHAR in 014.htm and set = P1+PRESS0 If 089 is loaded when TEM1 or P1 are already absolute; appropriate subtractions of TEMP0 or PRESSO must be made prior to loading CASE :fluid_name: OF ------------------- *** IDEAL GASES *** ------------------- WHEN Air,3 Air rho_expression = :ABP: / (286.7 * :ABT:) cp_expression=POL6((:abt:),1083.21,-.71465,2.08832e-3,$ -2.31029E-6,1.35272E-9,-4.12726E-13,5.27347E-17) enu_expression=1.e-7*POL6(:abt:,4.12235,.72111,-4.28123e-4,5.86$ 408E-8,1.25726E-10,-6.6876E-14,1.05134E-17)/rho1 cond_expression=.001*POL6(:abt:,.81363,.08083,6.13727E-5,$ -2.10686E-7,2.06091E-10,-8.42775E-14,1.27292E-17) WHEN Ammonia,7 ------------- Ammonia [NH3] ------------- rho_expression = :ABP: /(488 * :ABT:) cp_expression=POL3(:abt:,2666.4,-5.17638,0.01455,-9.98984e-6) emu_expression=1.e-7*POL3(:abt:,96.51,-.31668,1.49445e-3, -1.11027E$ -6) cond_expression=.001*POL3(:abt:,44.8,-.28639,9.27401E-4,-6.47476E-7) WHEN Argon,5 ---------- Argon [Ar] ---------- rho_expression = :ABP: /(208 * :ABT:) cp_expression=10^POL3(log10(:abt:),4.39149,-1.89717,0.71455, -.0894$ 8) emu_expression=1.e-7*POL3(:abt:,-4.28101,.92898,-5.76258e-4, 2.0328$ 6E-7) cond_expression=.001*POL3(:abt:,.30374,.06705,-3.38561E-5,1.01005E-$ 8) WHEN Carbon_Dioxide,14 -------------------- Carbon Dioxide [CO2] -------------------- rho_expression = :ABP: /(188.9 * :ABT:) cp_expression=POL3(:abt:,438.38,1.78238,-1.5182e-3,5.42004e-7) emu_expression=1.e-7*POL3(:abt:,-11.88865,.611933,-2.68668e-4, 6.44$ 312E-8) cond_expression=.001*POL3(:abt:,-.14761,.03134,1.02456E-4,-6.92754E$ -8) WHEN Carbon_Monoxide,15 -------------------- Carbon Monoxide [CO] -------------------- rho_expression = :ABP: /(297 * :ABT:) cp_expression=POL3(:abt:,1098.51,-.43638,9.56328e-4,-4.31702e-7) emu_expression=1.e-7*POL3(:abt:,12.4382,.66417,-4.36692e-4, 1.55269$ E-7) cond_expression=.001*POL3(:abt:,-1.11133,.10419,-6.56303E-5,2.98195$ E-8) WHEN Helium,6 ----------- Helium [He] ----------- rho_expression = :ABP: /(2079 * :ABT:) emu_expression=1.e-7*POL3(:abt:,53.934,.51628,-1.53395e-4, 3.09242E$ -8) cond_expression=.001*POL3(:abt:,37.147,.40379,-1.00316E-4,1.76137E-$ 8) WHEN Hydrogen,8 ------------- Hydrogen [H2] ------------- rho_expression = :ABP: /(4124 * :ABT:) cp_expression=POL6(:abt:,7359.,52.768,-.14776,2.01558e-4, -1.4161$ 6E-7,4.94995E-11,-6.78504E-15)) emu_expression=1.e-7*POL3(:abt:,21.297,.2458,-8.31987e-5, 1.74024E-$ 8) cond_expression=.001*POL3(:abt:,27.535,.54203,-1.7311E-4,5.85159E-8) WHEN Methane,7 ------------- Methane [CH4] ------------- rho_expression = :ABP: /(518.3 * :ABT:) cp_expression=POL3(:abt:,2632.,-6.29939,.02128,-1.53535e-5) emu_expression=1.e-7*POL3(:abt:,-4.26905,0.46482,-2.95931e-4, 1.030$ 3E-7) cond_expression=.001*POL3(:abt:,5.34048,.04799,1.77229E-4,-5.45455E$ -8) WHEN Nitrogen,8 ------------- Nitrogen [N2] ------------- rho_expression = :ABP: /(296.8 * :ABT:) cp_expression=POL3(:abt:,1093.84,-.37167,7.41563e-4,-2.94906e-7) emu_expression=1.e-7*POL3(:abt:,8.82353,.66875,-3.99709e-4, 1.20754$ E-7) cond_expression=.001*POL3(:abt:,-.02104,.10218,-6.02143E-5,2.28582E$ -8) WHEN Oxygen,6 ----------- Oxygen [O2] ----------- rho_expression = :ABP: /(259.8 * :ABT:) cp_expression=POL3(:abt:,1068.85,-1.59151,5.38758e-3,-7.08515e-8) emu_expression=1.e-7*POL3(:abt:,4.70598,.7951,-4.61161e-4, 1.37102E$ -7) cond_expression=.001*POL3(:abt:,.19674,.09693,-3.85719E-5,1.19573E-$ 8) WHEN superheated_water_vapour,24 ----------------------------------------------------- superheated water vapor [H2O] ----------------------------------------------------- rho_expression = :ABP: /(461.5 * :ABT:) cp_expression=POL4(:abt:,4691.7,-16.279,.03485,-3.16183e-5, 1.06647$ E-8) emu_expression=1.e-7*POL3(:abt:,-50.12,.50668,-1.59076e-4, 7.93564E$ -8) cond_expression=.001*POL3(:abt:,9.99947,-6.9857E-3,1.39211E-4,-4.43$ 512E-8) WHEN saturated_water_vapor,21 --------------------------- saturated_water_vapor [H2O] --------------------------- rho_expression=10^POL5(:abt:,-33.82056,.27917,-9.99909e-4, 1.941$ 68E-6,-1.96966E-9,8.23222E-13)) cp_expression=10^POL6(:abt:,59.30704,-.8509,5.31125e-3, -1.74424E$ -5,3.17764E-8,-3.04367E-11,1.19865E-14) enu_expression=1.e-7*POL5(:abt:,-2415.7,30.87097,-.15357, 3.8219E$ -4,-4.69789E-7,2.2805E-10)/(RHO1) cond_expression=.001*POL6(:abt:,2724.6,-41.81772,.26527,-8.86314E-4,$ 1.64831E-6,-1.61913e-9,6.57863e-13) WHEN saturated_vapor_Refrigerant-12,30 --------------------------------------- saturated vapor Refrigerant-12 [CCl2F2] --------------------------------------- rho_expression=POL4(:abt:,992.07,-15.84475,.09623,-2.66144e-4, 2.87$ 072E-7) cp_expression=POL3(:abt:,-3057.96,40.79543,-.15865,2.151e-4) emu_expression=1.e-6*POL3(:abt:,-44.055,.54351,-1.91416e-3, 2.47345$ E-6) cond_expression=.001*POL3(:abt:,-21.428,.25441,-8.54618E-4,1.1914E-$ 6) WHEN saturated_vapor_Refrigerant-134a,32 ------------------------------------------ saturated vapor Refrigerant-134a [CH2FCF3] ------------------------------------------ rho_expression=((.020213*:abt:-3.87325)/(1.-2.08494e-3*:abt:))^2 cp_expression=POL4(:abt:,44886.3,-676.907,3.85902,-9.73043e-3, 9.21$ 174E-6) cond_expression=.001*POL3(:abt:,-97.45182,1.04946,-3.6058E-3,4.4959$ 2E-6) --------------- *** LIQUIDS *** --------------- WHEN saturated_water,15 ------------------------------------- saturated_water ------------------------------------- rho_expression=POL5(:abt:,2446.,-20.6741,.11576,-3.12895e-4,4.0505$ E-7,-2.05460E-10) cp_expression=exp((8.29041-.012557*:abt:)/(1.-1.52373e-3*:abt:)) emu_expression=1.e-7*exp((1.12646-.039638*:abt:)/(1.-7.29769E-3*$ :abt:)) real(resenu) resenu=1.e-7*exp((1.12646-.039638*293.17)/(1.-7.29769E-3*$ 293.17))/1000. resenu cond_expression=.001*POL5(:abt:,62.282,-1.768417,.03499,-1.15706E-4,$ 1.53599E-7,-7.7477E-11) WHEN SAE_5W-30_engine_oil,20 ------------------------------------- SAE 5W-30 engine oil saturated liquid ------------------------------------- rho_expression=1052.3-0.6420*:abt: cp_expression=753.7+3.65*:abt: emu_expression=10.^(POL4(:abt:,58.2987,-.53817,1.92827e-3,-3.16448$ E-6,1.97922E-9)-2) cond_expression=0.1447-2.3073E-5*:abt: WHEN SAE_10W-30_engine_oil,21 -------------------------------------- SAE 10W-30 engine oil saturated liquid -------------------------------------- rho_expression=1064.2-0.6420*:abt: cp_expression=748.6+3.63*:abt: emu_expression=10.^(POL4(:abt:,63.8534,-.59352,2.13831e-3, -3.52257$ E-6,2.20961E-9)-2) cond_expression=0.1427-2.2763E-5*:abt: WHEN SAE_20W-20_engine_oil,21 -------------------------------------- SAE 20W-20 engine oil saturated liquid -------------------------------------- rho_expression=1084.4-0.665*:abt: cp_expression=743.1+3.6*:abt: emu_expression=10.^(POL4(:abt:,83.2643,-.79412,2.91605e-3, -4.87078$ E-6,3.08743E-9)-2) cond_expression=0.1406-2.2426E-5*:abt: WHEN Ethylene_Glycol,15 -------------------------------------------- Ethylene Glycol [C2H4(OH)2] saturated liquid -------------------------------------------- rho_expression=POL3(:abt:,1082.22,1.34031,-5.63885e-3, 4.97135E-6) cp_expression=POL3(:abt:,889.59,5.27817,-9.16533e-4, 8.80859E-7) emu_expression=10.^(POL3(:abt:,21.5824,-.149,3.40615e-4, -2.71336E-$ 7)-2) cond_expression=.001*POL3(:abt:,526.92,-.88347,3.1519E-4,-2.42873E-$ 7) WHEN Ethylene_Glycol_50%_by_volume_aqueous_solution,46 --------------------------------------------------------------- Ethylene Glycol 50% by volume aqueous solution saturated liquid --------------------------------------------------------------- rho_expression=POL3(:abt:,984.39,1.06485,-2.67192e-3, 2.53477E-7) cp_expression=POL3(:abt:,2177.05,3.605,7.77104e-4, -7.75096E-7) emu_expression=10.^(POL4(:abt:,65.677,-.73037,3.09175e-3, -5.93326E$ -6,4.30591E-9)-2) cond_expression=.001*POL3(:abt:,-239.67,3.70172,-6.14252E-3,2.49288$ E-6) WHEN Gasoline,8 ------------------------- Gasoline saturated liquid ------------------------- rho_expression=POL3(:abt:,898.4,-.24536,-7.65734e-4, -3.885E-7) cp_expression=POL3(:abt:,2302.,-7.784567,.03022, -2.2096E-5) emu_expression=1.e-6*POL3(:abt:,7733.,-50.393,.11339, -8.65385E-5) cond_expression=.001*POL3(:abt:,288.74,-1.15324,2.62273E-3,-2.36257$ E-6) WHEN Glycerin,8 ------------------------------------- Glycerin [C3H3(OH)3] saturated liquid ------------------------------------- rho_expression=POL2(:abt:,1124.6,1.33754,-2.99136e-3) cp_expression=POL2(:abt:,896.95,4.52584,1.43224e-3) emu_expression=10.^(POL3(:abt:,32.1741,-.20448,4.60137e-4, -3.81826$ E-7)-2) cond_expression=.001*POL2(:abt:,229.67,.1882,-7.71975E-5) WHEN Refrigerant-12,14 ---------------------------------------- Refrigerant-12 [CCl2F2] saturated liquid ---------------------------------------- rho_expression=POL2(:abt:,1728.1,.8676,-7.62887e-3) cp_expression=POL2(:abt:,1727.2,-7.9193,.01825) emu_expression=1.e-6*POL2(:abt:,2389.4,-12.67181,.01782) cond_expression=.001*POL2(:abt:,190.29,-.46193,1.73576E-4) WHEN Refrigerant-134a,16 ------------------------------------------- Refrigerant-134a [CH2FCF3] saturated liquid ------------------------------------------- rho_expression=POL3(:abt:,3287.84,-17.85245,.062308,-8.66841e-5) cp_expression=POL3(:abt:,-6077.85,84.02833,-.3266,4.32984e-4) emu_expression=1.e-6*POL3(:abt:,9274.56,-81.60132,.24932, -2.60198E$ -4) cond_expression=.001*POL3(:abt:,223.81,-.50536,1.41171E-4,-1.45688E$ -7) WHEN Therminol_59,12 ------------------------------------------- Therminol 59 saturated liquid ------------------------------------------- Therminol 59 is registered trademark of Monsanto Company. rho_expression=POL2(:abt:,1165.09,-.57939,-2.34264e-4) cp_expression=POL2(:abt:,784.12,2.90052,5.30178e-4) enu_expression=10^(POL6(:abt:,249.725,-3.60213,.0216, -6.87146E-$ 5,1.21942E-7,-1.14294E-10,4.41809E-14)-2)/(RHO1) cond_expression=.001*POL2(:abt:,130.35,9.10712E-3,-1.3588E-4) WHEN Therminol_66,12 ------------------------------------------- Therminol 66 saturated liquid ------------------------------------------- Therminol 66 is registered trademark of Monsanto Company. rho_expression=POL2(:abt:,1164.28,-.439676,-3.193e-4) cp_expression=POL2(:abt:,646.67,2.86896,8.57054e-4) ? emu_expression=10^(POL4(:abt:,308.86,-3.84808,.01998,-5.52177e-$ ? 5,8.53209E-8)-2) cond_expression=.001*POL2(:abt:,115.77,.05139,-1.54476E-4) WHEN Dowtherm_A,10 ------------------------------------------- Dowtherm A saturated liquid ------------------------------------------- Dowtherm is registered trademark of the Dow Chemical Company. rho_expression=POL2(:abt:,1196.98,-.2609,-7.40445e-4) cp_expression=POL2(:abt:,773.2,2.69482,-1.42544e-4) emu_expression=10^(POL4(:abt:,11.3823,-.08682,2.30885e-4, -2.8241E-$ 7,1.29978E-10)-2) cond_expression=.001*POL2(:abt:,174.61,-.11839,1.26737E-6) WHEN Syltherm_800,12 ------------------------------------------- Syltherm 800 saturated liquid ------------------------------------------- Syltherm 800 is registered trademark of Dow Corning Corporation. rho_expression=POL2(:abt:,1143.89,-.57526,-4.50374e-4) cp_expression=POL2(:abt:,1107.79,1.70918,-2.29585e-6) emu_expression=10^(POL4(:abt:,8.82459,-.06548,1.7918e-4, -2.32611E-$ 7,1.14055E-10)-2) cond_expression=.001*POL2(:abt:,189.98,-.18743,7.20079E-6) WHEN FC-72,5 ---------------------- FC-72 saturated liquid ---------------------- FC-72 is registered trademark of 3M. rho_expression=POL2(:abt:,2463.7,-2.69081,1.4985e-4) cp_expression=POL2(:abt:,584.75,1.58946,-6.49351e-5) emu_expression=10^(POL2(:abt:,4.01018,-.02904,3.87596e-5)-2) cond_expression=.001*POL2(:abt:,97.495,-.16494,9.99001E-5) WHEN HFE-7100,8 ------------------------- HFE-7100 saturated liquid ------------------------- HFE-7100 is registered trademark of 3M. rho_expression=2256.6-2.55934*:abt: cp_expression=587.+2.*:abt: emu_expression=10^(POL2(:abt:,3.37646,-.02537,3.3034e-5)-2) ? cond_expression=.001*POL2(:abt:,138.51,-.28616,1.74825) WHEN Mercury,7 --------------------------------------- Mercury [Hg] liquid metal (Tmp = 234 K) --------------------------------------- rho_expression=POL3(:abt:,14.293,-2.68226,5.3957e-4,-3.16674e-7) cp_expression=POL3(:abt:,159.54,-.10108,1.23163e-4,-3.60116e-8) emu_expression=POL3(:abt:,5.47854,-.02372,4.35299e-5,-2.$ 79475E-8) cond_expression=POL3(:abt:,3.90003,.01799,-8.20701E-6,1.52734E-9) ENDCASE unconditional property settings if(:condtn:.eq.t) then ! if 'condition' is simply 'true' IF(:rho_expression:.NE.XXX) THEN (property rho1 is :rho_expression:) ENDIF IF(:emu_expression:.NE.XXX) THEN (property enul is :emu_expression:/rho1) enu_expression=:emu_expression:/rho1 ENDIF IF(:enu_expression:.NE.XXX) THEN (property enul is :enu_expression:) ENDIF IF(:cp_expression:.NE.XXX) THEN (property cp1 is :cp_expression:) ENDIF IF(:cond_expression:.NE.XXX) THEN (stored cond is :cond_expression:) (property prndtl(tem1) is cond/CP1) ENDIF else ! otherwise conditional property settings IF(:rho_expression:.NE.XXX) THEN (property rho1 is :rho_expression: with :condtn:) ENDIF IF(:emu_expression:.NE.XXX) THEN (property enul is :emu_expression:/rho1 with :condtn:) enu_expression=:emu_expression:/rho1 ENDIF IF(:enu_expression:.NE.XXX) THEN (property enul is :enu_expression: with :condtn:) ENDIF IF(:cp_expression:.NE.XXX) THEN (property cp1 is :cp_expression: with :condtn:) ENDIF IF(:cond_expression:.NE.XXX) THEN (stored cond is :cond_expression:) (property prndtl(tem1) is -cond with :condtn:) Note the essential minus above so that prndtl(tem) is recognised as a conductivity ENDIF endif