Ray/others, In order to test the power loss minimization, I have tried modeling the generators with equal linear costs simulating substation busbars.
Please have a look at the attached m-file, it contains a simple test case consisting of two generators connected to a load, with different branch impedance values. Theoretically, minimizing losses would mean that the generator (substation busbar) connected to the load through the branch of lowest impedance would be serving the load at full power (Pmax), right? Can you please explain why it is not? Thanks, Guilherme
function mpc = SStests
mpc.version = '2';
mpc.baseMVA = 100;
%% bus data
% b t P Q G B a V V
b z V V
% u y d d s s r m a a o m i
% s p e s n a n
% _ e a e e x
% i K
% V
mpc.bus = [
1 3 0 0 0 0 1 1 0
13.8 1 1.1 0.9;
2 2 0 0 0 0 1 1 0 13.8 1 1.1 0.9;
3 1 50 20 0 0 1 1 0 13.8 1 1.1 0.9;
];
%% generator data
% b P Q Q Q V m s P P
P P Q Q Q Q r r r r
a
% u g g m m g B t m m c c c c c c a a a
a p
% s a i a a a i 1 2 1 1 2 2 m m m
m f
% x m s t x n m m m m p p p
p
% e u i a i a _ _ _
_
% s n x n x a 1 3
q
% g 0 0
% c
mpc.gen = [
1 30 10 30 -30 1 100 1 35
0 0 0 0 0 0 0 0 0 0
0 0;
2 20 10 30 -30 1 100 1 35 0
0 0 0 0 0 0 0 0 0 0
0;
];
%% branch data
% f t r x b r r r r
a s a a
% b b a a a a n t n n
% u u t t t t g a g g
% s s e e e i l t m m
% A B C o e u i a
% s n x
mpc.branch = [
1 3 0.001 0.01 0 100 100 100 0 0
1 -360 360;
2 3 0.90 0.01 0 100 100 100 0 0 1 -360 360;
];
%%----- OPF Data -----%%
%% area data
% a r
% r e
% e f
% a b
% u
% s
mpc.areas = [
1 1;
];
%% generator cost data
% 1 startup shutdown n x1 y1 ... xn
yn
% 2 startup shutdown n c(n-1) ... c0
mpc.gencost = [
1 0 0 2 0 10 600 10;
1 0 0 2 0 10 600 10;
];
