def pipeflow(net, sol_vec=None, **kwargs):
"""
The main method used to start the solver to calculate the veatity, pressure and temperature\
distribution for a given net. Different options can be entered for \\**kwargs, which control\
the solver behaviour (see function init constants for more information).
:param net: The pandapipes net for which to perform the pipeflow
:type net: pandapipesNet
:param sol_vec:
:type sol_vec:
:param kwargs: A list of options controlling the solver behaviour
:return: No output
EXAMPLE:
pipeflow(net, mode="hydraulic")
"""
local_params = dict(locals())
# Inputs & initialization of variables
# ------------------------------------------------------------------------------------------
# Init physical constants and options
init_options(net, local_params)
create_lookups(net, NodeComponent, BranchComponent, BranchWInternalsComponent)
node_pit, branch_pit = initialize_pit(net, Junction.table_name(),
NodeComponent, NodeElementComponent,
BranchComponent, BranchWInternalsComponent)
calculation_mode = get_net_option(net, "mode")
if get_net_option(net, "check_connectivity"):
nodes_connected, branches_connected = check_connectivity(
net, branch_pit, node_pit, check_heat=calculation_mode in ["heat", "all"])
else:
nodes_connected = node_pit[:, ACTIVE_ND].astype(np.bool)
branches_connected = branch_pit[:, ACTIVE_BR].astype(np.bool)
reduce_pit(net, node_pit, branch_pit, nodes_connected, branches_connected)
if calculation_mode == "hydraulics":
niter = hydraulics(net)
elif calculation_mode == "heat":
if net.user_pf_options["hyd_flag"]:
node_pit = net["_active_pit"]["node"]
node_pit[:, PINIT] = sol_vec[:len(node_pit)]
branch_pit = net["_active_pit"]["branch"]
branch_pit[:, VINIT] = sol_vec[len(node_pit):]
niter = heat_transfer(net)
else:
logger.warning("Converged flag not set. Make sure that hydraulic calculation results "
"are available.")
elif calculation_mode == "all":
niter = hydraulics(net)
niter = heat_transfer(net)
else:
logger.warning("No proper calculation mode chosen.")
extract_results_active_pit(net, node_pit, branch_pit, nodes_connected, branches_connected)
extract_all_results(net, Junction.table_name())
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)
pandapipes.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"][:, 0] == 0)
pipe_p_data = pipe_results["PINIT"][pipe_p_data_idx, 1]
pipe_v_data = pipe_results["VINIT"][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_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)