powerio-dist 0.11.7

Multiconductor distribution network model and lossless converters for OpenDSS, PMD JSON, BMOPF JSON, and distribution CIM.
Documentation 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
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

use super::*;

fn geometry() -> Geometry {
    Geometry {
        wires: (0..4)
            .map(|i| {
                let j = f64::from(i);
                super::super::Wire {
                    x: j * 0.1,
                    height: -1.0 - 0.05 * f64::from(i % 2),
                    radius: 0.008 + j * 0.0005,
                    capacitance_radius: 0.01 + j * 0.0005,
                    gmr: 0.006 + j * 0.0003,
                    rdc: 0.0002 * (1.0 + j / 5.0),
                    rac: 0.00021 * (1.0 + j / 5.0),
                    cable: (i < 3).then_some(Data {
                        k: (13 + i) as usize,
                        diastrand: 0.0015,
                        gmrstrand: 0.0006,
                        rstrand: 0.0008,
                        epsr: 2.3 + j * 0.1,
                        inslayer: 0.004,
                        diains: 0.026,
                        diacable: 0.031,
                    }),
                }
            })
            .collect(),
        phases: 3,
        reduce: true,
        reduction_order: None,
        earth_model: EarthModel::Carson,
        rho: 37.0,
    }
}

#[test]
fn cable_and_bare_neutral_reduction_matches_pinned_native_values() {
    // DSS C-API 87d85c2, dss_cn_geometry_oracle.py carson-60.0-4-3-True.
    let matrices = calculate(&geometry(), 60.0).unwrap();
    let r = [
        0.000_265_503_608_882_421_4,
        3.519_501_778_096_733e-6,
        6.498_786_398_968_169e-7,
        3.519_501_778_096_733e-6,
        0.000_302_151_139_117_561_83,
        2.374_709_548_094_71e-6,
        6.498_786_398_968_169e-7,
        2.374_709_548_094_71e-6,
        0.000_341_845_252_252_206_93,
    ];
    let x = [
        8.428_299_081_358_696e-5,
        -5.959_678_010_096_394e-6,
        -3.215_202_146_906_68e-6,
        -5.959_678_010_096_394e-6,
        8.041_560_706_820_995e-5,
        -5.162_644_883_651_848e-6,
        -3.215_202_146_906_68e-6,
        -5.162_644_883_651_848e-6,
        7.349_568_921_535_918e-5,
    ];
    for (actual, expected) in matrices
        .r
        .iter()
        .flatten()
        .zip(r)
        .chain(matrices.x.iter().flatten().zip(x))
    {
        assert!((actual - expected).abs() < 1e-15);
    }
    assert!((matrices.b[0][0] - 1.311_763_202_020_904_5e-7).abs() < 1e-17);
    assert_eq!(matrices.b[0][1].to_bits(), 0.0_f64.to_bits());
    let mut unreduced = geometry();
    unreduced.reduce = false;
    let matrices = calculate(&unreduced, 60.0).unwrap();
    assert_eq!(matrices.r.len(), 4);
    assert_eq!(matrices.b[3], vec![0.0; 4]);
}

#[test]
fn incomplete_and_mixed_phase_cables_are_explicitly_refused() {
    let mut input = geometry();
    input.wires[0].cable.as_mut().unwrap().diains = 0.001;
    assert!(calculate(&input, 60.0).unwrap_err().contains("incomplete"));
    input = geometry();
    input.wires[0].cable = None;
    assert!(
        calculate(&input, 60.0)
            .unwrap_err()
            .contains("all phase positions")
    );
    input = geometry();
    input.wires[0].height = f64::NAN;
    assert!(calculate(&input, 60.0).is_err());
}