powerflow/series_reactor.cpp

Go to the documentation of this file.

#include <stdlib.h>
#include <stdio.h>
#include <errno.h>
#include <math.h>

#include "series_reactor.h"
#include "node.h"

// series_reactor CLASS FUNCTIONS
CLASS* series_reactor::oclass = NULL;
CLASS* series_reactor::pclass = NULL;

series_reactor::series_reactor(MODULE *mod) : link(mod)
{
    if(oclass == NULL)
    {
        pclass = link::oclass;

        oclass = gl_register_class(mod,"series_reactor",sizeof(series_reactor),PC_PRETOPDOWN|PC_BOTTOMUP|PC_POSTTOPDOWN|PC_UNSAFE_OVERRIDE_OMIT);
        if(oclass == NULL)
            GL_THROW("unable to register object class implemented by %s",__FILE__);

        if(gl_publish_variable(oclass,
            PT_INHERIT, "link",
            PT_complex, "phase_A_impedance[Ohm]",PADDR(phase_A_impedance),PT_DESCRIPTION,"Series impedance of reactor on phase A",
            PT_double, "phase_A_resistance[Ohm]",PADDR(phase_A_impedance.Re()),PT_DESCRIPTION,"Resistive portion of phase A's impedance",
            PT_double, "phase_A_reactance[Ohm]",PADDR(phase_A_impedance.Im()),PT_DESCRIPTION,"Reactive portion of phase A's impedance",
            PT_complex, "phase_B_impedance[Ohm]",PADDR(phase_B_impedance),PT_DESCRIPTION,"Series impedance of reactor on phase B",
            PT_double, "phase_B_resistance[Ohm]",PADDR(phase_B_impedance.Re()),PT_DESCRIPTION,"Resistive portion of phase B's impedance",
            PT_double, "phase_B_reactance[Ohm]",PADDR(phase_B_impedance.Im()),PT_DESCRIPTION,"Reactive portion of phase B's impedance",
            PT_complex, "phase_C_impedance[Ohm]",PADDR(phase_C_impedance),PT_DESCRIPTION,"Series impedance of reactor on phase C",
            PT_double, "phase_C_resistance[Ohm]",PADDR(phase_C_impedance.Re()),PT_DESCRIPTION,"Resistive portion of phase C's impedance",
            PT_double, "phase_C_reactance[Ohm]",PADDR(phase_C_impedance.Im()),PT_DESCRIPTION,"Reactive portion of phase C's impedance",
            PT_double, "rated_current_limit[A]",PADDR(rated_current_limit),PT_DESCRIPTION,"Rated current limit for the reactor",
            NULL) < 1) GL_THROW("unable to publish properties in %s",__FILE__);
    }
}

int series_reactor::isa(char *classname)
{
    return strcmp(classname,"series_reactor")==0 || link::isa(classname);
}

int series_reactor::create()
{
    int result = link::create();
    phase_A_impedance = phase_B_impedance = phase_C_impedance = 0.0;

    return result;
}

