DP_Ph1_Switch.h

Go to the documentation of this file.

/* 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_Switch.h>
#include <dpsim-models/Definitions.h>
#include <dpsim-models/Logger.h>
#include <dpsim-models/MNASimPowerComp.h>
#include <dpsim-models/Solver/MNASwitchInterface.h>
#include <dpsim-models/Solver/MNAVariableCompInterface.h>
 
namespace CPS {
namespace DP {
namespace Ph1 {
class Switch : public MNASimPowerComp<Complex>,
               public Base::Ph1::Switch,
               public SharedFactory<Switch>,
               public MNASwitchInterface,
               public MNAVariableCompInterface {
public:
  enum class SwitchingMode {
    Ideal = 0,
    CurrentZero = 1,
    ExponentialZCSEmulation = 2,
  };
 
  const Attribute<Bool>::Ptr mOpeningRequested;
 
  const Attribute<Bool>::Ptr mPoleClosed;

  const Attribute<Real>::Ptr mInstantCurrent;
 
  const Attribute<Real>::Ptr mZeroCrossingTime;

  const Attribute<Bool>::Ptr mExponentialTransitionActive;
  const Attribute<Bool>::Ptr mExponentialTransitionClosing;
  const Attribute<Real>::Ptr mExponentialProgress;
  const Attribute<Real>::Ptr mExponentialTransitionStartTime;
  const Attribute<Real>::Ptr mExponentialTransitionEndTime;

  const Attribute<Real>::Ptr mEffectiveResistance;

  Switch(String uid, String name, Logger::Level loglevel = Logger::Level::off);
  Switch(String name, Logger::Level logLevel = Logger::Level::off)
      : Switch(name, name, logLevel) {}
 
  SimPowerComp<Complex>::Ptr clone(String name) override;
 
  void setSwitchingMode(SwitchingMode mode);
  SwitchingMode switchingMode() const { return mSwitchingMode; }
 
  void setZeroCrossingTolerance(Real tolerance);
 
  void setExponentialSwitchingTime(Real switchingTime);
  Real exponentialSwitchingTime() const { return mExponentialSwitchingTime; }
 
  void openSwitch() override;
  void closeSwitch() override;
 
  // #### General ####
  void initializeFromNodesAndTerminals(Real frequency) override;
 
  // #### General 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 mnaCompUpdateCurrent(const Matrix &leftVector) override;
  void mnaCompPostStep(Real time, Int timeStepCount,
                       Attribute<Matrix>::Ptr &leftVector) override;
  void
  mnaCompAddPostStepDependencies(AttributeBase::List &prevStepDependencies,
                                 AttributeBase::List &attributeDependencies,
                                 AttributeBase::List &modifiedAttributes,
                                 Attribute<Matrix>::Ptr &leftVector) override;
 
  // #### MNA section for switch ####
  Bool mnaIsClosed() override;
  void mnaCompApplySwitchSystemMatrixStamp(Bool closed,
                                           SparseMatrixRow &systemMatrix,
                                           Int freqIdx) override;
 
  Bool supportsPrecomputedSystemMatrices() const override;
 
  // #### MNA section for variable component ####
  Bool hasParameterChanged() override;
 
private:
  SwitchingMode mSwitchingMode = SwitchingMode::Ideal;
 
  Real mZeroCrossingTolerance = 1e-6;
  Real mExponentialSwitchingTime = 0.01;
  Real mTimeStep = 0.0;

  Real mShiftOmega = 0.0;
 
  Bool mPoleClosedPrev = false;
  Real mEffectiveResistancePrev = 0.0;
 
  Real mPreviousInstantaneousCurrent = 0.0;
  Real mPreviousCurrentTime = 0.0;
  Bool mPreviousCurrentValid = false;
 
  Bool mResetZeroCrossingHistory = false;
  Bool mExponentialTransitionStarted = false;
 
  void setPoleClosed(Bool closed);
 
  void resetZeroCrossingTime();
 
  void synchronizeEffectiveResistance(Bool closed);
 
  void resetExponentialTransition();
  void validateExponentialResistanceParameters() const;
  Real exponentialResistance(Real alpha) const;
  void updateExponentialTransition(Real time);
 
  Complex currentConductance() const;
 
  Real reconstructInstantaneousCurrent(Real time) const;
 
  Bool currentCrossedZero(Real previousCurrent, Real current) const;
 
  Real interpolateZeroCrossingTime(Real previousCurrent, Real current,
                                   Real previousTime, Real currentTime) const;
 
  void updateZeroCrossingState(Real time);
};
} // namespace Ph1
} // namespace DP
} // namespace CPS