initial commit / add all books

1 files changed, 265 insertions, 0 deletions

diff --git a/1328/CH13/EX13.6a/13_6a.sce b/1328/CH13/EX13.6a/13_6a.sce
new file mode 100644
index 000000000..1dde88f6f
--- /dev/null
+++ b/1328/CH13/EX13.6a/13_6a.sce
@@ -0,0 +1,265 @@
+printf("\t example 13.6a \n");

+printf("\t approximate values are mentioned in the book \n");

+

+ds=[0 10 20 30 40 50 60 70 80 90 100];

+tmp=[90 145 180 208 234 260 286 312 338 367 400];

+clf();

+subplot(3,2,1);

+plot2d(ds,tmp,style=2,rect=[0,80,100,400]);

+xtitle("Plot of ASTM curve",boxed=1);

+xlabel("Per cent distilled off");

+ylabel("Temperature °F");

+

+//From the plotted ASTM curve and reference line

+s = (312-145)/60; // (70% - 10%)/60%

+printf("\tSlope of ASTm = %.2f °F\n",s);

+ap = (180+260+338)/3; // (20% +50% +80%)/3

+printf("\tAverage 50prcnt point = %.1f °F\n",ap);

+

+fc = 38; //°F, from Fig.13.8

+printf("\t50prcnt point ASTM = 50prcnt point flash curve = %.0f °F\n",fc);

+fc1 = ap - fc; //°F, fixing first point on EFC

+printf("\t50prcnt on EFC = %.0f °F\n",fc1);

+

+s1 = 1.65; // (°F/%) from fig 13.10, upper curve

+ten = 221 - 40*s1; //

+printf("\t10prcnt on EFC = 50prcnt - 40prcnt = %.0f °F\n",ten);

+sty = 221 + 20*s1; //

+printf("\t70prcnt on EFC = 50prcnt + 20prcnt %.0f °F\n",sty);

+

+//Draw this line as a reference through the 50% point. Calculate the flash curve for different percentages off

+

+//0% off

+printf("\n\t0 prcnt off:\n");

+dela = 90 - 117; // Step (8)

+printf("\t\tDelT ASTM = %.0f °F\n",dela); 

+delE = dela * 0.50; // Step (9)

+printf("\t\tDelT EFC = %.1f °F\n",delE);

+FE = 139 - delE; // Step (10)

+printf("\t\t°F EFC = %.1f\n",FE);

+//end

+ov=13300; //lb/hr

+ng=90;//lb/hr

+mng=50;// mol. wt

+st=370;//lb/hr

+avG=50;//°F API

+//For 80%

+ouc=ov*0.80;//lb/hr

+printf("\toil uncondensed = %.0f lb/hr\n",ouc);

+avB=269;//°F,from Fig. 13.13

+printf("\tAverage boiling point from the EFC at 1 atm = %.0f°F\n",avB);

+avB1=avB+17;//°F,from Fig. 13.13

+printf("\tAverage boiling point from the EFC at 19.7 psia = %.0f°F\n",avB1);

+mwt=113;//mol. wt

+mtoc=ouc/mwt;

+printf("\tThe moles of oil still to be condensed = %.1f\n",mtoc);

+mg1=ng/mng;

+ms1=st/18;

+tm=mg1+ms1+mtoc;

+printf("\t\tMols gas = %.2f\n\t\tMols steam = %.1f\n",mg1,ms1);

+printf("\t\t\t    -----\n\t\tMols total = %.1f\n",tm);

+tp=19.7;//psia

+poil=(mtoc/tm)*tp;//psia

+printf("\tPartial pressure of oil = %.1f psia\n",poil);

+pgas=(mg1/tm)*tp;//psia

+printf("\tPartial pressure of NC gas = %.3f psia\n",pgas);

+tm(1)=95;//°F

+tm(2)=127;//°F

+tm(3)=163;//°F

+tm(4)=205;//°F

+tm(5)=240;//°F

+pp(1)=6.73;

+pp(2)=9.40;

