> restart;#'therm_eq6.mws'.13-04-2004.

Chapter 6: Set of pre-written equations available.

Basic Clapeyron's equation

> eq6_1:=Diff(p,T)[sat]=h[lv]/(T*v[lv]);

`:=`(eq6_1, (Diff(p, T))[sat] = `/`(`*`(h[lv]), `*`(T, `*`(v[lv]))))

Integrated Clapeyron's equation (p0,T0 is a point in the curve)

> eq6_4:=p=p[0]*exp(-(h[lv]/R)*(1/T-1/T[0]));

`:=`(eq6_4, p = `*`(p[0], `*`(exp(`+`(`-`(`/`(`*`(h[lv], `*`(`+`(`/`(1, `*`(T)), `-`(`/`(1, `*`(T[0])))))), `*`(R))))))))

Antoine's eqaution (pu,Tu are dimensional units)

> eq6_41:=p=p[u]*exp(A-B/(C+T/T[u]));eq6_42:=T=T[u]*(B/(A-ln(p/p[u]))-C);

`:=`(eq6_41, p = `*`(p[u], `*`(exp(`+`(A, `-`(`/`(`*`(B), `*`(`+`(C, `/`(`*`(T), `*`(T[u])))))))))))

`:=`(eq6_42, T = `*`(T[u], `*`(`+`(`/`(`*`(B), `*`(`+`(A, `-`(ln(`/`(`*`(p), `*`(p[u]))))))), `-`(C)))))

Vapour mass fraction (vapour quality)

> eq6_7:=x=m[v]/(m[l]+m[v]);

`:=`(eq6_7, x = `/`(`*`(m[v]), `*`(`+`(m[l], m[v]))))

Heat input needed to sustain a steady mass outflow

> eq6_8:=Q=m*h[lv]*(1+v[l]/v[lv]);

`:=`(eq6_8, Q = `*`(m, `*`(h[lv], `*`(`+`(1, `/`(`*`(v[l]), `*`(v[lv])))))))

Mean specific volume (lever rule for a biphasic state)

> eq6_9:=v=(1-x)*v[l]+x*v[v];

`:=`(eq6_9, v = `+`(`*`(`+`(1, `-`(x)), `*`(v[l])), `*`(x, `*`(v[v]))))

WARNING. A list of all variables follows, to copy and paste to the save command (after Maple7 there is no saving all). CAUTION: After pasting all, I must MANUALLY delete system variables: like RealRange_...

> sort([anames()]);

[TypeTools, eq6_1, eq6_4, eq6_41, eq6_42, eq6_7, eq6_8, eq6_9, exp, ln, `print/Diff`, `tools/membertype`, `type/SymbolicInfinity`, `type/arctrig`]
[TypeTools, eq6_1, eq6_4, eq6_41, eq6_42, eq6_7, eq6_8, eq6_9, exp, ln, `print/Diff`, `tools/membertype`, `type/SymbolicInfinity`, `type/arctrig`]

> save eq6_1, eq6_4, eq6_41, eq6_42, eq6_7, eq6_8, eq6_9, "../therm_eq6.m":

>