8 files changed, 226 insertions, 0 deletions
diff --git a/3523/CH5/EX5.6.10/Ex5_10.sce b/3523/CH5/EX5.6.10/Ex5_10.sce
new file mode 100644
index 000000000..003adc9cc
--- /dev/null
+++ b/3523/CH5/EX5.6.10/Ex5_10.sce
@@ -0,0 +1,29 @@
+//Example 10// Ch 5
+clc;
+clear;
+close;
+// given data
+m1=0.92;//smoothness coefficient
+m2=0.95;//weather coefficient
+Deq=600;//mean geometric distance b/w conductors in cm
+V = 275;//line operating at voltage V in kV
+p=75;//pressure in cm Hg
+t = 35;//in degree C
+r=1;//radius of conductors in cm
+delta=3.92*p/(273+t);//relative air density
+printf("relative air density %f",delta)
+E0=30*delta*(1+0.3/sqrt(delta*r))*m1*m2;//corona onset field
+printf("corona onset field %f kVpeak/cm",E0)
+V0 = E0*log(Deq);//onset voltage in kVpeak
+printf("onset voltage %f kVpeak",V0)
+V0rms=V0/sqrt(2);//rms onset voltage
+printf("rms onset voltage %f kV",V0rms)
+V0ll=V0rms*sqrt(3);//onset voltage line to line
+printf("line to line onset voltage %f kV line to line",V0ll)
+K= 0.05;
+f=50;//in Hz
+Vph=(V*10^3)/sqrt(3);
+Pc=3.73*K*f*(Vph^2)*10^-5/(Deq/r)^2;
+printf("corona power loss %f kW/(cond.km)",Pc)
+Ic=Pc/Vph;
+printf("corona current %e A/km",Ic)
diff --git a/3523/CH5/EX5.6.11/Ex5_11.sce b/3523/CH5/EX5.6.11/Ex5_11.sce
new file mode 100644
index 000000000..ef34232e8
--- /dev/null
+++ b/3523/CH5/EX5.6.11/Ex5_11.sce
@@ -0,0 +1,26 @@
+//Example 10// Ch 5
+clc;
+clear;
+close;
+// given data
+m1=0.9;//smoothness coefficient
+m2=0.9;//weather coefficient
+r=3.175;//radius of conductor in cm
+V=525;//rated voltage in kV where no corona is present
+delta=1;//relative air-density factor
+Deq=112.63;//in cm
+E0=30*delta*(1+0.3/sqrt(delta*r))*m1*m2;//corona onset field
+printf("corona onset field %f kVpeak/cm",E0)
+E0rms=E0/sqrt(2);
+printf("rms corona onset field %f kV/cm",E0rms)
+V0=E0*r*log(Deq);
+printf("corona onset voltage %f kV",V0)
+V0ll=V0*sqrt(3);
+printf("corona onset voltage lin to line %f kV",V0ll)
+V1=2.5*V;//line to line voltage higher than V0 so corona is present on the conductor
+re=5;//effective radius of corona envelope in cm
+printf("envelope radius %f cm",re)
+
+
+
+
diff --git a/3523/CH5/EX5.6.2/Ex5_2.sce b/3523/CH5/EX5.6.2/Ex5_2.sce
new file mode 100644
index 000000000..7a075768e
--- /dev/null
+++ b/3523/CH5/EX5.6.2/Ex5_2.sce
@@ -0,0 +1,22 @@
+//Example 2// Ch 5
+clc;
+clear;
+close;
+// given data
+d = 0.001;//in meters
+p = 101.3; //gas pressure in kPa
+C = -2400.4;//constant value
+A = 0.027;//constant value
+As = 10^8;//avalanche size
+//secondary ionization coefficient is much smaller than unity therefore ionization coefficient (alpha) is equal to electron attachment coefficient 
+E = (2400.4*p)/0.027; //alpha is equal to e- attachment coeff occurs at this eq
+Vs1 = E*d;//breakdown voltage in V
+printf("electric field %e V/m \n",E)
+printf("breakdown voltage %f V \n",Vs1)
+Vs2 = (log(As)-C*p*d)/A; //in V
+printf("breakdown voltage corresponding to an avalanche size %f V \n",Vs2)
+//as the avalanche self-space charge is neglected the breakdown voltage will be same irrespective of the polarity of the stressed plate acc. to eq (5.4) N2>=N1;
+Vspos = 9.4;//in kV when N2>=N1 in which no of e- in second avalanche is greater than equal to no of e- in first avalanche
+printf("positive voltage breakdown %f kV \n",Vspos)
+Vsneg = 9.2;//in kV when Neph >= 1 where Neph is no of e-photoemitted from the cathode
+printf("negative voltage breakdown %f kV \n",Vsneg)
diff --git a/3523/CH5/EX5.6.3/Ex5_3.sce b/3523/CH5/EX5.6.3/Ex5_3.sce
new file mode 100644
index 000000000..a3e5542b5
--- /dev/null
+++ b/3523/CH5/EX5.6.3/Ex5_3.sce
@@ -0,0 +1,28 @@
+//Example 3// Ch 5
+clc;
+clear;
+close;
+// given data
+d = 0.001;//in meters
+p1 = 3*101.3; //gas pressure of 3 atmp in kPa
+p2 = 5*101.3; //gas pressure of 5 atmp in kPa
+C =  2400.4;//constant value
+A = 0.027;//constant value
+As = 10^8;//avalanche size
+Vs1 = C*p1*d/A;//breakdown voltage at 3 atm
+Vs2 = C*p2*d/A;//breakdown voltage at 5 atm
+Vs3 = (log(As)+C*p1*d)/A;//breakdown voltage at 3 atm corresponding to an avalanche size
+Vs4 = (log(As)+C*p2*d)/A;//breakdown voltage at 5 atm corresponding to an avalanche size
+printf("breakdown voltage at 3 atm %f kV \n",Vs1*10^-3)
+printf("breakdown voltage at 5 atm %f kV \n",Vs2*10^-3)
+printf("breakdown voltage at 3 atm corresponding to an avalanche size %f kV \n",Vs3*10^-3)
+printf("breakdown voltage at 5 atm corresponding to an avalanche size %f kV \n",Vs4*10^-3)
+//acc. to eq N2>=N1 and Neph>=1 with increase of gas pressure improves the dielectric strength of the gas since breakdown voltage increses with gas pressure
+Vs1pos = 27.5;//postive breakdown voltage at 3 atm in kV
+Vs1neg = 27.73;//negative breakdown voltage at 3 atm in kV
+Vs2pos = 45.2;//postive breakdown voltage at 5 atm in kV
+Vs2neg = 45.5;//negative breakdown voltage at 5 atm in kV
+printf("positive breakdown voltage at 3 atm %f kV \n",Vs1pos)
+printf("negative breakdown voltage at 3 atm %f kV \n",Vs1neg)
+printf("positive breakdown voltage at 5 atm %f kV \n",Vs2pos)
+printf("negative breakdown voltage at 5 atm %f kV \n",Vs2neg)
diff --git a/3523/CH5/EX5.6.5/Ex5_5.sce b/3523/CH5/EX5.6.5/Ex5_5.sce
new file mode 100644
index 000000000..761be0065
--- /dev/null
+++ b/3523/CH5/EX5.6.5/Ex5_5.sce
@@ -0,0 +1,23 @@
+//Example 5// Ch 5
+clc;
+clear;
+close;
+// given data
+d=0.001;
+a = 0.1*10^-2;//radii of concentric circle in meters
+b = 2.1*10^-2;//radii of concentric circle in meters
+p = 101.3;//gas pressure in kPa
+p1=3*p;
+p2=5*p;
+C =  -2400.4;//constant value
+A = 0.027;//constant value
+As = 10^8;//avalanche size
+ri = 0.0772;//in m
+V0 = [log(10^8)-{(C*p)*(ri-a)}]*(b-a)/[A*{(1/a)-(1/ri)}];
+printf("corona onset voltage is %f kV \n",V0)
+V0pos = 13.1;//in kV
+V0neg = 13.7;//in kV
+printf("positive corona onset voltage %f kV \n",V0pos)
+printf("negative corona onset voltage %f kV \n",V0neg)
+
+//acc. to eq N2>=N1 and Neph>=1 with increase of gas pressure improves the dielectric strength of the gas since breakdown voltage increses with gas pressure
diff --git a/3523/CH5/EX5.6.6/Ex5_6.sce b/3523/CH5/EX5.6.6/Ex5_6.sce
new file mode 100644
index 000000000..13ef66fe7
--- /dev/null
+++ b/3523/CH5/EX5.6.6/Ex5_6.sce
@@ -0,0 +1,30 @@
+//Example 6// Ch 5
+clc;
+clear;
+close;
+// given data
+d=0.001;
+a = 0.1*10^-2;//radii of concentric circle in meters
+b = 2.1*10^-2;//radii of concentric circle in meters
+p = 101.3;//gas pressure in kPa
+p1=3*p;
+p2=5*p;
+C =  -2400.4;//constant value
+A = 0.027;//constant value
+As = 10^8;//avalanche size
+ri = 0.0772;//in m
+V01 = [log(10^8)-{(C*p1)*(ri-a)}]*(b-a)/[A*{(1/a)-(1/ri)}];
+V02 = [log(10^8)-{(C*p2)*(ri-a)}]*(b-a)/[A*{(1/a)-(1/ri)}];
+printf("corona onset voltage at 3atmp is %f kV \n",V01)
+printf("corona onset voltage at 5atmp is %f kV \n",V02)
+V01pos = 41.9;//in kV at 3 atmp
+V01neg = 42.2;//in kV at 3 atmp
+V02pos = 69.2;//in kV at 5 atmp
+V02neg = 69.8;//in kV at 5 atmp
+printf("positive corona onset voltage %f kV \n",V01pos)
+printf("negative corona onset voltage %f kV \n",V01neg)
+printf("positive corona onset voltage %f kV \n",V02pos)
+printf("negative corona onset voltage %f kV \n",V02neg)
+//answer given in the book is wrong
+
+//acc. to eq N2>=N1 and Neph>=1 with increase of gas pressure improves the dielectric strength of the gas since breakdown voltage increses with gas pressure
diff --git a/3523/CH5/EX5.6.8/Ex5_8.sce b/3523/CH5/EX5.6.8/Ex5_8.sce
new file mode 100644
index 000000000..06c969de8
--- /dev/null
+++ b/3523/CH5/EX5.6.8/Ex5_8.sce
@@ -0,0 +1,54 @@
+//Example 8// Ch 5
+clc;
+clear;
+close;
+// given data
+
+delta=1;//at standard temp and pressure
+r=1;//radius of conductors in cm
+s=40;//subconductor to subconductor spacing in cm
+D=500; //phase to phase spacing in cm
+E0=30*delta*(1+(0.3/sqrt(delta*r)));//corona onset field in kVpeak/cm
+printf("corona onset field %f kVpeak/cm",E0)
+
+V01=E0*log(D/r);//corona onset voltage using single conductor
+printf("corona onset voltage V01 is %f kVpeak",V01)
+V01rms=V01/sqrt(2);//rms onset voltage in kV
+printf("corona rms onset voltage V01rms %f kV",V01rms)
+
+x2 = log(D /(sqrt(s*r)));
+y2 = (1+((2*r)/s));
+
+V02=2*E0*r*(x2/y2);//corona onset voltage using bundle-2 conductor arranged horizontally and vertically
+printf("corona onset voltage V02 is %f kVpeak",V02)
+V02rms=V02/sqrt(2);//rms onset voltage in kV
+printf("corona rms onset voltage V02rms is %f kV",V02rms)
+
+
+x3 = log(D /((sqrt(2)*(s)^2*r)^0.3));
+y3 = (1+((3*sqrt(3)*r)/s));
+
+V03=3*E0*r*(x3/y3);//corona onset voltage using bundle-3 conductor arranged at vertices of an upright or inverted triangle
+printf("corona onset voltage V03 is %f kVpeak",V03)
+V03rms=V03/sqrt(2);//rms onset voltage in kV
+printf("corona rms onset voltage V03rms is %f kV",V03rms)
+
+
+x4 = log(D /((sqrt(2)*(s)^3*r)^0.25));
+y4 = (1+((4*sqrt(2)*r)/s));
+
+V04=4*E0*r*(x4/y4);//corona onset voltage using bundle-4 conductor arranged at vertices of a square
+printf("corona onset voltage V04 is %f kVpeak",V04)
+V04rms=V04/sqrt(2);//rms onset voltage in kV
+printf("corona rms onset voltage V04rms is %f kV",V04rms)
+
+
+x5 = log(D /((sqrt(2)*(s)^3*r)^0.25));
+y5 = (1+((3*sqrt(2)*r)/s));
+
+V05=4*E0*r*(x5/y5);//corona onset voltage using bundle-4 conductor arranged at vertices of a diamond form square
+printf("corona onset voltage V05 is %f kVpeak",V05)
+V05rms=V05/sqrt(2);//rms onset voltage in kV
+printf("corona rms onset voltage V05rms is %f kV",V05rms)
+
+//acc. to eq 5.18 in question 7 corona onset voltage is calculated
diff --git a/3523/CH5/EX5.6.9/Ex5_9.sce b/3523/CH5/EX5.6.9/Ex5_9.sce
new file mode 100644
index 000000000..e8f7c9910
--- /dev/null
+++ b/3523/CH5/EX5.6.9/Ex5_9.sce
@@ -0,0 +1,14 @@
+//Example 9// Ch 5
+clc;
+clear;
+close;
+// given data
+Deq=600;//mean geometric distance b/w conductors in cm
+delta=1;//at standard temp and pressure
+r=1;//radius of conductors in cm
+E0=30*delta*(1+(0.3/sqrt(delta*r)));//corona onset field in kVpeak/cm
+printf("corona onset field %f kVpeak/cm",E0)
+V0=E0*log(Deq);//corona onset voltage
+printf("corona onset voltage %f kVpeak",V0)
+V0rms=V0/sqrt(2);//rms onset voltage in kV
+printf("corona rms onset voltage %f kV",V0rms)