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

Re: Question regarding Matpower's continuation power flow

2022-11-01 Thread Abhyankar, Shrirang G
Hi Carlos, Event: See MATPOWER manual section 5.5 for description of ‘events’ in MATPOWER. Basically, an event is a point on the continuation curve where some sort of limit or threshold is met. For example, a generator hitting its active/reactive power limit is an event. The full list of

Re: Base Case CPF - Reg

2021-09-21 Thread Abhyankar, Shrirang G
It is hard to say anything without knowing exactly what commands you ran. Please refer to the script on page 59 of the manual. Thanks, Shri From: on behalf of "A.K Samy" Reply-To: MATPOWER discussion forum Date: Tuesday, September 21, 2021 at

Re: Request for technical data - Reg

2021-09-16 Thread Abhyankar, Shrirang G
The load real power (Pd) for a lambda value (lam) is given by Pd = (1-lam)*Pd_b + lam*Pd_t Where Pd_b and Pd_t are the load real powers for the base and target case, respectively. Similar calculation applies for reactive power. Shri From: on behalf of "A.K Samy" Reply-To: MATPOWER

Re: Continuous power flow question

2021-06-01 Thread Abhyankar, Shrirang G
Hi Simone, See the FAQ on power flow divergence. It lists some useful tips to debug power flow divergence when using MATPOWER and how to run a continuation power flow to check if a power flow solution exists (point v). Thanks, Shri From: on

Re: CPF power flow did not converge

2019-07-11 Thread Abhyankar, Shrirang G
Cynthia, Please see this MATPOWER FAQ on debugging power flow divergence. Also, you forgot to attach your case file. Thanks, Shri From: on behalf of W Ruoyu Reply-To: MATPOWER discussion forum Date: Thursday, July 11, 2019 at 10:54 AM To:

Re: Question about CPF

2019-07-02 Thread Abhyankar, Shrirang G
See below. Thanks, Shri From: on behalf of W Ruoyu Reply-To: MATPOWER discussion forum Date: Monday, July 1, 2019 at 5:35 PM To: MATPOWER discussion forum Subject: Re: Question about CPF Hi, Thanks for the reply ! That really helps. But I still have a little questions. To output the

Re: Maximum load limit for radial distribution cases.

2019-06-13 Thread Abhyankar, Shrirang G
For maximum loadability calculations, it is best to run a continuation power flow (CPF) rather than power flow with successive load increments. CPF avoids ill-conditioning of the matrix near collapse point which a simple power flow cannot. One can run a CPF in MATPOWER using the runcpf()

Re: Question about power flow for a system with faulted bus

2019-06-05 Thread Abhyankar, Shrirang G
, CPF is widely used to solve low voltage problem. Any idea? Thank you once again. Ismael On Jun 1, 2019, at 7:38 AM, Abhyankar, Shrirang G mailto:shrirang.abhyan...@pnnl.gov>> wrote: Ismael, Power flow with loads modeled as “constant power” faces convergence issue when the voltages a

Re: Question about power flow for a system with faulted bus

2019-06-01 Thread Abhyankar, Shrirang G
Ismael, Power flow with loads modeled as “constant power” faces convergence issue when the voltages are low. A common trick is to convert constant power loads to constant impedance loads for low voltages. In MATPOWER, you can convert all the loads to constant impedance using the options

Re: Continuation power flow_lambda

2019-05-12 Thread Abhyankar, Shrirang G
Lambda, a scalar quantity, represents the fraction of additional load/generation (transfer) on the network. The continuation power flow finds solutions starting from lambda = 0 (base loading/generation) to lambda_max (nose-point, max. loading/generation). You can translate the lambda value to

Re: Isolated Bus - Power flow solution

2019-04-29 Thread Abhyankar, Shrirang G
John, MATPOWER’s power flow does not do a topology check to find isolated buses and islands. So, any isolated bus i needs to have its bus type set to None (4) -- mpc.bus(i, BUS_TYPE) = 4. If you want to find out which buses get isolated when the status of one or more lines changes then use

Re: Removing a Bus and Evaluating the New Network

2019-04-29 Thread Abhyankar, Shrirang G
Mitchell, See below. Shri From: on behalf of "Dennis, Mitchell T." Reply-To: MATPOWER discussion forum Date: Friday, April 26, 2019 at 8:38 PM To: "MATPOWER-L@cornell.edu" Subject: Removing a Bus and Evaluating the New Network Dear All, I have a question regarding convergence. I am

Re: CPU time

