initial commit / add all books

16 files changed, 255 insertions, 0 deletions

diff --git a/3872/CH12/EX12.1/Ex12_1.JPG b/3872/CH12/EX12.1/Ex12_1.JPG
new file mode 100644
index 000000000..407efe386
--- /dev/null
+++ b/3872/CH12/EX12.1/Ex12_1.JPG
diff --git a/3872/CH12/EX12.1/Ex12_1.sce b/3872/CH12/EX12.1/Ex12_1.sce
new file mode 100644
index 000000000..2240ec1c0
--- /dev/null
+++ b/3872/CH12/EX12.1/Ex12_1.sce
@@ -0,0 +1,30 @@
+//Book - Power System: Analysis & Design 5th Edition

+//Authors - J. Duncan Glover, Mulukutla S. Sarma, and Thomas J. Overbye

+//Chapter - 12 ; Example 12.1

+//Scilab Version - 6.0.0 ; OS - Windows

+

+clc;

+clear;

+

+Tr=0;

+Ka=100;

+Ta=0.05;

+Vrmax=5;

+Vrmin=-5;

+Ke=1;

+Te=0.26;

+Kf=0.01;

+Tf=1;

+

+Efd=2.9135;             //Value taken from Example 11.10

+Vt=1.0946;              //Value taken from Example 11.10

+

+Vr=Ke*Efd;              //Initial value of Vr

+Vf=0;                   //Initial value of vf

+Vref=(Vr/Ka)+Vt+Vf;     //Initial value of Vref

+

+printf('The initial value of Vr is %.4f\n',Vr)

+printf('The initial value of Vf is %.4f\n',Vf)

+printf('The initial value of Vref is %.4f\n',Vref)

+

+//Section 'b' of this problem cannot be simulated using current version of Scilab

diff --git a/3872/CH12/EX12.2/Ex12_2.JPG b/3872/CH12/EX12.2/Ex12_2.JPG
new file mode 100644
index 000000000..b1e885d20
--- /dev/null
+++ b/3872/CH12/EX12.2/Ex12_2.JPG
diff --git a/3872/CH12/EX12.2/Ex12_2.sce b/3872/CH12/EX12.2/Ex12_2.sce
new file mode 100644
index 000000000..edf79be4b
--- /dev/null
+++ b/3872/CH12/EX12.2/Ex12_2.sce
@@ -0,0 +1,35 @@
+//Book - Power System: Analysis & Design 5th Edition

+//Authors - J. Duncan Glover, Mulukutla S. Sarma, and Thomas J. Overbye

+//Chapter - 12 ; Example 12.2

+//Scilab Version - 6.0.0 ; OS - Windows

+

+clc;

+clear;

+

+KQi=0.4;

+KVi=40;

+XIqmax=1.45;

+XIqmin=0.5;

+Vmax=1.1;

+Vmin=0.9;

+

+Tr=0;

+Ka=100;

+Ta=0.05;

+Vrmax=5;

+Vrmin=-5;

+Ke=1;

+Te=0.26;

+Kf=0.01;

+Tf=1;

+vt=0.5;

+Vf=0;                   //Initial value of vf

+Vref=1.0239;            //Initial value of Vref from Example 11.12

+Isorq=-XIqmax/0.8;      // Reactive component of Isorc

+Qcmd=0.22;              //Obtained from Example 11.12

+Qnet=(vt)*abs(Isorq)-(vt)^2/0.8;   // Net reactive power injection in pu

+

+printf('The initial value of reference voltage is %.4f pu\n',Vref)

+printf('The initial value of reactive power Qcmd = %.4f pu = %.4f Mvar\n',Qcmd,Qcmd*100)

+printf('The maxximum net reactive power Qnet = %.4f pu = %.4f Mvar\n',Qnet,Qnet*100)

+

diff --git a/3872/CH12/EX12.3/Ex12_3.JPG b/3872/CH12/EX12.3/Ex12_3.JPG
new file mode 100644
index 000000000..310cac9cb
--- /dev/null
+++ b/3872/CH12/EX12.3/Ex12_3.JPG
diff --git a/3872/CH12/EX12.3/Ex12_3.sce b/3872/CH12/EX12.3/Ex12_3.sce
new file mode 100644
index 000000000..7b71197ab
--- /dev/null
+++ b/3872/CH12/EX12.3/Ex12_3.sce
@@ -0,0 +1,19 @@
+//Book - Power System: Analysis & Design 5th Edition

