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

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

Test T F If Test T F if :

PR.TOL 105 -^:

n T P PR ν n 11 substr : o n strlen 1 > n 21 substr if : 12 mat κ 0.37464 1.54226 ω * + 0.26992 ω 2^* - : α 1 κ 1 T T.c / sqrt - (* + : 12 mat a 0.45724 R.m T.c * α * (2^P.c / * : b 0.07780 R.m T.c * P.c / * : 12 mat A a P * R.m 2^T 2^* / : B b P * R.m T * / : 12 mat Z B 3^B 2^+ A B * - A 3 B 2^* - 2 B * - B 1 - 1 stack polyroots : ZR 0 : 12 mat k 13 range Z k el abs Z k el ≡ ZR k el Z k el : ZR k el 0 : if for Z ZR min PR.TOL < ZR max ZR min if ZR max stack : 12 mat L 122 sqrt * / Z 12 sqrt B * + (+ Z 12 sqrt B * - (+ / ln * vectorize : ϕ Z 1 - Z B - ln - A B / L * - vectorize : 12 mat ι o f ≡ 1 o g ≡ 2 ϕ 1 el ϕ 2 el < 12 If If If : Z Z ι el : L L ι el : 12 mat ν Z ≡ Z ν f ≡ P e L^* ν v ≡ R.m Z * T P / * da 0.45724 κ * α * R.m 2^P.c / * T.c 3^T / sqrt * - : C CpIG : 12 mat ν h ≡ h R.m T * Z 1 - (* T da * a - b / 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 21 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 101 line :

In my version, takes a lot of time:

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 :

Propane

Ideal Gas Cp

CpIG 4.224 - 0.3063 0.0001586 - 3.215 108 -^* stack :

Critical Values

T.c 96.74 °C * : P.c 4251.2 kPa * : ω 0.1521 : 13 mat

MW 44.096 kg kmol / * :

Ref State

T.o 25 °C * : P.o 101 kPa * : 12 mat

h.o ho T.o P.o PR : kJ kmol / 77.0576 -

h.o MW / kJ kg / 1.7475 -

s.o so T.o P.o PR : kJ kmol K * / 0.2 -

s.o MW / kJ kg K * / 0.0045 -

ho and so are stored here, or can choose any value for ho, so

Some TP points

T 3055016041 mat °C * :

P 51004004041 mat kPa * :

Z Z T P PR vectorize : 0.9992 0.98 0.9457 0.9978 41 mat

f f T P PR vectorize : kPa 5.0003 100.123 401.7264 40.0125 41 mat

v v T P PR vectorize : m 3^kmol / 503.7115 22.6632 6.3521 89.8364 41 mat

h h T P PR vectorize : kJ kmol / 445.3855 1440.4083 - 1715.8657 11934.2988 41 mat

s s T P PR vectorize : kJ kmol K * / 26.4187 4.9234 - 5.7599 - 40.4004 41 mat

v.f vf T P PR vectorize : m 3^kmol / 0.0891 0.0802 0.1009 89.8364 41 mat

v.g vg T P PR vectorize : m 3^kmol / 503.7115 22.6632 6.3521 89.8364 41 mat

Vapor and liquid p hases

h.f hf T P PR vectorize : kJ kmol / 2074.4351 - 3736.1376 - 903.5222 - 11934.2988 41 mat

h.g hg T P PR vectorize : kJ kmol / 445.3855 1440.4083 - 1715.8657 11934.2988 41 mat

s.f sf T P PR vectorize : kJ kmol K * / 53.7259 - 61.9248 - 47.0439 - 40.4004 41 mat

s.g sg T P PR vectorize : kJ kmol K * / 26.4187 4.9234 - 5.7599 - 40.4004 41 mat

Calculus for mixture of phases

Given

x 0.4 :

x v.g v.f - (* v.f + m 3^kmol / 201.5381 9.1134 2.6014 89.8364 41 mat

x h.g h.f - (* h.f + kJ kmol / 1066.5068 - 2817.8459 - 144.2329 11934.2988 41 mat

x s.g s.f - (* s.f + kJ kmol K * / 21.6681 - 39.1242 - 30.5303 - 40.4004 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 / 11.4231 0.514 0.1441 2.0373 41 mat

h h MW / : kJ kg / 10.1004 32.6653 - 38.912 270.6436 41 mat

s s MW / : kJ kg K * / 0.5991 0.1117 - 0.1306 - 0.9162 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 / 57.0239 -

u h R.m MW / Z * T * - u.o - : kJ kg / 180.4982 - 223.2638 - 151.6865 - 80.045 41 mat

Alvaro