SP_Ph1_NetworkInjection.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 <dpsim-models/CompositePowerComp.h>
#include <dpsim-models/SP/SP_Ph1_VoltageSource.h>
#include <dpsim-models/Solver/DAEInterface.h>
#include <dpsim-models/Solver/PFSolverInterfaceBus.h>
 
namespace CPS {
namespace SP {
namespace Ph1 {
class NetworkInjection : public CompositePowerComp<Complex>,
                         public SharedFactory<NetworkInjection>,
                         public PFSolverInterfaceBus,
                         public DAEInterface {
 
private:
  // ### Electrical Subcomponents ###
  std::shared_ptr<SP::Ph1::VoltageSource> mSubVoltageSource;
 
  // #### solver ####
  std::vector<const Matrix *> mRightVectorStamps;
 
  // #### Powerflow section ####
  Real mBaseVoltage;
 
public:
  const Attribute<Complex>::Ptr mVoltageRef;
  const Attribute<Real>::Ptr mSrcFreq;
 
  // #### Powerflow section ####
  const Attribute<Real>::Ptr mVoltageSetPoint;
  const Attribute<Real>::Ptr mVoltageSetPointPerUnit;
  const Attribute<Real>::Ptr mActivePowerInjection;
  const Attribute<Real>::Ptr mReactivePowerInjection;

  NetworkInjection(String uid, String name,
                   Logger::Level logLevel = Logger::Level::off);
  NetworkInjection(String name, Logger::Level logLevel = Logger::Level::off)
      : NetworkInjection(name, name, logLevel) {}

  SimPowerComp<Complex>::Ptr clone(String name) override;
 
  // #### General ####
  void initializeFromNodesAndTerminals(Real frequency) override;
 
  // #### Powerflow section ####
  void setParameters(Real vSetPointPerUnit);
  // Get base voltage
  Real getBaseVoltage() const;
  void setBaseVoltage(Real baseVoltage);
  void calculatePerUnitParameters(Real baseApparentPower, Real baseOmega);
  void modifyPowerFlowBusType(PowerflowBusType powerflowBusType) override;
  void updatePowerInjection(Complex powerInj);
 
  // #### MNA Section ####
  void setParameters(Complex voltageRef, Real srcFreq = 0.0);
  void setParameters(Complex initialPhasor, Real freqStart, Real rocof,
                     Real timeStart, Real duration, bool smoothRamp = true);
  void setParameters(Complex initialPhasor, Real modulationFrequency,
                     Real modulationAmplitude, Real baseFrequency = 0.0,
                     bool zigzag = false);

  void mnaParentApplyRightSideVectorStamp(Matrix &rightVector) override;
  void mnaCompUpdateCurrent(const Matrix &leftVector) override;
  void mnaCompUpdateVoltage(const Matrix &leftVector) override;
  void mnaParentPreStep(Real time, Int timeStepCount) override;
  void mnaParentPostStep(Real time, Int timeStepCount,
                         Attribute<Matrix>::Ptr &leftVector) override;
  void mnaParentAddPreStepDependencies(
      AttributeBase::List &prevStepDependencies,
      AttributeBase::List &attributeDependencies,
      AttributeBase::List &modifiedAttributes) override;
  void
  mnaParentAddPostStepDependencies(AttributeBase::List &prevStepDependencies,
                                   AttributeBase::List &attributeDependencies,
                                   AttributeBase::List &modifiedAttributes,
                                   Attribute<Matrix>::Ptr &leftVector) override;
 
  // #### DAE Section ####
  void daeResidual(double ttime, const double state[], const double dstate_dt[],
                   double resid[], std::vector<int> &off) override;
  Complex daeInitialize() override;
};
} // namespace Ph1
} // namespace SP
} // namespace CPS