Maple Solutions to the Chemical Engineering Problem Set
Maple Solutions to the Chemical Engineering Problem Set
Maple Solutions to the Chemical Engineering Problem Set
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
0.02<br />
0<br />
-0.02<br />
-0.04<br />
-0.06<br />
-0.08<br />
-0.1<br />
0.2<br />
z<br />
0.4<br />
0.6<br />
0.8<br />
1<br />
This illustration clearly shows <strong>the</strong> only real root; <strong>the</strong> two complex roots are near <strong>the</strong> maximum in<br />
<strong>the</strong> curve. Under slightly different conditions <strong>the</strong> curve will cross <strong>the</strong> z axis in three places and we<br />
would have three real roots.<br />
><br />
><br />
Example 2 and 3<br />
Calculate <strong>the</strong> compressibility of ammonia at a reduced pressures of 1, 2, 4, 10, and 20 at a<br />
temperature of 450 K.<br />
The critical properties of ammonia are<br />
> CriticalProps[NH3]:={T[c] = 405.649994, P[c] = 11280000.0}: CriticalProps[NH3];<br />
{ T = 405.649994 ,<br />
P = .112800000 10 }<br />
c<br />
c<br />
8<br />
where <strong>the</strong> temperature is in kelvin and <strong>the</strong> pressure in pascals.<br />
The reduced pressure is related <strong>to</strong> <strong>the</strong> actual pressure by<br />
> prdef:= P[c]=P[r]/P: prdef;<br />
P<br />
r<br />
P =<br />
c P<br />
We will use this relation in <strong>the</strong> substitution <strong>to</strong> follow<br />
The temperature and pressure of interest (in <strong>the</strong> same units) are<br />
> Tspec:=T=450: Tspec;<br />
T = 450<br />
The desired dimensionless pressures are given in a list as<br />
> prlist := [1,2,4,10,20]: prlist;<br />
[ 1, 2, 4, 10, 20]<br />
We compute <strong>the</strong> compressibility at all reduced pressures in one loop<br />
><br />
> for Pr in prlist do<br />
Page 6