Re: Trying to solve runopf for t_case30_userfcns involving toggle_iflims

2024-06-29 Thread Ray Daniel Zimmerman
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

Trying to solve runopf for t_case30_userfcns involving toggle_iflims

2024-06-29 Thread Dr. Aamir Nawaz
Dear Dr. Ray, I am trying to do a very simple task as given in your manual. I am using matpower 7.1 version. I want to perform the following: mpc = loadcase('t_case30_userfcns'); mpc = toggle_iflims(mpc, 'on'); results = runopf(mpc); but I am constantly getting the error message as below:

Re: objective function for renewable energy generation

2024-06-10 Thread Ari N
Many thanks On Tue, Jun 11, 2024 at 7:07 AM Ray Daniel Zimmerman wrote: > It is certainly a valid option to model net load and incorporate the > renewables into the load values in the bus matrix. That is equivalent to > modeling it separately as a generator with a very low cost that forces it >

Re: Base kV for 14 bus system

2024-06-10 Thread Ray Daniel Zimmerman
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

Re: objective function for renewable energy generation

2024-06-10 Thread Ray Daniel Zimmerman
It is certainly a valid option to model net load and incorporate the renewables into the load values in the bus matrix. That is equivalent to modeling it separately as a generator with a very low cost that forces it to be dispatched. Ray On Jun 8, 2024, at 4:00 AM, Ari N wrote: Thanks

Base kV for 14 bus system

2024-06-10 Thread Muhammad Junaid
Hi Sir, I hope you are doing well. I am writing to ask about the base kV (kilovolts) for the 14-bus system in MATPOWER. Specifically, I am working on a project that requires accurate voltage information for the 14-bus system, and understanding the base kV is essential for my analysis. Could you

Re: objective function for renewable energy generation

2024-06-08 Thread Ari N
Thanks 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 wrote: > I suppose it depends

Re: objective function for renewable energy generation

2024-06-07 Thread Ray Daniel Zimmerman
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.

objective function for renewable energy generation

2024-05-31 Thread Ari N
Healthy greetings, In generator scheduling, what is the objective function for renewable energy generation? I hope someone can help, thank you. Regards, Ari

Re: CPF stopping Criteria

2024-05-21 Thread Abhyankar, Shrirang G
Yes, CPF can stop at a given lambda value. For example, mpopt = mpoption; mpopt.cpf.stop=1.0 runcpf(‘case9’,’case9target’,mpopt); will stop the CPF when lambda reaches 1.0. Shri From: bounce-128220963-83436...@list.cornell.edu on behalf of Ahmad Sadiq Abubakar Date: Tuesday, May 21, 2024

CPF stopping Criteria

2024-05-21 Thread Ahmad Sadiq Abubakar
Dear all, Please can the CPF stop at a target lambda value? If yes, What will be the field of this cpf_stop_at? I know for "FULL" and "NOSE" are to complete the cpf and stop at the nose point respectively. In anticipation of your response -- *AHMAD Abubakar Sadiq (PhD)* *Faculty Member

Re: UPFC modelling in MATPOWER

2024-05-20 Thread Muhammad Junaid
Thank you Sir for the clarification. @Ahmad Sadiq Sir, I am up for discussion. Thanks On Wed, May 15, 2024 at 3:32 AM Ahmad Sadiq Abubakar < ahmad.abuba...@futminna.edu.ng> wrote: > I have work around UPFC with the fictitious buses in MATPOWER. You may > send me a message to schedule discussion.

Re: UPFC modelling in MATPOWER

2024-05-14 Thread Ahmad Sadiq Abubakar
I have work around UPFC with the fictitious buses in MATPOWER. You may send me a message to schedule discussion. On Fri, May 10, 2024, 10:35 AM Muhammad Junaid wrote: > Dear Sir, > I hope you are doing well. Can you guide how to create two > fictitious buses in Matpower for modeling UPFC? or is

Re: UPFC modelling in MATPOWER

2024-05-14 Thread Ray Daniel Zimmerman
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

UPFC modelling in MATPOWER

2024-05-10 Thread Muhammad Junaid
Dear Sir, I hope you are doing well. Can you guide how to create two fictitious buses in Matpower for modeling UPFC? or is there any guide for modeling UPFC provided by Matpower? Thanks and regards

Penalty Factor Calculations in MATPOWER

2024-04-05 Thread Mohammed Alhajri
Dear Team Could you guide me on implementing Penalty Factor Calculations in MATPOWER? Regards

