Re: objective function for renewable energy generation

I suppose it depends on the market rules. But a cost of $0/MW is one option that is fairly typical. In some cases, there are subsidies involved which would give the renewables an incentive not to curtail even if the price goes a bit negative, in which case you could even use a negative cost.

Re: UPFC modelling in MATPOWER

Unfortunately, there is no such guide provided by MATPOWER. Modeling of various FACTS devices is something planned for version 9. If you are interested in building your own model, you may find How to Create a New Element Type and the next two how to

Re: Inquiry about Bus Voltage Variation in MOST Examples

resolution? Best, Lichen From: bounce-128084230-96962...@list.cornell.edu<mailto:bounce-128084230-96962...@list.cornell.edu> mailto:bounce-128084230-96962...@list.cornell.edu>> on behalf of Ray Daniel Zimmerman mailto:r...@cornell.edu>> Sent:

Re: Inquiry about Bus Voltage Variation in MOST Examples

Hello Lichen, I’m afraid MOST uses a DC network model. That is, bus voltage magnitudes are not modeled at all and instead are assumed to be approximated by 1 p.u. See Section 3.7 in the MATPOWER User’s Manual for more details on the DC

Re: Distributed Generators Model

(column 4) of branch data? Is this the correct way? On Thu, Nov 23, 2023 at 5:31 PM Muhammad Junaid mailto:gaylethunder...@gmail.com>> wrote: Thank you so much for the detailed answer. On Wed, Nov 22, 2023 at 10:11 PM Ray Daniel Zimmerman mailto:r...@cornell.edu>> wrote: Adding them

Re: Proposed 533-bus system based on real data

edu> mailto:bounce-127939752-95792...@list.cornell.edu>> För Ray Daniel Zimmerman Skickat: den 12 december 2023 19:10 Till: MATPOWER-DEV-L mailto:matpower-de...@list.cornell.edu>> Ämne: Re: Proposed 533-bus system based on real data Hi Gabriel, I prefer to avoid anything that require

Re: Proposed 533-bus system based on real data

beneficial to perform similar calculations for the actual cables in the test system to validate if corrections to the rated current are necessary. Best regards, Mirko On 12/6/23 18:28, Ray Daniel Zimmerman wrote: Hello Gabriel, I’d be happy to include this in the next MATPOWER release. I’ve attache

Re: Proposed 533-bus system based on real data

icity market environment. sir, If you have any basic demo source file, please send, that could be help in my research. for learning purpose please respond sir On Wed, Dec 6, 2023 at 10:58 PM Ray Daniel Zimmerman mailto:r...@cornell.edu>> wrote: Hello Gabriel, I’d be happy to include this in

Re: Proposed 533-bus system based on real data

, please go ahead and create a pull request on GitHub to add the file. Thanks again for this contribution! Ray On Nov 30, 2023, at 10:38 AM, Ray Daniel Zimmerman wrote: Absolutely, we would be interested in including this. Let me have a look at it and get back to you about any suggested

Re: Power Flow Calculation Error

tp://radial.pf/> separately? On Wed, Nov 22, 2023 at 10:16 PM Ray Daniel Zimmerman mailto:r...@cornell.edu>> wrote: What you show is a warning (not an error) coming from the Power Summation method. As long is the method converges with success = 1, the result should be valid. Also, if

Re: Proposed 533-bus system based on real data

Absolutely, we would be interested in including this. Let me have a look at it and get back to you about any suggested modifications. Ray On Nov 29, 2023, at 9:13 AM, Gabriel Malmer wrote: Dear Matpower-developers, My name is Gabriel Malmer and I’m a PhD student in power systems at Lund

Re: Overload Condition

You can certainly write a little function to return the quantities (power flows, currents, voltages) of interest based on the results of a PF or OPF, and compare these as you add generation. If the size of the matrices are not changing you can also use

Re: Power Flow Calculation Error

What you show is a warning (not an error) coming from the Power Summation method. As long is the method converges with success = 1, the result should be valid. Also, if the Newton power flow converges, its result should also be fine. Sometimes the Newton power flow has more difficulty on