+//Authors - J. Duncan Glover, Mulukutla S. Sarma, and Thomas J. Overbye

+//Chapter - 12 ; Example 12.3

+//Scilab Version - 6.0.0 ; OS - Windows

+

+clc;

+clear;

+

+M=500;               //MVA rating of the generator

+f=60;                //requency in Hertz

+R=0.05;              //Regulation constant in pu

+delF=0.01;           //Increase in frequency in Hertz

+

+delFpu=delF/f;       //Frequency increase in pu

+delPref=0;           //Since fixed reference power setting is assumed

+delPmpu=delPref-(1/R)*delFpu    //Change in mechanical power in pu

+delPm=delPmpu*M;                //Actual value of mechanical power in MW

+

+printf('The turbine mechanical power output decreases by %.3f MW.',abs(delPm))

diff --git a/3872/CH12/EX12.4/Ex12_4.JPG b/3872/CH12/EX12.4/Ex12_4.JPG
new file mode 100644
index 000000000..d4245e76a
--- /dev/null
+++ b/3872/CH12/EX12.4/Ex12_4.JPG
diff --git a/3872/CH12/EX12.4/Ex12_4.sce b/3872/CH12/EX12.4/Ex12_4.sce
new file mode 100644
index 000000000..602114e9e
--- /dev/null
+++ b/3872/CH12/EX12.4/Ex12_4.sce
@@ -0,0 +1,31 @@
+//Book - Power System: Analysis & Design 5th Edition

+//Authors - J. Duncan Glover, Mulukutla S. Sarma, and Thomas J. Overbye

+//Chapter - 12 ; Example 12.4

+//Scilab Version - 6.0.0 ; OS - Windows

+

+clc;

+clear;

+funcprot(0);

+f=60;                          //Frequency in Hertz

+G=[1000 750 500];              //Rating of unit 1,2 &3 respectively in MVA

+R=0.05;                        //Regulation constant of each unit in pu

+delP=200;                      //Load increment in MW

+SBnew=1000;                    //New MVA base of the entire system

+

+Rnew=R*(SBnew./G);             //Regulation of each generators with common base

+beta=sum(1 ./Rnew);             //area frequency response characteristic, beta

+

+printf('The area frequency response characteristic beta is %.2f per unit\n',beta)

+

+delPpu=delP/SBnew;             //Load increment in per unit

+delFpu=(-1/beta)*delPpu        //Frequency drop in per unit

+delF=delFpu*f;                 //Frequency drop in Hertz

+

+printf('The steady-state frequency drop is %.4f Hz\n',abs(delF))

+

+delPm=delFpu*(-1 ./Rnew);

+delPmact=SBnew*delPm;

+

+printf('The increase in turbine mechanical power output of unit1=%.4f pu = %.4f MW\n',delPm(1),delPmact(1))

+printf('The increase in turbine mechanical power output of unit2=%.4f pu = %.4f MW\n',delPm(2),delPmact(2))

+printf('The increase in turbine mechanical power output of unit3=%.4f pu = %.4f MW',delPm(3),delPmact(3))

diff --git a/3872/CH12/EX12.5/Ex12_5.JPG b/3872/CH12/EX12.5/Ex12_5.JPG
new file mode 100644
index 000000000..bdd340193
--- /dev/null
+++ b/3872/CH12/EX12.5/Ex12_5.JPG
diff --git a/3872/CH12/EX12.5/Ex12_5.sce b/3872/CH12/EX12.5/Ex12_5.sce
new file mode 100644
index 000000000..054733001
--- /dev/null
+++ b/3872/CH12/EX12.5/Ex12_5.sce
@@ -0,0 +1,40 @@
+//Book - Power System: Analysis & Design 5th Edition

+//Authors - J. Duncan Glover, Mulukutla S. Sarma, and Thomas J. Overbye

+//Chapter - 12 ; Example 12.5

+//Scilab Version - 6.0.0 ; OS - Windows

+

+clc;

+clear;

+

+f=60;                    //Frequency of the system in Hertz

+G1=2000;                 //Total generation of area 1 in MW

+G2=4000;                 //Total generation of area 2 in MW

+beta1=700;               //Area frequency response characteristic of area 1 in MW/Hz

