def test_random_net_and_one_node_net(create_test_net): """ :param create_test_net: :type create_test_net: :return: :rtype: """ net = copy.deepcopy(create_test_net) pp.create_fluid_from_lib(net, "water") j = create_junction(net, 1, 298.15) create_ext_grid(net, j, 1, 298.15) create_sink(net, j, 0.01) create_source(net, j, 0.02) pp.pipeflow(net) net = copy.deepcopy(create_test_net) pp.create_fluid_from_lib(net, "lgas") j = create_junction(net, 1, 298.15) create_ext_grid(net, j, 1, 298.15) create_sink(net, j, 0.01) create_source(net, j, 0.02) pp.pipeflow(net) assert np.isclose( net.res_ext_grid.values[-1] + net.res_sink.values[-1] - net.res_source.values[-1], 0)
def test_temperature_internal_nodes_tee_2ab_1zu(): """ :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_ext_grid(net, j0, p_bar=5, t_k=350, type="pt") pandapipes.create_sink(net, j2, mdot_kg_per_s=1) pandapipes.create_sink(net, j3, mdot_kg_per_s=1) pandapipes.create_pipe_from_parameters(net, j0, j1, 2.5, d, k_mm=.1, alpha_w_per_m2k=5) pandapipes.create_pipe_from_parameters(net, j1, j2, 2.5, d, k_mm=.1, alpha_w_per_m2k=5) pandapipes.create_pipe_from_parameters(net, j1, j3, 2.5, d, k_mm=.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", 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_tee_2ab_1zu_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_pump_from_regression_parameteres(): """ :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_ext_grid(net, j1, 5, 283.15, type="p") pandapipes.create_pump_from_parameters(net, j2, j3, 'P1', regression_parameters=[ -1.48620799e-04, -1.29656785e-02, 6.10000000e+00 ]) pandapipes.create_sink(net, j4, 0.02333) pandapipes.create_fluid_from_lib(net, "lgas", overwrite=True) pandapipes.pipeflow(net, stop_condition="tol", iter=3, 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_pump.csv"), sep=';') res_junction = net.res_junction.p_bar.values res_pipe = net.res_pipe.v_mean_m_per_s.values p_diff = np.abs(1 - res_junction / data['p'].dropna().values) v_diff = np.abs(1 - res_pipe / data['v'].dropna().values) assert np.all(p_diff < 0.01) assert np.all(v_diff < 0.01)
def get_gas_example(): net_gas = g_nw.gas_meshed_square() pandapipes.create_fluid_from_lib(net_gas, "hgas", overwrite=True) net_gas.sink.drop(index=0, inplace=True) net_gas.junction.pn_bar = 30 net_gas.ext_grid.p_bar = 30 net_gas.pipe.diameter_m = 0.8 return net_gas
def test_set_user_pf_options(create_test_net): """ :param create_test_net: :type create_test_net: :return: :rtype: """ net = copy.deepcopy(create_test_net) pandapipes.create_fluid_from_lib(net, "lgas") necessary_options = {'mode': 'hydraulics'} pandapipes.pipeflow(net, **necessary_options) old_options = net._options.copy() test_options = {key: i for i, key in enumerate(old_options.keys())} pandapipes.pipeflow_setup.set_user_pf_options(net, hello='bye', **test_options) test_options.update({'hello': 'bye'}) test_options.update({'hyd_flag': True}) assert net.user_pf_options == test_options # remove what is in user_pf_options and add hello=world pandapipes.pipeflow_setup.set_user_pf_options(net, reset=True, hello='world') assert net.user_pf_options == {'hello': 'world'} # check if 'hello' is added to net._options, but other options are untouched pandapipes.pipeflow(net, **necessary_options) assert 'hello' in net._options.keys() and net._options['hello'] == 'world' net._options.pop('hello') assert net._options == old_options # check if user_pf_options can be deleted and net._options is as it was before pandapipes.pipeflow_setup.set_user_pf_options(net, reset=True, hello='world') pandapipes.pipeflow_setup.set_user_pf_options(net, reset=True) assert net.user_pf_options == {} pandapipes.pipeflow(net, **necessary_options) assert 'hello' not in net._options.keys() # see if user arguments overrule user_pf_options, but other user_pf_options still have the # priority pandapipes.pipeflow_setup.set_user_pf_options(net, reset=True, tol_p=1e-6, tol_v=1e-6) pandapipes.pipeflow(net, tol_p=1e-8, **necessary_options) assert net.user_pf_options['tol_p'] == 1e-6 assert net._options['tol_p'] == 1e-8 assert net._options['tol_v'] == 1e-6
def test_g2g_multiple(get_gas_example): """ coupling of multiple elements in two gas grids with one MulitEnergyController gas-to-gas = e.g. hgas (methane) to Hydrogen conversion (SMR)""" # get the nets fluid1 = {"name": "hgas", "cal_value": 14.62197} net_gas1 = copy.deepcopy(get_gas_example) pandapipes.create_fluid_from_lib(net_gas1, fluid1["name"], overwrite=True) fluid2 = {"name": "hydrogen", "cal_value": 38.38024} net_gas2 = copy.deepcopy(get_gas_example) pandapipes.create_fluid_from_lib(net_gas2, fluid2["name"], overwrite=True) # set up multinet mn = create_empty_multinet("test_g2g") add_nets_to_multinet(mn, hgas_net=net_gas1, hydrogen_net=net_gas2) # dummy component for offset in sink/source indices: _ = pandapipes.create_sink(net_gas1, 0, mdot_kg_per_s=0.001) no_g2g = pandapipes.create_sources(net_gas2, [0, 3], mdot_kg_per_s=0.0314) # add components to represent G2P unit gas1_cons_kg_per_s = 0.05 g2g_ids_cons = pandapipes.create_sinks(net_gas1, range(1, 4), mdot_kg_per_s=gas1_cons_kg_per_s, name="SMR consumption") g2g_ids_prod = pandapipes.create_sources(net_gas2, [0, 2, 5], 0, name="SMR production") # add coupling controller eta = 0.65 GasToGasConversion(mn, g2g_ids_cons, g2g_ids_prod, efficiency=eta, name_gas_net_from='hgas_net', name_gas_net_to='hydrogen_net') run_control(mn) # nets must not be changed assert mn.nets["hgas_net"] == net_gas1 assert mn.nets["hydrogen_net"] == net_gas2 # check G2G result assert np.all(net_gas1.sink.loc[g2g_ids_cons, "mdot_kg_per_s"] == \ net_gas1.res_sink.loc[g2g_ids_cons, "mdot_kg_per_s"]) assert np.all(net_gas1.sink.loc[g2g_ids_cons, "mdot_kg_per_s"] == gas1_cons_kg_per_s) assert np.all(net_gas2.source.loc[no_g2g, "mdot_kg_per_s"] == 0.0314) assert np.allclose( net_gas2.source.loc[g2g_ids_prod, "mdot_kg_per_s"], (gas1_cons_kg_per_s * fluid1["cal_value"] / fluid2["cal_value"]) * eta)
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_g2g_single(get_gas_example): """gas-to-gas = hgas (methane) to Hydrogen conversion""" # get the nets fluid1 = {"name": "hgas", "cal_value": 14.62197} net_gas1 = copy.deepcopy(get_gas_example) pandapipes.create_fluid_from_lib(net_gas1, fluid1["name"], overwrite=True) fluid2 = {"name": "hydrogen", "cal_value": 38.38024} net_gas2 = copy.deepcopy(get_gas_example) pandapipes.create_fluid_from_lib(net_gas2, fluid2["name"], overwrite=True) # set up multinet mn = create_empty_multinet("test_g2g") add_nets_to_multinet(mn, hgas_net=net_gas1, hydrogen_net=net_gas2) # add components to represent G2P unit gas1_cons_kg_per_s = 0.5 g2g_id_cons = pandapipes.create_sink(net_gas1, 1, mdot_kg_per_s=gas1_cons_kg_per_s, name="SMR consumption") g2g_id_prod = pandapipes.create_source(net_gas2, 1, 0, name="SMR production") # add coupling controller eta = 0.65 GasToGasConversion(mn, g2g_id_cons, g2g_id_prod, efficiency=eta, name_gas_net_from='hgas_net', name_gas_net_to='hydrogen_net') run_control(mn) # nets must not be changed assert mn.nets["hgas_net"] == net_gas1 assert mn.nets["hydrogen_net"] == net_gas2 # check G2G result assert net_gas1.sink.at[g2g_id_cons, "mdot_kg_per_s"] == \ net_gas1.res_sink.at[g2g_id_cons, "mdot_kg_per_s"] assert net_gas1.sink.at[g2g_id_cons, "mdot_kg_per_s"] == gas1_cons_kg_per_s assert np.isclose( net_gas2.source.at[g2g_id_prod, "mdot_kg_per_s"], (gas1_cons_kg_per_s * fluid1["cal_value"] / fluid2["cal_value"]) * eta)
def test_pump_bypass_high_vdot(): """ High flow: pressure lift not <0, always >=0 :return: :rtype: """ net = pandapipes.create_empty_network("net", add_stdtypes=True) 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(net, j1, j2, std_type='2000_ST<16', k_mm=0.1, length_km=0.1) pandapipes.create_pipe(net, j3, j4, std_type='2000_ST<16', k_mm=0.1, length_km=0.1) pandapipes.create_ext_grid(net, j1, 5, 283.15, type="p") pandapipes.create_pump(net, j2, j3, std_type='P1') pandapipes.create_sink(net, j4, 1000) pandapipes.create_fluid_from_lib(net, "hgas", overwrite=True) pandapipes.pipeflow(net, stop_condition="tol", iter=30, friction_model="nikuradse", mode="hydraulics", transient=False, nonlinear_method="automatic", tol_p=1e-4, tol_v=1e-4) assert net.res_pump.deltap_bar.isin([0]).all() assert np.isclose(net.res_junction.loc[1, "p_bar"], net.res_junction.loc[2, "p_bar"])
def test_pump_bypass_on_reverse_flow(): """ reverse flow = no pressure lift :return: :rtype: """ net = pandapipes.create_empty_network("net", add_stdtypes=True) 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(net, j1, j2, std_type='125_PE_80_SDR_11', k_mm=1., length_km=10) pandapipes.create_pipe(net, j3, j4, std_type='125_PE_80_SDR_11', k_mm=1., length_km=12) pandapipes.create_ext_grid(net, j1, 5, 283.15, type="p") pandapipes.create_pump(net, j2, j3, std_type='P1') pandapipes.create_source(net, j4, 0.02333) pandapipes.create_fluid_from_lib(net, "hgas", overwrite=True) pandapipes.pipeflow(net, stop_condition="tol", iter=3, friction_model="nikuradse", mode="hydraulics", transient=False, nonlinear_method="automatic", tol_p=1e-4, tol_v=1e-4) assert net.res_pump.deltap_bar.isin([0]).all() assert np.isclose(net.res_junction.loc[1, "p_bar"], net.res_junction.loc[2, "p_bar"])
def test_connectivity_hydraulic(create_test_net): """ :param create_test_net: :type create_test_net: :return: :rtype: """ net = copy.deepcopy(create_test_net) net.junction.in_service = True net.pipe.in_service = [True, False, False, True, True] pandapipes.create_fluid_from_lib(net, "water") pandapipes.pipeflow(net, iter=100, tol_p=1e-7, tol_v=1e-7, friction_model="nikuradse") assert np.all(np.isnan(net.res_junction.loc[[2, 5, 6], :].values)) assert np.all(np.isnan(net.res_pipe.loc[[1, 2, 3], :].values)) assert not np.any(np.isnan(net.res_junction.loc[[0, 1, 3, 4], :].values)) assert not np.any( np.isnan(net.res_pipe.loc[ [0, 4], ["v_mean_m_per_s", "p_from_bar", "p_to_bar"]].values)) assert not np.any( np.isnan(net.res_sink.loc[[0, 2], "mdot_kg_per_s"].values)) assert np.all(np.isnan(net.res_sink.loc[[1, 3, 4], "mdot_kg_per_s"].values)) assert np.allclose(net.res_ext_grid.mdot_kg_per_s.sum(), -net.res_sink.mdot_kg_per_s.sum(), rtol=1e-10, atol=0) active_branches = get_lookup(net, "branch", "active") active_nodes = get_lookup(net, "node", "active") assert np.all(active_nodes == np.array( [True, True, False, True, True, False, False, False, False, True])) assert np.all(active_branches == np.array( [True, False, False, False, False, False, True, True, True, False]))
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 test_connectivity_hydraulic2(create_test_net): net = copy.deepcopy(create_test_net) net.junction.in_service = True net.pipe.in_service = True net.valve.opened = True pandapipes.create_fluid_from_lib(net, "water") pandapipes.pipeflow(net, iter=100, tol_p=1e-7, tol_v=1e-7, friction_model="nikuradse") pandapipes.create_junction(net, 1., 293.15) pandapipes.create_junction(net, 1., 293.15) j = pandapipes.create_junction(net, 1., 293.15) pandapipes.create_junction(net, 1, 293.15) pandapipes.create_junction(net, 1, 293.15) pandapipes.create_pipe_from_parameters(net, 0, j, 0.1, 0.1) pandapipes.create_sink(net, j, 0.1) pandapipes.pipeflow(net, iter=100, tol_p=1e-7, tol_v=1e-7, friction_model="nikuradse") active_branches = get_lookup(net, "branch", "active") active_nodes = get_lookup(net, "node", "active") assert np.all(active_nodes == np.array([ True, True, True, True, True, True, True, False, False, True, False, False, True, True, True ])) assert np.all(active_branches) assert not np.all( np.isnan(net.res_junction.loc[[0, 1, 2, 3, 4, 5, 9], :].values)) assert not np.all(np.isnan(net.res_pipe.values)) assert np.any(np.isnan(net.res_junction.loc[[7, 8, 10, 11], :].values))
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)
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_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_inservice_water(create_test_net): """ :param create_test_net: :type create_test_net: :return: :rtype: """ net = copy.deepcopy(create_test_net) pandapipes.create_fluid_from_lib(net, "water") pandapipes.pipeflow(net, iter=100, tol_p=1e-7, tol_v=1e-7, friction_model="nikuradse") assert np.all(np.isnan(net.res_pipe.loc[~net.pipe.in_service, :].values)) assert np.all(np.isnan(net.res_valve.loc[~net.valve.opened, :].values)) assert np.all( np.isnan(net.res_junction.loc[~net.junction.in_service, :].values)) oos_sinks = np.isin(net.sink.junction.values, net.junction.index[~net.junction.in_service]) \ | ~net.sink.in_service.values assert np.all(np.isnan(net.res_sink.loc[oos_sinks, :].values)) assert not any( np.isnan(net.res_pipe.v_mean_m_per_s.loc[net.pipe.in_service].values)) assert not any( np.isnan(net.res_valve.v_mean_m_per_s.loc[net.valve.opened].values)) assert not any( np.isnan(net.res_junction.p_bar.loc[net.junction.in_service].values)) assert not np.any(np.isnan(net.res_sink.loc[~oos_sinks, :].values)) assert np.allclose(net.res_ext_grid.mdot_kg_per_s.sum(), -net.res_sink.mdot_kg_per_s.sum(), rtol=1e-10, atol=0)
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_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 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 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 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)