def test_p_type():
"""
:return:
:rtype:
"""
net = pandapipes.create_empty_network("net")
d = 75e-3
pandapipes.create_junction(net, pn_bar=5, tfluid_k=293.15)
pandapipes.create_junction(net, pn_bar=5, tfluid_k=293.15)
pandapipes.create_pipe_from_parameters(net,
0,
1,
10,
diameter_m=d,
k_mm=0.1,
sections=1)
pandapipes.create_ext_grid(net, 0, p_bar=5, t_k=285.15, type="p")
pandapipes.create_sink(net, 1, mdot_kg_per_s=1)
pandapipes.create_fluid_from_lib(net, name="water")
pandapipes.pipeflow(net,
stop_condition="tol",
iter=70,
friction_model="nikuradse",
transient=False,
nonlinear_method="automatic",
tol_p=1e-4,
tol_v=1e-4)
data = pd.read_csv(os.path.join(internals_data_path, "ext_grid_p.csv"),
sep=';',
header=0,
keep_default_na=False)
p_comp = data["p"]
p_pandapipes = net.res_junction["p_bar"][0]
p_diff = np.abs(1 - p_pandapipes / p_comp.loc[0])
assert np.all(p_diff < 0.01)
def test_create_pump(create_empty_net):
net = copy.deepcopy(create_empty_net)
pandapipes.create_junction(net, 1, 293, index=8, geodata=(0, 1))
pandapipes.create_junction(net, 1, 293, index=9, geodata=(2, 2))
pandapipes.create_pump(net, 8, 9, "P1", index=2)
try:
pandapipes.create_pump(net, 8, 9, "P1", index=2)
assert False, "Shouldn't make pump, index exists!"
except UserWarning:
assert True
try:
pandapipes.create_pump(net, 8, 10, "P1")
assert False, "Shouldn't make pump, non-existent to_junction!"
except UserWarning:
assert True
try:
pandapipes.create_pump(net, 8, 9, "blah")
assert False, "Shouldn't make pump, non-existent std_type!"
except UserWarning:
assert True
try:
pandapipes.create_pump(net,
8,
9,
"P1",
geodata=[(0, 1), (1, 1), (2, 2)])
assert False, "Shouldn't make pump, geodata not possible!"
except ValueError:
assert True
net2 = pandapipes.create_empty_network(fluid="hgas", add_stdtypes=False)
pandapipes.create_junction(net2, 1, 293, index=8, geodata=(0, 1))
pandapipes.create_junction(net2, 1, 293, index=9, geodata=(2, 2))
try:
pandapipes.create_pump(net2, 8, 9, "P1")
assert False, "Shouldn't make pump, no std_types in net!"
except UserWarning:
assert True
def test_fluid_exceptions():
net = pandapipes.create_empty_network(fluid="hgas")
fluid = pandapipes.get_fluid(net)
try:
fluid.get_property("xyz", 100)
assert False, "Shouldn't find property xyz!"
except UserWarning:
pass
prop = pandapipes.FluidProperty()
try:
prop.get_property(100)
assert False, "Shouldn't have property defined!"
except NotImplementedError:
pass
try:
pandapipes.call_lib("natural_gas")
assert False, "Shouldn't have fluid natural_gas defined in lib!"
except AttributeError:
pass
def test_heat_exchanger():
"""
:return:
:rtype:
"""
net = pandapipes.create_empty_network("net", add_stdtypes=False)
d = 75e-3
pandapipes.create_junction(net, pn_bar=5, tfluid_k=283)
pandapipes.create_junction(net, pn_bar=5, tfluid_k=283)
pandapipes.create_heat_exchanger(net, 0, 1, d, qext_w=20000)
pandapipes.create_ext_grid(net, 0, p_bar=5, t_k=330, type="pt")
pandapipes.create_sink(net, 1, mdot_kg_per_s=1)
pandapipes.create_fluid_from_lib(net, "water", overwrite=True)
pandapipes.pipeflow(net,
stop_condition="tol",
iter=3,
friction_model="nikuradse",
mode="all",
transient=False,
nonlinear_method="automatic",
tol_p=1e-4,
tol_v=1e-4)
data = pd.read_csv(os.path.join(internals_data_path,
"heat_exchanger_test.csv"),
sep=';',
header=0,
keep_default_na=False)
temp_an = data["T1"]
t_pan = net.res_junction.t_k
temp_diff = np.abs(1 - t_pan / temp_an)
assert np.all(temp_diff < 0.01)
def test_temperature_internal_nodes_masche_1load():
"""
:return:
:rtype:
"""
net = pandapipes.create_empty_network("net", add_stdtypes=False)
d = 75e-3
j0 = pandapipes.create_junction(net, pn_bar=5, tfluid_k=283)
j1 = pandapipes.create_junction(net, pn_bar=5, tfluid_k=283)
j2 = pandapipes.create_junction(net, pn_bar=5, tfluid_k=283)
j3 = pandapipes.create_junction(net, pn_bar=5, tfluid_k=283)
pandapipes.create_pipe_from_parameters(net, j0, j1, 2.5, d, k_mm=.1, sections=6,
alpha_w_per_m2k=5)
pandapipes.create_pipe_from_parameters(net, j1, j2, 2.5, d, k_mm=.1, sections=6,
alpha_w_per_m2k=5)
pandapipes.create_pipe_from_parameters(net, j1, j3, 2.5, d, k_mm=.1, sections=6,
alpha_w_per_m2k=5)
pandapipes.create_pipe_from_parameters(net, j3, j2, 2.5, d, k_mm=.1, sections=6,
alpha_w_per_m2k=5)
pandapipes.create_ext_grid(net, j0, p_bar=5, t_k=350, type="pt")
pandapipes.create_sink(net, j2, mdot_kg_per_s=1)
pandapipes.create_fluid_from_lib(net, "water", overwrite=True)
pandapipes.pipeflow(net, stop_condition="tol", iter=70, friction_model="nikuradse",
mode='all', transient=False, nonlinear_method="automatic", tol_p=1e-4,
tol_v=1e-4)
data = pd.read_csv(os.path.join(internals_data_path, "Temperature_masche_1load_an.csv"),
sep=';', header=0, keep_default_na=False)
temp_an = data["T"]
temp_pandapipes = net.res_junction["t_k"]
temp_diff = np.abs(1 - temp_pandapipes / temp_an)
assert np.all(temp_diff < 0.01)
def test_hydraulic_only():
"""
:return:
:rtype:
"""
net = pp.create_empty_network("net")
d = 75e-3
pp.create_junction(net, pn_bar=5, tfluid_k=283)
pp.create_junction(net, pn_bar=5, tfluid_k=283)
pp.create_pipe_from_parameters(net, 0, 1, 6, diameter_m=d, k_mm=.1, sections=1,
alpha_w_per_m2k=5)
pp.create_ext_grid(net, 0, p_bar=5, t_k=330, type="pt")
pp.create_sink(net, 1, mdot_kg_per_s=1)
pp.create_fluid_from_lib(net, "water", overwrite=True)
pp.pipeflow(net, stop_condition="tol", iter=70, friction_model="nikuradse",
transient=False, nonlinear_method="automatic", tol_p=1e-4,
tol_v=1e-4)
data = pd.read_csv(os.path.join(data_path, "hydraulics.csv"), sep=';', header=0,
keep_default_na=False)
node_pit = net["_pit"]["node"]
branch_pit = net["_pit"]["branch"]
v_an = data.loc[0,"pv"]
p_an = data.loc[1:3, "pv"]
p_pandapipes = node_pit[:, PINIT]
v_pandapipes = branch_pit[:, VINIT]
p_diff = np.abs(1 - p_pandapipes / p_an)
v_diff = np.abs(v_pandapipes-v_an)
assert np.all(p_diff < 0.01)
assert(np.all(v_diff<0.05))
def test_create_pipe(create_empty_net):
net = copy.deepcopy(create_empty_net)
pandapipes.create_junction(net, 1, 293, index=8, geodata=(0, 1))
pandapipes.create_junction(net, 1, 293, index=9, geodata=(2, 2))
pandapipes.create_pipe(net,
8,
9,
"80_GGG",
0.3,
index=2,
geodata=[(0, 1), (1, 1), (2, 2)])
with pytest.raises(UserWarning):
pandapipes.create_pipe(net, 8, 9, "80_GGG", 0.3, index=2)
with pytest.raises(UserWarning):
pandapipes.create_pipe(net, 8, 10, "80_GGG", 0.3)
with pytest.raises(UserWarning):
pandapipes.create_pipe(net, 8, 9, "blah", 0.3)
net2 = pandapipes.create_empty_network(fluid="hgas", add_stdtypes=False)
pandapipes.create_junction(net2, 1, 293, index=8, geodata=(0, 1))
pandapipes.create_junction(net2, 1, 293, index=9, geodata=(2, 2))
with pytest.raises(UserWarning):
pandapipes.create_pipe(net2, 8, 9, "80_GGG", 0.3)
def test_create_pipe(create_empty_net):
net = copy.deepcopy(create_empty_net)
pandapipes.create_junction(net, 1, 293, index=8, geodata=(0, 1))
pandapipes.create_junction(net, 1, 293, index=9, geodata=(2, 2))
pandapipes.create_pipe(net,
8,
9,
"80_GGG",
0.3,
index=2,
geodata=[(0, 1), (1, 1), (2, 2)])
try:
pandapipes.create_pipe(net, 8, 9, "80_GGG", 0.3, index=2)
assert False, "Shouldn't make pipe, index exists!"
except UserWarning:
assert True
try:
pandapipes.create_pipe(net, 8, 10, "80_GGG", 0.3)
assert False, "Shouldn't make pipe, non-existent to_junction!"
except UserWarning:
assert True
try:
pandapipes.create_pipe(net, 8, 9, "blah", 0.3)
assert False, "Shouldn't make pipe, non-existent std_type!"
except UserWarning:
assert True
net2 = pandapipes.create_empty_network(fluid="hgas", add_stdtypes=False)
pandapipes.create_junction(net2, 1, 293, index=8, geodata=(0, 1))
pandapipes.create_junction(net2, 1, 293, index=9, geodata=(2, 2))
try:
pandapipes.create_pipe(net2, 8, 9, "80_GGG", 0.3)
assert False, "Shouldn't make pipe, no std_types in net!"
except UserWarning:
assert True
def test_t_type_single_pipe():
"""
:return:
:rtype:
"""
net = pandapipes.create_empty_network("net")
d = 75e-3
j0 = pandapipes.create_junction(net, pn_bar=5, tfluid_k=283)
j1 = pandapipes.create_junction(net, pn_bar=5, tfluid_k=283)
pandapipes.create_ext_grid(net, j0, 5, 645, type="pt")
pandapipes.create_sink(net, j1, 1)
pandapipes.create_pipe_from_parameters(net,
j0,
j1,
6,
diameter_m=d,
k_mm=.1,
sections=1,
alpha_w_per_m2k=5)
pandapipes.create_fluid_from_lib(net, "water", overwrite=True)
pandapipes.pipeflow(net,
stop_condition="tol",
iter=70,
friction_model="nikuradse",
transient=False,
nonlinear_method="automatic",
tol_p=1e-4,
tol_v=1e-4,
mode="all")
temp = net.res_junction.t_k.values
net2 = pandapipes.create_empty_network("net")
d = 75e-3
j0 = pandapipes.create_junction(net2, pn_bar=5, tfluid_k=283)
j1 = pandapipes.create_junction(net2, pn_bar=5, tfluid_k=283)
pandapipes.create_ext_grid(net2, j0, 5, 645, type="p")
pandapipes.create_ext_grid(net2, j1, 100, 323.15, type="t")
pandapipes.create_sink(net2, j1, 1)
pandapipes.create_pipe_from_parameters(net2,
j0,
j1,
6,
diameter_m=d,
k_mm=.1,
sections=1,
alpha_w_per_m2k=5)
pandapipes.create_fluid_from_lib(net2, "water", overwrite=True)
pandapipes.pipeflow(net2,
stop_condition="tol",
iter=70,
friction_model="nikuradse",
transient=False,
nonlinear_method="automatic",
tol_p=1e-4,
tol_v=1e-4,
mode="all")
temp2 = net2.res_junction.t_k.values
temp_diff = np.abs(1 - temp / temp2)
assert np.all(temp_diff < 0.01)
def test_heat_only():
net = pp.create_empty_network("net")
d = 75e-3
pp.create_junction(net, pn_bar=5, tfluid_k=283)
pp.create_junction(net, pn_bar=5, tfluid_k=283)
pp.create_pipe_from_parameters(net,
0,
1,
6,
diameter_m=d,
k_mm=.1,
sections=6,
alpha_w_per_m2k=5)
pp.create_ext_grid(net, 0, p_bar=5, t_k=330, type="pt")
pp.create_sink(net, 1, mdot_kg_per_s=1)
pp.create_fluid_from_lib(net, "water", overwrite=True)
pp.pipeflow(net,
stop_condition="tol",
iter=70,
friction_model="nikuradse",
nonlinear_method="automatic",
mode="all")
ntw = pp.create_empty_network("net")
d = 75e-3
pp.create_junction(ntw, pn_bar=5, tfluid_k=283)
pp.create_junction(ntw, pn_bar=5, tfluid_k=283)
pp.create_pipe_from_parameters(ntw,
0,
1,
6,
diameter_m=d,
k_mm=.1,
sections=6,
alpha_w_per_m2k=5)
pp.create_ext_grid(ntw, 0, p_bar=5, t_k=330, type="pt")
pp.create_sink(ntw, 1, mdot_kg_per_s=1)
pp.create_fluid_from_lib(ntw, "water", overwrite=True)
pp.pipeflow(ntw,
stop_condition="tol",
iter=50,
friction_model="nikuradse",
nonlinear_method="automatic",
mode="hydraulics")
p = ntw._pit["node"][:, 5]
v = ntw._pit["branch"][:, 12]
u = np.concatenate((p, v))
pp.pipeflow(ntw,
sol_vec=u,
stop_condition="tol",
iter=50,
friction_model="nikuradse",
nonlinear_method="automatic",
mode="heat")
T_net = net.res_junction.t_k
T_ntw = ntw.res_junction.t_k
T_diff = np.abs(1 - T_net / T_ntw)
assert np.all(T_diff < 0.01)
def release_control_test_network():
# empty net
net = pp.create_empty_network("net", add_stdtypes=False)
# fluid
pp.create_fluid_from_lib(net, "water", overwrite=True)
# junctions
pp.create_junction(net, pn_bar=3, tfluid_k=293, height_m=0, name="Junction 0", index=None, in_service=True,
type="junction", geodata=None)
pp.create_junction(net, pn_bar=3, tfluid_k=293, height_m=0, name="Junction 1")
pp.create_junction(net, pn_bar=3, tfluid_k=293, height_m=0, name="Junction 2")
pp.create_junction(net, pn_bar=3, tfluid_k=293, height_m=0, name="Junction 3")
pp.create_junction(net, pn_bar=3, tfluid_k=293, height_m=0, name="Junction 4")
pp.create_junction(net, pn_bar=3, tfluid_k=293, height_m=0, name="Junction 5")
pp.create_junction(net, pn_bar=3, tfluid_k=293, height_m=0, name="Junction 6")
pp.create_junction(net, pn_bar=3, tfluid_k=293, height_m=0, name="Junction 7")
pp.create_junction(net, pn_bar=3, tfluid_k=293, height_m=0, name="Junction 8")
pp.create_junction(net, pn_bar=3, tfluid_k=293, height_m=0, name="Junction 9")
pp.create_junction(net, pn_bar=3, tfluid_k=293, height_m=0, name="Junction 10")
pp.create_junction(net, pn_bar=3, tfluid_k=293, height_m=0, name="Junction 11")
pp.create_junction(net, pn_bar=3, tfluid_k=293, height_m=0, name="Junction 12")
# pipes
pp.create_pipe_from_parameters(net, from_junction=0, to_junction=8, length_km=3, diameter_m=0.01, k_mm=1,
loss_coefficient=0, sections=10, alpha_w_per_m2k=10, text_k=293,
qext_w=0., name="Pipe 0", index=None, geodata=None, in_service=True, type="pipe")
pp.create_pipe_from_parameters(net, 9, 2, length_km=6, diameter_m=0.075, k_mm=.1, sections=10,
alpha_w_per_m2k=3,
name="Pipe 1")
pp.create_pipe_from_parameters(net, 2, 12, length_km=5, diameter_m=0.06, k_mm=.1, sections=10,
alpha_w_per_m2k=20,
name="Pipe 2")
pp.create_pipe_from_parameters(net, 4, 12, length_km=0.1, diameter_m=0.07, k_mm=.1, sections=10,
alpha_w_per_m2k=2,
name="Pipe 3")
pp.create_pipe_from_parameters(net, 5, 3, length_km=1, diameter_m=0.09, k_mm=.1, sections=10, alpha_w_per_m2k=3,
name="Pipe 4")
pp.create_pipe_from_parameters(net, 4, 11, length_km=2.5, diameter_m=0.08, k_mm=.1, sections=10,
alpha_w_per_m2k=15,
name="Pipe 5")
pp.create_pipe_from_parameters(net, 7, 6, length_km=4.5, diameter_m=0.085, k_mm=.1, sections=10,
alpha_w_per_m2k=2.5, name="Pipe 6")
pp.create_pipe_from_parameters(net, 1, 7, length_km=4, diameter_m=0.03, k_mm=.1, sections=10, alpha_w_per_m2k=1,
name="Pipe 7")
# external grids
pp.create_ext_grid(net, junction=0, p_bar=3, t_k=300, name="External Grid 0", in_service=True, index=None,
type="pt")
pp.create_ext_grid(net, 1, p_bar=5, t_k=350, name="External Grid 1", type="pt")
# sinks
pp.create_sink(net, junction=2, mdot_kg_per_s=0.2, scaling=1., name="Sink 0", index=None, in_service=True,
type="sink")
pp.create_sink(net, 3, mdot_kg_per_s=0.1, name="Sink 1")
pp.create_sink(net, 4, mdot_kg_per_s=0.5, name="Sink 2")
pp.create_sink(net, 5, mdot_kg_per_s=0.07, name="Sink 3")
pp.create_sink(net, 6, mdot_kg_per_s=0.09, name="Sink 4")
pp.create_sink(net, 7, mdot_kg_per_s=0.1, name="Sink 5")
# sources
pp.create_source(net, junction=8, mdot_kg_per_s=0.1, scaling=1., name="Source 0", index=None, in_service=True,
type="source")
pp.create_source(net, junction=9, mdot_kg_per_s=0.03, name="Source 1")
pp.create_source(net, junction=10, mdot_kg_per_s=0.04, name="Source 2")
pp.create_source(net, junction=11, mdot_kg_per_s=0.09, name="Source 3")
# valves
pp.create_valve(net, from_junction=8, to_junction=9, diameter_m=0.1, opened=True, loss_coefficient=0,
name="Valve 0", index=None, type="valve")
pp.create_valve(net, 9, 4, diameter_m=0.05, opened=True, name="Valve 1")
# pump
pp.create_pump_from_parameters(net, from_junction=8, to_junction=3, new_std_type_name="Pump",
pressure_list=[6.1, 5.8, 4],
flowrate_list=[0, 19, 83], reg_polynomial_degree=2,
poly_coefficents=None, name=None, index=None, in_service=True,
type="pump")
# circulation pump mass
pp.create_circ_pump_const_mass_flow(net, from_junction=3, to_junction=4, p_bar=6, mdot_kg_per_s=1,
t_k=290, name="Circ. Pump Mass", index=None, in_service=True,
type="pt")
# circulation pump pressure
pp.create_circ_pump_const_pressure(net, from_junction=11, to_junction=5, p_bar=5, plift_bar=2,
t_k=290, name="Circ. Pump Pressure", index=None, in_service=True, type="pt")
# heat exchanger
pp.create_heat_exchanger(net, from_junction=10, to_junction=6, diameter_m=0.08, qext_w=50, loss_coefficient=0,
name="Heat Exchanger 0", index=None, in_service=True, type="heat_exchanger")
pp.create_heat_exchanger(net, from_junction=4, to_junction=10, diameter_m=0.08, qext_w=28000,
loss_coefficient=0,
name="Heat Exchanger 1", index=None, in_service=True, type="heat_exchanger")
# time series
profiles_sink = pd.read_csv(os.path.join(pp_dir, 'test', 'api', 'release_cycle',
'release_control_test_sink_profiles.csv'), index_col=0)
profiles_source = pd.read_csv(os.path.join(pp_dir, 'test', 'api', 'release_cycle',
'release_control_test_source_profiles.csv'), index_col=0)
ds_sink = DFData(profiles_sink)
ds_source = DFData(profiles_source)
const_sink = control.ConstControl(net, element='sink', variable='mdot_kg_per_s',
element_index=net.sink.index.values, data_source=ds_sink,
profile_name=net.sink.index.values.astype(str))
const_source = control.ConstControl(net, element='source', variable='mdot_kg_per_s',
element_index=net.source.index.values,
data_source=ds_source,
profile_name=net.source.index.values.astype(str))
const_sink.initial_run = False
const_source.initial_run = False
pp.pipeflow(net)
pp.to_json(net, os.path.join(path, 'example_%s.json' % pp.__version__))
return net
def test_create_network():
net = pandapipes.create_empty_network(fluid=3)
with pytest.raises(AttributeError):
pandapipes.get_fluid(net)
def test_temperature_internal_nodes_single_pipe():
"""
:return:
:rtype:
"""
net = pandapipes.create_empty_network("net", add_stdtypes=False)
d = 75e-3
pandapipes.create_junction(net, pn_bar=5, tfluid_k=283)
pandapipes.create_junction(net, pn_bar=5, tfluid_k=283)
pandapipes.create_pipe_from_parameters(net,
0,
1,
6,
d,
k_mm=.1,
sections=6,
alpha_w_per_m2k=5)
pandapipes.create_ext_grid(net, 0, p_bar=5, t_k=330, type="pt")
pandapipes.create_sink(net, 1, mdot_kg_per_s=1)
pandapipes.create_fluid_from_lib(net, "water", overwrite=True)
pandapipes.pipeflow(net,
stop_condition="tol",
iter=3,
friction_model="nikuradse",
mode="all",
transient=False,
nonlinear_method="automatic",
tol_p=1e-4,
tol_v=1e-4)
pipe_results = Pipe.get_internal_results(net, [0])
data = pd.read_csv(os.path.join(internals_data_path,
"Temperature_one_pipe_an.csv"),
sep=';',
header=0,
keep_default_na=False)
temp_an = data["T"]
pipe_temp_data_idx = np.where(pipe_results["TINIT"][:, 0] == 0)
pipe_temp_data = pipe_results["TINIT"][pipe_temp_data_idx, 1]
node_pit = net["_pit"]["node"]
junction_idx_lookup = get_lookup(net, "node",
"index")[Junction.table_name()]
from_junction_nodes = junction_idx_lookup[net["pipe"]
["from_junction"].values]
to_junction_nodes = junction_idx_lookup[net["pipe"]["to_junction"].values]
temp_pandapipes = np.zeros(len(pipe_temp_data[0]) + 2)
temp_pandapipes[0] = node_pit[from_junction_nodes[0], TINIT]
temp_pandapipes[1:-1] = pipe_temp_data[:]
temp_pandapipes[-1] = node_pit[to_junction_nodes[0], TINIT]
temp_diff = np.abs(1 - temp_pandapipes / temp_an)
assert np.all(temp_diff < 0.01)
import pandapipes as pp
net = pp.create_empty_network(fluid='lgas')
j1 = pp.create_junction(net, pn_bar=1.05, tfluid_k=293.15, name="Junction 1")
j2 = pp.create_junction(net, pn_bar=1.05, tfluid_k=293.15, name="Junction 2")
j3 = pp.create_junction(net, pn_bar=1.05, tfluid_k=293.15, name="Junction 3")
ext_grid = pp.create_ext_grid(net,
junction=j1,
p_bar=1.1,
t_k=293.15,
name="Grid Connection")
sink = pp.create_sink(net, junction=j3, mdot_kg_per_s=0.045, name="Sink")
pipe = pp.create_pipe(net,
from_junction=j1,
to_junction=j2,
length_km=0.1,
diameter_m=0.05,
name="Pipe 1")
valve = pp.create_valve(net,
from_junction=j2,
to_junction=j3,
diameter_m=0.05,
opened=True,
name="Valve 1")
pp.pipeflow(net)
print(net.res_junction)
print(net.res_pip)
def load_net():
# create test network
net = pandapipes.create_empty_network("test_net", fluid="lgas")
j1 = pandapipes.create_junction(net,
pn_bar=1.05,
tfluid_k=293.15,
name="Connection to External Grid",
geodata=(0, 0))
j2 = pandapipes.create_junction(net,
pn_bar=1.05,
tfluid_k=293.15,
name="Junction 2",
geodata=(2, 0))
j3 = pandapipes.create_junction(net,
pn_bar=1.05,
tfluid_k=293.15,
name="Junction 3",
geodata=(7, 4))
j4 = pandapipes.create_junction(net,
pn_bar=1.05,
tfluid_k=293.15,
name="Junction 4",
geodata=(7, -4))
j5 = pandapipes.create_junction(net,
pn_bar=1.05,
tfluid_k=293.15,
name="Junction 5",
geodata=(5, 3))
j6 = pandapipes.create_junction(net,
pn_bar=1.05,
tfluid_k=293.15,
name="Junction 6",
geodata=(5, -3))
pandapipes.create_ext_grid(net,
junction=j1,
p_bar=1.1,
t_k=293.15,
name="Grid Connection")
pandapipes.create_pipe_from_parameters(net,
from_junction=j1,
to_junction=j2,
length_km=10,
diameter_m=0.05,
name="Pipe 1",
geodata=[(0, 0), (2, 0)])
pandapipes.create_pipe_from_parameters(net,
from_junction=j2,
to_junction=j3,
length_km=2,
diameter_m=0.05,
name="Pipe 2",
geodata=[(2, 0), (2, 4), (7, 4)])
pandapipes.create_pipe_from_parameters(net,
from_junction=j2,
to_junction=j4,
length_km=2.5,
diameter_m=0.05,
name="Pipe 3",
geodata=[(2, 0), (2, -4), (7, -4)])
pandapipes.create_pipe_from_parameters(net,
from_junction=j3,
to_junction=j5,
length_km=1,
diameter_m=0.05,
name="Pipe 4",
geodata=[(7, 4), (7, 3), (5, 3)])
pandapipes.create_pipe_from_parameters(net,
from_junction=j4,
to_junction=j6,
length_km=1,
diameter_m=0.05,
name="Pipe 5",
geodata=[(7, -4), (7, -3), (5, -3)])
pandapipes.create_valve(net,
from_junction=j5,
to_junction=j6,
diameter_m=0.05,
opened=True)
pandapipes.create_sink(net,
junction=j4,
mdot_kg_per_s=5.45e-5,
name="Sink 1")
pandapipes.create_source(net, junction=j3, mdot_kg_per_s=3.45e-5)
return net
def create_base_net(oos, additional_pumps=True):
net = pandapipes.create_empty_network(fluid="lgas")
# create network elements, such as junctions, external grid, pipes, valves, sinks and sources
junction1 = pandapipes.create_junction(net, pn_bar=1.05, tfluid_k=293.15, in_service=not oos,
name="Connection to External Grid", geodata=(0, 0))
junction2 = pandapipes.create_junction(net, pn_bar=1.05, tfluid_k=293.15, name="Junction 2",
geodata=(2, 0), in_service=not oos)
junction3 = pandapipes.create_junction(net, pn_bar=1.05, tfluid_k=293.15, name="Junction 3",
geodata=(7, 4), in_service=not oos)
junction4 = pandapipes.create_junction(net, pn_bar=1.05, tfluid_k=293.15, name="Junction 4",
geodata=(7, -4), in_service=not oos)
junction5 = pandapipes.create_junction(net, pn_bar=1.05, tfluid_k=293.15, name="Junction 5",
geodata=(5, 3), in_service=not oos)
junction6 = pandapipes.create_junction(net, pn_bar=1.05, tfluid_k=293.15, name="Junction 6",
geodata=(5, -3))
junction7 = pandapipes.create_junction(net, pn_bar=1.05, tfluid_k=293.15, name="Junction 7",
geodata=(9, -4))
junction8 = pandapipes.create_junction(net, pn_bar=1.05, tfluid_k=293.15, name="Junction 8",
geodata=(9, 4))
junction9 = pandapipes.create_junction(net, pn_bar=1.05, tfluid_k=293.15, name="Junction 9",
geodata=(9, 0))
pandapipes.create_ext_grid(net, junction=junction1, p_bar=1.1, t_k=293.15,
name="Grid Connection")
pandapipes.create_pipe_from_parameters(net, from_junction=junction1, to_junction=junction2,
length_km=10, diameter_m=0.3, name="Pipe 1",
geodata=[(0, 0), (2, 0)], in_service=not oos)
pandapipes.create_pipe_from_parameters(net, from_junction=junction2, to_junction=junction3,
length_km=2, diameter_m=0.3, name="Pipe 2",
geodata=[(2, 0), (2, 4), (7, 4)], in_service=not oos)
pandapipes.create_pipe_from_parameters(net, from_junction=junction2, to_junction=junction4,
length_km=2.5, diameter_m=0.3, name="Pipe 3",
geodata=[(2, 0), (2, -4), (7, -4)], in_service=not oos)
pandapipes.create_pipe_from_parameters(net, from_junction=junction3, to_junction=junction5,
length_km=1, diameter_m=0.3, name="Pipe 4",
geodata=[(7, 4), (7, 3), (5, 3)])
pandapipes.create_pipe_from_parameters(net, from_junction=junction4, to_junction=junction6,
length_km=1, diameter_m=0.3, name="Pipe 5",
geodata=[(7, -4), (7, -3), (5, -3)])
pandapipes.create_pipe_from_parameters(net, from_junction=junction7, to_junction=junction8,
length_km=1, diameter_m=0.3, name="Pipe 6",
geodata=[(9, -4), (9, 0)])
pandapipes.create_pipe_from_parameters(net, from_junction=junction7, to_junction=junction8,
length_km=1, diameter_m=0.3, name="Pipe 7",
geodata=[(9, 0), (9, 4)])
pandapipes.create_valve(net, from_junction=junction5, to_junction=junction6, diameter_m=0.05,
opened=True)
pandapipes.create_heat_exchanger(net, junction3, junction8, diameter_m=0.3, qext_w=20000)
pandapipes.create_sink(net, junction=junction4, mdot_kg_per_s=0.545, name="Sink 1")
pandapipes.create_source(net, junction=junction3, mdot_kg_per_s=0.234)
pandapipes.create_pump_from_parameters(net, junction4, junction7, 'P1')
if additional_pumps:
pandapipes.create_circ_pump_const_mass_flow(net, junction9, junction5, 1.05, 1)
pandapipes.create_circ_pump_const_pressure(net, junction9, junction6, 1.05, 0.5)
if oos:
pandapipes.create_ext_grid(net, junction=junction1, p_bar=1.1, t_k=293.15,
name="Grid Connection", in_service=False)
pandapipes.create_heat_exchanger(net, junction3, junction8, diameter_m=0.3, qext_w=20000,
in_service=False)
pandapipes.create_sink(net, junction=junction4, mdot_kg_per_s=0.545, name="Sink 2",
in_service=False)
pandapipes.create_pump_from_parameters(net, junction4, junction7, 'P1', in_service=False)
pandapipes.create_source(net, junction=junction3, mdot_kg_per_s=0.234, in_service=False)
pandapipes.create_circ_pump_const_mass_flow(net, junction9, junction5, 1.05, 1,
in_service=False)
pandapipes.create_circ_pump_const_pressure(net, junction9, junction6, 1.05, 0.5,
in_service=False)
return net
def test_create_network():
net = pandapipes.create_empty_network(fluid=3)
with pytest.raises(UserWarning):
pandapipes.get_fluid(net)
import pandapipes as pp
# create empty net
net = pp.create_empty_network(fluid="lgas")
# create junction
j1 = pp.create_junction(net, pn_bar=1.05, tfluid_k=293.15, name="Junction 1")
j2 = pp.create_junction(net, pn_bar=1.05, tfluid_k=293.15, name="Junction 2")
j3 = pp.create_junction(net, pn_bar=1.05, tfluid_k=293.15, name="Junction 3")
# create junction elements
ext_grid = pp.create_ext_grid(net,
junction=j1,
p_bar=1.1,
t_k=293.15,
name="Grid Connection")
sink = pp.create_sink(net, junction=j3, mdot_kg_per_s=0.045, name="Sink")
# create branch elements
pipe = pp.create_pipe_from_parameters(net,
from_junction=j1,
to_junction=j2,
length_km=0.1,
diameter_m=0.05,
name="Pipe")
valve = pp.create_valve(net,
from_junction=j2,
to_junction=j3,
diameter_m=0.05,
opened=True,
name="Valve")
def test_pipe_velocity_results():
"""
This test verifies the entries in the results table for a pipe network with pipes consisting of
more than one section. The basic idea is that the computation is first done with only one section per pipe.
Afterwards, the same network is calculated with more nodes. if everything works correctly, the entries
in the result table for the velocities (from, to) should be the same.
A T-junction is used for the test setup
:return:
:rtype:
"""
net = pandapipes.create_empty_network("net", add_stdtypes=False)
d = 209.1e-3
pandapipes.create_junction(net, pn_bar=51, tfluid_k=285.15)
pandapipes.create_junction(net, pn_bar=51, tfluid_k=285.15)
pandapipes.create_junction(net, pn_bar=51, tfluid_k=285.15)
pandapipes.create_junction(net, pn_bar=51, tfluid_k=285.15)
pandapipes.create_pipe_from_parameters(net,
0,
1,
6.0,
d,
k_mm=.5,
sections=1)
pandapipes.create_pipe_from_parameters(net,
1,
2,
6.0,
d,
k_mm=.5,
sections=1)
pandapipes.create_pipe_from_parameters(net,
1,
3,
6.0,
d,
k_mm=.5,
sections=1)
pandapipes.create_ext_grid(net,
0,
p_bar=51 - 1.01325,
t_k=285.15,
type="pt")
pandapipes.create_sink(net, 2, mdot_kg_per_s=0.82752 * 45000 / 3600 / 3)
pandapipes.create_sink(net, 3, mdot_kg_per_s=0.82752 * 45000 / 3600 / 2)
_add_fluid_to_net(
net,
pandapipes.create_constant_fluid(name="natural_gas",
fluid_type="gas",
viscosity=11.93e-6,
heat_capacity=2185,
compressibility=1,
der_compressibility=0,
density=0.82752))
pandapipes.pipeflow(net,
stop_condition="tol",
iter=70,
friction_model="nikuradse",
transient=False,
nonlinear_method="automatic",
tol_p=1e-4,
tol_v=1e-4)
v_1_sec_from = net.res_pipe.v_from_m_per_s
v_1_sec_to = net.res_pipe.v_from_m_per_s
net = pandapipes.create_empty_network("net", add_stdtypes=False)
d = 209.1e-3
pandapipes.create_junction(net, pn_bar=51, tfluid_k=285.15)
pandapipes.create_junction(net, pn_bar=51, tfluid_k=285.15)
pandapipes.create_junction(net, pn_bar=51, tfluid_k=285.15)
pandapipes.create_junction(net, pn_bar=51, tfluid_k=285.15)
pandapipes.create_pipe_from_parameters(net,
0,
1,
6.0,
d,
k_mm=.5,
sections=3)
pandapipes.create_pipe_from_parameters(net,
1,
2,
6.0,
d,
k_mm=.5,
sections=4)
pandapipes.create_pipe_from_parameters(net,
1,
3,
6.0,
d,
k_mm=.5,
sections=2)
pandapipes.create_ext_grid(net,
0,
p_bar=51 - 1.01325,
t_k=285.15,
type="pt")
pandapipes.create_sink(net, 2, mdot_kg_per_s=0.82752 * 45000 / 3600 / 3)
pandapipes.create_sink(net, 3, mdot_kg_per_s=0.82752 * 45000 / 3600 / 2)
_add_fluid_to_net(
net,
pandapipes.create_constant_fluid(name="natural_gas",
fluid_type="gas",
viscosity=11.93e-6,
heat_capacity=2185,
compressibility=1,
der_compressibility=0,
density=0.82752))
pandapipes.pipeflow(net,
stop_condition="tol",
iter=70,
friction_model="nikuradse",
transient=False,
nonlinear_method="automatic",
tol_p=1e-4,
tol_v=1e-4)
v_n_sec_from = net.res_pipe.v_from_m_per_s
v_n_sec_to = net.res_pipe.v_from_m_per_s
diff_from = v_1_sec_from - v_n_sec_from
diff_to = v_1_sec_to - v_n_sec_to
assert np.all(np.abs(diff_from) < 1e-9)
assert np.all(np.abs(diff_to) < 1e-9)
def create_test_net():
"""
:return:
:rtype:
"""
net = pandapipes.create_empty_network()
j1 = pandapipes.create_junction(net, 1, 293.15)
j2 = pandapipes.create_junction(net, 1, 293.15)
j3 = pandapipes.create_junction(net, 1, 293.15, in_service=False)
j4 = pandapipes.create_junction(net, 1, 293.15)
j5 = pandapipes.create_junction(net, 1, 293.15)
j6 = pandapipes.create_junction(net, 1, 293.15)
j7 = pandapipes.create_junction(net, 1, 293.15, in_service=False)
pandapipes.create_ext_grid(net, j1, 1, 285.15, type="pt")
pandapipes.create_pipe_from_parameters(net,
j1,
j2,
0.1,
0.1,
sections=1,
alpha_w_per_m2k=5)
pandapipes.create_pipe_from_parameters(net,
j2,
j3,
0.1,
0.1,
sections=2,
alpha_w_per_m2k=5,
in_service=False)
pandapipes.create_pipe_from_parameters(net,
j4,
j6,
0.1,
0.1,
sections=2,
alpha_w_per_m2k=5)
pandapipes.create_pipe_from_parameters(net,
j6,
j7,
0.1,
0.1,
sections=1,
alpha_w_per_m2k=5,
in_service=False)
pandapipes.create_pipe_from_parameters(net,
j1,
j5,
0.1,
0.1,
sections=2,
alpha_w_per_m2k=5)
pandapipes.create_valve(net, j1, j4, 0.1)
pandapipes.create_valve(net, j4, j5, 0.1, opened=False)
pandapipes.create_sink(net, j2, mdot_kg_per_s=0.1)
pandapipes.create_sink(net, j3, mdot_kg_per_s=0.1)
pandapipes.create_sink(net, j5, mdot_kg_per_s=0.1)
pandapipes.create_sink(net, j6, mdot_kg_per_s=0.1)
pandapipes.create_sink(net, j7, mdot_kg_per_s=0.1)
return net
def net_plotting():
net = pandapipes.create_empty_network(fluid="lgas")
# create network elements, such as junctions, external grid, pipes, valves, sinks and sources
junction1 = pandapipes.create_junction(net,
pn_bar=1.05,
tfluid_k=293.15,
name="Connection to External Grid",
geodata=(0, 0))
junction2 = pandapipes.create_junction(net,
pn_bar=1.05,
tfluid_k=293.15,
name="Junction 2",
geodata=(2, 0))
junction3 = pandapipes.create_junction(net,
pn_bar=1.05,
tfluid_k=293.15,
name="Junction 3",
geodata=(7, 4))
junction4 = pandapipes.create_junction(net,
pn_bar=1.05,
tfluid_k=293.15,
name="Junction 4",
geodata=(7, -4))
junction5 = pandapipes.create_junction(net,
pn_bar=1.05,
tfluid_k=293.15,
name="Junction 5",
geodata=(5, 3))
junction6 = pandapipes.create_junction(net,
pn_bar=1.05,
tfluid_k=293.15,
name="Junction 6",
geodata=(5, -3))
junction7 = pandapipes.create_junction(net,
pn_bar=1.05,
tfluid_k=293.15,
name="Junction 7",
geodata=(9, -4))
junction8 = pandapipes.create_junction(net,
pn_bar=1.05,
tfluid_k=293.15,
name="Junction 7",
geodata=(9, 4))
pandapipes.create_ext_grid(net,
junction=junction1,
p_bar=1.1,
t_k=293.15,
name="Grid Connection")
pandapipes.create_pipe_from_parameters(net,
from_junction=junction1,
to_junction=junction2,
length_km=10,
diameter_m=0.3,
name="Pipe 1",
geodata=[(0, 0), (2, 0)])
pandapipes.create_pipe_from_parameters(net,
from_junction=junction2,
to_junction=junction3,
length_km=2,
diameter_m=0.3,
name="Pipe 2",
geodata=[(2, 0), (2, 4), (7, 4)])
pandapipes.create_pipe_from_parameters(net,
from_junction=junction2,
to_junction=junction4,
length_km=2.5,
diameter_m=0.3,
name="Pipe 3",
geodata=[(2, 0), (2, -4), (7, -4)])
pandapipes.create_pipe_from_parameters(net,
from_junction=junction3,
to_junction=junction5,
length_km=1,
diameter_m=0.3,
name="Pipe 4",
geodata=[(7, 4), (7, 3), (5, 3)])
pandapipes.create_pipe_from_parameters(net,
from_junction=junction4,
to_junction=junction6,
length_km=1,
diameter_m=0.3,
name="Pipe 5",
geodata=[(7, -4), (7, -3), (5, -3)])
pandapipes.create_pipe_from_parameters(net,
from_junction=junction7,
to_junction=junction8,
length_km=1,
diameter_m=0.3,
name="Pipe 6",
geodata=[(9, -4), (9, 4)])
pandapipes.create_valve(net,
from_junction=junction5,
to_junction=junction6,
diameter_m=0.05,
opened=True)
pandapipes.create_heat_exchanger(net,
junction3,
junction8,
diameter_m=0.3,
qext_w=20000)
pandapipes.create_sink(net,
junction=junction4,
mdot_kg_per_s=0.545,
name="Sink 1")
pandapipes.create_source(net, junction=junction3, mdot_kg_per_s=0.234)
pandapipes.create_pump_from_parameters(net, junction4, junction7, 'P1')
return net
def complex_heat_connectivity_grid():
net = pandapipes.create_empty_network()
pandapipes.create_fluid_from_lib(net, "water")
j1 = pandapipes.create_junction(net, 1, 320.15, index=1)
j2 = pandapipes.create_junction(net, 1, 320.15, index=2)
j3 = pandapipes.create_junction(net, 1, 320.15, index=3)
j4 = pandapipes.create_junction(net, 1, 320.15, index=4, in_service=False)
j5 = pandapipes.create_junction(net, 1, 320.15, index=5)
j6 = pandapipes.create_junction(net, 1, 320.15, index=6)
j7 = pandapipes.create_junction(net, 1, 320.15, index=7)
j8 = pandapipes.create_junction(net, 1, 320.15, index=8)
j9 = pandapipes.create_junction(net, 1, 320.15, index=9)
j10 = pandapipes.create_junction(net, 1, 320.15, index=10)
pandapipes.create_ext_grid(net, j1, 1, 320.15, type="p", index=5)
pandapipes.create_ext_grid(net, j7, 1, 320.15, type="t", index=2)
pandapipes.create_ext_grid(net, j10, 1, 320.15, type="pt", index=1)
pandapipes.create_pipe_from_parameters(net,
j1,
j2,
0.1,
0.1,
alpha_w_per_m2k=5,
index=3)
pandapipes.create_pipe_from_parameters(net,
j1,
j3,
0.1,
0.1,
alpha_w_per_m2k=5,
index=4)
pandapipes.create_pipe_from_parameters(net,
j2,
j4,
0.1,
0.1,
alpha_w_per_m2k=5,
in_service=False,
index=5)
pandapipes.create_pipe_from_parameters(net,
j3,
j5,
0.1,
0.1,
alpha_w_per_m2k=5,
in_service=False,
index=7)
pandapipes.create_pipe_from_parameters(net,
j6,
j7,
0.1,
0.1,
alpha_w_per_m2k=5,
index=9)
pandapipes.create_pipe_from_parameters(net,
j5,
8,
0.1,
0.1,
alpha_w_per_m2k=5,
in_service=False,
index=8)
pandapipes.create_pipe_from_parameters(net,
j8,
j10,
0.1,
0.1,
alpha_w_per_m2k=5,
index=1)
pandapipes.create_pipe_from_parameters(net,
j9,
j10,
0.1,
0.1,
alpha_w_per_m2k=5,
index=2)
pandapipes.create_valve(net, j5, j6, 0.1, index=10)
pandapipes.create_valve(net, j4, j5, 0.1, opened=False, index=12)
pandapipes.create_sink(net, j3, mdot_kg_per_s=0.1, index=3)
pandapipes.create_sink(net, j4, mdot_kg_per_s=0.1, index=4)
pandapipes.create_sink(net, j7, mdot_kg_per_s=0.2, index=5)
pandapipes.create_sink(net, j9, mdot_kg_per_s=0.1, index=7)
pandapipes.create_sink(net, j8, mdot_kg_per_s=0.1, index=1)
pandapipes.create_source(net, j5, mdot_kg_per_s=0.1, index=7)
pandapipes.create_source(net, j2, mdot_kg_per_s=0.05, index=2)
return net
def create_empty_net():
return pandapipes.create_empty_network()
def test_valve():
"""
:return:
:rtype:
"""
net = pandapipes.create_empty_network("net", add_stdtypes=True)
j0 = pandapipes.create_junction(net, pn_bar=5, tfluid_k=283.15, index=5)
j1 = pandapipes.create_junction(net, pn_bar=5, tfluid_k=283.15, index=3)
j2 = pandapipes.create_junction(net, pn_bar=5, tfluid_k=283.15, index=6)
j3 = pandapipes.create_junction(net, pn_bar=5, tfluid_k=283.15, index=9)
j4 = pandapipes.create_junction(net, pn_bar=5, tfluid_k=283.15, index=20)
j5 = pandapipes.create_junction(net, pn_bar=5, tfluid_k=283.15, index=45)
j6 = pandapipes.create_junction(net, pn_bar=5, tfluid_k=283.15, index=4)
j7 = pandapipes.create_junction(net, pn_bar=5, tfluid_k=283.15, index=8)
pandapipes.create_ext_grid(net, j0, 5, 283.15, type="p")
pandapipes.create_pipe_from_parameters(net,
j0,
j1,
diameter_m=.1,
k=0.1,
length_km=1.)
pandapipes.create_pipe_from_parameters(net,
j3,
j4,
diameter_m=.1,
k=0.1,
length_km=.5)
pandapipes.create_pipe_from_parameters(net,
j2,
j4,
diameter_m=.1,
k=0.1,
length_km=.5)
pandapipes.create_pipe_from_parameters(net,
j5,
j4,
diameter_m=.1,
k=0.1,
length_km=.35)
pandapipes.create_pipe_from_parameters(net,
j1,
j6,
diameter_m=.1,
k=0.1,
length_km=.1,
loss_coefficient=9000)
pandapipes.create_pipe_from_parameters(net,
j1,
j7,
diameter_m=.1,
k=0.1,
length_km=.1,
loss_coefficient=9000)
pandapipes.create_valve(net, j6, j2, diameter_m=0.1, opened=False)
pandapipes.create_valve(net, j7, j3, diameter_m=0.1, opened=True)
pandapipes.create_sink(net, j5, 0.11667)
pandapipes.create_fluid_from_lib(net, "lgas", overwrite=True)
pandapipes.pipeflow(net,
stop_condition="tol",
iter=10,
friction_model="nikuradse",
mode="hydraulics",
transient=False,
nonlinear_method="automatic",
tol_p=1e-4,
tol_v=1e-4)
data = pd.read_csv(os.path.join(internals_data_path, "test_valve.csv"),
sep=';')
data_p = data['p'].dropna(inplace=False)
data_v = data['v'].dropna(inplace=False)
res_junction = net.res_junction.p_bar.values
res_pipe = net.res_pipe.v_mean_m_per_s.values
zeros = res_pipe == 0
test_zeros = data_v.values == 0
check_zeros = zeros == test_zeros
assert np.all(check_zeros)
p_diff = np.abs(1 - res_junction / data_p[data_p != 0].values)
v_diff = np.abs(1 - res_pipe[res_pipe != 0] / data_v[data_v != 0].values)
assert np.all(p_diff < 0.01)
assert np.all(v_diff < 0.01)
def test_gas_internal_nodes():
"""
:return:
:rtype:
"""
net = pandapipes.create_empty_network("net", add_stdtypes=False)
d = 209.1e-3
pandapipes.create_junction(net, pn_bar=51, tfluid_k=285.15)
pandapipes.create_junction(net, pn_bar=51, tfluid_k=285.15)
pandapipes.create_pipe_from_parameters(net,
0,
1,
12.0,
d,
k_mm=.5,
sections=12)
pandapipes.create_ext_grid(net,
0,
p_bar=51 - 1.01325,
t_k=285.15,
type="pt")
pandapipes.create_sink(net, 1, mdot_kg_per_s=0.82752 * 45000 / 3600)
_add_fluid_to_net(
net,
pandapipes.create_constant_fluid(name="natural_gas",
fluid_type="gas",
viscosity=11.93e-6,
heat_capacity=2185,
compressibility=1,
der_compressibility=0,
density=0.82752))
pandapipes.pipeflow(net,
stop_condition="tol",
iter=70,
friction_model="nikuradse",
transient=False,
nonlinear_method="automatic",
tol_p=1e-4,
tol_v=1e-4)
pipe_results = Pipe.get_internal_results(net, [0])
data = pd.read_csv(os.path.join(internals_data_path,
"gas_sections_an.csv"),
sep=';',
header=0,
keep_default_na=False)
p_an = data["p1"] / 1e5
v_an = data["v"]
v_an = v_an.drop([0])
pipe_p_data_idx = np.where(pipe_results["PINIT"][:, 0] == 0)
pipe_v_data_idx = np.where(pipe_results["VINIT_MEAN"][:, 0] == 0)
pipe_p_data = pipe_results["PINIT"][pipe_p_data_idx, 1]
pipe_v_data = pipe_results["VINIT_MEAN"][pipe_v_data_idx, 1]
node_pit = net["_pit"]["node"]
junction_idx_lookup = get_lookup(net, "node",
"index")[Junction.table_name()]
from_junction_nodes = junction_idx_lookup[net["pipe"]
["from_junction"].values]
to_junction_nodes = junction_idx_lookup[net["pipe"]["to_junction"].values]
p_pandapipes = np.zeros(len(pipe_p_data[0]) + 2)
p_pandapipes[0] = node_pit[from_junction_nodes[0], PINIT]
p_pandapipes[1:-1] = pipe_p_data[:]
p_pandapipes[-1] = node_pit[to_junction_nodes[0], PINIT]
p_pandapipes = p_pandapipes + 1.01325
v_pandapipes = pipe_v_data[0, :]
p_diff = np.abs(1 - p_pandapipes / p_an)
v_diff = np.abs(v_an - v_pandapipes)
assert np.all(p_diff < 0.01)
assert np.all(v_diff < 0.4)
def test_collection_lengths():
net = pandapipes.create_empty_network(add_stdtypes=False)
d = 40e-3
j1 = pandapipes.create_junction(net,
pn_bar=5,
tfluid_k=293.15,
geodata=(0, 0))
j2 = pandapipes.create_junction(net,
pn_bar=5,
tfluid_k=293.15,
geodata=(2, 0))
j3 = pandapipes.create_junction(net,
pn_bar=5,
tfluid_k=293.15,
geodata=(4, 2))
j4 = pandapipes.create_junction(net,
pn_bar=5,
tfluid_k=293.15,
geodata=(4, 0))
j5 = pandapipes.create_junction(net,
pn_bar=5,
tfluid_k=293.15,
geodata=(4, -2))
j6 = pandapipes.create_junction(net,
pn_bar=5,
tfluid_k=293.15,
geodata=(6, 0))
j7 = pandapipes.create_junction(net,
pn_bar=5,
tfluid_k=293.15,
geodata=(8, 0))
j8 = pandapipes.create_junction(net,
pn_bar=5,
tfluid_k=293.15,
geodata=(10, 0))
pandapipes.create_ext_grid(net, j1, p_bar=5, t_k=293.15)
pandapipes.create_sink(net, j5, mdot_kg_per_s=0.5)
pandapipes.create_sink(net, j6, mdot_kg_per_s=0.5)
pandapipes.create_sink(net, j8, mdot_kg_per_s=0.5)
pandapipes.create_pipe_from_parameters(net,
j1,
j2,
0.1,
diameter_m=d,
k_mm=0.1,
geodata=[(0, 0), (2, 0)])
pandapipes.create_valve(net,
j2,
j3,
diameter_m=d,
opened=True,
loss_coefficient=5e3)
pandapipes.create_valve(net,
j2,
j4,
diameter_m=d,
opened=False,
loss_coefficient=0.01)
pandapipes.create_valve(net,
j2,
j5,
diameter_m=d,
opened=False,
loss_coefficient=0.01)
pandapipes.create_pipe_from_parameters(net,
j3,
j6,
0.1,
diameter_m=d,
k_mm=0.1,
geodata=[(4, 2), (6, 2), (6, 0)])
pandapipes.create_pipe_from_parameters(net,
j4,
j6,
0.2,
diameter_m=d,
k_mm=0.1,
geodata=[(4, 0), (6, 0)])
pandapipes.create_pipe_from_parameters(net,
j5,
j6,
0.2,
diameter_m=d,
k_mm=0.1,
geodata=[(4, -2), (6, -2), (6, 0)])
pandapipes.create_heat_exchanger(net, j6, j7, d, qext_w=20000)
pandapipes.create_pump_from_parameters(net, j7, j8, 'P1')
pipe_coll_direct = plot.create_pipe_collection(net,
use_junction_geodata=True)
pipe_coll_real = plot.create_pipe_collection(net)
assert np.all(
pipe_coll_direct.indices.astype(np.int64) ==
net.pipe.index.values.astype(np.int64))
assert len(pipe_coll_direct.get_paths()) == len(net.pipe)
assert all([len(p) == 2 for p in pipe_coll_direct.get_paths()])
assert all([
len(p) == len(net.pipe_geodata.coords.iloc[i])
for i, p in enumerate(pipe_coll_real.get_paths())
])
junction_coll = plot.create_junction_collection(net)
assert np.all(
junction_coll.node_indices.astype(np.int64) ==
net.junction.index.values.astype(np.int64))
assert len(junction_coll.get_paths()) == len(net.junction)
valve_coll_patches, valve_coll_lines = plot.create_valve_collection(net)
assert len(valve_coll_patches.get_paths()) == 2 * len(net.valve)
assert len(valve_coll_lines.get_paths()) == 2 * len(net.valve)
pump_coll_patches, pump_coll_lines = plot.create_pump_collection(net)
assert len(pump_coll_patches.get_paths()) == len(net.pump)
assert len(pump_coll_lines.get_paths()) == 4 * len(net.pump)
hex_coll_patches, hex_coll_lines = plot.create_heat_exchanger_collection(
net)
assert len(hex_coll_patches.get_paths()) == 2 * len(net.heat_exchanger)
assert len(hex_coll_lines.get_paths()) == 2 * len(net.heat_exchanger)
def test_circulation_pump_constant_mass():
"""
:return:
:rtype:
"""
net = pandapipes.create_empty_network("net", add_stdtypes=False)
j1 = pandapipes.create_junction(net, pn_bar=5, tfluid_k=283.15)
j2 = pandapipes.create_junction(net, pn_bar=5, tfluid_k=283.15)
j3 = pandapipes.create_junction(net, pn_bar=5, tfluid_k=283.15)
j4 = pandapipes.create_junction(net, pn_bar=5, tfluid_k=283.15)
pandapipes.create_pipe_from_parameters(net,
j1,
j2,
k_mm=1.,
length_km=0.43380,
diameter_m=0.1022)
pandapipes.create_pipe_from_parameters(net,
j3,
j4,
k_mm=1.,
length_km=0.26370,
diameter_m=0.1022)
pandapipes.create_circ_pump_const_mass_flow(net,
j1,
j4,
5,
5,
300,
type='pt')
pandapipes.create_heat_exchanger(net, j2, j3, 0.1, qext_w=200000)
pandapipes.create_sink(net, j1, 2)
pandapipes.create_source(net, j4, 2)
pandapipes.create_fluid_from_lib(net, "water", overwrite=True)
pandapipes.pipeflow(net,
stop_condition="tol",
iter=10,
friction_model="nikuradse",
mode="all",
transient=False,
nonlinear_method="automatic",
tol_p=1e-4,
tol_v=1e-4)
data = pd.read_csv(os.path.join(internals_data_path,
"test_circ_pump_mass.csv"),
sep=';')
res_junction = net.res_junction
res_pipe = net.res_pipe.v_mean_m_per_s.values
res_pump = net.res_circ_pump_mass
p_diff = np.abs(1 - res_junction.p_bar.values / data['p'].dropna().values)
t_diff = np.abs(1 - res_junction.t_k.values / data['t'].dropna().values)
v_diff = np.abs(1 - res_pipe / data['v'].dropna().values)
mdot_diff = np.abs(1 - res_pump['mdot_kg_per_s'].values /
data['mdot'].dropna().values)
deltap_diff = np.abs(1 - res_pump['deltap_bar'].values /
data['deltap'].dropna().values)
assert np.all(p_diff < 0.01)
assert np.all(t_diff < 0.01)
assert np.all(v_diff < 0.01)
assert np.all(mdot_diff < 0.01)
assert np.all(deltap_diff < 0.01)