minimal_example_lopf.py

# -*- coding: utf-8 -*-
## Minimal 3-node example of PyPSA linear optimal power flow
#
# Available as a Jupyter notebook at <https://pypsa.readthedocs.io/en/latest/examples/minimal_example_lopf.ipynb>.

import numpy as np

import pypsa

network = pypsa.Network()

# add three buses
for i in range(3):
    network.add("Bus", "My bus {}".format(i))

print(network.buses)

# add three lines in a ring
for i in range(3):
    network.add(
        "Line",
        "My line {}".format(i),
        bus0="My bus {}".format(i),
        bus1="My bus {}".format((i + 1) % 3),
        x=0.0001,
        s_nom=60,
    )

print(network.lines)

# add a generator at bus 0
network.add("Generator", "My gen 0", bus="My bus 0", p_nom=100, marginal_cost=50)

# add a generator at bus 1
network.add("Generator", "My gen 1", bus="My bus 1", p_nom=100, marginal_cost=25)

print(network.generators)

print(network.generators.p_set)

# add a load at bus 2
network.add("Load", "My load", bus="My bus 2", p_set=100)

print(network.loads)

print(network.loads.p_set)

# Do a linear OPF


def my_f(network, snapshots):
    print(snapshots)


network.lopf(extra_functionality=my_f)

# Cheap generator 1 cannot be fully dispatched because of network constraints,
# so expensive generator 0 also has to dispatch
print(network.generators_t.p)

# network flows
print(network.lines_t.p0)

# Line 1 is congested
print(abs(network.lines_t.p0) / network.lines.s_nom)

# Power flows towards lower voltage angles
print(network.buses_t.v_ang * 180 / np.pi)

# In linear approximation, all voltage magnitudes are nominal, i.e. 1 per unit
print(network.buses_t.v_mag_pu)

# At bus 2 the price is set above any marginal generation costs in the model, because to dispatch to
# it from expensive generator 0, also some dispatch from cheap generator 1 has to be substituted from generator0
# to avoid overloading line 1.
print(network.buses_t.marginal_price)