diff options
Diffstat (limited to 'Thermodynamics:_A_Core_Course/CH6.ipynb')
| -rw-r--r-- | Thermodynamics:_A_Core_Course/CH6.ipynb | 82 |
1 files changed, 0 insertions, 82 deletions
diff --git a/Thermodynamics:_A_Core_Course/CH6.ipynb b/Thermodynamics:_A_Core_Course/CH6.ipynb index b2d73f0e..4326009b 100644 --- a/Thermodynamics:_A_Core_Course/CH6.ipynb +++ b/Thermodynamics:_A_Core_Course/CH6.ipynb @@ -27,25 +27,16 @@ "cell_type": "code",
"collapsed": false,
"input": [
- "#Find molar volume of CO2\n",
"\n",
"from scipy.optimize import fsolve\n",
"\n",
- "#Variable declaration:\n",
- "#For CO2:\n",
"a=3.61 #atm L**2 mol**-2\n",
"b=4.29*10**-2 #L mol**-1\n",
"R=0.082 #L atm K**-1 mol**-1\n",
"T=500 #K\n",
"P=100 #atm\n",
"\n",
- "#CALCULATION\n",
"\n",
- "#(P+a/Vm^2)(Vm-b)=RT\n",
- "#or Vm^2(P+a/Vm^2)(Vm-b)=RTVm^2\n",
- "#or P*Vm**3+a*Vm-Pb*Vm**2-a*b=Rt*Vm**2\n",
- "#or Vm^^3-(b+RT/P)Vm**2+a/P*Vm-(a*b)/P=0\n",
- "#Let\n",
"C1=b+(R*T/P) #L mol**-1 [aSsume]\n",
"C2=a/P #L^2 mol^-2 [assume]\n",
"C3=C2*b #L^3mol**-3\n",
@@ -53,7 +44,6 @@ " return(x**3-C1*x**2+C2*x-C3)\n",
"x=fsolve(f,0.3)\n",
"\n",
- "#RESULT\n",
"\n",
"print \"x=\",round(x,3)\n",
"Vm=round(x,3)\n",
@@ -85,23 +75,18 @@ "cell_type": "code",
"collapsed": false,
"input": [
- "#example 6.5\n",
- "#To find the fugacity and fugacity coefficient\n",
- "#Variable declaration\n",
"b=0.0391 \t\t\t#Van der waals constant[dm3/mol]\n",
"R=0.082 \t\t\t#Universal gas constant[dm3*atm/mol]\n",
"P2=1000 \t\t\t#pressure [atm]\n",
"P1=0 \t\t\t\t#pressure [atm]\n",
"T=1273\t\t \t\t#Temperature [K]\n",
"import math\n",
- "#Calculation\n",
"\n",
"x=b*(P2-P1) \n",
"y=R*T \n",
"fc=math.exp(x/y) \t\t#fugacity coefficient\n",
"\n",
"f=P2*fc #fugacity[atm]\n",
- "#Result\n",
"print\"The fugacity coefficient is\",round(fc,3) \n",
"print\"The fugacity is\",round(f),\"atm\"\n"
],
@@ -131,22 +116,17 @@ "cell_type": "code",
"collapsed": false,
"input": [
- "#example 6.10\n",
- "#To find the partial pressure of CO2 gas \n",
- "#Variable declaration\n",
"m1=0.03 #mass of CO2(g)[gm]\n",
"w1=44.01 #molecular weight of CO2(g)[gm/mol]\n",
"m2=250 #mass of water[gm]\n",
"w2=18.02 #molecular weight of water[gm/mol]\n",
"k=1.25*10**6 #Henry's law constant[Torr]\n",
"T=298 #Temperature[K]\n",
- "#Calculation\n",
"\n",
"n1=m1/w1 #no. of moles of CO2\n",
"n2=m2/w2 #no. of moles of water\n",
"x1=n1/(n1+n2) #mole fraction of CO2\n",
"Pco2=k*x1 #Partial pressure of CO2[Torr]\n",
- "#Result\n",
"print\"The partial pressure of CO2 gas is\",round(Pco2,2),\"Torr\""
],
"language": "python",
@@ -174,10 +154,7 @@ "cell_type": "code",
"collapsed": false,
"input": [
- "#example 6.11\n",
"\n",
- "#To find the Volume of the solution\n",
- "#Variable declaration\n",
"W=1000 \t\t\t#Total mass of a solution[gm]\n",
"x1=0.5 \t\t\t#mole fraction of Chloroform\n",
"x2=0.5 \t\t\t#mole fraction of Acetone\n",
@@ -185,11 +162,9 @@ "V2m=74.166 \t\t#Partial molar volume of Acetone[cm3/mol]\n",
"M1=119.59 \t\t#molecular weight of chloroform[gm/mol]\n",
"M2=58 \t\t\t#molecular weight of Acetone[gm/mol]\n",
- "#Calculation\n",
"\n",
"nT=W/(x1*M1+x2*M2) \t#Total no. of moles\n",
"V=nT*(x1*V1m+x2*V2m) \t#Total volume[cm3]\n",
- "#Result\n",
"print\"The volume of the solution is\",round(V,1),\"cm^3 (approx)\""
],
"language": "python",
@@ -217,16 +192,11 @@ "cell_type": "code",
"collapsed": false,
"input": [
- "#example 6.12\n",
- "#to find the excess volume \n",
- "#Variable declaration\n",
"x1=0.5 #mole fraction of chloroform\n",
"x2=0.5 #mole fraction of p-xylene\n",
"T=298 #Temperature[K]\n",
- "#Calculation\n",
"\n",
"Ve=x1*x2*(0.585+0.085*(x1-x2)-0.165*(x1-x2)**2) #Excess volume measured by using a dilatometer\n",
- "#Result\n",
"print\"Ve/(cm3.mol**-1) = \",round(Ve,3) "
],
"language": "python",
@@ -254,10 +224,7 @@ "cell_type": "code",
"collapsed": false,
"input": [
- "#example 6.14\n",
"\n",
- "#To find the activity , molality of the electrolytes\n",
- "#Variable declaration\n",
"m1=0.01 \t\t#molality[m]\n",
"v11=1.0 \n",
"v12=2.0 \n",
@@ -267,7 +234,6 @@ "v22=1.0 \n",
"Y2=0.53 \n",
"\n",
- "#Calculation\n",
"\n",
"v1=(v11)+(v12) \n",
"v2=(v21)+(v22) \n",
@@ -279,12 +245,9 @@ "y=1.0/v2 \n",
"m2m=m2*(v21*v21*v22**v22)**y #molality[m]\n",
"a2m=a2**y \n",
- "#Result\n",
"print\"The activity of the electrolyte ZnCl2 is\",round(a1,8)\n",
"print\"The activity of the electrolyte CuSO4 is\",round(a2,8)\n",
- "#print\"The mean activity of ZnCl2 is\",a1m\n",
"print\"The mean molality of ZnCl2 in [m]\",round(m1m,4)\n",
- "#print\"The mean activity of CuSO4 is\",a2m \n",
"print\"The mean molality of CuSO4 in [m]\",m2m \n"
],
"language": "python",
@@ -315,20 +278,15 @@ "cell_type": "code",
"collapsed": false,
"input": [
- "#example 6.15\n",
"\n",
- "#To find the molecular weight of sucrose\n",
- "#Variable declaration\n",
"m2=3 \t\t\t#mass of the sucrose[gm]\n",
"m1=0.1 \t\t\t#mass of water [Kg]\n",
"Kf=1.86 \t\t#cryoscopic constant of water[K*Kg/mol]\n",
"dTf=0.16 \t\t#Lowering in freezing point[K]\n",
- "#Calculation\n",
"\t\n",
"a=m1*dTf \n",
"b=Kf*m2 \n",
"M2=b/a \t\t\t#molecular weight\n",
- "#Result\n",
"print\"M2=molecular weight , then M2=\",M2 "
],
"language": "python",
@@ -356,20 +314,15 @@ "cell_type": "code",
"collapsed": false,
"input": [
- "#example 5.16\n",
"\n",
- "#To find the molecular formula of sulphur\n",
- "#Variable declaration\n",
"dTf=0.088 \t\t\t#Lowering in freezing point[K]\n",
"m2=0.45 \t\t\t#mass of sulphur[gm]\n",
"m1=0.09955 \t\t\t#mass of benzene[gm]\n",
"Kf=5.07 \t\t\t#cryoscopic constant for benzene[K*Kg/mol]\n",
- "#Calculation\n",
"\n",
"a=m1*dTf \n",
"b=Kf*m2 \n",
"M2=b/a \t\t\t\t#molecular weight of sulphur\n",
- "#Result\n",
"print\"The molecular weight of sulphur is\",round(M2,1) \n",
"x=M2/32 \t\t\t#no. of sulphur atoms\n",
"print\"\\n The molecular formula of sulphur is S\",round(x) "
@@ -401,9 +354,6 @@ "cell_type": "code",
"collapsed": false,
"input": [
- "#example 6.17\n",
- "#To find the molar mass of macromolecule\n",
- "#Variable declaration\n",
"m2=1.35 \t\t\t#mass of a macromolecule[gm]\n",
"V=100\t \t\t\t#volume of solution[cm^3]\n",
"R=82 \t\t\t\t#Universal gas constant[atm.cm^3.K^-1]\n",
@@ -412,13 +362,11 @@ "d=1 \t\t\t\t#density\n",
"p=1013250 \t\t\t#Atmospheric pressure\n",
"g=980.67 \t\t\t#gravitational field\n",
- "#Calculation\n",
"\n",
"\n",
"a=m2*R*T*p \n",
"b=V*9.9*d*g \n",
"M2=a/b #molar mass of macromolecule\n",
- "#Result\n",
"print\" M2 = molar mass of macromolecule , therefore M2 = \",round(M2),\"g.mol^-1\""
],
"language": "python",
@@ -446,21 +394,16 @@ "cell_type": "code",
"collapsed": false,
"input": [
- "#example 6.18\n",
"\n",
- "#To find the osmotic pressure of a solution\n",
- "#Variable declaration\n",
"R=82 \t\t\t#Universal gas constant[atm.ml.K^-1.mol^-1]\n",
"T=298 \t\t\t#Temperature[K]\n",
"V=250 \t\t\t#volume of water[ml]\n",
"m2=2.6 \t\t\t#mass of the protein\n",
"M2=85000 \t\t#molar mass of protein[g.mol^-1]\n",
"\n",
- "#Calculation\n",
"\t\n",
"n2=m2/M2 \t\t\t#no. of moles of protein\n",
"II=(n2*R*T)/V \t\t\t#Osmotic pressure of a solution[atm]\n",
- "#Result\n",
"print\"The osmotic pressure is\",round(II,5),\"atm \"\n"
],
"language": "python",
@@ -488,19 +431,14 @@ "cell_type": "code",
"collapsed": false,
"input": [
- "#example 6.19\n",
- "#To find the Ebullioscopic constant of water\n",
- "#Variable declaration\n",
"R=8.314 \t\t\t#Universal gas constant[J.K**-1.mol**-1]\n",
"Tb=373.15 \t\t\t#Boiling point temperature[K]\n",
"M1=0.018 \t\t\t# mass of water[kg]\n",
"Hvap=40.7 \t\t\t#Enthalpy of vaporization[KJ.mol**-1]\n",
- "#Calculation\n",
"\n",
"a=R*0.001*Tb**2*M1 \n",
"b=Hvap \n",
"Kb=a/b \t\t\t\t#Ebullioscopic constant of water[K.Kg.mol**-1]\n",
- "#Result\n",
"print\"The Ebullioscopic constant of water is\",round(Kb,2),\"K.Kg.mol-1\""
],
"language": "python",
@@ -528,19 +466,14 @@ "cell_type": "code",
"collapsed": false,
"input": [
- "#example 6.20\n",
- "#To find the value of activity coefficient\n",
"print\"CaF2(s)<=>CaF2(aq)<=>Ca+2(aq) + 2F-(aq)\"\n",
"\n",
- "#Variable declaration\n",
"Ksp=4.0*(10**-11) \t#Solubility product of sparingly soluble salt CaF2\n",
- "#Calculation\n",
"\n",
"x=Ksp/4.0 \n",
"Cs=x**(1.0/3.0) \t\t#Solubility \n",
"y=Cs**2 \n",
"Y=(x/y)**(1.0/3.0) \t\t#activity coefficient\n",
- "#Result\n",
"print\"The activity coefficient is\",Y \n",
"print\"NOTE:please note that the value of Cs is wrongly calculated as 4.64*10^-11 in book\""
],
@@ -571,10 +504,7 @@ "cell_type": "code",
"collapsed": false,
"input": [
- "#example 6.21\n",
"\n",
- "#To find the mean activity coefficient of ZnCl2 solution\n",
- "#Variable declaration\n",
"R=8.314 \t\t\t#Universal gas constant[J/K/mol]\n",
"T=298 \t\t\t\t#Temperature[K]\n",
"F=96500 \t\t\t#Faraday's constant\n",
@@ -582,13 +512,11 @@ "E=1.16 \t\t\t\t#e.m.f of the cell[Volts]\n",
"m=0.01 \n",
"import math\n",
- "#Calculation\n",
"\n",
"a=R*T \n",
"b=2*F \n",
"x=a/b \n",
"Y=math.exp((Eo-E-(x*math.log(4*m*m*m)))/(3*x)) #mean activity coefficient\n",
- "#Result\n",
"print\"The mean activity coefficient is\",round(Y,2) "
],
"language": "python",
@@ -616,9 +544,6 @@ "cell_type": "code",
"collapsed": false,
"input": [
- "#example 6.22\n",
- "#To find the ionic strength in a solution\n",
- "#Variable declaration\n",
"M1=0.01 \t\t\t#no. of moles of KCl\n",
"M2=0.005 \t\t\t#no. of moles of MgCl2\n",
"M3=0.002 \t\t\t#no. of moles of MgSO4\n",
@@ -629,14 +554,12 @@ "z22=1 \n",
"z31=2 \n",
"z32=2 \n",
- "#Calculation\n",
"\t\n",
"m1=M1/M \t\t\t#molality of KCL[m]\n",
"m2=M2/M \t\t\t#molality of MgCl2[m]\n",
"m3=M3/M \t\t\t#molality of MgSO4[m]\n",
"\n",
"I=0.5*((m1*z11**2+m1*z12**2+m2*z21**2+2*m2*z22**2+m3*z31**2+m3*z32**2)) #[mol/Kg]\n",
- "#Result\n",
"print\"The Ionic strength of a solution is\",I,\"mol/Kg\""
],
"language": "python",
@@ -664,22 +587,17 @@ "cell_type": "code",
"collapsed": false,
"input": [
- "#example 6.23\n",
- "#To find the mean activity coefficient\n",
- "#Variable declaration\n",
"T=298 \t\t\t\t#Temperature[K]\n",
"P=1 \t\t\t\t#pressure [atm]\n",
"m=0.02\t \t\t\t#Ionic strength of HCl solution in CH3OH[mol/Kg]\n",
"E=32.6 \t\t\t\t#Di-electric constant\n",
"d=0.787 \t\t\t#Density[gm/cm3]\n",
- "#Calculation\n",
"\t\n",
"I=0.5*(0.02*1*1+0.02*1*1) \t#Ionic strength of HCl solution[mol/Kg]\n",
"a=I*d \n",
"b=(E**3)*(298**3) \n",
"x=(a/b)**0.5 \n",
"Y=10**(-1.825*1000000*1*1*x) \t#mean activity coefficient\n",
- "#Result\n",
"print\"The mean activity coefficient is\",round(Y,2) "
],
"language": "python",
|