initial commit / add all books

1 files changed, 32 insertions, 0 deletions

diff --git a/611/CH3/EX3.14/Chap3_Ex14_R1.sce b/611/CH3/EX3.14/Chap3_Ex14_R1.sce
new file mode 100755
index 000000000..4024868e3
--- /dev/null
+++ b/611/CH3/EX3.14/Chap3_Ex14_R1.sce
@@ -0,0 +1,32 @@
+// Y.V.C.Rao ,1997.Chemical Engineering Thermodynamics.Universities Press,Hyderabad,India.

+

+//Chapter-3,Example 14,Page 75

+//Title:Volume using generalized form of the Redlich-Kwong equation of state

+//================================================================================================================

+clear 

+clc

+

+//INPUT

+T=427.85;//temperature of n-octane vapour in K

+P=0.215;//pressure of n-ocatne vapour in MPa

+Tc=569.4;//critical temperature of n-octane in K

+Pc=2.497;//critical pressure of n-octane in MPa

+R=8.314;//universal gas constant in (Pa m^3)/(mol K)

+

+//CALCULATION

+Tr=T/Tc;//calculation of reduced temperature (no unit)

+Pr=P/Pc;//calculation of reduced pressure (no unit)

+z_init=1;//taking a guess value of z (compressibilty factor) to get a value of h for solving the system

+h=(0.08664*Pr)/(z_init*Tr);//calculation of h using Eq.(3.68)

+tolerance=1e-6;//Framing the tolerance limit for the convergence of the equation

+function[fn]=solver_func(zi)

+    fn=zi-((1/(1-h))-((h/(1+h))*(4.93398/(Tr^(3/2)))));//Function defined for solving the system using Eq.(3.67)

+endfunction

+[z]=fsolve(h,solver_func,tolerance)//using inbuilt function fsolve for solving the system of equations

+V=(z*R*T)/(P*10^6);//calculation of volume in m^3/mol using Eq.(3.59)

+

+//OUTPUT

+mprintf('\n The volume occupied by n-octane vapour obtained by the generalized form of Redlich-Kwong equation of state= %f m^3/mol\n',V);

+

+

+//===============================================END OF PROGRAM===================================================