search

sogno-platform / dpsim

Real-time power system simulator including powerflow, (dynamic) phasors and EMT
https://sogno.energy/dpsim/
Mozilla Public License 2.0
70 stars 51 forks source link

Use smaller snubber components in PowerTransformer #192

Open martinmoraga opened 1 year ago

martinmoraga commented 1 year ago

Running powerflow to initialize the dynamic simulation of the IEEE14 network I could observe that the use of snubber components leads to an accuracy decreceament. The snubber components are calculated in such a way that they consume only 0.1% of the base power. In IEEE14 the base power is equal to 100MW and then is the power consumption of the snubber resistors and capacitors is not insignificantly. As we can see the difference between PSAT (and matpower) and DPsim is not to small:

Bus Vmag [V] Vangle [deg] P [W] Q [VAr]
N1 0.000000e+00 0.000000e+00 -8.947953e+05 2.245238e+05
N2 4.627425e-07 3.960691e-06 2.775558e-09 -7.304798e+02
N3 -8.826473e-08 4.446239e-07 1.110223e-08 2.544786e+02
N4 -1.595693e+00 1.580482e-06 0.000000e+00 0.000000e+00
N5 3.644548e-05 4.840754e-06 0.000000e+00 0.000000e+00
N6 -4.451754e-07 2.406449e-06 -1.387779e-09 -1.110082e+03
N7 5.017682e-05 -2.275554e-07 0.000000e+00 0.000000e+00
N8 3.519465e-07 -2.275805e-07 0.000000e+00 -1.914669e+02
N9 7.713768e-01 -2.769656e-06 0.000000e+00 1.375890e+03
N10 -3.631533e-01 9.733956e-08 0.000000e+00 -6.938894e-10
N11 2.615480e-01 -3.855694e-06 0.000000e+00 0.000000e+00
N12 -2.757599e-01 9.607901e-07 0.000000e+00 0.000000e+00
N13 -5.834347e-01 3.610134e-06 0.000000e+00 -6.938894e-10
N14 6.442746e-06 2.247867e-06 0.000000e+00 0.000000e+00

Using $R=10^{12}$ for the snubber resistors and $C=10^{-12}$ for the snubber capacitors, the accuracy can be considerately improved:

Bus Vmag [V] Vangle [deg] P [W] Q [VAr]
N1 0.000000e+00 0.000000e+00 4.084656e+03 -1.959172e+02
N2 4.627425e-07 3.960691e-06 2.775558e-09 -7.304798e+02
N3 -8.826473e-08 4.446239e-07 1.110223e-08 2.544786e+02
N4 -1.595693e+00 1.580482e-06 0.000000e+00 0.000000e+00
N5 3.644548e-05 4.840754e-06 0.000000e+00 0.000000e+00
N6 -4.451754e-07 2.406449e-06 -1.387779e-09 -1.110082e+03
N7 5.017682e-05 -2.275554e-07 0.000000e+00 0.000000e+00
N8 3.519465e-07 -2.275805e-07 0.000000e+00 -1.914669e+02
N9 7.713768e-01 -2.769656e-06 0.000000e+00 1.375890e+03
N10 -3.631533e-01 9.733956e-08 0.000000e+00 -6.938894e-10
N11 2.615480e-01 -3.855694e-06 0.000000e+00 0.000000e+00
N12 -2.757599e-01 9.607901e-07 0.000000e+00 0.000000e+00
N13 -5.834347e-01 3.610134e-06 0.000000e+00 -6.938894e-10
N14 6.442746e-06 2.247867e-06 0.000000e+00 0.000000e+00
martinmoraga commented 1 year ago

another option would be try a new implementation of the PT without using snubber components. I think an implementation without virtual nodes would not need snubber components and we could also improve the computation time of the ieee14 bus (and bigger networks).