import numpy as np
from pydgrid.pydgrid import grid
from pydgrid.pf import pf_eval,time_serie

Kersting book example 6.1

Data

data = {
"lines":[
        {"bus_j": "Bus_1",  "bus_k": "Bus_2",  "code": "Kersting", "m": 1609.34}
        ],
"buses":[
            {"bus": "Bus_1",  "pos_x": 10, "pos_y":  0, "units": "m", "U_kV":12.47},
            {"bus": "Bus_2",  "pos_x": 200, "pos_y": 0, "units": "m", "U_kV":12.47}
            ],
"grid_formers":[
            {"bus": "Bus_1","bus_nodes": [1, 2, 3],
                    "kV": [7.53869, 7.45125, 7.48512],
                    "deg": [1.57248, -118.30047, 121.93184]
            }
            ],
"loads":[
            {"bus": "Bus_2" , "kVA": 6000.0, "pf": 0.9,"type":"3P"}
        ],
"line_codes":
            {"Kersting":
            {"R":[
                                    [0.8667, 0.2955, 0.2907],
                                    [0.2955, 0.8837, 0.2992],
                                    [0.2907, 0.2992, 0.8741]
                                    ],
             "X":[
                        [2.0417,0.9502, 0.7290],
          [0.9502,1.9852, 0.8023],
          [0.7290,0.8023, 2.0172]
                              ],
     "B_mu":[
                                     [10.7409, -3.4777, -1.3322],
                                     [-3.4777, 11.3208, -2.2140],
                                     [ -1.3322, -2.2140, 10.2104]
                             ],
"unit":"miles"
            }
            }
}

Execute power flow

grid_1 = grid()
grid_1.read(data)  # Load data
grid_1.pf_solver = 2
grid_1.pf()  # solve power flow

Graph with obtained results

from pydgrid.plot_bokeh import plot_results
plot_results(grid_1)

Loading BokehJS ...

Figure(

id = '1003', …)

above = [],

aspect_scale = 1,

background_fill_alpha = {'value': 1.0},

background_fill_color = {'value': '#ffffff'},

below = [LinearAxis(id='1013', ...)],

border_fill_alpha = {'value': 1.0},

border_fill_color = {'value': '#ffffff'},

css_classes = [],

disabled = False,

extra_x_ranges = {},

extra_y_ranges = {},

h_symmetry = True,

height = None,

hidpi = True,

js_event_callbacks = {},

js_property_callbacks = {},

left = [LinearAxis(id='1018', ...)],

lod_factor = 10,

lod_interval = 300,

lod_threshold = 2000,

lod_timeout = 500,

match_aspect = False,

min_border = 5,

min_border_bottom = None,

min_border_left = None,

min_border_right = None,

min_border_top = None,

name = None,

outline_line_alpha = {'value': 1.0},

outline_line_cap = 'butt',

outline_line_color = {'value': '#e5e5e5'},

outline_line_dash = [],

outline_line_dash_offset = 0,

outline_line_join = 'bevel',

outline_line_width = {'value': 1},

output_backend = 'canvas',

plot_height = 400,

plot_width = 600,

renderers = [LinearAxis(id='1013', ...), Grid(id='1017', ...), LinearAxis(id='1018', ...), Grid(id='1022', ...), BoxAnnotation(id='1031', ...), GlyphRenderer(id='1042', ...), GlyphRenderer(id='1048', ...), GlyphRenderer(id='1054', ...)],

right = [],

sizing_mode = 'fixed',

subscribed_events = [],

tags = [],

title = Title(id='1002', ...),

title_location = 'above',

toolbar = Toolbar(id='1029', ...),

toolbar_location = 'right',

toolbar_sticky = True,

v_symmetry = False,

width = None,

x_range = DataRange1d(id='1005', ...),

x_scale = LinearScale(id='1009', ...),

y_range = DataRange1d(id='1007', ...),

y_scale = LinearScale(id='1011', ...))

Get element transfomers results

grid_1.buses

[{'bus': 'Bus_1',
  'pos_x': 10,
  'pos_y': 0,
  'units': 'm',
  'U_kV': 12.47,
  'N_nodes': 3,
  'v_an': 7538.69,
  'v_bn': 7451.25,
  'v_cn': 7485.12,
  'v_ng': 0.0,
  'deg_an': 1.57248,
  'deg_bn': -118.30046999999999,
  'deg_cn': 121.93184,
  'deg_ng': 0.0,
  'v_ab': 12973.42498856987,
  'v_bc': 12920.120310382239,
  'v_ca': 13034.52171783379,
  'p_a': 1858935.7393203387,
  'p_b': 1835419.9515336645,
  'p_c': 1839830.0481158614,
  'q_a': 963491.4173624129,
  'q_b': 956338.0192537266,
  'q_c': 968270.9001076091},
 {'bus': 'Bus_2',
  'pos_x': 200,
  'pos_y': 0,
  'units': 'm',
  'U_kV': 12.47,
  'N_nodes': 3,
  'v_an': 7199.551320077657,
  'v_bn': 7199.557391530307,
  'v_cn': 7199.55602482274,
  'v_ng': 0.0,
  'deg_an': -3.170444377827802e-05,
  'deg_bn': -120.00002917542324,
  'deg_cn': 119.9999628259816,
  'deg_ng': 0.0,
  'v_ab': 12469.993777213469,
  'v_bc': 12469.998513070677,
  'v_ca': 12469.992408860844,
  'p_a': -1799998.8481163539,
  'p_b': -1800000.327593055,
  'p_c': -1800000.107596966,
  'q_a': -871778.0011666914,
  'q_b': -871778.8157972179,
  'q_c': -871778.399022273}]