It appears to me that the versions included in the 
<https://github.com/power-grid-lib/pglib-opf> PGLib-OPF 
c<https://github.com/power-grid-lib/pglib-opf>ollection are possibly older 
versions of the corresponding cases included in MATPOWER. A quick look at the 
output of case_info() for the two cases you mention reveals a rather large 
discrepancy between the total generation and total load. This is not an issue 
for an OPF, which dispatches generation to meet the load, but for a power flow, 
the specified generation must be close enough that any difference can be 
handled by the slack bus.

So, I would say that, with the generation as specified, those are simply not 
valid power flow cases.

    Ray



On May 25, 2021, at 12:44 PM, Abdelrahman Aldik 
<[email protected]<mailto:[email protected]>> wrote:

Thanks Ray however, my concern is about rte systems taken from the OPF 
taskforce collection. An example to that is the systems 
"pglib_opf_case2848_rte" and "pglib_opf_case2868_rte". Matpower always diverges 
when I try to use it for solving the AC PF problem. I would appreciate your 
feedback in this regard.

Regards

On Tue, May 18, 2021 at 12:56 PM Ray Daniel Zimmerman 
<[email protected]<mailto:[email protected]>> wrote:
All of the included RTE systems should solve without problems.

>> mpopt = mpoption('verbose', 1, 'out.all', 0);
>> runpf('case1888rte', mpopt)

MATPOWER Version 7.1, 08-Oct-2020 -- AC Power Flow (Newton)

Newton's method power flow (power balance, polar) converged in 2 iterations.
>> runpf('case1951rte', mpopt)

MATPOWER Version 7.1, 08-Oct-2020 -- AC Power Flow (Newton)

Newton's method power flow (power balance, polar) converged in 3 iterations.
>> runpf('case2848rte', mpopt)

MATPOWER Version 7.1, 08-Oct-2020 -- AC Power Flow (Newton)

Newton's method power flow (power balance, polar) converged in 2 iterations.
>> runpf('case2868rte', mpopt)

MATPOWER Version 7.1, 08-Oct-2020 -- AC Power Flow (Newton)

Newton's method power flow (power balance, polar) converged in 5 iterations.
>> runpf('case6468rte', mpopt)

MATPOWER Version 7.1, 08-Oct-2020 -- AC Power Flow (Newton)

Newton's method power flow (power balance, polar) converged in 3 iterations.
>> runpf('case6470rte', mpopt)

MATPOWER Version 7.1, 08-Oct-2020 -- AC Power Flow (Newton)

Newton's method power flow (power balance, polar) converged in 3 iterations.
>> runpf('case6495rte', mpopt)

MATPOWER Version 7.1, 08-Oct-2020 -- AC Power Flow (Newton)

Newton's method power flow (power balance, polar) converged in 2 iterations.
>> runpf('case6515rte', mpopt)

MATPOWER Version 7.1, 08-Oct-2020 -- AC Power Flow (Newton)

Newton's method power flow (power balance, polar) converged in 3 iterations.

If you get a different result, I would suggest reinstalling MATPOWER … maybe 
something got corrupted with your install.

   Ray


On May 14, 2021, at 12:54 AM, Abdelrahman Aldik 
<[email protected]<mailto:[email protected]>> wrote:

Dears

I would like to ask you if you have encountered any issues regarding the 
convergence of the AC power flow for the RTE french systems. I have tried 
different starting points with no avail. Any ideas about the reason why?

Regards

--
Abdelrahman Aldik
PhD Candidate
Department of Electrical and Computer Engineering and Centre of Urban Energy
Ryerson University
CUI-030, 44 Gerrard Street East, Toronto, ON, Canada
Email: [email protected]<mailto:[email protected]>



--
Abdelrahman Aldik
PhD Candidate
Department of Electrical and Computer Engineering and Centre of Urban Energy
Ryerson University
CUI-030, 44 Gerrard Street East, Toronto, ON, Canada
Email: [email protected]<mailto:[email protected]>


Reply via email to