+pp(3)=12.25;

+pp(4)=14.64;

+pp(5)=15.65;

+psat(1)=0.815;//From steam table

+psat(2)=2.050;//From steam table

+psat(3)=5.09;//From steam table

+psat(4)=12.77;//From steam table

+psat(5)=24.97;//From steam table

+printf("\n\t\tCALCULATION OF DEW POINT OF THE STEAM\n");

+printf("\tT,°F\t[pt  -  (poil+pNC)]  =  psteam\tpsat(steam tables)\n");

+i=1;

+while(i<6)

+    ps=tp-pp(i);

+    printf("\t"+string(tm(i))+"\t%.1f\t %.2f\t\t%.2f\t%.3f\n",tp,pp(i),ps,psat(i));

+    i=i+1;

+end

+subplot(3,2,2);

+plot2d(psat,tm,style=3,rect=[0,25,90,250]);

+xtitle("Computed pressure of steam",boxed=1);

+xlabel("Pressure of steam, psi");

+ylabel("Temperature °F");

+

+ds=6.417;//psia,at 173°F,

+printf("\tAt 173°F, the dew point of steam, psat = %.3f psia\n",ds);

+pd1=tp-ds;//psia

+printf("\tpoil + pNC = %.2f psia\n",pd1);

+x=((tp*ms1)/ds)-(ms1+mg1);// mols oil

+printf("\tOil = %.2f mols oil\n",x);

+mw=85;//From fig. 13.14

+printf("\tThe molecular weight of the vapors is %.0f\n",mw);

+lv=x*mw;//lb

+printf("\tLb/hr vapor = %.0f\n",lv);

+prc=((ov-lv)*100)/ov;//%

+printf("\tpercent Condensed = %.0f\n",prc);

+printf("\n\t\t\tOIL CONDENSING CURVE\n");

+printf("\tprcnt\tCondensables\t\tAv BP on EFC\t\t50° API\t\tMol oil\t\tMol NC gas\tMol steam\tMol total\tTotal pressure\tPartial pressure\tPartial pressure\tCond temp,°F\n\t\tlb.hr\t\t14.7 psia °F\t19.7 psia,°F\tmol.wt\t\t\t\t\t\t\t\t\t\tpsia\t\toil,psia\t\tNC gas, psia\n");

+mo(1)=107.5;

+mo(2)=94.3;

+mo(3)=77.7;

+mo(4)=57.4;

+mo(5)=31.8;

+mo(6)=17.1;

+mo(7)=8.9;

+i=1;

+while(i<8)

+    mt(i)=mo(i)+mg1+ms1;

+    ppo(i)=(mo(i)/mt(i))*tp;

+    ppg(i)=(mg1/mt(i))*tp;

+    i=i+1;

+end

+printf("\t---------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------\n");

+printf("\t100\t13330\t\t300\t\t317\t\t124\t\t%.1f\t\t1.8\t\t20.6\t\t%.1f\t\t19.7\t\t%.1f\t\t\t%.3f\t\t\t305\n",mo(1),mt(1),ppo(1),ppg(1));

+printf("\t80\t10664\t\t269\t\t286\t\t113\t\t%.1f\t\t1.8\t\t20.6\t\t%.1f\t\t19.7\t\t%.1f\t\t\t%.3f\t\t\t277\n",mo(2),mt(2),ppo(2),ppg(2));

+printf("\t60\t7998\t\t239\t\t256\t\t103\t\t%.1f\t\t1.8\t\t20.6\t\t%.1f\t\t19.7\t\t%.1f\t\t\t%.3f\t\t\t240\n",mo(3),mt(3),ppo(3),ppg(3));

+printf("\t40\t5332\t\t207\t\t224\t\t93\t\t%.1f\t\t1.8\t\t20.6\t\t%.1f\t\t19.7\t\t%.1f\t\t\t%.3f\t\t\t205\n",mo(4),mt(4),ppo(4),ppg(4));

