Bus 1 is the slack bus (1.0∠0∘) and bus 2 is a load (PQ) bus with P2sp=−1.0 pu and Q2sp=−0.5 pu. Unknowns: δ2 and ∣V2∣. Flat start: ∣V2∣0=1.0, δ20=0.
Y-bus
y12=1/j0.1=−j10 pu, so
Ybus=[−j10.0000j10.0000j10.0000−j10.0000]=[10∠−90∘10∠90∘10∠90∘10∠−90∘]
Power equations at bus 2
P2Q2=∣V2∣∣V1∣∣Y21∣cos(θ21−δ2+δ1)+∣V2∣2∣Y22∣cosθ22=10∣V2∣cos(90∘−δ2)=10∣V2∣sinδ2=−∣V2∣∣V1∣∣Y21∣sin(θ21−δ2+δ1)−∣V2∣2∣Y22∣sinθ22=−10∣V2∣cosδ2+10∣V2∣2
Initial power mismatch
At ∣V2∣=1, δ2=0: P20=10(1)sin0=0.0000 and Q20=−10+10=0.0000.
ΔP2ΔQ2=P2sp−P20=−1.0−0=−1.0000=Q2sp−Q20=−0.5−0=−0.5000
[ΔM0]=[ΔP2ΔQ2]=[−1.0000−0.5000]
Initial Jacobian
∂δ2∂P2∂∣V2∣∂P2∂δ2∂Q2∂∣V2∣∂Q2=10∣V2∣cosδ2=10=10sinδ2=0=10∣V2∣sinδ2=0=−10cosδ2+20∣V2∣=−10+20=10
[J0]=[10.00000.00000.000010.0000]
(For use in the first iteration: Δδ2=−1/10=−0.1000 rad =−5.730∘ and Δ∣V2∣=−0.5/10=−0.0500, giving V21=0.9500∠−5.730∘ pu.)
Answer: [ΔM0]=[−1.0000, −0.5000]T and [J0]=[100010].