Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Jump to bottom

Fix power flow initialisation of synchronous generators #238

Merged
merged 10 commits into from
Nov 30, 2023
Merged

Fix power flow initialisation of synchronous generators #238

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
10 commits
Select commit Hold shift + click to select a range
f18262f
fix initialization of SG power after powerflow
martinmoraga Jun 20, 2023
cef7439
adapt function initWithPowerflow to initialize SG power from power flow
martinmoraga Jun 20, 2023
c0fb8b7
remove unused function setSGPower
martinmoraga Jul 27, 2023
9201481
add function setPower<Complex> to SimNode
martinmoraga Aug 2, 2023
f2eec8e
correct power calculation of PQ buses after powerflow
martinmoraga Aug 2, 2023
0331903
add attributes mConductancePerUnit and mSusceptancePerUnit to SP_Shunt
martinmoraga Aug 3, 2023
981abc6
improve some comments
dinkelbachjan Nov 28, 2023
234b3f9
enforce setting of domain argument for initWithPowerflow
dinkelbachjan Nov 29, 2023
7fb0981
fix pybind double definition of components_at_node
dinkelbachjan Nov 29, 2023
9405f52
add comments on power solution setters in powerflow solver
dinkelbachjan Nov 29, 2023
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
9 changes: 4 additions & 5 deletions dpsim-models/include/dpsim-models/SP/SP_Ph1_Shunt.h
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -20,15 +20,14 @@ namespace SP {namespace Ph1 {
const Attribute<Real>::Ptr mConductance;
/// Susceptance [S]
const Attribute<Real>::Ptr mSusceptance;
/// Conductance [pu]
const Attribute<Real>::Ptr mConductancePerUnit;
/// Susceptance [pu]
const Attribute<Real>::Ptr mSusceptancePerUnit;
private:
/// Base voltage [V]
Real mBaseVoltage;

/// Conductance [pu]
Real mConductancePerUnit;
/// Susceptance [pu]
Real mSusceptancePerUnit;

public:
/// Defines UID, name, component parameters and logging level
Shunt(String uid, String name, Logger::Level logLevel = Logger::Level::off);
Expand Down
3 changes: 3 additions & 0 deletions dpsim-models/include/dpsim-models/SimNode.h
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -127,4 +127,7 @@ namespace CPS {

template<>
void SimNode<Complex>::setVoltage(Complex newVoltage);

template<>
void SimNode<Complex>::setPower(Complex newPower);
}
5 changes: 2 additions & 3 deletions dpsim-models/include/dpsim-models/SystemTopology.h
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -14,7 +14,6 @@
#include <dpsim-models/TopologicalPowerComp.h>
#include <dpsim-models/SimPowerComp.h>
#include <dpsim-models/SimNode.h>
#include <dpsim-models/DP/DP_Ph1_Resistor.h>

#ifdef WITH_GRAPHVIZ
#include <dpsim-models/Graph.h>
Expand Down Expand Up @@ -107,8 +106,8 @@ namespace CPS {
/// Add multiple components
void addComponents(const IdentifiedObject::List& components);

/// Initialize nodes from PowerFlow
void initWithPowerflow(const SystemTopology& systemPF);
/// Initialize nodes and SG power from PowerFlow
void initWithPowerflow(const SystemTopology& systemPF, CPS::Domain domain);

/// Adds component and initializes frequencies
void addTearComponent(IdentifiedObject::Ptr component);
Expand Down
12 changes: 7 additions & 5 deletions dpsim-models/src/SP/SP_Ph1_Shunt.cpp
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -13,7 +13,9 @@ using namespace CPS;
SP::Ph1::Shunt::Shunt(String uid, String name, Logger::Level logLevel)
: SimPowerComp<Complex>(uid, name, logLevel),
mConductance(mAttributes->create<Real>("G")),
mSusceptance(mAttributes->create<Real>("B")) {
mSusceptance(mAttributes->create<Real>("B")),
mConductancePerUnit(mAttributes->create<Real>("Gpu")),
mSusceptancePerUnit(mAttributes->create<Real>("Bpu")) {

SPDLOG_LOGGER_INFO(mSLog, "Create {} of type {}", this->type(), name);
setTerminalNumber(1);
Expand Down Expand Up @@ -42,15 +44,15 @@ void SP::Ph1::Shunt::calculatePerUnitParameters(Real baseApparentPower, Real bas
auto baseAdmittance = 1.0 / baseImpedance;
SPDLOG_LOGGER_INFO(mSLog, "Base Voltage={} [V] Base Admittance={} [S]", mBaseVoltage, baseAdmittance);

mConductancePerUnit = **mConductance / baseAdmittance;
mSusceptancePerUnit = **mSusceptance / baseAdmittance;
SPDLOG_LOGGER_INFO(mSLog, "Susceptance={} [pu] Conductance={} [pu]", mSusceptancePerUnit, mConductancePerUnit);
**mConductancePerUnit = **mConductance / baseAdmittance;
**mSusceptancePerUnit = **mSusceptance / baseAdmittance;
SPDLOG_LOGGER_INFO(mSLog, "Susceptance={} [pu] Conductance={} [pu]", **mSusceptancePerUnit, **mConductancePerUnit);
};


void SP::Ph1::Shunt::pfApplyAdmittanceMatrixStamp(SparseMatrixCompRow & Y) {
int bus1 = this->matrixNodeIndex(0);
Complex Y_element = Complex(mConductancePerUnit, mSusceptancePerUnit);
Complex Y_element = Complex(**mConductancePerUnit, **mSusceptancePerUnit);

if (std::isinf(Y_element.real()) || std::isinf(Y_element.imag())) {
std::cout << "Y:" << Y_element << std::endl;
Expand Down
20 changes: 19 additions & 1 deletion dpsim-models/src/SystemTopology.cpp
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -12,6 +12,7 @@
#include <unordered_map>

#include <dpsim-models/SystemTopology.h>
#include <dpsim-models/SP/SP_Ph1_SynchronGenerator.h>

using namespace CPS;

Expand Down Expand Up @@ -72,7 +73,8 @@ void SystemTopology::addComponents(const IdentifiedObject::List& components) {
addComponent(comp);
}

void SystemTopology::initWithPowerflow(const SystemTopology& systemPF) {
void SystemTopology::initWithPowerflow(const SystemTopology& systemPF, CPS::Domain domain) {

for (auto nodePF : systemPF.mNodes) {
if (auto node = this->node<TopologicalNode>(nodePF->name())) {
//SPDLOG_LOGGER_INFO(mSLog, "Updating initial voltage of {} according to powerflow", node->name());
Expand All @@ -81,6 +83,22 @@ void SystemTopology::initWithPowerflow(const SystemTopology& systemPF) {
//SPDLOG_LOGGER_INFO(mSLog, "Updated initial voltage: {}", node->initialSingleVoltage());
}
}

// set initial power of SG
for (auto compPF : systemPF.mComponents) {
if (auto genPF = std::dynamic_pointer_cast<CPS::SP::Ph1::SynchronGenerator>(compPF)) {
if (domain==CPS::Domain::DP || domain==CPS::Domain::SP) {
auto comp = this->component<SimPowerComp<Complex>>(compPF->name());
auto terminal = comp->terminals()[0];
terminal->setPower(-genPF->getApparentPower());
} else if (domain==CPS::Domain::EMT) {
auto comp = this->component<SimPowerComp<Real>>(compPF->name());
auto terminal = comp->terminals()[0];
terminal->setPower(-genPF->getApparentPower());
}
//SPDLOG_LOGGER_INFO(mSLog, "Updated initial power of gen {}: {}", compPF->name(), genPF->getApparentPower());
}
}
}

void SystemTopology::addTearComponent(IdentifiedObject::Ptr component) {
Expand Down
2 changes: 1 addition & 1 deletion dpsim/examples/cxx/CIM/DP_CIGRE_MV_withDG.cpp
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -79,7 +79,7 @@ int main(int argc, char** argv){
CIM::Reader reader2(simName, Logger::Level::debug, Logger::Level::debug);
SystemTopology systemDP = reader2.loadCIM(scenario.systemFrequency, filenames, CPS::Domain::DP);
Examples::Grids::CIGREMV::addInvertersToCIGREMV(systemDP, scenario, Domain::DP);
systemDP.initWithPowerflow(systemPF);
systemDP.initWithPowerflow(systemPF, CPS::Domain::DP);

auto logger = DPsim::DataLogger::make(simName);

Expand Down
4 changes: 1 addition & 3 deletions dpsim/examples/cxx/CIM/DP_CIGRE_MV_withDG_withLoadStep.cpp
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -75,15 +75,13 @@ int main(int argc, char** argv){
CIM::Reader reader2(simName, Logger::Level::info, Logger::Level::debug);
SystemTopology systemDP = reader2.loadCIM(scenario.systemFrequency, filenames, CPS::Domain::DP);
Examples::Grids::CIGREMV::addInvertersToCIGREMV(systemDP, scenario, Domain::DP);
systemDP.initWithPowerflow(systemPF);
systemDP.initWithPowerflow(systemPF, CPS::Domain::DP);

auto logger = DPsim::DataLogger::make(simName);

// log node voltages
for (auto node : systemDP.mNodes)
{
logger->logAttribute(node->name() + ".V", node->attribute("v"));
}

// log line currents
for (auto comp : systemDP.mComponents) {
Expand Down
2 changes: 1 addition & 1 deletion dpsim/examples/cxx/CIM/DP_CIGRE_MV_withoutDG.cpp
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -71,7 +71,7 @@ int main(int argc, char** argv){
Logger::setLogDir("logs/" + simName);
CIM::Reader reader2(simName, Logger::Level::debug, Logger::Level::debug);
SystemTopology systemDP = reader2.loadCIM(systemFrequency, filenames, CPS::Domain::DP);
systemDP.initWithPowerflow(systemPF);
systemDP.initWithPowerflow(systemPF, CPS::Domain::DP);

auto logger = DPsim::DataLogger::make(simName);

Expand Down
4 changes: 2 additions & 2 deletions dpsim/examples/cxx/CIM/DP_WSCC-9bus_IdealVS.cpp
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -47,7 +47,7 @@ int main(int argc, char *argv[]) {
CPS::CIM::Reader reader(simNamePF, Logger::Level::debug, Logger::Level::debug);
SystemTopology systemPF = reader.loadCIM(60, filenames, Domain::SP, PhaseType::Single, CPS::GeneratorType::PVNode);
systemPF.component<CPS::SP::Ph1::SynchronGenerator>("GEN1")->modifyPowerFlowBusType(CPS::PowerflowBusType::VD);

// define logging
auto loggerPF = DPsim::DataLogger::make(simNamePF);
for (auto node : systemPF.mNodes)
Expand All @@ -72,7 +72,7 @@ int main(int argc, char *argv[]) {

CPS::CIM::Reader reader2(simName, Logger::Level::debug, Logger::Level::debug);
SystemTopology sys = reader2.loadCIM(60, filenames, Domain::DP, PhaseType::Single, CPS::GeneratorType::IdealVoltageSource);
sys.initWithPowerflow(systemPF);
sys.initWithPowerflow(systemPF, CPS::Domain::DP);

// Logging
auto logger = DataLogger::make(simName);
Expand Down
7 changes: 2 additions & 5 deletions dpsim/examples/cxx/CIM/DP_WSCC9bus_SGReducedOrderIter.cpp
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -139,13 +139,10 @@ int main(int argc, char *argv[]) {
sys.addComponent(faultDP);
}

sys.initWithPowerflow(systemPF);
sys.initWithPowerflow(systemPF, Domain::DP);
for (auto comp : sys.mComponents) {
if (auto genReducedOrder = std::dynamic_pointer_cast<CPS::Base::ReducedOrderSynchronGenerator<Complex>>(comp)) {
auto genPF = systemPF.component<CPS::SP::Ph1::SynchronGenerator>(comp->name());
genReducedOrder->terminal(0)->setPower(-genPF->getApparentPower());
if (auto genReducedOrder = std::dynamic_pointer_cast<CPS::Base::ReducedOrderSynchronGenerator<Complex>>(comp))
genReducedOrder->scaleInertiaConstant(inertiaScalingFactor);
}
}

// Logging
Expand Down
4 changes: 1 addition & 3 deletions dpsim/examples/cxx/CIM/DP_WSCC9bus_SGReducedOrderVBR.cpp
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -135,11 +135,9 @@ int main(int argc, char *argv[]) {
sys.addComponent(faultDP);
}

sys.initWithPowerflow(systemPF);
sys.initWithPowerflow(systemPF, Domain::DP);
for (auto comp : sys.mComponents) {
if (auto genReducedOrder = std::dynamic_pointer_cast<CPS::Base::ReducedOrderSynchronGenerator<Complex>>(comp)) {
auto genPF = systemPF.component<CPS::SP::Ph1::SynchronGenerator>(comp->name());
genReducedOrder->terminal(0)->setPower(-genPF->getApparentPower());
genReducedOrder->scaleInertiaConstant(inertiaScalingFactor);
genReducedOrder->setModelAsNortonSource(false);
}
Expand Down
2 changes: 1 addition & 1 deletion dpsim/examples/cxx/CIM/EMT_CIGRE_MV_withDG.cpp
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -76,7 +76,7 @@ int main(int argc, char** argv){
CIM::Reader reader2(simName, Logger::Level::debug, Logger::Level::debug);
SystemTopology systemEMT = reader2.loadCIM(scenario.systemFrequency, filenames, CPS::Domain::EMT, PhaseType::ABC);
Examples::Grids::CIGREMV::addInvertersToCIGREMV(systemEMT, scenario, Domain::EMT);
systemEMT.initWithPowerflow(systemPF);
systemEMT.initWithPowerflow(systemPF, CPS::Domain::EMT);

auto logger = DPsim::DataLogger::make(simName);

Expand Down
2 changes: 1 addition & 1 deletion dpsim/examples/cxx/CIM/EMT_CIGRE_MV_withDG_withLoadStep.cpp
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -76,7 +76,7 @@ int main(int argc, char** argv){
CIM::Reader reader2(simName, Logger::Level::debug, Logger::Level::debug);
SystemTopology systemEMT = reader2.loadCIM(scenario.systemFrequency, filenames, CPS::Domain::EMT, PhaseType::ABC);
Examples::Grids::CIGREMV::addInvertersToCIGREMV(systemEMT, scenario, Domain::EMT);
systemEMT.initWithPowerflow(systemPF);
systemEMT.initWithPowerflow(systemPF, CPS::Domain::EMT);

auto logger = DPsim::DataLogger::make(simName);

Expand Down
2 changes: 1 addition & 1 deletion dpsim/examples/cxx/CIM/EMT_CIGRE_MV_withoutDG.cpp
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -69,7 +69,7 @@ int main(int argc, char** argv){
Logger::setLogDir("logs/" + simName);
CIM::Reader reader2(simName, Logger::Level::debug, Logger::Level::debug);
SystemTopology systemEMT = reader2.loadCIM(systemFrequency, filenames, CPS::Domain::EMT, PhaseType::ABC);
systemEMT.initWithPowerflow(systemPF);
systemEMT.initWithPowerflow(systemPF, CPS::Domain::EMT);

auto logger = DPsim::DataLogger::make(simName);

Expand Down
8 changes: 1 addition & 7 deletions dpsim/examples/cxx/CIM/EMT_WSCC-9bus_FullOrderSG.cpp
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -73,13 +73,7 @@ int main(int argc, char *argv[]) {
CPS::CIM::Reader reader2(simName, Logger::Level::debug, Logger::Level::debug);
SystemTopology sys = reader2.loadCIM(60, filenames, Domain::EMT, PhaseType::ABC, CPS::GeneratorType::FullOrder);

sys.initWithPowerflow(systemPF);
for (auto comp : sys.mComponents) {
if (auto genEMT = std::dynamic_pointer_cast<CPS::EMT::Ph3::SynchronGeneratorDQTrapez>(comp)) {
auto genPF = systemPF.component<CPS::SP::Ph1::SynchronGenerator>(comp->name());
genEMT->terminal(0)->setPower(-genPF->getApparentPower());
}
}
sys.initWithPowerflow(systemPF, CPS::Domain::EMT);

// Logging
auto logger = DataLogger::make(simName);
Expand Down
8 changes: 1 addition & 7 deletions dpsim/examples/cxx/CIM/EMT_WSCC-9bus_IdealCS.cpp
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -73,13 +73,7 @@ int main(int argc, char *argv[]) {
CPS::CIM::Reader reader2(simName, Logger::Level::debug, Logger::Level::debug);
SystemTopology sys = reader2.loadCIM(60, filenames, Domain::EMT, PhaseType::ABC, CPS::GeneratorType::IdealCurrentSource);

sys.initWithPowerflow(systemPF);
for (auto comp : sys.mComponents) {
if (auto genEMT = std::dynamic_pointer_cast<CPS::EMT::Ph3::SynchronGeneratorIdeal>(comp)) {
auto genPF = systemPF.component<CPS::SP::Ph1::SynchronGenerator>(comp->name());
genEMT->terminal(0)->setPower(-genPF->getApparentPower());
}
}
sys.initWithPowerflow(systemPF, CPS::Domain::EMT);

// Logging
auto logger = DataLogger::make(simName);
Expand Down
2 changes: 1 addition & 1 deletion dpsim/examples/cxx/CIM/EMT_WSCC-9bus_IdealVS.cpp
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -72,7 +72,7 @@ int main(int argc, char *argv[]) {

CPS::CIM::Reader reader2(simName, Logger::Level::debug, Logger::Level::debug);
SystemTopology sys = reader2.loadCIM(60, filenames, Domain::EMT, PhaseType::ABC, CPS::GeneratorType::IdealVoltageSource);
sys.initWithPowerflow(systemPF);
sys.initWithPowerflow(systemPF, CPS::Domain::EMT);

// Logging
auto logger = DataLogger::make(simName);
Expand Down
8 changes: 1 addition & 7 deletions dpsim/examples/cxx/CIM/EMT_WSCC-9bus_VBR.cpp
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -76,13 +76,7 @@ int main(int argc, char *argv[]) {
CPS::CIM::Reader reader2(simName, Logger::Level::debug, Logger::Level::debug);
SystemTopology sys = reader2.loadCIM(60, filenames, Domain::EMT, PhaseType::ABC, CPS::GeneratorType::FullOrderVBR);

sys.initWithPowerflow(systemPF);
for (auto comp : sys.mComponents) {
if (auto genEMT = std::dynamic_pointer_cast<CPS::EMT::Ph3::SynchronGeneratorVBR>(comp)) {
auto genPF = systemPF.component<CPS::SP::Ph1::SynchronGenerator>(comp->name());
genEMT->terminal(0)->setPower(-genPF->getApparentPower());
}
}
sys.initWithPowerflow(systemPF, CPS::Domain::EMT);

// Logging
auto logger = DataLogger::make(simName);
Expand Down
4 changes: 1 addition & 3 deletions dpsim/examples/cxx/CIM/EMT_WSCC9bus_SGReducedOrderVBR.cpp
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -132,11 +132,9 @@ int main(int argc, char *argv[]) {
sys.addComponent(sw);
}

sys.initWithPowerflow(systemPF);
sys.initWithPowerflow(systemPF, CPS::Domain::EMT);
for (auto comp : sys.mComponents) {
if (auto genReducedOrder = std::dynamic_pointer_cast<CPS::Base::ReducedOrderSynchronGenerator<Real>>(comp)) {
auto genPF = systemPF.component<CPS::SP::Ph1::SynchronGenerator>(comp->name());
genReducedOrder->terminal(0)->setPower(-genPF->getApparentPower());
genReducedOrder->scaleInertiaConstant(inertiaScalingFactor);
genReducedOrder->setModelAsNortonSource(false);
}
Expand Down
7 changes: 2 additions & 5 deletions dpsim/examples/cxx/CIM/SP_WSCC9bus_SGReducedOrderVBR.cpp
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -129,13 +129,10 @@ int main(int argc, char *argv[]) {
sys.addComponent(faultSP);
}

sys.initWithPowerflow(systemPF);
sys.initWithPowerflow(systemPF, Domain::SP);
for (auto comp : sys.mComponents) {
if (auto genReducedOrder = std::dynamic_pointer_cast<CPS::Base::ReducedOrderSynchronGenerator<Complex>>(comp)) {
auto genPF = systemPF.component<CPS::SP::Ph1::SynchronGenerator>(comp->name());
genReducedOrder->terminal(0)->setPower(-genPF->getApparentPower());
if (auto genReducedOrder = std::dynamic_pointer_cast<CPS::Base::ReducedOrderSynchronGenerator<Complex>>(comp))
genReducedOrder->scaleInertiaConstant(inertiaScalingFactor);
}
}

// Logging
Expand Down
2 changes: 1 addition & 1 deletion dpsim/examples/cxx/CIM/WSCC-9bus_CIM.cpp
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -72,7 +72,7 @@ int main(int argc, char *argv[]) {

CPS::CIM::Reader reader2(simName, Logger::Level::debug, Logger::Level::off);
SystemTopology sys = reader2.loadCIM(60, filenames, Domain::DP, PhaseType::Single, CPS::GeneratorType::IdealVoltageSource);
sys.initWithPowerflow(systemPF);
sys.initWithPowerflow(systemPF, CPS::Domain::EMT);

// Logging
auto logger = DataLogger::make(simName);
Expand Down
2 changes: 1 addition & 1 deletion dpsim/examples/cxx/Circuits/DP_Slack_PiLine_PQLoad_with_PF_Init.cpp
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -111,7 +111,7 @@ int main(int argc, char* argv[]) {
SystemComponentList{extnetDP, lineDP, loadDP});

// Initialization of dynamic topology
systemDP.initWithPowerflow(systemPF);
systemDP.initWithPowerflow(systemPF, Domain::DP);

// Logging
auto loggerDP = DataLogger::make(simNameDP);
Expand Down
2 changes: 1 addition & 1 deletion dpsim/examples/cxx/Circuits/DP_Slack_PiLine_VSI_Ramp_with_PF_Init.cpp
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -132,7 +132,7 @@ int main(int argc, char* argv[]) {
SystemComponentList{extnetDP, lineDP, pv});

// Initialization of dynamic topology with values from powerflow
systemDP.initWithPowerflow(systemPF);
systemDP.initWithPowerflow(systemPF, Domain::DP);

// Logging
auto loggerDP = DataLogger::make(simName);
Expand Down
2 changes: 1 addition & 1 deletion dpsim/examples/cxx/Circuits/DP_Slack_PiLine_VSI_with_PF_Init.cpp
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -120,7 +120,7 @@ int main(int argc, char* argv[]) {
SystemComponentList{extnetDP, lineDP, pv});

// Initialization of dynamic topology with values from powerflow
systemDP.initWithPowerflow(systemPF);
systemDP.initWithPowerflow(systemPF, Domain::DP);

// Logging
auto loggerDP = DataLogger::make(simNameDP);
Expand Down
6 changes: 3 additions & 3 deletions dpsim/examples/cxx/Circuits/DP_SynGenTrStab_3Bus_Fault.cpp
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -32,13 +32,13 @@ void DP_SynGenTrStab_3Bus_Fault(String simName, Real timeStep, Real finalTime, b
auto n3PF = SimNode<Complex>::make("n3", PhaseType::Single);

//Synchronous generator 1
auto gen1PF = SP::Ph1::SynchronGenerator::make("Generator", Logger::Level::debug);
auto gen1PF = SP::Ph1::SynchronGenerator::make("SynGen1", Logger::Level::debug);
// setPointVoltage is defined as the voltage at the transfomer primary side and should be transformed to network side
gen1PF->setParameters(ThreeBus.nomPower_G1, ThreeBus.nomPhPhVoltRMS_G1, ThreeBus.initActivePower_G1, ThreeBus.setPointVoltage_G1*ThreeBus.t1_ratio, PowerflowBusType::VD);
gen1PF->setBaseVoltage(ThreeBus.Vnom);

//Synchronous generator 2
auto gen2PF = SP::Ph1::SynchronGenerator::make("Generator", Logger::Level::debug);
auto gen2PF = SP::Ph1::SynchronGenerator::make("SynGen2", Logger::Level::debug);
// setPointVoltage is defined as the voltage at the transfomer primary side and should be transformed to network side
gen2PF->setParameters(ThreeBus.nomPower_G2, ThreeBus.nomPhPhVoltRMS_G2, ThreeBus.initActivePower_G2, ThreeBus.setPointVoltage_G2*ThreeBus.t2_ratio, PowerflowBusType::PV);
gen2PF->setBaseVoltage(ThreeBus.Vnom);
Expand Down Expand Up @@ -166,7 +166,7 @@ void DP_SynGenTrStab_3Bus_Fault(String simName, Real timeStep, Real finalTime, b
SystemComponentList{gen1DP, gen2DP, loadDP, line12DP, line13DP, line23DP, faultDP});

// Initialization of dynamic topology
systemDP.initWithPowerflow(systemPF);
systemDP.initWithPowerflow(systemPF, Domain::DP);

// Logging
auto loggerDP = DataLogger::make(simNameDP);
Expand Down
Loading
Loading