+printf("\t20\t2666\t\t178\t\t195\t\t84\t\t%.1f\t\t1.8\t\t20.6\t\t%.1f\t\t19.7\t\t%.1f\t\t\t%.3f\t\t\t163\n",mo(5),mt(5),ppo(5),ppg(5));

+printf("\t10\t1333\t\t155\t\t172\t\t78\t\t%.1f\t\t1.8\t\t20.6\t\t%.1f\t\t19.7\t\t%.1f\t\t\t%.3f\t\t\t127\n",mo(6),mt(6),ppo(6),ppg(6));

+printf("\t5\t667\t\t141\t\t158\t\t75\t\t%.1f\t\t1.8\t\t20.6\t\t%.1f\t\t19.7\t\t%.1f\t\t\t%.3f\t\t\t95\n",mo(7),mt(7),ppo(7),ppg(7));

+

+//Trail 1:

+m=78;//50° API mol. wt. for condesables 1333

+vap=(ov*0.10)/78;//Mol/hr

+printf("\n\t\t\t\tMol/hr\n\tOil vapor\t\t%.1f\n\tNC gas\t\t\t%.1f\n\tSteam\t\t\tX\n\tTotal\t\t\t18.9+X\n",vap,mg1);

+vap1=vap+mg1;//Mol/hr

+psteam=5.09;//psia, For 163°F

+x1=(psteam*vap1)/(tp-psteam);//mols steam

+printf("\tX = %.2f mols steam\n",x1);

+tv=vap1+x1;

+printf("\n\t\t\tMol/hr\tmf\tmf*pt = p-partial\n");

+mf1=vap/(tv);

+ppar1=mf1*tp;

+printf("\tOil vapor\t%.1f\t%.3f\t%.2f\n",vap,mf1,ppar1);

+mf2=mg1/tv;

+ppar2=mf2*tp;

+printf("\tNC gas\t\t%.1f\t%.3f\t%.2f\n",mg1,mf2,ppar2);

+mf3=x1/tv;

+ppar3=mf3*tp;

+printf("\tSteam\t\t%.2f\t%.3f\t%.2f\n",x1,mf3,ppar3);

+tot1=vap+mg1+x1;

+tot2=mf1+mf2+mf3;

+tot3=ppar1+ppar2+ppar3;

+printf("\tTotal\t\t%.2f\t%.3f\t%.2f\n",tot1,tot2,tot3);

+//Error was found. So trail 2 is done in a similar way

+printf("\n\tSimilarly,\n\tT,°F\tOil cond, prcnt\tOil cond, lb\tSteam cond,lb\n");

+printf("\t173\t74\t\t9863\t\t0\n\t163\t85\t\t11350\t\t204\n\t127\t97.5\t\t13000\t\t357\n\t95\t100\t\t13330\t\t370\n");

+//Condensing curve

+printf("\n\t\t\tOil\t\t\t\tSteam\n\t-----------------------------------------------------------------\n\tTc,°F\tHv,vapor\tHl,liquid\tHg or Hv,\tHl,liquid\n\t\t\t\t\t\tgas or vapor\n");

+printf("\t-----------------------------------------------------------------\n")

+printf("\t305\t368\t\t242\t\t1197.0\t\tSuperheated\n\t277\t359\t\t225\t\t1184.1\t\tSuperheated\n\t240\t337\t\t204\t\t1167.0\t\tSuperheated\n\t205\t322\t\t185\t\t1150.6\t\tSuperheated\n\t173\t310\t\t168\t\t1135.4\t\t140.9\n");//From fig.11 in Appendix and steam tables

+//Heat load

+//305°F:

+hvv=368;

+hvg=1197.0;

+olv=ov*hvv;

+stm=st*hvg;

+ncg=ng*(0.46*273);

+thh=olv+stm+ncg;

+printf("\n\t\t\t\tH\t\tq\n");

+printf("\tOil vapor\t\t%.2e\n\tSteam\t\t\t%.2e\n\tNC gas\t\t\t%.2e\n\t\t\t\t--------\n\t\t\t\t%.4e\t0\n",olv,stm,ncg,thh);