Re: Power Flow Calculation Error

I didn’t write the code for the power summation method, so I wasn’t aware that it did not handle multiple generators at a single bus. I consider this a bug. Would you mind creating a bug report for it here? Thanks, Ray On Nov 18, 2023, at 6:37

Re: Distributed Generators Model

Adding them to columns 3 and 4 of the bus matrix will be equivalent to adding them to columns 2 and 3 of the gen matrix, assuming the bus is set as a PQ bus and you are doing a simple power flow. Adding it to bus 5 and 6 of the bus matrix is not equivalent, since the power (as a constant

Re: Finding the maximum positive (and negative) power injection at a bus

it be used as negative generator? From: Ray Daniel Zimmerman Sent: Wednesday, 8 November 2023 9:14 PM To: MATPOWER-L Reply To: MATPOWER discussion forum Subject: Re: Finding the maximum positive (and negative) power injection at a bus See Section 6.4.2 on Dispatchable Loads in the MATPOWER User’s

Re: Piecewise Linear Generator Cost Function

1. Yes. 2. Yes. It is simply a sum of the y variables, but the y variables are constrained to lie above c(x). 3. Yes. On Nov 10, 2023, at 4:32 AM, Ronald Cabaoig wrote: Dear MATPOWER team, I'd like to clarify some things about the piecewise linear generator cost function. Question 1:

Re: I need help regarding MATPOWER implementation for my research

Hello Naresh, FYI, the MATPOWER-DEV-L is for discussion of the development of MATPOWER itself, not really for discussing implementation ideas for a particular application of MATPOWER, such as you describe. With regards to the use of MATPOWER for your problem, I’m afraid that’s beyond the

Re: Convertion of CIM CGMES format into a MATPOWER case struct

Unfortunately, I am not aware of any such converter and I am not in a position to create one. However, if someone else is, I’d be happy to include it in the MATPOWER distribution. Ray On Nov 9, 2023, at 1:43 AM, oliver.zehn...@elcom.admin.ch wrote: Dear all, ENTSO-E, the European Network

Re: How to run single side auction in MATPOWER

If you use all fixed loads (not dispatchable loads), it will be a single-sided market with just offers. If you use negative fixed loads to represent the generation, and make the loads dispatchable, you will have a single-sided market with just bids. Ray On Nov 8, 2023, at 4:33 AM, Dr. D.

Re: Finding the maximum positive (and negative) power injection at a bus

cenario 2: minimize the negative power injection in bus X, with a power factor constraint of 0.80 leading pf. Cheers, Ronald On Sat, Oct 14, 2023 at 6:29 AM Ronald Cabaoig mailto:rrcaba...@up.edu.ph>> wrote: Thank you very much! Cheers, Ronald On Sat, Oct 14, 2023 at 1:33 AM Ray Daniel Zimm

Re: Question on running MATPOWER runoff in LabVIEW mathscript

any suggestion. Best regards, Kevin From: bounce-127852508-95296...@list.cornell.edu<mailto:bounce-127852508-95296...@list.cornell.edu> mailto:bounce-127852508-95296...@list.cornell.edu>> on behalf of Ray Daniel Zimmerman mailto:r...@cornell.edu>> Date: Saturday, 28 Octob

Re: Question on running MATPOWER runoff in LabVIEW mathscript

I'm afraid I do not have experience with LabVIEW's MathScript, but if it can be made to run a MATLAB script, then you can certainly make it call MATPOWER's runpf(). You may have to write your own wrapper script to handle the inputs from LabVIEW and set up the call to runpf() and then package

Re: What the objective function denotes in the smart market

edu> mailto:bounce-127850109-88225...@list.cornell.edu>> on behalf of Ray Daniel Zimmerman mailto:r...@cornell.edu>> Sent: Thursday, October 26, 2023 11:31 PM To: MATPOWER-L mailto:matpowe...@list.cornell.edu>> Subject: Re: What the objective function denotes in the smart market I see wh

Re: What the objective function denotes in the smart market

127836391-88225...@list.cornell.edu> mailto:bounce-127836391-88225...@list.cornell.edu>> on behalf of Ray Daniel Zimmerman mailto:r...@cornell.edu>> Sent: Friday, October 20, 2023 2:53 AM To: MATPOWER-L mailto:matpowe...@list.cornell.edu>> Subject: Re: What the objecti

Re: What the objective function denotes in the smart market

2985.78 represents. | From: bounce-127828667-88225...@list.cornell.edu<mailto:bounce-127828667-88225...@list.cornell.edu> mailto:bounce-127828667-88225...@list.cornell.edu>> on behalf of Ray Daniel Zimmerman mailto:r...@cornell.edu>> Sen

Re: What the objective function denotes in the smart market

The smartmarket code simply sets up the cost functions and generator limits for the OPF problem. The OPF minimizes total generation cost. In the case that includes dispatchable loads, modeled as negative generation, this cost is essential the cost of the normal generators plus the negative cost

Re: Finding the maximum positive (and negative) power injection at a bus

+ 0Pg2 + ...) 4. Retain all other information Thanks in advance! Cheers, Ronald On Tue, Oct 10, 2023 at 5:23 AM Ray Daniel Zimmerman mailto:r...@cornell.edu>> wrote: You can do this by placing a dummy generator at the bus of interest and using a very large negative (to maximize t

Re: Finding the maximum positive (and negative) power injection at a bus

You can do this by placing a dummy generator at the bus of interest and using a very large negative (to maximize the injection) or positive (to minimize the injection) generation cost. You can either zero out all other costs, or simply make the magnitude of the cost large enough to dominate all

Re: Plot SE gain matrix

That is correct … PV buses come first then PQ buses. Ray > On Sep 28, 2023, at 1:33 PM, Themistoklis Xygkis > wrote: > > Hello collegues. > > I have a small question regarding the SE function; > in doSE, J, that is, the gain matrix of the problem is calculated. > The diagonal entries of

Re: MOST Optimization Issue with PV Generation Less Than Half of Load

This mailing list (or a GitHub issue) is still my preferred way to respond to these sorts of questions, but for some reason I saw and responded first here this time. Ray On Sep 25,

Re: loss optimization. Case with pv bus.

Could you please clarify? I’m afraid I don’t follow you. For case118, as with any other case, you can minimize losses by setting the costs of all generators to an equal positive cost (e.g. $1/MW) and running an OPF. The concept of a PV bus is not relevant for the OPF, only for power flow.

Re: runopf. how do I find the constraint. which removes cost calculation

en out limit check_index=3 - VG out limit check_index=4 -. .. 22 вересня 2023, 22:02:25, від "Ray Daniel Zimmerman" mailto:r...@cornell.edu>>: I assume you are asking whether it’s possible to have checklimits() return a 0 or

Re: Understanding the Profile Struct in MOST

For future reference, it looks like I saw and answered it here first. Ray On Sep 22, 2023, at 8:29 AM, keir steegstra wrote: Hello everyone, I'm working on a project where I'm

Re: Lagrangian multiplier and kuhn tucker multiplier for loss reduction

Your constant and equal generator cost is the conversion factor you need to convert between changes in cost cost and changes in loss, which are equal to changes in total generation. So if all of your generator costs are $1, then shadow prices on inequality constraints are equal to the

Re: runopf. how do I find the constraint. which removes cost calculation

I assume you are asking whether it’s possible to have checklimits() return a 0 or 1 depending on whether or not there are limit violations. I don’t really have the time to work on this now, but I shouldn’t be too hard to add. You’d simply check to see if all of the following outputs are empty

Re: runopf. how do I find the constraint. which removes cost calculation

Hello Igor, First, by default, runopf() does not consider the VG column of the gen matrix, but rather optimizes the generator voltage setpoint, keeping it between the VMIN and VMAX for the corresponding bus. So changing VG in the input will have no effect. You can set the opf.use_vg option to

Re: Battery storages

Hi Heller, MATPOWER’s OPF does not have an explicit storage model. A battery is by nature a device that links operation across time, whereas an OPF is a single instant snapshot. So for the purposes of an OPF, you have a couple of options. * Determine the storage dispatch by some process

Re: Jacobian Matrix

(:, [PD QD])*1.0; mpci = ext2int(mpct); Jacc = makeJac(mpci); Jac=full(Jacc); Jac_inv= inv(Jac);% inverse of J sen_vpp=Jac_inv(1:117,118:181); % sensitivity of voltages to real power injection On Wednesday, September 6, 2023 at 10:46:30 AM EDT, Ray Daniel Zimmerman wrote: Can you describe

Re: fix the generator capacities in runpf

Remember the definition of bus types for the power flow problem. PV means that active power and voltage magnitude are specified and we are solving for voltage angle and reactive power. PQ means that active and reactive power are specified, and we are solving for voltage magnitude and angle.

Re: Jacobian Matrix

Can you describe exactly what you are doing in a bit more detail (or maybe provide a few lines of code to illustrate)? I’m not sure where you are getting a matrix that is 117x64. Ray On Sep 3, 2023, at 7:10 AM, Kaled Zareer wrote: Dear MATPOWER team, I used inv(Jac) to obtain the

Re: OPF with Multiple Objectives

I’m pretty confident that the answer to the question you are asking is “no”. But I’m still curious what you mean by “run the OPF with two or three objectives”. What would characterize the kind of solution you are looking for if it is not optimizing a weighted sum of the multiple objectives?

Re: OPF with Distribution Network

It will only be useful if you have some dispatchable elements (generation or loads), but it is not fundamentally different. Whether or not the load flow converges will depend on the parameters of the system, the algorithm used, the starting point, etc. Ray > On Aug 25, 2023, at 8:46 AM,

Re: MATPOWER 8.0b1 Installation problem

See https://github.com/MATPOWER/matpower/issues/201 Ray On Aug 18, 2023, at 4:55 AM, Jessica Korkmaz wrote: Hello, I am using MATLAB R2023a, and I am installing matpower 8.0b1. After running 'install_matpower.m' and running the required test, I got the following error message: Would you

Re: Dispatchable component that acts as a load or a generator

t;generator" mode, P is injected) or a negative ("load" mode, P is drawn) value depending on the system scenario? Cheers, Ronald On Tue, Jul 18, 2023 at 11:52 PM Ray Daniel Zimmerman mailto:r...@cornell.edu>> wrote: Yes. A generator is simply a complex power injection and

Re: Dispatchable component that acts as a load or a generator

Yes. A generator is simply a complex power injection and its active and reactive limits can be anything. Obviously, feasibility requires that PMIN <= PMAX and QMIN <= QMAX. The only special case is that of a dispatchable load, defined as a generator with PMAX=0 and PMIN < 0, in which case an

Re: derive intermediate steps of each variable

That depends on the solver being used. The default MIPS solver does not output intermediate results, but I suppose you create a modified version of mips.m where you collect the variable x at each iteration and add it to the hist struct that gets returned in results.raw.output.hist. Ray

Re: voltage difference of generators runpf runopf case9

The generator bus voltage for the power flow is set by the VG column of the gen matrix. For the OPF it is an optimization variable must remain between the bounds set by the VMIN and VMAX columns of the bus matrix. You can also fix the voltages at the VG setpoint in the OPF by setting the

Re: On using Piecewise or Nonlinear Cost Functions

There is not an easy way to do piecewise non-linear cost functions for a single generator. The easiest approach is probably to split the single generator into two, each with a quadratic cost. If this is not practical for some reason, you’ll probably have to do it via an extension using

Re: Maximize the distance from several feasible points in AC OPF problem

ror, I found that the gradient DF should be a column vector, like the requirement of MIPS. Would you please kindly clarify this point for me? Thank you. At 2023-04-08 01:09:29, "Ray Daniel Zimmerman" wrote: Another great question. The way to do it is to define a function handle that

Re: Maximize the distance from several feasible points in AC OPF problem

variable. Would you please kindly advise me how can I use args inside fcn more elegantly? Thank you very much. Best, Yuanxi Wu At 2023-04-07 23:22:44, "Ray Daniel Zimmerman" wrote: 1) Yes, that is correct. 2) Yes, the individual cost terms are simply summed together. And

