powerflow/powerflow.h

Go to the documentation of this file.

#ifndef _POWERFLOW_H
#define _POWERFLOW_H

#include "gridlabd.h"
#include "solver_nr.h"

#ifdef _POWERFLOW_CPP
#define GLOBAL
#define INIT(A) = (A)
#else
#define GLOBAL extern
#define INIT(A)
#endif

#ifdef _DEBUG
void print_matrix(complex mat[3][3]);
#endif

#define GETOBJECT(obj) ((OBJECT *) obj - 1)
#define IMPORT_CLASS(name) extern CLASS *name##_class

//Deltamode use
#define TSNVRDBL 9223372036854775808.0

typedef enum {SM_FBS=0, SM_GS=1, SM_NR=2} SOLVERMETHOD;     
typedef enum {MM_SUPERLU=0, MM_EXTERN=1} MATRIXSOLVERMETHOD;    
typedef enum {
    MD_NONE=0,          
    MD_ONCE=1,          
    MD_PERCALL=2,       
    MD_ALL=3            
} MATRIXDUMPMETHOD;

typedef enum {
    LS_OPEN=0,          
    LS_CLOSED=1         
} LINESTATUS;   //Line/link status - made at powerflow level for reusability

typedef enum {
    IRM_NONE=0,         
    IRM_UNDEFINED=1,    
    IRM_TRAPEZOIDAL=2,  
    IRM_BACKEULER=3     
} INRUSHINTMETHOD;  //Selection of the integration method for the inrush calculations

typedef enum {
    FMM_NONE=0,     
    FMM_SIMPLE=1,   
    FMM_PLL=2       
} FREQMEASDEFAULT;  //Base default method for frequency measurement

//Structure to hold external LU solver calls
typedef struct s_ext_fxn {
    void *dllLink;
    void *ext_init;
    void *ext_alloc;
    void *ext_solve;
    void *ext_destroy;
} EXT_LU_FXN_CALLS;

GLOBAL char256 LUSolverName INIT("");               
GLOBAL EXT_LU_FXN_CALLS LUSolverFcns;               
GLOBAL SOLVERMETHOD solver_method INIT(SM_FBS);     
GLOBAL char256 MDFileName INIT("");                 
GLOBAL MATRIXDUMPMETHOD NRMatDumpMethod INIT(MD_NONE);  
GLOBAL bool NRMatReferences INIT(false);            
GLOBAL bool use_line_cap INIT(false);               
GLOBAL bool use_link_limits INIT(true);             
GLOBAL MATRIXSOLVERMETHOD matrix_solver_method INIT(MM_SUPERLU);    
GLOBAL unsigned int NR_bus_count INIT(0);           
GLOBAL unsigned int NR_branch_count INIT(0);        
GLOBAL BUSDATA *NR_busdata INIT(NULL);              
GLOBAL BRANCHDATA *NR_branchdata INIT(NULL);        
GLOBAL NR_SOLVER_STRUCT NR_powerflow;               
GLOBAL int NR_islands_detected INIT(0);             
GLOBAL bool NR_island_fail_method INIT(false);      
GLOBAL bool NR_solver_working INIT(false);          
GLOBAL int NR_curr_bus INIT(-1);                    
GLOBAL int NR_curr_branch INIT(-1);                 
GLOBAL int64 NR_iteration_limit INIT(500);          
GLOBAL bool NR_dyn_first_run INIT(true);            
GLOBAL bool NR_admit_change INIT(true);             
GLOBAL int NR_superLU_procs INIT(1);                
GLOBAL TIMESTAMP NR_retval INIT(TS_NEVER);          
GLOBAL OBJECT *NR_swing_bus INIT(NULL);             
GLOBAL int NR_swing_bus_reference INIT(-1);         
GLOBAL int64 NR_delta_iteration_limit INIT(10);     
GLOBAL bool FBS_swing_set INIT(false);              
GLOBAL bool show_matrix_values INIT(false);         
GLOBAL double primary_voltage_ratio INIT(60.0);     
GLOBAL double nominal_frequency INIT(60.0);         
GLOBAL double warning_underfrequency INIT(55.0);    
GLOBAL double warning_overfrequency INIT(65.0);     
GLOBAL double warning_undervoltage INIT(0.8);       
GLOBAL double warning_overvoltage INIT(1.2);        
GLOBAL double warning_voltageangle INIT(2.0);       
GLOBAL bool require_voltage_control INIT(false);    
GLOBAL double geographic_degree INIT(0.0);          
GLOBAL complex fault_Z INIT(complex(1e-6,0));       
GLOBAL complex ground_Z INIT(complex(1e-6,0));      
GLOBAL double default_maximum_voltage_error INIT(1e-6); 
GLOBAL double default_maximum_power_error INIT(0.0001); 
GLOBAL OBJECT *restoration_object INIT(NULL);       
GLOBAL OBJECT *fault_check_object INIT(NULL);       
GLOBAL bool fault_check_override_mode INIT(false);  
GLOBAL bool meshed_fault_checking_enabled INIT(false);  /*** fault_check object flag for possible meshing -- adjusts how reliability-related code runs */
GLOBAL bool restoration_checks_active INIT(false);  /***< Overall flag for when reconfigurations are occurring - special actions in devices */
GLOBAL double reliability_metrics_recloser_counts INIT(0.0);    /***< Recloser counts for momentary outages for power_metrics -- primarily to do backwards compatibility*/