+beta2=1400;              //Area frequency response characteristic of area 2 in MW/Hz

+delPD1=100;              //Load increment of area 1 in MW

+delPD2=0;

+

+//WITHOUT LFC

+delF=(delPD1+delPD2)/(-(beta1+beta2));    //Frequency Change in Hertz

+delPm1=-beta1*delF;                       //Change in power of area 1

+delPm2=-beta2*delF;                       //Change in power of area 2

+delPtie1=-delPm2                          //Tie line power flow from area 1 to 2

+delPtie2=delPm2                          //Tie line power flow from area 2 to 1

+

+disp('RESULTS WITHOUT LFC')

+printf('\nThe steady state frequency error is %.4f Hz',delF)

+printf('\nThe tie-line power flow from area 1 is %.4f MW',delPtie1)

+printf('\nThe tie-line power flow from area 2 is %.4f MW\n',delPtie2)

+

+//WITH LFC

+

+delFl=0/(beta1+beta2);              //Frequency Change in Hertz (as ACE1+ACE2=0)

+delPm1=-beta1*delFl;                //Change in power of area 1

+delPm2=-beta2*delFl;                //Change in power of area 2

+delPtie1=-delPm2                    //Tie line power flow from area 1 to 2

+delPtie2=delPm2                     //Tie line power flow from area 2 to 1

+

+disp('RESULTS WITH LFC')

+printf('\nThe steady state frequency error is %.4f Hz',delFl)

+printf('\nThe tie-line power flow from area 1 is %.4f MW',delPtie1)

+printf('\nThe tie-line power flow from area 2 is %.4f MW',delPtie2)

diff --git a/3872/CH12/EX12.6/Ex12_6.JPG b/3872/CH12/EX12.6/Ex12_6.JPG
new file mode 100644
index 000000000..7956950a9
--- /dev/null
+++ b/3872/CH12/EX12.6/Ex12_6.JPG
diff --git a/3872/CH12/EX12.6/Ex12_6.sce b/3872/CH12/EX12.6/Ex12_6.sce
new file mode 100644
index 000000000..19fe7085a
--- /dev/null
+++ b/3872/CH12/EX12.6/Ex12_6.sce
@@ -0,0 +1,24 @@
+//Book - Power System: Analysis & Design 5th Edition

+//Authors - J. Duncan Glover, Mulukutla S. Sarma, and Thomas J. Overbye

+//Chapter - 12 ; Example 12.6

+//Scilab Version - 6.0.0 ; OS - Windows

+

+clc;

+clear;

+

+C1=[8e-3 10 0]                  //Coefficients of cost equation for unit 1

+C2=[9e-3 8  0]                  //Coefficients of cost equation for unit 2

+

+dC1=[2*C1(1) C1(2)]        //Coefficients of incremental cost equation for unit 1

+dC2=[2*C2(1) C2(2)]        //Coefficients of incremental cost equation for unit 2

+result=[];

+for PT=500:100:1500

+    P1=(dC2(1)*PT+(dC2(2)-dC1(2)))/(dC2(1)+dC1(1));

+    P2=PT-P1;

+    dC1value=dC1(1)*P1+dC1(2);

+    dC2value=dC2(1)*P2+dC2(2);

+    CT=C1(1)*P1^2+C1(2)*P1+C1(3)+C2(1)*P2^2+C2(2)*P2+C2(3);

+    result=[result;PT P1 P2 dC1value CT]

+end

+

+disp(result,' PT(MW)     P1(MW)      P2(MW)  dC1=dC2($/MWhr) CT($/hr)');

diff --git a/3872/CH12/EX12.7/Ex12_7.JPG b/3872/CH12/EX12.7/Ex12_7.JPG
new file mode 100644
index 000000000..1ecd121a2
--- /dev/null
+++ b/3872/CH12/EX12.7/Ex12_7.JPG
diff --git a/3872/CH12/EX12.7/Ex12_7.sce b/3872/CH12/EX12.7/Ex12_7.sce
new file mode 100644
index 000000000..6bf46aba1
--- /dev/null
+++ b/3872/CH12/EX12.7/Ex12_7.sce
@@ -0,0 +1,46 @@
+//Book - Power System: Analysis & Design 5th Edition

+//Authors - J. Duncan Glover, Mulukutla S. Sarma, and Thomas J. Overbye

