4d9d4542-be88-4fae-b951-efd9003a2f43 213 4 5 decimal

&[DATE] &[TIME] - &[FILENAME]

Test T F If Test T F if :

Π.3TOL 105 -^:

ZR k el Z k el abs Z k el ≡ (Z k el (* :

ZR k el Z k el Im 0 ≡ Z k el 0 If :

Cool Prop Wrapper Π Prop Π Fluid Prop : Π F Prop F strtoupper Π strtoupper CoolProp_Props1 : Π a b Prop Π a b Fluid Prop : Π a b F Prop U K Pa kg m 3^/ J kg / J kg K * / mol m 3^/ J mol / J mol K * / 1 T P d h s dmolar hmolar smolar q 29 mat transpose : A U a UnitsOf 1 findrows 12 el : B U b UnitsOf 1 findrows 12 el : Π A a B b F Prop 41 line : Π P1 V1 P2 V2 F Prop Π strtoupper P1 strtoupper V1 P2 strtoupper V2 F strtoupper CoolProp_Props : End 61 line

n T P Π.3 ν n 11 substr : o n strlen 1 > n 21 substr if : 12 mat Tr T T.c / : κ 0 : α 0 : a R.m T.c * (2^P.c / : b R.m T.c * P.c / : 15 mat Π.3EOS SRK ≡ κ 0.480 1.574 ω * + 0.176 ω 2^* - : α 1 κ 1 Tr sqrt - (* + : 12 mat Π.3EOS PR ≡ κ 0.37464 1.54226 ω * + 0.26992 ω 2^* - : α 1 κ 1 Tr sqrt - (* + : 12 mat κ 0.134 0.508 ω * + 0.0467 ω 2^* - : α e 2 0.836 Tr * +^Tr κ^* sqrt : 12 mat if if Π.3EOS VdW ≡ a 2764 / a * : b b 8 / : 12 mat Π.3EOS RK ≡ a 0.45724 a * T.c K / sqrt * : b 0.07780 b * : 12 mat a 0.45724 α 2^* a * : b 0.08664 b * : 12 mat if if A a P * R.m T * (2^/ : B b P * R.m T * / : 12 mat Π.3EOS VdW ≡ C A B * - A 1 - B - stack : Π.3EOS RK ≡ A A K T / sqrt * : C A B * - A B - B 2^- 1 - stack : 12 mat Π.3EOS SRK ≡ C A B * - A B - B 2^- 1 - stack : C B 3^B 2^+ A B * - A 3 B 2^* - 2 B * - B 1 - stack : if if if ZP C 1 stack polyroots : ZR 0 : Z 0 : 13 mat k 13 range ZP k el abs ZP k el ≡ ZR k el ZP k el : ZR k el 0 : if for ZR min Π.3TOL < Z ZR max ZR max stack : Z ZR min ZR max stack : if 12 mat Π.3EOS VdW ≡ L B Z * vectorize : Π.3EOS RK ≡ (Π.3EOS SRK ≡ (| L 1 B Z / + ln vectorize : L 122 sqrt * / Z 12 sqrt + (B * + Z 12 sqrt - (B * + / ln * vectorize : if if φ Z 1 - Z B - ln - A B / L * - vectorize : ν 0 ≡ Z 1 el B - Z 2 el B - / ln A B / L 2 el L 1 el - Z 2 el Z 1 el - / (* + ι Z rows 1 ≡ 1 o f ≡ 1 o g ≡ 2 φ 1 el φ 2 el < 12 If If If If : Z Z ι el : L L ι el : φ φ ι el : 13 mat ν Z ≡ Z ν φ ≡ e φ^ν f ≡ P e φ^* ν v ≡ R.m Z * T P / * ν d ≡ P R.m Z * T * / da 0.45724 κ * α * R.m 2^* T.c P.c / * 1 Tr / sqrt * - : C CpIG : 12 mat ν h ≡ Π.3EOS VdW ≡ h Z 2 -^: Π.3EOS RK ≡ h 1.5 : Π.3EOS SRK ≡ h 1 κ Tr sqrt * α / + : Π.3EOS PR ≡ h κ Tr sqrt * α / - 1 - : h 1 : if if if if h R.m T * Z 1 - A B / h * L * + (* : o o ≡ h h h.o - kJ kmol / C k el T k^T.o k^- (* k K k^* / (k 1 C rows sum * + if 31 line ν s ≡ s R.m Z B - abs ln * R.m P P.o / ln * - L b / da * + : o o ≡ s s s.o - kJ kmol K * / C 1 el T T.o / ln * C k 1 + el T k^T.o k^- (* k K k^* / (k 1 C rows 1 - sum + (* + if 21 line Not Implemented. if if 21 line if if if if if 31 line if 111 line (:

Z ZR max 11 mat :

X 0.5 11 mat :

X 1 + vectorize 1.5 11 mat

L 122 sqrt * / X 12 sqrt + (1 * + X 12 sqrt - (1 * + / ln * vectorize : 1.2465

BUG

2a3 − 9ab + 27c

α e 2 0.836 Tr * +^Tr κ^* sqrt :

X Π.3 X UnitsOf 1 ≡ f.0 5.92714 6.09648 X / - 1.28862 X ln * - 0.16935 X 6^* + : f.1 15.2518 15.6875 X / - 13.4721 X ln * - 0.43577 X 6^* + : f.0 ω f.1 * + 31 line X UnitsOf K ≡ P P.c e X T.c / Π.3^* : 0 X P.c Π.3 21 line 1 if if 11 line :

ho T.o P.o Π.3

http://people.ds.cam.ac.uk/pjb10/thermo/pure.html

CoolProp.CoolProp.set_reference_state set_reference_state(signatures, args, kwargs, defaults) Accepts one of two signatures: Type #1: set_reference_state(FluidName,reference_state) FluidName The name of the fluid param reference_state The reference state to use, one of IIR (h=200 kJ/kg, s=1 kJ/kg/K at 0C sat. liq.) ASHRAE (h=0,s=0 @ -40C sat liq) NBP (h=0,s=0 @ 1.0 bar sat liq.) Type #2: set_reference_state(FluidName,T0,rho0,h0,s0) FluidName The name of the fluid T0 The temperature at the reference point [K] rho0 The density at the reference point [kg/m^3] h0 The enthalpy at the reference point [J/kg] s0 The entropy at the reference point [J/kg]

Z 1 el B - Z 2 el B - / ln A B / L 2 el L 1 el - Z 2 el Z 1 el - / (* +

B Z 2 el Z 1 el - (* Z 1 el B - det ln Z 2 el B - det ln - (* A L 2 el L 1 el - (* +

10 °C * Π.3 kPa #

Fluid Propane : Π.3EOS PR : 12 mat CpIG 4.224 - 0.3063 0.0001586 - 3.215 108 -^* stack : T.c 96.74 °C * : P.c 4251.2 kPa * : ω 0.1521 : MW 44.096 kg kmol / * : 14 mat T.o 0 °C * : P.o 101 kPa * : h.o MW 200 * kJ kg / * : s.o so T.o P.o Π.3 : 14 mat 41 line

h.o 8819.2 J mol / *

Fluid RESET CoolProp_set_reference_stateS 1

Fluid T.o d T.o P.o Π.3 h.o s.o CoolProp_set_reference_stateD 1

Some TP points

T 3055016041 mat °C * :

P 51004004041 mat kPa * :

Z Z T P Π.3 vectorize : 0.9992 0.9803 0.9469 0.9979 41 mat

T K 303.15 278.15 323.15 433.15 41 mat

P bar 0.05 14 0.4 41 mat

Z T P Prop vectorize 0.9993 0.98 0.949 0.9981 41 mat

0 T 1 el P 1 el Π.3 11.9167 -

φ φ T P Π.3 vectorize : 0.9992 0.9806 0.949 0.9979 41 mat

d d T P Π.3 vectorize : mol m 3^/ 1.9852 44.1108 157.2269 11.1305 41 mat

d T P Prop MW / vectorize mol m 3^/ 1.9852 44.1212 156.8819 11.128 41 mat

h h T P Π.3 vectorize : kJ kmol / 6675.7512 - 8598.1804 - 5534.5364 - 4806.4774 41 mat

s s T P Π.3 vectorize : kJ kmol K * / 32.7319 1.3462 0.4116 46.7074 41 mat

MW h T P Prop * vectorize kJ kmol / 28284.9611 26347.5841 29382.5459 39610.6162 41 mat

MW s T P Prop * vectorize kJ kmol K * / 151.9525 120.5127 119.4759 165.5019 41 mat

h.o ho T P Π.3 vectorize : kJ kmol / 5.3643 - 123.1294 - 406.6307 - 24.9403 - 41 mat

s.o so T P Π.3 vectorize : kJ kmol K * / 24.9794 0.1969 - 12.2669 - 7.6614 41 mat

Fluid T.o d T.o P.o Π.3 h.o s.o CoolProp_set_reference_stateD MW H T T P P Fluid strtoupper CoolProp_Props * vectorize 21 line kJ kmol / 28284.9611 26347.5841 29382.5459 39610.6162 41 mat

Departures

v.f vf T P Π.3 vectorize : m 3^kmol / 0.1162 0.097 6.3602 89.8428 41 mat

v.g vg T P Π.3 vectorize : m 3^kmol / 503.7179 22.6702 6.3602 89.8428 41 mat

Vapor and liquid p hases

h.f hf T P Π.3 vectorize : kJ kmol / 19758.5238 - 23211.3873 - 5534.5364 - 4806.4774 41 mat

h.g hg T P Π.3 vectorize : kJ kmol / 6675.7512 - 8598.1804 - 5534.5364 - 4806.4774 41 mat

s.f sf T P Π.3 vectorize : kJ kmol K * / 55.2119 - 67.1057 - 0.4116 46.7074 41 mat

s.g sg T P Π.3 vectorize : kJ kmol K * / 32.7319 1.3462 0.4116 46.7074 41 mat

Calculus for mixture of phases

Given

x 0.4 :

x v.g v.f - (* v.f + m 3^kmol / 201.5569 9.1263 6.3602 89.8428 41 mat

x h.g h.f - (* h.f + kJ kmol / 14525.4148 - 17366.1045 - 5534.5364 - 4806.4774 41 mat

x s.g s.f - (* s.f + kJ kmol K * / 20.0343 - 39.725 - 0.4116 46.7074 41 mat

volume, enthalpy and entropy of the phase mixture, for a quality of 40%.

Notice that

T T.c < vectorize 111041 mat

and

P P.c < vectorize 111141 mat

so isn´t valid for the last row (T>Tc)

To mass basis

v v MW / : m 3^kg / #

h h MW / : kJ kg / 151.3913 - 194.9878 - 125.5111 - 109.0003 41 mat

s s MW / : kJ kg K * / 0.7423 0.0305 0.0093 1.0592 41 mat

ToDo: Saturation lines, Cp, quality and internal energy:

u.o h.o MW / R.m MW / Z T.o P.o PR * T.o * - : kJ kg / #

u h R.m MW / Z * T * - u.o - : kJ kg / #

Alvaro