The numerical divergence could also possibly be due to the infeasibility of the problem. Can you try relaxing some of the constraints (especially the voltage magnitude constraints) and see if the OPF converges?
Switching to another solver is also an option. See the list of available optimization packages that MATPOWER supports here http://www.pserc.cornell.edu//matpower/#optionalpackages These archive threads may be relevant for your problem. http://www.mail-archive.com/[email protected]/msg02915.html http://www.mail-archive.com/[email protected]/msg01298.html Shri From: Angelina <[email protected]<mailto:[email protected]>> Reply-To: MATPOWER discussion forum <[email protected]<mailto:[email protected]>> Date: Wed, 29 Jan 2014 15:35:23 +0000 To: MATPOWER discussion forum <[email protected]<mailto:[email protected]>> Subject: numerically failed opf Dear Dr Zimmerman, I am trying to solve an ac opf with the mips solver (560) in a lightly meshed 6.6 kV distribution system of 32 branches. The range for the reactances is [0.05-1.9 pu], whereas some of them have almost negligible values. When I am solving the simple ac power flow the problem converges with no problems to be appeared. However, when I am trying to solve the same problem with optimal power flow, assigning some costs to the generators, the ac opf fails numerically or the message ‘ the matrix is singular to working precision, the results may be inaccurate’ appears. I suspect that the problem is the solver. Could you recommend me another solver to use? Any guidance would be appreciated. Kind Regards, Angelina
