PFSolverPowerPolar.h

/* Copyright 2017-2021 Institute for Automation of Complex Power Systems,
 *                     EONERC, RWTH Aachen University
 *
 * This Source Code Form is subject to the terms of the Mozilla Public
 * License, v. 2.0. If a copy of the MPL was not distributed with this
 * file, You can obtain one at https://mozilla.org/MPL/2.0/.
 *********************************************************************************/
 
#pragma once
 
#include <map>
 
#include <dpsim/PFSolver.h>
 
namespace DPsim {
class PFSolverPowerPolar : public PFSolver {
protected:
  CPS::Vector sol_P;
  CPS::Vector sol_Q;
  CPS::Vector sol_V;
  CPS::Vector sol_D;
  CPS::VectorComp sol_V_complex;
  CPS::VectorComp sol_S_complex;
 
  CPS::Vector Pesp;
  CPS::Vector Qesp;
 
  // Cache for last converged solution
  CPS::Vector mLastConvergedV;
  CPS::Vector mLastConvergedD;
  bool mHasLastConvergedSolution = false;
 
  // Core methods
  void generateInitialSolution(Real time,
                               bool keep_last_solution = false) override;
  void calculateJacobian() override;
  void updateSolution() override;
  void setSolution() override;
  void calculateMismatch() override;
 
  // Helper methods
  void resize_sol(CPS::Int n);
  void resize_complex_sol(CPS::Int n);
  CPS::Real sol_Vr(CPS::UInt k);
  CPS::Real sol_Vi(CPS::UInt k);
  CPS::Complex sol_Vcx(CPS::UInt k);
  CPS::Real P(CPS::UInt k);
  CPS::Real Q(CPS::UInt k);
  void calculatePAndQAtSlackBus();
  void calculateQAtPVBuses();
  CPS::Real generatorReactivePowerPerUnit(CPS::TopologicalNode::Ptr node);
  CPS::Real loadReactivePowerPerUnit(CPS::TopologicalNode::Ptr node);
  CPS::Bool enforceReactiveLimits() override;
  void clearReactiveLimitState() override;
  std::map<CPS::TopologicalNode::Ptr, bool> mQLimitConvertedAtMax;
  std::map<CPS::TopologicalNode::Ptr, CPS::UInt> mQLimitSwitchCount;
  void calculatePAndQInjectionPQBuses();
  void calculateBranchFlow();
  void calculateNodalInjection();
 
public:
  PFSolverPowerPolar(CPS::String name, const CPS::SystemTopology &system,
                     CPS::Real timeStep, CPS::Logger::Level logLevel);
  virtual ~PFSolverPowerPolar(){};
};
} // namespace DPsim