DPsim::PFSolverPowerPolarSparse Class Reference

#include <PFSolverPowerPolarSparse.h>

Inheritance diagram for DPsim::PFSolverPowerPolarSparse:

Collaboration diagram for DPsim::PFSolverPowerPolarSparse:

Public Member Functions
	PFSolverPowerPolarSparse (CPS::String name, const CPS::SystemTopology &system, CPS::Real timeStep, CPS::Logger::Level logLevel)
	Constructor to be used in simulation examples.
Public Member Functions inherited from DPsim::PFSolverPowerPolar
	PFSolverPowerPolar (CPS::String name, const CPS::SystemTopology &system, CPS::Real timeStep, CPS::Logger::Level logLevel)
	Constructor to be used in simulation examples.
Public Member Functions inherited from DPsim::PFSolver
	PFSolver (CPS::String name, CPS::SystemTopology system, Real timeStep, CPS::Logger::Level logLevel)
	Constructor to be used in simulation examples.
void	setVDNode (CPS::String name)
	Set a node to VD using its name.
void	modifyPowerFlowBusComponent (CPS::String name, CPS::PowerflowBusType powerFlowBusType)
	Allows to modify the powerflow bus type of a specific component.
void	setSolverAndComponentBehaviour (Solver::Behaviour behaviour) override
	set solver and component to initialization or simulation behaviour
void	setKeepLastSolution (CPS::Bool keepLastSolution)
void	setEnforceReactiveLimits (CPS::Bool value)
	Enable generator reactive-limit enforcement (PV<->PQ outer loop)
void	setBaseVoltageLooseTolerance (CPS::Real tolerance)
	Raise for grids with legitimate large voltage drop (e.g. untapped feeders)
void	setBaseVoltageStrictTolerance (CPS::Real tolerance)
	Override the tolerance between authoritative base-voltage sources within a zone.
CPS::Bool	getKeepLastSolution () const
void	setBaseApparentPowerFallback (CPS::Real baseApparentPowerFallback)
CPS::Real	getBaseApparentPowerFallback () const
void	setMaxIterations (CPS::UInt maxIterations)
CPS::UInt	getMaxIterations () const
Public Member Functions inherited from DPsim::Solver
	Solver (String name, CPS::Logger::Level logLevel)
void	setTimeStep (Real timeStep)
void	doFrequencyParallelization (Bool freqParallel)
virtual void	setSystem (const CPS::SystemTopology &system)
void	doSystemMatrixRecomputation (Bool value)
void	setLogSolveTimes (Bool value)
void	doSteadyStateInit (Bool f)
	activate steady state initialization
void	setSteadStIniTimeLimit (Real v)
	set steady state initialization time limit
void	setSteadStIniAccLimit (Real v)
	set steady state initialization accuracy limit
void	doInitFromNodesAndTerminals (Bool f)
	activate powerflow initialization
virtual void	setDirectLinearSolverConfiguration (DirectLinearSolverConfiguration &)
	set direct linear solver configuration (only available in MNA for now)
virtual void	logLUTimes ()
	log LU decomposition times, if applicable
virtual void	log (Real time, Int timeStepCount)
	Log results.
void	setMaxNumberOfIterations (int maxIterations)

Protected Member Functions
void	setUpJacobianStorage () override
	Build the fixed sparsity pattern, create the linear solver and analyze the pattern once.
void	calculateJacobian () override
	Fill the sparse Jacobian values in place (pattern unchanged)
void	solveJacobianSystem () override
	Factorize (first iteration of a run) or refactorize, then solve.
void	buildJacobianPattern ()
	Construct the structural sparsity pattern of the Jacobian from the bus admittance matrix.
bool	isConnected (CPS::UInt k, CPS::UInt j)
	Whether buses k and j are connected (off-diagonal Jacobian entry exists)
Protected Member Functions inherited from DPsim::PFSolverPowerPolar
void	generateInitialSolution (Real time, bool keep_last_solution=false) override
	Generate initial solution for current time step.
void	updateSolution () override
	Update solution in each iteration.
void	setSolution () override
	Set final solution.
void	calculateMismatch () override
	Calculate mismatch.
void	resize_sol (CPS::Int n)
	Resize solution vector.
void	resize_complex_sol (CPS::Int n)
	Resize complex solution vector.
CPS::Real	sol_Vr (CPS::UInt k)
	Calculate real part of voltage from sol_V and sol_D.
CPS::Real	sol_Vi (CPS::UInt k)
	Calculate imaginary part of voltage from sol_V and sol_D.
CPS::Complex	sol_Vcx (CPS::UInt k)
	Calculate complex voltage from sol_V and sol_D.
