Skip to content

Add sectors

Add basic EV loads and chargers to a PyPSA network.

add_carrier_if_missing(n, carrier_name)

Add a carrier to the network if it doesn't already exist.

Parameters:

Name Type Description Default
n Network

PyPSA network to modify.

 required
carrier_name str

Name of the carrier to add.

 required
Source code in workflow/scripts/add_sectors.py 
12
13
14
15
16
17
18
19
20
21
def add_carrier_if_missing(n: pypsa.Network, carrier_name: str):
    """Add a carrier to the network if it doesn't already exist.

    Args:
        n (pypsa.Network): PyPSA network to modify.
        carrier_name (str): Name of the carrier to add.
    """
    if carrier_name not in n.carriers.index:
        n.add("Carrier", carrier_name)
        logger.debug("Carrier '%s' added to network.", carrier_name)

attach_simple_ev(n, p_set, nodes, options, ev_type)

Attach electric vehicle demand and charging links to the PyPSA network.

This function implements an EV demand model with a fixed charging profile.
For each node, it creates: • an EV load bus, • a Load component representing the EV charging demand time series, • and a charger Link connecting the AC bus to the EV load bus.

Parameters:

Name Type Description Default
n Network

Network to modify in place.

 required
p_set DataFrame

EV charging demand time series (MW), indexed by snapshots and with nodes as columns.

 required
nodes Index

AC buses where EV components are added.

 required
options dict

EV configuration parameters, including: - annual_consumption (float): annual energy demand per vehicle (MWh/year) - charge_rate (float): charger power rating per vehicle (MW) - share_charger (float): fraction of vehicles that can charge simultaneously

 required
ev_type str

EV category for component naming (e.g., "passenger", "freight").

 required

Returns:

Name Type Description
None

Modifies the network in place.
Source code in workflow/scripts/add_sectors.py 
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
def attach_simple_ev(
    n: pypsa.Network, p_set: pd.DataFrame, nodes: pd.Index, options: dict, ev_type: str
):
    """
    Attach electric vehicle demand and charging links to the PyPSA network.

    This function implements an EV demand model with a fixed charging profile.  
    For each node, it creates:
      • an EV load bus,
      • a Load component representing the EV charging demand time series,
      • and a charger Link connecting the AC bus to the EV load bus.

    Args:
        n (pypsa.Network): Network to modify in place.
        p_set (pd.DataFrame): EV charging demand time series (MW), indexed by
            snapshots and with nodes as columns.
        nodes (pd.Index): AC buses where EV components are added.
        options (dict): EV configuration parameters, including:
            - annual_consumption (float): annual energy demand per vehicle (MWh/year)
            - charge_rate (float): charger power rating per vehicle (MW)
            - share_charger (float): fraction of vehicles that can charge simultaneously
        ev_type (str): EV category for component naming (e.g., "passenger", "freight").

    Returns:
        None: Modifies the network in place.
    """

    total_energy = p_set.sum().sum()
    total_number_evs = total_energy / max(options["annual_consumption"], 1e-6)
    node_ratio = p_set.sum() / max(total_energy, 1e-6)
    number_evs = node_ratio * total_number_evs
    charge_power = (
        number_evs * options["charge_rate"] * options["share_charger"]
    ).clip(lower=0.001)

    logger.info("EV %s: %s vehicles (direct charging)", ev_type, f"{int(total_number_evs):,}")
    logger.debug(
        "EV %s: total energy %.2f MWh, annual consumption %.3f MWh/veh",
        ev_type,
        total_energy,
        options["annual_consumption"],
    )
    logger.debug(
        "EV %s: sample node power (MW)\n%s",
        ev_type,
        charge_power.head().to_string(),
    )

    ev_load_carrier = f"EV_{ev_type}_load"
    add_carrier_if_missing(n, ev_load_carrier)
    ev_load_bus = nodes + f" EV_{ev_type}_load"
    n.add("Bus", nodes, suffix=f" EV_{ev_type}_load", carrier=ev_load_carrier)

    n.add(
        "Load",
        nodes,
        suffix=f" EV_{ev_type}_load",
        bus=ev_load_bus,
        carrier=ev_load_carrier,
        p_set=p_set.loc[n.snapshots, nodes],
    )

    charger_carrier = f"EV_{ev_type}_charger"
    add_carrier_if_missing(n, charger_carrier)
    n.add(
        "Link",
        nodes,
        suffix=f" EV_{ev_type}_charger",
        bus0=nodes,
        bus1=ev_load_bus,
        carrier=charger_carrier,
        p_nom=charge_power,
        efficiency=1.0,
    )