Short circuits in MV line

In [12]:

import numpy as np
from pydgrid import grid
from pydgrid.plot_bokeh import plot_results

System

In [2]:

data = {
"transformers":[
        {"bus_j": "B0",  "bus_k": "B1",  "S_n_kVA": 250000.0, "U_j_kV":66.0, "U_k_kV":20.0,
        "R_cc_pu": 0.01, "X_cc_pu":0.08, "connection": "Dyn1",
        "conductors_1": 3, "conductors_2": 4}
        ],
"lines":[
        {"bus_j": "B1",  "bus_k": "B2",  "code": "mv_al_120", "m": 200.0}
        ],
"buses":[
        {"bus": "B0",  "pos_x": 0,  "pos_y": 0, "units": "m", "U_kV":66.0},
        {"bus": "B1",  "pos_x":30,  "pos_y": 0, "units": "m", "U_kV":20.0},
        {"bus": "B2",  "pos_x":200, "pos_y": 0, "units": "m", "U_kV":20.0}
        ],
"grid_formers":[
        {"bus": "B0",
        "bus_nodes": [1, 2, 3], "deg": [0, -120, -240],
        "kV": [38.105, 38.105, 38.105]}
        ],
"grid_feeders":[{"bus": "B2","bus_nodes": [1, 2, 3],
                "kW": [0, 0, 0], "kvar": [0,0,0],
                "kA": [0,0,0], "phi_deg":[30, 30, 30]}
                ],
    "line_codes":
        {"mv_al_50":  {"R1":0.8,    "X1":   0.148, "R0":0.8,   "X0":        0.148},
         "mv_al_95":  {"R1":0.403,  "X1":   0.129, "R0":0.403, "X0":        0.129},
         "mv_al_120": {"R1":0.321,  "X1":   0.123, "R0":0.321, "X0":        0.321},
         "mv_al_185": {"R1":0.209,  "X1":   0.113, "R0":0.209, "X0":        0.209},
         "mv_al_300": {"R1":0.128,  "X1":   0.105, "R0":0.128, "X0":        0.128}
         },
"shunts":[
        {"bus": "B1" , "R": 1e12, "X": 0.0, "bus_nodes": [1,0]},
        {"bus": "B1" , "R": 1e-8, "X": 0.0, "bus_nodes": [1,2]}, # applied fault
        {"bus": "B1" , "R": 1e-8, "X": 0.0, "bus_nodes": [2,3]}, # applied fault
        ]
}

Three phase short circuit

In [3]:

grid_1 = grid()
grid_1.read(data)  # Load data
grid_1.pf()  # solve power flow
p=plot_results(grid_1)

Phase ground fault (isolated neutral)

In [4]:

data['shunts'] = [
        {"bus": "B1" , "R": 1e12, "X": 0.0, "bus_nodes": [1,0]},
        {"bus": "B1" , "R": 1e-8, "X": 0.0, "bus_nodes": [1,0]}, # applied fault
]

In [5]:

grid_ph_g = grid()
grid_ph_g.read(data)  # Load data
grid_ph_g.pf()  # solve power flow
p=plot_results(grid_ph_g)

Phase ground fault (grounded neutral)

In [6]:

data['shunts'] = [
        {"bus": "B1" , "R": 1e-12, "X": 0.0, "bus_nodes": [4,0]},
        {"bus": "B1" , "R": 1e-12, "X": 0.0, "bus_nodes": [1,0]}, # applied fault
]

In [7]:

grid_ph_n = grid()
grid_ph_n.read(data)  # Load data
grid_ph_n.pf()  # solve power flow
p=plot_results(grid_ph_n)

Phase ground fault (grounded neutral with impedance)

In [10]:

data['shunts'] = [
        {"bus": "B1" , "R": 12.0, "X": 0.0, "bus_nodes": [4,0]},
        {"bus": "B1" , "R": 1e-12, "X": 0.0, "bus_nodes": [1,0]}, # applied fault
]

In [11]:

grid_ph_n = grid()
grid_ph_n.read(data)  # Load data
grid_ph_n.pf()  # solve power flow
p=plot_results(grid_ph_n)

Run this notebook in binder:

Documentation Status