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() function. Besides the
maximum loadability calculation, CPF includes several features including
step-control, real and reactive power limits, voltage and flow limits,
user-defined callbacks, termination criteria, and others. For more info, see
section 5 of the MATPOWER manual.
You can run a continuation power flow to find the maximum loading limit on the
case33bw system with the following code.
mpcbase = runpf('case33bw');
mpctarget = mpcbase;
define_constants;
mpctarget.bus(:,PD) = mpctarget.bus(:,PD)*2;
mpctarget.bus(:,QD) = mpctarget.bus(:,QD)*2;
mpopt = mpoption;
mpopt.cpf.plot.level = 2;
results = runcpf(mpcbase,mpctarget,mpopt);
Thanks,
Shri
Note that CPF uses Newton’s method, not forward-backward sweep.
From: <[email protected]> on behalf of akhil
<[email protected]>
Reply-To: MATPOWER discussion forum <[email protected]>
Date: Thursday, June 13, 2019 at 9:36 AM
To: "[email protected]" <[email protected]>
Subject: Fwd: Maximum load limit for radial distribution cases.
---------- Forwarded message ---------
From: akhil <[email protected]<mailto:[email protected]>>
Date: Thu, 13 Jun 2019 at 19:52
Subject: Maximum load limit for radial distribution cases.
To: <[email protected]<mailto:[email protected]>>
I am running the radial distribution power flow which is based on backward
forward sweep algorithm. for the network data which is provided with matpower
the radial load flow converges successfully. But to get the maximum load limit
of the system, when I am running the load flow with increased load at each bus
in a step of 1 % the radial load flow does not converge after some point.
The change in Load is as per Pnew=Lambda*P0 and Qnew=Lambda*Q0 with 1 % change
in load at all buses simultaneously.
I am Trying to get the minimum bus voltage of the system after each change in
load and running corresponding radial flow successively until the minimum
voltage of the system (bus with minimum voltage) getting reduced smoothly.
The code which I run is attached herewith for your reference.
Thank you
Yours
Akhilesh K Barnwal
Ph.D. Research Scholar
IIT BHU Varanasi.
mail: [email protected]<mailto:[email protected]>
--
@khil
