powerflow/solver_nr.h

/* $Id
 * Newton-Raphson solver
 */

#ifndef _SOLVER_NR
#define _SOLVER_NR

#include "complex.h"
#include "object.h"

typedef struct  {
    int type;               
    unsigned char phases;   
    unsigned char origphases;   
    set *busflag;           
    complex *V;             
    complex *S;             
    complex *Y;             
    complex *I;             
    complex *prerot_I;      
    complex *S_dy;          
    complex *Y_dy;          
    complex *I_dy;          
    complex *full_Y;        
    complex *full_Y_all;    
    complex *full_Y_load;   
    complex *extra_var;     
    complex *house_var;     
    int *Link_Table;        
    unsigned int Link_Table_Size;   
    double PL[3];           
    double QL[3];           
    bool *dynamics_enabled; 
    bool swing_functions_enabled;   
    complex *PGenTotal;     
    complex *DynCurrent;    
    complex *BusHistTerm;   
    complex *BusSatTerm;    
    double volt_base;       
    double mva_base;        
    double Jacob_A[3];      // Element a in equation (37), which is used to update the Jacobian matrix at each iteration
    double Jacob_B[3];      // Element b in equation (38), which is used to update the Jacobian matrix at each iteration
    double Jacob_C[3];      // Element c in equation (39), which is used to update the Jacobian matrix at each iteration
    double Jacob_D[3];      // Element d in equation (40), which is used to update the Jacobian matrix at each iteration
    unsigned int Matrix_Loc;// Starting index of this object's place in all matrices/equations
    double max_volt_error;  
    char *name;             
    OBJECT *obj;            
    FUNCTIONADDR ExtraCurrentInjFunc;   
    OBJECT *ExtraCurrentInjFuncObject;  
    int island_number;      
} BUSDATA;

typedef struct {
    complex *Yfrom;         
    complex *Yto;           
    complex *YSfrom;        
    complex *YSto;          
    unsigned char phases;   
    unsigned char origphases;   
    unsigned char faultphases;  
    int from;               
    int to;                 
    int fault_link_below;    
    enumeration *status;    
    unsigned char lnk_type; 
    double v_ratio;         
    char *name;             
    OBJECT *obj;            
    complex *If_from;       
    complex *If_to;         
    FUNCTIONADDR limit_check;   
    FUNCTIONADDR ExtraDeltaModeFunc;    
    int island_number;      
} BRANCHDATA;

typedef struct Y_NR{
    int row_ind;  
    int col_ind;  
    double Y_value; 
} Y_NR;

typedef struct {
    int row_ind;  
    int col_ind;  
    complex Y[3][3]; 
    char size;      
} Bus_admit;

typedef struct {
    double *a_LU;
    double *rhs_LU;
    int *cols_LU;
    int *rows_LU;
} NR_SOLVER_VARS;

typedef enum {
    PF_NORMAL=0,    
    PF_DYNINIT=1,   
    PF_DYNCALC=2    
} NRSOLVERMODE;

// Sparse element
typedef struct SP_E {
    struct SP_E* next;
    int row_ind;  
    double value; 
} SP_E;

// Sparse matrix
typedef struct {
    SP_E** cols;
    SP_E* llheap;
    unsigned int llptr;
    unsigned int ncols;
} SPARSE;

typedef struct {
    double *deltaI_NR;                  
    unsigned int size_offdiag_PQ;       
    unsigned int size_diag_fixed;       
    unsigned int total_variables;       
    unsigned int size_diag_update;      
    unsigned int size_Amatrix;          
    unsigned int max_size_offdiag_PQ;   
    unsigned int max_size_diag_fixed;   
    unsigned int max_total_variables;   
    unsigned int max_size_diag_update;  
    unsigned int prev_m;                
    unsigned int index_count;           
    unsigned int bus_count;             
    unsigned int indexer;               
    bool NR_realloc_needed;             
    Y_NR *Y_offdiag_PQ;                 
    Y_NR *Y_diag_fixed;                 
    Y_NR *Y_diag_update;                
    SPARSE *Y_Amatrix;                  
    NR_SOLVER_VARS matrices_LU;         
    void *LU_solver_vars;               
    int64 iteration_count;              
    bool new_iteration_required;        
    bool swing_converged;               
    bool swing_is_a_swing;              
    bool SaturationMismatchPresent;     
    int solver_info;                    
    int64 return_code;                  
    double max_mismatch_converge;       
} NR_MATRIX_CONSTRUCTION;

typedef struct {
    NR_MATRIX_CONSTRUCTION *island_matrix_values;   
    Bus_admit *BA_diag;                 
} NR_SOLVER_STRUCT;

//Mesh-fault-related structure - passing information
typedef struct {
    complex *z_matrix;          
    int NodeRefNum;             
    int return_code;            
} NR_MESHFAULT_IMPEDANCE;

//Function prototypes for external solver interface
//void *ext_solver_init(void *ext_array);
//void ext_solver_alloc(void *ext_array, unsigned int rowcount, unsigned int colcount, bool admittance_change);
//int ext_solver_solve(void *ext_array, NR_SOLVER_VARS *system_info_vars, unsigned int rowcount, unsigned int colcount);
//void ext_solver_destroy(void *ext_array, bool new_iteration);

int64 solver_nr(unsigned int bus_count, BUSDATA *bus, unsigned int branch_count, BRANCHDATA *branch, NR_SOLVER_STRUCT *powerflow_values, NRSOLVERMODE powerflow_type , NR_MESHFAULT_IMPEDANCE *mesh_imped_vals, bool *bad_computations);
void compute_load_values(unsigned int bus_count, BUSDATA *bus, NR_SOLVER_STRUCT *powerflow_values, bool jacobian_pass, int island_number);

//Newton-Raphson solver array handlers
STATUS NR_array_structure_free(NR_SOLVER_STRUCT *struct_of_interest,int number_of_islands);     /* Handles freeing NR_SOLVER_STRUCT arrays */
STATUS NR_array_structure_allocate(NR_SOLVER_STRUCT *struct_of_interest,int number_of_islands); /* Allocates NR_SOLVER_STRUCT item for a given number of entries/islands */

#endif

---

GridLAB-D™ Version 4.1

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