2019-04-12 Thread Abhyankar, Shrirang G
Ahmad, Can you be more specific with your question. How big of a system are you planning to solve? Do you need to perform the CPF run within some stipulated time? CPF basically solves a series of power flow solutions given some transfer direction. Roughly, its computational complexity

Re: case300 problem

2019-04-12 Thread Abhyankar, Shrirang G
Naime, The bus numbers given in the casefile, “external bus numbers”, need not be consecutive integers from 1..nbus. These can be any integers. Internally, MATPOWER converts them to “internal numbering” (consecutive integers from 1..nbus). If you want to operate with internal numbers then

Re: AC power flow not converge

2018-11-14 Thread Abhyankar, Shrirang G.
Hi, The AC power flow (runpf) is converging for your system, while the AC optimal power flow (runopf) does not. The reason being the line ratings (rateA, B, C) are less than the load they are supplying to. For instance, the load on bus 2 is 2.8 MW, but there is a single line 1-2 supplying

Re: Power Flow in Matpower

2018-03-16 Thread Abhyankar, Shrirang G.
See FAQ #5 Thanks, Shri Ph: (630) 252 0219 www.mcs.anl.gov/~abhyshr From: on behalf of Mohammed Alhajri Reply-To: MATPOWER discussion forum

Re: AC OPF Convergence Issue

2018-03-08 Thread Abhyankar, Shrirang G.
ractical solution. Shri sincerely, Aakanksha Dubey On Wed, Mar 7, 2018 at 7:25 PM, Abhyankar, Shrirang G. <abhy...@anl.gov<mailto:abhy...@anl.gov>> wrote: The voltage magnitudes in MATPOWER data files, (bus(:,VM)), are given in per unit and are generally close to 1.0. You have very

Re: AC OPF Convergence Issue

2018-03-07 Thread Abhyankar, Shrirang G.
The voltage magnitudes in MATPOWER data files, (bus(:,VM)), are given in per unit and are generally close to 1.0. You have very low voltages for many of the buses. Similarly, the VMAX and VMIN values are also in per unit and generally they are 1.0 and 0.95, respectively. The following code

Re: Introducing contingency in the network

2018-01-17 Thread Abhyankar, Shrirang G.
MATPOWER’s power flow does not enforce generator active power limits. To do so, run a dc optimal power flow (rundcopf). Thanks, Shri Ph: (630) 252 0219 www.mcs.anl.gov/~abhyshr From: on behalf of Zain Anwer Memon

Re: The question about extend OPF and the optimization tool

2018-01-17 Thread Abhyankar, Shrirang G.
Not sure if this is what you want, but here’s one way. https://www.mail-archive.com/matpower-l@cornell.edu/msg00613.html Thanks, Shri Ph: (630) 252 0219 www.mcs.anl.gov/~abhyshr From: on behalf of

Re: Question about the continuation power flow

2017-12-28 Thread Abhyankar, Shrirang G.
The base case (with PD = QD = PG = 0) is not converging because some of the generators have high Qmin limits. You’ll also need to relax the generator Qmin limits for your base case. Thanks, Shri Ph: (630) 252 0219 www.mcs.anl.gov/~abhyshr From:

Re: PF did not converge

2017-12-14 Thread Abhyankar, Shrirang G.
Try the suggestions in FAQ #5 Thanks, Shri Ph: (630) 252 0219 www.mcs.anl.gov/~abhyshr From: on behalf of Mohammed Alhajri Reply-To:

Re: How to find max-load ability for each bus of power system

2017-11-30 Thread Abhyankar, Shrirang G.
bus is critical for voltage stability studies? On Mon, Nov 27, 2017 at 11:22 PM, Abhyankar, Shrirang G. <abhy...@anl.gov<mailto:abhy...@anl.gov>> wrote: Hello, MATPOWER’s CPF requires two data files; the first one having the base loading/generation, called basecasedata, and the ot

Re: How to find max-load ability for each bus of power system

2017-11-27 Thread Abhyankar, Shrirang G.
Hello, MATPOWER’s CPF requires two data files; the first one having the base loading/generation, called basecasedata, and the other that has the target loading/generation, called targetcasedata. The CPF traces the trajectory of solutions along increasing lambda direction, with lambda = 0 at

Re: PF and OPF did not converge

2017-10-03 Thread Abhyankar, Shrirang G.
http://www.pserc.cornell.edu/matpower/#pfconvergence From: on behalf of Charalambos Ioannou Reply-To: MATPOWER discussion forum Date: Tuesday, October 3, 2017 at 9:12 AM To:

