Re: Unbalanced 3 phase load flow and IEEE Test Systems

I don’t believe so. Ray On Sep 23, 2024, at 9:54 AM, Oguzhan Ceylan wrote: Dear All, I was wondering if any of the IEEE distribution network test systems (except IEEE 4 node) has been implemented in Matpower? I have checked this link: How to Run a Three-Phase Power Flow — MATPOWER Docu

Re: Running load flow on three phase and single phase system.

MATPOWER 8 does have some prototype 3-phase modeling (see How to Run a Three-Phase Power Flow). While it does allow for connecting a 3-phase unbalanced network (e.g. distribution) to a balanced single-phase representation of a 3-phase network (e.g

Re: Inquiry About Dynamic Load Flow Analysis in MATPOWER

MATPOWER does not directly support dynamic load flow analysis. You may want to check out MatDyn, though I have no personal experience with it. Ray On Sep 5, 2024, at 7:17 AM, Muhammad Junaid wrote: Dear Sir, I hope this message fin

Re: Trying to solve runopf for t_case30_userfcns involving toggle_iflims

Those command work fine on 7.1 for me. I suspect that maybe something got messed up with your installation. You might try reinstalling a fresh cope of MATPOWER. Does test_matpower pass all of the tests for you? And if you are going to reinstall, I suggest trying MATPOWER 8

Re: Base kV for 14 bus system

I’m afraid the only information I have about the system is what is available here … https://labs.ece.uw.edu/pstca/pf14/pg_tca14bus.htm That was the source of the data. Ray On Jun 10, 2024, at 8:38 AM, Muhammad Junaid wrote: Hi Sir, I hope you are doing well. I am writing to ask ab

Re: objective function for renewable energy generation

for the answer. Is it possible for renewables to be used to reduce the load? Renewable power can be obtained from weather forecasting (temperature & irradiance for PV or wind speed for wind turbines) Ari On Fri, Jun 7, 2024 at 11:42 PM Ray Daniel Zimmerman mailto:r...@cornell.edu>>

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. Als

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 modelin

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: 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: 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 compare_case()

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 distri

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 AM

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 impeda

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

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: Suppose

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. K

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

= -9Pdummy) and zero out all other costs, turning the objective function to: min (-9Pdummy + 0Pg1 + 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

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 October

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 th

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

du<mailto:bounce-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:

Re: What the objective function denotes in the smart market

the value 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>&

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 the i

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 in

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, 202

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. R

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

- Qgen 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 o

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 modeli

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 sensitivi

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 (n

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 1

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 pr

Re: Jacobian Matrix

= mpct.bus(:, [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 yo

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. So,

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 voltage

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, A

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 ad

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 On

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 opf.us

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 callback

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

with is to define args as a global 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 term

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 m

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

t regards, Yuanxi Wu 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

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, 2023,

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 User’s

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 cont

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 capa

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), wher

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 to

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 exce

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 us

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 colum

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>> wrote:

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 sur

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 mailto:

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 M

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 co

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, 2022,

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 po

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 th

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 t

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 insta

Re: Load Flexibility

es 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 Zim

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

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

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 can

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); Here

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 poin

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 ahmad

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 firs

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 how

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 value

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 understanding,

Re: LODF but for reconnection of a open line

e computing 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

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

Re: mpc.branch

Get Outlook for Android<https://aka.ms/AAb9ysg> From: bounce-126765030-91948...@list.cornell.edu<mailto:bounce-126765030-91948...@list.cornell.edu> mailto:bounce-126765030-91948...@list.cornell.edu>> on behalf of Ray Daniel Zimmerman mailto:r..

Re: mpc.branch

Hi Scott, The explanation is in the fine print … namely the 3rd footnote under Table B-3. Those extra columns are included in the output but are typically not included (or else they are ignored) on input. Ray On Aug 22, 2022, at 3:16 PM, Scott seedell mailto:scottseed...@hotmail.com>>

Re: Active and reactive power losses - Question

Hi Fábio, Your Qloss includes both the series reactive losses returned by get_losses() in the imaginary part of the loss return value as well as the line charging capacitance returned in the chg (or fchg and tchg) return value. See help get_losses() for details about the line charging part.

Re: [EXTERNAL] TSPOPF PDIPM AC solver not a valid value

pt); Error in runopf (line 75) [r, success] = opf(casedata, mpopt); Which I don’t know how to interpret. Thanks, Ivy From: bounce-126635453-91654...@list.cornell.edu<mailto:bounce-126635453-91654...@list.cornell.edu> mailto:bounce-126635453-91654...@list.cornell.edu>> On Behalf O

Re: TSPOPF PDIPM AC solver not a valid value

That error message is very strange. Are you sure you are using straight, single-quotes around the string arguments to mpoption? Like this ... mpopt = mpoption('opf.ac.solver','PDIPM') If so, I’d try reinstalling MATPOWER and TSPOPF and trying again. Ray On Jun 21, 2022, at 4:22 PM, Rachae

  1   2   3   4   >