CPS::Real	P (CPS::UInt k)
	Calculate active power at a bus from current solution.
CPS::Real	Q (CPS::UInt k)
	Calculate the reactive power at a bus from current solution.
void	calculatePAndQAtSlackBus ()
	Calculate P and Q at slack bus from current solution.
void	calculateQAtPVBuses ()
	Calculate the reactive power at all PV buses from current solution.
CPS::Real	generatorReactivePowerPerUnit (CPS::TopologicalNode::Ptr node)
	Generator reactive-power injection at a bus [pu], used by the Q-limit check.
CPS::Real	loadReactivePowerPerUnit (CPS::TopologicalNode::Ptr node)
	Total local load reactive power at a bus [pu].
CPS::Bool	enforceReactiveLimits () override
	Q-limit PV<->PQ switching pass (overrides the base no-op)
void	clearReactiveLimitState () override
	Clear the Q-limit conversion bookkeeping between solves.
void	calculatePAndQInjectionPQBuses ()
	Calculate complex power flowing from this node to the other nodes.
void	calculateBranchFlow ()
	Calculate branch flows from current solution and store them in line and transformer components.
void	calculateNodalInjection ()
	Calculate nodal power injections and store them in first line or transformer (in case no line is connected), so that lower level classes (Node, TopologicalTerminal) can stay untouched.
Protected Member Functions inherited from DPsim::PFSolver
void	initialize () override
	Initialization of the solver.
void	initializeComponents ()
	Initialization of individual components.
void	assignMatrixNodeIndices ()
	Assignment of matrix indices for nodes.
void	setBaseApparentPower ()
	Set apparent base power of per-unit system.
void	determinePFBusType ()
	Determine bus type for all buses.
void	rebuildBusIndexAggregates ()
	Rebuild index vectors + counts from the PQ/PV/VD node lists.
void	reclassifyBuses ()
	Re-derive index vectors and resize the system after PV<->PQ switching.
void	resetToOriginalClassification ()
	Restore the original PV/PQ classification before a fresh solve.
CPS::Real	componentBaseVoltage (CPS::TopologicalPowerComp::Ptr comp, CPS::TopologicalNode::Ptr node)
	Base voltage a single component reports for node, or 0 if unknown.
void	propagateAndVerifyBaseVoltage ()
	Determine, verify and propagate each node's base voltage per electrical zone.
void	composeAdmittanceMatrix ()
	Compose admittance matrix.
CPS::Real	G (int i, int j)
	Gets the real part of admittance matrix element.
CPS::Real	B (int i, int j)
	Gets the imaginary part of admittance matrix element.
Bool	solvePowerflow ()
	Solves the powerflow problem.
Bool	runNewtonRaphson ()
	Run a single Newton-Raphson solve with the current bus classification.
CPS::Bool	checkConvergence ()
	Check whether below tolerance.
CPS::String	logVector (std::vector< CPS::UInt > indexVector)
	Logging for integer vectors.
CPS::Task::List	getTasks () override
	Get tasks for scheduler.

Protected Attributes
CPS::SparseMatrixRow	mJsparse
	Sparse Jacobian (row-major, fixed pattern, values updated in place each iteration)
std::shared_ptr< DirectLinearSolver >	mLinearSolver
	Direct linear solver reused across iterations (symbolic factorization computed once)
std::vector< std::pair< CPS::UInt, CPS::UInt > >	mVariableSystemMatrixEntries
	Empty list: PF uses full refactorization, not partial refactorization.
Protected Attributes inherited from DPsim::PFSolverPowerPolar
CPS::Vector	sol_P
	Solution vector of active power.
CPS::Vector	sol_Q
	Solution vector of reactive power.
CPS::Vector	sol_V
	Solution vector of voltage magnitude.
CPS::Vector	sol_D
	Solution vector of voltage angle.
CPS::VectorComp	sol_V_complex
	Solution vector of representing sol_V and sol_D as complex quantity.
CPS::VectorComp	sol_S_complex
	Solution vector of representing sol_P and sol_Q as complex quantity.
CPS::Vector	Pesp
CPS::Vector	Qesp
CPS::Vector	mLastConvergedV
CPS::Vector	mLastConvergedD
bool	mHasLastConvergedSolution = false
std::map< CPS::TopologicalNode::Ptr, bool >	mQLimitConvertedAtMax
	Buses switched PV->PQ by the Q-limit loop -> pinned at Qmax (true) or Qmin (false)
std::map< CPS::TopologicalNode::Ptr, CPS::UInt >	mQLimitSwitchCount
	Total PV<->PQ switches per bus over the solve (anti-oscillation cap)