Re: Calculation of r, x, b from PF results

2017-08-03 Thread Abhyankar, Shrirang G.
This is an inverse problem. You'll need to use a nonlinear least squares solver (lsqnonlin in MATLAB) to compute the line parameters given the different measurements. MATPOWER does not have an explicit function to do what you want, but you can use its different functions (or extract code from

Re: Objective Functions

2017-07-15 Thread Abhyankar, Shrirang G.
MATPOWER supports two types of generator cost curves: quadratic and piece-wise linear. I don’t think there is support for user-defined cost curves, Ray would be able to provide a more definitive answer. So, I think there is an no easy way to implement the objective you need. You’ll need to

Re: Get_losses

2017-07-15 Thread Abhyankar, Shrirang G.
dloss_dv is the partial derivative of the losses w.r.t. voltage angle. It is an optional output argument. If you do not need it then don’t use it. FYI: MATPOWER has an online function reference where you can see the usage of all of its

Re: Compare case

2017-07-15 Thread Abhyankar, Shrirang G.
compare_case() cannot be used on two different networks. It compares the runpf() (or runopf()) outputs for two different scenarios on the “same” network. For example, mpc1 = runpf(‘case9’);% run power flow on 9-bus network mpc2 = loadcase(‘case9’);% load 9-bus network

Re: multi-period AC Power Flow (not OPF)

2017-07-15 Thread Abhyankar, Shrirang G.
Using a loop would be the easiest way. Alternately, you could concatenate all the mpc bus, gen, etc. arrays into one giant mpc case and then solve it via a single runpf(). However, this is more complicated (changing of bus,gen, branch numbers for the later time-steps) compared to using a loop.

Re: CPF Callback Functions

2017-07-05 Thread Abhyankar, Shrirang G.
Here’s an example of a CPF call back function. The callback function can be set through MATPOWER options struct mpopt =

Re: Load increase in one busbar using runcpf

2017-05-14 Thread Abhyankar, Shrirang G.
Assuming 101 is the row number in the mpc.bus array for bus 101: define_constants; mpopt = mpoption('out.all', 1, 'verbose', 2); mpopt = mpoption(mpopt, 'cpf.stop_at', 'FULL', 'cpf.step', 0.2,’cpf.plot.bus’,101); mpopt = mpoption(mpopt, 'cpf.plot.level', 2); mpcb = loadcase('Case_04'); % load

Re: CPF termination

2017-05-10 Thread Abhyankar, Shrirang G.
m_event as template, do I also need to write a callback function for my Vlimts and Llimits event? or modify the default callback function for cpf termination? Kindly excuse my too basic questions, I am novice in coding. In anticipation of your reply On Apr 25, 2017 6:45 PM, "Abhyankar, S

FW: MATPOWER, 9 bus System, Need help solving

2017-04-29 Thread Abhyankar, Shrirang G.
Dear Saad, It is always best to send your querries on the MATPOWER mailing list than to direct it to individual members. You will multiple and much better responses by mailing the list. I don’t think MATPOWER has a function that can do the state estimation as you’ve described in the document.

Re: CPF termination

2017-04-25 Thread Abhyankar, Shrirang G.
Ahmad, You need to implement an event function (see section 5.5 of the manual) for cpf with logic on how you want cpf to terminate. There are a few event functions currently implemented — nose point detection, target limit, active/reactive power limits — but there is no user-defined event

Re: Line loading limit

2017-04-19 Thread Abhyankar, Shrirang G.
of 'rundcopf' On Wed, Apr 19, 2017 at 11:39 AM, Abhyankar, Shrirang G. <abhy...@anl.gov<mailto:abhy...@anl.gov>> wrote: runpf does not enforce any line limits. You need to run optimal power flow (runopf). Shri From: Akash Tyagi <akashtyagi0...@gmail.com<mailto:akashtyagi0.

Re: Line loading limit

2017-04-19 Thread Abhyankar, Shrirang G.
runpf does not enforce any line limits. You need to run optimal power flow (runopf). Shri From: Akash Tyagi > Reply-To: MATPOWER discussion forum > Date: Tuesday, April 18,

Re: regarding CPF

2017-04-13 Thread Abhyankar, Shrirang G.
I think the base case has the generator at bus 8 turned off, while the target case has it ON. My guess is this is causing problems in the CPF when creating the transfer vectors due to differences in sizes of the Pg/Qg. Can you keep the generator ON at bus 8 in the base case with zero real and

Re: False Data Detection in State Estimation