GLOBAL bool enable_subsecond_models INIT(false);        /* normally not operating in delta mode */
GLOBAL bool all_powerflow_delta INIT(false);            /* Flag to make all powerflow objects participate in deltamode -- otherwise is individually flagged per object */
GLOBAL FREQMEASDEFAULT all_powerflow_freq_measure_method INIT(FMM_NONE);        /* Flag to enable all capable powerflow objects to do frequency measurements */
GLOBAL unsigned long deltamode_timestep INIT(10000000); /* deltamode timestep value - 10 ms timestep, at first - intermnal */
GLOBAL double deltamode_timestep_publish INIT(10000000.0); /* deltamode module-published 10 ms timestep, at first -- module property version, to be converted*/
GLOBAL OBJECT **delta_objects INIT(NULL);               /* Array pointer objects that need deltamode interupdate calls */
GLOBAL FUNCTIONADDR *delta_functions INIT(NULL);    /* Array pointer functions for objects that need deltamode interupdate calls */
GLOBAL FUNCTIONADDR *post_delta_functions INIT(NULL);       /* Array pointer functions for objects that need deltamode postupdate calls */
GLOBAL int pwr_object_count INIT(0);                /* deltamode object count */
GLOBAL int pwr_object_current INIT(-1);             /* Index of current deltamode object */
GLOBAL TIMESTAMP deltamode_starttime INIT(TS_NEVER);    /* Tracking variable for next desired instance of deltamode */
GLOBAL TIMESTAMP deltamode_endtime INIT(TS_NEVER);      /* Tracking variable to see when deltamode ended - so differential calculations don't get messed up */
GLOBAL double deltamode_endtime_dbl INIT(TSNVRDBL);     /* Tracking variable to see when deltamode ended - double valued for explicit movement calculations */
GLOBAL TIMESTAMP deltamode_supersec_endtime INIT(TS_NEVER); /* Tracking variable to indicate the "floored" time of detamode_endtime */
GLOBAL double current_frequency INIT(60.0);         
GLOBAL bool master_frequency_update INIT(false);    
GLOBAL bool enable_frequency_dependence INIT(false);    
GLOBAL double default_resistance INIT(1e-4);        
//In-rush deltamode stuff
GLOBAL bool enable_inrush_calculations INIT(false); 
GLOBAL INRUSHINTMETHOD default_inrush_integration_method INIT(IRM_TRAPEZOIDAL); 
GLOBAL double impedance_conversion_low_pu INIT(0.7);    
GLOBAL double deltatimestep_running INIT(-1.0);         
//Mesh fault current stuff
GLOBAL bool enable_mesh_fault_current INIT(false);  
//Market-based item -- moved from triplex_meter
GLOBAL char1024 market_price_name INIT("current_market.clearing_price");

// Deltamode stuff
void schedule_deltamode_start(TIMESTAMP tstart);    /* Anticipated time for a deltamode start, even if it is now */

/* used by many powerflow enums */
#define UNKNOWN 0
#define ROUNDOFF 1e-6           // numerical accuracy for zero in float comparisons


#include "powerflow_object.h"

#endif // _POWERFLOW_H

---

GridLAB-D™ Version 4.1

An open-source smart grid simulator created by PNNL for the US Department of Energy Office of Electricity