EMT_Ph3_PQLoad.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/MNASimPowerComp.h>
#include <dpsim-models/Solver/MNASyncGenInterface.h>
 
namespace CPS {
namespace EMT {
namespace Ph3 {

class PQLoad : public MNASimPowerComp<Real>,
               public MNASyncGenInterface,
               public SharedFactory<PQLoad> {
protected:
  static Matrix quadratureVoltage(const Matrix &voltage);
  Matrix calculateCurrent(const Matrix &voltage) const;
  void updateMeasuredPower();
 
public:
  const Attribute<Real>::Ptr mActivePower;
  const Attribute<Real>::Ptr mReactivePower;
  const Attribute<Real>::Ptr mNomVoltage;
  const Attribute<Real>::Ptr mMinimumVoltagePerUnit;
  const Attribute<Real>::Ptr mMeasuredActivePower;
  const Attribute<Real>::Ptr mMeasuredReactivePower;

  PQLoad(String uid, String name, Logger::Level logLevel = Logger::Level::off);
  PQLoad(String name, Logger::Level logLevel = Logger::Level::off);
  PQLoad(String name, Real activePower, Real reactivePower, Real nominalVoltage,
         Logger::Level logLevel = Logger::Level::off);
  PQLoad(String uid, String name, Real activePower, Real reactivePower,
         Real nominalVoltage, Logger::Level logLevel = Logger::Level::off);

  void setParameters(Real activePower, Real reactivePower, Real nominalVoltage,
                     Real minimumVoltagePerUnit = 0.1);
  SimPowerComp<Real>::Ptr clone(String name) override;
 
  String description() override {
    return fmt::format("Active: {}MW, Reactive: {}MVAr, Voltage: {}kV",
                       **mActivePower / 1e6, **mReactivePower / 1e6,
                       **mNomVoltage / 1e3);
  }

  void initializeFromNodesAndTerminals(Real frequency) override;
 
  // #### MNA section ####
  void mnaCompInitialize(Real omega, Real timeStep,
                         Attribute<Matrix>::Ptr leftVector) override;
  void mnaCompApplySystemMatrixStamp(SparseMatrixRow &systemMatrix) override;
  void mnaCompApplyRightSideVectorStamp(Matrix &rightVector) override;
  void mnaCompUpdateVoltage(const Matrix &leftVector) override;
  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;
 
  // #### Same-time-step fixed-point correction ####
  void correctorStep() override;
  void updateVoltage(const Matrix &leftVector) override;
  bool requiresIteration() override;
};
 
} // namespace Ph3
} // namespace EMT
} // namespace CPS