Working Example of Solving RLCircuit and Microgrid with Jacobians

CMake/FindSuiteSparse.cmake

@@ -69,11 +69,11 @@ Author(s):
 set(SUITESPARSE_MODULES
   amd
   colamd
-  klu)
+  klu
+  suitesparseconfig)
 
 find_library(SUITESPARSE_LIBRARY
   NAMES
-  suitesparseconfig
   ${SUITESPARSE_MODULES}
   PATHS
   ${SUITESPARSE_DIR} $ENV{SUITESPARSE_DIR} ${SUITESPARSE_ROOT_DIR}
@@ -97,6 +97,7 @@ find_path(SUITESPARSE_INCLUDE_DIR
   amd.h
   colamd.h
   klu.h
+  SuiteSparse_config.h
   PATHS
   ${SUITESPARSE_DIR} $ENV{SUITESPARSE_DIR} ${SUITESPARSE_ROOT_DIR} ${SUITESPARSE_LIBRARY_DIR}/..
   PATH_SUFFIXES

CMakeLists.txt

@@ -81,7 +81,7 @@ option(GRIDKIT_ENABLE_IPOPT "Enable Ipopt support" ON)
 option(GRIDKIT_ENABLE_SUNDIALS "Enable SUNDIALS support" ON)
 
 # Enable KLU
-option(GRIDKIT_ENABLE_SUNDIALS_SPARSE "Enable SUNDIALS sparse linear solvers" OFF)
+option(GRIDKIT_ENABLE_SUNDIALS_SPARSE "Enable SUNDIALS sparse linear solvers" ON)
 
 
 set(CMAKE_MACOSX_RPATH 1)
@@ -110,7 +110,7 @@ endif("${isSystemDir}" STREQUAL "-1")
 
 
 # TODO: Probably beter to set a debug interface target
-set(CMAKE_CXX_FLAGS_DEBUG "-Wall -O0 -g")
+set(CMAKE_CXX_FLAGS_DEBUG "-Wall -O0 -g -DDEBUG")
 
 set(CMAKE_CXX_STANDARD 17)
 
@@ -143,6 +143,9 @@ add_subdirectory(ComponentLib)
 # General Utilities and File IO
 add_subdirectory(Utilities)
 
+#Local Sparse matrix operations
+add_subdirectory(SparseMatrix)
+
 # Create solvers
 add_subdirectory(Solver)
 

CircuitGraph.hpp

@@ -0,0 +1,119 @@
+
+
+#include <iostream>
+#include <cmath>
+#include <vector>
+#include <unordered_set>
+#include <map>
+
+/**
+ * @brief A very basic hypergraph setup for circuit representation. This forms the hypergraph as a bipartite graph. Doesn't allow removing. Can only grab sets of connections to nodes
+ * 
+ * @todo should replace with something better and more efficent. Should replace with a libraries setup instead. This would allow fast and easy partitioning of circuits
+ * 
+ * @todo should replace N and E with Node and Component classes respectively
+ * 
+ * @tparam IdxT 
+ * @tparam Label 
+ */
+template <typename N, typename E>
+class CircuitGraph
+{
+private:
+	std::set<N> hypernodes;
+	std::set<E> hyperedges;
+	std::map<E, std::set<N>> edgestonodes;
+public:
+	CircuitGraph();
+	~CircuitGraph();
+	bool addHyperEdge(E he);
+	bool addHyperNode(N hn);
+	bool addConnection(N hn, E he);
+	std::set<N> getHyperEdgeConnections(E he);
+	size_t amountHyperNodes();
+	size_t amountHyperEdges();
+	void printBiPartiteGraph(bool verbose = false);
+};
+
+template <typename N, typename E>
+CircuitGraph<N, E>::CircuitGraph()
+{
+}
+
+template <typename N, typename E>
+CircuitGraph<N, E>::~CircuitGraph()
+{
+}
+
+template <typename N, typename E>
+bool  CircuitGraph<N, E>::addHyperNode(N hn)
+{
+	return this->hypernodes.insert(hn).second;
+}
+
+template <typename N, typename E>
+bool  CircuitGraph<N, E>::addHyperEdge(E he)
+{
+	return this->hyperedges.insert(he).second;
+}
+
+
+template <typename N, typename E>
+bool  CircuitGraph<N, E>::addConnection(N hn, E he)
+{
+	if(this->hyperedges.count(he) == 0 || this->hypernodes.count(hn) == 0)
+	{
+		return false;
+	} 
+	return this->edgestonodes[he].insert(hn).second;
+}
+
+
+template <typename N, typename E>
+std::set<N>  CircuitGraph<N, E>::getHyperEdgeConnections(E he)
+{
+	return this->edgestonodes[he];
+}
+
+
+template <typename N, typename E>
+size_t  CircuitGraph<N, E>::amountHyperNodes()
+{
+	return this->hypernodes.size();
+}
+
+
+template <typename N, typename E>
+size_t  CircuitGraph<N, E>::amountHyperEdges()
+{
+	return this->hyperedges.size();
+}
+
+/**
+ * @brief Printing
+ * 
+ * @todo need to add verbose printing for connections display
+ * 
+ * @tparam IdxT 
+ * @param verbose 
+ */
+
+template <typename N, typename E>
+void CircuitGraph<N, E>::printBiPartiteGraph(bool verbose)
+{
+
+	std::cout << "Amount of HyperNodes: " << this->amountHyperNodes() << std::endl;
+	std::cout << "Amount of HyperEdges: " << this->amountHyperEdges() << std::endl;
+	std::cout << "Connections per Edge:" << std::endl;
+	for (auto i : this->edgestonodes)
+	{
+		std::cout << i.first << " : {";
+		for (auto j : i.second){
+			std::cout << j << ", ";
+		}
+		std::cout << "}\n";
+
+	}
+
+
+}

ComponentLib/Bus/BusPV.cpp

@@ -74,7 +74,7 @@ namespace ModelLib {
  */
 template <class ScalarT, typename IdxT>
 BusPV<ScalarT, IdxT>::BusPV()
-  : BaseBus<ScalarT, IdxT>(0), V_(0.0), theta0_(0.0), Pg_(0.0)
+  : BaseBus<ScalarT, IdxT>(0), V_(0.0), theta0_(0.0)
 {
     //std::cout << "Create BusPV..." << std::endl;
     //std::cout << "Number of equations is " << size_ << std::endl;
@@ -83,9 +83,7 @@ BusPV<ScalarT, IdxT>::BusPV()
 }
 
 /*!
- * @brief BusPV constructor.
- *
- * This constructor sets initial values for voltage and phase angle.
+ * @brief Constructor for a PV bus
  *
  * @todo Arguments that should be passed to ModelEvaluatorImpl constructor:
  * - Number of equations = 1 (size_)
@@ -94,10 +92,10 @@ BusPV<ScalarT, IdxT>::BusPV()
  * - Number of optimization parameters = 0
  */
 template <class ScalarT, typename IdxT>
-BusPV<ScalarT, IdxT>::BusPV(ScalarT V, ScalarT theta0, ScalarT Pg)
-  : BaseBus<ScalarT, IdxT>(0), V_(V), theta0_(theta0), Pg_(Pg)
+BusPV<ScalarT, IdxT>::BusPV(ScalarT V, ScalarT theta0)
+  : BaseBus<ScalarT, IdxT>(0), V_(V), theta0_(theta0)
 {
-    //std::cout << "Create BusPV ..." << std::endl;
+    //std::cout << "Create BusPV..." << std::endl;
     //std::cout << "Number of equations is " << size_ << std::endl;
 
     size_ = 1;
@@ -171,8 +169,8 @@ template <class ScalarT, typename IdxT>
 int BusPV<ScalarT, IdxT>::evaluateResidual()
 {
     // std::cout << "Evaluating residual of a PV bus ...\n";
-    P() = Pg_;
-    Q() = 0.0;
+    P() = 0.0; // <-- Residual P
+    Q() = 0.0; // <-- Output Qg, the reactive power generator needs to supply
 
     return 0;
 }

ComponentLib/Bus/BusPV.hpp

@@ -91,7 +91,7 @@ namespace ModelLib
         using BusData = GridKit::PowerSystemData::BusData<real_type, IdxT>;
 
         BusPV();
-        BusPV(ScalarT V, ScalarT theta0, ScalarT P);
+        BusPV(ScalarT V, ScalarT theta0);
         BusPV(BusData& data);
         virtual ~BusPV();
 
@@ -197,10 +197,9 @@ namespace ModelLib
         }
 
     private:
-        ScalarT V_;
-        ScalarT theta0_;  ///< Default initial value for phase
-        ScalarT Pg_;      ///< Generator injection
-        ScalarT Q_;
+        ScalarT V_;      ///< Bus voltage magnitude
+        ScalarT theta0_; ///< Default initial value for phase
+        ScalarT Q_;      ///< Reactive power that generator needs to provide
 
         ScalarT VB_;
         ScalarT thetaB_;

ComponentLib/Bus/BusSlack.cpp

@@ -97,6 +97,8 @@ BusSlack<ScalarT, IdxT>::BusSlack(ScalarT V, ScalarT theta)
 {
     //std::cout << "Create BusSlack..." << std::endl;
     //std::cout << "Number of equations is " << size_ << std::endl;
+    P() = 0.0;
+    Q() = 0.0;
     size_ = 0;
 }
 
@@ -106,6 +108,8 @@ BusSlack<ScalarT, IdxT>::BusSlack(BusData& data)
 {
     //std::cout << "Create BusSlack..." << std::endl;
     //std::cout << "Number of equations is " << size_ << std::endl;
+    P() = 0.0;
+    Q() = 0.0;
     size_ = 0;
 }
 

ComponentLib/CMakeLists.txt

@@ -69,3 +69,5 @@ add_subdirectory(Generator4Governor)
 add_subdirectory(Generator4Param)
 add_subdirectory(Load)
 add_subdirectory(MiniGrid)
+add_subdirectory(PowerElectronicsComponents)
+

ComponentLib/Generator/GeneratorPV.hpp

@@ -124,14 +124,18 @@ namespace ModelLib
             return P_;
         }
 
+        /// @brief  Reactive power excess on PV bus
+        /// @return reference to negative PV generator reactive power
         virtual ScalarT& Q()
         {
-            return bus_->Q();
+            return (bus_->Q());
         }
 
+        /// @brief  Reactive power excess on PV bus
+        /// @return const reference to negative PV generator reactive power
         virtual const ScalarT& Q() const
         {
-            return bus_->Q();
+            return (bus_->Q());
         }
 
     private:

ComponentLib/PowerElectronicsComponents/CMakeLists.txt

@@ -0,0 +1,14 @@
+
+
+add_subdirectory(Capacitor)
+add_subdirectory(Resistor)
+add_subdirectory(VoltageSource)
+add_subdirectory(Inductor)
+add_subdirectory(LinearTransformer)
+add_subdirectory(InductionMotor)
+add_subdirectory(SynchronousMachine)
+add_subdirectory(DistributedGenerator)
+add_subdirectory(TransmissionLine)
+add_subdirectory(MicrogridLoad)
+add_subdirectory(MicrogridLine)
+add_subdirectory(MicrogridBusDQ)

ComponentLib/PowerElectronicsComponents/Capacitor/CMakeLists.txt

@@ -0,0 +1,8 @@
+
+
+
+gridkit_add_library(powerelec_capacitor
+  SOURCES
+    Capacitor.cpp
+  OUTPUT_NAME
+    gridkit_powerelec_capacitor)