DP_Ph1_Capacitor.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/Base/Base_Ph1_Capacitor.h>
#include <dpsim-models/MNASimPowerComp.h>
 
namespace CPS {
namespace DP {
namespace Ph1 {
class Capacitor : public MNASimPowerComp<Complex>,
                  public Base::Ph1::Capacitor,
                  public SharedFactory<Capacitor> {
protected:
  MatrixComp mEquivCurrent;
  MatrixComp mEquivCond;
  MatrixComp mPrevVoltCoeff;
 
public:
  Capacitor(String uid, String name,
            Logger::Level logLevel = Logger::Level::off);
  Capacitor(String name, Logger::Level logLevel = Logger::Level::off)
      : Capacitor(name, name, logLevel) {}
 
  SimPowerComp<Complex>::Ptr clone(String name) override;
 
  // #### General ####
  void initializeFromNodesAndTerminals(Real frequency) override;
  void initialize(Matrix frequencies) override;
 
  // #### MNA section ####
  Complex getMNAConductance() const;
  void mnaCompInitialize(Real omega, Real timeStep,
                         Attribute<Matrix>::Ptr leftVector) override;
  void mnaCompInitializeHarm(
      Real omega, Real timeStep,
      std::vector<Attribute<Matrix>::Ptr> leftVector) override;
  void mnaCompApplySystemMatrixStamp(SparseMatrixRow &systemMatrix) override;
  void mnaCompApplySystemMatrixStampHarm(SparseMatrixRow &systemMatrix,
                                         Int freqIdx) override;
  void mnaCompApplyRightSideVectorStamp(Matrix &rightVector) override;
  void mnaCompApplyRightSideVectorStampHarm(Matrix &rightVector) override;
  void mnaCompUpdateVoltage(const Matrix &leftVector) override;
  void mnaCompUpdateVoltageHarm(const Matrix &leftVector, Int freqIdx);
  void mnaCompApplyRightSideVectorStampHarm(Matrix &sourceVector,
                                            Int freqIdx) override;
  void mnaCompUpdateCurrent(const Matrix &leftVector) override;
  void mnaCompUpdateCurrentHarm();
  void mnaCompPreStep(Real time, Int timeStepCount) override;
  void mnaCompPostStep(Real time, Int timeStepCount,
                       Attribute<Matrix>::Ptr &leftVector) override;
  void mnaCompAddPreStepDependencies(
      AttributeBase::List &prevStepDependencies,
      AttributeBase::List &attributeDependencies,
      AttributeBase::List &modifiedAttributes) override;
  void
  mnaCompAddPostStepDependencies(AttributeBase::List &prevStepDependencies,
                                 AttributeBase::List &attributeDependencies,
                                 AttributeBase::List &modifiedAttributes,
                                 Attribute<Matrix>::Ptr &leftVector) override;
 
  class MnaPreStepHarm : public CPS::Task {
  public:
    MnaPreStepHarm(Capacitor &capacitor)
        : Task(**capacitor.mName + ".MnaPreStepHarm"), mCapacitor(capacitor) {
      // actually depends on C, but then we'd have to modify the system matrix anyway
      mModifiedAttributes.push_back(capacitor.attribute("right_vector"));
      mPrevStepDependencies.push_back(capacitor.attribute("i_intf"));
      mPrevStepDependencies.push_back(capacitor.attribute("v_intf"));
    }
    void execute(Real time, Int timeStepCount);
 
  private:
    Capacitor &mCapacitor;
  };
 
  class MnaPostStepHarm : public CPS::Task {
  public:
    MnaPostStepHarm(Capacitor &capacitor,
                    const std::vector<Attribute<Matrix>::Ptr> &leftVectors)
        : Task(**capacitor.mName + ".MnaPostStepHarm"), mCapacitor(capacitor),
          mLeftVectors(leftVectors) {
      for (UInt i = 0; i < mLeftVectors.size(); i++)
        mAttributeDependencies.push_back(mLeftVectors[i]);
      mModifiedAttributes.push_back(mCapacitor.attribute("v_intf"));
      mModifiedAttributes.push_back(mCapacitor.attribute("i_intf"));
    }
    void execute(Real time, Int timeStepCount);
 
  private:
    Capacitor &mCapacitor;
    std::vector<Attribute<Matrix>::Ptr> mLeftVectors;
  };
};
} // namespace Ph1
} // namespace DP
} // namespace CPS