Re: Inquiry about Bus Voltage Variation in MOST Examples

2024-03-18 Thread Ray Daniel Zimmerman
The time resolution is completely up to you when you are running your own sequence of snapshots using runopf() as in your example. That is, a MATPOWER OPF is for a specific instant in time, so it knows nothing of any time resolution. You are free to define whatever tilme step you like for

Re: Inquiry about Bus Voltage Variation in MOST Examples

2024-03-15 Thread Lichen Wu
Thank you, Ray! It really helps. for i = 1:length(profiles) mpc.bus(2, PD) = 150 * profiles(i); results = runopf(mpc); vm_results(i) = results.bus(2, VM); end Could you please clarify the defined time resolution, such as timestep, microsecond, second, or hour? If so, is it feasible to

Re: Inquiry about Bus Voltage Variation in MOST Examples

2024-03-12 Thread Ray Daniel Zimmerman
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

2024-03-12 Thread amir ali Hosseini
Hello sir, I hope you are well and healthy! If you need to have a parallel load in that branch, you can use column five and six. On Tue, Mar 5, 2024 at 8:30 PM Ray Daniel Zimmerman wrote: > That is correct … assuming you meant to say *column* 5 and 6 of the bus > matrix. > >Ray > > On Mar

Inquiry about Bus Voltage Variation in MOST Examples

2024-03-07 Thread Lichen Wu
Dear All, I hope this message finds you well. I am currently working on replicating a specific scenario from research paper, https://doi.org/10.1016/j.apenergy.2020.115523, voltage smoothing control. I attempted to reproduce the bus voltage variation every second (shown in fig1.fig19.png)

Re: Distributed Generators Model

2024-03-05 Thread Ray Daniel Zimmerman
That is correct … assuming you meant to say column 5 and 6 of the bus matrix. Ray On Mar 4, 2024, at 2:50 AM, Muhammad Junaid wrote: For shunt compensation, do we have to add the shunt compensation value at bus 5 and 6 of bus data, and for series compensation do we have to add it to X

Re: Distributed Generators Model

2024-03-04 Thread Muhammad Junaid
For shunt compensation, do we have to add the shunt compensation value at bus 5 and 6 of bus data, and for series compensation do we have to add it to X (column 4) of branch data? Is this the correct way? On Thu, Nov 23, 2023 at 5:31 PM Muhammad Junaid wrote: > Thank you so much for the

Sv: Proposed 533-bus system based on real data

2023-12-18 Thread Gabriel Malmer
Good, the case files are now uploaded to GitHub as pull requests (case533mt_lo.m and case533mt_hi.m). Best, Gabriel Från: bounce-127945200-95792...@list.cornell.edu För Ray Daniel Zimmerman Skickat: den 15 december 2023 17:15 Till: MATPOWER-DEV-L Kopia: lovisa.tho...@svk.se Ämne: Re:

Re: Proposed 533-bus system based on real data

2023-12-15 Thread Ray Daniel Zimmerman
Sounds good … you might also something like case533mt_lo.m and case533mt_hi.m for the names. Not a big deal either way to me. Ray On Dec 15, 2023, at 7:33 AM, Gabriel Malmer wrote: Hi again, Ray, thank you for suggestion, I was thinking about that as well. I will then upload them as two

Sv: Proposed 533-bus system based on real data

2023-12-15 Thread Gabriel Malmer
Hi again, Ray, thank you for suggestion, I was thinking about that as well. I will then upload them as two separate case files (case533mt.m and case533mt2.m), with description of the loads in the help text. Mirko, I investigated the rated currents when merging parallel branches according to

Re: Proposed 533-bus system based on real data

2023-12-12 Thread Ray Daniel Zimmerman
Hi Gabriel, I prefer to avoid anything that requires modifying a file in the official distribution, even if just to uncomment a single line. So, there are two options for including high load data without putting it in the same case file. 1. Provide the high load case as a separate case file

Sv: Proposed 533-bus system based on real data

2023-12-11 Thread Gabriel Malmer
Hello, Ray, thank you for looking into our case file. We are very glad that you want to include this in the next MATPOWER release. The comments and copyright + license lines are perfectly fine. As an appendix to the file, I added the opportunity to also run the maximum net load hour of 2022.

Re: Proposed 533-bus system based on real data

2023-12-07 Thread Mirko Todorovski
Hello, I'm glad there is new test system based on real data. While browsing through http://lup.lub.lu.se/student-papers/record/9124806, I noticed something worth discussing. Figure 3.3 (p. 20) illustrates merging parallel branches into a single equivalent branch by simply adding the