int series_reactor::init(OBJECT *parent)
{
    int result = link::init(parent);

    a_mat[0][0] = d_mat[0][0] = A_mat[0][0] = (has_phase(PHASE_A)) ? 1.0 : 0.0;
    a_mat[1][1] = d_mat[1][1] = A_mat[1][1] = (has_phase(PHASE_B)) ? 1.0 : 0.0;
    a_mat[2][2] = d_mat[2][2] = A_mat[2][2] = (has_phase(PHASE_C)) ? 1.0 : 0.0;

    if (solver_method==SM_FBS)
    {
        c_mat[0][0] = 0.0;
        c_mat[1][1] = 0.0;
        c_mat[2][2] = 0.0;

        //Pass in impedance values
        b_mat[0][0] = B_mat[0][0] = phase_A_impedance;
        b_mat[1][1] = B_mat[1][1] = phase_B_impedance;
        b_mat[2][2] = B_mat[2][2] = phase_C_impedance;

    }
    else
    {
        //Flag it as special (we'll forgo inversion processes on this)
        SpecialLnk = SWITCH;

        //Initialize off-diagonals just in case
        From_Y[0][1] = From_Y[0][2] = From_Y[1][0] = 0.0;
        From_Y[1][2] = From_Y[2][0] = From_Y[2][1] = 0.0;

        //See if it has a particular phase, if so populate it.  If not
        //and has zero impedance, put a "default" value in its place
        if (has_phase(PHASE_A))
        {
            if (phase_A_impedance==0.0)
                From_Y[0][0] = complex(1e4,1e4);
            else
                From_Y[0][0] = complex(1.0,0.0)/phase_A_impedance;
        }
        else
            From_Y[0][0] = 0.0;     //Should already be 0, but let's be paranoid

        if (has_phase(PHASE_B))
        {
            if (phase_B_impedance==0.0)
                From_Y[1][1] = complex(1e4,1e4);
            else
                From_Y[1][1] = complex(1.0,0.0)/phase_B_impedance;
        }
        else
            From_Y[1][1] = 0.0;     //Should already be 0, but let's be paranoid

        if (has_phase(PHASE_C))
        {
            if (phase_C_impedance==0.0)
                From_Y[2][2] = complex(1e4,1e4);
            else
                From_Y[2][2] = complex(1.0,0.0)/phase_C_impedance;
        }
        else
            From_Y[2][2] = 0.0;     //Should already be 0, but let's be paranoid
    }

    return result;
}

// IMPLEMENTATION OF CORE LINKAGE: series_reactor

EXPORT int commit_series_reactor(OBJECT *obj)
{
    if (solver_method==SM_FBS)
    {
        series_reactor *plink = OBJECTDATA(obj,series_reactor);
        plink->calculate_power();
    }
    return 1;
}
EXPORT int create_series_reactor(OBJECT **obj, OBJECT *parent)
{
    try
    {
        *obj = gl_create_object(series_reactor::oclass);
        if (*obj!=NULL)
        {
            series_reactor *my = OBJECTDATA(*obj,series_reactor);
            gl_set_parent(*obj,parent);
            return my->create();
        }
    }
    catch (const char *msg)
    {
        gl_error("create_series_reactor: %s", msg);
    }
    return 0;
}

EXPORT int init_series_reactor(OBJECT *obj)
{
    series_reactor *my = OBJECTDATA(obj,series_reactor);
    try {
        return my->init(obj->parent);
    }
    catch (const char *msg)
    {
        GL_THROW("%s (series_reactor:%d): %s", my->get_name(), my->get_id(), msg);
        return 0;
    }
}

EXPORT TIMESTAMP sync_series_reactor(OBJECT *obj, TIMESTAMP t0, PASSCONFIG pass)
{
    series_reactor *pObj = OBJECTDATA(obj,series_reactor);
    try {
        TIMESTAMP t1 = TS_NEVER;
        switch (pass) {
        case PC_PRETOPDOWN:
            return pObj->presync(t0);
        case PC_BOTTOMUP:
            return pObj->sync(t0);
        case PC_POSTTOPDOWN:
            t1 = pObj->postsync(t0);
            obj->clock = t0;
            return t1;
        default:
            throw "invalid pass request";
        }
    } catch (const char *error) {
        GL_THROW("%s (series_reactor:%d): %s", pObj->get_name(), pObj->get_id(), error);
        return 0;
    } catch (...) {
        GL_THROW("%s (series_reactor:%d): %s", pObj->get_name(), pObj->get_id(), "unknown exception");
        return 0;
    }
}

EXPORT int isa_series_reactor(OBJECT *obj, char *classname)
{
    return OBJECTDATA(obj,series_reactor)->isa(classname);
}

---

GridLAB-D^TM Version 2.0

An open-source project initiated by the US Department of Energy