SP_Ph1_PiLine.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_PiLine.h>
#include <dpsim-models/CompositePowerComp.h>
#include <dpsim-models/SP/SP_Ph1_Capacitor.h>
#include <dpsim-models/SP/SP_Ph1_Inductor.h>
#include <dpsim-models/SP/SP_Ph1_Resistor.h>
#include <dpsim-models/Solver/MNATearInterface.h>
#include <dpsim-models/Solver/PFSolverInterfaceBranch.h>
 
namespace CPS {
namespace SP {
namespace Ph1 {
class PiLine : public CompositePowerComp<Complex>,
               public Base::Ph1::PiLine,
               public MNATearInterface,
               public SharedFactory<PiLine>,
               public PFSolverInterfaceBranch {
private:
  Real mBaseVoltage;
 
public:
  // #### Power flow results ####
  const Attribute<MatrixComp>::Ptr mCurrent;

  const Attribute<Matrix>::Ptr mActivePowerBranch;

  const Attribute<Matrix>::Ptr mReactivePowerBranch;
 
protected:
  Real mBaseCurrent;
  Real mBaseApparentPower;
  Real mBaseOmega;
  Real mBaseImpedance;
  Real mBaseAdmittance;
  Real mBaseInductance;
  Real mBaseCapacitance;

  Real mSeriesResPerUnit;
  Real mParallelCapPerUnit;
  Real mSeriesIndPerUnit;
  Real mParallelCondPerUnit;
 
  // #### Admittance matrix stamp ####
  MatrixComp mY_element;
  std::shared_ptr<Inductor> mSubSeriesInductor;
  std::shared_ptr<Resistor> mSubSeriesResistor;
  std::shared_ptr<Resistor> mSubParallelResistor0;
  // Parallel Capacitor submodel at Terminal 0
  std::shared_ptr<Capacitor> mSubParallelCapacitor0;
  std::shared_ptr<Resistor> mSubParallelResistor1;
  std::shared_ptr<Capacitor> mSubParallelCapacitor1;
  std::vector<const Matrix *> mRightVectorStamps;
 
public:
  const Attribute<Real>::Ptr mActivePowerInjection;
  const Attribute<Real>::Ptr mReactivePowerInjection;
 
  // #### General ####
  PiLine(String uid, String name, Logger::Level logLevel = Logger::Level::off);
  PiLine(String name, Logger::Level logLevel = Logger::Level::off)
      : PiLine(name, name, logLevel) {}

  SimPowerComp<Complex>::Ptr clone(String copySuffix) override;
  void setParameters(Real resistance, Real inductance, Real capacitance = -1,
                     Real conductance = -1);
  void createSubComponents() override;
  void initializeParentFromNodesAndTerminals(Real frequency) override;
 
  // #### Powerflow section ####
  Real getBaseVoltage() const;
  void setBaseVoltage(Real baseVoltage);
  void calculatePerUnitParameters(Real baseApparentPower, Real baseOmega);
  void pfApplyAdmittanceMatrixStamp(SparseMatrixCompRow &Y) override;
  void updateBranchFlow(VectorComp &current, VectorComp &powerflow);
  void storeNodalInjection(Complex powerInjection);
 
  // #### Getter ####
  MatrixComp Y_element();
 
  // #### MNA section ####
  void mnaCompUpdateCurrent(const Matrix &leftVector) override;
  void mnaCompUpdateVoltage(const Matrix &leftVector) override;
  void mnaParentPostStep(Real time, Int timeStepCount,
                         Attribute<Matrix>::Ptr &leftVector) override;
  void
  mnaParentAddPostStepDependencies(AttributeBase::List &prevStepDependencies,
                                   AttributeBase::List &attributeDependencies,
                                   AttributeBase::List &modifiedAttributes,
                                   Attribute<Matrix>::Ptr &leftVector) override;
 
  MNAInterface::List mnaTearGroundComponents() override;
  void mnaTearInitialize(Real omega, Real timeStep) override;
  void mnaTearApplyMatrixStamp(SparseMatrixRow &tearMatrix) override;
  void mnaTearApplyVoltageStamp(Matrix &voltageVector) override;
  void mnaTearPostStep(Complex voltage, Complex current) override;
};
} // namespace Ph1
} // namespace SP
} // namespace CPS