コード例 #1
0
def test_add_fluid():
    net = pandapipes.create_empty_network()
    fluid_old = pandapipes.call_lib("air")

    with pytest.raises(AttributeError, match="no fluid"):
        pandapipes.get_fluid(net)

    _add_fluid_to_net(net, fluid_old)
    fluid_new = pandapipes.create_constant_fluid("arbitrary_gas2",
                                                 "gas",
                                                 density=2,
                                                 compressibility=2)
    _add_fluid_to_net(net, fluid_new, overwrite=False)
    assert pandapipes.get_fluid(net) == fluid_old

    _add_fluid_to_net(net, fluid_new)
    assert pandapipes.get_fluid(net) == fluid_new

    net["fluid"] = "Hello"

    _add_fluid_to_net(net, fluid_new, overwrite=False)
    assert pandapipes.get_fluid(net) == "Hello"

    _add_fluid_to_net(net, fluid_new)
    assert pandapipes.get_fluid(net) == fluid_new
コード例 #2
0
def test_add_fluid():
    net = pandapipes.create_empty_network()
    fluid_old = pandapipes.call_lib("air")

    try:
        pandapipes.get_fluid(net)
        assert False, "should'nt get here"
    except UserWarning:
        pass

    _add_fluid_to_net(net, fluid_old)
    fluid_new = pandapipes.create_constant_fluid("arbitrary_gas2",
                                                 "gas",
                                                 density=2,
                                                 compressibility=2)
    _add_fluid_to_net(net, fluid_new, overwrite=False)
    assert pandapipes.get_fluid(net) == fluid_old

    _add_fluid_to_net(net, fluid_new)
    assert pandapipes.get_fluid(net) == fluid_new

    net["fluid"] = "Hello"

    _add_fluid_to_net(net, fluid_new, overwrite=False)
    assert pandapipes.get_fluid(net) == "Hello"

    _add_fluid_to_net(net, fluid_new)
    assert pandapipes.get_fluid(net) == fluid_new
コード例 #3
0
def create_fluid_from_lib(net, name, overwrite=True):
    """
    Creates a fluid from library (if there is an entry) and sets net["fluid"] to this value.
    Currently existing fluids in the library are: "hgas", "lgas", "hydrogen", "water", "air".

    :param net: The net for which this fluid should be created
    :type net: pandapipesNet
    :param name: The name of the fluid that shall be extracted from the fluid lib
    :type name: str
    :param overwrite: Flag if a possibly existing fluid in the net shall be overwritten
    :type overwrite: bool, default True
    :return: No output

    :Example:
        >>> pp.create_fluid_from_lib(net, name="water")

    """
    _add_fluid_to_net(net, call_lib(name), overwrite=overwrite)
コード例 #4
0
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)
コード例 #5
0
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)