Merge pull request #1 from efantnu/new-inertia

Correct inertia of turbo generators

README.md

@@ -6,12 +6,12 @@ This repository contains the data set and simulation files of the paper "Sizing
 # File organization
 - energy_storage.slx: Simulink file containing the surrogate model of the case study presented in the section "Sizing validation"
 - energy_storage_pre.m: MATLAB script that should be executed before running the Simulink model. Contains the parameters of all equipment and simulation options.
-- energy_storage_post.m: MATLAB script that should be executed after running the Simulink model. It produces the datasets required for Figures 8 and 9. It also calculates the energy supplied by the battery system.  
+- energy_storage_post.m: MATLAB script that should be executed after running the Simulink model. It produces the datasets required for Figures 9 and 10. It also calculates the energy supplied by the battery system.  
 - load-pdf.txt: dataset used to produce Figure 6. 
-- results-step3-noess.txt: dataset from case 1 used to produce Figure 8. 
-- results-step3.txt: dataset from case 2 used to produce Figure 8. 
-- results-step12-noess.txt: dataset from case 3 used to produce Figure 9. 
-- results-step12.txt: dataset from case 4 used to produce Figure 9. 
+- results-step3-noess.txt: dataset from case 1 used to produce Figure 9. 
+- results-step3.txt: dataset from case 2 used to produce Figure 9. 
+- results-step12-noess.txt: dataset from case 3 used to produce Figure 10. 
+- results-step12.txt: dataset from case 4 used to produce Figure 10. 
 - VSCdesignpu.m: MATLAB function that designs the AC-side LC filter and DC-link capacitor of a grid-connected three-phase voltage source converter (VSC). Based on this data, it also defines optimal controller parameters for a current controller on the dq0 reference frame.
 
 All files were created using MATLAB Simulink R2018a, the Simscape Electrical Specialized Power Systems Toolbox and the [NTNU Power Systems Library (pwrsys-matlab)](https://github.com/efantnu/pwrsys-matlab).

energy_storage_post.m

@@ -1,7 +1,8 @@
 %% Script Parameters
 tStart = 0;
-tEnd = 340;
-saveResults = 1;
+tEnd = simu.tTotal;
+saveResults = 0;
+saveTxt = 0;
 
 switch(caseNumber)
     case 1
@@ -66,24 +67,24 @@
 PlotTable.PWT = Plot.PWT.Data;
 
 
-if saveResults
+if saveTxt
     filename = strcat('results',caseDescription,'.txt');
     writetable(PlotTable,filename,'Delimiter','\t');
 end
 
 %% Energy calculations
 rtr = 0.05;     % Maximum transient frequency deviation [pu]
-rss = 0.02;     % Maximum steady-state frequency deviation [pu]
-Tarr = 8; %25;      % Arrest period [s]
-Treb = 120; %15;      % Rebound period [s]
-Trec = 165; %120;     % Recovery period [s]
+rss = 0.02; %0.03;     % Maximum steady-state frequency deviation [pu]
+Tarr = 11; %3.5;       % Arrest period [s]
+Treb = 17; %120;       % Rebound period [s]
+Trec = 150; %230;      % Recovery period [s]
 
 % Proposed
 EES1arr = 1/12 * Tarr * (2*rtr^2 + 3*rtr);
 EES1reb = 1/6 * Treb * (2*rtr^3 + 3*rtr^2 - 2*rss^3 - 3*rss^2) / (rtr - rss);
 EES1rec = 1/12 * Trec * (2*rss^2 + 3*rss);
 
-ES1calc = (es1.Prated/rss)/3600 * (EES1arr + EES1reb + EES1rec);
+EES1calc = (es1.Prated/rss)/3600 * (EES1arr + EES1reb + EES1rec);
 
 % Simulation
 EES1sim = trapz(0.1/3600,Plot.PES1.Data)*grid.Sb;

energy_storage_pre.m

@@ -118,7 +118,7 @@
 gt1.SM.TdppSC = 0.022;          % d-axis subtransient time - short-circuit [s]
 gt1.SM.TqppSC = 0.0334;         % q-axis subtransient time - short-circuit [s]
 gt1.SM.Rs = 0.0242;             % Stator (armature) resistance [pu]
-gt1.SM.M = 2*5.1*70/88;         % Inertia constant [s]
+gt1.SM.M = 5.1*88/70;           % Inertia constant [s]
 gt1.SM.p = 2;                   % Pairs of poles
 
 % Turbine