Stage 1.0.05 - Solver NR development. Part 7. Stable NR power flow.

ealux · May 6, 2022 · 4fbb0e8 · 4fbb0e8
2 parents 1a62c98 + d8e0a81
commit 4fbb0e8
Show file tree

Hide file tree

Showing 5 changed files with 386 additions and 162 deletions.
diff --git a/PowerFlowCore/Engine/Engine.cs b/PowerFlowCore/Engine/Engine.cs
@@ -71,12 +71,9 @@ public void Calculate()
         {
             //this.Grid.SolverGS(this.Grid.Uinit, this.Options);
 
-            //this.Grid.SolverGS(this.Grid.Uinit, new CalculationOptions() { IterationsCount = 10 }).SolverNR(this.Grid.Ucalc, this.Options);
+            //this.Grid.SolverGS(this.Grid.Uinit, new CalculationOptions() { IterationsCount = 15 }).SolverNR(this.Grid.Ucalc, this.Options);
             this.Grid.SolverNR(this.Grid.Uinit, this.Options);
 
-            //this.Grid.SolverNR2(this.Grid.Uinit, this.Options);
-            //this.Grid.SolverGS(this.Grid.Uinit, new CalculationOptions() { IterationsCount = 5}).SolverNR2(this.Grid.Ucalc, this.Options);
-
 
             this.Grid.CalculatePowerMatrix();
             this.NeedsToCalc = false;

diff --git a/PowerFlowCore/Options/CalculationOptions.cs b/PowerFlowCore/Options/CalculationOptions.cs
@@ -14,7 +14,7 @@ public sealed class CalculationOptions
         ///   Gauss-Seidel:   Voltage step defference;
         ///   Newton-Raphson: Power residual tolerance.
         /// </summary>
-        public double Accuracy { get; set; } = 1e-6;
+        public double Accuracy { get; set; } = 1e-8;
 
 
         /// <summary>
@@ -32,7 +32,7 @@ public sealed class CalculationOptions
         /// <summary>
         /// Voltage convergence criteria (Newton-Raphson method only)
         /// </summary>
-        public double VoltageConvergence { get; set; } = 1e-6;
+        public double VoltageConvergence { get; set; } = 1e-8;
 
 
         /// <summary>

diff --git a/PowerFlowCore/Solvers/SolverGS.cs b/PowerFlowCore/Solvers/SolverGS.cs
@@ -83,6 +83,14 @@ public static Grid SolverGS(this Grid grid,
             for (int n = 0; n < grid.Nodes.Count; n++) 
                 grid.Nodes[n].U = U[n];
 
+
+            // !!!!!
+            // TODO:
+            // 1. (Not)Success logic
+            // 2. Voltage restrictions
+            // 3. Angle restrictions in branches
+            // !!!!!
+
             return grid;
         }