1 files changed, 57 insertions, 0 deletions
diff --git a/3472/CH15/EX15.4/Example15_4.sce b/3472/CH15/EX15.4/Example15_4.sce
new file mode 100644
index 000000000..adf1ede3f
--- /dev/null
+++ b/3472/CH15/EX15.4/Example15_4.sce
@@ -0,0 +1,57 @@
+// A Texbook on POWER SYSTEM ENGINEERING
+// A.Chakrabarti, M.L.Soni, P.V.Gupta, U.S.Bhatnagar
+// DHANPAT RAI & Co.
+// SECOND EDITION 
+
+// PART II : TRANSMISSION AND DISTRIBUTION
+// CHAPTER 8: CORONA
+
+// EXAMPLE : 8.4 :
+// Page number 228-229
+clear ; clc ; close ; // Clear the work space and console
+
+// Given data
+V = 110.0            // Operating line voltage(kV)
+f = 50.0             // Frequency(Hz)
+l = 175.0            // Line length(km)
+d = 1.0              // Diameter of conductor(cm)
+D = 300.0            // Distance b/w conductor(cm)
+t = 26.0             // Temperature(°C)
+b = 74.0             // Barometric pressure(cm)
+m = 0.85             // Irregularity factor
+m_v_local = 0.72     // Roughness factor for local corona
+m_v_gen = 0.82       // Roughness factor for general corona
+
+// Calculations
+delta = 3.92*b/(273.0+t)                                                  // Air density factor
+r = d/2.0                                                                 // Radius of conductor(cm)
+E_0 = 21.1*m*delta*r*log(D/r)                                             // Critical disruptive voltage(kV) rms
+E_v_local = 21.1*m_v_local*delta*r*(1+(0.3/(delta*r)**0.5))*log(D/r)      // Critical disruptive voltage for local corona(kV) rms
+E_v_gen = 21.1*m_v_gen*delta*r*(1+(0.3/(delta*r)**0.5))*log(D/r)          // Critical disruptive voltage for general corona(kV) rms
+E = V/3**0.5                                                              // Phase voltage(kV)
+// Case(i)
+P_c_i = 244.0*10**-5*(f+25)/delta*(r/D)**0.5*(E-E_0)**2                   // Peek"s formula for fair weather condition(kW/km/phase)
+P_c_i_total = P_c_i*l*3                                                   // Total power loss(kW)
+// Case(ii)
+P_c_ii = 244.0*10**-5*(f+25)/delta*(r/D)**0.5*(E-0.8*E_0)**2              // Peek"s formula for stormy condition(kW/km/phase)
+P_c_ii_total = P_c_ii*l*3                                                 // Total power loss(kW)
+// Case(iii)
+F_iii = 0.0713                                                            // From text depending on E/E_0
+P_c_iii = 21.0*10**-6*f*E**2*F_iii/(log10(D/r))**2                        // Peterson"s formula for fair condition(kW/km/phase)
+P_c_iii_total = P_c_iii*l*3                                               // Total power loss(kW)
+// Case(iv)
+F_iv = 0.3945                                                             // From text depending on E/E_0
+P_c_iv = 21.0*10**-6*f*E**2*F_iv/(log10(D/r))**2                          // Peterson"s formula for stormy condition(kW/km/phase)
+P_c_iv_total = P_c_iv*l*3                                                 // Total power loss(kW)
+    
+// Results
+disp("PART II - EXAMPLE : 8.4 : SOLUTION :-")
+printf("\nCase(i)  : Power loss due to corona using Peek formula for fair weather condition, P_c  = %.3f kW/km/phase", P_c_i)
+printf("\n           Total corona loss in fair weather condition using Peek formula = %.1f kW", P_c_i_total)
+printf("\nCase(ii) : Power loss due to corona using Peek formula for stormy weather condition, P_c  = %.2f kW/km/phase", P_c_ii)
+printf("\n           Total corona loss in stormy condition using Peek formula = %.f kW", P_c_ii_total)
+printf("\nCase(iii): Power loss due to corona using Peterson formula for fair weather condition, P_c  = %.4f kW/km/phase", P_c_iii)
+printf("\n           Total corona loss in fair condition using Peterson formula = %.2f kW",P_c_iii_total)
+printf("\nCase(iii): Power loss due to corona using Peterson formula for fair weather condition, P_c  = %.4f kW/km/phase", P_c_iv)
+printf("\n           Total corona loss in stormy condition using Peterson formula = %.1f kW \n",P_c_iv_total)
+printf("\nNOTE: ERROR: Calculation mistake in the final answer in textbook")