EMT_Ph3_HalfDecouplingLine.h

Go to the documentation of this file.

// SPDX-FileCopyrightText: 2026 Institute for Automation of Complex Power Systems, EONERC, RWTH Aachen University
// SPDX-License-Identifier: MPL-2.0
#pragma once
#include <dpsim-models/CompositePowerComp.h>
#include <dpsim-models/Definitions.h>
#include <dpsim-models/EMT/EMT_Ph3_ControlledCurrentSource.h>
#include <dpsim-models/EMT/EMT_Ph3_Resistor.h>
#include <dpsim-models/Solver/MNAInterface.h>
 
namespace CPS {
namespace EMT {
namespace Ph3 {
class HalfDecouplingLine : public CompositePowerComp<Real>,
                           public SharedFactory<HalfDecouplingLine> {
protected:
  Real mDelay;
  Matrix mResistance = Matrix::Zero(3, 3);
  Matrix mInductance = Matrix::Zero(3, 3);
  Matrix mCapacitance = Matrix::Zero(3, 3);
  Matrix mSurgeImpedance = Matrix::Zero(3, 3);
 
  // ### Electrical Subcomponents ###
  std::shared_ptr<EMT::Ph3::ControlledCurrentSource> mSubCtrledCurrentSource;
  std::shared_ptr<EMT::Ph3::Resistor> mSubRes;

  Matrix mVoltBuf, mCurBuf;
  UInt mBufIdx = 0;
  UInt mBufSize;
  Real mAlpha;

  Matrix interpolate(Matrix &data);
 
public:
  typedef std::shared_ptr<HalfDecouplingLine> Ptr;

  const Attribute<Matrix>::Ptr mSrcCtrledCurrent;
  const Attribute<Matrix>::Ptr mSrcRes;
  const Attribute<Matrix>::Ptr mReceivingVolt;
  const Attribute<Matrix>::Ptr mReceivingCur;
  const Attribute<Matrix>::Ptr mSendingVolt;
  const Attribute<Matrix>::Ptr mSendingCur;

  HalfDecouplingLine(String name, Logger::Level logLevel = Logger::Level::off)
      : HalfDecouplingLine(name, name, logLevel) {}

  HalfDecouplingLine(String uid, String name,
                     Logger::Level logLevel = Logger::Level::off);
 
  // #### General ####
  void setParameters(Matrix resistance, Matrix inductance, Matrix capacitance);
  void setCouplingSource(Attribute<Matrix>::Ptr receivingVolt,
                         Attribute<Matrix>::Ptr receivingCur);
  void createSubComponents() override;
  void initializeParentFromNodesAndTerminals(Real frequency) override;

  void step(Real time, Int timeStepCount);
  void postStep();
 
  // #### MNA section ####
  void mnaParentInitialize(Real omega, Real timeStep,
                           Attribute<Matrix>::Ptr leftVector) 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;
};
} // namespace Ph3
} // namespace EMT
} // namespace CPS