diff options
author | priyanka | 2015-06-24 15:03:17 +0530 |
---|---|---|
committer | priyanka | 2015-06-24 15:03:17 +0530 |
commit | b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b (patch) | |
tree | ab291cffc65280e58ac82470ba63fbcca7805165 /3417/CH13/EX13.4.3/Ex13_4_3.sce | |
download | Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.tar.gz Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.tar.bz2 Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.zip |
initial commit / add all books
Diffstat (limited to '3417/CH13/EX13.4.3/Ex13_4_3.sce')
-rw-r--r-- | 3417/CH13/EX13.4.3/Ex13_4_3.sce | 58 |
1 files changed, 58 insertions, 0 deletions
diff --git a/3417/CH13/EX13.4.3/Ex13_4_3.sce b/3417/CH13/EX13.4.3/Ex13_4_3.sce new file mode 100644 index 000000000..e42ddea2b --- /dev/null +++ b/3417/CH13/EX13.4.3/Ex13_4_3.sce @@ -0,0 +1,58 @@ +//Ex.13.4.3;maximum efficiency,no. of thermocouple in series,open ckt voltage,heat i/p and reject at full load.
+
+kA=0.02;//unit=watt/cm degree kelvin
+kB=0.03;//unit=watt/cm degree kelvin
+pA=0.01;//unit=ohm cm
+pB=0.012;//unit=ohm cm
+TH=1500;//unit=degree kelvin
+TC=1000;//unit=degree kelvin
+AA=43.5;//unit=cm^2
+AB=48.6;//unit=cm^2
+LA=0.49;//unit=cm
+LB=0.49;//unit=cm
+I=20*48.6;//Current density in the element limited to,I=20 amp/cm^2
+output=100;//unit=kW
+//alpha_SAB at 1250 degree kelvin=0.0012 volt/degree kelvin=alpha_SA-alpha_SB
+alpha_SAB=0.0012;//unit=volt/degree kelvin
+//let
+b=(pA*kA);
+c=(pB*kB);
+A=sqrt(b);
+B=sqrt(c);
+C=(A+B);
+//figure of merit
+Z=(alpha_SAB/C)^2;
+printf(" Z=%f degree k^-1",Z);
+M=(1+((Z/2)*(TH+TC)))^0.5;
+printf("\n M=%f",M);
+//let
+aa=((TH-TC)/TH);
+bb=(M-1)/(M+(TC/TH));
+//1] MAx. efficiency of a thermoelectric converter is given by n_max=((TH-TC)/TH)*[(M-1)/(M+(TC/TH))]*100;
+n_max=aa*bb*100;
+printf("\n Maximum efficiency n_max=%f persent",n_max);
+//2] No. of thermocouple in series
+V=alpha_SAB*(TH-TC);
+printf("\n V=%f volt",V);
+R=((pA*LA)/AA)+((pB*LB)/AB);//since R=RA+RB=((pA*LA)/AA)+((pB*LB)/AB);
+printf("\n R=%f ohm",R);
+VL=V-(R*I);
+printf("\n VL=%f volt",VL);
+//NTCS=total voltage required/voltage required by one couple
+NTCS=115/VL;
+printf("\n No. of thermocouple in series=%f",NTCS);
+//3] Open circuit voltage
+OCV=V*309;
+printf("\n Open circuit voltage=%f volt",OCV)
+//4] Heat input and reject at full load.
+//Heat input at full load.=output/efficency=100/0.091
+HIFL=output/(n_max/100);
+printf("\n Heat input at full load=%f kW",HIFL)
+// Heat reject at full load. =Heat input-Work output
+HRFL=HIFL-output;
+printf("\n Heat reject at full load=%f kW",HRFL)
+
+
+
+//The value of "pB" is misprinted
+//The values are taken in the text book is approximately equal to calculated values
|