Hi Shri I tried using a smaller continuation step size as your suggestion (step size=0.05). And I tried decreasing step size to 0.002, I still got the same result, the PV curve is not full nose curve.
Best Regards, Natakorn Thasnas Ph.D Student Khonkaen University, Thailand 2015-02-12 20:41 GMT+07:00 Abhyankar, Shrirang G. <[email protected]>: > Try using a smaller continuation step size. > mpopt = mpoption(mpopt, 'cpf.stop_at', 'full', 'cpf.step', 0.05); > > Shri > > > On Feb 12, 2015, at 4:03 AM, "Natakorn Thasnas" <[email protected]> > wrote: > > > > Hi All > > > > I am trying to run a continuation power flow with the following code > > > > define_constants; > > mpopt = mpoption('out.all', 0, 'verbose', 2); > > mpopt = mpoption(mpopt, 'cpf.stop_at', 'full', 'cpf.step', 0.2); > > mpopt = mpoption(mpopt, 'cpf.plot.level', 2 ,'cpf.plot.bus', 14); > > mpcb = loadcase('case14result'); % load base case > > mpct = mpcb; % set up target case with > > mpct.gen(:, [PG QG]) = mpcb.gen(:, [PG QG])*2.5; > > mpct.bus(:,[PD QG]) = mpcb.bus(:,[PD QG])*2.5; > > results = runcpf(mpcb, mpct, mpopt); > > > > After running it, I got the warning message "corrector did not converge > in 1000 iterations" and the PV curve is not full nose curve as showing in > attached file (c14.jpg) > > > > <c14.jpg> > > > > Then I increased "max_it" ​from 1000 to 10000. I got the warning > messages as following > > > > Warning: Matrix is close to singular or badly scaled. > > Results may be inaccurate. RCOND = 1.944067e-018. > > > In cpf_corrector at 143 > > In runcpf at 268 > > In cpfrun at 55 > > step 48 : lambda = -7.770, corrector did not converge in 10000 > iterations > > > > > > Could you please suggest me how to fix it to get "the full nose curve"? > > > > Thank you in advance. > > > > Best Regards, > > > > Natakorn Thasnas > > Ph.D Student > > Khonkaen University, Thailand > >
