Fix power flow initialisation of synchronous generators (#238)

In powerflow simulations, the power of synchronous generators is not being always correctly calculated, since the power of the loads connected to nodes is not considered in the calculation. This PR fixes the calculation of the synchronous generators' power in powerflow simulations and adapts the funcion `initWithPowerFlow` to initialialize the terminal power of the synchronous generators after the powerflow calculation. Closes #223
sogno-platform · Nov 30, 2023 · 4b9ab5f · 4b9ab5f
2 parents dfd811d + 9405f52
commit 4b9ab5f
Show file tree

Hide file tree

Showing 85 changed files with 1,436 additions and 1,124 deletions.
diff --git a/dpsim-models/include/dpsim-models/SP/SP_Ph1_Shunt.h b/dpsim-models/include/dpsim-models/SP/SP_Ph1_Shunt.h
@@ -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);

diff --git a/dpsim-models/include/dpsim-models/SimNode.h b/dpsim-models/include/dpsim-models/SimNode.h
@@ -127,4 +127,7 @@ namespace CPS {
 
 	template<>
 	void SimNode<Complex>::setVoltage(Complex newVoltage);
+
+	template<>
+	void SimNode<Complex>::setPower(Complex newPower);
 }
diff --git a/dpsim-models/include/dpsim-models/SystemTopology.h b/dpsim-models/include/dpsim-models/SystemTopology.h
@@ -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>
@@ -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);

diff --git a/dpsim-models/src/SP/SP_Ph1_Shunt.cpp b/dpsim-models/src/SP/SP_Ph1_Shunt.cpp
@@ -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);
@@ -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;

diff --git a/dpsim-models/src/SystemTopology.cpp b/dpsim-models/src/SystemTopology.cpp
@@ -12,6 +12,7 @@
 #include <unordered_map>
 
 #include <dpsim-models/SystemTopology.h>
+#include <dpsim-models/SP/SP_Ph1_SynchronGenerator.h>
 
 using namespace CPS;
 
@@ -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());
@@ -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) {

diff --git a/dpsim/examples/cxx/CIM/DP_CIGRE_MV_withDG.cpp b/dpsim/examples/cxx/CIM/DP_CIGRE_MV_withDG.cpp
@@ -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);
 

diff --git a/dpsim/examples/cxx/CIM/DP_CIGRE_MV_withDG_withLoadStep.cpp b/dpsim/examples/cxx/CIM/DP_CIGRE_MV_withDG_withLoadStep.cpp
@@ -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) {

diff --git a/dpsim/examples/cxx/CIM/DP_CIGRE_MV_withoutDG.cpp b/dpsim/examples/cxx/CIM/DP_CIGRE_MV_withoutDG.cpp
@@ -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);
 

diff --git a/dpsim/examples/cxx/CIM/DP_WSCC-9bus_IdealVS.cpp b/dpsim/examples/cxx/CIM/DP_WSCC-9bus_IdealVS.cpp
@@ -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)
@@ -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);

diff --git a/dpsim/examples/cxx/CIM/DP_WSCC9bus_SGReducedOrderIter.cpp b/dpsim/examples/cxx/CIM/DP_WSCC9bus_SGReducedOrderIter.cpp
@@ -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

diff --git a/dpsim/examples/cxx/CIM/DP_WSCC9bus_SGReducedOrderVBR.cpp b/dpsim/examples/cxx/CIM/DP_WSCC9bus_SGReducedOrderVBR.cpp
@@ -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);
 		}

diff --git a/dpsim/examples/cxx/CIM/EMT_CIGRE_MV_withDG.cpp b/dpsim/examples/cxx/CIM/EMT_CIGRE_MV_withDG.cpp
@@ -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);
 

diff --git a/dpsim/examples/cxx/CIM/EMT_CIGRE_MV_withDG_withLoadStep.cpp b/dpsim/examples/cxx/CIM/EMT_CIGRE_MV_withDG_withLoadStep.cpp
@@ -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);
 

diff --git a/dpsim/examples/cxx/CIM/EMT_CIGRE_MV_withoutDG.cpp b/dpsim/examples/cxx/CIM/EMT_CIGRE_MV_withoutDG.cpp
@@ -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);
 

diff --git a/dpsim/examples/cxx/CIM/EMT_WSCC-9bus_FullOrderSG.cpp b/dpsim/examples/cxx/CIM/EMT_WSCC-9bus_FullOrderSG.cpp
@@ -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);

diff --git a/dpsim/examples/cxx/CIM/EMT_WSCC-9bus_IdealCS.cpp b/dpsim/examples/cxx/CIM/EMT_WSCC-9bus_IdealCS.cpp
@@ -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);

diff --git a/dpsim/examples/cxx/CIM/EMT_WSCC-9bus_IdealVS.cpp b/dpsim/examples/cxx/CIM/EMT_WSCC-9bus_IdealVS.cpp
@@ -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);

diff --git a/dpsim/examples/cxx/CIM/EMT_WSCC-9bus_VBR.cpp b/dpsim/examples/cxx/CIM/EMT_WSCC-9bus_VBR.cpp
@@ -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);

diff --git a/dpsim/examples/cxx/CIM/EMT_WSCC9bus_SGReducedOrderVBR.cpp b/dpsim/examples/cxx/CIM/EMT_WSCC9bus_SGReducedOrderVBR.cpp
@@ -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);
 		}

diff --git a/dpsim/examples/cxx/CIM/SP_WSCC9bus_SGReducedOrderVBR.cpp b/dpsim/examples/cxx/CIM/SP_WSCC9bus_SGReducedOrderVBR.cpp
@@ -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

diff --git a/dpsim/examples/cxx/CIM/WSCC-9bus_CIM.cpp b/dpsim/examples/cxx/CIM/WSCC-9bus_CIM.cpp
@@ -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);

diff --git a/dpsim/examples/cxx/Circuits/DP_Slack_PiLine_PQLoad_with_PF_Init.cpp b/dpsim/examples/cxx/Circuits/DP_Slack_PiLine_PQLoad_with_PF_Init.cpp
@@ -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);

diff --git a/dpsim/examples/cxx/Circuits/DP_Slack_PiLine_VSI_Ramp_with_PF_Init.cpp b/dpsim/examples/cxx/Circuits/DP_Slack_PiLine_VSI_Ramp_with_PF_Init.cpp
@@ -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);

diff --git a/dpsim/examples/cxx/Circuits/DP_Slack_PiLine_VSI_with_PF_Init.cpp b/dpsim/examples/cxx/Circuits/DP_Slack_PiLine_VSI_with_PF_Init.cpp
@@ -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);

diff --git a/dpsim/examples/cxx/Circuits/DP_SynGenTrStab_3Bus_Fault.cpp b/dpsim/examples/cxx/Circuits/DP_SynGenTrStab_3Bus_Fault.cpp
@@ -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);
@@ -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);