Protected Attributes inherited from DPsim::PFSolver
UInt	mNumPQBuses = 0
	Number of PQ nodes.
UInt	mNumPVBuses = 0
	Number of PV nodes.
UInt	mNumVDBuses = 0
	Number of PV nodes.
UInt	mNumUnknowns = 0
	Number of unknowns, defining system dimension.
CPS::TopologicalNode::List	mPQBuses
	Vector of nodes characterized as PQ buses.
CPS::TopologicalNode::List	mPVBuses
	Vector of nodes characterized as PV buses.
CPS::TopologicalNode::List	mVDBuses
	Vector of nodes characterized as VD buses.
CPS::TopologicalNode::List	mPQBusesOrig
	Original PQ/PV classification (snapshot before Q-limit switching)
CPS::TopologicalNode::List	mPVBusesOrig
std::vector< CPS::UInt >	mPQBusIndices
	Vector with indices of PQ buses.
std::vector< CPS::UInt >	mPVBusIndices
	Vector with indices of PV buses.
std::vector< CPS::UInt >	mVDBusIndices
	Vector with indices of VD buses.
std::vector< CPS::UInt >	mPQPVBusIndices
	Vector with indices of both PQ and PV buses.
CPS::SparseMatrixCompRow	mY
	Admittance matrix.
CPS::Matrix	mJ
	Jacobian matrix.
CPS::Vector	mX
	Solution vector.
CPS::Vector	mF
	Vector of mismatch values.
CPS::SystemTopology	mSystem
	System list.
std::vector< std::shared_ptr< CPS::SP::Ph1::Transformer > >	mTransformers
	Vector of transformer components.
std::vector< std::shared_ptr< CPS::SP::Ph1::SolidStateTransformer > >	mSolidStateTransformers
	Vector of solid state transformer components.
std::vector< std::shared_ptr< CPS::SP::Ph1::SynchronGenerator > >	mSynchronGenerators
	Vector of synchronous generator components.
std::vector< std::shared_ptr< CPS::SP::Ph1::Load > >	mLoads
	Vector of load components.
std::vector< std::shared_ptr< CPS::SP::Ph1::PiLine > >	mLines
	Vector of line components.
std::vector< std::shared_ptr< CPS::SP::Ph1::Shunt > >	mShunts
	Vector of shunt components.
std::vector< std::shared_ptr< CPS::SP::Ph1::NetworkInjection > >	mExternalGrids
	Vector of external grid components.
std::vector< std::shared_ptr< CPS::SP::Ph1::AvVoltageSourceInverterDQ > >	mAverageVoltageSourceInverters
	Vector of average voltage source inverters.
std::map< CPS::TopologicalNode::Ptr, CPS::Real >	mBaseVoltageAtNode
	Map providing determined base voltages for each node.
Real	mTolerance = 1e-8
	Solver tolerance.
CPS::UInt	mMaxIterations = 20
	Maximum number of iterations.
CPS::UInt	mIterations
	Actual number of iterations.
CPS::Bool	mEnforceReactiveLimits = false
	Enforce generator reactive-power limits via PV<->PQ outer-loop switching.
CPS::UInt	mMaxOuterIterations = 10
	Maximum number of Q-limit outer iterations.
CPS::UInt	mMaxQLimitSwitchesPerBus = 2
	Maximum number of PV<->PQ switches per bus before it is frozen (anti-oscillation)
CPS::Real	mBaseVoltageLooseTolerance = 0.1
	Relative tolerance for non-authoritative (e.g. Load) base-voltage candidates vs. the zone's rating.
CPS::Real	mBaseVoltageStrictTolerance = 0.01
	Relative tolerance between authoritative base-voltage sources (generator/transformer/network-injection/VSI) within a zone.
CPS::Real	mBaseApparentPower
	Base power of per-unit system.
CPS::Real	mBaseApparentPowerFallback = 100e6
	Fallback base apparent power if no generator or transformer rating is found.
CPS::Bool	isConverged = false
	Convergence flag.
CPS::Bool	solutionInitialized = false
	Flag whether solution vectors are initialized.
CPS::Bool	solutionComplexInitialized = false
	Flag whether complex solution vectors are initialized.
CPS::Bool	mKeepLastSolution = false
	Use last converged solution as initial guess.
Protected Attributes inherited from DPsim::Solver
String	mName
	Name for logging.
CPS::Logger::Level	mLogLevel
	Logging level.
Bool	mLogSolveTimes = true
	Collect step time for logging.
CPS::Logger::Log	mSLog
	Logger.
Real	mTimeStep
	Time step for fixed step solvers.
