Hi,

The answer below given by Ilias is correct. However, you do not have to call radial_pf directly. It is better to call it via runpf like this

opt = mpoption('pf.alg', 'PQSUM');

result = runpf('case33bw', opt);

In this case runpf makes all input data preparation and indexing. Before running any power flow solution method in runpf there is a call to ext2int which removes all isolated buses, off-line generators and branches, as well as any generators or branches connected to isolated buses.

Best regards,

Mirko

On 10/23/20 11:03 AM, Ilias Sarantakos wrote:
Hi Paulo,

Could you please share your network (case69teste1)? Because, I think it runs for (just) case69.

Thanks, and regards,

Ilias

Στις Παρ, 23 Οκτ 2020 στις 9:42 π.μ., ο/η Paulo Moisés <[email protected] <mailto:[email protected]>> έγραψε:

    Dear Ilias,

    Thank you very much for your answer to my question.
    It was very helpful but I still having some problems with my case
    study.
    In fact, after run the code you send me, I obtain the following
    errors.
    Can you, please, help me in understanding the problem?

    Thank you again
    Paulo

    image.png

    Ilias Sarantakos <[email protected]
    <mailto:[email protected]>> escreveu no dia sexta,
    23/10/2020 à(s) 09:02:

        Hi Paulo,

        Please see the code below:

        mpc = case33bw;
        c = mpc.branch(:,11); % 0,1 branch status
        mpc2 = mpc;
        mpc2.branch = mpc.branch(c~=0,:);
        mpopt = mpoption('pf.alg', 'PQSUM');
        result = radial_pf(mpc2,mpopt);

        I remember that in order to use radial_pf, the network should
        be radial, which means that even if there are open branches
        which create a loop, they should be removed.

        Hope this helps.

        Kind regards,

        Ilias

        *Dr Ilias Sarantakos*

        Research Associate

        School of Engineering

        Newcastle University

        Στις Παρ, 23 Οκτ 2020 στις 6:48 π.μ., ο/η Paulo Moisés
        <[email protected]
        <mailto:[email protected]>> έγραψε:

            Hi,

            I have a difficulty in running a radial power flow by
            using the radia_pf function in the matpower.
            Can you help me in this regard?

            Thanks
            Paulo

Reply via email to