2017-04-02 Thread Abhyankar, Shrirang G.
MATPOWER assumes the network is balanced so only a positive sequence power flow calculation is performed for transmission and distribution. If you’d like to do an unbalanced distribution power flow calculation then you’ll need to use other tools like OpenDSS or GridLab-D. Shri From: Akash

Re: CPF Question

2017-03-06 Thread Abhyankar, Shrirang G.
(FULL) of any arbitrary bus using the CPF in MATPOWER. So ideally what I want is bus voltage plotted with mu. Can i get additional insights in to how can I modify the default callback function keeping in view my requirements? Thanks, Aaqib On Sun, Mar 5, 2017 at 12:22 PM, Abhyankar, Shrirang G

Re: CPF Question

2017-03-05 Thread Abhyankar, Shrirang G.
Hello Aaqib, There is no in-built function in MATPOWER to plot what you need. However, you can do so easily in two ways: 1. Modify the default callback function, cpf_default_callback.m, to plot the active power - voltage curve. Currently, it is set to plot lambda-voltage curve (both

Re: CPF with Matpower

2017-03-02 Thread Abhyankar, Shrirang G.
ile. Last time I hyperlinked with my google drive. bests- Mahi On Mon, Feb 27, 2017 at 8:03 PM, Abhyankar, Shrirang G. <abhy...@anl.gov<mailto:abhy...@anl.gov>> wrote: Hi Mahi, No file was attached with your email. Shri From: Mahbubur R <tanim.3...@gmail.com<mailto:tani

Re: CPF with Matpower

2017-02-27 Thread Abhyankar, Shrirang G.
teps for the cpf? Thanks Mahi On Wed, Feb 22, 2017 at 11:16 AM, Abhyankar, Shrirang G. <abhy...@anl.gov<mailto:abhy...@anl.gov>> wrote: Hello Mahi, basecasedata is a MATPOWER data struct that has the base case loading/generation, while targetcasedata has that for the target

Re: CPF with Matpower

2017-02-22 Thread Abhyankar, Shrirang G.
Hello Mahi, basecasedata is a MATPOWER data struct that has the base case loading/generation, while targetcasedata has that for the target loading/generation. A continuation power flow traces power flow solutions starting from the base case to the target case, if it can be reached. Please

Re: load reduction in continuation power flow (CPF)

2017-02-16 Thread Abhyankar, Shrirang G.
I agree with Ray. The base case dispatch being updated by OPF causes a reduction in generation, commensurate to the load reduction, and hence you are seeing a decrease in the load margin. What happens when you decrease the load at some other buses? Do you get the same behavior? Shri From: Ray

Re: runpf help

2016-10-05 Thread Abhyankar, Shrirang G.
results=runpf(this_case,mpopt); J=full(makeJac(results)); > On Oct 5, 2016, at 6:05 PM, Aaqib Peerzada wrote: > > results=runpf(this_case,mpopt); > J=full(makeJac(this_case));

Re: Problem with runcpf

2016-09-02 Thread Abhyankar, Shrirang G.
I suspect you have isolated buses or islands that do not have a swing bus. Have you tried I-iv in MatPower faq on power flow divergence. http://www.pserc.cornell.edu/matpower/#pfconvergence Shri On Sep 2, 2016, at 6:51 AM, Gamze Dogan >

Re: Admittance formula for branches with transformers

2016-08-26 Thread Abhyankar, Shrirang G.
Please take a look at the function makeYbus(). > On Aug 26, 2016, at 10:05 AM, davor sutic wrote: > > I'm trying to understand the effects of having a transformer connected to a > branch on the implicated admittance matrix. In the MatPower manual, section > 3.2,

Re: How do I print the reduced Ybus matrix

2016-07-27 Thread Abhyankar, Shrirang G.
1. You need to have the MATPOWER directory on your MATLAB path. http://www.mathworks.com/help/matlab/matlab_env/add-remove-or-reorder-folders-on-the-search-path.html 2. YBUS, YF, YT] = makeYbus(mpc) should be YBUS, YF, YT] = makeYbus(mpcreduced); Shri From: Philip Adedotun Olaniyi

Re: How do I print the reduced Ybus matrix

2016-07-27 Thread Abhyankar, Shrirang G.
Use makeYbus() to create the admittance matrix for the reduced network. http://www.pserc.cornell.edu/matpower/docs/ref/matpower5.0/makeYbus.html Shri On Jul 26, 2016, at 10:19 PM, Philip Adedotun Olaniyi > wrote: Good day sir, I am trying to extract

Re: Query about LMP function

2016-06-29 Thread Abhyankar, Shrirang G.
http://www.pserc.cornell.edu//matpower/docs/ref/matpower6.0b1/scale_load.ht ml FYI: The list of available MATPOWER functions can be found here http://www.pserc.cornell.edu//matpower/docs/ref/ Shri From: Priyanka Shinde Reply-To: MATPOWER discussion forum

Re: AW: CPF - Functionality

2016-06-10 Thread Abhyankar, Shrirang G.
t matter how big the Sxfr is. Important is only the relation between the different rows of the vector? Thomas Von: bounce-120551126-75986...@list.cornell.edu<mailto:bounce-120551126-75986...@list.cornell.edu> [mailto:bounce-120551126-75986...@list.cornell.edu] Im Auftrag von Abhy

Re: CPF - Functionality

2016-06-09 Thread Abhyankar, Shrirang G.
1. Yes. 2. CPF with generator active and reactive power limits should be hopefully available in the upcoming version of MATPOWER. Checking for min. or max. bus voltages is not available in CPF. With some modifications in the CPF code, you should be able to check voltages. We were working on

Re: Wind generation

2016-05-22 Thread Abhyankar, Shrirang G.
No. > On May 22, 2016, at 5:25 PM, Mohamed Taleb wrote: > > Dears, > > Greetings, > I there any modeling for the wind turbine that simulates its behavior? > > Regards, > Mohamed Taleb

Re: modification on Matpower case.m file for power flow of distribution system

2016-05-04 Thread Abhyankar, Shrirang G.
http://ewh.ieee.org/soc/pes/dsacom/testfeeders/index.html However, as Jovan mentioned, MATPOWER is for solving balanced three-phase transmission systems. If you are interested in solving unbalanced three-phase power flows for distribution networks then refer to either OpenDSS or GridLAB-D.

Re: scaleload

2016-02-23 Thread Abhyankar, Shrirang G.
mpc=loadcase('case14o’); for t=1:2 mpc(t)=mpc; mpc(t).bus = scale_load(0.5, mpc(t).bus); The error is in the third line of your code. You are resizing the mpc struct created in the first line and hence getting the structure resizing error. Use the following code.

Re: power flow question

2016-02-18 Thread Abhyankar, Shrirang G.
I agree. Having the warnings for all generators, voltages, etc. would be helpful for the cases when users are also varying the generation on PV buses. I would even like MATPOWER fixing the active power at PV buses to the max. or min. limit, if the generation exceeds these limits, and informing

Re: power flow question

2016-02-17 Thread Abhyankar, Shrirang G.
Wenlei, The power flow routine in MATPOWER does not enforce generation limits on real power. Exceeding the capacity does not necessarily imply a divergence of the power flow. Shri From: , Wenlei > Reply-To: MATPOWER discussion forum

Re: Generating Larger Test Cases from existing Matpower Cases

2016-02-15 Thread Abhyankar, Shrirang G.
Paul, To add to this discussion, I recall Carlos mentioning that one needs to be careful when trying to ³stitch² together different systems due to the large angle differences at the interconnection buses, which might lead to solver infeasibility. The following email threads might be helpful in

Re: demand response

2016-02-14 Thread Abhyankar, Shrirang G.
Please ALWAYS send the code you were trying to run and the entire error message. Shri On Feb 13, 2016, at 11:16 PM, Mounika Vanjarapu > wrote: when i tried for case9 or case5 ,its getting error.what is the reason On Wed, Feb 10,

Re: save losses

2016-02-13 Thread Abhyankar, Shrirang G.
The two error messages you’ve sent look unrelated to me. Are these from two different runs? In any case, most likely, you are using an older version of MATPOWER. Download the latest version from here http://www.pserc.cornell.edu//matpower/. If you are still getting errors, please send the

Re: How to model a substation in Matpower

2016-02-12 Thread Abhyankar, Shrirang G.
Dear Matpowers, Hope you are doing great. We are trying to run PF for a distribution network in Matpower. Since the R/X ratio is acceptable we tried NR. It did converge, however, the loss is not correct. The loss formulae used in MATPOWER are given in section 9.2.4 of the manual

Re: demand reponse

2016-02-02 Thread Abhyankar, Shrirang G.
cult.in<http://difficult.in> which way i can know which is most preferrable outage by contingency analysis using matpower. On Tue, Feb 2, 2016 at 10:18 PM, Abhyankar, Shrirang G. <abhy...@anl.gov<mailto:abhy...@anl.gov>> wrote: Each generator, branch has a status that you can change t

Re: demand reponse

2016-02-02 Thread Abhyankar, Shrirang G.
Each generator, branch has a status that you can change to simulate a contingency. See the functions idx_brch() and idx_gen() for the appropriate columns for the status flag in the branch and generator data, respectively. Shri From: Mounika Vanjarapu

Re: Running Available Transfer Capability (ATC-Power System)

2016-01-04 Thread Abhyankar, Shrirang G.
The number of input arguments to runpf() is incorrect. Please send the error log (output in the command window). Shri From: abde wahab naser > Reply-To: MATPOWER discussion forum

Re: ask in matpower 3.2

2015-12-31 Thread Abhyankar, Shrirang G.
You can calculate the branch currents (from and to end) using the real and reactive power flows (from the branch data) and the bus voltages (in the bus data). See the function idx_brch() and idx_bus() for the location of the flows and voltages. Shri From: Mahmoud Abdallah

Re: Difference between System and Load Losses

2015-12-29 Thread Abhyankar, Shrirang G.
Arun, I don’t think there is a “load losses” function in MATPOWER that returns how much load is lost akin to a LOLE calculation. If you want to calculate loss of load for reliability or restoration calculations, then you can run an optimal power flow with dispatchable loads. (set very high

Re: 39-bus system

2015-12-20 Thread Abhyankar, Shrirang G.
Yes. However, note that each area needs to have enough generation capacity to supply the load. Shri > On Dec 19, 2015, at 12:42 PM, "Owusuboateng, Prince Charles" > wrote: > > Hello, > > I have the 39 bus system divided into 4 independent areas. Can I

Re: Continuation Power Flow

2015-12-02 Thread Abhyankar, Shrirang G.
You need to create a callback routine for CPF that plots what you need and set it in the options struct. mpopt = mpoption(‘cpf.user_callback’,’cpf_user_callback’); where cpf_user_callback is the name of your callback routine. Read pages 38 - 40 of the manual. Shri From: Anjali Sharma

Re: ACOPF constraints

2015-11-30 Thread Abhyankar, Shrirang G.
Rectangular. > On Nov 30, 2015, at 5:10 AM, "Vladimir Frolov" > wrote: > > Dear Matpower community, > > Do you know if opportunity cost constraint is used in the ACOPF formulation > of Matpower? > In other words, does an allowed P,Q area for generator is just

Re: VSC-Multiple Terminals

2015-10-20 Thread Abhyankar, Shrirang G.
You can run a load flow with MATPOWER, but it does not have small signal analysis capability yet. Shri On Oct 20, 2015, at 12:08 PM, "Rajesh Vaid" > wrote: Hello Matpower community, I am working on DC grid system where in the system is

Re: Question about sparsity-based implementation in MATPower

2015-10-18 Thread Abhyankar, Shrirang G.
Hi Shruti, The direct linear solver used by MATLAB depends on the symmetry of the Jacobian matrix. For MATPOWER test cases that have symmetric Jacobians (due to inactive taps), a Cholesky factorization is used (LL^T = A). For cases that lead to non-symmetric Jacobian, MATLAB uses UMFPACK for

Re: Question about sparsity-based implementation in MATPower

2015-10-18 Thread Abhyankar, Shrirang G.
ion followed by forward,backward substitutions. From your response I assume that it is true that MATpower uses "\" right? Thank you for your response. On Sun, Oct 18, 2015 at 6:27 PM, Abhyankar, Shrirang G. <abhy...@anl.gov<mailto:abhy...@anl.gov>> wrote: Hi Shrut

Re: Computing "A" matrix for Eigenvalue analysis

2015-09-29 Thread Abhyankar, Shrirang G.
I think you want to calculate the eigen values the Jacobian matrix for the AC-DC power flow. You need to get in touch with MATACDC developers to find out what API is provided for obtaining the Jacobian matrix. MATPOWER has the function makeJac()

Re: Constraint types and its effect on OPF convergence

2015-09-10 Thread Abhyankar, Shrirang G.
I am not sure there is a clear obvious way to tell which constraints pose the most problem for the optimization. Perhaps doing some sensitivity analysis might help, but with different operating conditions you will get different sensitivities. There are a few papers published on convergence

Re: Jacobian matrix

2015-09-07 Thread Abhyankar, Shrirang G.
The size of the square Jacobian matrix, returned by makeJac(), in MATPOWER is 2*npv + npq, where npv is the number of PV buses and npq is number of PQ buses. You can create appropriate permutation vectors to reorganize the matrix entries. However, from what I understand, you want to put the

Re: REQUEST

2015-09-02 Thread Abhyankar, Shrirang G.
What do you mean by 'simple (direct or undirect) network'? Shri On Aug 31, 2015, at 9:06 PM, "Mohamed-Larbi Rebaiaia" > wrote: Hi... I would like to proceed to the transformation of an electric network ( e.g.

Re: REQUEST

2015-08-31 Thread Abhyankar, Shrirang G.
Please take a look at the continuation power flow chapter in the MAtPOWER manual. Shri On Aug 31, 2015, at 11:05 AM, "Owusuboateng, Prince Charles" > wrote: Hi All, I have a divided the IEE 39bus system in 4

Re: REQUEST

2015-08-31 Thread Abhyankar, Shrirang G.
ateng, Prince Charles" <powusuboateng8...@mytu.tuskegee.edu<mailto:powusuboateng8...@mytu.tuskegee.edu>> wrote: Please it has not been of help. I am looking to load each of the 4 regions independently and perform power flows on each one of them. On Mon, Aug 31, 2015 at 2:58 PM,

Re: transformers and tap coefficients

2015-08-26 Thread Abhyankar, Shrirang G.
Yes, it is a typo. It should be r =b =0. Shri On Aug 26, 2015, at 6:47 AM, Vladimir Frolov vladimir.fro...@skolkovotech.ru wrote: Dear Matpower community, My question is about tap coefficients. It is said in the data description that branch includes transformer if r=x=0. But actually

Re: About three winding transformers and contingency analysis

2015-08-23 Thread Abhyankar, Shrirang G.
It would be hard to find which branches in the MATPOWER data constitute a three-winding transformer if you don’t have the original cases. I guess you’ll just have to choose sets of three branches and assume they constitute three-winding transformers. Shri From: Jesus Maria Lopez Lezama

Re: Add a load as a variable to the extended optimal power

2015-08-21 Thread Abhyankar, Shrirang G.
See section 6.4.2 on ‘dispatchable loads’ in the MATPOWER manual. Shri From: Mahdi Kiaee meh19...@yahoo.commailto:meh19...@yahoo.com Reply-To: MATPOWER discussion forum matpowe...@list.cornell.edumailto:matpowe...@list.cornell.edu Date: Friday, August 21, 2015 at 1:35 PM To:

Re: convergence problem in runpf.

2015-08-16 Thread Abhyankar, Shrirang G.
, Abhyankar, Shrirang G. abhy...@anl.govmailto:abhy...@anl.gov wrote: Jose, Would it be fine with you if the steps you’ve mentioned below are added to MATPOWER FAQ#5 http://www.pserc.cornell.edu//matpower/#pfconvergence Many a times, useful and detailed suggestions, such as what you’ve enumerated

Re: convergence problem in runpf.

2015-08-14 Thread Abhyankar, Shrirang G.
,grid_transformer) which is similar to matpower test cases.). So could you please suggest me what necessary steps I should follow? Thank you for your time. Regards Mirish Thakur KIT, University. On Mon, Aug 10, 2015 at 7:14 PM, Abhyankar, Shrirang G. abhy...@anl.govmailto:abhy...@anl.gov wrote: I would

Re: PV curve using CPF

2015-08-13 Thread Abhyankar, Shrirang G.
but that is for that bus only for which p-v curve is plotted. My question is How to find Voltage Stability Margin for whole Network using P-V curve ? I mean how to find maximum lamda for whole network using p-v curve? Thanks. From: Abhyankar, Shrirang G. abhy...@anl.govmailto:abhy...@anl.gov Sent: Fri

Re: convergence problem in runpf.

2015-08-11 Thread Abhyankar, Shrirang G.
, Shrirang G. abhy...@anl.govmailto:abhy...@anl.gov wrote: I would suggest trying the following: 1. Use the solution of a fast decoupled power flow or an optimal power flow (with line limits and voltage limits relaxed) as the initial guess for the power flow. 2. Follow step 5 in http

Re: convergence problem in runpf.

2015-08-10 Thread Abhyankar, Shrirang G.
I would suggest trying the following: 1. Use the solution of a fast decoupled power flow or an optimal power flow (with line limits and voltage limits relaxed) as the initial guess for the power flow. 2. Follow step 5 in http://www.pserc.cornell.edu/matpower/#pfconvergence making CPF to

Re: PF Not converging

2015-08-08 Thread Abhyankar, Shrirang G.
See Faq#5 Why does MATPOWER power flow not converge?http://www.pserc.cornell.edu//matpower/#pfconvergence Shri From: Syed, Waqas Ahmad Raazi waqas.s...@rwth-aachen.demailto:waqas.s...@rwth-aachen.de Reply-To: MATPOWER discussion forum

Re: PV curve using CPF

2015-08-07 Thread Abhyankar, Shrirang G.
I don’t quite understand your question, can you please elaborate. The maximum value of loading scaling parameter ‘lambda’ gives a measure of how much power can be transferred for a given transfer direction. So, lambda is also a measure of the nose point for the whole network. Shri From:

Re: LMP calculation using MATPOWER

2015-07-31 Thread Abhyankar, Shrirang G.
I don’t think there exist built-in MATPOWER functions to calculate loss and delivery factors. However, you can easily compute them, and possibly think about contributing, using the available MATPOWER functions. The full list of available functions can be found here

Re: Questions about dynamic model parameter of IEEE 39bus system

2015-07-27 Thread Abhyankar, Shrirang G.
You'll have to contact MatDYN developers to figure out what's wrong. Shri On Jul 25, 2015, at 6:22 AM, 陈春宇 chunyu...@gmail.com wrote: Hi everyone, I try to build the dynamic model of IEEE 39 bus in MatDyn, however, everytime I run the simulation, it reports errors which say the

Re: Equvalence of IEEE 118 bus system

2015-07-21 Thread Abhyankar, Shrirang G.
The latest version of MATPOWER has code for creating reduced equivalents of larger networks. The code is in the ‘reduction’ folder in the extras directory. See page 148 of the MATPOWER manual. Shri From: Nandu Balachandran nbala...@my.uno.edumailto:nbala...@my.uno.edu Reply-To: MATPOWER

Re: Problem in finding inverse

2015-06-26 Thread Abhyankar, Shrirang G.
Niharika, There could be multitude of reasons causing matrix singularity for power flow or optimal power flow problems. See FAQ 5 http://www.pserc.cornell.edu//matpower/#pfconvergence for some helpful tips on debugging this issue. Shri From: Niharika Iter

Re: Jacobian matrix values

2015-06-03 Thread Abhyankar, Shrirang G.
. Shri Hiranya Sent from my iPhone On 3 Jun 2015, at 02:26, Abhyankar, Shrirang G. abhy...@anl.gov wrote: The Jacobian is stored in a sparse matrix format. The Jacobian values not displayed are all zeros. MATPOWER orders the Jacobian as follows: J = [dP_dVa dP_dVm dQ_dVa

Re: Jacobian matrix values

2015-06-02 Thread Abhyankar, Shrirang G.
could please help me out as my project deadline is next week. Thank you Hiranya On 2 Jun 2015, at 17:08, Abhyankar, Shrirang G. abhy...@anl.gov wrote: fullJac is a flag that you can set to either access the full Jacobian matrix or a partial Jacobian. With fullJac = 1, the routine returns

Re: Jacobian matrix values

2015-06-02 Thread Abhyankar, Shrirang G.
http://www.pserc.cornell.edu//matpower/docs/ref/matpower5.1/makeJac.html -Original Message- From: Hiranya Seneviratne hiranya.seneviratn...@gmail.com Reply-To: MATPOWER discussion forum matpowe...@list.cornell.edu Date: Tuesday, June 2, 2015 at 4:04 AM To: MATPOWER discussion forum

Re: Power factor

2015-06-02 Thread Abhyankar, Shrirang G.
Before running the power flow, you can change the load P and Q values in the bus data. Shri -Original Message- From: Hiranya Seneviratne hiranya.seneviratn...@gmail.com Reply-To: MATPOWER discussion forum matpowe...@list.cornell.edu Date: Tuesday, June 2, 2015 at 6:17 AM To: MATPOWER

Re: Power factor

2015-06-02 Thread Abhyankar, Shrirang G.
:31, Abhyankar, Shrirang G. abhy...@anl.gov wrote: Before running the power flow, you can change the load P and Q values in the bus data. Shri -Original Message- From: Hiranya Seneviratne hiranya.seneviratn...@gmail.com Reply-To: MATPOWER discussion forum matpowe

Re: Jacobian matrix values

2015-06-02 Thread Abhyankar, Shrirang G.
Sent from my iPhone On 2 Jun 2015, at 16:29, Abhyankar, Shrirang G. abhy...@anl.gov wrote: http://www.pserc.cornell.edu//matpower/docs/ref/matpower5.1/makeJac.html -Original Message- From: Hiranya Seneviratne hiranya.seneviratn...@gmail.com Reply-To: MATPOWER discussion forum

  1   2   3   >