Dear Igor, please see comments below.
Igor Verbruggen wrote:
Dear all,
I am attempting to use
MOST to formulate a DC security constrained OPF, with
preventive and corrective redisp
naime ahmadi wrote:
Hello,
I want to calculate shunt capacitance (c) in the
attached figure.
In the branch data of MATPOWER we have b.
Is it correct to use this formulation to calculate c
c=b/(2*pi*f)
Slight correction:
Sf = V( mpc.branch(:, F_BUS) ) .* conj(Yf * V);
St = V( mpc.branch(:, T_BUS) ) .* conj(Yf * V);
carlos.
Ray Daniel Zimmerman wrote:
In matrix form, you can do the corresponding calculations as
Ahmad Bariq Al Fahri wrote:
Thank you Carlos for the information.
As in Case 9, all bus voltage is specified to be 345 kV. If
the transformer defaults ratio is 1:1, does it mean that the
generator is connected in 345 kV
Ahmad Bariq Al Fahri wrote:
I would like to ask about case9.m (9 bus, 3 generators
case) system configuration. As I know that the parameters of
'TAP' (column 9) in the branch data are set to zero, does it
mean that the system does no
naime ahmadi wrote:
Hello,
How can I change
the power
injection
to flow
measurement?
For
example, in the IEEE 14-bus system, bus 2 is connected to
buses 1, 3
Harun Hökelek wrote:
Dear all,
I’m trying to obtain bus impedance matrix of IEEE 33-Bus test system by
inversing admittance matrix.
However inv(makeYbus(case-data)) gives badly scaled or singular warning.
Is there any way to get the bus impedance matrix of this system?
Sincerely,
Harun
Dear Ha
naime ahmadi wrote:
Dear Professor Ray,
I have a 3 bus system and want to calculate the DC Jacobian. I
attached the result that I want. My problem is that I do not know
which function in MATPOWER does exactly this calculation.
I also attached the 3 bus system in MATPOWER format.
Look at the f
Müller, Tobias wrote:
Hello everyone,
I have problems when I
have a piecewise linear cost function on some generators
which increases slower after a break point as the power
increas
Matthias:
You are specifying piece-wise linear costs with corner points at
(100, 1000) and (110,1000). That's a single flat linear segment,
which of course will have a marginal cost of zero. In general, if
these flexible generators are not against a limit
yangyang wrote:
Thank you Carlos, I have also found it; the tap should be the
ratio of per unit voltages instead of nominal voltages. But
after changing the tap to 1, the power flow shows that all buses
are 1 pu 0 degree, and it still differ
yangyang wrote:
Hi Carlos, I think there is a bus labeled 1 in the bus field,
next to the "mpc.bus" in the file, and it is a slack bus. If
there was not a slack bus, runpf() will report an error. So this
may not be the problem of my case...
yangyang wrote:
Dear all,
I am trying to construct the
power flow of IEEE 2nd benchmark for SSR study; and my
implementation is as attached. The capacity compensation is
55% and the other data can be found in the paper
Hi all;
Actually, what you need to do is to set equal linear costs for
all generators. If all linear cost coefficients are 1.0, for
example, then you get to minimize the sum of all generation, which
is a proxy for minimizing losses.
carlos.
Russ Patterson wrote:
Hi
- I am still trying to hand calculate the flow
into branch 2 from bus 1 to bus 2. I can’t get my results
to match MATPOWER.
I
get Q into the banks fro
Russ Patterson wrote:
Hi friends,
Per Table 4-2 in the manual, I
want to increase the number of iterations to solve my power
flow. So, I added a line to my .m file as you see below.
But,
Russ Patterson wrote:
Hi
Carlos,
Thank
you. That makes sense. Based on that, I would expect to be
able to set the bus voltage, Vm, and it stick, but it
doesn’t. If I set Vm=
Russ Patterson wrote:
Hi
Uriel,
Is
there documentation (apologies if I missed it) that explain
how the power flow handles input data? For example, from
the case you ju
Rene Kruessmann wrote:
Hello everyone,
I have a question regarding how storages can get connected at different
positions within the grid. I want to add 2 storages, one at the top of the
grid and the other one at the bottom, that are connected with each other
to reduce the costs of redispatch. Th
Carlos Ferrandon Cervantes wrote:
Hello everyone:
I’ve been working quite a while with a UC problem
in MOST, using the IEEE RTS. I realised that when I compare the
mdo.UC.CommitSched results with mdo.results.ExpectedDispatch
the
Dear Yang:
If I understand correctly what you have in mind, then it seems to me
that you have two choices: 1) approximate the exponential in a limited
range with a higher order polynomial (which MATPOWER can easily handle),
or 2) write matlab code to implement user-defined cost functions (i.e
#FOO CEJUN, JOEL# wrote:
I am a student working on an OPF model for a 15 bus power
system.
I discovered an issue whereby the converged voltage magnitude
output values do not keep within the volta
Pfeifer, Lars Philipp wrote:
Hello everybody,
I have a question about the usage of the transformer in the branch matrix.
I have been reading the manual a few times and also took many looks in the
example cases.
But I still don’t really understand how I have to fill in the matrix for the
transf
.
Please use it, break it, and tell us about any issues that you find.
Sincerely,
Carlos E. Murillo-Sanchez, Universidad Nacional de Colombia - Manizales
Sergio García Marín, Universidad Nacional de Colombia - Manizales
Wilson González Vanegas, Universidad Tecnológica de Pereira
To expand a little bit on Ray's
response: In power systems, we tend to have labels for very
specific kinds of problems;
one such label is "power flow" or "load flow". This refers to a
problem in which given some input parameters (values
of loads, generator
Reinhold Bertram wrote:
Dear all,
is it possible to speed up the repeated solving of a deterministic
multi-period dispatch problem with MOST when only the loadprofile changes?
At the moment SetupTime is about 2 seconds for each iteration. Is it
possible to skip/shorten this time?
Best rega
Actually, since the quadratic costs
that Mariusz wants are diagonal, it
is possible to use the normal polynomial cost mechanism in
mpc.gencost (including the cost on the reactive generations; see
the footnote in the manual at page 152).
carl
Ismael K Abdulrahman wrote:
Dear Matpower community,
I have a system with one load bus faulted to ground (v=0). The system with
normal operation works fine. The system with the fault, however, cannot be
solved owing to low voltage problem. I tried to deal the faulted load bus as a
generator b
Jubeyer Rahman wrote:
Hi,
I have several questions regarding the AC power flow in matpower.
a. I would like to know whether the AC Power Flow command ('runpf')
,when called, does it have any constraint to constrain the transformer
power rating violation which is related to the real and reacti
amir ali Hosseini wrote:
Dear Sir,
i trying make a new network. 37 bus
but ! after runpf .i have warning!
can you help me ?
thanks a lot.
my warning is :
Warning: Matrix is singular to working precision.
> In newtonpf (line 89)
In runpf (line 204)
find attachment .
new network with 37 node.
aheen wrote:
No, I don't
On Sun, Feb 10, 2019, 1:24 AM Carlos E
Murillo-Sanchez <ce.murillosanc...@gmail.com
wrote:
Do you
include reactive generation in the decision
On Sat, 9 Feb 2019 at 15:42,
Carlos E Murillo-Sanchez <ce.murillosanc...@gmail.com>
wrote:
Rather,
it seems that your optimization method should take into
account any constraint viol
Rather, it seems that your optimization
method should take into account any constraint violation and fix
the issue by itself. In some methods, such as GA, one way to do
this is by using penalty functions on constraint violations.
carlos.
naime ahmadi wrote:
Hi,
I have question about branch data in case 57 bus system in matpower.
When I look at diagram of IEEE 57 bus system, I can see it have some
crossing in the network between lines.
Are there any crossing in power system?
Thanks,
Naime
Hi Naime: Short answer is yes. Longer
sidering
the ramping up/down limits of the generation?
Kind regards,
Carlos Ferrandon
El mar., 4 dic. 2018 a las 16:30, Carlos E
Murillo-Sanchez (<ce.murillosanc...@gmail.c
Carlos: MOST models post-contingency
states explicitly, and these states influence the possible base
case states directly so that the base cases are secure in light of
the specific contingencies considered.
carlos.
Carlos Ferrandon Cervantes w
Shiva Moshtagh wrote:
Hello
Can I use rectangular coordinates in matpower ? for example I know
that there are real and imaginary parts of power injections and power
flows as well as bus voltages, but for branch current there is not
rectangular coordinate. Suppose that I want to add derivativ
To clarify a little further on Ray's
advice:
An SVC can be approximated by a generator under voltage control
that can sink QMIN or source QMAX MVArs, but this would be
accurate for an SVC with a slope of infinity (i.e., constant
voltage). If you want
The equations are generally available
in the OPF literature. For specific GAMS modeling, you may want
to look at
https://www.gams.com/latest/psoptlib_ml/libhtml/index.html
carlos.
Jane Cheung wrote:
Dear all,
Actually, if it is just a linear cost
on the absolute value of the reactive output of a generator, you
can define a "V"-shaped piecewise linear cost for it using the
standard piecewise linear cost capability.
carlos.
Ray Zimmerman wrote:
e of the branch current. Thanks.
On Sat, Mar 31, 2018 at 12:41 PM,
Carlos E Murillo-Sanchez <ce.murillosanc...@gmail.com>
wrote:
Assuming
that the system data is in structure 'mpc', and that th
Carlos.
Mostafa Mohammadpourfard wrote:
Thank
you, Carlos. May I ask you to write the code here? I am not a
power system student. Thanks
On Sat, Mar 31, 2018 at 9:41 AM, Carlos
E Murillo-Sanchez <c
You have the bus voltages; now you just need to multiply the system's
nodal admittance matrix by the complex voltage vector to get the complex
currents vector. See the makeYbus function in MATTPOWER.
Carlos.
Arkan Arkan wrote:
Thank you, Ilias. But, unfortunately, there is not any direct
inf
The Ybus matrix computed from the data
in your file has NaN's and Inf's because branch # 69 from bus 16
to bus 17 has zero series impedance. You must collapse buses 16
and 17 into a single bus before applying any algorithm to the
system because the electrical "di
Let's see, Mohammed: I tried the data
in MATPOWER format that you posted originally, and applied the
case_info(0) function to it, yielding the following results:
>> case_info(mpc)
Checking connectivity ... 3 connected groups, 1 isolated bus
Elapsed
In the actual optimization problem, everything is in per unit.
Therefore, you probably want to specify the right hand side of your
inequality as 70/mpc.baseMVA .
Carlos.
Teiji Ponishi wrote:
Dear all,
I am trying to add the following simple constraints by direct
specification.
0 <= P
Matt: I have proof-of-concept code that implements line switching with
disjunctive constraints. Email me separately.
Carlos.
Matt Roveto wrote:
Hello,
I hope you are all well. I am trying to implement a transmission
switching DCOPF model to MATPOWER for use in a research project
according
slack bus
?
Thanks
On 4 Oct 2017, at 18:12, Carlos E
Murillo-Sanchez <ce.murillosanc...@gmail.com>
wrote:
Dear Charalambos:
Dear Charalambos:
To begin with, you have a load of 894.31 MW (computed with
sum(mpc.bus(:,PD)) and a total generation capacity of 30
(computed with sum(mpc.gen(:, PMAX)) . Your slack generator has
PMAX specified as zero. And it seems to me that many of
Stephen:
For hydro plants with reservoir, the way I have done this is by
assigning a daily (or weekly, depending on the length of the planning
horizon) amount of energy as the initial storage state, to be used
throughout the horizon optimally. As for the cost, it could even be
zero; the pr
Ashkan: the profit depends on the
commercial structure of the market and the type of auction
implemented, if it is an auction at all. There are many different
ways of choosing "payout functions" in an auction, and the profit
will depend on the chosen type.
Saeed: you need to add the sdp_pf
sub-folder to MATLAB's execution path.
carlos.
Saeed Ahmed wrote:
Hi
Community,
For
Making the he connection or bus incidence matrix i use
You can get an approximation to this on
the active power flow using the angular difference limits in
specified in the branch table. Columns ANGMAX and ANGMIN tell the
opf solver to implement
ANGMIN <= \theta_f - \theta_t <= ANGMAX
Carlos.
Oh, and by the way, do try IPOPT if you
still have not done so.
Carlos
Carlos E Murillo-Sanchez wrote:
Enrico:
I would certainly experiment with the values of the reactive
component of loads. In fact, an extreme
Best regards,
Enrico
2017-02-25 22:01 GMT+01:00 Carlos E
Murillo-Sanchez <ce.murillosanc...@gmail.com>:
Perhaps
you can try to run an AC OPF with branches that are
represented only by the seri
Perhaps you can try to run an AC OPF with branches that are represented
only by the series reactance (no line charging capacitance nor
resistance); one of the AC solvers might be able to solve the resulting
problem.
Carlos.
Enrico Vaccariello wrote:
Dear all MATPOWER developers and users,
I
Brenda: That is exactly what the slack generator is expected to do in a
traditional load flow. I suggest that you look up any classic textbook
(Bergen, Wood & Wollenberg, etc.) and learn about the load flow problem
and its standard assumptions.
carlos.
BRENDA ROJAS DELGADO wrote:
Hi, everyo
With the loading pattern in your
network, it is simply not possible to attain the voltage levels
that you want without using some sort of voltage compensation at
one or more locations. Even setting the voltage setpoint of the
source at 1.1 p.u., the resulting vol
Hi;
Looking at the diagram of the 118 bus system, it seems that the
branches in question are a combination of a transformer and a
line. Evidently, the pi model in the case data is an
approximation of this combined setup. But that is how the
origina
...@list.cornell.edu
[mailto:bounce-120671035-76197...@list.cornell.edu]
On Behalf Of Carlos E Murillo-Sanchez
Sent: Friday, August 05, 2016 11:10 AM
To: MATPOWER discussion forum
;
MATPOWER-L@cornell.edu
The LMP are shadow prices (Lagrange multipliers) on the active nodal
balance. Even if there is no load at a bus, its balance equation must
still be fulfilled and it carries with it a multiplier in the lagrangian
function.
Mounika Vanjarapu wrote:
I have a doubt regarding the LMPs in an IEEE-1
When you simplified away the PQ buses,
you lost the ability to specify loads in them. So, theoretically,
you could perform a load flow in which you can only have loads at
the 54 generator buses. You cannot use runpf, though, as runpf
starts by building the Ybus
we hope that it will
be useful for most people and most problems.
carlos.
Carlos E Murillo-Sanchez wrote:
Because modeling as a dispatchable
(up to Pmax) generator allows the optimization to curtail the
wind intake for security
From:
bounce-120538094-71172...@list.cornell.edu
[mailto:bounce-120538094-71172...@list.cornell.edu]
On Behalf Of Carlos E Murillo-Sanchez
Sent: Friday, June 03, 2016 11:29 AM
To: MATPOWER discussion
For each time period, there can be
several wind scenarios. The idea is to capture the variability of
wind, given a forecast, using those plausible wind realizations.
In each scenario, the Pmax value of the wind generator varies to
indicate the availability of win
Dear Mounika:
We need to understand the nature of your question. It seems as if you
might be trying to force a two-stage decision problem onto the framework
of a single-stage one. Hence Ray's answer. Please ellucidate.
carlos.
Mounika Vanjarapu wrote:
sir
can we pose penalties for th
Dear Victor:
Derivation of the lagrangian function with respect to Pg1 yields
the first order optimality condition
f'(Pg1) + mu_+ - mu_- - \lambda_P1 = 0
or \lambda_P1 = f'(Pg1) + mu_+ - mu_-
where \lambda_P1 is the multiplier
the 32
bit version with the command
: matlab -maci &
If you are running on windows, you might have the 32 bit version
somewhere in the Program Files-x86 folder.
Carlos.
Malcolm Barnacle wrote:
Dear Carlos E. Murillo-Sanchez,
I am trying to use the bpmpd solver for a dcopf problem. I
Hi;
You have used the TAP column of the branch table to indicate two
transformers with a 55:1 ratio, and then you used the LV
impedances, I think, still in PU with a base voltage of 22KV; they
are pretty large. The TAP column should be used to specify smal
Dear Dominic:
The PF solver takes the generator voltage setpoint from the gen(:,
VG), thus changing the voltage values set by the OPF at generator
buses. Try
resultsOPF=runopf(mpc)
%calculate opf
We use
gP(Theta,V,Pg) = gp1(Theta,V) + Pd - Pg = 0
One advantage of writing the mismatches like this is that the
lagrange multipliers come out naturally positive out of the OPF.
Abhyankar, Shrirang G. wrote:
It
Actually, the capacitive susceptance
seems to be way too high (at least an order of magnitude, perhaps
even more considering it is an MV network) for the given R and X.
Carlos.
Marco Barbetta wrote:
Yes, i report a little part of mpc.
2014-11-28 20:34 GMT+01:00 Carlos E
Murillo-Sanchez <ce.murillosanc...@gmail.com>:
If there
is only one point of connection from the MV network to the
HV one, why can't you model this particular bus as a
reference bus wit
If there is only one point of connection from the MV network to the HV
one, why can't you model this particular bus as a reference bus with a
generator that can have both positive and negative Pg? Just set Pmin to
something negative and Pmax to some positive value.
Carlos.
Marco Barbetta wro
Isabel:
Try creating a small band of acceptable voltage limits at the
generator bus
say VMIN = Vcenter - epsilon <= V <= Vcenter + epsilon =
VMAX, then start slowly
moving Vcenter from the initial optimization problem's solution
towards th
4...@list.cornell.edu
[bounce-117292143-69334...@list.cornell.edu]" im Auftrag von "Carlos
E Murillo-Sanchez [ce.murillosanc...@gmail.com]
Gesendet: Montag, 14. Juli 2014 19:47
An: MATPOWER discussion forum
Betreff: Re: AW: Convergence issues: OPF simulation of several interconnect
Tell me a little more about how you stitched together the different
systems. Did you have separate data files for different sections of
the system and then you actually interconnected them? What do you mean
by "realistic" at this stage? Or did you have one single file with all
sections in i
Are there any phase shifters in the system?
Eser Patrick wrote:
Dear MATPOWER Community,
I am trying to solve a relatively large (>1000 buses) AC OPF problem. It
consists of several individual countries (200-300 buses each). Please allow me to
describe my procedures, to paint the full picture
[ref pv, pq] = bustypes(bus, gen) ( or [ref pv, pq] =
bustypes(results.bus, results.gen) )
and then
results.bus(pq, VM)
carlos.
Kusi, Samuel A wrote:
Is there any way in MatPower to call out just the load voltages?
results.bus(:, VM) calls all buss voltages. Is there any code that can
Branch admittances are additive in the
Ybus matrix construction procedure. makeYbus uses the sparse()
function to build Ybus, and automatically adds any admittances
resulting from parallel branches in the right places. You should
have no problem with parallel li
It seems to me that you have an
instance of a data file which is aimed at detailed modeling of all
of the system's elements. This kind of data is useful for
graphical representation, but it poses a major problem: it
assigns different bus numbers to elements that
No formula is used; it comes directly
from the optimizer algorithm, which is trying to solve the first
order optimality conditions of the optimization problem. These
conditions include the Lagrange multipliers for the power balance
equations and that's where the
Dear Jiashen:
The model used in Matpower is the standard pi one-line equivalent of a
three phase transmission line, with quantities in per-unit. This is
constructed from a base voltage equal to sqrt(3) times the phase voltage
and a base power equal to three times the per-phase power. With thi
One possible way, using system replication to build an augmented system
with coupling constraints among the different islands, is what we
originally employed a few years ago for this:
http://www.sciencedirect.com/science/article/pii/S0167923613001152
linhomeperson wrote:
> Dear MATPOWER users,
>
Mr. Acker:
The mathematical problem that is solved behind the scenes is simply
singular in light of such modeling detail. Any robust solver would
actually look for such occurences in the data and collapse the busbars
into a single one, and then, as a post-solution-processing strategy,
split
As well as adding coordination costs in
separation and coordination optimization schemes :-) .
Ray Zimmerman wrote:
The basekV column of the bus matrix has no effect on the
solution.
I am not aware of any real cost ass
Unfortunately, this problem is of an
intrinsic nature. If it is an OPF that you are running, you might
try using additional linear constraints to equate the two bus bars
magnitudes and angles, and hope that the underlying linear algebra
algorithms will deal with
Ahmad:
If baseMVA is 100, then Bs in pu should be 0.19 .
carlos.
ahmad rezaee wrote:
Dear
Dr Zimmerman and all friends
I
have written one load flow code with Gauss-Sidel and one with Newton-Raphson.
Both of them lead to the same results. For validating my codes I compared the
results with
Santiago:
You can use a piecewise linear cost model with two segments, with vertex
points
at (-Qnegmax, Qnegprice*Qnegmax), (0,0), (Qposmax, Qposprice*Qposmax) .
Santiago Torres wrote:
Dear Ray and Matpower friends
I want to use the reactive power cost of AC optimal power flow. I
understa
You can specify bus shunt conductances and capacitive susceptances in
the GS and BS columns of the bus matrix as MW and Mvars at nominal
voltage. And yes, the assumption that loads are modeled as PQ quantities
is standard in the load flow problem. Using a different model for loads
is a non-stan
Dear Bela:
Running MINOPF on the system, it seems that both active and reactive
balances are hard to fulfill in bus 706.
Carlos.
Bela Genge wrote:
Dear all,
I am having some problems running an AC power flow with the attached
case system. This is an approximate EU Power Grid model that was
Richard:
One easy way to minimize losses using the opf solver is to assign to all
generators the same linear "fuel" cost. If all sources have the same
cost, then the optimizer will try to minimize losses. However, you have
to provide sensible VMIN, VMAX limits for all buses in the bus table,
I think that this is actually a good review of the issues involved in
attempting a true probabilistic power flow:
http://ieeexplore.ieee.org/xpl/login.jsp?tp=&arnumber=4523658&url=http%3A%2F%2Fieeexplore.ieee.org%2Fiel5%2F4511470%2F4523365%2F04523658.pdf%3Farnumber%3D4523658
iman wrote:
Dear A
Iman: two methods are compared in this paper. The 2-point estimate
tries to capture moments of the distribution using just two
_deterministic_ OPFs. The cumulant method assumes linear propagation
(i.e., based on a sensitivity analysis of the KKT conditions) of the
pdfs of the quantities deem
Hi;
Unfortunately, the original data does not include those
parameters.
A while ago, based on information about base transmission line
voltage, standard conductor sizes, and information about whether
certain lines were single or double circuit,
Richard:
Are you sure that you want to specify the power injections of the
photovoltaic generators as zero in a power flow?
The trapezoidal model for (P,Q) feasibility is intended for use in an
optimal power flow setting, not a power flow setting. If you have all
generators with Pg = 0 then
Hi;
There are diagrams for the ieee 30 bus system and the 118 bus system
around, may also for others.
I will send you these two in a separate mail.
carlos.
Shri wrote:
On Aug 20, 2012, at 11:11 AM, Carlos E Murillo-Sanchez
wrote:
Power networks are usually almost planar, with a few
my question is what is MATPOWER doing for solving the DC power flow
that causes this.
It would be of great help to me if you answer this.
All the best,
Hanie
On Tue, Aug 21, 2012 at 2:46 PM, Carlos E Murillo-Sanchez
mailto:ce.murillosanc...@gmail.com>> wrote:
Let me see if I ge
Let me see if I get this rightt:
1) You define real power dispatches of non-slack generators in a test
case, perhaps also change the amount of load, and then
2) you run a power flow, ostensibly with runpf. But then,
3) you find that the dispatch values that you had defined in step (1)
for gene
Hi;
I am not sure that I completely follow your second question. With
regards to the first one, we really have just the ones that are
distributed with MATPOWER, plus some others that we use in research and
which may contain propietary data.
Your second question... if the solvers find a sol
1 - 100 of 127 matches
Mail list logo
