Dear Dr Zimmerman, Thanks for your reply. I would appreciate it if you let me know how to find the code for running Implicit Zbus Gauss method. Regards Zeinab Ghofrani
________________________________ From: Ray Zimmerman <[email protected]> To: MATPOWER discussion forum <[email protected]> Sent: Friday, October 12, 2012 7:28 AM Subject: Re: convergence of 32-bus distribution network PF That is not the Implicit Zbus Gauss method described in my dissertation. That is the standard Gauss-Seidel method that is already implemented in MATPOWER. And by the way, the easy way to run that one is … opt = mpoption('PF_ALG', 4); runpf(case_name, opt); -- Ray Zimmerman Senior Research Associate 419A Warren Hall, Cornell University, Ithaca, NY 14853 phone: (607) 255-9645 On Oct 10, 2012, at 3:26 AM, Zeinab Ghofrani <[email protected]> wrote: Dear Dr Zimmerman > >I implemented the Gauss method using the following codes, but still there is >convergence problem: > > >case_name='case_dist_32bus';mpc = loadcase(case_name); > >ref=[33]; pv=[]; >pq=1:32; pq=pq'; >V0=ones(size(mpc.bus,1),1);Sbus=makeSbus(mpc.baseMVA, mpc.bus, mpc.gen); >Ybus=makeYbus(mpc);results=gausspf(Ybus, Sbus, V0, ref, pv, pq); > >The system matpower data is attached to this email. > >Thanks >Zeinab Ghofrani > > >________________________________ >From: Ray Zimmerman <[email protected]> >To: MATPOWER discussion forum <[email protected]> >Sent: Monday, October 8, 2012 7:49 AM >Subject: Re: convergence of 32-bus distribution network PF > > >I wish I could. Unfortunately, the code I used for my thesis was tightly tied >to some proprietary data formats and it also does not run on current versions >of Matlab. I you already have the data in MATPOWER format, it would probably >be rather straightforward to implement the Implicit Zbus Gauss method, since >MATPOWER can easily build the Ybus matrix for you. That's probably the >approach I would take if I were you. > > >-- >Ray Zimmerman >Senior Research Associate >419A Warren Hall, Cornell University, Ithaca, NY 14853 >phone: (607) 255-9645 > > > > >On Oct 6, 2012, at 1:27 AM, Zeinab Ghofrani <[email protected]> wrote: > >Dear Dr Zimmerman; >> >>Thanks for sending me your thesis, which was very useful. Could you please >>provide me with a reliable MATLAB code for running backward-forward sweep >>load flow on a radial distribution matpower case? >> >> >> >> >>Zeinab Ghofrani >> >> >> >>________________________________ >>From: Ray Zimmerman <[email protected]> >>To: MATPOWER discussion forum <[email protected]> >>Sent: Friday, October 5, 2012 6:21 AM >>Subject: Re: convergence of 32-bus distribution network PF >> >> >>My guess is that it may require a different power flow formulation, namely >>one that is specifically designed for radial distribution systems. Feel free >>to look at my dissertation and the references in it for some ideas. >> >> >>-- >>Ray Zimmerman >>Senior Research Associate >>419A Warren Hall, Cornell University, Ithaca, NY 14853 >>phone: (607) 255-9645 >> >> >> >> >>On Oct 5, 2012, at 2:37 AM, Zeinab Ghofrani <[email protected]> wrote: >> >>Dear all; >>> >>>I'm tryning to run power flow on the 32 bus distribution system described in >>>the following paper: >>>M. E. Baran and F. F. Wu, “Network Reconfiguration in Distribution Systems >>>for Loss Reduction and Load Balancing”, IEEE Trans on Power Delivery, 1989, >>>Vol 4(2), pp >>>1401-1407. >>>I have chnaged the place of the tie switches so that the modified system has >>>a long latteral with high voltage drop. The power flow doesn't converge, >>>unless we set Vg of mpc.gen matrix to 1.3. How can I have a converged load >>>flow when Vg is set to 1? >>>The data I use is attached to this email. Waiting for your reply. >>> >>>Zeinab Ghofrani<case_dist_32bus.m> >> >> >> > > ><case_dist_32bus.m>