+//Similarily at other temperatures

+ttp(1)=305;//°F

+ttp(2)=277;//°F

+ttp(3)=240//°F

+ttp(4)=205;//°F

+ttp(5)=173;//°F, Dew point of steam

+ttp(6)=163;//°F

+ttp(7)=127;//°F

+ttp(8)=95;//°F

+

+hld(1)=0;//million Btu

+hld(2)=0.55;//milllion Btu

+hld(3)=1.2;//million Btu

+hld(4)=1.75;//million Btu

+hld(5)=2.3;//million Btu

+hld(6)=2.73;//million Btu

+hld(7)=3.3;//million Btu

+hld(8)=3.66;//million Btu

+subplot(2,2,3);

+plot2d(hld,ttp,style=6,rect=[0,60,3.8,320]);

+xtitle("Condensation of mixed hydrocarbons with gas and steam",boxed=1);

+xlabel("Heat load, million Btu");

+ylabel("Temperature °F");

+//summary

+dp=3042800;//Btu/hr

+ttt=3638400;//Btu/hr

+i2s=thh-dp;//Btu/hr

+printf("\tInlet to steam dew point = %.4eBtu/hr\n",i2s);

+so=dp-1735900;//Btu/hr

+printf("\tSteam dew point to outlet = %.4e Btu/hr\n",so);

+totl=i2s+so;//Btu/hr

+printf("\tTotal\t\t\t= %.4e Btu/hr\n",totl);

+twa=ttt/(120-85);

+printf("\tTotal water = %.2e lb/hr\n",twa);

+wt=85+((1306900/ttt)*35);//°F

+printf("\tWater temperature at dew point of steam = %.0f°F\n",wt);

+//Weighted true temperature difference, delT:

+    //Inlet to dew point of steam:

+delq=2331500;

+delt1=122.2;

+UA1=delq/delt1;

+printf("\tUA = %.0f\n",UA1);

+printf("\n\tDew point of steam to oulet\n");

+printf("\tq\tdelq\tTc\ttw\tdelTav\t(delq/delTav) = UA\n");

+printf("\t----------------------------------------------------------\n");

+q(1)=2331500;

+q(2)=2500000;

+q(3)=2750000;

+q(4)=3000000;

+q(5)=3250000;

+q(6)=3500000;

+q(7)=3638000;

+i=1;

+while(i<7)

+    dq(i)=q(i+1)-q(i);

+    i=i+1;

+end

+dpt(1)=173;

+dpt(2)=169;

+dpt(3)=161;

+dpt(4)=149;

+dpt(5)=134;

+dpt(6)=112;

+dpt(7)=95;

+dtw(1)=97.5;

+dtw(2)=96;

+dtw(3)=93;

+dtw(4)=91;

+dtw(5)=89;

+dtw(6)=86;

+dtw(7)=85;

+i=1;

+tua=0;

+while(i<7)

+    dpdelt(i)=((dpt(i+1)-dtw(i+1))+(dpt(i)-dtw(i)))/2;

+    UA(i)=dq(i)/dpdelt(i);

+    tua=tua+UA(i);

+    i=i+1;

+end

+printf("\t2331500\t......\t173\t173\t97.5\n");

+i=1;

+while(i<7)

+    printf("\t"+string(q(i+1))+"\t"+string(dq(i))+"\t"+string(dpt(i+1))+"\t"+string(dtw(i+1))+"\t"+string(dpdelt(i))+"\t%.0f\n",UA(i));//from Fig. 13.16

+i=i+1;

+end

+

+printf("\t\t\t\t\t\t%.0f\tUA = sigma{delq/delt}\n",tua);

+wdt=1306900/tua;//°F

+printf("\tWeighted delt = %.1f°F\n",wdt);

+owdt=ttt/(tua+UA1);//°F

+printf("\tOverall weighted temperature difference = %.1f °F\n",owdt);

+printf("\tThe uncorrected LMTD is 60.1°F\n");

+//end

+