Re: Proposed 533-bus system based on real data

2023-12-07 Thread Ray Daniel Zimmerman
For general MATPOWER help, please post (a new thread, not just reply to an existing one) on the MATPOWER Discussion e-mail list, not this developer’s list. For getting started with using MATPOWER for electricity markets, check out Appendix F in

Re: Proposed 533-bus system based on real data

2023-12-06 Thread Naresh Boda
Thank you so much sir. Sir, I need help regarding the MATPOWER toolbox, I am the new to use MATPOWER, how to use the matpower in electricity 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

Re: Proposed 533-bus system based on real data

2023-12-06 Thread Ray Daniel Zimmerman
Hello Gabriel, I’d be happy to include this in the next MATPOWER release. I’ve attached a version with a few tweaks to the comments at the beginning of the file (no data changes). Please check that they are ok with you and the providers of the original data. I added some of the description of

Re: Power Flow Calculation Error

2023-12-06 Thread Ray Daniel Zimmerman
It will be included in the next MATPOWER release (planned to be 8.0). But, if you want to use the updated version now, you can either … 1) clone/download and install the current MATPOWER master branch (to download, click the green “Code” button, then select

Re: Power Flow Calculation Error

2023-12-05 Thread Muhammad Junaid
With reference to issue 210 , how can we utilize the updated file? Will it be integrated into MATPOWER 7.1, or do we need to download the updated radial.pf separately? On Wed, Nov 22, 2023 at 10:16 PM Ray Daniel Zimmerman wrote: > What you show

Re: Proposed 533-bus system based on real data

2023-11-30 Thread Ray Daniel Zimmerman
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

2023-11-30 Thread Ray Daniel Zimmerman
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

Proposed 533-bus system based on real data

2023-11-29 Thread Gabriel Malmer
Dear Matpower-developers, My name is Gabriel Malmer and I'm a PhD student in power systems at Lund University. For my master's thesis, me and my colleague Lovisa Thorin developed a 533-bus distribution system in Matpower, based on real system data from the local DSO Kraftringen in southern

Overload Condition

2023-11-28 Thread Scott seedell
Hi, I'd like to be able add some new generation and monitor the branches for overcurrent/overvoltage. Is there a relatively easy way to extract this info, perhaps by comparing Branch Data before and after adding additional generation? Regards

Re: Distributed Generators Model

2023-11-23 Thread Muhammad Junaid
Thank you so much for the detailed answer. On Wed, Nov 22, 2023 at 10:11 PM Ray Daniel Zimmerman wrote: > 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

Re: Power Flow Calculation Error

2023-11-22 Thread Ray Daniel Zimmerman
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

2023-11-22 Thread Ray Daniel Zimmerman
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

2023-11-22 Thread Ray Daniel Zimmerman
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: Power Flow Calculation Error

2023-11-18 Thread Muhammad Junaid
When the generator matrix has the same bus numbers for new rows (other generators placed on the same bus) it gives an error for the power summation method. However, if we add those P and Q's to the same bus for multiple DG's P and Q and represent it as a single bus then it doesn't give an error.

Power Flow Calculation Error

2023-11-18 Thread Muhammad Junaid
Hi Sir, I hope you are doing well. Sir, when we use this data for the generator matrix, it gets solved using Newton's Power Flow but gives this error: > MATPOWER Version 7.1, 08-Oct-2020 -- AC Power Flow (Power Summation) > Warning: Matrix is singular to working precision. > > In calc_v_pq_sum

Distributed Generators Model

2023-11-17 Thread Muhammad Junaid
> > Hi Sir. I hope you are doing well. If EVs are to be added as loads on > buses that ave P and Q. Should the load be added in bus data or should I be > added as shunt element or can it be used as negative generator? You replied > If they are considered as fixed P and Q they can be added to

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

2023-11-10 Thread Ray Daniel Zimmerman
FYI, please start a new thread for a new question, as opposed to just replying to an existing thread. If they are considered as fixed P and Q they can be added to any existing load in the bus matrix or they can be added as a negative generator. If you add it as a shunt (in the bus matrix) it

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

2023-11-10 Thread Muhammad Junaid
Hi Sir. I hope you are doing well. If EVs are to be added as loads on buses that ave P and Q. Should the load be added in bus data or should I be added as shunt element or can it be used as negative generator?

Re: Piecewise Linear Generator Cost Function

2023-11-10 Thread Ray Daniel Zimmerman
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:

Piecewise Linear Generator Cost Function

2023-11-10 Thread Ronald Cabaoig
Dear MATPOWER team, I'd like to clarify some things about the piecewise linear generator cost function. *Question 1*: Suppose I have a cost function as shown: [image: image.png] Is the following the right gencost for the graph above: mpc.gencost = [1 0 0 3 0 0 20 400 30 750]; ?

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

2023-11-09 Thread Jose Luis Marín
I know you want to convert in the direction CIM CGMES ==> MATPOWER, but just in case, the PowSyBl project has converters in the other direction: *

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

2023-11-09 Thread oliver.zehnder
Dear Ray, thank you very much for your prompt feedback. Consequently, we will have to find another way. Best regards, Oliver Von: bounce-127879544-95614...@list.cornell.edu Im Auftrag von Ray Daniel Zimmerman Gesendet: Donnerstag, 9. November 2023 17:25 An: MATPOWER-L Betreff: Re:

Re: I need help regarding MATPOWER implementation for my research

2023-11-09 Thread Ray Daniel Zimmerman
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

2023-11-09 Thread Ray Daniel Zimmerman
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

I need help regarding MATPOWER implementation for my research

2023-11-09 Thread Naresh Boda
dear sir , Please check the doc below. . Thanks & Regards Naresh Boda Research Scholar, EED MNNIT Allahabad, Prayagraj Help.docx Description: MS-Word 2007 document

Convertion of CIM CGMES format into a MATPOWER case struct

2023-11-09 Thread oliver.zehnder
Dear all, ENTSO-E, the European Network of Transmission System Operators for Electricity, publishes datasets representing the power systems of the ENTSO-E members in the CIM CGMES format. Is there a converter available for MATPOWER, to convert data in the CIM CGMES format into a MATPOWER case

Re: How to run single side auction in MATPOWER

2023-11-08 Thread Ray Daniel Zimmerman
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

2023-11-08 Thread Ray Daniel Zimmerman
See Section 6.4.2 on Dispatchable Loads in the MATPOWER User’s

How to run single side auction in MATPOWER

2023-11-08 Thread Dr. D. Karthikaikannan .
Dear sir, how to run a single side auction in MATPOWER where we submit either offer or bid not both.

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

2023-11-05 Thread Ronald Cabaoig
Dear MATPOWER team, In connection with placing a dummy generator to maximize/minimize the power injection at the bus of interest, is it also possible in MATPOWER to impose a power factor constraint on the said dummy generator? Examples: - Scenario 1: minimize the negative power injection in bus

Re: Question on running MATPOWER runoff in LabVIEW mathscript

2023-11-02 Thread Ray Daniel Zimmerman
If it’s a MATPOWER error, I may be able to help, but would need the details of where in MATPOWER the error is happening. Unfortunately, if it is a LabVIEW error, I’m afraid I have no way to help you, since I do not have access to LabVIEW. Ray On Nov 1, 2023, at 2:12 AM, Kevin Yong

Re: Question on running MATPOWER runoff in LabVIEW mathscript

2023-11-01 Thread Kevin Yong
Dear Ray, I have tried to load runpf.m in LabVIEW MATLAB script. However, I am running into error saying that “input must be a MAT-file or an ASCII file containing numeric data with same number of columns in each row” I am stuck and have no idea how to proceed. Would appreciate if you have any

Re: Question on running MATPOWER runoff in LabVIEW mathscript

2023-10-27 Thread Ray Daniel Zimmerman
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

2023-10-27 Thread Ray Daniel Zimmerman
No. It is the total cost (based on offers) of generation, minus the total value (based on bids) of the dispatchable loads. If you have a case where all of the loads are dispatchable (not true for t_auction_case) and have bids that accurately reflect their value, and generators offer their true

Re: What the objective function denotes in the smart market

2023-10-27 Thread Dr. D. Karthikaikannan
thank you sir. is that the objective function value 2985.78 $/hr represents the social welfare (supplier profit for the given offers +consumer net surplus for the given bid). From: bounce-127850109-88225...@list.cornell.edu on behalf of Ray Daniel Zimmerman

Question on running MATPOWER runoff in LabVIEW mathscript

2023-10-26 Thread Kevin Yong
Hello Everyone, I was wondering if it is possible to run MATPOWER power flow m file in LabVIEW math script. A bit of background of my work. I am currently constructing a digital twin of a microgrid in LabVIEW and would like to power flow in LabVIEW. I have configured a case file specifically

Re: What the objective function denotes in the smart market

2023-10-26 Thread Ray Daniel Zimmerman
I see what’s going on. Turns out the explanation is simple. There are 2 sets of “costs” that are in play here. 1. The original generator costs in mpc.gencost. These represent the actual costs of the units. 2, The “costs” derived from the bids and offers. These are used for the actual

Re: What the objective function denotes in the smart market

2023-10-26 Thread Dr. D. Karthikaikannan .
Dear sir, this is the coding and results clc clear mkt.OPF='AC'; mkt.auction_type=1; mpc=loadcase('t_auction_case'); define_constants offers.P.qty=[... 12 24 24 12 24 24 12 24 24 12 24 24 12 24 24 12 24 24]; offers.P.prc=[... 20 50 60; 20 40 70; 20

Re: What the objective function denotes in the smart market

2023-10-19 Thread Ray Daniel Zimmerman
Can you send me the details I need to reproduce your result? The simple OPF on that case does show consistency between the objective function value and the sum of the generator costs. >> r = runopf('t_auction_case'); >> sum(totcost(r.gencost, r.gen(:, PG))) ans = -245.2078 >> r.f ans

Re: What the objective function denotes in the smart market

2023-10-17 Thread Dr. D. Karthikaikannan
this is the result i got it when i run for t_auction_case | Market Summary |

Re: What the objective function denotes in the smart market

2023-10-16 Thread Ray Daniel Zimmerman
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

What the objective function denotes in the smart market

2023-10-14 Thread Dr. D. Karthikaikannan
Dear sir, In optimal power flow the objective function is sum of the cost function and benefit function of generator and load respectively. I run smart market and cross checked the objective function displayed with the corresponding cost function and generation using total cost function and it

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

2023-10-13 Thread Ronald Cabaoig
Thank you very much! Cheers, Ronald On Sat, Oct 14, 2023 at 1:33 AM Ray Daniel Zimmerman wrote: > That’s correct. Except, if you are zeroing out all other costs, there is > no need to use such a huge magnitude for the cost. Any reasonable non-zero > number would do (like 10 or 100). When the

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

2023-10-13 Thread Ray Daniel Zimmerman
That’s correct. Except, if you are zeroing out all other costs, there is no need to use such a huge magnitude for the cost. Any reasonable non-zero number would do (like 10 or 100). When the problem includes coefficients with huge ranges it can potentially result in numerical issues. Ray

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

2023-10-12 Thread Ronald Cabaoig
Hi sir, Thanks for the suggestion! Let me confirm my understanding: 1. Create a dummy generator with the following limits: -9 MW <= Pdummy <= 9 MW 2. To minimize the injection, set a very large positive cost coefficient (C = 9Pdummy) and zero out all other costs, turning the

自动回复: Re: Finding the maximum positive (and negative) power injection at a bus

2023-10-09 Thread 1097997121
这是来自QQ邮箱的假期自动回复邮件。   您好,我最近正在休假中,无法亲自回复您的邮件。我将在假期结束后,尽快给您回复。

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

2023-10-09 Thread Ray Daniel Zimmerman
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

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

2023-10-05 Thread Ronald Cabaoig
Dear MATPOWER team, Greetings! In the current MATPOWER version, is it possible to reformulate the objective function to find the maximum positive (and negative) power injection in a certain bus: Maximize the positive power injection in bus X (or could also be: minimize the negative power

自动回复: Re: Plot SE gain matrix

2023-10-02 Thread 1097997121
这是来自QQ邮箱的假期自动回复邮件。   您好,我最近正在休假中,无法亲自回复您的邮件。我将在假期结束后,尽快给您回复。

Re: Plot SE gain matrix

2023-10-02 Thread Ray Daniel Zimmerman
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

Plot SE gain matrix

2023-09-28 Thread Themistoklis Xygkis
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 invG are the variances of the estimates of the state variables. So, I want to plot them and I found a small detail; since the index

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

2023-09-27 Thread Ray Daniel Zimmerman
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.

2023-09-27 Thread Ray Daniel Zimmerman
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

2023-09-27 Thread Ray Daniel Zimmerman
Thanks for the suggestion. The best way to register this sort of suggestion in a way that I won’t lose track of it is to submit it as a new issue here. Ray On Sep 24, 2023, at 6:04 AM, Ігор Ткач wrote: checklimits is a very useful

Re: Understanding the Profile Struct in MOST

2023-09-27 Thread Ray Daniel Zimmerman
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

MOST Optimization Issue with PV Generation Less Than Half of Load

2023-09-25 Thread keir steegstra
Hi all, Background: I am working on a project that involves simulating a system using MATPOWER's MOST function. The system consists of a single bus with a PV generator, a load, and a battery storage. The objective is to optimize the operation of the battery based on the PV generation profile

loss optimization. Case with pv bus.

2023-09-24 Thread Ігор Ткач
Dear all, As loss optimization for case 118. Case118 with pv bus. costcase(1,6) = 1 is not enough. It is necessary to have the total active power for generators and loads. Igor.

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

2023-09-24 Thread Ігор Ткач
checklimits is a very useful feature. For ease of transition from text to numbers, I humbly suggest the following addition to the checklimits function for a future version of Matpower. Indexes for out of limits : [check_index.]=checklimits where [check_index = 0 - all o'key check_index=1 - Pgen

Re: Lagrangian multiplier and kuhn tucker multiplier for loss reduction

2023-09-22 Thread Ray Daniel Zimmerman
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

2023-09-22 Thread Ray Daniel Zimmerman
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

Understanding the Profile Struct in MOST

2023-09-22 Thread keir steegstra
Hello everyone, I'm working on a project where I'm modeling two isolated towns as a two-bus system. Each town (represented as a bus) will have a storage system, load, and a PV generator. My objective is to determine the feasibility of connecting these towns with a lossy link (initially using a

Re: Lagrangian multiplier and kuhn tucker multiplier for loss reduction

2023-09-22 Thread Dr. Jora Gonda , M
Dear Dr. AK, Please bring the full model. On Fri, 22 Sep 2023, 14:11 Dr. D. Karthikaikannan ., < karthikaikan...@eee.sastra.edu> wrote: > Dear sir, > When we do loss reduction in opf by equating all gencost same, the lam_p > and others represent cost change with respect bus powers etc. How to

自动回复: Lagrangian multiplier and kuhn tucker multiplier for loss reduction

2023-09-22 Thread 1097997121
这是来自QQ邮箱的假期自动回复邮件。   您好,我最近正在休假中,无法亲自回复您的邮件。我将在假期结束后,尽快给您回复。

Lagrangian multiplier and kuhn tucker multiplier for loss reduction

2023-09-22 Thread Dr. D. Karthikaikannan .
Dear sir, When we do loss reduction in opf by equating all gencost same, the lam_p and others represent cost change with respect bus powers etc. How to find lagrangian and kuhn tucker multiplier that represent change in loss with respect to bus powers etc

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

2023-09-22 Thread Ігор Ткач
Is it possible to [Index_check]=checklimits(results, 1) output Index_check = 1 or 0 when there is an output of constrains or not out of limits? 20 вересня 2023, 16:45:03, від "Ігор Ткач" : Dear all, Question. how do I find the constraint. which removes my cost calculation when calculating the

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

2023-09-20 Thread Ray Daniel Zimmerman
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

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

2023-09-20 Thread Ігор Ткач
Dear all, Question. how do I find the constraint. which removes my cost calculation when calculating the runopf cost? For case30, I study the dependence of cost when adjusting the voltage of   generators. For example, I change the voltage of the 4th generator: cost=function('VM')  ==  

Re: Battery storages

2023-09-13 Thread Ray Daniel Zimmerman
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

Battery storages

2023-09-12 Thread Heller , Leopold Bernd Jürgen
Hello everyone, Recently, I carried out an AC OPF for a high-voltage power grid. With roughly 70% of the installed net capacity from renewables, I'm curious if Matpower offers a battery storage model. I'm considering a basic model where it behaves like a generator when injecting power into

Re: Jacobian Matrix

2023-09-12 Thread Ray Daniel Zimmerman
If I understand, you are attempting to find the sensitivity of bus voltage magnitudes to active power injections at load (PQ) buses, under standard power flow asusmptions, right? In that case, I think you want … J = makeJac(mpci); invJ = inv(full(J)); dVdP = invJ(118:end, 54:117); The rows of

Re: Jacobian Matrix

2023-09-09 Thread Kaled Zareer
Thank toy for your response. Here is the code.  define_constants;mpcb = loadcase('case118');mpct = mpcb;mpct.bus(:, [PD QD]) = mpct.bus(:, [PD QD])*1.0;mpci = ext2int(mpct);Jacc = makeJac(mpci);Jac=full(Jacc); Jac_inv= inv(Jac);% inverse of Jsen_vpp=Jac_inv(1:117,118:181); % sensitivity of

  1   2   3   4   5   6   7   8   9   10   >