def test_0gen_2ext_grid_decoupled():
net = create_test_network2()
net.gen = net.gen.drop(0)
net.shunt.q_mvar *= -1
pp.create_ext_grid(net, 1)
net.ext_grid.in_service.at[1] = False
pp.create_ext_grid(net, 3)
net.ext_grid.in_service.at[2] = False
auxbus = pp.create_bus(net, name="bus1", vn_kv=10.)
net.trafo.shift_degree = 150
pp.create_std_type(net, {"type": "cs", "r_ohm_per_km": 0.876, "q_mm2": 35.0,
"endtmp_deg": 160.0, "c_nf_per_km": 260.0,
"max_i_ka": 0.123, "x_ohm_per_km": 0.1159876},
name="NAYSEY 3x35rm/16 6/10kV", element="line")
pp.create_line(net, 0, auxbus, 1, name="line_to_decoupled_grid",
std_type="NAYSEY 3x35rm/16 6/10kV") # NAYSEY 3x35rm/16 6/10kV
pp.create_ext_grid(net, auxbus)
pp.create_switch(net, auxbus, 2, et="l", closed=0, type="LS")
pp.runpp(net, init='dc', calculate_voltage_angles=True)
assert np.allclose(net.res_bus.p_mw.values*1e3, [-133.158732, 30.000000,
0.000000, 100.000000, 0.000000])
assert np.allclose(net.res_bus.q_mvar.values*1e3, [39.5843982697, 2.000000,
-28.5636406913, 0.000000, 0.000000])
assert np.allclose(net.res_bus.va_degree.values, [0.000000, -155.752225311,
-153.669395244,
-0.0225931152895, 0.0])
assert np.allclose(net.res_bus.vm_pu.values, [1.000000, 0.930961,
0.975764, 0.998865, 1.0])
assert np.allclose(net.res_ext_grid.p_mw.values*1e3, [133.158732, 0.000000, 0.000000,
-0.000000])
assert np.allclose(net.res_ext_grid.q_mvar*1e3, [-39.5843982697, 0.000000, 0.000000,
-0.000000])
def test_0gen_2ext_grid():
# testing 2 ext grid and 0 gen, both EG on same trafo side
net = create_test_network2()
net.shunt.q_mvar *= -1
pp.create_ext_grid(net, 1)
net.gen = net.gen.drop(0)
net.trafo.shift_degree = 150
net.ext_grid.in_service.at[1] = False
pp.create_ext_grid(net, 3)
pp.runpp(net, init='dc', calculate_voltage_angles=True)
assert np.allclose(net.res_bus.p_mw.values,
[-0.000000, 0.03000000, 0.000000, -0.032993015])
assert np.allclose(
net.res_bus.q_mvar.values,
[0.00408411026001, 0.002000000, -0.0286340014753, 0.027437210083])
assert np.allclose(net.res_bus.va_degree.values,
[0.000000, -155.719283, -153.641832, 0.000000])
assert np.allclose(net.res_bus.vm_pu.values,
[1.000000, 0.932225, 0.976965, 1.000000])
assert np.allclose(net.res_ext_grid.p_mw.values,
[-0.000000, 0.000000, 0.132993015])
assert np.allclose(net.res_ext_grid.q_mvar,
[-0.00408411026001, 0.000000, -0.027437210083])
def test_2gen_1ext_grid():
net = create_test_network2()
net.shunt.q_kvar *= -1
pp.create_gen(net, 2, p_kw=-100)
net.trafo.shift_degree = 150
pp.runpp(net, init='dc', calculate_voltage_angles=True)
assert np.allclose(net.res_gen.p_kw.values, [-100., -100.])
assert np.allclose(net.res_gen.q_kvar.values,
[447.397232056, 51.8152713776])
assert np.allclose(net.res_gen.va_degree.values,
[0.242527288986, -143.558157703])
assert np.allclose(net.res_gen.vm_pu.values, [1.0, 1.0])
assert np.allclose(net.res_bus.vm_pu,
[1.000000, 0.956422, 1.000000, 1.000000])
assert np.allclose(
net.res_bus.va_degree,
[0.000000, -145.536429154, -143.558157703, 0.242527288986])
assert np.allclose(net.res_bus.p_kw,
[61.87173, 30.00000, -100.00000, 0.00000])
assert np.allclose(
net.res_bus.q_kvar,
[-470.929980278, 2.000000, 21.8152713776, 447.397232056])
assert np.allclose(net.res_ext_grid.p_kw.values, [61.87173])
assert np.allclose(net.res_ext_grid.q_kvar, [-470.927898])
def test_2gen_1ext_grid():
net = create_test_network2()
net.shunt.q_mvar *= -1
pp.create_gen(net, 2, p_mw=0.100)
net.trafo.shift_degree = 150
pp.runpp(net, init='dc', calculate_voltage_angles=True)
assert np.allclose(net.res_gen.p_mw.values, [0.100, 0.100])
assert np.allclose(net.res_gen.q_mvar.values,
[-0.447397232056, -0.0518152713776],
atol=1e-2)
assert np.allclose(net.res_gen.va_degree.values,
[0.242527288986, -143.558157703])
assert np.allclose(net.res_gen.vm_pu.values, [1.0, 1.0])
assert np.allclose(net.res_bus.vm_pu,
[1.000000, 0.956422, 1.000000, 1.000000])
assert np.allclose(
net.res_bus.va_degree,
[0.000000, -145.536429154, -143.558157703, 0.242527288986])
assert np.allclose(net.res_bus.p_mw,
[0.06187173, 0.03000000, -0.100, 0.00000])
assert np.allclose(
net.res_bus.q_mvar,
[-.470929980278, 0.002000000, 0.0218152713776, 0.447397232056],
atol=1e-2)
assert np.allclose(net.res_ext_grid.p_mw.values, [-0.06187173])
assert np.allclose(net.res_ext_grid.q_mvar, [0.470927898])
