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 /551/CH4 | |
| 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 '551/CH4')
62 files changed, 1427 insertions, 0 deletions
diff --git a/551/CH4/EX4.1/1.sce b/551/CH4/EX4.1/1.sce new file mode 100755 index 000000000..82b6ad3e1 --- /dev/null +++ b/551/CH4/EX4.1/1.sce @@ -0,0 +1,7 @@ +clc
+Q=-50; //kJ/kg
+W=-100; //kJ/kg
+dU=Q-W;
+disp("gain in internal energy = ")
+disp(dU)
+disp("kJ/kg")
\ No newline at end of file diff --git a/551/CH4/EX4.10/10.sce b/551/CH4/EX4.10/10.sce new file mode 100755 index 000000000..825365ffd --- /dev/null +++ b/551/CH4/EX4.10/10.sce @@ -0,0 +1,26 @@ +clc
+Q_lqm=168; //kJ
+W_lqm=64; //kJ
+dU_lm=Q_lqm - W_lqm;
+W_lnm=21; //kJ
+W_ml=-42; //kJ
+
+Q_lnm=dU_lm + W_lnm;
+disp("(i)Q_lnm=")
+disp(Q_lnm)
+disp("kJ")
+
+
+Q_ml=W_ml - dU_lm;
+disp("(ii)Q_ml = ")
+disp(Q_ml)
+disp("kJ")
+
+
+W_ln=21; //kJ
+dU_ln=84; //kJ
+Q_ln=dU_ln + W_ln;
+Q_nm=Q_lnm-Q_ln;
+disp("Q_nm = ")
+disp(Q_nm)
+disp("kJ")
\ No newline at end of file diff --git a/551/CH4/EX4.11/11.sce b/551/CH4/EX4.11/11.sce new file mode 100755 index 000000000..9ba496bc0 --- /dev/null +++ b/551/CH4/EX4.11/11.sce @@ -0,0 +1,9 @@ +clc
+T1=55; //0C
+T2=95; //0C
+W=integrate('200', 'T', T1, T2);
+Q=integrate('160', 'T', T1, T2);
+dU=Q-W;
+disp("change in internal energy=")
+disp(dU/10^3)
+disp("kJ")
\ No newline at end of file diff --git a/551/CH4/EX4.12/12.sce b/551/CH4/EX4.12/12.sce new file mode 100755 index 000000000..684761e69 --- /dev/null +++ b/551/CH4/EX4.12/12.sce @@ -0,0 +1,49 @@ +clc
+Q=-340; //kJ
+n=200; //cycles/min
+
+//For Process 1-2
+
+W_12=4340; //kJ/min
+Q_12=0;
+
+dE_12=Q_12-W_12;
+disp("dE_12 =")
+disp(dE_12)
+disp("kJ/min")
+
+//For process 2-3
+
+Q_23=42000; //kJ/min
+W_23=0;
+
+dE_23=Q_23-W_23;
+disp("dE_23 =")
+disp(dE_23)
+disp("kJ/min")
+
+//For process 3-4
+
+Q_34=-4200; //kJ/min
+dE_34=-73200; //kJ/min
+
+W_34=Q_34-dE_34;
+disp("W_34 =")
+disp(W_34)
+disp("kJ/min")
+
+//For process 4-1
+Q_41=Q*n-Q_12-Q_23-Q_34;
+disp("Q_41 =")
+disp(Q_41)
+disp("kJ/min")
+
+dE_41=0-dE_12-dE_23-dE_34;
+disp("dE_41 =")
+disp(dE_41)
+disp("kJ/min")
+
+W_41=Q_41-dE_41;
+disp("W_41 =")
+disp(W_41)
+disp("kJ/min")
\ No newline at end of file diff --git a/551/CH4/EX4.13/13.sce b/551/CH4/EX4.13/13.sce new file mode 100755 index 000000000..87156207b --- /dev/null +++ b/551/CH4/EX4.13/13.sce @@ -0,0 +1,11 @@ +clc
+P=1200; //kW
+Qin=3360; //kJ/kg
+Qout=2520; //kJ/kg
+F=6; //kW
+dQ=Qin - Qout;
+dW=P-F; //kJ/s
+m = dW/dQ;
+disp("Steam flow round the cycle")
+disp(m)
+disp("kg/s")
\ No newline at end of file diff --git a/551/CH4/EX4.14/14.sce b/551/CH4/EX4.14/14.sce new file mode 100755 index 000000000..c7ad70abd --- /dev/null +++ b/551/CH4/EX4.14/14.sce @@ -0,0 +1,16 @@ +clc
+dT=25; //0C
+Q=30; //kJ
+cv=1.2; //kJ/kg.0C
+m=2.5; //kg
+dU=m*cv*dT;
+
+
+disp("change in internal energy = ")
+disp(dU)
+disp("kJ")
+
+W=Q - dU;
+disp("Work done = ")
+disp(W)
+disp("kJ")
\ No newline at end of file diff --git a/551/CH4/EX4.15/15.sce b/551/CH4/EX4.15/15.sce new file mode 100755 index 000000000..f46b046b1 --- /dev/null +++ b/551/CH4/EX4.15/15.sce @@ -0,0 +1,31 @@ +clc
+Q=50; //kJ
+dV=0.14; //m^3
+p=1.2*10^5; //N/m^2
+m=90; //kg
+d=5.5; //m
+g=9.8; //m/s^2
+W_adb = -110; //kJ
+Wnet=m*g*d/1000; //kJ
+
+
+disp("(i)Change in internal energy")
+W=p*dV/1000 + Wnet; //kJ
+dE=Q-W;
+disp("dE=")
+disp(dE)
+disp("kJ")
+
+disp("(ii) Adiabatic process")
+Q=0;
+dE=-W_adb;
+disp("dE=")
+disp(dE)
+disp("kJ")
+
+disp("(iii)Change in internal energy")
+Q=50; //kJ
+dE=Q - [W_adb+W];
+disp("dE=")
+disp(dE)
+disp("kJ")
\ No newline at end of file diff --git a/551/CH4/EX4.16/16.sce b/551/CH4/EX4.16/16.sce new file mode 100755 index 000000000..5a54bf7c9 --- /dev/null +++ b/551/CH4/EX4.16/16.sce @@ -0,0 +1,21 @@ +clc
+V1=0.15; //m^3
+V2=0.05; //m^3
+Q=-45; //kJ
+p1=(5/V1+1.5)*10^5; //N/m^2
+p2=(5/V2+1.5)*10^5; //N/m^2
+W=integrate('(5/V+1.5)*10^2', 'V', V1, V2);
+
+
+disp("(i)Change in internal energy = ")
+dU=Q-W;
+disp("dU=")
+disp(dU)
+disp("kJ")
+
+
+disp("(ii) Change in enthalpy")
+dH=(dU*10^3+(p2*V2-p1*V1))/10^3;
+disp("dH=")
+disp(dH)
+disp("kJ")
\ No newline at end of file diff --git a/551/CH4/EX4.17/17.sce b/551/CH4/EX4.17/17.sce new file mode 100755 index 000000000..4b53e4a70 --- /dev/null +++ b/551/CH4/EX4.17/17.sce @@ -0,0 +1,27 @@ +clc
+V1=0.25; //m^3
+p1=500; //kPa
+p2=100; //kPa
+V2=V1*(p1/p2)^(1/1.25)
+n=1.25
+dU=3.64*(p2*V2 - p1*V1)
+
+
+disp("(i) If the expansion is quasi-static")
+W=(p1*V1-p2*V2)/(n-1);
+Q=dU+W
+disp("Heat transfered=")
+disp(Q)
+disp("kJ")
+
+
+disp("(ii) In another process")
+Q=32; //kJ
+W=Q-dU;
+disp("Work done=")
+disp(W)
+disp("kJ")
+
+
+disp("(iii)The difference")
+disp("(iii) The work in (ii) is not equal to ∫ p dV since the process is not quasi-static.")
\ No newline at end of file diff --git a/551/CH4/EX4.18/18.sce b/551/CH4/EX4.18/18.sce new file mode 100755 index 000000000..3f3709894 --- /dev/null +++ b/551/CH4/EX4.18/18.sce @@ -0,0 +1,33 @@ +clc
+v1=0.3; //m^3/kg
+T1=20; //0C
+v2=0.55; //m^3/kg
+T2=260; //0C
+p=1.6*10^5; //Pa
+
+disp("(i)Heat added per kg = ")
+Q=integrate('1.5 + 75/(T+45)', 'T', T1,T2);
+disp("Q=")
+disp(Q)
+disp("kJ/kg")
+
+
+disp("(ii)The work done per kg of fluid")
+W=p*(v2-v1)/1000; //kJ/kg
+disp("W=")
+disp(W)
+disp("kJ/kg")
+
+
+disp("(iii)Change in internal energy")
+dU=Q-W;
+disp("dU=")
+disp(dU)
+disp("kJ/kg")
+
+
+disp("(iv)Change in enthalpy")
+dH=Q;
+disp("dH=")
+disp(dH)
+disp("kJ/kg")
\ No newline at end of file diff --git a/551/CH4/EX4.19/19.sce b/551/CH4/EX4.19/19.sce new file mode 100755 index 000000000..3f14a2dcf --- /dev/null +++ b/551/CH4/EX4.19/19.sce @@ -0,0 +1,11 @@ +clc
+m=1; //kg
+du=-42000; //J
+cp=840; //J/kg.0C
+cv=600; //J/kg.0C
+dT=du/m/cv;
+Q=m*cp*dT;
+W=(Q-du)/10^3;
+disp("Work done=")
+disp(W)
+disp("kJ")
\ No newline at end of file diff --git a/551/CH4/EX4.2/2.sce b/551/CH4/EX4.2/2.sce new file mode 100755 index 000000000..49ac83698 --- /dev/null +++ b/551/CH4/EX4.2/2.sce @@ -0,0 +1,8 @@ +clc
+u1=450; //kJ/kg
+u2=220; //kJ/kg
+W=120; //kJ/kg
+Q=(u2-u1) + W;
+disp("Heat rejected by air =")
+disp(-Q)
+disp("kJ/kg")
\ No newline at end of file diff --git a/551/CH4/EX4.20/20.sce b/551/CH4/EX4.20/20.sce new file mode 100755 index 000000000..f6fe8e3d8 --- /dev/null +++ b/551/CH4/EX4.20/20.sce @@ -0,0 +1,25 @@ +clc
+p1=190; //kPa
+V1=0.035; //m^3
+p2=420; //kPa
+V2=0.07; //m^3
+dU=3.6*(p2*V2-p1*V1);
+// p=a+b*v
+// a+0.035*b=190
+// a+0.07*b=420
+// solving the two equations
+p=[1,0.035;1,0.07];
+q=[190;420];
+X=inv(p)*q;
+a=X(1,1);
+b=X(2,1);
+W=integrate('a+b*V', 'V', V1, V2);
+disp("Work done by the system = ")
+disp(W)
+disp("kJ")
+
+
+Q=dU+W;
+disp("Heat transfer into the system = ")
+disp(Q)
+disp("kJ")
\ No newline at end of file diff --git a/551/CH4/EX4.21/21.sce b/551/CH4/EX4.21/21.sce new file mode 100755 index 000000000..7704f51a5 --- /dev/null +++ b/551/CH4/EX4.21/21.sce @@ -0,0 +1,21 @@ +clc
+Qv=90; //kJ
+Qp=-95; //kJ
+W=-18; //kJ
+U_l=105; //kJ
+W_lm=0;
+Q_lm=90;
+U_m=U_l+90;
+dU_mn=Qp-W;
+U_n=U_m+dU_mn;
+dQ=Qv+Qp;
+dW=dQ;
+W_nl=dW-W;
+disp("W_nl(in kJ)=")
+disp(W_nl)
+disp("U_l in kJ =")
+disp(U_l)
+disp("U_m in kJ =")
+disp(U_m)
+disp("U_n in kJ")
+disp(U_n)
\ No newline at end of file diff --git a/551/CH4/EX4.23/23.sce b/551/CH4/EX4.23/23.sce new file mode 100755 index 000000000..2396a1171 --- /dev/null +++ b/551/CH4/EX4.23/23.sce @@ -0,0 +1,29 @@ +clc
+V1=0.2; //m^3
+p1=4*10^5; //N/m^2
+T1=403; //K
+p2=1.02*10^5; //N/m^2
+dH=72.5; //kJ
+Q_23=dH;
+cp=1; //kJ/kg
+cv=0.714; //kJ/kg
+y=1.4;
+V2=V1*(p1/p2)^(1/y);
+T2=T1*((p2/p1)^((y-1)/y));
+R=(cp-cv)*1000; //J/kg.K
+m=p1*V1/R/T1;
+T3=Q_23/(m*cp) +T2;
+V3=V2*T3/T2;
+W_12=(p1*V1 - p2*V2)/(y-1);
+W_23=p2*(V3-V2);
+W_123=W_12+W_23;
+disp("Total work done = ")
+disp(W_123)
+disp("J")
+
+
+disp("(ii) Index of expansion, n")
+p3=p2;
+n=(p1*V1-p3*V3)/W_123 + 1;
+disp("value of index = ")
+disp(n)
\ No newline at end of file diff --git a/551/CH4/EX4.25/25.sce b/551/CH4/EX4.25/25.sce new file mode 100755 index 000000000..278dce86c --- /dev/null +++ b/551/CH4/EX4.25/25.sce @@ -0,0 +1,20 @@ +clc
+d=0.15; //m
+T=303; //K
+p=3*10^5; //N/m^2
+l=0.085; //m
+Q=-4000; //J
+
+disp("(i) Workdone by the system")
+dv=%pi/4*d^2*l;
+W=p*dv;
+disp("W=")
+disp(W/10^3)
+disp("kJ")
+
+
+disp("(ii) Decrease in internal energy of the system")
+dU=(Q-W)/10^3;
+disp("Decrease in internal energy = ")
+disp(-dU)
+disp("kJ")
\ No newline at end of file diff --git a/551/CH4/EX4.27/27.sce b/551/CH4/EX4.27/27.sce new file mode 100755 index 000000000..0ca9be4c4 --- /dev/null +++ b/551/CH4/EX4.27/27.sce @@ -0,0 +1,37 @@ +clc
+y=1.4
+R=294.2; //J/kg.0C
+p1=1*10^5; //N/m^2
+T1=353; //K
+V1=0.45; //m^3
+V2=0.13; //m^3
+p2=5*10^5; //N/m^2
+cv=R/(y-1);
+
+disp("(i) The mass of gas")
+m=p1*V1/R/T1;
+disp("m=")
+disp(m)
+disp("kg")
+
+
+disp("(ii) The value of index ‘n’ for compression")
+n=log(p2/p1)/log(V1/V2);
+disp("n=")
+disp(n)
+
+
+disp("(iii) The increase in internal energy of the gas")
+T2=T1*(V1/V2)^(n-1);
+dU=m*cv*(T2-T1)/10^3;
+disp("dU=")
+disp(dU)
+disp("kJ")
+
+
+disp("(iv) The heat received or rejected by the gas during compression.")
+W=m*R*(T1-T2)/(n-1)/10^3;
+Q=dU+W;
+disp("Q=")
+disp(Q)
+disp("kJ")
\ No newline at end of file diff --git a/551/CH4/EX4.28/28.sce b/551/CH4/EX4.28/28.sce new file mode 100755 index 000000000..9e9f9002c --- /dev/null +++ b/551/CH4/EX4.28/28.sce @@ -0,0 +1,30 @@ +clc
+p1=1.02*10^5; //Pa
+T1=295; //K
+V1=0.015; //m^3
+p2=6.8*10^5; //Pa
+y=1.4;
+
+
+disp("(i) Final temperature")
+T2=T1*(p2/p1)^((y-1)/y);
+t2=T2-273;
+disp("t2=")
+disp(t2)
+disp("°C")
+
+
+disp("(ii) Final volume :")
+V2=V1*(p1/p2)^(1/y);
+disp("V2=")
+disp(V2)
+disp("m^3")
+
+
+disp("(iii)Work done")
+R=287;
+m=p1*V1/R/T1;
+W=m*R*(T1-T2)/(y-1)/10^3;
+disp("W=")
+disp(W)
+disp("kJ")
\ No newline at end of file diff --git a/551/CH4/EX4.29/29.sce b/551/CH4/EX4.29/29.sce new file mode 100755 index 000000000..c94f4d613 --- /dev/null +++ b/551/CH4/EX4.29/29.sce @@ -0,0 +1,24 @@ +clc
+m=0.44; //kg
+T1=453; //K
+ratio=3; //ratio=V2/V1
+T2=288; //K
+W_12=52.5; //kJ
+y=log(T2/T1)/ log(1/ratio) + 1;
+R=W_12*(y-1)/m/(T1-T2);
+// We have got two equations
+// cp-cv=R
+// cp-y*cv=0
+M=[1,-1;1,-y];
+N=[R;0];
+X=inv(M)*N;
+cp=X(1,1);
+cv=X(2,1);
+disp("cp=")
+disp(cp)
+disp("kJ/kg.K")
+
+
+disp("cv=")
+disp(cv)
+disp("kJ/kg.K")
\ No newline at end of file diff --git a/551/CH4/EX4.3/3.sce b/551/CH4/EX4.3/3.sce new file mode 100755 index 000000000..b88baa3a9 --- /dev/null +++ b/551/CH4/EX4.3/3.sce @@ -0,0 +1,11 @@ +clc
+m=0.3; //kg
+cv=0.75; //kJ/kg.K
+T1=313; //K
+T2=433; //K
+W=-30; //kJ
+dU=m*cv*(T2-T1);
+Q=dU + W;
+disp("Heat rejected during the process=")
+disp(-Q)
+disp("kJ")
\ No newline at end of file diff --git a/551/CH4/EX4.30/30.sce b/551/CH4/EX4.30/30.sce new file mode 100755 index 000000000..1499b08a1 --- /dev/null +++ b/551/CH4/EX4.30/30.sce @@ -0,0 +1,18 @@ +clc
+n=1.3;
+m=1; //kg
+p1=1.1; //bar
+T1=300; //K
+p2=6.6; //bar
+R0=8314;
+M=30;
+cp=1.75; //kJ/kg.K
+R=R0/M/1000; //kJ/kg.K
+cv=cp - R;
+y=cp/cv;
+T2=T1 *(p2/p1)^((n-1)/n);
+W=R*(T1-T2)/(n-1);
+Q=((y-n)/(y-1))*W;
+disp("Heat supplied = ")
+disp(Q)
+disp("kJ/kg")
\ No newline at end of file diff --git a/551/CH4/EX4.31/31.sce b/551/CH4/EX4.31/31.sce new file mode 100755 index 000000000..7ecaad192 --- /dev/null +++ b/551/CH4/EX4.31/31.sce @@ -0,0 +1,42 @@ +clc
+cp=14.3; //kJ/kg.K
+cv=10.2; //kJ/kg.K
+V1=0.1; //m^3
+T1=300; //K
+p1=1; //bar
+p2=8; //bar
+y=cp/cv;
+R=cp-cv;
+V2=V1*(p1/p2)^(1/y);
+V3=V2;
+T2=T1*(p2/p1)^((y-1)/y);
+p3=p1*V1/V3;
+T3=300; //K
+
+
+disp("(i) Pressure at the end of constant volume cooling = ")
+disp(p3)
+disp("bar")
+
+
+disp("(ii) Change in internal energy during constant volume process")
+m=p1*V1/R/T1*10^2; //kg
+
+dU_23=m*cv*(T3-T2);
+disp("dU_23 = ")
+disp(dU_23)
+disp("kJ")
+
+
+disp("(iii) Net work done and heat transferred during the cycle")
+W_12=m*R*(T1-T2)/(y-1);
+W_23=0;
+W_31=p3*V3*log(V1/V3)*10^2; //kJ
+disp("Net work done=")
+Wnet=W_12+W_23+W_31;
+disp(Wnet)
+disp("kJ")
+Qnet=Wnet;
+disp("Heat transferred during the complete cycle = ")
+disp(Qnet)
+disp("kJ")
\ No newline at end of file diff --git a/551/CH4/EX4.32/32.sce b/551/CH4/EX4.32/32.sce new file mode 100755 index 000000000..b97de6dd6 --- /dev/null +++ b/551/CH4/EX4.32/32.sce @@ -0,0 +1,24 @@ +clc
+V1=0.15; //m^3
+p1=15; //bar
+T1=550; //K
+T2=T1;
+r=4; //r=V2/V1
+V2=r*V1;
+T3=290; //K
+p2=p1*V1/V2;
+W_12=p1*V1*log(V2/V1)*10^2; //kJ
+V3=V2;
+p3=p2*T3/T2;
+W_23=0;
+n=log(p1/p3)/log(V3/V1);
+W_31=(p3*V3-p1*V1)/(n-1)*10^2; //kJ
+
+disp("net work done = ")
+Wnet=W_12+W_23+W_31
+disp("kJ")
+
+Qnet=Wnet;
+disp("Heat transferred during the cycle = ")
+disp(Qnet)
+disp("kJ")
\ No newline at end of file diff --git a/551/CH4/EX4.33/33.sce b/551/CH4/EX4.33/33.sce new file mode 100755 index 000000000..7ad63df02 --- /dev/null +++ b/551/CH4/EX4.33/33.sce @@ -0,0 +1,71 @@ +clc
+m=1; //kg
+p1=5; //bar
+V1=0.02; //m^3
+V2=0.08; //m^3
+p2=1.5; //bar
+
+function p=f(V)
+ p=a+b*V;
+endfunction
+
+// 5=a+0.02*b
+// 1.5=a+0.08*b
+// Solving above two equations
+
+A=[1,0.02;1,0.08];
+B=[5;1.5];
+X=inv(A)*B;
+a=X(1,1);
+b=X(2,1);
+
+disp("(i) p-V diagram")
+
+V=0.02:0.001:0.08;
+p=a+b*V;
+plot(V,p,'b')
+
+V=[0.0667 0.08];
+p=[1.5 1.5];
+plot(V,p,'g')
+
+V=0.02:0.001:0.0667;
+function p=fa(V)
+ p=0.1/V;
+endfunction
+
+plot(V,fa,'r')
+
+V=[0.0667 0.0667];
+p=[1.5 0];
+plot(V,p,'--')
+
+xtitle("p-V diagram", "V(m^3)", "p(bar)");
+legend("p=a+b*V","p=constant","pv=constant")
+
+
+disp("(ii) Work done and heat transfer")
+
+W_12=integrate('(a+b*V)*10^2','V',V1,V2);
+disp("Work done by the system =")
+disp(W_12)
+disp("kJ")
+
+p3=p2;
+V3=p1*V1/p3;
+W_23=p2*(V3-V2)*10^2; //kJ
+
+W_31=p3*V3*log(V1/V3)*10^2; //kJ
+disp("Work done on the system =")
+disp(W_31)
+disp("kJ")
+
+W_net=W_12+W_23+W_31;
+disp("Net work done =")
+disp(W_net)
+disp("kJ")
+
+Q_net=W_net;
+disp("Heat transferred during the complete cycle =")
+disp(Q_net)
+disp("kJ")
\ No newline at end of file diff --git a/551/CH4/EX4.34/34.sce b/551/CH4/EX4.34/34.sce new file mode 100755 index 000000000..a60573e82 --- /dev/null +++ b/551/CH4/EX4.34/34.sce @@ -0,0 +1,24 @@ +clc
+cv=0.71; //kJ/kg.K
+R=0.287; //kJ/kg.K
+d=8; //cm
+l=3.5; //cm
+S=150; //N/cm
+p1=30; //N/cm
+V1=45; //cm^3
+T1=293; //K
+cv=0.71; //kJ/kg.K
+R=0.287; //kJ/kg.K
+A=%pi/4*d^2;
+C=p1-S/A^2*V1;
+dV=l*A;
+V2=V1+dV;
+p2=S/A^2*V2 + C;
+W=integrate('A^2/S*p/100', 'p', p1, p2);
+T2=p2*V2*T1/p1/V1;
+m=p1*V1/R/T1/10^5; //kg
+dU=m*cv*(T2-T1);
+Q_12=dU + W*10^(-3);
+disp("Amount of heat added to the system = ")
+disp(Q_12)
+disp("kJ")
\ No newline at end of file diff --git a/551/CH4/EX4.35/35.sce b/551/CH4/EX4.35/35.sce new file mode 100755 index 000000000..5f8edaff2 --- /dev/null +++ b/551/CH4/EX4.35/35.sce @@ -0,0 +1,36 @@ +clc
+r=10; //kg/min
+p1=1.5*10^5; //N/m^2
+rho1=26; //kg/m^3
+C1=110; //m/s
+u1=910; //kJ/kg
+p2=5.5*10^5; //N/m^2
+rho2=5.5; //kg/m^3
+C2=190; //m/s
+u2=710; //kJ/kg
+Q=55; //kJ/s
+h=55; //m
+g=9.81; //m/s^2
+v2=1/rho2;
+v1=1/rho1;
+
+disp("(i) Change in enthalpy")
+dh=u2-u1+ (p2*v2-p1*v1)/10^3;
+disp(dh)
+disp("kJ/kg")
+
+disp("(ii) Work done during the process (W).")
+
+Q=330; //kJ/kg
+KE=(C2^2-C1^2)/2/10^3; //kJ
+PE=g*h/10^3; //kJ
+W=-Q-KE-PE-dh;
+disp("Work done = ")
+disp(W)
+disp("kJ")
+
+
+disp("Work done per second = ")
+P=W*10/60;
+disp(P)
+disp("kW")
\ No newline at end of file diff --git a/551/CH4/EX4.36/36.sce b/551/CH4/EX4.36/36.sce new file mode 100755 index 000000000..74c5a1352 --- /dev/null +++ b/551/CH4/EX4.36/36.sce @@ -0,0 +1,22 @@ +clc
+m=15; //kg/s
+v=0.45; //m^3/kg
+P=12000; //kW
+W=P/m; //kJ/kg
+h1=1260; //kJ/kg
+h2=400; //kJ/kg
+C1=50; //m/s
+C2=110; //m/s
+
+disp("(i) Heat rejected = ")
+Q=h2-h1+(C2^2-C1^2)/2/10^3 +W;
+Qnet=m*Q;
+disp("Qnet=")
+disp(-Qnet)
+disp("kW")
+
+disp("(ii) Inlet area")
+A=v*m/C1;
+disp("A=")
+disp(A)
+disp("m^2")
\ No newline at end of file diff --git a/551/CH4/EX4.37/37.sce b/551/CH4/EX4.37/37.sce new file mode 100755 index 000000000..79e8ef68b --- /dev/null +++ b/551/CH4/EX4.37/37.sce @@ -0,0 +1,30 @@ +clc
+m=0.5; //kg/s
+C1=6; //m/s
+C2=5; //m/s
+p1=1; //bar
+p2=7; //bar
+v1=0.85; //m^3/kg
+v2=0.16; //m^3/kg
+du=90; //kJ/kg
+Q=-120; //kJ/kg
+
+
+disp("(i) Power required to drive the compressor")
+W=-du+(C1^2-C2^2)/2/1000 + (p1*v1 - p2*v2)*10^2 + Q;
+Power=m*W;
+disp("Power=")
+disp(-Power)
+disp("kW")
+
+
+disp("(ii) Inlet and outlet pipe cross-sectional areas")
+A1=m*v1/C1;
+A2=m*v2/C2;
+disp("Inlet crosssectional area = ")
+disp(A1)
+disp("m^2")
+
+disp("Outlet crossectional area=")
+disp(A2)
+disp("m^2")
\ No newline at end of file diff --git a/551/CH4/EX4.38/38.sce b/551/CH4/EX4.38/38.sce new file mode 100755 index 000000000..0d979f9e2 --- /dev/null +++ b/551/CH4/EX4.38/38.sce @@ -0,0 +1,12 @@ +clc
+h1=800; //kJ/kg
+C1=5; //m/s
+h2=2520; //kJ/kg
+C2=50; //m/s
+dZ=4; //m
+g=9.81; //m/s^2
+Q=2180; //kJ/kg
+W=h1-h2+(C1^2 - C2^2)/2/1000 +dZ*g/1000+Q;
+disp("Power developed = ")
+disp(W)
+disp("kW")
\ No newline at end of file diff --git a/551/CH4/EX4.39/39.sce b/551/CH4/EX4.39/39.sce new file mode 100755 index 000000000..e49581e80 --- /dev/null +++ b/551/CH4/EX4.39/39.sce @@ -0,0 +1,14 @@ +clc
+g=9.8; //m/s^2
+m=4500/3600; //kg/s
+C1=2800/60; //m/s
+Z1=5.5; //m
+h1=2800; //kJ/g
+C2=5600/60; //m/s
+Z2=1.5; //m
+h2=2300; //kJ/kg
+Q=-16000/3600; //kJ/s
+W=Q-m*[(h1-h2) + (C2^2 - C1^2)/2/1000 + (Z2-Z1)*g/1000];
+disp("Power output of the turbine = ")
+disp(-W)
+disp("kW")
\ No newline at end of file diff --git a/551/CH4/EX4.4/4.sce b/551/CH4/EX4.4/4.sce new file mode 100755 index 000000000..21efac61b --- /dev/null +++ b/551/CH4/EX4.4/4.sce @@ -0,0 +1,11 @@ +clc
+p1=105; //kPa
+V1=0.4; //m^3
+p2=p1;
+V2=0.20; //m^3
+Q=-42.5; //kJ
+W=p1*(V2-V1);
+dU=Q-W;
+disp("change in internal energy = ")
+disp(dU)
+disp("kJ")
\ No newline at end of file diff --git a/551/CH4/EX4.40/40.sce b/551/CH4/EX4.40/40.sce new file mode 100755 index 000000000..4690e366e --- /dev/null +++ b/551/CH4/EX4.40/40.sce @@ -0,0 +1,11 @@ +clc
+p1=6.87; //bar
+C1=50; //m/s
+p2=1.37; //bar
+C2=500; //m/s
+disp("From steam table corresponding to p1")
+h1=2850; //kJ/kg
+h2=h1 - (C2^2-C1^2)/2/1000;
+disp("Final enthalpy of steam = ")
+disp(h2)
+disp("kJ")
\ No newline at end of file diff --git a/551/CH4/EX4.41/41.sce b/551/CH4/EX4.41/41.sce new file mode 100755 index 000000000..48163c0bb --- /dev/null +++ b/551/CH4/EX4.41/41.sce @@ -0,0 +1,15 @@ +clc
+m=220/60; //kg/s
+C1=320; //m/s
+p1=6*10^5; //N/m^2
+u1=2000*10^3; //J/kg
+v1=0.36; //m^3/kg
+C2=140; //m/s
+p2=1.2*10^5; //N/m^2
+u2=1400*10^3; //J/kg
+v2=1.3; //m^3/kg
+Q=100*10^3; //J/s
+W=(m*[(u1-u2)+ (p1*v1 - p2*v2) + (C1^2-C2^2)/2] -Q)/10^6;
+disp("power capacity of the system = ")
+disp(W)
+disp("MW")
\ No newline at end of file diff --git a/551/CH4/EX4.42/42.sce b/551/CH4/EX4.42/42.sce new file mode 100755 index 000000000..5feee9359 --- /dev/null +++ b/551/CH4/EX4.42/42.sce @@ -0,0 +1,13 @@ +clc
+p1=7.5*10^5; //N/m^2
+C1=140; //m/s
+h1=950*10^3; //J/kg
+p2=2*10^5; //N/m^2
+C2=280; //m/s
+h2=650*10^3; //J/kg
+m=5; //kg/s
+W=(h1-h2)+(C1^2-C2^2)/2
+Power=m*W/1000;
+disp("Power capacity of turbine = ")
+disp(Power)
+disp("kW")
\ No newline at end of file diff --git a/551/CH4/EX4.43/43.sce b/551/CH4/EX4.43/43.sce new file mode 100755 index 000000000..acba2a3e9 --- /dev/null +++ b/551/CH4/EX4.43/43.sce @@ -0,0 +1,22 @@ +clc
+C1=12; //m/s
+p1=1*10^5; //N/m^2
+v1=0.5; //m^3/kg
+C2=90; //m/s
+p2=8*10^5; //N/m^2
+v2=0.14; //m^3/kg
+dh=150; //kJ/kg
+Q=-11.67; //kJ/s
+m=0.2; //kg/s
+
+disp("(i) Motor power required to drive the compressor")
+W=m*[-dh + (C1^2-C2^2)/2/1000] +Q;
+disp("Power=")
+disp(-W)
+disp("kW")
+
+
+disp("(ii)Ratio of inlet to outlet pipi diameter")
+ratio=sqrt(C2/C1*v1/v2);
+disp("ratio=")
+disp(ratio)
\ No newline at end of file diff --git a/551/CH4/EX4.44/44.sce b/551/CH4/EX4.44/44.sce new file mode 100755 index 000000000..ea1304d75 --- /dev/null +++ b/551/CH4/EX4.44/44.sce @@ -0,0 +1,9 @@ +clc
+W=-175; //kJ/kg
+dh=70; //kJ/kg
+Q_water=-92; //kJ/kg
+Q=dh+W;
+Q_atm=Q-Q_water;
+disp("Heat transferred to the atmosphere = ")
+disp(-Q_atm)
+disp("kJ/kg")
\ No newline at end of file diff --git a/551/CH4/EX4.45/45.sce b/551/CH4/EX4.45/45.sce new file mode 100755 index 000000000..24c2bd6e7 --- /dev/null +++ b/551/CH4/EX4.45/45.sce @@ -0,0 +1,27 @@ +clc
+h1=2800*10^3; //J/kg
+C1=50; //m/s
+A1=900*10^(-4); //m^2
+v1=0.187; //m^3/kg
+h2=2600*10^3; //J/kg
+v2=0.498; //m^3/kJ
+disp("(i) Velocity at exit of the nozzle")
+C2=sqrt(2*[(h1-h2) + C1^2/2]);
+
+disp("C2=")
+disp(C2)
+disp("m/s")
+
+
+disp("(ii) Mass flow rate")
+m=A1*C1/v1;
+disp("m=")
+disp(m)
+disp("kg/s")
+
+
+disp("(iii) Area at the exit")
+A2=m*v2/C2*10^4;
+disp("A2=")
+disp(A2)
+disp("cm^2")
\ No newline at end of file diff --git a/551/CH4/EX4.46/46.sce b/551/CH4/EX4.46/46.sce new file mode 100755 index 000000000..ed10ee594 --- /dev/null +++ b/551/CH4/EX4.46/46.sce @@ -0,0 +1,10 @@ +clc
+h1=240; //kJ/kg
+h2=192; //kJ/kg
+dZ=20; //m
+g=9.81; //m/s^2
+Q=(h2-h1)+dZ*g/1000;
+
+disp("heat transfer = ")
+disp(-Q)
+disp("kJ/kg")
\ No newline at end of file diff --git a/551/CH4/EX4.47/47.sce b/551/CH4/EX4.47/47.sce new file mode 100755 index 000000000..6e419a580 --- /dev/null +++ b/551/CH4/EX4.47/47.sce @@ -0,0 +1,10 @@ +clc
+p1=2; //bar
+C1=300; //m/s
+Q=0;
+h1=915*10^3; //J/kg
+h2=800*10^3; //J/kg
+C2=sqrt(2*[h1-h2 + C1^2/2]);
+disp("Relative velocity of gas leaving the pipe=")
+disp(C2)
+disp("m/s")
\ No newline at end of file diff --git a/551/CH4/EX4.48/48.sce b/551/CH4/EX4.48/48.sce new file mode 100755 index 000000000..515ae81b2 --- /dev/null +++ b/551/CH4/EX4.48/48.sce @@ -0,0 +1,17 @@ +clc
+mw=50; //kg/s
+p1=10^5; //N/m^2
+p2=4.2*10^5; //N/m^2
+h=10.7; //m
+d1=0.2; //m
+d2=0.1; //m
+v1=1/1000;
+v2=1/1000;
+g=9.81; //m/s^2
+
+C1=mw*4/%pi/d1^2*v1;
+C2=mw*4/%pi/d2^2*v2;
+W=mw*[(p1*v1-p2*v2) + (g*(0-h))+(C1^2-C2^2)/2]/10^3;
+disp("Capacity of electric motor")
+disp(-W)
+disp("kW")
\ No newline at end of file diff --git a/551/CH4/EX4.49/49.sce b/551/CH4/EX4.49/49.sce new file mode 100755 index 000000000..240864beb --- /dev/null +++ b/551/CH4/EX4.49/49.sce @@ -0,0 +1,16 @@ +clc
+Ca=250; //m/s
+t=-14; //0C
+ha=250; //kJ/kg
+hg=900; //kJ/kg
+ratio=0.0180;
+Ef=45*10^3; //kJ/kg
+Q=-21; //kJ/kg
+ma=1; //kg
+mg=1.018; //kg
+mf=0.018; //kg
+Eg=0.06*mf/mg*Ef;
+Cg=sqrt(2000*([ma*(ha+Ca^2/2/1000) + mf*Ef + Q]/mg -hg-Eg));
+disp("velocity of exhaust gas jet = ")
+disp(Cg)
+disp("m/s")
\ No newline at end of file diff --git a/551/CH4/EX4.50/50.sce b/551/CH4/EX4.50/50.sce new file mode 100755 index 000000000..a8f701e3a --- /dev/null +++ b/551/CH4/EX4.50/50.sce @@ -0,0 +1,31 @@ +clc
+t1=20; //0C
+C1=40; //m/s
+t2=820; //0C
+C2=40; //m/s
+t3=620; //0C
+C3=55; //m/s
+t4=510; //0C
+m=2.5; //kg/s
+cp=1.005; //kJ/kg.0C
+
+
+disp("(i) Heat exchanger")
+Q_12=m*cp*(t2-t1);
+disp("rate of heat transfer=")
+disp(Q_12)
+disp("kJ/s")
+
+
+disp("(ii) Turbine")
+W_23=m*[(cp*(t2-t3))+(C2^2-C3^2)/2/1000];
+disp("Power output of turbine=")
+disp(W_23)
+disp("kW")
+
+
+disp("(iii) Nozzle")
+C4=sqrt(2*1000*(cp*(t3-t4)+C3^2/2/1000));
+disp("Velocity at exit from the nozzle= ")
+disp(C4)
+disp("m/s")
\ No newline at end of file diff --git a/551/CH4/EX4.51/51.sce b/551/CH4/EX4.51/51.sce new file mode 100755 index 000000000..816c8187b --- /dev/null +++ b/551/CH4/EX4.51/51.sce @@ -0,0 +1,28 @@ +clc
+V=0.028; //m^3
+p1=80; //bar
+t=350; //0C
+p2=50; //bar
+v1=0.02995; //m^3/kg
+h1=2987.3; //kJ/kg
+v2=0.02995; //m^3/kg
+vg2=0.0394; //m^3/kg
+uf2=1149; //kJ/kg
+ug2=2597; //kJ/kg
+
+m=V/v1;
+u1=h1 - (p1*v1*10^2); //kJ/kg
+
+
+disp("(i) State of steam after cooling")
+x2=v2/vg2;
+disp("dryness fraction = ")
+disp(x2)
+
+
+disp("(ii) Heat rejected by the steam")
+u2=(1-x2)*uf2 + x2*ug2;
+Q=m*(u2-u1);
+disp("Heat rejected = ")
+disp(-Q)
+disp("kJ")
\ No newline at end of file diff --git a/551/CH4/EX4.52/52.sce b/551/CH4/EX4.52/52.sce new file mode 100755 index 000000000..18da48a90 --- /dev/null +++ b/551/CH4/EX4.52/52.sce @@ -0,0 +1,22 @@ +clc
+m=0.08; //kg
+p=2*10^5; //Pa
+V=0.10528; //m^3
+h1=2706.3; //kJ/kg
+h2=3071.8; //kJ/kg
+v1=0.885; //m^3/kg
+v2=V/m; //m^3/kg
+
+
+disp("(i) Heat supplied")
+Q=m*(h2-h1);
+disp("Q=")
+disp(Q)
+disp("kJ")
+
+disp("(ii) Work done")
+W=p*(v2-v1);
+W_total=m*W/10^3;
+disp("Total work done = ")
+disp(W_total)
+disp("kJ")
\ No newline at end of file diff --git a/551/CH4/EX4.53/53.sce b/551/CH4/EX4.53/53.sce new file mode 100755 index 000000000..ab3be111f --- /dev/null +++ b/551/CH4/EX4.53/53.sce @@ -0,0 +1,46 @@ +clc
+m=1; //kg
+p=8; //bar
+s1=6.55; //kJ/kg.K
+T=200; //0C
+s_f1=2.0457; //kJ/kg.K
+s_fg1=4.6139; //kJ/kg.K
+h_f1=720.9; //kJ/kg
+h_fg1=2046.5; //kJ/kg
+h2=2839.3; //kJ/kg
+
+x1=(s1-s_f1)/s_fg1;
+h1=h_f1+x1*h_fg1;
+Q=h2-h1;
+disp("Heat supplied=")
+disp(Q)
+disp("kJ/kg")
+
+// For T-s diagram
+
+s=0:0.01:10;
+T=(-(s-5)^2+298);
+plot(s,T)
+
+T=[295.44 295.44];
+s=[6.6 3.45];
+plot(s,T,'g')
+
+s=[6.6 7];
+T=[295.44 300];
+plot(s,T,'g')
+
+s=[6.55 6.55];
+T=[270 295.44];
+plot(s,T,'r')
+
+s=[6.6 6.6];
+T=[270 295.44];
+plot(s,T,'--r')
+
+s=[6.66 6.66];
+T=[270 295.44];
+plot(s,T,'r')
+
+xtitle("T-s diagram", "s(kJ/kg K)", "T(K)")
+//The area in red represents the heat flow and it goes upto x-axis
\ No newline at end of file diff --git a/551/CH4/EX4.54/54.sce b/551/CH4/EX4.54/54.sce new file mode 100755 index 000000000..d22eac08a --- /dev/null +++ b/551/CH4/EX4.54/54.sce @@ -0,0 +1,39 @@ +clc
+p1=7*10^5; //Pa
+p2=1.5*10^5; //Pa
+Q=420; //kJ/kg
+uf=696; //kJ/kg
+x=0.95;
+ug=2573; //kJ/kg
+u_f2=2580; //kJ/kg
+u_g2=2856; //kJ/kg
+x2=15/50;
+h_f1=697.1; //kJ/kg
+h_fg1=2064.9; //kJ.kg
+h_f2=2772.6; //kJ/kg
+h_g2=2872.9; //kJ/kg
+
+
+disp("(i) Change of internal energy")
+u1=(1-x)*uf + x*ug;
+u2=2602.8; //kJ/kg
+du=u2-u1;
+disp("du=")
+disp(du)
+disp("kJ/kg")
+
+
+disp("(ii) Change in enthalpy")
+h1=h_f1+x*h_fg1;
+h2=h_f2+x2*(h_g2-h_f2);
+dh=h2-h1;
+disp("dh=")
+disp(dh)
+disp("kJ/kg")
+
+
+disp("(iii) Work done = ")
+W=Q-du;
+disp("W=")
+disp(W)
+disp("kJ/kg")
\ No newline at end of file diff --git a/551/CH4/EX4.55/55.sce b/551/CH4/EX4.55/55.sce new file mode 100755 index 000000000..a6e182780 --- /dev/null +++ b/551/CH4/EX4.55/55.sce @@ -0,0 +1,21 @@ +clc
+p1=5.5*10^5; //Pa
+x1=1;
+p2=0.75*10^5; //Pa
+v1=0.3427; //m^3/kg
+v2=p1*v1/p2;
+// Since v2 > vg (at 0.75 bar), therefore, the steam is superheated at state 2.
+u2=2567.25; //kJ/kg
+u1=2565; //kJ/kg
+du=u2-u1; //kJ/kg
+C=p1*v1;
+
+disp("Work done = ")
+W=integrate('C/v', 'v', v1,v2)
+disp("N-m/kg")
+
+
+disp("Heat supplied = ")
+Q=du+W/10^3;
+disp(Q)
+disp("kJ/kg")
\ No newline at end of file diff --git a/551/CH4/EX4.56/56.sce b/551/CH4/EX4.56/56.sce new file mode 100755 index 000000000..0e257fd03 --- /dev/null +++ b/551/CH4/EX4.56/56.sce @@ -0,0 +1,21 @@ +p1=100; //bar
+p2=10; //bar
+s1=5.619; //kJ/kg.K
+T=584; //K
+s2=7.163; //kJ/kg.K
+u1=2545; //kJ/kg
+u2=2811.8; //kJ/kg
+
+
+disp("(i)Heat supplied ")
+Q=T*(s2-s1);
+disp("Q=")
+disp(Q)
+disp("kJ/kg")
+
+
+disp("(ii) Work done")
+W=Q-(u2-u1);
+disp("W=")
+disp(W)
+disp("kJ/kg")
\ No newline at end of file diff --git a/551/CH4/EX4.57/57.sce b/551/CH4/EX4.57/57.sce new file mode 100755 index 000000000..2b673337b --- /dev/null +++ b/551/CH4/EX4.57/57.sce @@ -0,0 +1,13 @@ +clc
+m=1; //kg
+p1=120*10^5; //N/m^2
+t1=400; //0C
+p2=38; //bar
+h1=3051.3; //kJ/kg
+v1=0.02108; //m^3/kg
+u1=h1-p1*v1/10^3; //kJ/kg
+u2=2602; //kJ/kg
+disp("WOrk done = ")
+W=u1-u2;
+disp(W)
+disp("kJ/kg")
\ No newline at end of file diff --git a/551/CH4/EX4.58/58.sce b/551/CH4/EX4.58/58.sce new file mode 100755 index 000000000..15bb06c02 --- /dev/null +++ b/551/CH4/EX4.58/58.sce @@ -0,0 +1,32 @@ +clc
+p1=7*10^5; //N/m^2
+x1=0.98;
+p2=0.34*10^5; //N/m^2
+vg=0.273; //m^3/kg
+n=1.1;
+v_g2=4.65; //m^3/kg
+u_f1=696; //kJ/kg
+u_g1=2573; //kJ/kg
+u_f2=302; //kJ/kg
+u_g2=2472; //kJ/kg
+
+
+v1=x1*vg;
+v2=v1*(p1/p2)^(1/n);
+x2=v2/v_g2;
+
+
+disp("(i) Work done by the steam during the process")
+W=(p1*v1-p2*v2)/(n-1)/10^3; //kJ/kg
+disp("W=")
+disp(W)
+disp("kJ/kg")
+
+
+disp("(ii) Heat transferred")
+u1=(1-x1)*u_f1+x1*u_g1;
+u2=(1-x2)*u_f2+x2*u_g2;
+Q=u2-u1 + W;
+disp("Q=")
+disp(Q)
+disp("kJ/kg")
\ No newline at end of file diff --git a/551/CH4/EX4.59/59.sce b/551/CH4/EX4.59/59.sce new file mode 100755 index 000000000..2d3adc890 --- /dev/null +++ b/551/CH4/EX4.59/59.sce @@ -0,0 +1,21 @@ +clc
+p1=15; //bar
+t1=350; //0C
+C1=60; //m/s
+p2=1.2; //bar
+C2=180; //m/s
+s1=7.102; //kJ/kg
+s_f2=1.3609; //kJ/kg
+s_g2=7.2884; //kJ/kg
+h_f2=439.4; //kJ/kg
+h_fg2=2241.1; //kJ/kg
+h1=3147.5; //kJ/kg
+
+
+x2=(s1 - s_f2)/(s_g2-s_f2);
+h2=h_f2+x2*h_fg2;
+
+W=(h1-h2) + (C1^2 - C2^2)/2/1000;
+disp("Work done = ")
+disp(W)
+disp("kJ/kg")
\ No newline at end of file diff --git a/551/CH4/EX4.6/25.sce b/551/CH4/EX4.6/25.sce new file mode 100755 index 000000000..3de31da00 --- /dev/null +++ b/551/CH4/EX4.6/25.sce @@ -0,0 +1,15 @@ +clc
+p=1.1; //bar
+x=0.95;
+c_pw=4.18;
+m1=90; //kg
+m2=5.25; //kg
+T1=25; //0C
+T2=40; //0C
+m=m1+m2;
+h_f=428.8; //kJ/kg
+h_fg= 2250.8; //kJ/kg
+m_s= [m*c_pw*(T2-T1)]/[(h_f + x*h_fg) - c_pw*T2];
+disp("Mass of steam condensed=")
+disp(m_s)
+disp("kg")
\ No newline at end of file diff --git a/551/CH4/EX4.60/60.sce b/551/CH4/EX4.60/60.sce new file mode 100755 index 000000000..430146318 --- /dev/null +++ b/551/CH4/EX4.60/60.sce @@ -0,0 +1,13 @@ +clc
+p1=10; //bar
+t1=200; //0C
+C1=60; //m/s^2
+c2=650; //m/s
+p2=1.5; //bar
+h1=2827.9; //kJ/kg
+h_f2=467.1; //kJ/kg
+h2=2618.45; //kJ/kg
+h_g2=2693.4; //kJ/kg
+x2=(h2-h_f2)/(h_g2-h_f2);
+disp("quality of steam leaving the nozzle=")
+disp(x2)
\ No newline at end of file diff --git a/551/CH4/EX4.61/61.sce b/551/CH4/EX4.61/61.sce new file mode 100755 index 000000000..e201c6b9b --- /dev/null +++ b/551/CH4/EX4.61/61.sce @@ -0,0 +1,8 @@ +clc
+h1=2776.4; //kJ/kg
+h2=h1;
+h_f1=884.6; //kJ/kg
+h_fg1=1910.3; //kJ/kg
+x1=(h1-h_f1)/h_fg1;
+disp("Initial dryness fraction = ")
+disp(x1)
\ No newline at end of file diff --git a/551/CH4/EX4.62/62.sce b/551/CH4/EX4.62/62.sce new file mode 100755 index 000000000..0b979f644 --- /dev/null +++ b/551/CH4/EX4.62/62.sce @@ -0,0 +1,9 @@ +p1=10; //bar
+x1=0.9; //bar
+p2=2; //bar
+
+
+// Using Mollier chart, we get
+x2=0.94;
+disp("x2 =")
+disp(x2)
\ No newline at end of file diff --git a/551/CH4/EX4.63/63.sce b/551/CH4/EX4.63/63.sce new file mode 100755 index 000000000..c92ae57b4 --- /dev/null +++ b/551/CH4/EX4.63/63.sce @@ -0,0 +1,74 @@ +clc
+
+disp("(a)From steam tables")
+
+p1=15*10^5; //Pa
+p2=7.5*10^5; //Pa
+h_f1=844.7; //kJ/kg
+ts1=198.3; //0C
+s_f1=2.3145; //kJ/kg.K
+s_g1=6.4406; //kJ/kg.K
+v_g1=0.132; //m^3/kg
+h_fg1=1945.2; //kJ/kg
+x1=0.95;
+h_f2=709.3; //kJ/kg
+h_fg2=2055.55; //kJ/kg
+s_f2=2.0195; //kJ/kg
+s_g2=6.6816; //kJ/kg.K
+v_g2=0.255; //m^3/kg
+x2=0.9;
+x3=1;
+s_f3=0.521; //kJ/kg K
+s_g3=8.330; //kJ/kg K
+
+h2=h_f2+x2*h_fg2;
+h1=h_f1 + x1*h_fg1;
+s1=s_f1 + x1*(s_g1-s_f1);
+s2=s1;
+ds_12=s2-s1;
+
+s3=s_f3+x3*(s_g3-s_f3);
+ds_23=s3-s2;
+
+ds=ds_12+ds_23;
+disp("(i) Change in entropy =")
+disp(ds)
+disp("kJ/kg K")
+
+h3=h2;
+
+disp("(ii) Change in enthalpy")
+dh=h2-h1;
+disp(dh)
+disp("kJ/kg")
+
+
+disp("(iii) Change in internal energy")
+u1=h1-p1*x1*v_g1/10^3;
+u2=h2-p2*x2*v_g2/10^3;
+du=u2-u1;
+disp("du=")
+disp(du)
+disp("kJ/kg")
+
+
+// Only the expansion of steam from point 1 to 2 (i.e., isentropic expansion) is reversible because of unresisted flow whereas the expansion from point 2 to point 3 (i.e., throttling expansion) is irreversible because of frictional resistance to flow. Increase of entropy also shows that expansion from point 2 to point 3 is irreversible.
+
+
+disp("(b) Using Mollier chart")
+h1=2692; //kJ/kg
+h2=2560; //kJ/kg
+s1=6.23; //kJ/kg K
+s2=s1;
+s3=8.3; //kJ/kg K
+
+disp("(i) Change in entropy =")
+ds=s3-s1;
+disp(ds)
+disp("kJ/kg K")
+
+
+disp("(ii) Change in enthalpy =")
+dh=h2-h1;
+disp(dh)
+disp("kJ/kg")
\ No newline at end of file diff --git a/551/CH4/EX4.64/64.sce b/551/CH4/EX4.64/64.sce new file mode 100755 index 000000000..8a566ef3a --- /dev/null +++ b/551/CH4/EX4.64/64.sce @@ -0,0 +1,17 @@ +clc
+V1=5.5; //m^3
+p1=16*10^5; //Pa
+T1=315; //K
+V2=V1;
+p2=12*10^5; //Pa
+R=0.287*10^3;
+y=1.4;
+
+m1=p1*V1/R/T1;
+T2=T1*(p2/p1)^((y-1)/y);
+
+m2=p2*V2/R/T2;
+disp("Mass of air which left the receiver=")
+m=m1-m2;
+disp(m)
+disp("kg")
\ No newline at end of file diff --git a/551/CH4/EX4.65/65.sce b/551/CH4/EX4.65/65.sce new file mode 100755 index 000000000..df54e6221 --- /dev/null +++ b/551/CH4/EX4.65/65.sce @@ -0,0 +1,20 @@ +clc
+cp=1; //kJ/kg.K
+cv=0.711; //kJ/kg.K
+V1=1.6; //m^3
+V2=V1;
+p1=5*10^5; //Pa
+T1=373; //K
+p2=1*10^5; //Pa
+R=287;
+y=1.4;
+
+m1=p1*V1/R/T1;
+T2=T1*(p2/p1)^((y-1)/y);
+m2=p2*V2/R/T2;
+KE=(m1*cv*T1)-(m2*cv*T2)-(m1-m2)*cp*T2;
+disp("Kinetic energy of discharge air =")
+disp(KE)
+disp("kJ")
+
+disp("This is the exact answer when using proper value of cv")
\ No newline at end of file diff --git a/551/CH4/EX4.66/66.sce b/551/CH4/EX4.66/66.sce new file mode 100755 index 000000000..12a8c7662 --- /dev/null +++ b/551/CH4/EX4.66/66.sce @@ -0,0 +1,52 @@ +clc
+//For oxygen
+cpa=0.88; //kJ/kg K
+Ra=0.24; //kJ/kg K
+V1a=0.035; //m^3
+p1a=4.5; //bar
+T1a=333; //K
+V2a=0.07; //m^3
+
+//For methane
+V1b=0.07; //m^3
+V2b=0.035; //m^3
+p1b=4.5; //bar
+T1b=261; //K
+cpb=1.92; //kJ/kg K
+Rb=0.496; //kJ/kg K
+
+yb=cpb/(cpb-Rb); //for methane
+cva=cpa-Ra; //for oxygen
+
+disp("(i) Final state condition")
+
+p2b=p1b*(V1b/V2b)^yb;
+disp("p2 for methane =")
+disp(p2b)
+disp("bar")
+
+T2b=p2b*V2b*T1b/p1b/V1b;
+disp("T2 for methane =")
+disp(T2b)
+disp("K")
+
+p2a=p2b;
+
+T2a=p2a*V2a/p1a/V1a*T1a;
+disp("T2 for oxygen =")
+disp(T2a)
+disp("K")
+
+Wb=(p1b*V1b - p2b*V2b)/(yb-1)*100; //kJ
+
+disp("(ii) Work done by the piston ")
+disp("The piston will be in virtual equilibrium and hence zero work is effected by the piston.")
+
+Wa=-Wb;
+
+ma=p1a*V1a/Ra/T1a*10^2;
+
+Q=ma*cva*(T2a-T1a) + Wa;
+disp("(iii) Heat transferred to oxygen =")
+disp(Q)
+disp("kJ")
\ No newline at end of file diff --git a/551/CH4/EX4.7/7.sce b/551/CH4/EX4.7/7.sce new file mode 100755 index 000000000..0a9620e1e --- /dev/null +++ b/551/CH4/EX4.7/7.sce @@ -0,0 +1,10 @@ +clc
+W_12=-82; //kJ
+Q_12=-45; //kJ
+dU_12=Q_12 - W_12;
+W_21=100; //kJ
+dU_21=-dU_12;
+Q_21=dU_21 + W_21;
+disp("Heat added to the system = ")
+disp(Q_21)
+disp("kJ")
\ No newline at end of file diff --git a/551/CH4/EX4.8/8.sce b/551/CH4/EX4.8/8.sce new file mode 100755 index 000000000..430763545 --- /dev/null +++ b/551/CH4/EX4.8/8.sce @@ -0,0 +1,12 @@ +clc
+Q2=9000; //kJ
+Q1=3000; //kJ
+Q=Q1-Q2;
+W=0;
+dU=W-Q;
+disp("Work done = ")
+disp(W)
+
+disp("Change in internal energy = ")
+disp(dU)
+disp("kJ")
\ No newline at end of file diff --git a/551/CH4/EX4.9/9.sce b/551/CH4/EX4.9/9.sce new file mode 100755 index 000000000..640186906 --- /dev/null +++ b/551/CH4/EX4.9/9.sce @@ -0,0 +1,23 @@ +m=20; //kg
+g=9.81; //m/s^2
+z2=0;
+z1=15;
+
+disp("(i) When the stone is about to enter the water")
+Q=0
+W=0
+dU=0
+PE=m*g*(z2-z1)
+KE=-PE
+
+disp("(ii) When the stone dips into the tank and comes to rest")
+Q=0
+W=0
+KE=0
+PE=m*g*(z2-z1)
+dU=-PE
+
+disp("(iii) When the water and stone come to their initial temperature")
+W=0
+KE=0
+Q=-dU
\ No newline at end of file |