SP_Ph1_RXLine.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/Solver/PFSolverInterfaceBranch.h"
#include <dpsim-models/Base/Base_Ph1_PiLine.h>
#include <dpsim-models/CompositePowerComp.h>
#include <dpsim-models/SP/SP_Ph1_Inductor.h>
#include <dpsim-models/SP/SP_Ph1_Resistor.h>
 
namespace CPS {
namespace SP {
namespace Ph1 {
 
class RXLine : public CompositePowerComp<Complex>,
               public SharedFactory<RXLine>,
               public PFSolverInterfaceBranch,
               public Base::Ph1::PiLine {
private:
  Real mBaseVoltage;
 
protected:
  Real mCapacitance;
  Real mBaseCurrent;
  Real mBaseApparentPower;
  Real mBaseOmega;
  Real mBaseImpedance;
  Real mBaseAdmittance;
  Real mBaseInductance;
  Real mBaseCapacitance;

  Real mSeriesResPerUnit;
  Real mParallelCapPerUnit;
  Real mSeriesIndPerUnit;
 
  // #### Admittance matrix stamp ####
  MatrixComp mY_element;

  std::shared_ptr<Inductor> mSubInductor;
  std::shared_ptr<Resistor> mSubResistor;
  std::shared_ptr<Resistor> mInitialResistor;
 
public:
  const Attribute<Real>::Ptr mInductance;

  const Attribute<Real>::Ptr mActivePowerInjection;
  const Attribute<Real>::Ptr mReactivePowerInjection;
 
  // #### Power flow results ####
  const Attribute<MatrixComp>::Ptr mCurrent;

  const Attribute<Matrix>::Ptr mActivePowerBranch;

  const Attribute<Matrix>::Ptr mReactivePowerBranch;
 
  // #### constructors ####
  // power flow

  RXLine(String uid, String name, Real baseVoltage, Real resistance,
         Real inductance, Logger::Level logLevel = Logger::Level::off);
  // MNA
  RXLine(String uid, String name, Logger::Level logLevel = Logger::Level::off);
  RXLine(String name, Logger::Level logLevel = Logger::Level::off)
      : RXLine(name, name, logLevel) {}
 
  // #### General ####
  void setPerUnitSystem(Real baseApparentPower, Real baseOmega);
  void transformParametersToPerUnitSystem();
 
  // #### Powerflow section ####
  Real getBaseVoltage() const;
  void pfApplyAdmittanceMatrixStamp(SparseMatrixCompRow &Y) override;
  void updateBranchFlow(VectorComp &current, VectorComp &powerflow);
  void storeNodalInjection(Complex powerInjection);
 
  // #### Getter ####
  MatrixComp Y_element();
 
  // #### MNA Section ####
 
  SimPowerComp<Complex>::Ptr clone(String name) override;
 
  // #### General ####
  void initializeFromNodesAndTerminals(Real frequency) override;
 
  // #### MNA section ####
  void mnaCompUpdateVoltage(const Matrix &leftVector) override;
  void mnaCompUpdateCurrent(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;
};
} // namespace Ph1
} // namespace SP
} // namespace CPS