Bool	mFrequencyParallel = false
	Activates parallelized computation of frequencies.
Real	mSteadStIniTimeLimit = 10
	steady state initialization time limit
Real	mSteadStIniAccLimit = 0.0001
	steady state initialization accuracy limit
Bool	mSteadyStateInit = false
	Activates steady state initialization.
Bool	mIsInInitialization = false
	Determines if solver is in initialization phase, which requires different behavior.
Bool	mInitFromNodesAndTerminals = true
Bool	mSystemMatrixRecomputation = false
	Enable recomputation of system matrix during simulation.
Behaviour	mBehaviour = Solver::Behaviour::Simulation
	Solver behaviour initialization or simulation.

Additional Inherited Members
Public Types inherited from DPsim::Solver
enum	Behaviour { Initialization , Simulation }
enum class	Type { MNA , DAE , NRP }
typedef std::shared_ptr< Solver >	Ptr
typedef std::vector< Ptr >	List
Public Attributes inherited from DPsim::Solver
int	mMaxIterations = 10

Detailed Description

Polar-coordinate powerflow solver that assembles a sparse Jacobian and reuses the symbolic factorization across Newton-Raphson iterations (scales to large grids).

Definition at line 12 of file PFSolverPowerPolarSparse.h.

Constructor & Destructor Documentation

PFSolverPowerPolarSparse()

PFSolverPowerPolarSparse::PFSolverPowerPolarSparse	(	CPS::String	name,
		const CPS::SystemTopology &	system,
		CPS::Real	timeStep,
		CPS::Logger::Level	logLevel )

Constructor to be used in simulation examples.

Definition at line 17 of file PFSolverPowerPolarSparse.cpp.

~PFSolverPowerPolarSparse()

virtual DPsim::PFSolverPowerPolarSparse::~PFSolverPowerPolarSparse ( )

inlinevirtual

Definition at line 38 of file PFSolverPowerPolarSparse.h.

Member Function Documentation

buildJacobianPattern()

void PFSolverPowerPolarSparse::buildJacobianPattern ( )

protected

Construct the structural sparsity pattern of the Jacobian from the bus admittance matrix.

Definition at line 26 of file PFSolverPowerPolarSparse.cpp.

calculateJacobian()

void PFSolverPowerPolarSparse::calculateJacobian ( )

overrideprotectedvirtual

Fill the sparse Jacobian values in place (pattern unchanged)

Reimplemented from DPsim::PFSolverPowerPolar.

Definition at line 84 of file PFSolverPowerPolarSparse.cpp.

isConnected()

bool PFSolverPowerPolarSparse::isConnected ( CPS::UInt k, CPS::UInt j )

protected

Whether buses k and j are connected (off-diagonal Jacobian entry exists)

Definition at line 22 of file PFSolverPowerPolarSparse.cpp.

setUpJacobianStorage()

void PFSolverPowerPolarSparse::setUpJacobianStorage ( )

overrideprotectedvirtual

Build the fixed sparsity pattern, create the linear solver and analyze the pattern once.

Reimplemented from DPsim::PFSolver.

Definition at line 66 of file PFSolverPowerPolarSparse.cpp.

solveJacobianSystem()

void PFSolverPowerPolarSparse::solveJacobianSystem ( )

overrideprotectedvirtual

Factorize (first iteration of a run) or refactorize, then solve.

Reimplemented from DPsim::PFSolver.

Definition at line 165 of file PFSolverPowerPolarSparse.cpp.

Member Data Documentation

mJsparse

CPS::SparseMatrixRow DPsim::PFSolverPowerPolarSparse::mJsparse

protected

Sparse Jacobian (row-major, fixed pattern, values updated in place each iteration)

Definition at line 15 of file PFSolverPowerPolarSparse.h.

mLinearSolver

std::shared_ptr<DirectLinearSolver> DPsim::PFSolverPowerPolarSparse::mLinearSolver

protected

Direct linear solver reused across iterations (symbolic factorization computed once)

Definition at line 17 of file PFSolverPowerPolarSparse.h.

mVariableSystemMatrixEntries

std::vector<std::pair<CPS::UInt, CPS::UInt> > DPsim::PFSolverPowerPolarSparse::mVariableSystemMatrixEntries

protected

Empty list: PF uses full refactorization, not partial refactorization.

Definition at line 19 of file PFSolverPowerPolarSparse.h.

---

The documentation for this class was generated from the following files:

• /__w/dpsim/dpsim/dpsim/include/dpsim/PFSolverPowerPolarSparse.h
• /__w/dpsim/dpsim/dpsim/src/PFSolverPowerPolarSparse.cpp