You have not included the reactive injections from line charging capacitance.
For WB5 …
–7.0 (Q gen) – 50.0 (Q load) – 35.3 (Q series losses) + 92.3 (Q branch
charging injection) = 0
--Ray
> On Apr 30, 2017, at 8:56 AM, piero reis <[email protected]> wrote:
>
> Dear MATPOWER developers,
>
> My name is Piero Carlo dos Reis and I am the MAVT student who has
> previously asked you for material on the courses in electrical engineering
> you hold at ETH.
> Since my current semester project involves an optimization code for AC
> newton Raphson procedure, I have been asked by my professor to validate it. I
> have followed the procedure illustrated by Prof. Shahidehpour in his paper
> "Security-Constrained Unit Commitment With AC Constraints" in order to model
> the Jacobian matrix (in which I also included transformer variables such as
> tap and transformer angle). A validation method advised to me by my Professor
> and which I have tried to pursue is the use of AC opf with matpower to solve
> a case study from internet (and to try to get same results with my code). I
> have download matpower version 6. However, when solving case studies (such as
> case WB5 which I am annexing together with a text file of the matpower
> solution), I get that the sum of Q generation is different from the sum of Q
> load + sum of Q losses (in WB5 matpower ac opf solution, it is evident due to
> the fact that the result is sum Q generation is -7 p.u., sum Q load is 50
> p.u. and sum Q losses is 35.30 p.u.). Would you have an explanation to this
> results and to how matpower solves the Q nodal balance of power/Q mismatch?
> In another case study I choose to run, case5, I got the also that this
> equality was not true, although the difference is not as evident as in WB5. I
> am also annexing the results from my code and from matpower in text files (I
> evaluated the losses as sum of Gen - sum of load for simplicity.)
>
> I thank you for your attention, sincerely, Piero carlo dos Reis
> <WB5.m><sol_WB5_matpower.txt><sol_case5_matpower.txt><my results-case
> 5.txt><case5.m>
