def test_const_pq_out_of_service(simple_test_net):
# allows to use recycle = {"bus_pq"} and fast output read
net = simple_test_net
for i in range(3):
b = pp.create_bus(net, 20., in_service=False)
pp.create_line(net,
2,
b,
std_type="149-AL1/24-ST1A 20.0",
length_km=1.,
in_service=False)
pp.create_transformer(net,
2,
b,
std_type="25 MVA 110/20 kV",
in_service=False)
pp.create_transformer3w(net,
1,
2,
b,
std_type="63/25/38 MVA 110/20/10 kV",
in_service=False)
_, ds = create_data_source(n_timesteps)
# 1load
c = add_const(net, ds, recycle=None)
vm_pu, ll = _run_recycle(net)
del c
# calculate the same results without recycle
c = add_const(net, ds, recycle=False)
ow = _run_normal(net)
in_service = net.bus.loc[net.bus.in_service].index
assert np.allclose(vm_pu.loc[:, in_service],
ow.output["res_bus.vm_pu"].loc[:, in_service])
in_service = net.line.loc[net.line.in_service].index
assert np.allclose(
ll.loc[:, in_service],
ow.output["res_line.loading_percent"].loc[:, in_service])
def test_const_pq_gen_trafo_tap(simple_test_net):
# allows to use recycle = {"bus_pq", "gen", "trafo"}
net = simple_test_net
profiles, _ = create_data_source(n_timesteps)
profiles['ext_grid'] = np.ones(n_timesteps) + np.arange(0, n_timesteps) * 1e-2
ds = DFData(profiles)
# 1load
c1 = add_const(net, ds, recycle=None)
c2 = add_const(net, ds, recycle=None, profile_name="ext_grid", variable="vm_pu", element_index=0,
element="ext_grid")
c3 = ContinuousTapControl(net, 0, 1.01, tol=1e-9)
vm_pu, ll = _run_recycle(net)
del c1, c2, c3
# calculate the same results without recycle
c = add_const(net, ds, recycle=False)
c2 = add_const(net, ds, recycle=False, profile_name="ext_grid", variable="vm_pu", element_index=0,
element="ext_grid")
c3 = ContinuousTapControl(net, 0, 1.01, recycle=False, tol=1e-9)
ow = _run_normal(net)
assert np.allclose(vm_pu, ow.output["res_bus.vm_pu"])
assert np.allclose(ll, ow.output["res_line.loading_percent"])
def test_const_ext_grid(simple_test_net):
# allows to use recycle = {"gen"} and fast output read
net = simple_test_net
profiles, _ = create_data_source(n_timesteps)
profiles['ext_grid'] = np.ones(n_timesteps) + np.arange(0,
n_timesteps) * 1e-2
ds = DFData(profiles)
# 1load
c1 = add_const(net, ds, recycle=None)
c2 = add_const(net,
ds,
recycle=None,
profile_name="ext_grid",
variable="vm_pu",
element_index=0,
element="ext_grid")
vm_pu, ll = _run_recycle(net)
del c1, c2
# calculate the same results without recycle
c = add_const(net, ds, recycle=False)
c2 = add_const(net,
ds,
recycle=False,
profile_name="ext_grid",
variable="vm_pu",
element_index=0,
element="ext_grid")
ow = _run_normal(net)
assert np.allclose(vm_pu, ow.output["res_bus.vm_pu"])
assert np.allclose(ll, ow.output["res_line.loading_percent"])
assert np.allclose(ow.output["res_bus.vm_pu"].values[:, 0],
profiles["ext_grid"])
assert np.allclose(vm_pu.values[:, 0], profiles["ext_grid"])
def test_batch_output_reader(simple_test_net):
net = simple_test_net
_, ds = create_data_source(n_timesteps)
# 1load
c = add_const(net, ds, recycle=True)
ow = OutputWriter(net,
output_path=tempfile.gettempdir(),
output_file_type=".json",
log_variables=list())
ow.log_variable('ppc_bus', 'vm')
ow.log_variable('ppc_bus', 'va')
recycle = dict(trafo=False, gen=False, bus_pq=True)
only_v = True
pp.runpp(net, only_v_results=only_v, recycle=recycle)
run_timeseries(net,
time_steps,
recycle=recycle,
only_v_results=only_v,
verbose=False)
vm, va = ow.output["ppc_bus.vm"], ow.output["ppc_bus.va"]
# ppc was removed (set to None) in cleanup() after the timeseries was over. now we have to remake the net._ppc by running the power flow again
pp.runpp(net, only_v_results=only_v, recycle=recycle)
s, s_abs, i_abs = v_to_i_s(net, vm, va)
del net, ow, c
net = simple_test_net
net.output_writer.drop(index=net.output_writer.index, inplace=True)
net.controller.drop(index=net.controller.index, inplace=True)
add_const(net, ds, recycle=False)
ow = OutputWriter(net,
output_path=tempfile.gettempdir(),
output_file_type=".json",
log_variables=list())
ow.log_variable('res_line', 'i_from_ka')
ow.log_variable('res_line', 'loading_percent')
ow.log_variable('res_line', 'i_to_ka')
ow.log_variable('res_line', 'loading_percent')
ow.log_variable('res_trafo', 'loading_percent')
ow.log_variable('res_trafo', 'i_hv_ka')
ow.log_variable('res_trafo', 'i_lv_ka')
ow.log_variable('res_trafo3w', 'loading_percent')
ow.log_variable('res_trafo3w', 'i_hv_ka')
ow.log_variable('res_trafo3w', 'i_mv_ka')
ow.log_variable('res_trafo3w', 'i_lv_ka')
ow.log_variable('res_bus', 'vm_pu')
ow.log_variable('res_bus', 'va_degree')
pp.runpp(net)
run_timeseries(net, time_steps, trafo_loading="current", verbose=False)
vm_pu, va_degree = ow.output["res_bus.vm_pu"], ow.output[
"res_bus.va_degree"]
i_from_ka_normal, i_to_ka_normal = ow.output[
"res_line.i_from_ka"], ow.output["res_line.i_to_ka"]
i_hv_ka_normal, i_lv_ka_normal = ow.output["res_trafo.i_hv_ka"], ow.output[
"res_trafo.i_lv_ka"]
i3_hv_ka, i3_lm_ka, i3_lv_ka = ow.output["res_trafo3w.i_hv_ka"], ow.output["res_trafo3w.i_mv_ka"], \
ow.output["res_trafo3w.i_lv_ka"]
t3_loading_percent_normal = ow.output["res_trafo3w.loading_percent"]
loading_percent_normal = ow.output["res_line.loading_percent"]
t_loading_percent_normal = ow.output["res_trafo.loading_percent"]
v_normal = polar_to_rad(vm_pu, va_degree)
v_batch = polar_to_rad(vm, va)
i_ka, i_from_ka, i_to_ka, loading_percent = get_batch_line_results(
net, i_abs)
i_ka, i_hv_ka, i_lv_ka, s_mva, ld_trafo = get_batch_trafo_results(
net, i_abs, s_abs)
i_h, i_m, i_l, ld3_trafo = get_batch_trafo3w_results(net, i_abs, s_abs)
# v batch contains aux buses of trafo3w. we can only compare non aux bus voltages
assert np.allclose(v_normal, v_batch[:, :v_normal.shape[1]])
# line loading
assert np.allclose(i_from_ka_normal, i_from_ka)
assert np.allclose(i_to_ka_normal, i_to_ka)
assert np.allclose(loading_percent_normal, loading_percent)
# trafo
assert np.allclose(i_hv_ka_normal, i_hv_ka)
assert np.allclose(i_lv_ka_normal, i_lv_ka)
assert np.allclose(t_loading_percent_normal, ld_trafo)
# trafo3w
assert np.allclose(i3_hv_ka, i_h)
assert np.allclose(i3_lm_ka, i_m)
assert np.allclose(i3_lv_ka, i_l)
assert np.allclose(t3_loading_percent_normal, ld3_trafo)