Re: Maximize the distance from several feasible points in AC OPF problem

1) Yes, that is correct. 2) Yes, the individual cost terms are simply summed together. And I believe a negative definite Q should be fine. 3) This is a much harder question. I’m afraid I don’t have a definitive answer to why the matrix for the Newton step in the MIPS primal-dual interior point

Re: How to rename variable set names of optimization vector in AC OPF

At 2023-04-04 03:04:06, "Ray Daniel Zimmerman" wrote: Good questions … 1) It’s the same as the order in the mpc.bus matrix. 2) You cannot change the names of the existing variable sets, but in the functions where you calculate the constraints, gradients and hessians, you can

Re: How to rename variable set names of optimization vector in AC OPF

Good questions … 1) It’s the same as the order in the mpc.bus matrix. 2) You cannot change the names of the existing variable sets, but in the functions where you calculate the constraints, gradients and hessians, you can pull out a subset of the full vector for your calculations. For this you

Re: FMINCON - constraint violations

The easiest solution is probably to enable soft constraints using toggle_softlims(). See Section 7.6.4 in the MATPOWER User’s Manual for more info. Ray On Mar 29,

Re: runpf and generation limits constraints

By default the runpf() command does not enforce any of the limits you mention. The purpose of power flow is to compute the unknown voltages and flows for a set of specified generator voltages and active power injections. There is an option (pf.enforce_q_lims, see Table 4-2 in the MATPOWER

Re: Matpower runopf cost function

Hello Muhammad, Table B-4 in the MATPOWER User’s Manual describes the gencost matrix which provides the generator costs that are minimized by runopf() by default. Specifically, as it says in the footnote … ... the first ng rows of gencost

Re: Using Matpower to find capability curves of distribution grid using OPF

Sorry, I was dealing with e-mail issues that caused me to miss messages to this list for a week or so. See my response here. Ray On Mar 6, 2023, at 2:49 AM, Akhtar Hussain Javed wrote: Dear All, I am trying to find the reactive power

Re: Adding PV to bus 5

Did you add a row to gencost when you added a row to gen? If they don’t have the same number of rows, that’s what’s causing the error. Ray On Feb 26, 2023, at 10:05 PM, naime ahmadi wrote: Dear Mr. Hosseini, I appreciate your guidance on changing the bus type of bus 5 from "PQ" to "PV".

Re: Using matpower to draw PV Curves.

What’s described in the current User’s Manual would be the way to go. The other option appears to be from MATPOWER 4.0 (from 2011), which had continuation power flow included as contributed code in extras. It was cleaned up and moved into the main MATPOWER distribution for version 5 (2014),

Re: Regarding storing of the result from optimal power flow

See Section 6.6 in the MATPOWER User’s Manual, especially Table 6-1. There are a few things, like losses, that are displayed in the console output, that are not included in the results struct. These can be computed using get_losses(). Ray

Re: Regarding to the ill-condition problem of Newton's method

I’m afraid I don’t have anything very concrete to offer in response to your questions … but in my experience also I’ve been surprised at how well the primal-dual interior point method converges even in the presence of fairly ill-conditioned matrices. My only suggestion in terms of approaches

Re: MIPS - Matpower Numerical error for arbitrary initial values

Hi Àlex, Sorry for the delay in responding. Unfortunately, I’m not sure what to suggest. All I can tell you is that the alphap is used to scale the size of the Newton step in order to prevent the slack variable Z from becoming negative (see Section 3.3.4 in the MIPS User’s

Re: Exclude control variables

Simply constrain the variables that you want to “exclude” by setting the upper and lower limits to the desired value. Ray On Dec 23, 2022, at 4:42 AM, karthikaikan...@eee.sastra.edu wrote: Dear sir, In opf, I want to exclude all control variables

Re: Best Sitting and Sizing of Distributed Generators on MOST

What is the error message you are getting? Is it coming from MATPOWER, or from your code? Ray > On Dec 7, 2022, at 7:01 AM, feyyaz fatih aydın wrote: > > Hi, > > I am trying to apply Best Sitting and Sizing of Distrubuted Generators > problem to find the minimum LMP by using PSO. > I

Re: Contingency Result

The probabilities are used to weight the costs of the scenarios, but the corresponding contingency scenario will always reflect the specified outages, regardless of the probabilities assigned. Ray On Dec 5, 2022, at 8:00 AM, Yusuf Susilo Wijoyo mailto:yusu...@ugm.ac.id>> wrote: Dear All,

Re: Adding wind turbine and getting results based on transmat

I assume you are trying to set up a single-period problem in MOST, given that you mention a transmat input that appears to be a single row vector. In that case, results of each scenario j should be available in mdo.flow(1, j, 1).mpc. Be careful how you interpret the LMP values in the LAM_P

Re: Can the runopf control the switched shunts?

MATPOWER does not currently optimize switched shunt settings. The closest you could get would be to treat it as a (continuously) dispatchable reactive power source, that is a generator with PMIN=PMAX=0. Ray On Nov 18, 2022, at 8:44 AM, Jianqiao Huang mailto:jhuan...@hawk.iit.edu>>

Re: Fixed Reserve and Contingency on Stocastic DC OPF with MOST

As described in Section 4.2.2 of the MOST User’s Manual, the zonal reserves and the contingency constraints are considered as 2 alternative ways to handle security. MOST was not designed (or tested) for applying both simultaneously … so I’m not

Re: Minimize the curtailment of renewable energy production

You can do it simply by setting the gencost appropriately. For example, you can use a negative marginal cost on your renewable units. For example, a cost of -$10/MW means the unit would need to be paid $10 for every megawatt of curtailment. Ray On Nov 20, 2022, at 10:00 PM, seuyxw

Re: Probabilistic AC-OPF with Matpower

I believe Carlos Murillo and some of his students have experience with what is essentially an AC OPF version of MOST, but I’m not of the status of that work. Ray On Nov 18, 2022, at 8:24 AM, Santiago Torres mailto:santiago.i...@gmail.com>> wrote: Hi Leopold, maybe you can do it using a

Re: Gurobi and MOST Interaction

MATPOWER checks for the presence of solvers, like Gurobi, only once and then caches the availability information. So after you install Gurobi, you have to either restart MATLAB or type: have_feature('all', 'clear_cache') to clear the cached information so it will detect the presence of Gurobi

Re: Query in power flow analysis of radial distribution system with two voltage levels

If you convert everything into per-unit correctly, using a different voltage base for the different buses, it shouldn’t be a problem. I suggest finding a good textbook on converting parameters to per-unit to make sure you get it right, especially for the transformer. Ray > On Nov 7,

MATPOWER manual

ist.cornell.edu> mailto:bounce-126942489-88411...@list.cornell.edu>> On Behalf Of Ray Daniel Zimmerman Sent: Wednesday, November 2, 2022 9:02 AM To: MATPOWER-L mailto:matpowe...@list.cornell.edu>> Subject: FEEDBACK REQUEST: Octave requirements for new MATPOWER Hi MATPOWER Users, I need you

Re: FEEDBACK REQUEST: Octave requirements for new MATPOWER

y savy enough to run simultaneous versions without trouble. Best regards, russ From: bounce-126942489-88411...@list.cornell.edu<mailto:bounce-126942489-88411...@list.cornell.edu> mailto:bounce-126942489-88411...@list.cornell.edu>> On Behalf Of Ray Daniel Zimmerman Sent: Wednesday, Nov

Re: nominal value

Without more context, I’m afraid I’m not sure what you mean by “nominal value” and “power injection” in this case? Injection per generator? per load? per bus? Ray > On Nov 2, 2022, at 10:31 PM, naime ahmadi wrote: > > Hello, > How can I calculate the nominal value of active and reactive

Re: FEEDBACK REQUEST: Octave requirements for new MATPOWER

option to switch between HTML and PDF or remove the links altogether. Richard On Wed, 2 Nov 2022 at 16:00, Ray Daniel Zimmerman mailto:r...@cornell.edu>> wrote: Thanks, Richard. Btw, I’m still impressed with https://matpower.app<https://matpower.app/> and the technology that make it possib

Re: Adding DG and Multiobjective For Power Transmission Systems

You can model a DG as you would any other generator, just make sure the capacities and costs reflect their actual capabilities. In terms of optimizing energy cost and emissions, simply combine the two into a single cost that prioritizes / weights them according to your requirements, then use

Re: FEEDBACK REQUEST: Octave requirements for new MATPOWER

something related to OpenMP or message passing or multiple precision arithmetic. Richard On Mon, 31 Oct 2022 at 23:59, Ray Daniel Zimmerman mailto:r...@cornell.edu>> wrote: Hi MATPOWER Developers, I need your feedback on a quick question. I’m working on finalizing things for a beta

Re: FEEDBACK REQUEST: Octave requirements for new MATPOWER

something related to OpenMP or message passing or multiple precision arithmetic. Richard On Mon, 31 Oct 2022 at 23:59, Ray Daniel Zimmerman mailto:r...@cornell.edu>> wrote: Hi MATPOWER Developers, I need your feedback on a quick question. I’m working on finalizing things for a beta

FEEDBACK REQUEST: Octave requirements for new MATPOWER

Hi MATPOWER Users, I need your feedback on a quick question. I’m working on finalizing things for a beta release of what amounts to a nearly complete re-write of MATPOWER for version 8.0. More on that soon. Since this new version defines tons of new classes, I thought it would be nice to put

FEEDBACK REQUEST: Octave requirements for new MATPOWER

Hi MATPOWER Developers, I need your feedback on a quick question. I’m working on finalizing things for a beta release of what amounts to a nearly complete re-write of MATPOWER for version 8.0. More on that soon. Since this new version defines tons of new classes, I thought it would be nice to

Re: Load Flexibility

e load Best regards, Lois On Wed, Oct 19, 2022 at 3:09 PM Ray Daniel Zimmerman mailto:r...@cornell.edu>> wrote: When using dispatchable loads, the range of load variation is determined by the limits on the associated “generator”. So for a load that can be dispatched between 50 and 100 M

Re: error in runpf function

I am not able to reproduce that error. That command works fine for me on MATPOWER 6.0. I suspect there may be something wrong with your MATPOWER installation. I’m guessing either one of the MATPOWER files has been inadvertently modified from the original, or MATPOWER has not been properly

Re: Load Flexibility

t does not increase. I tried using the if/else rule but it didn't work. Please how can I do this, Also, how can I add constraints in my OPF code, say I want to ensure that the load reduction doesn't exceed a particular value? Best Regards, Lois. On Sat, Sep 10, 2022 at 1:35 AM Ray Daniel Zimmerman ma

Re: Ybus

0.05*Q_bus(i,:); Q_inj_min(i,:) = Q_bus(i,:)+0.05*Q_bus(i,:); end end On Thu, Oct 13, 2022 at 1:33 AM Ray Daniel Zimmerman mailto:r...@cornell.edu>> wrote: I’m afraid I don’t understand what you are trying to calculate. Ray On Oct 11, 2022, at 11:22 AM, naime ahmadi mailto:nai

Re: Ybus

_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 follows (e.g. for case9): define_constants mpc = loadcase('case9'); [Ybus, Yf, Yt] = makeYbus(mpc); V = mpc.bus(:, VM) .* exp(

Re: Ybus

Good catch! On Oct 10, 2022, at 5:13 PM, Carlos E Murillo-Sanchez mailto:ce.murillosanc...@gmail.com>> wrote: 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

Re: Ybus

In matrix form, you can do the corresponding calculations as follows (e.g. for case9): define_constants mpc = loadcase('case9'); [Ybus, Yf, Yt] = makeYbus(mpc); V = mpc.bus(:, VM) .* exp(1j * mpc.bus(:, VA)*pi/180); Sbus = V .* conj(Ybus * V); Sf = V .* conj(Yf * V); St = V .* conj(Yt * V);

Re: Output

mpopt = mpoption('out.all', 0); Ray On Oct 5, 2022, at 2:45 PM, Scott seedell mailto:scottseed...@hotmail.com>> wrote: Hi, What is the best method to suppress the printing on the screen after each power flow? I’m running a monte carlo simulation with a large number of cycles and would be

Re: Base impedance inload flow

All p.u. parameters (e.g. branch impedances) are given with respect to a particular per-unit power base and should be specified (or converted) to match a common, system-wide base. If you do that, as long as it used consistently for all p.u. parameters, you should get the same solution. Ray

Re: Question about OPF problem

This is good question. It turns out that, in general, increasing the voltage tends to reduce losses, which in turn reduces total generation and generation cost. So minimizing generation cost also tends to dispatch reactive power so as to push voltages upward. Fiddling with the OPF starting

Re: Question about fmincopf_solver in MATPOWER

The fmincopf_solver() function is not a top-level user function. It is called indirectly by runopf() when you have the opf.ac.solver option set to 'FMINCON'. >> mpopt = mpoption('opf.ac.solver', 'FMINCON'); >> results = runopf(mpc, mpopt); Best, Ray On Sep 28, 2022, at 8:18 AM, naime

Re: Enter data for normally open ties

I believe treating them as branches is probably the best way in MATPOWER. Then you can simply use the BR_STATUS column to indicate whether they are open or closed. You’ll have to keep a separate list of which branches are ties. If it is convenient, you could put a list of the branch row indices

Re: Regarding to first order algorithms for AC-OPF

The AC OPF problem is solved by any of the nonlinear programming (NLP) solvers supported by MP-Opt-Model, which currently includes MIPS, fmincon, IPOPT, and Artelys Knitro. Some, like fmincon include multiple algorithms, but I don’t believe any are

Re: Load Flexibility

Sorry, no video about it, but you can use the load2disp() function to turn your loads into dispatchable loads, then they will be dispatched based on cost, just like a generator. Ray > On Sep 7, 2022, at 6:16 AM, Lois Efe wrote: > > Hi all, > > Can anyone recommend a link or video on

Re: Every steps CPF data - Reg

See Table 5-1 in the MATPOWER User’s Manual. The complex voltages are found in results.cpf.V. You can compute the real and reactive powers directly using the corresponding lambda values found in results.cpf.lam. Essentially use the lambda

Re: How to get the PMU data for IEEE 30 bus system

I’m sorry, but this does not sound like a MATPOWER question to me … maybe there is a forum somewhere for Simscape Electrical toolbox questions. Ray > On Sep 3, 2022, at 11:08 AM, naime ahmadi wrote: > > Hello, > I want to get PMU data for IEEE 30 bus system. According to my

Re: LODF but for reconnection of a open line

ing the LODF. I will try testing against results from PowerWorld and see if the -LODF approach works. Richard On Wed, 31 Aug 2022 at 22:57, Ray Daniel Zimmerman mailto:r...@cornell.edu>> wrote: Hello Richard, Do you know the precise definition of the LCDF? The LODF_ij is the chan

Re: LODF but for reconnection of a open line

Hello Richard, Do you know the precise definition of the LCDF? The LODF_ij is the change in flow in branch i resulting from a branch j outage as a fraction of the pre-outage flow on branch j. So would the LCDF_ij be the change in flow in branch i as a result of a branch j clousre as a fraction

  1   2   3   4   >