+//Chapter - 12 ; Example 12.7

+//Scilab Version - 6.0.0 ; OS - Windows

+

+clc;

+clear;

+

+C1=[8e-3 10 0]              //Coefficients of cost equation for unit 1

+C2=[9e-3 8  0]              //Coefficients of cost equation for unit 2

+

+dC1=[2*C1(1) C1(2)]        //Coefficients of incremental cost equation for unit 1

+dC2=[2*C2(1) C2(2)]        //Coefficients of incremental cost equation for unit 2

+

+P1lim=[100 600];           //Lower and upper generation limit for unit 1

+P2lim=[400 1000];          //Lower and upper generation limit for unit 2

+

+result=[];

+for PT=[500 600 700 725 800 900 1000 1100 1200 1244 1300 1400 1500]

+    P1=(dC2(1)*PT+(dC2(2)-dC1(2)))/(dC2(1)+dC1(1));

+    P2=PT-P1;

+    dC1value=dC1(1)*P1+dC1(2);

+    dC2value=dC2(1)*P2+dC2(2);

+            

+    if P1<P1lim(1) | P1>P1lim(2)       //Checking for limits of P1

+        if P1<P1lim(1)

+            P1=P1lim(1)

+        else

+            P1=P1lim(2)

+        end

+        P2=PT-P1;

+        dC1value=dC2(1)*P2+dC2(2);

+    elseif P2<P2lim(1) | P2>P2lim(2)   //Checking for limits of P2

+        if P2<P2lim(1)

+            P2=P2lim(1)

+        else

+            P2=P2lim(2)

+        end

+        P1=PT-P2;

+        dC1value=dC1(1)*P1+dC1(2);

+    end

+    

+    CT=C1(1)*P1^2+C1(2)*P1+C1(3)+C2(1)*P2^2+C2(2)*P2+C2(3);  //Total cost in $/hr

+    result=[result;PT P1 P2 dC1value CT]

+end

+disp(result,' PT(MW)     P1(MW)      P2(MW)  dC/dP($/MWhr) CT($/hr)');

diff --git a/3872/CH12/EX12.9/Ex12_9.JPG b/3872/CH12/EX12.9/Ex12_9.JPG
new file mode 100644
index 000000000..3a323446a
--- /dev/null
+++ b/3872/CH12/EX12.9/Ex12_9.JPG
diff --git a/3872/CH12/EX12.9/Ex12_9.sce b/3872/CH12/EX12.9/Ex12_9.sce
new file mode 100644
index 000000000..6cd8eba6e
--- /dev/null
+++ b/3872/CH12/EX12.9/Ex12_9.sce
@@ -0,0 +1,30 @@
+//Book - Power System: Analysis & Design 5th Edition

+//Authors - J. Duncan Glover, Mulukutla S. Sarma, and Thomas J. Overbye

+//Chapter - 12 ; Example 12.9

+//Scilab Version - 6.0.0 ; OS - Windows

+

+clc;

+clear;

+

+B11=1.5e-4; B12=2e-5; B22=3e-5;               // Loss coefficients

+lamda=16;                                     //Area incremental cost in $/MWhr

+

+e1=[20.8e-3 32e-5 6];       //Coefficients of incremental operating cost equation1

+e2=[32e-5 18.96e-3 8];      //Coefficients of incremental operating cost equation2

+

+P1=(e2(2)*e1(3)-e1(2)*e2(3))/(e2(2)*e1(1)-e1(2)*e2(1));   //Solution of P1 from incremental cost equations

+P2=(e2(1)*e1(3)-e1(1)*e2(3))/(e1(2)*e2(1)-e2(2)*e1(1));   //Solution of P2 from incremental cost equations

+

+Pl=B11*P1^2+B12*P1*P2+B22*P2^2;            //Total losses

+

+Pt=P1+P2-Pl;                               //Total demand

+

+CT=10*P1+8e-3*P1^2+8*P2+9e-3*P2^2;        // Cost equation taken from example 12.6

+

+printf('The output of each unit are given by P1=%d MW and P2=%d MW\n',P1,P2)

+printf('The total transmission loss is %.2f MW\n',Pl)

+printf('The total demand is %.2f MW\n',Pt)

+printf('The total operation cost is %.2f $/hr',CT)

+

+

+